Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed whole-genome sequencing study.

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions in lipid metabolism. Large-scale whole-genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess more associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with measurement of blood lipids and lipoproteins (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare-variant aggregate association tests using the STAAR (variant-set test for association using annotation information) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare-coding variants in nearby protein-coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500-kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variation and rare protein-coding variation at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNAs.

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.

Mendelian randomization in the multivariate general linear model framework.

Mendelian randomization (MR) is an application of instrumental variable (IV) methods to observational data in which the IV is a genetic variant. MR methods applicable to the general exponential family of distributions are currently not well characterized. We adapt a general linear model framework to the IV setting and propose a general MR method applicable to any full-rank distribution from the exponential family. Empirical bias and coverage are estimated via simulations. The proposed method is compared to several existing MR methods. Real data analyses are performed using data from the REGARDS study to estimate the potential causal effect of smoking frequency on stroke risk in African Americans. In simulations with binary variates and very weak instruments the proposed method had the lowest median [Q1 , Q3 ] bias (0.10 [-3.68 to 3.62]); compared with 2SPS (0.27 [-3.74 to 4.26]) and the Wald method (-0.69 [-1.72 to 0.35]). Low bias was observed throughout other simulation scenarios; as well as more than 90% coverage for the proposed method. In simulations with count variates, the proposed method performed comparably to 2SPS; the Wald method maintained the most consistent low bias; and 2SRI was biased towards the null. Real data analyses find no evidence for a causal effect of smoking frequency on stroke risk. The proposed MR method has low bias and acceptable coverage across a wide range of distributional scenarios and instrument strengths; and provides a more parsimonious framework for asymptotic hypothesis testing compared to existing two-stage procedures.

Alpha Globin Gene Copy Number Is Associated with Prevalent Chronic Kidney Disease and Incident End-Stage Kidney Disease among Black Americans.

BACKGROUND: α-Globin is expressed in endothelial cells of resistance arteries, where it limits endothelial nitric oxide signaling and enhances α-adrenergic-mediated vasoconstriction. α-Globin gene (HBA) copy number is variable in people of African descent and other populations worldwide. Given the protective effect of nitric oxide in the kidney, we hypothesized that HBA copy number would be associated with kidney disease risk. METHODS: Community-dwelling Black Americans aged ≥45 years old were enrolled in a national longitudinal cohort from 2003 through 2007. HBA copy number was measured using droplet digital PCR. The prevalence ratio (PR) of CKD and the relative risk (RR) of incident reduced eGFR were calculated using modified Poisson multivariable regression. The hazard ratio (HR) of incident ESKD was calculated using Cox proportional hazards multivariable regression. RESULTS: Among 9908 participants, HBA copy number varied from 2 to 6. In analyses adjusted for demographic, clinical, and genetic risk factors, a one-copy increase in HBA was associated with 14% greater prevalence of CKD (PR, 1.14; 95% CI, 1.07 to 1.21; P<0.0001). While HBA copy number was not associated with incident reduced eGFR (RR, 1.06; 95% CI, 0.94 to 1.19; P=0.38), the hazard of incident ESKD was 32% higher for each additional copy of HBA (HR, 1.32; 95% CI, 1.09 to 1.61; P=0.005). CONCLUSIONS: Increasing HBA copy number was associated with a greater prevalence of CKD and incidence of ESKD in a national longitudinal cohort of Black Americans.

Multi-ethnic GWAS and fine-mapping of glycaemic traits identify novel loci in the PAGE Study.

AIMS/HYPOTHESIS: Type 2 diabetes is a growing global public health challenge. Investigating quantitative traits, including fasting glucose, fasting insulin and HbA1c, that serve as early markers of type 2 diabetes progression may lead to a deeper understanding of the genetic aetiology of type 2 diabetes development. Previous genome-wide association studies (GWAS) have identified over 500 loci associated with type 2 diabetes, glycaemic traits and insulin-related traits. However, most of these findings were based only on populations of European ancestry. To address this research gap, we examined the genetic basis of fasting glucose, fasting insulin and HbA1c in participants of the diverse Population Architecture using Genomics and Epidemiology (PAGE) Study. METHODS: We conducted a GWAS of fasting glucose (n = 52,267), fasting insulin (n = 48,395) and HbA1c (n = 23,357) in participants without diabetes from the diverse PAGE Study (23% self-reported African American, 46% Hispanic/Latino, 40% European, 4% Asian, 3% Native Hawaiian, 0.8% Native American), performing transethnic and population-specific GWAS meta-analyses, followed by fine-mapping to identify and characterise novel loci and independent secondary signals in known loci. RESULTS: Four novel associations were identified (p < 5 × 10-9), including three loci associated with fasting insulin, and a novel, low-frequency African American-specific locus associated with fasting glucose. Additionally, seven secondary signals were identified, including novel independent secondary signals for fasting glucose at the known GCK locus and for fasting insulin at the known PPP1R3B locus in transethnic meta-analysis. CONCLUSIONS/INTERPRETATION: Our findings provide new insights into the genetic architecture of glycaemic traits and highlight the continued importance of conducting genetic studies in diverse populations. DATA AVAILABILITY: Full summary statistics from each of the population-specific and transethnic results are available at NHGRI-EBI GWAS catalog ( https://www.ebi.ac.uk/gwas/downloads/summary-statistics ).

Inflammation biomarkers and incident coronary heart disease: the Reasons for Geographic And Racial Differences in Stroke Study.

BACKGROUND: Individual inflammation biomarkers are associated with incident coronary heart disease (CHD) events. However, there is limited research on whether the risk for incident CHD is progressively higher with a higher number of inflammation biomarkers in abnormal levels. METHODS: We used data from 15,758 Reasons for Geographic and Racial Differences in Stroke (REGARDS) study participants aged ≥45 years without a history of CHD at baseline in 2003-2007. Abnormal levels of baseline high-sensitivity C-reactive protein, leukocyte count and serum albumin were defined as ≥3.8 mg/L (3rd tertile), ≥6.3 x 109 cells/L (3rd tertile), and <4.0 g/dL (1st tertile), respectively. The outcome was a composite of incident myocardial infarction or CHD death. RESULTS: Overall, 38.9% (n = 6,123) had 0, 36.6% (n = 5,774) had 1, 19.8% (n = 3,113) had 2 and 4.7% (n = 748) had 3 biomarkers of inflammation in abnormal levels. Over a median follow-up of 11.4 years, 954 (6.1%) participants had incident CHD. The rate of incident CHD per 1000 person-years for individuals with 0, 1, 2, and 3 biomarkers of inflammation in abnormal levels was 4.4 (95% confidence interval [CI]: 3.9-5.0), 6.3 (95% CI: 5.6-6.9), 8.8 (95% CI: 7.8-9.9), and 10.6 (95% CI: 8.1-13.1), respectively. Multi-variable adjusted hazard ratios for incident CHD associated with 1, 2 and 3 versus no inflammation biomarker in abnormal levels were 1.26 (95% CI: 1.07-1.49), 1.72 (95% CI: 1.43-2.07), and 1.84 (95% CI: 1.37-2.47), respectively (P-trend < .001). CONCLUSIONS: The number of inflammation markers in abnormal levels was associated with increased risk of incident CHD after multi-variable adjustment.

Increased chest CT derived bone and muscle measures capture markers of improved morbidity and mortality in COPD.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) is a disease of accelerated aging and is associated with comorbid conditions including osteoporosis and sarcopenia. These extrapulmonary conditions are highly prevalent yet frequently underdiagnosed and overlooked by pulmonologists in COPD treatment and management. There is evidence supporting a role for bone-muscle crosstalk which may compound osteoporosis and sarcopenia risk in COPD. Chest CT is commonly utilized in COPD management, and we evaluated its utility to identify low bone mineral density (BMD) and reduced pectoralis muscle area (PMA) as surrogates for osteoporosis and sarcopenia. We then tested whether BMD and PMA were associated with morbidity and mortality in COPD. METHODS: BMD and PMA were analyzed from chest CT scans of 8468 COPDGene participants with COPD and controls (smoking and non-smoking). Multivariable regression models tested the relationship of BMD and PMA with measures of function (6-min walk distance (6MWD), handgrip strength) and disease severity (percent emphysema and lung function). Multivariable Cox proportional hazards models were used to evaluate the relationship between sex-specific quartiles of BMD and/or PMA derived from non-smoking controls with all-cause mortality. RESULTS: COPD subjects had significantly lower BMD and PMA compared with controls. Higher BMD and PMA were associated with increased physical function and less disease severity. Participants with the highest BMD and PMA quartiles had a significantly reduced mortality risk (36% and 46%) compared to the lowest quartiles. CONCLUSIONS: These findings highlight the potential for CT-derived BMD and PMA to characterize osteoporosis and sarcopenia using equipment available in the pulmonary setting.

APOL1 Risk Variants Associated with Serum Albumin in a Population-Based Cohort Study.

INTRODUCTION: The association of apolipoprotein L1 (APOL1) nephropathy risk variants (APOL1), unique to African-ancestry (African-American [AA]) populations, with systemic inflammation, a contributor to chronic kidney disease (CKD) and end-stage kidney disease (ESKD) is ill-defined. This study aimed to describe the role of inflammatory markers in the relationship between APOL1 and incident kidney outcomes using a prospective cohort study. METHODS: APOL1 high-risk status under a recessive genetic model was studied in 10,605 AA adults aged ≥45 years from the Reasons for Geographic and Racial Differences in Stroke study. The primary variables of interest were inflammatory markers: C-reactive protein (mg/dL), white blood cell count (cells/mm3), and serum albumin (sALB) (mg/dL). High inflammation status was defined if at least one of these inflammatory markers exceeded clinical threshold. The association between APOL1 and biomarkers were assessed using regression models adjusting for age, sex, ancestry, hypertension, lipid medications, albumin-to-creatinine ratio, and estimated glomerular filtration rate (eGFR). Models were stratified by diabetes status. We identified incident ESKD using USRDS linkage, and we defined incident CKD as an eGFR <60 mL/min/1.73 m2 and ≥25% decline in the eGFR and normal baseline eGFR and tested for mediation of APOL1 and outcomes by biomarkers using the causal inference approach. RESULTS: Among 7,151 participants with data available on all inflammation markers, 4,479 participants had ≥1 marker meeting the clinical threshold. APOL1 high-risk status was associated with lower adjusted odds of reduced sALB {odds ratio (OR) (95% confidence interval [CI]): 0.59 [0.36, 0.96])}, and this association was significant in people with diabetes (OR [95% CI]: 0.40 [0.18, 0.89]) but not in those without diabetes. There was no association of APOL1 high-risk status with other markers or high inflammation status. APOL1 was independently associated with ESKD (OR [95% CI] = 1.78 [1.28, 2.48]) and CKD (OR [95% CI] = 1.38 [1.00, 1.91]). On mediation analysis, the direct effect between APOL1 and ESKD strengthened after accounting for sALB, but the estimated mediated effect was not statistically significant (OR [95% CI]: 0.98 [0.92, 1.05], p = 0.58). CONCLUSION: APOL1 high-risk variants were associated with sALB. However, sALB did not statistically mediate the association between APOL1 and incident ESKD.

