Performing highly parallelized and reproducible GWAS analysis on biobank-scale data

Genome-wide association studies (GWAS) are transforming genetic research and enable the detection of novel genotype-phenotype relationships. In the last two decades, over 60 000 genetic associations across thousands of traits have been discovered using a GWAS approach. Due to increasing sample sizes, researchers are increasingly faced with computational challenges. A reproducible, modular and extensible pipeline with a focus on parallelization is essential to simplify data analysis and to allow researchers to devote their time to other essential tasks. Here we present nf-gwas, a Nextflow pipeline to run biobank-scale GWAS analysis. The pipeline automatically performs numerous pre- and post-processing steps, integrates regression modeling from the REGENIE package and supports single-variant, gene-based and interaction testing. It includes an extensive reporting functionality that allows to inspect thousands of phenotypes and navigate interactive Manhattan plots directly in the web browser. The pipeline is tested using the unit-style testing framework nf-test, a crucial requirement in clinical and pharmaceutical settings. Furthermore, we validated the pipeline against published GWAS datasets and benchmarked the pipeline on high-performance computing and cloud infrastructures to provide cost estimations to end users. nf-gwas is a highly parallelized, scalable and well-tested Nextflow pipeline to perform GWAS analysis in a reproducible manner.

Penetrance, cancer incidence and survival in HFE haemochromatosis-A population-based cohort study

BACKGROUND AND AIMS

Haemochromatosis is characterized by progressive iron overload affecting the liver and can cause cirrhosis and hepatocellular carcinoma. Most haemochromatosis patients are homozygous for p.C282Y in HFE, but only a minority of individuals with this genotype will develop the disease. The aim was to assess the penetrance of iron overload, fibrosis, hepatocellular carcinoma and life expectancy.

METHODS

A total of 8839 individuals from the Austrian region of Tyrol were genotyped for the p.C282Y variant between 1997 and 2021. Demographic, laboratory parameters and causes of death were assessed from health records. Penetrance, survival, and cancer incidence were ascertained from diagnosed cases, insurance- and cancer registry data. Outcomes were compared with a propensity score-matched control population.

RESULTS

Median age at diagnosis in 542 p.C282Y homozygous individuals was 47.8 years (64% male). At genotyping, the prevalence of iron overload was 55%. The cumulative penetrance of haemochromatosis defined as the presence of provisional iron overload was 24.2% in males and 10.5% in females aged 60 years or younger. Among p.C282Y homozygotes of the same ages, the cumulative proportion of individuals without fibrosis (FIB-4 score < 1.3) was 92.8% in males and 96.7% in females. Median life expectancy was reduced by 6.8 years in individuals homozygous for p.C282Y when compared with population-matched controls (p = .001). Hepatocellular carcinoma incidence was not significantly higher in p.C282Y homozygotes than in controls matched for age and sex.

CONCLUSION

Reduced survival and the observed age-dependent increase in penetrance among p.C282Y homozygotes call for earlier diagnosis of haemochromatosis to prevent complications.

Nuclear and mitochondrial genetic variants associated with mitochondrial DNA copy number

Mitochondrial DNA copy number (mtDNA-CN) is a biomarker for mitochondrial dysfunction associated with several diseases. Previous genome-wide association studies (GWAS) have been performed to unravel underlying mechanisms of mtDNA-CN regulation. However, the identified gene regions explain only a small fraction of mtDNA-CN variability. Most of this data has been estimated from microarrays based on various pipelines. In the present study we aimed to (1) identify genetic loci for qPCR-measured mtDNA-CN from three studies (16,130 participants) using GWAS, (2) identify potential systematic differences between our qPCR derived mtDNA-CN measurements compared to the published microarray intensity-based estimates, and (3) disentangle the nuclear from mitochondrial regulation of the mtDNA-CN phenotype. We identified two genome-wide significant autosomal loci associated with qPCR-measured mtDNA-CN: at HBS1L (rs4895440, p = 3.39 × 10-13) and GSDMA (rs56030650, p = 4.85 × 10-08) genes. Moreover, 113/115 of the previously published SNPs identified by microarray-based analyses were significantly equivalent with our findings. In our study, the mitochondrial genome itself contributed only marginally to mtDNA-CN regulation as we only detected a single rare mitochondrial variant associated with mtDNA-CN. Furthermore, we incorporated mitochondrial haplogroups into our analyses to explore their potential impact on mtDNA-CN. However, our findings indicate that they do not exert any significant influence on our results.

Use of Analog and Human Insulin in a European Hemodialysis Cohort With Type 2 Diabetes: Associations With Mortality, Hospitalization, MACE, and Hypoglycemia

RATIONALE & OBJECTIVE

Poor glycemic control may contribute to the high mortality rate in patients with type 2 diabetes receiving hemodialysis. Insulin type may influence glycemic control, and its choice may be an opportunity to improve outcomes. This study assessed whether treatment with analog insulin compared with human insulin is associated with different outcomes in people with type 2 diabetes and kidney failure receiving hemodialysis.

STUDY DESIGN

Retrospective cohort study.

SETTING & PARTICIPANTS

People in the Analyzing Data, Recognizing Excellence and Optimizing Outcomes (AROii) study with kidney failure commencing hemodialysis and type 2 diabetes being treated with insulin within 288 dialysis facilities between 2007 and 2009 across 7 European countries. Study participants were followed for 3 years. People with type 1 diabetes were excluded using an established administrative data algorithm.

EXPOSURE

Treatment with an insulin analog or human insulin.

OUTCOME

All-cause mortality, major adverse cardiovascular events (MACE), all-cause hospitalization, and confirmed hypoglycemia (blood glucose<3.0mmol/L sampled during hemodialysis).

ANALYTICAL APPROACH

Inverse probability weighted Cox proportional hazards models to estimate hazard ratios for analog insulin compared with human insulin.

RESULTS

There were 713 insulin analog and 733 human insulin users. Significant variation in insulin type by country was observed. Comparing analog with human insulin at 3 years, the percentage of patients experiencing end points and adjusted hazard ratios (AHR) were 22.0% versus 31.4% (AHR, 0.808 [95% CI, 0.66-0.99], P=0.04) for all-cause mortality, 26.8% versus 35.9% (AHR, 0.817 [95% CI, 0.68-0.98], P=0.03) for MACE, and 58.2% versus 75.0% (AHR, 0.757 [95% CI, 0.67-0.86], P<0.001) for hospitalization. Hypoglycemia was comparable between insulin types at 14.1% versus 15.0% (AHR, 1.169 [95% CI, 0.80-1.72], P=0.4). Consistent strength and direction of the associations were observed across sensitivity analyses.

LIMITATIONS

Residual confounding, lack of more detailed glycemia data.

CONCLUSIONS

In this large multinational cohort of people with type 2 diabetes and kidney failure receiving maintenance hemodialysis, treatment with analog insulins was associated with better clinical outcomes when compared with human insulin.

PLAIN-LANGUAGE SUMMARY

People with diabetes who are receiving dialysis for kidney failure are at high risk of cardiovascular disease and death. This study uses information from 1,446 people with kidney failure from 7 European countries who are receiving dialysis, have type 2 diabetes, and are prescribed either insulin identical to that made in the body (human insulin) or insulins with engineered extra features (insulin analog). After 3 years, fewer participants receiving analog insulins had died, had been admitted to the hospital, or had a cardiovascular event (heart attack, stroke, heart failure, or peripheral vascular disease). These findings suggest that analog insulins should be further explored as a treatment leading to better outcomes for people with diabetes on dialysis.

Meta-GWAS on PCSK9 concentrations reveals associations of novel loci outside the PCSK9 locus in White populations

BACKGROUND AND AIMS

Proprotein convertase subtilisin/kexin type 9 (PCSK9) is a key regulator of lipid homeostasis. A few earlier genome-wide association studies (GWAS) investigated genetic variants associated with circulating PCSK9 concentrations. However, uncertainty remains about some of the genetic loci discovered beyond the PCSK9 locus. By conducting the largest PCSK9 meta-analysis of GWAS (meta-GWAS) so far, we aimed to identify novel loci and validate the previously reported loci that regulate PCSK9 concentrations.

METHODS

We performed GWAS for PCSK9 concentrations in two large cohorts (GCKD (n = 4,963) and KORA F3 (n = 2,895)). These were meta-analyzed with previously published data encompassing together 20,579 individuals. We further conducted a second meta-analysis in statin-naïve individuals (n = 15,390). A genetic risk score (GRS) was constructed on PCSK9-increasing SNPs and assessed its impact on the risk for coronary artery disease (CAD) in 394,943 statin-naïve participants (17,077 with events) of the UK Biobank by performing CAD-free survival analysis.

