Blood transcriptomic biomarkers of alcohol consumption and cardiovascular disease risk factors: the Framingham Heart Study.

BACKGROUND: The relations of alcohol consumption and gene expression remain to be elucidated. MATERIALS AND METHODS: We examined cross-sectional associations between alcohol consumption and whole blood derived gene expression levels and between alcohol-associated genes and obesity, hypertension, and diabetes in 5531 Framingham Heart Study (FHS) participants. RESULTS: We identified 25 alcohol-associated genes. We further showed cross-sectional associations of 16 alcohol-associated genes with obesity, nine genes with hypertension, and eight genes with diabetes at P < 0.002. For example, we observed decreased expression of PROK2 (ß = -0.0018; 95%CI: -0.0021, -0.0007; P = 6.5e - 5) and PAX5 (ß = -0.0014; 95%CI: -0.0021, -0.0007; P = 6.5e - 5) per 1 g/day increase in alcohol consumption. Consistent with our previous observation on the inverse association of alcohol consumption with obesity and positive association of alcohol consumption with hypertension, we found that PROK2 was positively associated with obesity (OR = 1.42; 95%CI: 1.17, 1.72; P = 4.5e - 4) and PAX5 was negatively associated with hypertension (OR = 0.73; 95%CI: 0.59, 0.89; P = 1.6e - 3). We also observed that alcohol consumption was positively associated with expression of ABCA13 (ß = 0.0012; 95%CI: 0.0007, 0.0017; P = 1.3e - 6) and ABCA13 was positively associated with diabetes (OR = 2.57; 95%CI: 1.73, 3.84; P = 3.5e - 06); this finding, however, was inconsistent with our observation of an inverse association between alcohol consumption and diabetes. CONCLUSIONS: We showed strong cross-sectional associations between alcohol consumption and expression levels of 25 genes in FHS participants. Nonetheless, complex relationships exist between alcohol-associated genes and CVD risk factors.

A naturally occurring variant of SHLP2 is a protective factor in Parkinson's disease.

Mitochondrial DNA single nucleotide polymorphisms (mtSNPs) have been associated with a reduced risk of developing Parkinson's disease (PD), yet the underlying mechanisms remain elusive. In this study, we investigate the functional role of a PD-associated mtSNP that impacts the mitochondrial-derived peptide (MDP) Small Humanin-like Peptide 2 (SHLP2). We identify m.2158 T > C, a mtSNP associated with reduced PD risk, within the small open reading frame encoding SHLP2. This mtSNP results in an alternative form of SHLP2 (lysine 4 replaced with arginine; K4R). Using targeted mass spectrometry, we detect specific tryptic fragments of SHLP2 in neuronal cells and demonstrate its binding to mitochondrial complex 1. Notably, we observe that the K4R variant, associated with reduced PD risk, exhibits increased stability compared to WT SHLP2. Additionally, both WT and K4R SHLP2 show enhanced protection against mitochondrial dysfunction in in vitro experiments and confer protection against a PD-inducing toxin, a mitochondrial complex 1 inhibitor, in a mouse model. This study sheds light on the functional consequences of the m.2158 T > C mtSNP on SHLP2 and provides insights into the potential mechanisms by which this mtSNP may reduce the risk of PD.

Epigenome-wide association study of mitochondrial genome copy number.

We conducted cohort- and race-specific epigenome-wide association analyses of mitochondrial deoxyribonucleic acid (mtDNA) copy number (mtDNA CN) measured in whole blood from participants of African and European origins in five cohorts (n = 6182, mean age = 57-67 years, 65% women). In the meta-analysis of all the participants, we discovered 21 mtDNA CN-associated DNA methylation sites (CpG) (P < 1 × 10-7), with a 0.7-3.0 standard deviation increase (3 CpGs) or decrease (18 CpGs) in mtDNA CN corresponding to a 1% increase in DNA methylation. Several significant CpGs have been reported to be associated with at least two risk factors (e.g. chronological age or smoking) for cardiovascular disease (CVD). Five genes [PR/SET domain 16, nuclear receptor subfamily 1 group H member 3 (NR1H3), DNA repair protein, DNA polymerase kappa and decaprenyl-diphosphate synthase subunit 2], which harbor nine significant CpGs, are known to be involved in mitochondrial biosynthesis and functions. For example, NR1H3 encodes a transcription factor that is differentially expressed during an adipose tissue transition. The methylation level of cg09548275 in NR1H3 was negatively associated with mtDNA CN (effect size = -1.71, P = 4 × 10-8) and was positively associated with the NR1H3 expression level (effect size = 0.43, P = 0.0003), which indicates that the methylation level in NR1H3 may underlie the relationship between mtDNA CN, the NR1H3 transcription factor and energy expenditure. In summary, the study results suggest that mtDNA CN variation in whole blood is associated with DNA methylation levels in genes that are involved in a wide range of mitochondrial activities. These findings will help reveal molecular mechanisms between mtDNA CN and CVD.