Use of >100,000 NHLBI Trans-Omics for Precision Medicine (TOPMed) Consortium whole genome sequences improves imputation quality and detection of rare variant associations in admixed African and Hispanic/Latino populations.

Most genome-wide association and fine-mapping studies to date have been conducted in individuals of European descent, and genetic studies of populations of Hispanic/Latino and African ancestry are limited. In addition, these populations have more complex linkage disequilibrium structure. In order to better define the genetic architecture of these understudied populations, we leveraged >100,000 phased sequences available from deep-coverage whole genome sequencing through the multi-ethnic NHLBI Trans-Omics for Precision Medicine (TOPMed) program to impute genotypes into admixed African and Hispanic/Latino samples with genome-wide genotyping array data. We demonstrated that using TOPMed sequencing data as the imputation reference panel improves genotype imputation quality in these populations, which subsequently enhanced gene-mapping power for complex traits. For rare variants with minor allele frequency (MAF) < 0.5%, we observed a 2.3- to 6.1-fold increase in the number of well-imputed variants, with 11-34% improvement in average imputation quality, compared to the state-of-the-art 1000 Genomes Project Phase 3 and Haplotype Reference Consortium reference panels. Impressively, even for extremely rare variants with minor allele count <10 (including singletons) in the imputation target samples, average information content rescued was >86%. Subsequent association analyses of TOPMed reference panel-imputed genotype data with hematological traits (hemoglobin (HGB), hematocrit (HCT), and white blood cell count (WBC)) in ~21,600 African-ancestry and ~21,700 Hispanic/Latino individuals identified associations with two rare variants in the HBB gene (rs33930165 with higher WBC [p = 8.8x10-15] in African populations, rs11549407 with lower HGB [p = 1.5x10-12] and HCT [p = 8.8x10-10] in Hispanics/Latinos). By comparison, neither variant would have been genome-wide significant if either 1000 Genomes Project Phase 3 or Haplotype Reference Consortium reference panels had been used for imputation. Our findings highlight the utility of the TOPMed imputation reference panel for identification of novel rare variant associations not previously detected in similarly sized genome-wide studies of under-represented African and Hispanic/Latino populations.

