Multi-ancestry polygenic risk scores for venous thromboembolism.

Venous thromboembolism (VTE) is a significant contributor to morbidity and mortality, with large disparities in incidence rates between Black and White Americans. Polygenic risk scores (PRSs) limited to variants discovered in genome-wide association studies in European-ancestry samples can identify European-ancestry individuals at high risk of VTE. However, there is limited evidence on whether high-dimensional PRS constructed using more sophisticated methods and more diverse training data can enhance the predictive ability and their utility across diverse populations. We developed PRSs for VTE using summary statistics from the International Network against Venous Thrombosis (INVENT) consortium genome-wide association studies meta-analyses of European- (71 771 cases and 1 059 740 controls) and African-ancestry samples (7482 cases and 129 975 controls). We used LDpred2 and PRS-CSx to construct ancestry-specific and multi-ancestry PRSs and evaluated their performance in an independent European- (6781 cases and 103 016 controls) and African-ancestry sample (1385 cases and 12 569 controls). Multi-ancestry PRSs with weights tuned in European-ancestry samples slightly outperformed ancestry-specific PRSs in European-ancestry test samples (e.g. the area under the receiver operating curve [AUC] was 0.609 for PRS-CSx_combinedEUR and 0.608 for PRS-CSxEUR [P = 0.00029]). Multi-ancestry PRSs with weights tuned in African-ancestry samples also outperformed ancestry-specific PRSs in African-ancestry test samples (PRS-CSxAFR: AUC = 0.58, PRS-CSx_combined AFR: AUC = 0.59), although this difference was not statistically significant (P = 0.34). The highest fifth percentile of the best-performing PRS was associated with 1.9-fold and 1.68-fold increased risk for VTE among European- and African-ancestry subjects, respectively, relative to those in the middle stratum. These findings suggest that the multi-ancestry PRS might be used to improve performance across diverse populations to identify individuals at highest risk for VTE.

Whole-exome sequencing study identifies four novel gene loci associated with diabetic kidney disease.

Diabetic kidney disease (DKD) is recognized as an important public health challenge. However, its genomic mechanisms are poorly understood. To identify rare variants for DKD, we conducted a whole-exome sequencing (WES) study leveraging large cohorts well-phenotyped for chronic kidney disease and diabetes. Our two-stage WES study included 4372 European and African ancestry participants from the Chronic Renal Insufficiency Cohort and Atherosclerosis Risk in Communities studies (stage 1) and 11 487 multi-ancestry Trans-Omics for Precision Medicine participants (stage 2). Generalized linear mixed models, which accounted for genetic relatedness and adjusted for age, sex and ancestry, were used to test associations between single variants and DKD. Gene-based aggregate rare variant analyses were conducted using an optimized sequence kernel association test implemented within our mixed model framework. We identified four novel exome-wide significant DKD-related loci through initiating diabetes. In single-variant analyses, participants carrying a rare, in-frame insertion in the DIS3L2 gene (rs141560952) exhibited a 193-fold increased odds [95% confidence interval (CI): 33.6, 1105] of DKD compared with noncarriers (P = 3.59 × 10-9). Likewise, each copy of a low-frequency KRT6B splice-site variant (rs425827) conferred a 5.31-fold higher odds (95% CI: 3.06, 9.21) of DKD (P = 2.72 × 10-9). Aggregate gene-based analyses further identified ERAP2 (P = 4.03 × 10-8) and NPEPPS (P = 1.51 × 10-7), which are both expressed in the kidney and implicated in renin-angiotensin-aldosterone system modulated immune response. In the largest WES study of DKD, we identified novel rare variant loci attaining exome-wide significance. These findings provide new insights into the molecular mechanisms underlying DKD.

Antithrombin, Protein C, and Protein S: Genome and Transcriptome-Wide Association Studies Identify 7 Novel Loci Regulating Plasma Levels.