Associations of Alcohol Consumption with Cardiovascular Disease-Related Proteomic Biomarkers: The Framingham Heart Study.

BACKGROUND: Alcohol consumption and cardiovascular disease (CVD) have a complex relation. OBJECTIVES: We examined the associations between alcohol consumption, fasting plasma proteins, and CVD risk. METHODS: We performed cross-sectional association analyses of alcohol consumption with 71 CVD-related plasma proteins, and also performed prospective association analyses of alcohol consumption and protein concentrations with 3 CVD risk factors (obesity, hypertension, and diabetes) in 6745 Framingham Heart Study (FHS) participants (mean age 49 y; 53% women). RESULTS: A unit increase in log10 transformed alcohol consumption (g/d) was associated with an increased risk of hypertension (HR = 1.14; 95% CI: 1.04, 1.26; P = 0.007), and decreased risks of obesity (HR = 0.80; 95% CI: 0.71, 0.91; P = 4.6 × 10-4) and diabetes (HR: 0.68; 95% CI: 0.58, 0.80; P = 5.1 × 10-6) in a median of 13-y (interquartile = 7, 14) of follow-up. We identified 43 alcohol-associated proteins in a discovery sample (n = 4348, false discovery rate <0.05) and 20 of them were significant (P <0.05/43) in an independent validation sample (n = 2397). Eighteen of the 20 proteins were inversely associated with alcohol consumption. Four of the 20 proteins demonstrated 3-way associations, as expected, with alcohol consumption and CVD risk factors. For example, a greater concentration of APOA1 was associated with higher alcohol consumption (P = 1.2 × 10-65), and it was also associated with a lower risk of diabetes (P = 8.5 × 10-6). However, several others showed unexpected 3-way associations. CONCLUSIONS: We identified 20 alcohol-associated proteins in 6745 FHS samples. These alcohol-associated proteins demonstrated complex relations with the 3 CVD risk factors. Future studies with integration of more proteomic markers and larger sample size are warranted to unravel the complex relation between alcohol consumption and CVD risk.

Association analysis of mitochondrial DNA heteroplasmic variants: methods and application.

We rigorously assessed a comprehensive association testing framework for heteroplasmy, employing both simulated and real-world data. This framework employed a variant allele fraction (VAF) threshold and harnessed multiple gene-based tests for robust identification and association testing of heteroplasmy. Our simulation studies demonstrated that gene-based tests maintained an appropriate type I error rate at α=0.001. Notably, when 5% or more heteroplasmic variants within a target region were linked to an outcome, burden-extension tests (including the adaptive burden test, variable threshold burden test, and z-score weighting burden test) outperformed the sequence kernel association test (SKAT) and the original burden test. Applying this framework, we conducted association analyses on whole-blood derived heteroplasmy in 17,507 individuals of African and European ancestries (31% of African Ancestry, mean age of 62, with 58% women) with whole genome sequencing data. We performed both cohort- and ancestry-specific association analyses, followed by meta-analysis on both pooled samples and within each ancestry group. Our results suggest that mtDNA-encoded genes/regions are likely to exhibit varying rates in somatic aging, with the notably strong associations observed between heteroplasmy in the RNR1 and RNR2 genes (p<0.001) and advance aging by the Original Burden test. In contrast, SKAT identified significant associations (p<0.001) between diabetes and the aggregated effects of heteroplasmy in several protein-coding genes. Further research is warranted to validate these findings. In summary, our proposed statistical framework represents a valuable tool for facilitating association testing of heteroplasmy with disease traits in large human populations.