Association of Transthyretin Val122Ile Variant With Incident Heart Failure Among Black Individuals.

Importance: A genetic variant in the TTR gene (rs76992529; Val122Ile), present more commonly in individuals with African ancestry (population frequency: 3%-4%), causes misfolding of the tetrameric transthyretin protein complex that accumulates as extracellular amyloid fibrils and results in hereditary transthyretin amyloidosis. Objective: To estimate the association of the amyloidogenic Val122Ile TTR variant with the risk of heart failure and mortality in a large, geographically diverse cohort of Black individuals. Design, Setting, and Participants: Retrospective population-based cohort study of 7514 self-identified Black individuals living in the US participating in the REGARDS (Reasons for Geographic and Racial Differences in Stroke) study with genetic data available and without heart failure at baseline. The participants were enrolled at the baseline visit (2003-2007). The end of follow-up for the majority of outcomes was on December 31, 2018. All-cause mortality data were available through December 31, 2020. Exposures: TTR Val122Ile (rs76992529) genotype. Main Outcome and Measures: The primary outcome was incident heart failure (first hospitalization for heart failure or death due to heart failure). The secondary outcomes were heart failure mortality, cardiovascular mortality, and all-cause mortality. The multivariable Cox proportional hazards regression analyses were adjusted for genetic ancestry and demographic, clinical, and social factors. Results: Among 7514 Black participants (median age, 64 years [IQR, 57-70 years]; 61% women), the population frequency of the TTR Val122Ile variant was 3.1% (232 variant carriers and 7282 noncarriers). During a median follow-up of 11.1 years (IQR, 5.9-13.5 years), incident heart failure occurred in 535 individuals (34 variant carriers and 501 noncarriers) and the incidence of heart failure was 15.64 per 1000 person-years among variant carriers vs 7.16 per 1000 person-years among noncarriers (adjusted hazard ratio [HR], 2.43 [95% CI, 1.71-3.46]; P < .001). Deaths due to heart failure occurred in 141 individuals (13 variant carriers and 128 noncarriers) and the incidence of heart failure mortality was 6.11 per 1000 person-years among variant carriers vs 1.85 per 1000 person-years among noncarriers (adjusted HR, 4.19 [95% CI, 2.33-7.54]; P < .001). Deaths due to cardiovascular causes occurred in 793 individuals (34 variant carriers and 759 noncarriers) and the incidence of cardiovascular death was 15.18 per 1000 person-years among variant carriers vs 10.61 per 1000 person-years among noncarriers (adjusted HR, 1.69 [95% CI, 1.19-2.39]; P = .003). Deaths due to any cause occurred in 2715 individuals (100 variant carriers and 2615 noncarriers) and the incidence of all-cause mortality was 41.46 per 1000 person-years among variant carriers vs 33.94 per 1000 person-years among noncarriers (adjusted HR, 1.46 [95% CI, 1.19-1.78]; P < .001). There was no significant interaction between TTR variant carrier status and sex on incident heart failure and the secondary outcomes. Conclusions and Relevance: Among a cohort of Black individuals living in the US, being a carrier of the TTR Val122Ile variant was significantly associated with an increased risk of heart failure.