BACKGROUND: Antithrombin, PC (protein C), and PS (protein S) are circulating natural anticoagulant proteins that regulate hemostasis and of which partial deficiencies are causes of venous thromboembolism. Previous genetic association studies involving antithrombin, PC, and PS were limited by modest sample sizes or by being restricted to candidate genes. In the setting of the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium, we meta-analyzed across ancestries the results from 10 genome-wide association studies of plasma levels of antithrombin, PC, PS free, and PS total. METHODS: Study participants were of European and African ancestries, and genotype data were imputed to TOPMed, a dense multiancestry reference panel. Each of the 10 studies conducted a genome-wide association studies for each phenotype and summary results were meta-analyzed, stratified by ancestry. Analysis of antithrombin included 25 243 European ancestry and 2688 African ancestry participants, PC analysis included 16 597 European ancestry and 2688 African ancestry participants, PSF and PST analysis included 4113 and 6409 European ancestry participants. We also conducted transcriptome-wide association analyses and multiphenotype analysis to discover additional associations. Novel genome-wide association studies and transcriptome-wide association analyses findings were validated by in vitro functional experiments. Mendelian randomization was performed to assess the causal relationship between these proteins and cardiovascular outcomes. RESULTS: Genome-wide association studies meta-analyses identified 4 newly associated loci: 3 with antithrombin levels (GCKR, BAZ1B, and HP-TXNL4B) and 1 with PS levels (ORM1-ORM2). transcriptome-wide association analyses identified 3 newly associated genes: 1 with antithrombin level (FCGRT), 1 with PC (GOLM2), and 1 with PS (MYL7). In addition, we replicated 7 independent loci reported in previous studies. Functional experiments provided evidence for the involvement of GCKR, SNX17, and HP genes in antithrombin regulation. CONCLUSIONS: The use of larger sample sizes, diverse populations, and a denser imputation reference panel allowed the detection of 7 novel genomic loci associated with plasma antithrombin, PC, and PS levels.

Associations of plasma sphingolipids with measures of insulin sensitivity, ß-cell function, and incident diabetes in Japanese Americans.

BACKGROUND AND AIMS: To prospectively investigate associations of plasma sphingolipids with insulin sensitivity, ß-cell function, and incident diabetes in the Japanese American Community Diabetes Study. METHODS AND RESULTS: Baseline plasma samples from adults without diabetes (n = 349; mean age 56.7 years, 51 % men) were assayed for circulating ceramide and sphingomyelin species. Adjusted regression models examined cross-sectional and longitudinal associations with insulin sensitivity (HOMA2-%S), ß-cell function (oral disposition index: DIo) and with incident diabetes over 5 years follow-up. Concentrations of four species (Ceramide C16:0, C18:0, C20:0, and C22:0) were inversely associated with HOMA2-%S at baseline (all P values < 0.05, Q values < 0.05) and change in HOMA2-%S over 5 years (all P values < 0.05, Q values < 0.05). No sphingolipids were associated with baseline or change in DIo. Of the four species associated with HOMA2-%S, only Ceramide C18:0 was significantly and positively associated with incident diabetes (RR/1SD 1.44, 95 % CI 1.10-1.80, P = 0.006, Q = 0.024). The association of plasma Ceramide C18:0 with the risk of diabetes was partially mediated by change in HOMA2-%S between baseline and 5 years (mediation proportion: 61.5 %, 95 % CI 21.1%-212.5 %). CONCLUSION: Plasma Ceramide C18:0 was associated with higher risk of incident diabetes which was partially mediated through a decrease in insulin sensitivity between baseline and five years. Circulating Ceramide C18:0 could be a potential biomarker for identifying those at risk of developing diabetes.

Comparison of adaptive multiple phenotype association tests using summary statistics in genome-wide association studies.

Genome-wide association studies have been successful mapping loci for individual phenotypes, but few studies have comprehensively interrogated evidence of shared genetic effects across multiple phenotypes simultaneously. Statistical methods have been proposed for analyzing multiple phenotypes using summary statistics, which enables studies of shared genetic effects while avoiding challenges associated with individual-level data sharing. Adaptive tests have been developed to maintain power against multiple alternative hypotheses because the most powerful single-alternative test depends on the underlying structure of the associations between the multiple phenotypes and a single nucleotide polymorphism (SNP). Here we compare the performance of six such adaptive tests: two adaptive sum of powered scores (aSPU) tests, the unified score association test (metaUSAT), the adaptive test in a mixed-models framework (mixAda) and two principal-component-based adaptive tests (PCAQ and PCO). Our simulations highlight practical challenges that arise when multivariate distributions of phenotypes do not satisfy assumptions of multivariate normality. Previous reports in this context focus on low minor allele count (MAC) and omit the aSPU test, which relies less than other methods on asymptotic and distributional assumptions. When these assumptions are not satisfied, particularly when MAC is low and/or phenotype covariance matrices are singular or nearly singular, aSPU better preserves type I error, sometimes at the cost of decreased power. We illustrate this trade-off with multiple phenotype analyses of six quantitative electrocardiogram traits in the Population Architecture using Genomics and Epidemiology (PAGE) study.