Association Between Whole Blood-Derived Mitochondrial DNA Copy Number, Low-Density Lipoprotein Cholesterol, and Cardiovascular Disease Risk.

Background The relationship between mitochondrial DNA copy number (mtDNA CN) and cardiovascular disease remains elusive. Methods and Results We performed cross-sectional and prospective association analyses of blood-derived mtDNA CN and cardiovascular disease outcomes in 27 316 participants in 8 cohorts of multiple racial and ethnic groups with whole-genome sequencing. We also performed Mendelian randomization to explore causal relationships of mtDNA CN with coronary heart disease (CHD) and cardiometabolic risk factors (obesity, diabetes, hypertension, and hyperlipidemia). P<0.01 was used for significance. We validated most of the previously reported associations between mtDNA CN and cardiovascular disease outcomes. For example, 1-SD unit lower level of mtDNA CN was associated with 1.08 (95% CI, 1.04-1.12; P<0.001) times the hazard for developing incident CHD, adjusting for covariates. Mendelian randomization analyses showed no causal effect from a lower level of mtDNA CN to a higher CHD risk (ß=0.091; P=0.11) or in the reverse direction (ß=-0.012; P=0.076). Additional bidirectional Mendelian randomization analyses revealed that low-density lipoprotein cholesterol had a causal effect on mtDNA CN (ß=-0.084; P<0.001), but the reverse direction was not significant (P=0.059). No causal associations were observed between mtDNA CN and obesity, diabetes, and hypertension, in either direction. Multivariable Mendelian randomization analyses showed no causal effect of CHD on mtDNA CN, controlling for low-density lipoprotein cholesterol level (P=0.52), whereas there was a strong direct causal effect of higher low-density lipoprotein cholesterol on lower mtDNA CN, adjusting for CHD status (ß=-0.092; P<0.001). Conclusions Our findings indicate that high low-density lipoprotein cholesterol may underlie the complex relationships between mtDNA CN and vascular atherosclerosis.

Comparison of mitochondrial DNA sequences from whole blood and lymphoblastoid cell lines.

Lymphoblastoid cell lines (LCLs) provide an unlimited source of genomic DNA for genetic studies. Here, we compared mtDNA sequence variants, heteroplasmic or homplasmic, between LCL (sequenced by mitoRCA-seq method) and whole blood samples (sequenced through whole genome sequencing approach) of the same 130 participants in the Framingham Heart Study. We applied harmonization of sequence coverages and consistent quality control to mtDNA sequences. We identified 866 variation sites in the 130 LCL samples and 666 sites in the 130 blood samples. More than 94% of the identified homoplasmies were present in both LCL and blood samples while more than 70% of heteroplasmic sites were uniquely present either in LCL or in blood samples. The LCL and whole blood samples carried a similar number of homoplasmic variants (p = 0.45) per sample while the LCL carried a greater number of heteroplasmic variants than whole blood per sample (p < 2.2e-16). Furthermore, the LCL samples tended to accumulate low level heteroplasmies (heteroplasmy level in 3-25%) than their paired blood samples (p = 0.001). These results suggest that cautions should be taken in the interpretation and comparison of findings when different tissues/cell types or different sequencing technologies are applied to obtain mtDNA sequences.

Association of Serum Low-Density Lipoprotein, High-Density Lipoprotein, and Total Cholesterol With Development of Knee Osteoarthritis.

OBJECTIVE: Studies suggest an association between elevated total serum cholesterol, particularly low-density lipoprotein (LDL), and osteoarthritis (OA). The present study was undertaken to evaluate the association between total cholesterol, LDL, and high-density lipoprotein (HDL) and risk of knee OA. METHODS: We studied participants from the Multicenter Osteoarthritis study (MOST) cohort at risk of developing knee OA. From baseline through 7 years, repeated knee radiographs and magnetic resonance images (MRIs) were obtained, and knee symptoms were queried. From baseline fasting blood samples, lipids and lipoproteins were analyzed using standard assays. After excluding participants with baseline OA, we defined 2 sets of patients: those developing radiographic OA, and those developing symptomatic OA (knee pain and radiographic OA). Controls did not develop these outcomes. Additionally, we examined worsening of cartilage loss and synovitis on MRI and of knee pain using the Western Ontario and McMaster Universities Osteoarthritis Index scale. We carried out logistic regression adjusting for age, sex, body mass index, education, baseline pain, and depressive symptoms, testing total cholesterol and lipoproteins as continuous measures, and we performed sensitivity analyses examining whether commonly used thresholds for high cholesterol, LDL, or low HDL increased risk. RESULTS: We studied 337 patients with incident symptomatic OA and 283 patients with incident radiographic OA. The mean age at baseline was 62 years (55% women). Neither total cholesterol, LDL, nor HDL showed a significant association with radiographic or symptomatic OA. Additionally, we found no association of these lipid measures with cartilage loss, worsening synovitis, or worsening knee pain. CONCLUSION: Our data do not support an association between total cholesterol, LDL, or HDL with OA outcomes.