RESULTS

Nine loci were genome-wide significantly associated with PCSK9 concentrations. These included the previously described PCSK9, APOB, KCNA1/KCNA5, and TM6SF2/SUGP1 loci. All imputed SNPs in the PCSK9 locus account for ∼15% of variance of PCSK9 concentrations. We further identified FADS2 as a novel locus that was also found in statin-naïve participants. All imputed SNPs within the FADS2 locus explain ∼1.2% of variance of PCSK9 concentrations. Additionally, four further loci (a region on chromosome 5, SDK1, SPATA16 and HPR) were genome-wide significant in either the main model or the statin-naïve subset. The linear increase in a PCSK9 genetic risk score was associated with 1.41-fold (95%CI 1.16-1.72, p < 0.001) higher risk for incident CAD.

CONCLUSIONS

We identified five novel loci (FADS2, SPATA16, SDK1, HPR and a region on chromosome 5) for PCSK9 concentrations that would require further research. Additionally, we confirm the genome-wide significant loci that were previously detected.

Systemic Evidence for Mitochondrial Dysfunction in Age-Related Macular Degeneration as Revealed by mtDNA Copy Number Measurements in Peripheral Blood

Mitochondrial dysfunction is a common occurrence in the aging process and is observed in diseases such as age-related macular degeneration (AMD). Increased levels of reactive oxygen species lead to damaged mitochondrial DNA (mtDNA), resulting in dysfunctional mitochondria, and, consequently, mtDNA causes further harm in the retinal tissue. However, it is unclear whether the effects are locally restricted to the high-energy-demanding retinal pigment epithelium or are also systematically present. Therefore, we measured mtDNA copy number (mtDNA-CN) in peripheral blood using a qPCR approach with plasmid normalization in elderly participants with and without AMD from the AugUR study (n = 2262). We found significantly lower mtDNA-CN in the blood of participants with early (n = 453) and late (n = 170) AMD compared to AMD-free participants (n = 1630). In regression analyses, we found lower mtDNA-CN to be associated with late AMD when compared with AMD-free participants. Each reduction of mtDNA-CN by one standard deviation increased the risk for late AMD by 24%. This association was most pronounced in geographic atrophy (OR = 1.76, 95% CI 1.19-2.60, p = 0.004), which has limited treatment options. These findings provide new insights into the relationship between mtDNA-CN in blood and AMD, suggesting that it may serve as a more accessible biomarker than mtDNA-CN in the retina.

A polygenic and family risk score are both independently associated with risk of type 2 diabetes in a population-based study

The availability of polygenic scores for type 2 diabetes (T2D) raises the question, whether assessing family history might become redundant. However, family history not only involves shared genetics, but also shared environment. It was the aim of this study to assess the independent and combined effects of one family risk score (FamRS) and a polygenic score (PGS) on prevalent and incident T2D risk in a population-based study from Germany (n = 3071). The study was conducted in 2004/2005 with up to 12 years of follow-up. The FamRS takes into account not only the number of diseased first grade relatives, but also age at onset of the relatives and age of participants. 256 prevalent and additional 163 incident T2D cases were recorded. Prevalence of T2D increased sharply for those within the top quantile of the PGS distribution resulting in an OR of 19.16 (p < 2 × 10^-16) for the top 20% compared to the remainder of the population, independent of age, sex, BMI, physical activity and FamRS. On the other hand, having a very strong family risk compared to average was still associated with an OR of 2.78 (p = 0.001), independent of the aforementioned factors and the PGS. The PGS and FamRS were only slightly correlated (r²_Spearman = 0.018). The combined contribution of both factors varied with varying age-groups, though, with decreasing influence of the PGS with increasing age. To conclude, both, genetic information and family history are relevant for the prediction of T2D risk and might be used for identification of high risk groups to personalize prevention measures.

Bariatric surgery prevents carotid wall thickness progression

BACKGROUND

Bariatric surgery is a treatment option for patients with severe obesity and improves parameters of cardiovascular and/or metabolic disease. Carotid intima media thickness (C-IMT) is a surrogate measure of subclinical atherosclerosis. Previous studies showed short to mid-term arrest and even regression of C‑IMT progression following bariatric surgery. We aimed to investigate the long-term effect of weight loss on C‑IMT progression 10 years after bariatric surgery in comparison to a population-based control cohort.

METHODS

In total, 21 eligible patients were examined preoperatively, at 5 and 10 years after bariatric surgery. Anthropometric parameters, plasma triglycerides, total cholesterol, high-density lipoprotein cholesterol (HDL-C), insulin, and glucose were assessed at all three study visits. C‑IMT was measured via B‑mode scans of the common carotid artery. C‑IMT progression was measured in an age-matched and BMI-matched cohort selected from the population-based Bruneck study to compare with changes in C‑IMT progression after bariatric surgery.

RESULTS

C‑IMT remained stable over the 10-year observation period after bariatric surgery. The control cohort showed a significant C‑IMT progression over 10 years. The difference in C‑IMT progression over 10 years was significant (p < 0.01) between both cohorts.

CONCLUSION

Weight loss induced by bariatric surgery halts the natural progression of C‑IMT over a 10-year observation period.

KLKB1 and CLSTN2 are associated with HDL-mediated cholesterol efflux capacity in a genome-wide association study

BACKGROUND AND AIMS

HDL-mediated cholesterol efflux capacity (CEC) may protect from cardiovascular disease. Thus, we aimed to identify its genetic and non-genetic determinants.

METHODS

We measured CEC to 2% apolipoprotein B-depleted serum using BODIPY-cholesterol and cAMP-stimulated J774A.1 macrophages using serum samples from 4,981 participants in the German Chronic Kidney Disease (GCKD) study. Variance of CEC explained by clinical and biochemical parameters in a multivariable linear regression model was calculated by proportional marginal variance decomposition. A genome-wide association study with 7,746,917 variants was performed based on an additive genetic model. The main model was adjusted for age, sex and principal components 1-10. Further models were selected for sensitivity analysis and to reduce residual variance by known CEC pathways.

RESULTS

Variables that explained 1% and more of the variance of CEC were concentrations of triglycerides (12.9%), HDL-cholesterol (11.8%), LDL-cholesterol (3.0%), apolipoprotein A-IV (2.8%), PCSK9 (1.0%), and eGFR (1.0%). The KLKB1 (chr4) and APOE/C1 (chr19) loci were genome-wide significantly (p < 5x10-8) associated with CEC in our main model (p = 8.8x10-10 and p = 3.3x10-10, respectively). KLKB1 remained significantly associated after additional adjustment for either kidney parameters, HDL-cholesterol, triglycerides or apolipoprotein A-IV concentrations, while the APOE/C1 locus was not significantly associated anymore after adjustment for triglycerides. Adjustment for triglycerides also revealed an association with the CLSTN2 locus (chr3; p = 6.0x10-9).

CONCLUSIONS

We identified HDL-cholesterol and triglycerides as the main determinants of CEC. Furthermore, we newly found a significant association of CEC with the KLKB1 and the CLSTN2 locus and confirmed the association with the APOE/C1 locus, likely mediated by triglycerides.

Implicating genes, pleiotropy, and sexual dimorphism at blood lipid loci through multi-ancestry meta-analysis

BACKGROUND

Genetic variants within nearly 1000 loci are known to contribute to modulation of blood lipid levels. However, the biological pathways underlying these associations are frequently unknown, limiting understanding of these findings and hindering downstream translational efforts such as drug target discovery.

RESULTS

To expand our understanding of the underlying biological pathways and mechanisms controlling blood lipid levels, we leverage a large multi-ancestry meta-analysis (N = 1,654,960) of blood lipids to prioritize putative causal genes for 2286 lipid associations using six gene prediction approaches. Using phenome-wide association (PheWAS) scans, we identify relationships of genetically predicted lipid levels to other diseases and conditions. We confirm known pleiotropic associations with cardiovascular phenotypes and determine novel associations, notably with cholelithiasis risk. We perform sex-stratified GWAS meta-analysis of lipid levels and show that 3-5% of autosomal lipid-associated loci demonstrate sex-biased effects. Finally, we report 21 novel lipid loci identified on the X chromosome. Many of the sex-biased autosomal and X chromosome lipid loci show pleiotropic associations with sex hormones, emphasizing the role of hormone regulation in lipid metabolism.

CONCLUSIONS

Taken together, our findings provide insights into the biological mechanisms through which associated variants lead to altered lipid levels and potentially cardiovascular disease risk.