Lipid Phenotypes and DNA Methylation: a Review of the Literature.

PURPOSE OF REVIEW: Epigenetic modifications via DNA methylation have previously been linked to blood lipid levels, dyslipidemias, and atherosclerosis. The purpose of this review is to discuss current literature on the role of DNA methylation on lipid traits and their associated pathologies. RECENT FINDINGS: Candidate gene and epigenome-wide approaches have identified differential methylation of genes associated with lipid traits (particularly CPT1A, ABCG1, SREBF1), and novel approaches are being implemented to further characterize these relationships. Moreover, studies on environmental factors have shown that methylation variations at lipid-related genes are associated with diet and pollution exposure. Further investigation is needed to elucidate the directionality of the associations between the environment, lipid traits, and epigenome. Future studies should also seek to increase the diversity of cohorts, as European and Asian ancestry populations are the predominant study populations in the current literature.

Age and sex are associated with the plasma lipidome: findings from the GOLDN study.

BACKGROUND: Developing an understanding of the biochemistry of aging in both sexes is critical for managing disease throughout the lifespan. Lipidomic associations with age and sex have been reported, but prior studies are limited by measurements in serum rather than plasma or by participants taking lipid-lowering medications. METHODS: Our study included lipidomic data from 980 participants aged 18-87 years old from the Genetics of Lipid-Lowering Drugs and Diet Network (GOLDN). Participants were off lipid-lowering medications for at least 4 weeks, and signal intensities of 413 known lipid species were measured in plasma. We examined linear age and sex associations with signal intensity of (a) 413 lipid species; (b) 6 lipid classes (glycerolipids, glycerophospholipids, sphingolipids, sterol lipids, fatty acids, and acylcarnitines); and (c) 15 lipid subclasses; as well as with the particle sizes of three lipoproteins. RESULTS: Significant age associations were identified in 4 classes, 11 subclasses, 147 species, and particle size of one lipoprotein while significant sex differences were identified in 5 classes, 12 subclasses, 248 species, and particle sizes of two lipoproteins. For many lipid species (n = 97), age-related associations were significantly different between males and females. Age*sex interaction effects were most prevalent among phosphatidylcholines, sphingomyelins, and triglycerides. CONCLUSION: We identified several lipid species, subclasses, and classes that differ by age and sex; these lipid phenotypes may serve as useful biomarkers for lipid changes and associated cardiovascular risk with aging in the future. Future studies of age-related changes throughout the adult lifespan of both sexes are warranted. TRIAL REGISTRATION: ClinicalTrials.gov NCT00083369 ; May 21, 2004.