Association of brain microbleeds with risk factors, cognition, and MRI markers in MESA.

INTRODUCTION: Little is known about the epidemiology of brain microbleeds in racially/ethnically diverse populations. METHODS: In the Multi-Ethnic Study of Atherosclerosis, brain microbleeds were identified from 3T magnetic resonance imaging susceptibility-weighted imaging sequences using deep learning models followed by radiologist review. RESULTS: Among 1016 participants without prior stroke (25% Black, 15% Chinese, 19% Hispanic, 41% White, mean age 72), microbleed prevalence was 20% at age 60 to 64.9 and 45% at ≥85 years. Deep microbleeds were associated with older age, hypertension, higher body mass index, and atrial fibrillation, and lobar microbleeds with male sex and atrial fibrillation. Overall, microbleeds were associated with greater white matter hyperintensity volume and lower total white matter fractional anisotropy. DISCUSSION: Results suggest differing associations for lobar versus deep locations. Sensitive microbleed quantification will facilitate future longitudinal studies of their potential role as an early indicator of vascular pathology.

Epigenetic Age and the Risk of Incident Atrial Fibrillation.

BACKGROUND: The most prominent risk factor for atrial fibrillation (AF) is chronological age; however, underlying mechanisms are unexplained. Algorithms using epigenetic modifications to the human genome effectively predict chronological age. Chronological and epigenetic predicted ages may diverge in a phenomenon referred to as epigenetic age acceleration (EAA), which may reflect accelerated biological aging. We sought to evaluate for associations between epigenetic age measures and incident AF. METHODS: Measures for 4 epigenetic clocks (Horvath, Hannum, DNA methylation [DNAm] PhenoAge, and DNAm GrimAge) and an epigenetic predictor of PAI-1 (plasminogen activator inhibitor-1) levels (ie, DNAm PAI-1) were determined for study participants from 3 population-based cohort studies. Cox models evaluated for associations with incident AF and results were combined via random-effects meta-analyses. Two-sample summary-level Mendelian randomization analyses evaluated for associations between genetic instruments of the EAA measures and AF. RESULTS: Among 5600 participants (mean age, 65.5 years; female, 60.1%; Black, 50.7%), there were 905 incident AF cases during a mean follow-up of 12.9 years. Unadjusted analyses revealed all 4 epigenetic clocks and the DNAm PAI-1 predictor were associated with statistically significant higher hazards of incident AF, though the magnitudes of their point estimates were smaller relative to the associations observed for chronological age. The pooled EAA estimates for each epigenetic measure, with the exception of Horvath EAA, were associated with incident AF in models adjusted for chronological age, race, sex, and smoking variables. After multivariable adjustment for additional known AF risk factors that could also potentially function as mediators, pooled EAA measures for 2 clocks remained statistically significant. Five-year increases in EAA measures for DNAm GrimAge and DNAm PhenoAge were associated with 19% (adjusted hazard ratio [HR], 1.19 [95% CI, 1.09-1.31]; P<0.01) and 15% (adjusted HR, 1.15 [95% CI, 1.05-1.25]; P<0.01) higher hazards of incident AF, respectively. Mendelian randomization analyses for the 5 EAA measures did not reveal statistically significant associations with AF. CONCLUSIONS: Our study identified adjusted associations between EAA measures and incident AF, suggesting that biological aging plays an important role independent of chronological age, though a potential underlying causal relationship remains unclear. These aging processes may be modifiable and not constrained by the immutable factor of time.

Circulating ceramides and sphingomyelins and the risk of incident cardiovascular disease among people with diabetes: the strong heart study.