Relative Telomere Length Is Associated With Age-Related Macular Degeneration in Women

Purpose

Relative telomere length (RTL) is a biomarker for physiological aging. Premature shortening of telomeres is associated with oxidative stress, which is one possible pathway that might contribute to age-related macular degeneration (AMD). We therefore aimed to investigate the association between RTL and AMD in a well-characterized group of elderly individuals.

Methods

We measured RTL in participants of the AugUR study using a multiplex quantitative PCR-based assay determining the ratio between the telomere product and a single-copy gene product (T/S ratio). AMD was assessed by manual grading of color fundus images using the Three Continent AMD Consortium Severity Scale.

Results

Among the 2262 individuals 70 to 95 years old (627 with AMD and 1635 without AMD), RTL was significantly shorter in individuals with AMD compared to AMD-free participants. In age- and sex-adjusted logistic regression analyses, we observed an 8% higher odds for AMD per 0.1 unit shorter RTL (odds ratio [OR] = 1.08; 95% confidence interval [CI], 1.02–1.14; P = 0.005). The estimates remained stable when adjusted for smoking, high-density lipoprotein cholesterol, cardiovascular disease, diabetes, and hypertension. Interestingly, this association was only present in women (OR = 1.14; 95% CI, 1.06–1.23; P < 0.001), but not in men (OR = 1.01; 95% CI, 0.93–1.10; P = 0.76). A significant sex-by-RTL interaction on AMD was detected (P = 0.043).

Conclusions

Our results show an association of RTL with AMD that was restricted to women. This is in line with altered reactive oxygen species levels and higher telomerase activity in women and provides an indication for a sex-differential pathway for oxidative stress and AMD.

Apolipoprotein A-IV concentrations and clinical outcomes in a large chronic kidney disease cohort: Results from the GCKD study

BACKGROUND

Chronic kidney disease (CKD) represents a chronic proinflammatory state and is associated with very high cardiovascular risk. Apolipoprotein A-IV (apoA-IV) has antiatherogenic, antioxidative, anti-inflammatory and antithrombotic properties and levels increase significantly during the course of CKD.

OBJECTIVES

We aimed to investigate the association between apoA-IV and all-cause mortality and cardiovascular outcomes in the German Chronic Kidney Disease study.

METHODS

This was a prospective cohort study including 5141 Caucasian patients with available apoA-IV measurements and CKD. The majority of the patients had an estimated glomerular filtration rate (eGFR) of 30-60 ml/min/1.73m² or an eGFR >60 ml/min/1.73m² in the presence of overt proteinuria. Median follow-up was 6.5 years. The association of apoA-IV with comorbidities at baseline and endpoints during follow-up was modelled adjusting for major confounders.

RESULTS

Mean apoA-IV concentrations of the entire cohort were 28.9 ± 9.8 mg/dl. Patients in the highest apoA-IV quartile had the lowest high-sensitivity C-reactive protein values despite the highest prevalence of diabetes, albuminuria and the lowest eGFR. Each 10 mg/dl higher apoA-IV translated into lower odds of prevalent cardiovascular disease (1289 cases, odds ratio = 0.80, 95% confidence interval [CI] 0.72-0.86, p = 0.0000003). During follow-up, each 10 mg/dl higher apoA-IV was significantly associated with a lower risk for all-cause mortality (600 cases, hazard ratio [HR] = 0.81, 95% CI 0.73-0.89, p = 0.00004), incident major adverse cardiovascular events (506 cases, HR = 0.88, 95% CI 0.79-0.99, p = 0.03) and death or hospitalizations due to heart failure (346 cases, HR = 0.84, 95% CI 0.73-0.96, p = 0.01).

CONCLUSIONS

These data support a link between elevated apoA-IV concentrations and reduced inflammation in moderate CKD. ApoA-IV appears to be an independent risk marker for reduced all-cause mortality, cardiovascular events and heart failure in a large cohort of patients with CKD.

The effect of LPA Thr3888Pro on lipoprotein(a) and coronary artery disease is modified by the LPA KIV-2 variant 4925G>A

BACKGROUND AND AIMS

High lipoprotein(a) [Lp(a)] concentrations are associated with increased coronary artery disease (CAD) risk. Lp(a) is regulated mainly genetically by the LPA gene but involved genetic variants have not been fully elucidated. Improved understanding of the entanglements of genetic Lp(a) regulation may enhance genetic prediction of Lp(a) and CAD risk. We investigated an interaction between the well-known LPA missense SNP rs41272110 (known as Thr3888Pro) and the frequent LPA splicing mutation KIV-2 4925G>A.

METHODS

Effects on Lp(a) concentrations were investigated by multiple quantile regression in the German Chronic Kidney Disease (GCKD) study, KORA-F3 and KORA-F4 (n_total = 10,405) as well as in the UK Biobank (UKB) 200k exome dataset (n = 173,878). The impact of the interaction on CAD risk was assessed by survival analysis in UKB.

RESULTS

We observed a significant SNP-SNP interaction in all studies (p = 1.26e-05 to 3.03e-04). In quantile regression analysis, rs41272110 as a predictor shows no impact on Lp(a) (β = -0.06 [-0.79; 0.68], p = 0.879), but in a joint model including both SNPs as predictors, rs41272110 is associated with markedly higher Lp(a) (β = +9.40 mg/dL [6.45; 12.34], p = 4.07e-10). Similarly, rs41272110 shows no effect on CAD in UKB (HR = 1.01 [0.97; 1.04], p = 0.731), while rs41272110 carriers not carrying 4925G>A show an increased CAD risk (HR = 1.10 [1.04; 1.16], p = 6.9e-04). This group corresponds to 4% of the population. Adjustment for apolipoprotein(a) isoforms further modified the effect estimates markedly.

CONCLUSIONS

This work emphasizes the complexity of the genetic regulation of Lp(a) and the importance to account for genetic subgroups in Lp(a) association studies and when interpreting genetic cardiovascular risk profiles.

A Family and a Genome-Wide Polygenic Risk Score Are Independently Associated With Stroke in a Population-Based Study

BACKGROUND

Positive family history and genetic risk scores have been shown to independently capture those individuals with high risk for stroke. The aim of our study was to evaluate the amount of shared information between family history and genetic risk and to investigate their combined effect on the association with prevalent and incident stroke cases.

METHODS

We obtained a family risk score (FamRS), weighted for disease onset and family size as well as genome-wide polygenic risk score (PGS) including over 3.2 million single-nucleotide polymorphisms in the population-based prospective KORA F3 (Cooperative Health Research in the Region of Augsburg) study (n=3071) from Southern Germany. FamRS and PGS were evaluated separately and combined. The measures were once treated as continuous variables but also divided in the highest 20%, 10%, 5%, and 1% percentiles. Odds ratios via logistic regression and hazard ratios via Cox regression were estimated. A stroke event was defined as a hospitalization for stroke that was self-reported in a standardized interview by certified and supervised personnel.

RESULTS

The FamRS outperformed other simplified family measures such as affected parents or number of affected family members. FamRS and PGS were not correlated, and no individuals were observed with both very high FamRS and very high PGS (top 1% percentile). In a combined model, both FamRS and PGS were independently from each other associated with risk of stroke, also independent of other traditional risk factors (p [FamRS]=0.02, p [PGS]=0.005). Individuals in the top 1% of either FamRS or PGS were found to have >5-fold risk for stroke (odds ratios, 5.82 [95% CI, 2.08-14]; P=0.0002). The results for incident stroke events showed the same trend but were not significant.

CONCLUSIONS

Our study shows that a family risk score and PGS capture different information concerning individual stroke risk. Combining the risk measures FamRS and PGS increases predictive power, as demonstrated in a population-based study.

PCSK9 and Cardiovascular Disease in Individuals with Moderately Decreased Kidney Function