DNA Methylation Analysis Identifies Loci for Blood Pressure Regulation.

Genome-wide association studies have identified hundreds of genetic variants associated with blood pressure (BP), but sequence variation accounts for a small fraction of the phenotypic variance. Epigenetic changes may alter the expression of genes involved in BP regulation and explain part of the missing heritability. We therefore conducted a two-stage meta-analysis of the cross-sectional associations of systolic and diastolic BP with blood-derived genome-wide DNA methylation measured on the Infinium HumanMethylation450 BeadChip in 17,010 individuals of European, African American, and Hispanic ancestry. Of 31 discovery-stage cytosine-phosphate-guanine (CpG) dinucleotides, 13 replicated after Bonferroni correction (discovery: N = 9,828, p < 1.0 × 10-7; replication: N = 7,182, p < 1.6 × 10-3). The replicated methylation sites are heritable (h2 > 30%) and independent of known BP genetic variants, explaining an additional 1.4% and 2.0% of the interindividual variation in systolic and diastolic BP, respectively. Bidirectional Mendelian randomization among up to 4,513 individuals of European ancestry from 4 cohorts suggested that methylation at cg08035323 (TAF1B-YWHAQ) influences BP, while BP influences methylation at cg00533891 (ZMIZ1), cg00574958 (CPT1A), and cg02711608 (SLC1A5). Gene expression analyses further identified six genes (TSPAN2, SLC7A11, UNC93B1, CPT1A, PTMS, and LPCAT3) with evidence of triangular associations between methylation, gene expression, and BP. Additional integrative Mendelian randomization analyses of gene expression and DNA methylation suggested that the expression of TSPAN2 is a putative mediator of association between DNA methylation at cg23999170 and BP. These findings suggest that heritable DNA methylation plays a role in regulating BP independently of previously known genetic variants.

Integrating hypertension phenotype and genotype with hybrid non-negative matrix factorization.