BACKGROUND: Plasma ceramides and sphingomyelins have been independently linked to diabetes risk, glucose and insulin levels, and the risk of several cardiovascular (CVD) outcomes. However, whether individual ceramide and sphingomyelin species contribute to CVD risk among people with type 2 diabetes is uncertain. Our goal was to evaluate associations of 4 ceramide and 4 sphingomyelin species with incident CVD in a longitudinal population-based study among American Indians with diabetes. METHODS: This analysis included participants with prevalent type 2 diabetes from two cohorts: a prospective cohort of 597 participants in the Strong Heart Family Study (116 incident CVD cases; mean age: 49 years; average length of follow-up: 14 years), and a nested case-control sample of 267 participants in the Strong Heart Study (78 cases of CVD and 189 controls; mean age: 61 years; average time until incident CVD in cases: 3.8 years). The average onset of diabetes was 7 years prior to sphingolipid measurement. Sphingolipid species were measured using liquid chromatography and mass spectrometry. Cox regression and logistic regression were used to assess associations of sphingolipid species with incident CVD; results were combined across cohorts using inverse-variance weighted meta-analysis. RESULTS: There were 194 cases of incident CVD in the two cohorts. In meta-analysis of the 2 cohort results, higher plasma levels of Cer-16 (ceramide with acylated palmitic acid) were associated with higher CVD risk (HR per two-fold higher Cer-16: 1.85; 95% CI 1.05-3.25), and higher plasma levels of sphingomyelin species with a very long chain saturated fatty acid were associated with lower CVD risk (HR per two-fold higher SM-22: 0.48; 95% CI 0.26-0.87), although none of the associations met our pre-specified threshold for statistical significance of p = 0.006. CONCLUSIONS: While replication of the findings from the SHS in other populations is warranted, our findings add to a growing body of research suggesting that ceramides, in particular Cer-16, not only are associated with higher diabetes risk, but may also be associated with higher CVD risk after diabetes onset. We also find support for the hypothesis that sphingomyelins with a very long chain saturated fatty acid are associated with lower CVD risk among adults with type 2 diabetes.

Nonclassical Monocytes (CD14dimCD16+) Are Associated With Carotid Intima-Media Thickness Progression for Men but Not Women: The Multi-Ethnic Study of Atherosclerosis-Brief Report.

[Figure: see text].

Proteomics and Population Biology in the Cardiovascular Health Study (CHS): design of a study with mentored access and active data sharing.

BACKGROUND: In the last decade, genomic studies have identified and replicated thousands of genetic associations with measures of health and disease and contributed to the understanding of the etiology of a variety of health conditions. Proteins are key biomarkers in clinical medicine and often drug-therapy targets. Like genomics, proteomics can advance our understanding of biology. METHODS AND RESULTS: In the setting of the Cardiovascular Health Study (CHS), a cohort study of older adults, an aptamer-based method that has high sensitivity for low-abundance proteins was used to assay 4979 proteins in frozen, stored plasma from 3188 participants (61% women, mean age 74 years). CHS provides active support, including central analysis, for seven phenotype-specific working groups (WGs). Each CHS WG is led by one or two senior investigators and includes 10 to 20 early or mid-career scientists. In this setting of mentored access, the proteomic data and analytic methods are widely shared with the WGs and investigators so that they may evaluate associations between baseline levels of circulating proteins and the incidence of a variety of health outcomes in prospective cohort analyses. We describe the design of CHS, the CHS Proteomics Study, characteristics of participants, quality control measures, and structural characteristics of the data provided to CHS WGs. We additionally highlight plans for validation and replication of novel proteomic associations. CONCLUSION: The CHS Proteomics Study offers an opportunity for collaborative data sharing to improve our understanding of the etiology of a variety of health conditions in older adults.

Plasma ceramides containing saturated fatty acids are associated with risk of type 2 diabetes.