Background and objectives: Proprotein convertase subtilisin/kexin type 9 (PCSK9) is a key regulator of lipid homeostasis. Studies investigating the association between PCSK9 and cardiovascular disease in large cohorts of CKD patients are limited. Design, setting, participants, and measurements: The association of PCSK9 concentrations with prevalent and incident cardiovascular disease was investigated in 5,138 Caucasians of the German Chronic Kidney Disease study with a median follow-up of 6.5 years. Inclusion criteria were eGFR of 30-60 mL/min/1.73 m² or an eGFR>60 mL/min/1.73m² in the presence of overt proteinuria (urine albumin-creatinine ratio >300mg/g or equivalent). Prevalent cardiovascular disease was defined as history of non-fatal myocardial infarction, coronary artery bypass grafting, percutaneous transluminal coronary angioplasty, carotid arteries interventions and stroke. Incident major adverse cardiovascular disease events included death from cardiovascular causes, acute non-fatal myocardial infarction and non-fatal stroke. Results: Median PCSK9 concentration in the cohort was 285 ng/mL (interquartile range 231-346 ng/mL). There was no association between PCSK9 concentrations and baseline eGFR and albuminuria. With each 100 ng/ml increment of PCSK9, the odds for prevalent cardiovascular disease (n=1,284) was 1.22-fold (95%CI 1.12-1.34, p<0.001) higher in a model with extended adjustment for major confounders. This association was stronger in non-statin than statin-users (p-value for interaction=0.009). During follow-up 474 individuals experienced a major adverse cardiovascular disease event and participants in PCSK9 quartiles 2 to 4 had a 32%-47% higher risk compared to quartile 1 (p<0.05). Subgroup analysis revealed this association was restricted to those participants who already had cardiovascular disease at baseline (all HR>1.75, p<0.05). In addition, PCSK9 showed a valuable gain in classification accuracy both for prevalent cardiovascular disease (NRI=0.27, 95%CI:0.20-0.33) and incident major adverse cardiovascular disease events during follow-up (NRI=0.10, 95%CI:0.008-0.21) when added to an extended adjustment model. Conclusions: Our findings reveal no relation of PCSK9 with baseline eGFR and albuminuria, but a significant association between higher PCSK9 concentrations and risk of cardiovascular disease independent of traditional risk factors including LDL-cholesterol levels.

Genome-Wide Characterization of a Highly Penetrant Form of Hyperlipoprotein(a)emia Associated With Genetically Elevated Cardiovascular Risk

Lp(a) (lipoprotein [a]) is a highly atherogenic lipoprotein strongly associated with coronary artery disease (CAD). Lp(a) concentrations are chiefly determined genetically. Investigation of large pedigrees with extreme Lp(a) using modern whole-genome approaches may unravel the genetic determinants underpinning this pathological phenotype.

Cis-epistasis at the LPA locus and risk of cardiovascular diseases

AIMS

Coronary artery disease (CAD) has a strong genetic predisposition. However, despite substantial discoveries made by genome-wide association studies (GWAS), a large proportion of heritability awaits identification. Non-additive genetic effects might be responsible for part of the unaccounted genetic variance. Here, we attempted a proof-of-concept study to identify non-additive genetic effects, namely epistatic interactions, associated with CAD.

METHODS AND RESULTS

We tested for epistatic interactions in 10 CAD case-control studies and UK Biobank with focus on 8068 SNPs at 56 loci with known associations with CAD risk. We identified a SNP pair located in cis at the LPA locus, rs1800769 and rs9458001, to be jointly associated with risk for CAD [odds ratio (OR) = 1.37, P = 1.07 × 10-11], peripheral arterial disease (OR = 1.22, P = 2.32 × 10-4), aortic stenosis (OR = 1.47, P = 6.95 × 10-7), hepatic lipoprotein(a) (Lp(a)) transcript levels (beta = 0.39, P = 1.41 × 10-8), and Lp(a) serum levels (beta = 0.58, P = 8.7 × 10-32), while individual SNPs displayed no association. Further exploration of the LPA locus revealed a strong dependency of these associations on a rare variant, rs140570886, that was previously associated with Lp(a) levels. We confirmed increased CAD risk for heterozygous (relative OR = 1.46, P = 9.97 × 10-32) and individuals homozygous for the minor allele (relative OR = 1.77, P = 0.09) of rs140570886. Using forward model selection, we also show that epistatic interactions between rs140570886, rs9458001, and rs1800769 modulate the effects of the rs140570886 risk allele.

CONCLUSIONS

These results demonstrate the feasibility of a large-scale knowledge-based epistasis scan and provide rare evidence of an epistatic interaction in a complex human disease. We were directed to a variant (rs140570886) influencing risk through additive genetic as well as epistatic effects. In summary, this study provides deeper insights into the genetic architecture of a locus important for cardiovascular diseases.

The Manchester Triage System in a Pediatric Emergency Department of an Austrian University Hospital: A Retrospective Analysis of Urgency Levels

OBJECTIVES

The Manchester Triage System (MTS) has entered widespread international use in emergency departments (EDs). This retrospective study analyzes urgency of patient visits (PV) at the ED of the Clinic for Pediatrics at the Medical University of Innsbruck.

METHODS

We collected demographic and outcome information, including PV urgency levels (UL) according to the MTS, for 3 years (2015-2018), separating PV during regular office hours (ROH; 8:00 am to 5:00 pm) from PV during afternoon and night hours (5:00 pm to 8:00 am), and PV on weekdays from PV on weekends and bank holidays (WE).

RESULTS

A total of 56,088 PV were registered with a UL. Most (68.4%) PV were classified as nonurgent. During ROH, more PV per hour (PV/h) were recorded than during afternoon and night hours (3.0 PV/h vs 1.6 PV/h), with a higher proportion of less urgent cases during ROH. On WE, the amount of PV/h was higher than on weekdays (3.6 PV/h vs 2.8 PV/h), with a higher proportion of nonurgent cases (74.6% vs 68.6%). Likelihoods of inpatient admission and hospital stay lengths increased in step with UL.

CONCLUSIONS

The MTS proved useful for delineating UL distributions. The MTS analyses may be of value in managing EDs. Prompted by the results of our study, a general practice pediatric care unit was established to support the ED during WE.

Lipoprotein(a) and SARS-CoV-2 infections: Susceptibility to infections, ischemic heart disease and thromboembolic events

BACKGROUND

Comorbidities including ischemic heart disease (IHD) worsen outcomes after SARS-CoV-2 infections. High lipoprotein(a) [Lp(a)] concentrations are a strong risk factor for IHD and possibly for thromboembolic events. We therefore evaluated whether SARS-CoV-2 infections modify the risk of high Lp(a) concentrations for IHD or thromboembolic events during the first 8.5 months follow-up of the pandemic.

METHOD

Cohort study using data from the UK Biobank during the SARS-CoV-2 pandemic. Baseline Lp(a) was compared between SARS-CoV-2 positive patients and the population controls.

RESULTS

SARS-CoV-2 positive patients had Lp(a) concentrations similar to the population controls. The risk for IHD increased with higher Lp(a) concentrations in both, the population controls (n = 435,104) and SARS-CoV-2 positive patients (n = 6937). The causality of the findings was supported by a genetic risk score for Lp(a). A SARS-CoV-2 infection modified the association with a steeper increase in risk for infected patients (interaction p-value = 0.03). Although SARS-CoV-2 positive patients had a five-times higher frequency of thromboembolic events compared to the population controls (1.53% vs. 0.31%), the risk was not influenced by Lp(a).

CONCLUSIONS

SARS-CoV-2 infections enforce the association between high Lp(a) and IHD but the risk for thromboembolic events is not influenced by Lp(a).

Afamin predicts the prevalence and incidence of nonalcoholic fatty liver disease

OBJECTIVES

In the general population, increased afamin concentrations are associated with the prevalence and incidence of metabolic syndrome as well as type 2 diabetes. Although metabolic syndrome is commonly associated with nonalcoholic fatty liver disease (NAFLD), there exist no information on afamin and NAFLD.

METHODS

Afamin concentrations were cross-sectionally measured in 146 Austrian patients with NAFLD, in 45 patients without NAFLD, and in 292 age- and sex-matched healthy controls. Furthermore, the feasibility of afamin to predict incident NAFLD was evaluated in 1,434 adult participants in the population-based Cardiovascular Risk in Young Finns Study during a 10-year follow-up.

RESULTS

Median afamin concentrations were significantly higher in NAFLD patients (83.6 mg/L) than in patients without NAFLD (61.6 mg/L, p<0.0001) or in healthy controls (63.9 mg/L, p<0.0001). In age- and sex-adjusted logistic regression analyses a 10 mg/L increase of afamin was associated with a 1.5-fold increase of having NAFLD as compared with patients without NAFLD and the risk was even two-fold when compared with healthy controls. In the population-based cohort, afamin concentrations at baseline were significantly lower in participants without NAFLD (n=1,195) than in 239 participants who developed NAFLD (56.5 vs. 66.9 mg/L, p<0.0001) during the 10-year follow up, with highest afamin values observed in individuals developing severe forms of NAFLD. After adjustment for several potentially confounding parameters, afamin remained an independent predictor for the development of NAFLD (OR=1.37 [95% CI 1.23-1.54] per 10 mg/L increase, p<0.0001).