MOTIVATION: Hypertension is a heterogeneous syndrome in need of improved subtyping using phenotypic and genetic measurements with the goal of identifying subtypes of patients who share similar pathophysiologic mechanisms and may respond more uniformly to targeted treatments. Existing machine learning approaches often face challenges in integrating phenotype and genotype information and presenting to clinicians an interpretable model. We aim to provide informed patient stratification based on phenotype and genotype features. RESULTS: In this article, we present a hybrid non-negative matrix factorization (HNMF) method to integrate phenotype and genotype information for patient stratification. HNMF simultaneously approximates the phenotypic and genetic feature matrices using different appropriate loss functions, and generates patient subtypes, phenotypic groups and genetic groups. Unlike previous methods, HNMF approximates phenotypic matrix under Frobenius loss, and genetic matrix under Kullback-Leibler (KL) loss. We propose an alternating projected gradient method to solve the approximation problem. Simulation shows HNMF converges fast and accurately to the true factor matrices. On a real-world clinical dataset, we used the patient factor matrix as features and examined the association of these features with indices of cardiac mechanics. We compared HNMF with six different models using phenotype or genotype features alone, with or without NMF, or using joint NMF with only one type of loss We also compared HNMF with 3 recently published methods for integrative clustering analysis, including iClusterBayes, Bayesian joint analysis and JIVE. HNMF significantly outperforms all comparison models. HNMF also reveals intuitive phenotype-genotype interactions that characterize cardiac abnormalities. AVAILABILITY AND IMPLEMENTATION: Our code is publicly available on github at https://github.com/yuanluo/hnmf. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

APOL1 Nephropathy Risk Alleles and Mortality in African American Adults: A Cohort Study.

RATIONALE & OBJECTIVE: APOL1 nephropathy risk alleles are associated with the development of chronic kidney disease (CKD) in African Americans. Although CKD is an established risk factor for mortality, associations of APOL1 risk alleles with mortality are uncertain. STUDY DESIGN: Prospective cohort. SETTINGS & PARTICIPANTS: 10,380 African American and 17,485 white American participants in the Reasons for Geographic and Racial Differences in Stroke (REGARDS) Study. EXPOSURES: APOL1 nephropathy risk alleles. OUTCOMES: All-cause and cause-specific mortality. ANALYTICAL APPROACH: Cox proportional hazards models were used to examine the association of APOL1 high-risk genotypes (2 risk alleles) versus APOL1 low-risk genotypes (0/1 risk allele) with all-cause and cause-specific mortality in African Americans and examine the risk for all-cause mortality in African Americans with high-risk genotypes versus African Americans with low-risk genotypes and white Americans. RESULTS: APOL1 high-risk participants were younger and had a higher prevalence of albuminuria than low-risk participants. There was no statistically significant association of APOL1 high- versus low-risk genotypes with all-cause mortality in models adjusted for sociodemographic variables, comorbid conditions, and kidney function (HR, 0.88; 95% CI, 0.77-1.01). After further adjustment for genetic ancestry in a subset with available data, a statistically significant association emerged (HR, 0.81; 95% CI, 0.69-0.96). Associations differed by CKD status (Pinteraction=0.04), with African Americans with high-risk genotypes having lower risk for mortality than those with low-risk genotypes in fully adjusted models (HR, 0.78; 95% CI, 0.62-0.99) among those with CKD, but not those without CKD (HR, 0.84; 95% CI, 0.66-1.05). Compared with white Americans, African Americans with high-risk genotypes had a similar rate of mortality, whereas African Americans with low-risk genotypes had a higher rate of mortality (HR, 1.07; 95% CI, 1.00-1.14) in fully adjusted models. LIMITATIONS: Lack of follow-up measures of kidney function. CONCLUSIONS: African Americans with high-risk APOL1 genotypes had lower mortality than those with low-risk genotypes in multivariable-adjusted models including genetic ancestry.

Genome-Wide Association Meta-Analysis of Individuals of European Ancestry Identifies Suggestive Loci for Sodium Intake, Potassium Intake, and Their Ratio Measured from 24-Hour or Half-Day Urine Samples.