Recent studies suggest that the type of saturated fatty acid bound to sphingolipids influences the biological activity of those sphingolipids. However, it is unknown whether associations of sphingolipids with diabetes may differ by the identity of bound lipid species. Here, we investigated associations of 15 ceramide (Cer) and SM species (i.e., all sphingolipids, measured with coefficient of variation less than 20%) with incident type 2 diabetes in the Cardiovascular Health Study (n = 3,645), a large cohort study of cardiovascular disease among elderly adults who were followed from 1989 to 2015. Diabetes incidence was defined as fasting glucose ≥126 mg/dl or nonfasting glucose ≥200 mg/dl; reported use of insulin or oral hypoglycemic medication; or documentation of diabetes diagnosis through the Centers for Medicare and Medicaid Services records. Associations of each sphingolipid with incident diabetes were assessed using a Cox proportional hazards regression model. We found that higher circulating levels of Cer with acylated palmitic acid (Cer-16), stearic acid containing Cer (Cer-18), arachidic acid containing Cer (Cer-20), and behenic acid containing Cer (Cer-22) were each associated with a higher risk of diabetes. The hazard ratios for incident diabetes per 1 SD higher log levels of each Cer species were as follows: 1.21 (95% CI: 1.09-1.34) for Cer-16, 1.23 (95% CI: 1.10-1.37) for Cer-18, 1.14 (95% CI: 1.02-1.26) for Cer-20, and 1.18 (95% CI: 1.06-1.32) for Cer-22. In conclusion, higher levels of Cer-16, Cer-18, Cer-20, and Cer-22 were associated with a higher risk of diabetes.

Multi-ethnic genome-wide association analyses of white blood cell and platelet traits in the Population Architecture using Genomics and Epidemiology (PAGE) study.

BACKGROUND: Circulating white blood cell and platelet traits are clinically linked to various disease outcomes and differ across individuals and ancestry groups. Genetic factors play an important role in determining these traits and many loci have been identified. However, most of these findings were identified in populations of European ancestry (EA), with African Americans (AA), Hispanics/Latinos (HL), and other races/ethnicities being severely underrepresented. RESULTS: We performed ancestry-combined and ancestry-specific genome-wide association studies (GWAS) for white blood cell and platelet traits in the ancestrally diverse Population Architecture using Genomics and Epidemiology (PAGE) Study, including 16,201 AA, 21,347 HL, and 27,236 EA participants. We identified six novel findings at suggestive significance (P < 5E-8), which need confirmation, and independent signals at six previously established regions at genome-wide significance (P < 2E-9). We confirmed multiple previously reported genome-wide significant variants in the single variant association analysis and multiple genes using PrediXcan. Evaluation of loci reported from a Euro-centric GWAS indicated attenuation of effect estimates in AA and HL compared to EA populations. CONCLUSIONS: Our results highlighted the potential to identify ancestry-specific and ancestry-agnostic variants in participants with diverse backgrounds and advocate for continued efforts in improving inclusion of racially/ethnically diverse populations in genetic association studies for complex traits.

A systematic analysis of protein-altering exonic variants in chronic obstructive pulmonary disease.

Genome-wide association studies (GWASs) have identified regions associated with chronic obstructive pulmonary disease (COPD). GWASs of other diseases have shown an approximately 10-fold overrepresentation of nonsynonymous variants, despite limited exonic coverage on genotyping arrays. We hypothesized that a large-scale analysis of coding variants could discover novel genetic associations with COPD, including rare variants with large effect sizes. We performed a meta-analysis of exome arrays from 218,399 controls and 33,851 moderate-to-severe COPD cases. All exome-wide significant associations were present in regions previously identified by GWAS. We did not identify any novel rare coding variants with large effect sizes. Within GWAS regions on chromosomes 5q, 6p, and 15q, four coding variants were conditionally significant (P < 0.00015) when adjusting for lead GWAS single-nucleotide polymorphisms A common gasdermin B (GSDMB) splice variant (rs11078928) previously associated with a decreased risk for asthma was nominally associated with a decreased risk for COPD [minor allele frequency (MAF) = 0.46, P = 1.8e-4]. Two stop variants in coiled-coil α-helical rod protein 1 (CCHCR1), a gene involved in regulating cell proliferation, were associated with COPD (both P < 0.0001). The SERPINA1 Z allele was associated with a random-effects odds ratio of 1.43 for COPD (95% confidence interval = 1.17-1.74), though with marked heterogeneity across studies. Overall, COPD-associated exonic variants were identified in genes involved in DNA methylation, cell-matrix interactions, cell proliferation, and cell death. In conclusion, we performed the largest exome array meta-analysis of COPD to date and identified potential functional coding variants. Future studies are needed to identify rarer variants and further define the role of coding variants in COPD pathogenesis.