CONCLUSIONS

Afamin concentrations are increased in patients with NAFLD and independently predict the development of NAFLD in a population-based cohort.

Association of the metabolic syndrome with mortality and major adverse cardiac events: A large chronic kidney disease cohort

BACKGROUND

Metabolic syndrome with its key components insulin resistance, central obesity, dyslipidaemia, and hypertension is associated with a high risk for cardiovascular events and all-cause mortality in the general population. However, evidence that these findings apply to patients with chronic kidney disease (CKD) with moderately reduced estimated glomerular filtration rate and/or albuminuria is limited.

OBJECTIVES

We aimed to investigate the association between metabolic syndrome and its components with all-cause mortality and cardiovascular outcomes in CKD patients.

METHODS

Prospective observation of a cohort of 5110 CKD patients from the German Chronic Kidney Disease study with 3284 (64.3%) of them having a metabolic syndrome at baseline.

RESULTS

During the follow-up of 6.5 years, 605 patients died and 650 patients experienced major cardiovascular events. After extended data adjustment, patients with a metabolic syndrome had a higher risk for all-cause mortality (hazard ratio [HR] = 1.26, 95% confidence interval [CI]: 1.04-1.54) and cardiovascular events (HR = 1.48, 95% CI: 1.22-1.79). The risk increased steadily with a growing number of metabolic syndrome components (increased waist circumference, glucose, triglycerides, hypertension and decreased HDL cholesterol): HR per component = 1.09 (95% CI: 1.02-1.17) for all-cause mortality and 1.23 (95% CI: 1.15-1.32) for cardiovascular events. This resulted in hazard ratios between 1.50 and 2.50 in the case when four or five components are present. An analysis of individual components of metabolic syndrome showed that the glucose component led to the highest increase in risk for all-cause mortality (HR = 1.68, 95% CI: 1.38-2.03) and cardiovascular events (HR = 1.81, 95% CI: 1.51-2.18), followed by the HDL cholesterol and triglyceride components.

CONCLUSIONS

We observed a high prevalence of metabolic syndrome among patients with moderate CKD. Metabolic syndrome increases the risk for all-cause mortality and cardiovascular events. The glucose and lipid components seem to be the main drivers for the association with outcomes.

The power of genetic diversity in genome-wide association studies of lipids

Increased blood lipid levels are heritable risk factors of cardiovascular disease with varied prevalence worldwide owing to different dietary patterns and medication use1. Despite advances in prevention and treatment, in particular through reducing low-density lipoprotein cholesterol levels2, heart disease remains the leading cause of death worldwide3. Genome-wideassociation studies (GWAS) of blood lipid levels have led to important biological and clinical insights, as well as new drug targets, for cardiovascular disease. However, most previous GWAS4-23 have been conducted in European ancestry populations and may have missed genetic variants that contribute to lipid-level variation in other ancestry groups. These include differences in allele frequencies, effect sizes and linkage-disequilibrium patterns24. Here we conduct a multi-ancestry, genome-wide genetic discovery meta-analysis of lipid levels in approximately 1.65 million individuals, including 350,000 of non-European ancestries. We quantify the gain in studying non-European ancestries and provide evidence to support the expansion of recruitment of additional ancestries, even with relatively small sample sizes. We find that increasing diversity rather than studying additional individuals of European ancestry results in substantial improvements in fine-mapping functional variants and portability of polygenic prediction (evaluated in approximately 295,000 individuals from 7 ancestry groupings). Modest gains in the number of discovered loci and ancestry-specific variants were also achieved. As GWAS expand emphasis beyond the identification of genes and fundamental biology towards the use of genetic variants for preventive and precision medicine25, we anticipate that increased diversity of participants will lead to more accurate and equitable26 application of polygenic scores in clinical practice.

Prognostic value of depressed cardiac index after STEMI: a phase-contrast magnetic resonance study

AIMS 

An invasively measured cardiac index (CI) of ≤2.2 L/min/m2 is one of the strongest prognostic indicators after ST-elevation myocardial infarction (STEMI), however, knowledge is mainly based on invasive evaluations performed in the pre-stent era. Velocity-encoded phase-contrast cardiac magnetic resonance (PC-CMR) allows non-invasive determination of CI.

METHODS AND RESULTS 

In this prospective study, CMR was performed in 406 stable and contemporarily revascularized patients a median of 3 days after STEMI. Forward stroke volume was assessed at the level of the ascending aorta by PC-CMR. Left ventricular ejection fraction (LVEF) and global longitudinal strain (GLS) were determined by cine CMR. Major adverse cardiac events (MACE) were defined as the composite of death, myocardial infarction, or hospitalization for heart failure. Median CI was 2.52 L/min/m2 and 27% of patients had ≤2.2 L/min/m2. Median LVEF was 53% and median GLS was -12.2%. During a median follow-up of 14.2 [95% confidence interval (95% CI) 13.6-14.7] months, 41 patients (10.1%) experienced a MACE. A depressed CI was significantly associated with MACE after adjustment for LVEF, GLS, Thrombolysis in Myocardial Infarction (TIMI) risk score, and infarct size [hazard ratio = 3.15 (95% CI 1.53-6.47); P = 0.002] and led to significant discrimination improvement [net reclassification improvement 0.61 (95% CI 0.25-0.97); P < 0.001].

CONCLUSIONS 

A CI of 2.2 L/min/m2 or less as measured by PC-CMR was present in 27% of clinically stable patients after STEMI and strongly and independently predicted medium-term MACE. The prognostic value of a depressed CI was superior and incremental to LVEF, GLS, TIMI risk score, and infarct size.

Time-dependent lipid profile inversely associates with mortality in hemodialysis patients - independent of inflammation/malnutrition

BACKGROUND

Patients with end-stage kidney disease have an extremely high cardiovascular mortality rate, but there is a paradoxical relationship between lipid profile and survival in haemodialysis patients. To investigate whether inflammation/malnutrition confounds the associations between lipids and mortality, we studied a full lipid profile comprising of five clinically well-established lipid parameters and its associations with mortality in a large, multinational European cohort with a median follow-up >3 years.

METHODS

The association between quartiles of total, high-density lipoprotein (HDL), non-HDL, low-density lipoprotein (LDL) cholesterol, as well as triglyceride, levels and the end-points of all-cause, cardiovascular and non-cardiovascular mortality was assessed in a cohort of 5,382 incident, adult haemodialysis patients from >250 Fresenius Medical Care dialysis centres out of 14 participating countries using baseline and time-dependent Cox models. Analyses were fully adjusted and stratified for inflammation/malnutrition status and other patient-level variables.

RESULTS

Time-dependent quartiles of total, HDL, non-HDL and LDL cholesterol were inversely associated with the hazard for all-cause, cardiovascular and non-cardiovascular mortality. Compared with the lowest quartile of the respective lipid parameter, hazard ratios of other quartiles were <0.86. Similar, albeit weaker, associations were found with baseline lipid profile and mortality. Neither time-dependent nor baseline associations between lipid profile and mortality were affected by inflammation/malnutrition, statin use or geography.

CONCLUSIONS

Baseline and time-dependent lipid profile are inversely associated with mortality in a large, multicentre cohort of incident haemodialysis patients. Inflammation/malnutrition is not a confounder nor effect modificator of the associations between lipid profile and mortality in European haemodialysis patients.

Lysis reagents, cell numbers, and calculation method influence high-throughput measurement of HDL-mediated cholesterol efflux capacity

HDL-mediated cholesterol efflux capacity (CEC) may protect against cardiovascular disease. However, CEC assays are not standardized, hampering their application in large cohorts and comparison between studies. To improve standardization, we systematically investigated technical differences between existing protocols that influence assay performance that have not been previously addressed. CEC was measured in 96-well plates using J774A.1 macrophages labeled with BODIPY-cholesterol and incubated for 4 h with 2% apolipoprotein B-depleted human serum. The time zero method, which calculates CEC using control wells, and the per-well method, which calculates CEC based on the actual content of BODIPY-cholesterol in each well, were compared in 506 samples. We showed that the per-well method had a considerably lower sample rejection rate (4.74% vs. 13.44%) and intra-assay (4.48% vs. 5.28%) and interassay coefficients of variation (two controls: 7.85%, 9.86% vs. 13.58%, 15.29%) compared with the time zero method. Correction for plate-to-plate differences using four controls on each plate also improved assay performance of both methods. In addition, we observed that the lysis reagent used had a significant effect. Compared with cholic acid, lysis with sodium hydroxide results in higher (P = 0.0082) and Triton X-100 in lower (P = 0.0028) CEC values. Furthermore, large cell seeding errors (30% variation) greatly biased CEC for both referencing methods (P < 0.0001) as measured by a resazurin assay. In conclusion, lysis reagents, cell numbers, and assay setup greatly impact the quality and reliability of CEC quantification and should be considered when this method is newly established in a laboratory.