Whole Genome DNA and RNA Sequencing of Whole Blood Elucidates the Genetic Architecture of Gene Expression Underlying a Wide Range of Diseases.

To create a scientific resource of expression quantitative trail loci (eQTL), we conducted a genome-wide association study (GWAS) using genotypes obtained from whole genome sequencing (WGS) of DNA and gene expression levels from RNA sequencing (RNA-seq) of whole blood in 2622 participants in Framingham Heart Study. We identified 6,778,286 cis -eQTL variant-gene transcript (eGene) pairs at p <5×10 -8 (2,855,111 unique cis -eQTL variants and 15,982 unique eGenes) and 1,469,754 trans -eQTL variant-eGene pairs at p <1e-12 (526,056 unique trans -eQTL variants and 7,233 unique eGenes). In addition, 442,379 cis -eQTL variants were associated with expression of 1518 long non-protein coding RNAs (lncRNAs). Gene Ontology (GO) analyses revealed that the top GO terms for cis- eGenes are enriched for immune functions (FDR <0.05). The cis -eQTL variants are enriched for SNPs reported to be associated with 815 traits in prior GWAS, including cardiovascular disease risk factors. As proof of concept, we used this eQTL resource in conjunction with genetic variants from public GWAS databases in causal inference testing (e.g., COVID-19 severity). After Bonferroni correction, Mendelian randomization analyses identified putative causal associations of 60 eGenes with systolic blood pressure, 13 genes with coronary artery disease, and seven genes with COVID-19 severity. This study created a comprehensive eQTL resource via BioData Catalyst that will be made available to the scientific community. This will advance understanding of the genetic architecture of gene expression underlying a wide range of diseases.

Whole Genome DNA and RNA Sequencing of Whole Blood Elucidates the Genetic Architecture of Gene Expression Underlying a Wide Range of Diseases.

To create a scientific resource of expression quantitative trail loci (eQTL), we conducted a genome-wide association study (GWAS) using genotypes obtained from whole genome sequencing (WGS) of DNA and gene expression levels from RNA sequencing (RNA-seq) of whole blood in 2622 participants in Framingham Heart Study. We identified 6,778,286 cis -eQTL variant-gene transcript (eGene) pairs at p < 5x10 - 8 (2,855,111 unique cis -eQTL variants and 15,982 unique eGenes) and 1,469,754 trans -eQTL variant-eGene pairs at p < 1e-12 (526,056 unique trans -eQTL variants and 7,233 unique eGenes). In addition, 442,379 cis -eQTL variants were associated with expression of 1518 long non-protein coding RNAs (lncRNAs). Gene Ontology (GO) analyses revealed that the top GO terms for cis- eGenes are enriched for immune functions (FDR < 0.05). The cis -eQTL variants are enriched for SNPs reported to be associated with 815 traits in prior GWAS, including cardiovascular disease risk factors. As proof of concept, we used this eQTL resource in conjunction with genetic variants from public GWAS databases in causal inference testing (e.g., COVID-19 severity). After Bonferroni correction, Mendelian randomization analyses identified putative causal associations of 60 eGenes with systolic blood pressure, 13 genes with coronary artery disease, and seven genes with COVID-19 severity. This study created a comprehensive eQTL resource via BioData Catalyst that will be made available to the scientific community. This will advance understanding of the genetic architecture of gene expression underlying a wide range of diseases.

Whole genome DNA and RNA sequencing of whole blood elucidates the genetic architecture of gene expression underlying a wide range of diseases.

