Search

While polygenic risk scores (PRSs) enable early identification of genetic risk for chronic obstructive pulmonary disease (COPD), predictive performance is limited when the discovery and target populations are not well matched. Hypothesizing that the biological mechanisms of disease are shared across ancestry groups, we introduce a PrediXcan-derived polygenic transcriptome risk score (PTRS) to improve cross-ethnic portability of risk prediction. We constructed the PTRS using summary statistics from application of PrediXcan on large-scale GWASs of lung function (forced expiratory volume in 1 s [FEV(1)] and its ratio to forced vital capacity [FEV(1)/FVC]) in the UK Biobank. We examined prediction performance and cross-ethnic portability of PTRS through smoking-stratified analyses both on 29,381 multi-ethnic participants from TOPMed population/family-based cohorts and on 11,771 multi-ethnic participants from TOPMed COPD-enriched studies. Analyses were carried out for two dichotomous COPD traits (moderate-to-severe and severe COPD) and two quantitative lung function traits (FEV(1) and FEV(1)/FVC). While the proposed PTRS showed weaker associations with disease than PRS for European ancestry, the PTRS showed stronger association with COPD than PRS for African Americans (e.g., odds ratio [OR] = 1.24 [95% confidence interval [CI]: 1.08–1.43] for PTRS versus 1.10 [0.96–1.26] for PRS among heavy smokers with ≥ 40 pack-years of smoking) for moderate-to-severe COPD. Cross-ethnic portability of the PTRS was significantly higher than the PRS (paired t test p

Background: Plasma proteins are critical mediators of cardiovascular processes and are the targets of many drugs. Previous efforts to characterize the genetic architecture of the plasma proteome have been limited by a focus on individuals of European descent and leveraged genotyping arrays and imputation. Here we describe whole genome sequence analysis of the plasma proteome in individuals with greater African ancestry, increasing our power to identify novel genetic determinants. Methods: Proteomic profiling of 1301 proteins was performed in 1852 Black adults from the Jackson Heart Study using aptamer-based proteomics (SomaScan). Whole genome sequencing association analysis was ascertained for all variants with minor allele count ≥5. Results were validated using an alternative, antibody-based, proteomic platform (Olink) as well as replicated in the Multi-Ethnic Study of Atherosclerosis and the HERITAGE Family Study (Health, Risk Factors, Exercise Training and Genetics). Results: We identify 569 genetic associations between 479 proteins and 438 unique genetic regions at a Bonferroni-adjusted significance level of 3.8×10 -11 . These associations include 114 novel locus-protein relationships and an additional 217 novel sentinel variant-protein relationships. Novel cardiovascular findings include new protein associations at the APOE gene locus including ZAP70 (sentinel single nucleotide polymorphism [SNP] rs7412-T, β=0.61±0.05, P =3.27×10 -30 ) and MMP-3 (β=-0.60±0.05, P =1.67×10 -32 ), as well as a completely novel pleiotropic locus at the HPX gene, associated with 9 proteins. Further, the associations suggest new mechanisms of genetically mediated cardiovascular disease linked to African ancestry; we identify a novel association between variants linked to APOL1-associated chronic kidney and heart disease and the protein CKAP2 (rs73885319-G, β=0.34±0.04, P =1.34×10 -17 ) as well as an association between ATTR amyloidosis and RBP4 levels in community-dwelling individuals without heart failure. Conclusions: Taken together, these results provide evidence for the functional importance of variants in non-European populations, and suggest new biological mechanisms for ancestry-specific determinants of lipids, coagulation, and myocardial function.

Background: Sleep problems (SP) are common in cancer patients but have not been previously assessed in patients receiving immune checkpoint inhibitors (ICI).Methods: We collected questionnaire data on sleep apnea risk, insomnia and general sleep patterns. We used an adjusted multivariate Poisson regression to calculate prevalence ratios (PRs) and associated 95% confidence intervals (CIs) for associations between these SP and metastatic versus localized cancer stage (M1 vs M0), and adjusted logistic regression models to calculate ORs for associations between SP with the number of ICI infusions completed (6+ vs. Results: Among 32 patients who received ICI treatment, the prevalence of low, intermediate, and high-risk OSA risk was 36%, 42%, and 21%, respectively. Overall, 58% of participants reported clinically significant insomnia. We did not find a significant association between intermediate or high risk OSA (vs. low risk) and metastatic cancer status [PR=1.01 (95% CI: 0.28, 3.67)]. Patients in the cohort who reported taking >15 minutes to fall asleep were 3.6 times more likely to be diagnosed with metastatic cancer compared to those reporting shorter sleep latency [95% CI (1.74, 7.35)]. We did not find a significant association between SP and number of ICI infusions completed.Conclusion: Our data associating sleep apnea risk, insomnia and sleep patterns with more advanced cancer encourages further exploration in larger-scale observational studies and suggests interventional clinical trials focused on sleep quality improvement could result in better outcomes for these patients.

Sex differences in Parkinson Disease (PD) risk are well-known. However, it is still unclear the role of sex chromosomes in the development and progression of PD. We performed the first X-chromosome Wide Association Study (XWAS) for PD risk in Latin American individuals. We used data from three admixed cohorts: (i) Latin American Research consortium on the GEnetics of Parkinson’s Disease (n=1,504) as discover cohort and (ii) Latino cohort from International Parkinson Disease Genomics Consortium (n = 155) and (iii) Bambui Aging cohort (n= 1,442) as replication cohorts. After developing a X-chromosome framework specifically designed for admixed populations, we identified eight linkage disequilibrium regions associated with PD. We fully replicated one of these regions (top variant rs525496; discovery OR [95%CI]: 0.60 [0.478 - 0.77], p = 3.13 ×10-5; replication OR: 0.60 [0.37-0.98], p = 0.04). rs525496 is an expression quantitative trait loci for several genes expressed in brain tissues, includingRAB9B, H2BFM, TSMB15BandGLRA4. We also replicated a previous XWAS finding (rs28602900), showing that this variant is associated with PD in non-European populations. Our results reinforce the importance of including X-chromosome and diverse populations in genetic studies.

Alzheimer disease (AD) is the most common form of senile dementia, with high incidence late in life in many populations including Caribbean Hispanic (CH) populations. Such admixed populations, descended from more than one ancestral population, can present challenges for genetic studies, including limited sample sizes and unique analytical constraints. Therefore, CH populations and other admixed populations have not been well represented in studies of AD, and much of the genetic variation contributing to AD risk in these populations remains unknown. Here, we conduct genome-wide analysis of AD in multiplex CH families from the Alzheimer Disease Sequencing Project (ADSP). We developed, validated, and applied an implementation of a logistic mixed model for admixture mapping with binary traits that leverages genetic ancestry to identify ancestry-of-origin loci contributing to AD. We identified three loci on chromosome 13q33.3 associated with reduced risk of AD, where associations were driven by Native American (NAM) ancestry. This AD admixture mapping signal spans the FAM155A, ABHD13, TNFSF13B, LIG4, and MYO16 genes and was supported by evidence for association in an independent sample from the Alzheimer's Genetics in Argentina-Alzheimer Argentina consortium (AGA-ALZAR) study with considerable NAM ancestry. We also provide evidence of NAM haplotypes and key variants within 13q33.3 that segregate with AD in the ADSP whole-genome sequencing data. Interestingly, the widely used genome-wide association study approach failed to identify associations in this region. Our findings underscore the potential of leveraging genetic ancestry diversity in recently admixed populations to improve genetic mapping, in this case for AD-relevant loci.

Many common and rare variants associated with hematologic traits have been discovered through imputation on large-scale reference panels. However, the majority of genome-wide association studies (GWASs) have been conducted in Europeans, and determining causal variants has proved challenging. We performed a GWAS of total leukocyte, neutrophil, lymphocyte, monocyte, eosinophil, and basophil counts generated from 109,563,748 variants in the autosomes and the X chromosome in the Trans-Omics for Precision Medicine (TOPMed) program, which included data from 61,802 individuals of diverse ancestry. We discovered and replicated 7 leukocyte trait associations, including (1) the association between a chromosome X, pseudo-autosomal region (PAR), noncoding variant located between cytokine receptor genes (CSF2RA and CLRF2) and lower eosinophil count; and (2) associations between single variants found predominantly among African Americans at the S1PR3 (9q22.1) and HBB (11p15.4) loci and monocyte and lymphocyte counts, respectively. We further provide evidence indicating that the newly discovered eosinophil-lowering chromosome X PAR variant might be associated with reduced susceptibility to common allergic diseases such as atopic dermatitis and asthma. Additionally, we found a burden of very rare FLT3 (13q12.2) variants associated with monocyte counts. Together, these results emphasize the utility of whole-genome sequencing in diverse samples in identifying associations missed by European-ancestry-driven GWASs.

Prevalence of smoking is higher in Alaska Native and American Indian (ANAI) populations living in Alaska than the general US population. Genetic factors contribute to smoking and cessation rates. The objective of this study was to compare CYP2A6 genetic variation and CYP2A6 enzyme activity toward nicotine in an ANAI population. ANAI (N = 151) people trying to quit smoking were recruited. DNA samples were genotyped for CYP2A6 variants *1X2A, *1B, *2, *4, *9, *10, *12, and *35. Multiple nicotine metabolites were measured in plasma and urine samples, including cotinine and 3′‐hydroxycotinine used to determine CYP2A6 activity (e.g., nicotine metabolite ratio [NMR]). We calculated summary statistics for all of the genotypes and metabolites and assigned CYP2A6 activity scores based on known information. We studied the association of CYP2A6 variants with the NMR and smoking histories. The overall frequency of the CYP2A6*1B gain of function allele was high in the ANAI versus non‐ANAI populations in other studies. Both *4 null and *9 decrease of function alleles had frequencies similar to previous studies of ANAI populations. In a multivariate analysis, the genotype‐inferred CYP2A6 activity score was associated with both plasma and urine NMR (p value = 8.56E‐08 and 4.08E‐13, respectively). Plasma NMR was also associated with duration of smoking (p value

Platelets play a key role in thrombosis and hemostasis. Platelet count (PLT) and mean platelet volume (MPV) are highly heritable quantitative traits, with hundreds of genetic signals previously identified, mostly in European ancestry populations. We here utilize whole genome sequencing (WGS) from NHLBI’s Trans-Omics for Precision Medicine initiative (TOPMed) in a large multi-ethnic sample to further explore common and rare variation contributing to PLT (n = 61 200) and MPV (n = 23 485). We identified and replicated secondary signals at MPL (rs532784633) and PECAM1 (rs73345162), both more common in African ancestry populations. We also observed rare variation in Mendelian platelet-related disorder genes influencing variation in platelet traits in TOPMed cohorts (not enriched for blood disorders). For example, association of GP9 with lower PLT and higher MPV was partly driven by a pathogenic Bernard-Soulier syndrome variant (rs5030764, p.Asn61Ser), and the signals at TUBB1 and CD36 were partly driven by loss of function variants not annotated as pathogenic in ClinVar (rs199948010 and rs571975065). However, residual signal remained for these gene-based signals after adjusting for lead variants, suggesting that additional variants in Mendelian genes with impacts in general population cohorts remain to be identified. Gene-based signals were also identified at several genome-wide association study identified loci for genes not annotated for Mendelian platelet disorders (PTPRH, TET2, CHEK2), with somatic variation driving the result at TET2. These results highlight the value of WGS in populations of diverse genetic ancestry to identify novel regulatory and coding signals, even for well-studied traits like platelet traits.

Alaska Native people are under-represented in genetic research but have unique gene variation that may critically impact their response to pharmacotherapy. Full resequencing of CYP2C9 in a cross-section of this population identified CYP2C9 Met1Leu (M1L), a novel, relatively common single nucleotide polymorphism hypothesized to confer CYP2C9 poor metabolizer phenotype by disrupting the start codon. M1L is present at a minor allele frequency of 6.3% in Yup’ik Alaska Native people and thus can contribute to the risk of an adverse drug response from narrow-therapeutic-index CYP2C9 substrates such as (S)-warfarin. This study’s objective was to characterize the catalytic efficiency of the Leu1 variant enzyme in vivo by evaluating the pharmacokinetic behavior of naproxen, a probe substrate for CYP2C9 activity, in genotyped Yup’ik participants. We first confirmed the selectivity of (S)-naproxen O-demethylation by CYP2C9 using activity-phenotyped human liver microsomes and selective cytochrome P450 inhibitors and then developed and validated a novel liquid chromatography mass spectrometry method for simultaneous quantification of (S)-naproxen, (S)-O-desmethylnaproxen, and naproxen acyl glucuronide in human urine. The average ratio of (S)-O-desmethylnaproxen to unchanged (S)-naproxen in urine was 18.0 ± 8.0 (n = 11) for the homozygous CYP2C9Met1 reference group and 10.3 ± 6.6 (n = 11) for the Leu1 variant carrier group (P = 0.011). The effect of M1L variation on CYP2C9 function and its potential to alter the pharmacokinetics of drugs metabolized by the enzyme has clinical implications and should be included in a variant screening panel when pharmacogenetic testing in the Alaska Native population is warranted. Significance Statement The novel CYP2C9 Met1Leu variant in Alaska Native people was recently identified. This study validated (S)-naproxen as a CYP2C9 probe substrate to characterize the in vivo functional activity of the CYP2C9 Leu1 variant. The results of this pharmacogenetic-pharmacokinetic study suggest that the CYP2C9 Leu1 variant exhibits loss of enzyme activity. This finding may be important to consider when administering narrow-therapeutic-index medications metabolized by CYP2C9 and also compels further investigation to characterize novel genetic variation in understudied populations.

The Trans-Omics for Precision Medicine (TOPMed) programme seeks to elucidate the genetic architecture and biology of heart, lung, blood and sleep disorders, with the ultimate goal of improving diagnosis, treatment and prevention of these diseases. The initial phases of the programme focused on whole-genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here we describe the TOPMed goals and design as well as the available resources and early insights obtained from the sequence data. The resources include a variant browser, a genotype imputation server, and genomic and phenotypic data that are available through dbGaP (Database of Genotypes and Phenotypes)1. In the first 53,831 TOPMed samples, we detected more than 400 million single-nucleotide and insertion or deletion variants after alignment with the reference genome. Additional previously undescribed variants were detected through assembly of unmapped reads and customized analysis in highly variable loci. Among the more than 400 million detected variants, 97% have frequencies of less than 1% and 46% are singletons that are present in only one individual (53% among unrelated individuals). These rare variants provide insights into mutational processes and recent human evolutionary history. The extensive catalogue of genetic variation in TOPMed studies provides unique opportunities for exploring the contributions of rare and noncoding sequence variants to phenotypic variation. Furthermore, combining TOPMed haplotypes with modern imputation methods improves the power and reach of genome-wide association studies to include variants down to a frequency of approximately 0.01%., The goals, resources and design of the NHLBI Trans-Omics for Precision Medicine (TOPMed) programme are described, and analyses of rare variants detected in the first 53,831 samples provide insights into mutational processes and recent human evolutionary history.

The frequencies of genetic variants in the CYP3A4 and CYP3A5 genes differ greatly across global populations, leading to profound differences in the metabolic activity of these enzymes and resulting drug metabolism rates, with important consequences for therapeutic safety and efficacy. Yet, the impact of genetic variants on enzyme activity are incompletely described, particularly in American Indian and Alaska Native (AIAN) populations. To characterize genetic variation in CYP3A4 and CYP3A5 and its effect on enzyme activity, we partnered with AIAN people living in two regions of Alaska: Yup’ik Alaska Native people living in the Yukon‐Kuskokwim Delta region of rural southwest Alaska and AIAN people receiving care at the Southcentral Foundation in Anchorage, Alaska. We identified low frequencies of novel and known variation in CYP3A4 and CYP3A5, including low frequencies of the CYP3A4*1G and CYP3A5*1 variants, and linkage disequilibrium patterns that differed from those we previously identified in an American Indian population in western Montana. We also identified increased activity of the CYP3A4*1G allele in vitro and in vivo. We demonstrated that the CYP3A4*1G allele confers increased protein content in human lymphoblastoid cells and both increased protein content and increased activity in human liver microsomes. We confirmed enhanced CYP3A4‐mediated 4β‐vitamin D hydroxylation activity in Yup’ik people with the CYP3A4*1G allele. AIAN people in Alaska and Montana who carry the CYP3A4*1G allele—coupled with low frequency of the functional CYP3A5*1 variant—may metabolize CYP3A substrates more rapidly than people with the reference CYP3A4 allele.

Background Prospective studies have suggested higher factor VIII (FVIII) levels are an independent risk factor for coronary heart disease (CHD) and stroke. However, limited information, including on genetic and epigenetic contributors to FVIII variation, is available specifically among African Americans (AAs), who have higher FVIII levels than Europeans. Objectives We measured FVIII levels in ~3400 AAs from the community-based Jackson Heart Study and assessed genetic, epigenetic, and epidemiological correlates of FVIII, as well as incident cardiovascular disease (CVD) associations. Methods We assessed cross-sectional associations of FVIII with CVD risk factors as well as incident CHD, stroke, heart failure, and mortality associations. We additionally assessed associations with TOPMed whole genome sequencing data and an epigenome-wide methylation array. Results Our results confirmed associations between FVIII and risk of incident CHD events and total mortality in AAs; mortality associations were largely independent of traditional risk factors. We also demonstrate an association of FVIII with incident heart failure, independent of B-type natriuretic peptide. Two genomic regions were strongly associated with FVIII (ABO and VWF). The index variant at VWF is specific to individuals of African descent and is distinct from the previously reported European VWF association signal. Epigenome-wide association analysis showed significant FVIII associations with several CpG sites in the ABO region. However, after adjusting for ABO genetic variants, ABO CpG sites were not significant. Conclusions Larger sample sizes of AAs will be required to discover additional genetic and epigenetic contributors to FVIII phenotypic variation, which may have consequences for CVD health disparities.

BACKGROUND: Asthma is a common respiratory disorder with a highly heterogeneous nature that remains poorly understood. The objective was to use whole genome sequencing (WGS) data to identify regions of common genetic variation contributing to lung function in individuals with a diagnosis of asthma. METHODS: WGS data were generated for 1,053 individuals from trios and extended pedigrees participating in the family-based Genetic Epidemiology of Asthma in Costa Rica study. Asthma affection status was defined through a physician’s diagnosis of asthma, and most participants with asthma also had airway hyperresponsiveness (AHR) to methacholine. Family-based association tests for single variants were performed to assess the associations with lung function phenotypes. RESULTS: A genome-wide significant association was identified between baseline FEV(1)/FVC ratio and a single-nucleotide polymorphism in the top hit cysteine-rich secretory protein LCCL domain-containing 2 (CRISPLD2) (rs12051168; P = 3.6 × 10(−8) in the unadjusted model) that retained suggestive significance in the covariate-adjusted model (P = 5.6 × 10(−6)). Rs12051168 was also nominally associated with other related phenotypes: baseline FEV(1) (P = 3.3 × 10(−3)), postbronchodilator (PB) FEV(1) (7.3 × 10(−3)), and PB FEV(1)/FVC ratio (P = 2.7 × 10(−3)). The identified baseline FEV(1)/FVC ratio and rs12051168 association was meta-analyzed and replicated in three independent cohorts in which most participants with asthma also had confirmed AHR (combined weighted z-score P = .015) but not in cohorts without information about AHR. CONCLUSIONS: These findings suggest that using specific asthma characteristics, such as AHR, can help identify more genetically homogeneous asthma subgroups with genotype-phenotype associations that may not be observed in all children with asthma. CRISPLD2 also may be important for baseline lung function in individuals with asthma who also may have AHR.

Background: suPAR (Soluble urokinase plasminogen activator receptor) has emerged as an important biomarker of coagulation, inflammation, and cardiovascular disease (CVD) risk. The contribution of suPAR to CVD risk and its genetic influence in Black populations have not been evaluated. Methods: We measured suPAR in 3492 Black adults from the prospective, community-based JHS (Jackson Heart Study). Cross-sectional associations of suPAR with lifestyle and CVD risk factors were assessed, whole-genome sequence data were used to evaluate genetic associations of suPAR, and relationships of suPAR with incident CVD outcomes and overall mortality were estimated over follow-up. Results: In Cox models adjusted for traditional CVD risk factors, estimated glomerular filtration rate, and CRP (C-reactive protein), each 1-SD higher suPAR was associated with a 21% to 31% increased risk of incident coronary heart disease, heart failure, stroke, and mortality. In the genome-wide association study, 2 missense (rs399145 encoding p.Thr86Ala, rs4760 encoding p.Phe272Leu) and 2 noncoding regulatory variants (rs73935023 within an enhancer element and rs4251805 within the promoter) of PLAUR on chromosome 19 were each independently associated with suPAR and together explained 14% of suPAR phenotypic variation. The allele frequencies of each of the four suPAR-associated genetic variants differ considerably across African and European populations. We further show that PLAUR rs73935023 can alter transcriptional activity in vitro. We did not find any association between genetically determined suPAR and CVD in JHS or a larger electronic medical record-based analyses of Blacks or Whites. Conclusions: Our results demonstrate the importance of ancestry-differentiated genetic variation on suPAR levels and indicate suPAR is a CVD biomarker in Black adults.

Linear mixed models are widely used in ecological and biological applications, especially in genetic studies. Reliable estimation of variance components is crucial for using linear mixed models. However, standard methods, such as the restricted maximum likelihood (REML), are computationally inefficient and may be unstable with small samples. Other commonly used methods, such as the Haseman-Elston (HE) regression, may yield negative estimates of variances. Utilizing regularized estimation strategies, we propose the restricted Haseman-Elston (REHE) regression and REHE with resampling (reREHE) estimators, along with an inference framework for REHE, as fast and robust alternatives that provide non-negative estimates with comparable accuracy to REML. The merits of REHE are illustrated using real data and benchmark simulation studies.

Background Genetic studies have primarily been conducted in European ancestry populations, identifying dozens of loci associated with late-onset Alzheimer’s disease (AD). However, much of AD’s heritability remains unexplained; as the prevalence of AD varies across populations, the genetic architecture of the disease may also vary by population with the presence of novel variants or loci. Methods We conducted genome-wide analyses of AD in a sample of 2565 Caribbean Hispanics to better understand the genetic contribution to AD in this population. Statistical analysis included both admixture mapping and association testing. Evidence for differential gene expression within regions of interest was collected from independent transcriptomic studies comparing AD cases and controls in samples with primarily European ancestry. Results Our genome-wide association study of AD identified no loci reaching genome-wide significance. However, a genome-wide admixture mapping analysis that tests for association between a haplotype’s ancestral origin and AD status detected a genome-wide significant association with chromosome 3q13.11 (103.7–107.7Mb, P = 8.76E−07), driven by a protective effect conferred by the Native American ancestry (OR = 0.58, 95%CI = 0.47−0.73). Within this region, two variants were significantly associated with AD after accounting for the number of independent tests (rs12494162, P = 2.33E−06; rs1731642, P = 6.36E−05). The significant admixture mapping signal is composed of 15 haplotype blocks spanning 5 protein-coding genes (ALCAM, BBX, CBLB, CCDC54, CD47) and four brain-derived topologically associated domains, and includes markers significantly associated with the expression of ALCAM, BBX, CBLB, and CD47 in the brain. ALCAM and BBX were also significantly differentially expressed in the brain between AD cases and controls with European ancestry. Conclusion These results provide multiethnic evidence for a relationship between AD and multiple genes at 3q13.11 and illustrate the utility of leveraging genetic ancestry diversity via admixture mapping for new insights into AD.

OBJECTIVE This work was undertaken in order to identify Parkinson's disease (PD) risk variants in a Latino cohort, to describe the overlap in the genetic architecture of PD in Latinos compared to European-ancestry subjects, and to increase the diversity in PD genome-wide association (GWAS) data. METHODS We genotyped and imputed 1,497 PD cases and controls recruited from nine clinical sites across South America. We performed a GWAS using logistic mixed models; variants with a p-value

The cytochrome P450 2D6 (CYP2D6) gene locus is challenging to accurately genotype due to numerous single nucleotide variants and complex structural variation. Our goal was to determine whether the CYP2D6 genotype-phenotype correlation is improved when diplotype assignments incorporate structural variation, identified by the bioinformatics tool Stargazer, with next-generation sequencing data. Using CYP2D6 activity measured with substrates dextromethorphan and metoprolol, activity score explained 40% and 34% of variability in metabolite formation rates, respectively, when diplotype calls incorporated structural variation, increasing from 36% and 31%, respectively, when diplotypes did not incorporate structural variation. We also investigated whether the revised Clinical Pharmacogenetics Implementation Consortium (CPIC) recommendations for translating genotype to phenotype improve CYP2D6 activity predictions over the current system. Although the revised recommendations do not improve the correlation between activity score and CYP2D6 activity, perhaps because of low frequency of the CYP2D6*10 allele, the correlation with metabolizer phenotype group was significantly improved for both substrates. We also measured the function of seven rare coding variants: one (A449D) exhibited decreased (44%) and another (R474Q) increased (127%) activity compared with reference CYP2D6.1 protein. Allele-specific analysis found that A449D is part of a novel CYP2D6*4 suballele, CYP2D6*4.028. The novel haplotype containing R474Q was designated CYP2D6*138 by PharmVar; another novel haplotype containing R365H was designated CYP2D6*139. Accuracy of CYP2D6 phenotype prediction is improved when the CYP2D6 gene locus is interrogated using next-generation sequencing coupled with structural variation analysis. Additionally, revised CPIC genotype to phenotype translation recommendations provides an improvement in assigning CYP2D6 activity.

Alaska Native and American Indian (AN/AI) people have unique pharmacogene variation that may affect warfarin disposition and therapeutic response. We performed targeted genotyping for cytochrome P450 (CYP)2C9, vitamin K epoxide oxidase reductase complex subunit 1 (VKORC1), CYP4F2, CYP4F11, and gamma‐glutamyl carboxylase (GGCX) variants in AN/AI people receiving warfarin. The primary outcome was stable warfarin dose, defined as one dose, and associated international normalized ratio within the target range, at least 6 months after starting therapy, with two matching doses at least 2 weeks apart. Genotype–phenotype relationships were assessed by multivariate regression analysis, adjusted for self‐reported heritage, age, gender, and concurrent statin use. VKORC1 genotype explained 34% of dose variability, with VKORC1 −1639G>A and 1173C>T associated with a 1.7 mg/day (P = 1.4e‐05) dose reduction. Additionally, CYP2C9 N218I was suggestively significant (P = 0.077), with heterozygotes requiring 1.1 mg/day less than reference individuals. Self‐reported heritage was significantly associated with dose, largely driven by differences in the diagnostic VKORC1 allele frequencies among AN/AI people.

Alaska Native and American Indian (AN/AI) people have unique pharmacogene variation that may affect warfarin disposition and therapeutic response. We performed targeted genotyping for cytochrome P450 (CYP)2C9, vitamin K epoxide oxidase reductase complex subunit 1 (VKORC1), CYP4F2, CYP4F11, and gamma-glutamyl carboxylase (GGCX) variants in AN/AI people receiving warfarin. The primary outcome was stable warfarin dose, defined as one dose, and associated international normalized ratio within the target range, at least 6 months after starting therapy, with two matching doses at least 2 weeks apart. Genotype–phenotype relationships were assessed by multivariate regression analysis, adjusted for self-reported heritage, age, gender, and concurrent statin use. VKORC1 genotype explained 34% of dose variability, with VKORC1 −1639G>A and 1173C>T associated with a 1.7 mg/day (P = 1.4e-05) dose reduction. Additionally, CYP2C9 N218I was suggestively significant (P = 0.077), with heterozygotes requiring 1.1 mg/day less than reference individuals. Self-reported heritage was significantly associated with dose, largely driven by differences in the diagnostic VKORC1 allele frequencies among AN/AI people.

Background Admixture mapping is a powerful approach for gene mapping of complex traits that leverages the diverse genetic ancestry in populations with recent admixture, such as Hispanic or Latino individuals in the United States. These individuals have an increased risk of chronic kidney disease (CKD). Methods We performed genome-wide admixture mapping for both CKD and estimated glomerular filtration rate (eGFR) in a sample of 12,601 participants from the Hispanic Community Health Study/Study of Latinos, with validation in a sample of 8191 African American participants from the Women's Health Initiative (WHI). We also compared the findings with those from a conventional genome-wide association study. Results Three novel ancestry-of-origin loci were identified on chromosomes 2, 14, and 15 for CKD and eGFR-loci harboring genes that have been implicated in kidney function. The chromosome 2 locus comprises two European ancestry regions encompassing the FSHR and NRXN1 genes, with European ancestry at this locus associated with increased CKD risk. The chromosome 14 locus, found within the DLK1-DIO3 imprinted domain, was associated with lower eGFR and driven by European ancestry. The eGFR-associated locus on chromosome 15 included intronic variants of RYR3 and was within an African-specific genomic region associated with higher eGFR. The genomewide association study failed to identify significant associations in these regions. We validated the chromosome 14 and 15 loci associated with eGFR in the WHI African American participants. Conclusions This study provides evidence of shared ancestry-specific genomic regions influencing eGFR in Hispanic or Latino individuals and African American individuals, and illustrates the potential for leveraging genetic ancestry in recently admixed populations for the discovery of novel candidate loci for kidney phenotypes.

BackgroundAdmixture mapping is a powerful approach for gene mapping of complex traits that leverages the diverse genetic ancestry in populations with recent admixture such as U.S. Hispanics/Latinos (HL), who have increased risk of chronic kidney disease (CKD).MethodsGenome-wide admixture mapping was performed for CKD and estimated glomerular filtration rate (eGFR) in a sample of 12,601 participants from the Hispanic Community Health Study/Study of Latinos, with validation in a sample of 8191 African Americans from the Women’s Health Initiative (WHI).ResultsThree novel ancestry-of-origin loci were identified on chromosomes 2, 14 and 15 for CKD and eGFR. The chromosome 2 locus (2p16.3) consisted of two European ancestry regions encompassing the FSHR and NRXN1 genes, with European ancestry at this locus associated with increased risk for CKD. The chromosome 14 locus (14q32.2) located within the DLK1-DIO3 imprinted domain was driven by European ancestry, and was associated with lower eGFR. The chromosome 15 locus (15q13.3-14) included intronic variants of RYR3 and was within an African-specific genomic region that was associated with higher eGFR. These findings were compared to the conventional genome-wide association study that failed to identify significant associations in these regions. We validated the chromosome 14 and 15 loci for eGFR in the WHI African Americans.ConclusionsThis study provides evidence of shared ancestry-specific genomic regions influencing eGFR in HL and African Americans, and illustrates the potential for leveraging genetic ancestry in recently admixed populations for novel discovery of kidney trait loci.

Objective: sCD14 (soluble CD14) is a circulating pattern recognition receptor involved in inflammatory signaling. Although it is known that sCD14 levels vary by race, information on the genetic and cardiovascular disease (CVD) risk relationships of sCD14 in Black participants is limited. Approach and Results: We measured sCD14 in plasma at the baseline exam from n=3492 Black participants from the JHS (Jackson Heart Study). We evaluated associations between sCD14 and subclinical CVD measures, incident CVD events, and mortality under 3 levels of covariate adjustment. We used whole-genome sequence data from the Trans-Omics for Personalized Medicine program to identify genetic associations with sCD14. Adjusting for CVD risk factors and C-reactive protein, higher sCD14 was significantly associated with increased risk of mortality (hazard ratio, 1.25 [95% CI, 1.17–1.32]), incident coronary heart disease (hazard ratio, 1.28 [95% CI, 1.11–1.47]), and incident heart failure (hazard ratio, 1.27 [95% CI, 1.15–1.41]), but not stroke (hazard ratio, 0.96 [95% CI, 0.84–1.09]). Some of these relationships differed by age or sex: the association between sCD14 and heart failure was only observed in females; there was an association between sCD14 and stroke only at younger ages (in the lowest tertile of age, CD14 region of chromosome 5q31. We additionally identified 2 novel statistically distinct genetic associations with sCD14 represented by index variants rs770147646 and rs57599368, also in the chromosome 5q31 region. Conclusions: sCD14 independently predicts CVD-related outcomes and mortality in Black participants.

Embracing Genetic Diversity to Improve Black Health Ideally, race will be replaced with genetic ancestry as a variable in medical research and practice. But until more ancestry data are available, ...

To date, over 90 Parkinson’s disease (PD) risk variants have been reported from genome-wide association studies (GWAS). However, these GWAS efforts have been limited to individuals of European and East Asian ancestry. We performed the first GWAS of Latino PD patients from South America, comparing 807 cases against 690 controls followed by association testing of suggestive loci in a replication cohort of 1,234 cases and 439,522 controls. We demonstrated that SNCA plays a significant role in PD etiology in a Latino cohort and identified a suggestive locus near NRROS on chromosome 3 that appeared to be driven by Peruvian subjects. We also characterized the overlap of PD genetic architecture between Europeans and Latinos with a replication of significant variants identified by Nalls et al. in their 2019 GWAS1, finding 80% concordance in direction of effect. We then leveraged the population history of Latinos via admixture mapping, identifying a significant locus on chromosome 14 in a joint test of ancestries, driven by the Native American ancestral background, and a significant locus on chromosome 6 in our test of African ancestry, containing the genes STXBP6 and RPS6KA2, respectively. Ultimately, our work reflects the most comprehensive characterization of PD genetic architecture in Latinos to date.

Alaska Native people are under-represented in genetic research but have unique gene variation that may critically impact their response to pharmacotherapy. Full resequencing of CYP2C9 in a cross-section of this population identified CYP2C9 Met1Leu (M1L), a novel, relatively common single nucleotide polymorphism hypothesized to confer CYP2C9 poor metabolizer phenotype by disrupting the start codon. M1L is present at a minor allele frequency of 6.3% in Yup’ik Alaska Native people and thus can contribute to the risk of an adverse drug response from narrow-therapeutic-index CYP2C9 substrates such as (S)-warfarin. This study’s objective was to characterize the catalytic efficiency of the Leu1 variant enzyme in vivo by evaluating the pharmacokinetic behavior of naproxen, a probe substrate for CYP2C9 activity, in genotyped Yup’ik participants. We first confirmed the selectivity of (S)-naproxen O-demethylation by CYP2C9 using activity-phenotyped human liver microsomes and selective cytochrome P450 inhibitors and then developed and validated a novel liquid chromatography mass spectrometry method for simultaneous quantification of (S)-naproxen, (S)-O-desmethylnaproxen, and naproxen acyl glucuronide in human urine. The average ratio of (S)-O-desmethylnaproxen to unchanged (S)-naproxen in urine was 18.0 ± 8.0 (n = 11) for the homozygous CYP2C9 Met1 reference group and 10.3 ± 6.6 (n = 11) for the Leu1 variant carrier group (P = 0.011). The effect of M1L variation on CYP2C9 function and its potential to alter the pharmacokinetics of drugs metabolized by the enzyme has clinical implications and should be included in a variant screening panel when pharmacogenetic testing in the Alaska Native population is warranted. Significance Statement The novel CYP2C9 Met1Leu variant in Alaska Native people was recently identified. This study validated (S)-naproxen as a CYP2C9 probe substrate to characterize the in vivo functional activity of the CYP2C9 Leu1 variant. The results of this pharmacogenetic-pharmacokinetic study suggest that the CYP2C9 Leu1 variant exhibits loss of enzyme activity. This finding may be important to consider when administering narrow-therapeutic-index medications metabolized by CYP2C9 and also compels further investigation to characterize novel genetic variation in understudied populations.

Maternal genetic variations, including variations in mitochondrial biogenesis (MB) and oxidative phosphorylation (OP), are associated with placental abruption (PA). However, the role of maternal-fetal genetic interactions (MFGI) and parent-of-origin (imprinting) effects in PA remain unknown. We investigated MFGI in MB-OP, and imprinting effects in relation to risk of PA. Among Peruvian mother-infant pairs (503 PA cases and 1,052 controls), independent single nucleotide polymorphisms (SNPs), with linkage-disequilibrium coefficient

Introduction Recent studies suggest that both sex-specific genetic risk factors and those shared between dementia and stroke are involved in dementia pathogenesis. Methods We performed both single-variant and gene-based genome-wide association studies of >11,000 whole genome sequences from the Women's Health Initiative cohort to discover loci associated with dementia, with adjustment for age, ethnicity, stroke, and venous thromboembolism status. Evidence for prior evidence of association and differential gene expression in dementia-related tissues and samples was gathered for each locus. Results Our multiethnic studies identified significant associations between variants within APOE, MYH11, FZD3, SORCS3, and GOLGA8B and risk of dementia. Ten genes implicated by these loci, including MYH11, FZD3, SORCS3, and GOLGA8B, were differentially expressed in the context of Alzheimer's disease. Discussion Our association of MYH11, FZD3, SORCS3, and GOLGA8B with dementia is supported by independent functional studies in human subjects, model systems, and associations with shared risk factors for stroke and dementia.

Background: Parkinson's disease is the second most common neurodegenerative disorder and affects people from all ethnic backgrounds, yet little is known about the genetics of Parkinson's disease in non-European populations. In addition, the overall identification of copy number variants at a genome-wide level has been understudied in Parkinson's disease patients. Objectives: To understand the genome-wide burden of copy number variants in Latinos and its association with Parkinson's disease. Methods: We used genome-wide genotyping data from 747 Parkinson's disease patients and 632 ancestry matched controls from the Latin American Research Consortium on the Genetics of Parkinson's disease. Results: Genome-wide copy number burden analysis showed no difference between patients vs. controls, whereas patients were significantly enriched for copy number variants overlapping known Parkinson's disease genes compared to controls (Odds Ratio: 3.97 [1.69 - 10.5], P = 0.018). PARK2 showed the strongest copy number burden, with 20 copy number variant carriers. These patients presented an earlier age of disease onset compared to patients with other copy number variants (median age at onset: 31 years vs. 57 years, P = 7.46 x 10-7). Conclusions: We found that Parkinson's disease patients are significantly enriched with copy number variants affecting known Parkinson's disease genes. We also identified that out of 250 patients with early-onset disease, 5.6% carried a copy number variant on PARK2 in our cohort. Our study is the first to analyze genome-wide copy number variants association in Latino Parkinson's disease patients and provides insights about this complex disease in this understudied population.

Biomarkers of chronic inflammation (such as C-reactive protein) have long been associated with cardiovascular disease and mortality; however, biomarkers involved in antiviral cytokine induction and adaptive immune system activation remain largely unexamined. We hypothesized the cytokine interferon gamma inducible protein 10 (IP-10) would be associated with clinical and subclinical cardiovascular disease and all-cause mortality in African Americans. We assessed these associations in the Jackson Heart Study (JHS) cohort and the REasons for Geographic and Racial Differences in Stroke (REGARDS) study. There was a modest association of IP-10 with higher odds of left ventricular hypertrophy (OR = 1.20 (95% confidence interval (CI) 1.03, 1.41) per standard deviation (SD) higher natural log-transformed IP-10 in JHS). We did not observe associations with ankle brachial index, intima-media thickness, or arterial calcification. Each SD higher increment of ln-transformed IP-10 concentration was associated with incident heart failure (hazard ratio (HR) 1.26; 95% CI 1.11, 1.42, p = 4x10-4) in JHS, and with overall mortality in both JHS (HR 1.12 per SD, 95% CI 1.03, 1.21, p = 7.5x10-3) and REGARDS (HR 1.31 per SD, 95% CI 1.10, 1.55, p = 2.0 x 10-3), adjusting for cardiovascular risk factors and C-reactive protein. However, we found no association between IP-10 and stroke or coronary heart disease. These results suggest a role of IP-10 in heart failure and mortality risk independent of C-reactive protein. Further research is needed to investigate how the body's response to chronic viral infection may mediate heart failure and overall mortality risk in African Americans.

Polygenic risk scores (PRSs) are weighted sums of risk allele counts of single-nucleotide polymorphisms (SNPs) associated with a disease or trait. PRSs are typically constructed based on published results from Genome-Wide Association Studies (GWASs), and the majority of which has been performed in large populations of European ancestry (EA) individuals. Although many genotype-trait associations have generalized across populations, the optimal choice of SNPs and weights for PRSs may differ between populations due to different linkage disequilibrium (LD) and allele frequency patterns. We compare various approaches for PRS construction, using GWAS results from both large EA studies and a smaller study in Hispanics/Latinos: The Hispanic Community Health Study/Study of Latinos (HCHS/SOL,mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"mml:min/mml:mimml:mo=/mml:momml:mn12/mml:mnmml:mo,/mml:momml:mn803/mml:mn/mml:math). We consider multiple approaches for selecting SNPs and for computing SNP weights. We study the performance of the resulting PRSs in an independent study of Hispanics/Latinos from the Women's Health Initiative (WHI,mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"mml:min/mml:mimml:mo=/mml:momml:mn3/mml:mnmml:mo,/mml:momml:mn582/mml:mn/mml:math). We support our investigation with simulation studies of potential genetic architectures in a single locus. We observed that selecting variants based on EA GWASs generally performs well, except for blood pressure trait. However, the use of EA GWASs for weight estimation was suboptimal. Using non-EA GWAS results to estimate weights improved results.

The Alzheimer’s Disease Sequencing Project (ADSP) undertook whole exome sequencing in 5,740 late-onset Alzheimer disease (AD) cases and 5,096 cognitively normal controls primarily of European ancestry (EA), among whom 218 cases and 177 controls were Caribbean Hispanic (CH). An age-, sex- and APOE based risk score and family history were used to select cases most likely to harbor novel AD risk variants and controls least likely to develop AD by age 85 years. We tested ~1.5 million single nucleotide variants (SNVs) and 50,000 insertion-deletion polymorphisms (indels) for association to AD, using multiple models considering individual variants as well as gene-based tests aggregating rare, predicted functional, and loss of function variants. Sixteen single variants and 19 genes that met criteria for significant or suggestive associations after multiple-testing correction were evaluated for replication in four independent samples; three with whole exome sequencing (2,778 cases, 7,262 controls) and one with genome-wide genotyping imputed to the Haplotype Reference Consortium panel (9,343 cases, 11,527 controls). The top findings in the discovery sample were also followed-up in the ADSP whole-genome sequenced family-based dataset (197 members of 42 EA families and 501 members of 157 CH families). We identified novel and predicted functional genetic variants in genes previously associated with AD. We also detected associations in three novel genes: IGHG3 (p = 9.8 × 10−7), an immunoglobulin gene whose antibodies interact with β-amyloid, a long non-coding RNA AC099552.4 (p = 1.2 × 10−7), and a zinc-finger protein ZNF655 (gene-based p = 5.0 × 10−6). The latter two suggest an important role for transcriptional regulation in AD pathogenesis.

Lipoprotein(a), Lp(a), is a modified low-density lipoprotein particle that contains apolipoprotein(a), encoded by LPA, and is a highly heritable, causal risk factor for cardiovascular diseases that varies in concentrations across ancestries. Here, we use deep-coverage whole genome sequencing in 8392 individuals of European and African ancestry to discover and interpret both single-nucleotide variants and copy number (CN) variation associated with Lp(a). We observe that genetic determinants between Europeans and Africans have several unique determinants. The common variant rs12740374 associated with Lp(a) cholesterol is an eQTL for SORT1 and independent of LDL cholesterol. Observed associations of aggregates of rare non-coding variants are largely explained by LPA structural variation, namely the LPA kringle IV 2 (KIV2)-CN. Finally, we find that LPA risk genotypes confer greater relative risk for incident atherosclerotic cardiovascular diseases compared to directly measured Lp(a), and are significantly associated with measures of subclinical atherosclerosis in African Americans., Circulating lipoprotein(a) is an important risk factor for cardiovascular disease and shows variability between different ethnic groups. Here, Zekavat et al. perform whole-genome sequencing in individuals of European and African ancestries and find ancestry-specific genetic determinants for lipoprotein(a) levels.

The major objective of this study was to investigate the association of genetic and nongenetic factors with variability in protein abundance and in vitro activity of the androgen-metabolizing enzyme UGT2B17 in human liver microsomes (n = 455). UGT2B17 abundance was quantified by liquid chromatography-tandem mass spectrometry proteomics, and enzyme activity was determined by using testosterone and dihydrotestosterone as in vitro probe substrates. Genotyping or gene resequencing and mRNA expression were also evaluated. Multivariate analysis was used to test the association of UGT2B17 copy number variation, single nucleotide polymorphisms (SNPs), age, and sex with its mRNA expression, abundance, and activity. UGT2B17 gene copy number and SNPs (rs7436962, rs9996186, rs28374627, and rs4860305) were associated with gene expression, protein levels, and androgen glucuronidation rates in a gene dose-dependent manner. UGT2B17 protein (mean ± S.D. picomoles per milligram of microsomal protein) is sparsely expressed in children younger than 9 years (0.12 ± 0.24 years) but profoundly increases from age 9 years to adults (∼10-fold) with ∼2.6-fold greater abundance in males than in females (1.2 vs. 0.47). Association of androgen glucuronidation with UGT2B15 abundance was observed only in the low UGT2B17 expressers. These data can be used to predict variability in the metabolism of UGT2B17 substrates. Drug companies should include UGT2B17 in early phenotyping assays during drug discovery to avoid late clinical failures.

Despite evidence that pharmacogenetics can improve tamoxifen pharmacotherapy, there are few studies with American Indian and Alaska Native (AIAN) people. We examined variation in cytochrome P450 (CYP) genes (CYP2D6, CYP3A4, CYP3A5, and CYP2C9) and tamoxifen biotransformation in AIAN patients with breast cancer (n = 42) from the Southcentral Foundation in Alaska and the Confederated Salish and Kootenai Tribes in Montana. We tested for associations between CYP diplotypes and plasma concentrations of tamoxifen and metabolites. Only the CYP2D6 variation was significantly associated with concentrations of endoxifen (P = 0.0008) and 4-hydroxytamoxifen (P = 0.0074), tamoxifen's principal active metabolites, as well as key metabolic ratios. The CYP2D6 was also the most significant predictor of active metabolites and metabolic ratios in a multivariate regression model, including all four genes as predictors, with minor roles for other CYP genes. In AIAN populations, CYP2D6 is the largest contributor to tamoxifen bioactivation, illustrating the importance of validating pharmacogenetic testing for therapy optimization in an understudied population.

Summary The Genomic Data Storage (GDS) format provides efficient storage and retrieval of genotypes measured by microarrays and sequencing. We developed GENESIS to perform various single- and aggregate-variant association tests using genotype data stored in GDS format. GENESIS implements highly flexible mixed models, allowing for different link functions, multiple variance components and phenotypic heteroskedasticity. GENESIS integrates cohesively with other R/Bioconductor packages to build a complete genomic analysis workflow entirely within the R environment. Availability and implementation https://bioconductor.org/packages/GENESIS; vignettes included. Supplementary information Supplementary data are available at Bioinformatics online.