Frequent LPA KIV-2 Variants Lower Lipoprotein(a) Concentrations and Protect Against Coronary Artery Disease

Background

Lipoprotein(a) (Lp(a)) concentrations are a major independent risk factor for coronary artery disease (CAD) and are mainly determined by variation in LPA. Up to 70% of the LPA coding sequence is located in the hyper-variable kringle IV type 2 (KIV-2) region. It is hardly accessible by conventional technologies, but may contain functional variants.

Objectives

This study sought to investigate the new, very frequent splicing variant KIV-2 4733G>A on Lp(a) and CAD.

Methods

We genotyped 4733G>A in the GCKD (German Chronic Kidney Disease) study (n = 4,673) by allele-specific polymerase chain reaction, performed minigene assays, identified proxy single nucleotide polymorphisms and used them to characterize its effect on CAD by survival analysis in UK Biobank (n = 440,234). Frequencies in ethnic groups were assessed in the 1000 Genomes Project.

Results

The 4733G>A variant (38.2% carrier frequency) was found in most isoform sizes. It reduces allelic expression without abolishing protein production, lowers Lp(a) by 13.6 mg/dL (95% CI: 12.5-14.7; P < 0.0001) and is the strongest variance-explaining factor after the smaller isoform. Splicing of minigenes was modified. Compound heterozygosity (4.6% of the population) for 4733G>A and 4925G>A, another KIV-2 splicing mutation, reduces Lp(a) by 31.8 mg/dL and most importantly narrows the interquartile range by 9-fold (from 42.1 to 4.6 mg/dL) when compared to the wild type. In UK Biobank 4733G>A alone and compound heterozygosity with 4925G>A reduced HR for CAD by 9% (95% CI: 7%-11%) and 12% (95% CI: 7%-16%) (both P < 0.001). Frequencies in ethnicities differ notably.

Conclusions

Functional variants in the previously inaccessible LPA KIV-2 region cooperate in determining Lp(a) variance and CAD risk. Even a moderate but lifelong genetic Lp(a) reduction translates to a noticeable CAD risk reduction. (J Am Coll Cardiol 2021;78:437–49)

The size of apolipoprotein (a) is an independent determinant of the reduction in lipoprotein (a) induced by PCSK9 inhibitors

AIMS

Lipoprotein (a) [Lp(a)] is a lipoprotein species causatively associated with atherosclerosis. Unlike statins, PCSK9 inhibitors (PCSK9i) reduce Lp(a), but this reduction is highly variable. Levels of Lp(a) are chiefly governed by the size of its signature protein, apolipoprotein (a) [apo(a)]. Whether this parameter determines some of the reduction in Lp(a) induced by PCSK9i remains unknown. We aimed to investigate if the Lp(a) lowering efficacy of PCSK9i is modulated by the size of apo(a), which is genetically determined by the variable number of KIV domains present on that protein.

METHODS AND RESULTS

The levels of Lp(a) and the size of apo(a) were assessed in plasma samples from 268 patients before and after treatment with PCSK9i. Patients were recruited at the Outpatient Lipid Clinic of the Charité Hospital (Berlin) between 2015 and 2020. They were hypercholesterolemic at very high CVD risk with LDL-cholesterol levels above therapeutic targets despite maximally tolerated lipid-lowering therapy. Patients received either Alirocumab (75 or 150 mg) or Evolocumab (140 mg) every 2 weeks. Apo(a), apoB100, and apoE concentrations as well as apoE major isoforms were determined by liquid chromatography high-resolution mass spectrometry. Apo(a) isoforms sizes were determined by Western Blot. PCSK9i sharply reduced LDL-cholesterol (-57%), apoB100 (-47%) and Lp(a) (-36%). There was a positive correlation between the size of apo(a) and the relative reduction in Lp(a) induced by PCSK9i (r = 0.363, p = 0.0001). The strength of this association remained unaltered after adjustment for baseline Lp(a) levels and all other potential confounding factors. In patients with two detectable apo(a) isoforms, there was also a positive correlation between the size of apo(a) and the reduction in Lp(a), separately for the smaller (r = 0.350, p = 0.0001) and larger (r = 0.324, p = 0.0003) isoforms. The relative contribution of the larger isoform to the total concentration of apo(a) was reduced from 29% to 15% (p < 0.0001).

CONCLUSIONS

The size of apo(a) is an independent determinant of the response to PCSK9i. Each additional kringle domain is associated with a 3% additional reduction in Lp(a). This explains in part the variable efficacy of PCSK9i and allows to identify patients who will benefit most from these therapies in terms of Lp(a) lowering.

TRANSLATIONAL PERSPECTIVE

Unlike statins, PCSK9 inhibitors reduce the circulating levels of the highly atherogenic Lipoprotein (a). The underlying mechanism remains a matter of considerable debate. The size of apo(a), the signature protein of Lp(a), is extremely variable (300 to more than 800 kDa) and depends on its number of kringle domains. We now show that each increase in apo(a) size by one kringle domain is associated with a 3% additional reduction in Lp(a) following PCSK9i treatment and that apo(a) size polymorphism is an independent predictor of the reduction in Lp(a) induced by these drugs. In an era of personalized medicine, this allows to identify patients who will benefit most from PCSK9i in terms of Lp(a) lowering.

Association of mitochondrial DNA copy number with metabolic syndrome and type 2 diabetes in 14 176 individuals

BACKGROUND

Mitochondria play an important role in cellular metabolism, and their dysfunction is postulated to be involved in metabolic disturbances. Mitochondrial DNA is present in multiple copies per cell. The quantification of mitochondrial DNA copy number (mtDNA-CN) might be used to assess mitochondrial dysfunction.

OBJECTIVES

We aimed to investigate the cross-sectional association of mtDNA-CN with type 2 diabetes and the potential mediating role of metabolic syndrome.

METHODS

We examined 4812 patients from the German Chronic Kidney Disease (GCKD) study and 9364 individuals from the Cooperative Health Research in South Tyrol (CHRIS) study. MtDNA-CN was measured in whole blood using a plasmid-normalized qPCR-based assay.

RESULTS

In both studies, mtDNA-CN showed a significant correlation with most metabolic syndrome parameters: mtDNA-CN decreased with increasing number of metabolic syndrome components. Furthermore, individuals with low mtDNA-CN had significantly higher odds of metabolic syndrome (OR = 1.025; 95% CI = 1.011-1.039, P = 3.19 × 10-4 , for each decrease of 10 mtDNA copies) and type 2 diabetes (OR = 1.027; 95% CI = 1.012-1.041; P = 2.84 × 10-4 ) in a model adjusted for age, sex, smoking and kidney function in the meta-analysis of both studies. Mediation analysis revealed that the association of mtDNA-CN with type 2 diabetes was mainly mediated by waist circumference in the GCKD study (66%) and by several metabolic syndrome parameters, especially body mass index and triglycerides, in the CHRIS study (41%).

CONCLUSIONS

Our data show an inverse association of mtDNA-CN with higher risk of metabolic syndrome and type 2 diabetes. A major part of the total effect of mtDNA-CN on type 2 diabetes is mediated by obesity parameters.

Elevated levels of serum PCSK9 in male patients with symptomatic peripheral artery disease: The CAVASIC study

BACKGROUND AND AIMS

Peripheral artery disease (PAD) affects more than 200 million people worldwide. Increased low-density lipoprotein cholesterol (LDL-C)levels are a risk factor for PAD and the concentrations are influenced by proprotein convertase subtilisin/kexin type 9 (PCSK9). PCSK9 regulates the recycling of the LDL receptors to the cell membrane surface. Only a limited number of mostly small studies investigated the association between serum PCSK9 concentrations and PAD of different definition, which revealed contrasting results.

METHODS

Serum PCSK9, lipoprotein(a) [Lp(a)] and other lipoprotein concentrations were measured in male participants of the CAVASIC study, a case-control study of 248 patients with intermittent claudication and 251 age and diabetes-matched controls.

RESULTS