To create a scientific resource of expression quantitative trail loci (eQTL), we conducted a genome-wide association study (GWAS) using genotypes obtained from whole genome sequencing (WGS) of DNA and gene expression levels from RNA sequencing (RNA-seq) of whole blood in 2622 participants in Framingham Heart Study. We identified 6,778,286 cis-eQTL variant-gene transcript (eGene) pairs at p < 5 × 10-8 (2,855,111 unique cis-eQTL variants and 15,982 unique eGenes) and 1,469,754 trans-eQTL variant-eGene pairs at p < 1e-12 (526,056 unique trans-eQTL variants and 7233 unique eGenes). In addition, 442,379 cis-eQTL variants were associated with expression of 1518 long non-protein coding RNAs (lncRNAs). Gene Ontology (GO) analyses revealed that the top GO terms for cis-eGenes are enriched for immune functions (FDR < 0.05). The cis-eQTL variants are enriched for SNPs reported to be associated with 815 traits in prior GWAS, including cardiovascular disease risk factors. As proof of concept, we used this eQTL resource in conjunction with genetic variants from public GWAS databases in causal inference testing (e.g., COVID-19 severity). After Bonferroni correction, Mendelian randomization analyses identified putative causal associations of 60 eGenes with systolic blood pressure, 13 genes with coronary artery disease, and seven genes with COVID-19 severity. This study created a comprehensive eQTL resource via BioData Catalyst that will be made available to the scientific community. This will advance understanding of the genetic architecture of gene expression underlying a wide range of diseases.

Phenylalanine Is a Novel Marker for Radiographic Knee Osteoarthritis Progression: The MOST Study.

OBJECTIVE: To identify plasma markers associated with an increased risk of radiographic knee osteoarthritis(OA) progression using a metabolomics approach. METHODS: Study participants were from the Multicenter Osteoarthritis Study (MOST) and were categorized into 2 groups based on the presence of baseline radiographic OA. Subjects in group 1 had unilateral knee OA and subjects in group 2 had bilateral knee OA. Progression was defined as a half-grade or greater worsening in joint space width at 30-month follow-up. For group 1, a participant progressed when their OA knee showed radiographic progression and the contralateral knee developed OA; for group 2, a participant progressed when both knees with OA showed radiographic progression. Metabolomic profiling was performed on plasma samples collected at baseline and logistic regression was performed to test the association between each metabolite and knee OA progression after adjustment for age, sex, BMI, and clinic site. Significance was defined as P ≤ 0.0003 in the combined analysis. RESULTS: There were 234 progressors (57 in group 1 and 177 in group 2) and 322 nonprogressors (206 in group 1 and 116 in group 2) included in the analyses. Among 157 metabolites studied, we found that odds of progression were 1.46 times higher per SD increase of phenylalanine level (95% CI 1.20-1.77, P = 0.0001) in the combined analysis. Sex-specific analysis showed that an association was seen in women (P = 0.0002) but not in men. CONCLUSION: Our data suggest that phenylalanine might be a novel plasma marker for higher risk of bilateral radiographic knee OA progression in women.

Presence and transmission of mitochondrial heteroplasmic mutations in human populations of European and African ancestry.

We investigated the concordance of mitochondrial DNA heteroplasmic mutations (heteroplasmies) in 6745 maternal pairs of European (EA, n = 4718 pairs) and African (AA, n = 2027 pairs) Americans in whole blood. Mother-offspring pairs displayed the highest concordance rate, followed by sibling-sibling and more distantly-related maternal pairs. The allele fractions of concordant heteroplasmies exhibited high correlation (R2 = 0.8) between paired individuals. Discordant heteroplasmies were more likely to be in coding regions, be nonsynonymous or nonsynonymous-deleterious (p < 0.001). The number of deleterious heteroplasmies was significantly correlated with advancing age (20-44, 45-64, and ≥65 years, p-trend = 0.01). One standard deviation increase in heteroplasmic burden (i.e., the number of heteroplasmies carried by an individual) was associated with 0.17 to 0.26 (p < 1e - 23) standard deviation decrease in mtDNA copy number, independent of age. White blood cell count and differential count jointly explained 0.5% to 1.3% (p ≤ 0.001) variance in heteroplasmic burden. A genome-wide association and meta-analysis identified a region at 11p11.12 (top signal rs779031139, p = 2.0e - 18, minor allele frequency = 0.38) associated with the heteroplasmic burden. However, the 11p11.12 region is adjacent to a nuclear mitochondrial DNA (NUMT) corresponding to a 542 bp area of the D-loop. This region was no longer significant after excluding heteroplasmies within the 542 bp from the heteroplasmic burden. The discovery that blood mtDNA heteroplasmies were both inherited and somatic origins and that an increase in heteroplasmic burden was strongly associated with a decrease in average number of mtDNA copy number in blood are important findings to be considered in association studies of mtDNA with disease traits.