BACKGROUND: Excess sodium intake and insufficient potassium intake are risk factors for hypertension, but there is limited knowledge regarding genetic factors that influence intake. Twenty-hour or half-day urine samples provide robust estimates of sodium and potassium intake, outperforming other measures such as spot urine samples and dietary self-reporting. OBJECTIVE: The aim of this study was to investigate genomic regions associated with sodium intake, potassium intake, and sodium-to-potassium ratio measured from 24-h or half-day urine samples. METHODS: Using samples of European ancestry (mean age: 54.2 y; 52.3% women), we conducted a meta-analysis of genome-wide association studies in 4 cohorts with 24-h or half-day urine samples (n = 6,519), followed by gene-based analysis. Suggestive loci (P < 10-6) were examined in additional European (n = 844), African (n = 1,246), and Asian (n = 2,475) ancestry samples. RESULTS: We found suggestive loci (P < 10-6) for all 3 traits, including 7 for 24-h sodium excretion, 4 for 24-h potassium excretion, and 4 for sodium-to-potassium ratio. The most significant locus was rs77958157 near cocaine- and amphetamine-regulated transcript prepropeptide (CARTPT) , a gene involved in eating behavior and appetite regulation (P = 2.3 × 10-8 with sodium-to-potassium ratio). Two suggestive loci were replicated in additional samples: for sodium excretion, rs12094702 near zinc finger SWIM-type containing 5 (ZSWIM5) was replicated in the Asian ancestry sample reaching Bonferroni-corrected significance (P = 0.007), and for potassium excretion rs34473523 near sodium leak channel (NALCN) was associated at a nominal P value with potassium excretion both in European (P = 0.043) and African (P = 0.043) ancestry cohorts. Gene-based tests identified 1 significant gene for sodium excretion, CDC42 small effector 1 (CDC42SE1), which is associated with blood pressure regulation. CONCLUSIONS: We identified multiple suggestive loci for sodium and potassium intake near genes associated with eating behavior, nervous system development and function, and blood pressure regulation in individuals of European ancestry. Further research is needed to replicate these findings and to provide insight into the underlying genetic mechanisms by which these genomic regions influence sodium and potassium intake.

A lipidome-wide association study of the lipoprotein insulin resistance index.

BACKGROUND: The lipoprotein insulin resistance (LPIR) score was shown to predict insulin resistance (IR) and type 2 diabetes (T2D) in healthy adults. However, the molecular basis underlying the LPIR utility for classification remains unclear. OBJECTIVE: To identify small molecule lipids associated with variation in the LPIR score, a weighted index of lipoproteins measured by nuclear magnetic resonance, in the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN) study (n = 980). METHODS: Linear mixed effects models were used to test the association between the LPIR score and 413 lipid species and their principal component analysis-derived groups. Significant associations were tested for replication with homeostatic model assessment-IR (HOMA-IR), a phenotype correlated with the LPIR score (r = 0.48, p <  0.001), in the Heredity and Phenotype Intervention (HAPI) Heart Study (n = 590). RESULTS: In GOLDN, 319 lipids were associated with the LPIR score (false discovery rate-adjusted p-values ranging from 4.59 × 10- 161 to 49.50 × 10- 3). Factors 1 (triglycerides and diglycerides/storage lipids) and 3 (mixed lipids) were positively (ß = 0.025, p = 4.52 × 10- 71 and ß = 0.021, p = 5.84 × 10- 41, respectively) and factor 2 (phospholipids/non-storage lipids) was inversely (ß = - 0.013, p = 2.28 × 10- 18) associated with the LPIR score. These findings were replicated for HOMA-IR in the HAPI Heart Study (ß = 0.10, p = 1.21 × 10- 02 for storage, ß = - 0.13, p = 3.14 × 10- 04 for non-storage, and ß = 0.19, p = 8.40 × 10- 07 for mixed lipids). CONCLUSIONS: Non-storage lipidomics species show a significant inverse association with the LPIR metabolic dysfunction score and present a promising focus for future therapeutic and prevention studies.

Transitions to family caregiving: enrolling incident caregivers and matched non-caregiving controls from a population-based study.

BACKGROUND AND AIM: Providing care to an older adult with a disability has been associated with increased risk to the caregiver's health, but most previous studies of caregiving and health compare persons who are already caregivers with poorly matched non-caregiving controls and are often based on convenience samples. In this report, we describe the enrollment of persons who transitioned into a family caregiving role while participating in a national epidemiological study. METHODS: Participants in the REasons for Geographic And Racial Differences in Stroke (REGARDS) study were asked on two occasions 9-14 years apart if they were providing care on an ongoing basis to a family member with a chronic illness or disability. Those who answered "no" and "yes", respectively, to this caregiving question and reported sufficient caregiving responsibilities after their transitions were enrolled in the present study as incident caregivers (N = 251). Participants matched on multiple demographic and health history variables and who reported no history of caregiving were enrolled as non-caregiving controls (N = 251). RESULTS: Among eligible participants, 84% agreed to participate, and 47% of caregivers reported caring for a person with dementia. Descriptive analyses confirmed the success of the matching procedures for balancing the groups on multiple demographic and pre-caregiving health variables. Depressive symptoms and perceived stress increased significantly after the transition to caregiving. CONCLUSION: Comparable, population-based samples of incident caregivers and matched non-caregivers have been enrolled. Future analyses will examine within-person changes in health and circulating biomarkers as a function of the transition to caregiving.