PAD patients had significantly higher PCSK9 concentrations when compared to controls (250 ± 77 vs. 222 ± 68 ng/mL, p < 0.001). Logistic regression analysis with adjustment for age revealed that an increase in PCSK9 concentrations of 100 ng/mL was associated with a 1.78-fold higher risk for PAD (95%CI 1.38-2.33, p = 1.43 × 10^-5). The association attenuated, but was still significant when adjusting additionally for age, Lp(a)-corrected LDL cholesterol, HDL cholesterol, high-sensitivity-CRP, statin treatment, hypertension, diabetes mellitus and smoking (OR = 1.49, 95%CI 1.03-2.18, p = 0.035). The strongest association was observed when both PCSK9 concentrations were above the median and Lp(a) concentrations were above 30 mg/dL (OR = 3.35, 95%CI 1.49-7.71, p = 0.0038).

CONCLUSIONS

Our findings suggest an association of higher PCSK9 concentrations with PAD, which was independent of other lipid parameters and classical cardiovascular risk factors.

HDL cholesterol efflux capacity is inversely associated with subclinical cardiovascular risk markers in young adults: The cardiovascular risk in Young Finns study

The atherogenic process begins already in childhood and progresses to symptomatic condition with age. We investigated the association of cholesterol efflux capacity (CEC) and vascular markers of subclinical atherosclerosis in healthy, young adults. CEC was determined in 2282 participants of the Young Finns study using cAMP treated ³H-cholesterol-labeled J774 cells. The CEC was correlated to baseline and 6-year follow-up data of cardiovascular risk factors and ultrasound measurements of arterial structure and function. CEC was higher in women, correlated with total cholesterol, HDL-C, and apolipoprotein A-I, but not with LDL-C or apolipoprotein B. Compared to the lowest CEC quartile, the highest CEC quartile was significantly associated with high CRP levels and inversely associated with adiponectin. At baseline, high CEC was associated with decreased flow-mediated dilation (FMD) and carotid artery distensibility, as well as an increased Young's modulus of elasticity, indicating adverse changes in arterial structure, and function. The association reversed with follow-up FMD data, indicating the interaction of preclinical parameters over time. A higher CEC was directly associated with a lower risk of subclinical atherosclerosis at follow-up. In young and healthy subjects, CEC was associated with important lipid risk parameters at baseline, as in older patients and CAD patients, but inversely with early risk markers for subclinical atherosclerosis.

The haemochromatosis gene Hfe and Kupffer cells control LDL cholesterol homeostasis and impact on atherosclerosis development

AIMS

Imbalances of iron metabolism have been linked to the development of atherosclerosis. However, subjects with hereditary haemochromatosis have a lower prevalence of cardiovascular disease. The aim of our study was to understand the underlying mechanisms by combining data from genome-wide association study analyses in humans, CRISPR/Cas9 genome editing, and loss-of-function studies in mice.

METHODS AND RESULTS

Our analysis of the Global Lipids Genetics Consortium (GLGC) dataset revealed that single nucleotide polymorphisms (SNPs) in the haemochromatosis gene HFE associate with reduced low-density lipoprotein cholesterol (LDL-C) in human plasma. The LDL-C lowering effect could be phenocopied in dyslipidaemic ApoE-/- mice lacking Hfe, which translated into reduced atherosclerosis burden. Mechanistically, we identified HFE as a negative regulator of LDL receptor expression in hepatocytes. Moreover, we uncovered liver-resident Kupffer cells (KCs) as central players in cholesterol homeostasis as they were found to acquire and transfer LDL-derived cholesterol to hepatocytes in an Abca1-dependent fashion, which is controlled by iron availability.

CONCLUSION

Our results disentangle novel regulatory interactions between iron metabolism, KC biology and cholesterol homeostasis which are promising targets for treating dyslipidaemia but also provide a mechanistic explanation for reduced cardiovascular morbidity in subjects with haemochromatosis.

Association of changes in bone mineral parameters with mortality in haemodialysis patients: insights from the ARO cohort

BACKGROUND

There is little information in haemodialysis (HD) patients on whether temporal changes in serum calcium, phosphate or intact parathyroid hormone (iPTH) are associated with mortality.

METHODS

We analysed associations of phosphate, total calcium and iPTH with all-cause and cardiovascular mortality in 8817 incident HD patients from the European second Analyzing Data, Recognizing Excellence and Optimizing Outcomes (AROii) cohort enrolled in 2007-09, which were prospectively followed for a median of 3 years, using time-dependent Cox proportional hazards models. We evaluated changes in risk over time depending on changes in phosphate, calcium or iPTH.

RESULTS

The association of phosphate and iPTH with all-cause mortality was U-shaped, with the lowest risk ranges between 1.20 and 1.89 mmol/L for phosphate and between 239 and 710 ng/L for iPTH. For total calcium, the associations were J-shaped, with an increased risk for all-cause mortality at levels >2.36 mmol/L. Lowest risk ranges for cardiovascular mortality did not change markedly for all three parameters. If iPTH was below the lowest risk range at baseline (iPTH <239 ng/L), a subsequent increase in levels was associated with improved survival. For phosphate, an increase or decrease out of the lowest risk range was associated with increased mortality risk. For calcium, this was only the case when the values increased above the lowest risk range.

CONCLUSION

In the AROii cohort, the ranges of bone mineral biomarkers associated with the lowest mortality ranges were largely consistent with the current Kidney Disease: Improving Global Outcomes chronic kidney disease-mineral and bone disorder guideline recommendations. Allowing a suppressed iPTH to increase was associated with a lower mortality, whereas shifts of phosphate or calcium outside the lowest risk range increased mortality.

Investigation of a nonsense mutation located in the complex KIV-2 copy number variation region of apolipoprotein(a) in 10,910 individuals

Background

The concentrations of the highly atherogenic lipoprotein(a) [Lp(a)] are mainly genetically determined by the LPA gene locus. However, up to 70% of the coding sequence is located in the complex so-called kringle IV type 2 (KIV-2) copy number variation, a region hardly accessible by common genotyping and sequencing technologies. Despite its size, little is known about genetic variants in this complex region. The R21X variant is a functional variant located in this region, but it has never been analyzed in large cohorts.

Methods

We typed R21X in 10,910 individuals from three European populations using a newly developed high-throughput allele-specific qPCR assay. R21X allelic location was determined by separating the LPA alleles using pulsed-field gel electrophoresis (PFGE) and typing them separately. Using GWAS data, we identified a proxy SNP located outside of the KIV-2. Linkage disequilibrium was determined both statistically and by long-range haplotyping using PFGE. Worldwide frequencies were determined by reanalyzing the sequencing data of the 1000 Genomes Project with a dedicated pipeline.

Results

R21X carriers (frequency 0.016–0.021) showed significantly lower mean Lp(a) concentrations (− 11.7 mg/dL [− 15.5; − 7.82], p = 3.39e−32). The variant is located mostly on medium-sized LPA alleles. In the 1000 Genome data, R21X mostly occurs in Europeans and South Asians, is absent in Africans, and shows varying frequencies in South American populations (0 to 0.022). Of note, the best proxy SNP was another LPA null mutation (rs41272114, D′ = 0.958, R² = 0.281). D′ was very high in all 1000G populations (0.986–0.996), although rs41272114 frequency varies considerably (0–0.182). Co-localization of both null mutations on the same allele was confirmed by PFGE-based long-range haplotyping.

Conclusions

We performed the largest epidemiological study on an LPA KIV-2 variant so far, showing that it is possible to assess LPA KIV-2 mutations on a large scale. Surprisingly, in all analyzed populations, R21X was located on the same haplotype as the splice mutation rs41272114, creating “double-null” LPA alleles. Despite being a nonsense variant, the R21X status does not provide additional information beyond the rs41272114 genotype. This has important implications for studies using LPA loss-of-function mutations as genetic instruments and emphasizes the complexity of LPA genetics.

The relationship between blood metabolites of the tryptophan pathway and kidney function: a bidirectional Mendelian randomization analysis

Blood metabolites of the tryptophan pathway were found to be associated with kidney function and disease in observational studies. In order to evaluate causal relationship and direction, we designed a study using a bidirectional Mendelian randomization approach. The analyses were based on published summary statistics with study sizes ranging from 1,960 to 133,413. After correction for multiple testing, results provided no evidence of an effect of metabolites of the tryptophan pathway on estimated glomerular filtration rate (eGFR). Conversely, lower eGFR was related to higher levels of four metabolites: C-glycosyltryptophan (effect estimate = − 0.16, 95% confidence interval [CI] (− 0.22; − 0.1); p = 9.2e−08), kynurenine (effect estimate = − 0.18, 95% CI (− 0.25; − 0.11); p = 1.1e−06), 3-indoxyl sulfate (effect estimate = − 0.25, 95% CI (− 0.4; − 0.11); p = 6.3e−04) and indole-3-lactate (effect estimate = − 0.26, 95% CI (− 0.38; − 0.13); p = 5.4e−05). Our study supports that lower eGFR causes higher blood metabolite levels of the tryptophan pathway including kynurenine, C-glycosyltryptophan, 3-indoxyl sulfate, and indole-3-lactate. These findings aid the notion that metabolites of the tryptophan pathway are a consequence rather than a cause of reduced eGFR. Further research is needed to specifically examine relationships with respect to chronic kidney disease (CKD) progression among patients with existing CKD.