Comparison of Ultrasound Features of Major Salivary Glands in Sarcoidosis, Amyloidosis, and Sjögren's Syndrome.

OBJECTIVE: While salivary gland ultrasound (SGU) has gained prominence for evaluating Sjögren's syndrome, little information exists on SGU appearance of sarcoidosis and amyloidosis, potential mimics of Sjögren's syndrome. Our goal in this study was to estimate the diagnostic accuracy of major SGU features in differentiating Sjögren's syndrome from sarcoidosis, amyloidosis, and controls. METHODS: We enrolled consecutive adult ambulatory patients with a clinical diagnosis of Sjögren's syndrome fulfilling the 2016 American College of Rheumatology (ACR) classification criteria; we also enrolled patients with a clinical diagnosis of sarcoidosis or systemic immunoglobulin light chain (AL) amyloidosis, with histologic confirmation from any tissue, and rheumatology outpatients without diagnoses affecting salivary glands. Subjects underwent major SGU using the Hocevar protocol, with resulting video clips reviewed blind to clinical diagnosis. RESULTS: Sjögren's syndrome SGU scores were greater than in patients from the other groups, but there were no distinguishing salivary gland features from AL amyloidosis or sarcoidosis. None of the patients in the control group scored higher than 17, a cutoff previously suggested for Sjögren's syndrome, but 27% of patients with AL amyloidosis and 19% with sarcoidosis scored higher than 17. Adding Hocevar SGU scores of ≥17 to the 2016 ACR/European League Against Rheumatism criteria in a parallel scheme increased the sensitivity for Sjögren's syndrome from 87% to 98%, while combining the 2 criteria in series increased specificity from 81% to 98%. CONCLUSION: Sjögren's syndrome, sarcoidosis, and AL amyloidosis share common SGU features that can help distinguish these conditions from patients without systemic rheumatologic disease. Clinicians should carefully consider these potential mimics when interpreting salivary gland US results.

Association of Visceral Adiposity With Pain but Not Structural Osteoarthritis.

OBJECTIVE: Osteoarthritis (OA) and pain are both made more severe by low-grade inflammation. This study was undertaken to examine whether visceral fat, a major source of inflammatory cytokines and adipokines, is associated with an increased risk of knee OA or musculoskeletal pain. METHODS: Subjects in the Multicenter Osteoarthritis Study cohort, who were age 50-79 years and had or were at high risk of knee OA, underwent whole-body dual x-ray absorptiometry (DXA) at baseline. At baseline, 30 months, and 60 months radiographs and magnetic resonance images (MRIs) of the knees were obtained, and patients were asked to score the severity of their knee pain and to identify sites of joint pain using a body homunculus. Baseline DXA scans were used to measure total body fat and visceral and subcutaneous fat in the torso. The association of fat depot size with structural outcomes (incident radiographic OA and cartilage loss and synovitis on MRI) and with pain outcomes (worsening knee pain, number of painful joints, and widespread pain) was assessed. Regression analyses were adjusted for age, sex, race, education level, smoking status, physical activity, body mass index (BMI), and depressive symptoms. RESULTS: Of the 2,961 participants at baseline, 60.7% were women. The mean age was 62.5 years and mean BMI was 30.5 kg/m2 . After adjustment for covariates, no fat measures were associated with any structural outcomes. However, total and visceral, but not subcutaneous, fat were positively associated with worsening knee pain (P = 0.0005 for total fat and P = 0.007 for visceral fat) and widespread pain (P = 0.001 for total fat and P = 0.02 for visceral fat), and the amount of visceral fat was associated with the number of painful joints (P = 0.07). CONCLUSION: Our findings indicate that visceral fat is associated with an increased risk of musculoskeletal and widespread pain.