Discovery and fine-mapping of adiposity loci using high density imputation of genome-wide association studies in individuals of African ancestry: African Ancestry Anthropometry Genetics Consortium.

Genome-wide association studies (GWAS) have identified >300 loci associated with measures of adiposity including body mass index (BMI) and waist-to-hip ratio (adjusted for BMI, WHRadjBMI), but few have been identified through screening of the African ancestry genomes. We performed large scale meta-analyses and replications in up to 52,895 individuals for BMI and up to 23,095 individuals for WHRadjBMI from the African Ancestry Anthropometry Genetics Consortium (AAAGC) using 1000 Genomes phase 1 imputed GWAS to improve coverage of both common and low frequency variants in the low linkage disequilibrium African ancestry genomes. In the sex-combined analyses, we identified one novel locus (TCF7L2/HABP2) for WHRadjBMI and eight previously established loci at P < 5×10-8: seven for BMI, and one for WHRadjBMI in African ancestry individuals. An additional novel locus (SPRYD7/DLEU2) was identified for WHRadjBMI when combined with European GWAS. In the sex-stratified analyses, we identified three novel loci for BMI (INTS10/LPL and MLC1 in men, IRX4/IRX2 in women) and four for WHRadjBMI (SSX2IP, CASC8, PDE3B and ZDHHC1/HSD11B2 in women) in individuals of African ancestry or both African and European ancestry. For four of the novel variants, the minor allele frequency was low (<5%). In the trans-ethnic fine mapping of 47 BMI loci and 27 WHRadjBMI loci that were locus-wide significant (P < 0.05 adjusted for effective number of variants per locus) from the African ancestry sex-combined and sex-stratified analyses, 26 BMI loci and 17 WHRadjBMI loci contained ≤ 20 variants in the credible sets that jointly account for 99% posterior probability of driving the associations. The lead variants in 13 of these loci had a high probability of being causal. As compared to our previous HapMap imputed GWAS for BMI and WHRadjBMI including up to 71,412 and 27,350 African ancestry individuals, respectively, our results suggest that 1000 Genomes imputation showed modest improvement in identifying GWAS loci including low frequency variants. Trans-ethnic meta-analyses further improved fine mapping of putative causal variants in loci shared between the African and European ancestry populations.

APOL1 Kidney Risk Variants and Cardiovascular Disease: An Individual Participant Data Meta-Analysis.

BACKGROUND: Two coding variants in the apo L1 gene (APOL1) are strongly associated with kidney disease in blacks. Kidney disease itself increases the risk of cardiovascular disease, but whether these variants have an independent direct effect on the risk of cardiovascular disease is unclear. Previous studies have had inconsistent results. METHODS: We conducted a two-stage individual participant data meta-analysis to assess the association of APOL1 kidney-risk variants with adjudicated cardiovascular disease events and death, independent of kidney measures. The analysis included 21,305 blacks from eight large cohorts. RESULTS: Over 8.9±5.0 years of follow-up, 2076 incident cardiovascular disease events occurred in the 16,216 participants who did not have cardiovascular disease at study enrollment. In fully-adjusted analyses, individuals possessing two APOL1 kidney-risk variants had similar risk of incident cardiovascular disease (coronary heart disease, myocardial infarction, stroke and heart failure; hazard ratio 1.11, 95% confidence interval, 0.96 to 1.28) compared to individuals with zero or one kidney-risk variant. The risk of coronary heart disease, myocardial infarction, stroke and heart failure considered individually was also comparable by APOL1 genotype. APOL1 genotype was also not associated with death. There was no difference in adjusted associations by level of kidney function, age, diabetes status, or body-mass index. CONCLUSIONS: In this large, two-stage individual participant data meta-analysis, APOL1 kidney-risk variants were not associated with incident cardiovascular disease or death independent of kidney measures.