Circulating Ceramides and Sphingomyelins and Risk of Mortality: The Cardiovascular Health Study.

BACKGROUND: Recent studies suggest that associations of ceramides (Cer) and sphingomyelins (SM) with health outcomes differ according to the fatty acid acylated to the sphingoid backbone. The purpose of this study was to assess associations of Cer and SM species with mortality. METHODS: The study population included participants from the Cardiovascular Health Study (CHS), a community-based cohort of adults aged ≥65 years who were followed from 1992-2015 (n = 4612). Associations of plasma Cer and SM species carrying long-chain (i.e., 16:0) and very-long-chain (i.e., 20:0, 22:0, 24:0) saturated fatty acids with mortality were assessed using Cox proportional hazards models. RESULTS: During a median follow-up of 10.2 years, 4099 deaths occurred. High concentrations of Cer and SM carrying fatty acid 16:0 were each associated with an increased risk of mortality. Conversely, high concentrations of several ceramide and sphingomyelin species carrying longer fatty acids were each associated with a decreased risk of mortality. The hazard ratios for total mortality per 2-fold difference in each Cer and SM species were: 1.89 (95% CI), 1.65-2.17 for Cer-16, 0.79 (95% CI, 0.70-0.88) for Cer-22, 0.74 (95% CI, 0.65-0.84) for Cer-24, 2.51 (95% CI, 2.01-3.14) for SM-16, 0.68 (95% CI, 0.58-0.79) for SM-20, 0.57 (95% CI, 0.49-0.67) for SM-22, and 0.66 (0.57-0.75) for SM-24. We found no association of Cer-20 with risk of death. CONCLUSIONS: Associations of Cer and SM with the risk of death differ according to the length of their acylated saturated fatty acid. Future studies are needed to explore mechanisms underlying these relationships.

Natural killer cells, gamma delta T cells and classical monocytes are associated with systolic blood pressure in the multi-ethnic study of atherosclerosis (MESA).

BACKGROUND: Hypertension is a major source of cardiovascular morbidity and mortality. Recent evidence from mouse models, genetic, and cross-sectional human studies suggest increased proportions of selected immune cell subsets may be associated with levels of systolic blood pressure (SBP). METHODS: We assayed immune cells from cryopreserved samples collected at the baseline examination (2000-2002) from 1195 participants from the multi-ethnic study of atherosclerosis (MESA). We used linear mixed models, with adjustment for age, sex, race/ethnicity, smoking, exercise, body mass index, education, diabetes, and cytomegalovirus titers, to estimate the associations between 30 immune cell subsets (4 of which were a priori hypotheses) and repeated measures of SBP (baseline and up to four follow-up measures) over 10 years. The analysis provides estimates of the association with blood pressure level. RESULTS: The mean age of the MESA participants at baseline was 64 ± 10 years and 53% were male. A one standard deviation (1-SD) increment in the proportion of γδ T cells was associated with 2.40 mmHg [95% confidence interval (CI) 1.34-3.42] higher average systolic blood pressure; and for natural killer cells, a 1-SD increment was associated with 1.88 mmHg (95% CI 0.82-2.94) higher average level of systolic blood pressure. A 1-SD increment in classical monocytes (CD14++CD16-) was associated with 2.01 mmHG (95% CI 0.79-3.24) lower average systolic blood pressure. There were no associations of CD4+ T helper cell subsets with average systolic blood pressure. CONCLUSION: These findings suggest that the innate immune system plays a role in levels of SBP whereas there were no associations with adaptive immune cells.

Ancestry-specific associations identified in genome-wide combined-phenotype study of red blood cell traits emphasize benefits of diversity in genomics.