Profiling of Mitochondrial DNA Heteroplasmy in a Prospective Oral Squamous Cell Carcinoma Study

While a shift in energy metabolism is essential to cancers, the knowledge about the involvement of the mitochondrial genome in tumorigenesis and progression in oral squamous cell carcinoma (OSCC) is still very limited. In this study, we evaluated 37 OSCC tumors and the corresponding benign mucosa tissue pairs by deep sequencing of the complete mitochondrial DNA (mtDNA). After extensive quality control, we identified 287 variants, 137 in tumor and 150 in benign samples exceeding the 1% threshold. Variant heteroplasmy levels were significantly increased in cancer compared to benign tissues (p = 0.0002). Furthermore, pairwise high heteroplasmy frequency difference variants (∆HF% > 20) with potential functional impact were increased in the cancer tissues (p = 0.024). Fourteen mutations were identified in the protein-coding region, out of which thirteen were detected in cancer and only one in benign tissue. After eight years of follow-up, the risk of mortality was higher for patients who harbored at least one ∆HF% > 20 variant in mtDNA protein-coding regions relative to those with no mutations (HR = 4.6, (95%CI = 1.3-17); p = 0.019 in primary tumor carriers). Haplogroup affiliation showed an impact on survival time, which however needs confirmation in a larger study. In conclusion, we observed a significantly higher accumulation of somatic mutations in the cancer tissues associated with a worse prognosis.

Estimation of the Required Lipoprotein(a)-Lowering Therapeutic Effect Size for Reduction in Coronary Heart Disease Outcomes: A Mendelian Randomization Analysis

Importance

Genetic and epidemiologic data suggest that lipoprotein(a) (Lp[a]) is one of the strongest genetically determined risk factors for coronary heart disease (CHD). Specific therapies to lower Lp(a) are on the horizon, but the required reduction of Lp(a) to translate into clinically relevant lowering of CHD outcomes is a matter of debate.

Objective

To estimate the required Lp(a)-lowering effect size that may be associated with a reduction of CHD outcomes compared with the effect size of low-density lipoprotein cholesterol (LDL-C)-lowering therapies.

Design, Setting, and Participants

Genetic epidemiologic study using a mendelian randomization analysis to estimate the required Lp(a)-lowering effect size for a clinically meaningful effect on outcomes. We used the effect estimates for Lp(a) from a genome-wide association study (GWAS) and meta-analysis on Lp(a) published in 2017 of 5 different primarily population-based studies of European ancestry. All Lp(a) measurements were performed in 1 laboratory. Genetic estimates for 27 single-nucleotide polymorphisms on Lp(a) concentrations were used. Odds ratios for these 27 single-nucleotide polymorphisms associated with CHD risk were retrieved from a subsample of the CHD Exome+ consortium.

Exposures

Genetic LPA score, plasma Lp(a) concentrations, and observations of statin therapies on CHD outcomes.

Main Outcomes and Measures

Coronary heart disease.

Results

The study included 13 781 individuals from the Lp(a)-GWAS-Consortium from 5 primarily population-based studies and 20 793 CHD cases and 27 540 controls from a subsample of the CHD Exome+ consortium. Four of the studies were similar in age distribution (means between 51 and 59 years), and 1 cohort was younger; mean age, 32 years. The frequency of women was similar between 51% and 55%. We estimated that the required reduction in Lp(a) effect size would be 65.7 mg/dL (95% CI, 46.3-88.3) to reach the same potential effect on clinical outcomes that can be reached by lowering LDL-C by 38.67 mg/dL (to convert to millimoles per liter, multiply by 0.0259).

Conclusions and Relevance

This mendelian randomization analysis estimated a required Lp(a)-lowering effect size of 65.7 mg/dL to reach the same effect as a 38.67-mg/dL lowering of LDL-C. However, this estimate is determined by the observed effect estimates of single-nucleotide polymorphisms on Lp(a) concentrations and is therefore influenced by the standardization of the Lp(a) assay used. As a consequence, calculations of the required Lp(a)-lowering potential of a drug to be clinically effective might have been overestimated in the past.

Lipoprotein(a) plasma levels are not associated with incident microvascular complications in type 2 diabetes mellitus

AIMS/HYPOTHESIS

Microvascular disease in type 2 diabetes is a significant cause of end-stage renal disease, blindness and peripheral neuropathy. The strict control of known risk factors, e.g. lifestyle, hyperglycaemia, hypertension and dyslipidaemia, reduces the incidence of microvascular complications, but a residual risk remains. Lipoprotein (a) [Lp(a)] is a strong risk factor for macrovascular disease in the general population. We hypothesised that plasma Lp(a) levels and the LPA gene SNPs rs10455872 and rs3798220 are associated with the incident development of microvascular complications in type 2 diabetes.

METHODS

Analyses were performed of data from the DiaGene study, a prospective study for complications of type 2 diabetes, collected in the city of Eindhoven, the Netherlands (n = 1886 individuals with type 2 diabetes, mean follow-up time = 6.97 years). To assess the relationship between plasma Lp(a) levels and the LPA SNPs with each newly developed microvascular complication (retinopathy n = 223, nephropathy n = 246, neuropathy n = 236), Cox proportional hazards models were applied and adjusted for risk factors for microvascular complications (age, sex, mean arterial pressure, non-HDL-cholesterol, HDL-cholesterol, BMI, duration of type 2 diabetes, HbA_1c and smoking).

RESULTS

No significant associations of Lp(a) plasma levels and the LPA SNPs rs10455872 and rs3798220 with prevalent or incident microvascular complications in type 2 diabetes were found. In line with previous observations the LPA SNPs rs10455872 and rs3798220 did influence the plasma Lp(a) levels.

CONCLUSIONS/INTERPRETATION

Our data show no association between Lp(a) plasma levels and the LPA SNPs with known effect on Lp(a) plasma levels with the development of microvascular complications in type 2 diabetes. This indicates that Lp(a) does not play a major role in the development of microvascular complications. However, larger studies are needed to exclude minimal effects of Lp(a) on the development of microvascular complications.

Mitochondrial DNA copy number is associated with all-cause mortality and cardiovascular events in patients with peripheral arterial disease

BACKGROUND

Dysfunctional mitochondria have an influence on inflammation and increased oxidative stress due to an excessive production of reactive oxygen species. The mitochondrial DNA copy number (mtDNA-CN) is a potential biomarker for mitochondrial dysfunction and has been associated with various diseases. However, results were partially contrasting which might have been caused by methodological difficulties to quantify mtDNA-CN.

OBJECTIVE

We aimed to investigate whether mtDNA-CN is associated with peripheral arterial disease (PAD) as well as all-cause mortality and cardiovascular events during seven years of follow-up.

METHODS

A total of 236 male patients with PAD from the Cardiovascular Disease in Intermittent Claudication (CAVASIC) study were compared with 249 age- and diabetes-matched controls. MtDNA-CN was measured with a well-standardized plasmid-normalized quantitative PCR-based assay determining the ratio between mtDNA-CN and nuclear DNA.

RESULTS

Individuals in the lowest quartile of mtDNA-CN had a twofold increased risk for PAD which, however, was no longer significant after adjusting for leukocytes and platelets. About 67 of the 236 patients had already experienced a cardiovascular event at baseline and those in the lowest mtDNA-CN quartile had a 2.34-fold increased risk for these events (95% CI 1.08-5.13). During follow-up, 37 PAD patients died and 66 patients experienced a cardiovascular event. Patients in the lowest mtDNA-CN quartile had hazard ratios of 2.66 (95% CI 1.27-5.58) for all-cause-mortality and 1.82 (95% CI 1.02-3.27) for cardiovascular events compared with the combined quartile 2-4 (adjusted for age, smoking, CRP, diabetes, prevalent cardiovascular disease, leukocytes and platelets).

CONCLUSION

This investigation supports the hypothesis of mitochondrial dysfunction in peripheral arterial disease and shows an association of low mtDNA-CNs with all-cause-mortality and prevalent and incident cardiovascular disease in PAD patients with intermittent claudication.