BACKGROUND: Quantitative red blood cell (RBC) traits are highly polygenic clinically relevant traits, with approximately 500 reported GWAS loci. The majority of RBC trait GWAS have been performed in European- or East Asian-ancestry populations, despite evidence that rare or ancestry-specific variation contributes substantially to RBC trait heritability. Recently developed combined-phenotype methods which leverage genetic trait correlation to improve statistical power have not yet been applied to these traits. Here we leveraged correlation of seven quantitative RBC traits in performing a combined-phenotype analysis in a multi-ethnic study population. RESULTS: We used the adaptive sum of powered scores (aSPU) test to assess combined-phenotype associations between ~ 21 million SNPs and seven RBC traits in a multi-ethnic population (maximum n = 67,885 participants; 24% African American, 30% Hispanic/Latino, and 43% European American; 76% female). Thirty-nine loci in our multi-ethnic population contained at least one significant association signal (p < 5E-9), with lead SNPs at nine loci significantly associated with three or more RBC traits. A majority of the lead SNPs were common (MAF > 5%) across all ancestral populations. Nineteen additional independent association signals were identified at seven known loci (HFE, KIT, HBS1L/MYB, CITED2/FILNC1, ABO, HBA1/2, and PLIN4/5). For example, the HBA1/2 locus contained 14 conditionally independent association signals, 11 of which were previously unreported and are specific to African and Amerindian ancestries. One variant in this region was common in all ancestries, but exhibited a narrower LD block in African Americans than European Americans or Hispanics/Latinos. GTEx eQTL analysis of all independent lead SNPs yielded 31 significant associations in relevant tissues, over half of which were not at the gene immediately proximal to the lead SNP. CONCLUSION: This work identified seven loci containing multiple independent association signals for RBC traits using a combined-phenotype approach, which may improve discovery in genetically correlated traits. Highly complex genetic architecture at the HBA1/2 locus was only revealed by the inclusion of African Americans and Hispanics/Latinos, underscoring the continued importance of expanding large GWAS to include ancestrally diverse populations.

Incorporating sampling weights into robust estimation of Cox proportional hazards regression model, with illustration in the Multi-Ethnic Study of Atherosclerosis.

BACKGROUND: Cox proportional hazards regression models are used to evaluate associations between exposures of interest and time-to-event outcomes in observational data. When exposures are measured on only a sample of participants, as they are in a case-cohort design, the sampling weights must be incorporated into the regression model to obtain unbiased estimating equations. METHODS: Robust Cox methods have been developed to better estimate associations when there are influential outliers in the exposure of interest, but these robust methods do not incorporate sampling weights. In this paper, we extend these robust methods, which already incorporate influence weights, so that they also accommodate sampling weights. RESULTS: Simulations illustrate that in the presence of influential outliers, the association estimate from the weighted robust method is closer to the true value than the estimate from traditional weighted Cox regression. As expected, in the absence of outliers, the use of robust methods yields a small loss of efficiency. Using data from a case-cohort study that is nested within the Multi-Ethnic Study of Atherosclerosis (MESA) longitudinal cohort study, we illustrate differences between traditional and robust weighted Cox association estimates for the relationships between immune cell traits and risk of stroke. CONCLUSIONS: Robust weighted Cox regression methods are a new tool to analyze time-to-event data with sampling, e.g. case-cohort data, when exposures of interest contain outliers.

Genome-wide meta-analysis of SNP-by9-ACEI/ARB and SNP-by-thiazide diuretic and effect on serum potassium in cohorts of European and African ancestry.

We evaluated interactions of SNP-by-ACE-I/ARB and SNP-by-TD on serum potassium (K+) among users of antihypertensive treatments (anti-HTN). Our study included seven European-ancestry (EA) (N = 4835) and four African-ancestry (AA) cohorts (N = 2016). We performed race-stratified, fixed-effect, inverse-variance-weighted meta-analyses of 2.5 million SNP-by-drug interaction estimates; race-combined meta-analysis; and trans-ethnic fine-mapping. Among EAs, we identified 11 significant SNPs (P < 5 × 10-8) for SNP-ACE-I/ARB interactions on serum K+ that were located between NR2F1-AS1 and ARRDC3-AS1 on chromosome 5 (top SNP rs6878413 P = 1.7 × 10-8; ratio of serum K+ in ACE-I/ARB exposed compared to unexposed is 1.0476, 1.0280, 1.0088 for the TT, AT, and AA genotypes, respectively). Trans-ethnic fine mapping identified the same group of SNPs on chromosome 5 as genome-wide significant for the ACE-I/ARB analysis. In conclusion, SNP-by-ACE-I /ARB interaction analyses uncovered loci that, if replicated, could have future implications for the prevention of arrhythmias due to anti-HTN treatment-related hyperkalemia. Before these loci can be identified as clinically relevant, future validation studies of equal or greater size in comparison to our discovery effort are needed.