DDN3.0: determining significant rewiring of biological network structure with differential dependency networks.

MOTIVATION: Complex diseases are often caused and characterized by misregulation of multiple biological pathways. Differential network analysis aims to detect significant rewiring of biological network structures under different conditions and has become an important tool for understanding the molecular etiology of disease progression and therapeutic response. With few exceptions, most existing differential network analysis tools perform differential tests on separately learned network structures that are computationally expensive and prone to collapse when grouped samples are limited or less consistent. RESULTS: We previously developed an accurate differential network analysis method-differential dependency networks (DDN), that enables joint learning of common and rewired network structures under different conditions. We now introduce the DDN3.0 tool that improves this framework with three new and highly efficient algorithms, namely, unbiased model estimation with a weighted error measure applicable to imbalance sample groups, multiple acceleration strategies to improve learning efficiency, and data-driven determination of proper hyperparameters. The comparative experimental results obtained from both realistic simulations and case studies show that DDN3.0 can help biologists more accurately identify, in a study-specific and often unknown conserved regulatory circuitry, a network of significantly rewired molecular players potentially responsible for phenotypic transitions. AVAILABILITY AND IMPLEMENTATION: The Python package of DDN3.0 is freely available at https://github.com/cbil-vt/DDN3. A user's guide and a vignette are provided at https://ddn-30.readthedocs.io/.

CAM3.0: determining cell type composition and expression from bulk tissues with fully unsupervised deconvolution.

MOTIVATION: Complex tissues are dynamic ecosystems consisting of molecularly distinct yet interacting cell types. Computational deconvolution aims to dissect bulk tissue data into cell type compositions and cell-specific expressions. With few exceptions, most existing deconvolution tools exploit supervised approaches requiring various types of references that may be unreliable or even unavailable for specific tissue microenvironments. RESULTS: We previously developed a fully unsupervised deconvolution method-Convex Analysis of Mixtures (CAM), that enables estimation of cell type composition and expression from bulk tissues. We now introduce CAM3.0 tool that improves this framework with three new and highly efficient algorithms, namely, radius-fixed clustering to identify reliable markers, linear programming to detect an initial scatter simplex, and a smart floating search for the optimum latent variable model. The comparative experimental results obtained from both realistic simulations and case studies show that the CAM3.0 tool can help biologists more accurately identify known or novel cell markers, determine cell proportions, and estimate cell-specific expressions, complementing the existing tools particularly when study- or datatype-specific references are unreliable or unavailable. AVAILABILITY AND IMPLEMENTATION: The open-source R Scripts of CAM3.0 is freely available at https://github.com/ChiungTingWu/CAM3/(https://github.com/Bioconductor/Contributions/issues/3205). A user's guide and a vignette are provided.

Adipose tissue-derived metabolite risk scores and risk for type 2 diabetes in South Asians.

BACKGROUND: South Asians are at higher risk for type 2 diabetes (T2D) than many other race/ethnic groups. Ectopic adiposity, specifically hepatic steatosis and visceral fat may partially explain this. Our objective was to derive metabolite risk scores for ectopic adiposity and assess associations with incident T2D in South Asians. METHODS: We examined 550 participants in the Mediators of Atherosclerosis in South Asians Living in America (MASALA) cohort study aged 40-84 years without known cardiovascular disease or T2D and with metabolomic data. Computed tomography scans at baseline assessed hepatic attenuation and visceral fat area, and fasting serum specimens at baseline and after 5 years assessed T2D. LC-MS-based untargeted metabolomic analysis was performed followed by targeted integration and reporting of known signals. Elastic net regularized linear regression analyses was used to derive risk scores for hepatic steatosis and visceral fat using weighted coefficients. Logistic regression models associated metabolite risk score and incident T2D, adjusting for age, gender, study site, BMI, physical activity, diet quality, energy intake and use of cholesterol-lowering medication. RESULTS: Average age of participants was 55 years, 36% women with an average body mass index (BMI) of 25 kg/m2 and 6% prevalence of hepatic steatosis, with 47 cases of incident T2D at 5 years. There were 445 metabolites of known identity. Of these, 313 metabolites were included in the MET-Visc score and 267 in the MET-Liver score. In most fully adjusted models, MET-Liver (OR 2.04 [95% CI 1.38, 3.03]) and MET-Visc (OR 2.80 [1.75, 4.46]) were associated with higher odds of T2D. These associations remained significant after adjustment for measured adiposity. CONCLUSIONS: Metabolite risk scores for intrahepatic fat and visceral fat were strongly related to incident T2D independent of measured adiposity. Use of these biomarkers to target risk stratification may help capture pre-clinical metabolic abnormalities.

Metabolomic Profiling of Cholesterol Efflux Capacity in a Multiethnic Population: Insights From MESA.

BACKGROUND: Impaired cholesterol efflux capacity (CEC) is a novel lipid metabolism trait associated with atherosclerotic cardiovascular disease. Mechanisms underlying CEC variation are unknown. We evaluated associations of circulating metabolites with CEC to advance understanding of metabolic pathways involved in cholesterol efflux regulation. METHODS: Participants enrolled in the MESA (Multi-Ethnic Study of Atherosclerosis) who underwent nuclear magnetic resonance metabolome profiling and CEC measurement (N=3543) at baseline were included. Metabolite associations with CEC were evaluated using standard linear regression analyses. Repeated ElasticNet and multilayer perceptron regression were used to assess metabolite profile predictive performance for CEC. Features important for CEC prediction were identified using Shapley Additive Explanations values. RESULTS: Greater CEC was significantly associated with metabolite clusters composed of the largest-sized particle subclasses of VLDL (very-low-density lipoprotein) and HDL (high-density lipoprotein), as well as their constituent apo A1, apo A2, phospholipid, and cholesterol components (ß=0.072-0.081; P<0.001). Metabolite profiles had poor accuracy for predicting in vitro CEC in linear and nonlinear analyses (R2<0.02; Spearman ρ<0.18). The most important feature for CEC prediction was race, with Black participants having significantly lower CEC compared with other races. CONCLUSIONS: We identified independent associations among CEC, the largest-sized particle subclasses of VLDL and HDL, and their constituent apolipoproteins and lipids. A large proportion of variation in CEC remained unexplained by metabolites and traditional clinical risk factors, supporting further investigation into genomic, proteomic, and phospholipidomic determinants of CEC.

AtheroSpectrum Reveals Novel Macrophage Foam Cell Gene Signatures Associated With Atherosclerotic Cardiovascular Disease Risk.

BACKGROUND: Whereas several interventions can effectively lower lipid levels in people at risk for atherosclerotic cardiovascular disease (ASCVD), cardiovascular event risks remain, suggesting an unmet medical need to identify factors contributing to cardiovascular event risk. Monocytes and macrophages play central roles in atherosclerosis, but studies have yet to provide a detailed view of macrophage populations involved in increased ASCVD risk. METHODS: A novel macrophage foaming analytics tool, AtheroSpectrum, was developed using 2 quantitative indices depicting lipid metabolism and the inflammatory status of macrophages. A machine learning algorithm was developed to analyze gene expression patterns in the peripheral monocyte transcriptome of MESA participants (Multi-Ethnic Study of Atherosclerosis; set 1; n=911). A list of 30 genes was generated and integrated with traditional risk factors to create an ASCVD risk prediction model (30-gene cardiovascular disease risk score [CR-30]), which was subsequently validated in the remaining MESA participants (set 2; n=228); performance of CR-30 was also tested in 2 independent human atherosclerotic tissue transcriptome data sets (GTEx [Genotype-Tissue Expression] and GSE43292). RESULTS: Using single-cell transcriptomic profiles (GSE97310, GSE116240, GSE97941, and FR-FCM-Z23S), AtheroSpectrum detected 2 distinct programs in plaque macrophages-homeostatic foaming and inflammatory pathogenic foaming-the latter of which was positively associated with severity of atherosclerosis in multiple studies. A pool of 2209 pathogenic foaming genes was extracted and screened to select a subset of 30 genes correlated with cardiovascular event in MESA set 1. A cardiovascular disease risk score model (CR-30) was then developed by incorporating this gene set with traditional variables sensitive to cardiovascular event in MESA set 1 after cross-validation generalizability analysis. The performance of CR-30 was then tested in MESA set 2 (P=2.60×10-4; area under the receiver operating characteristic curve, 0.742) and 2 independent data sets (GTEx: P=7.32×10-17; area under the receiver operating characteristic curve, 0.664; GSE43292: P=7.04×10-2; area under the receiver operating characteristic curve, 0.633). Model sensitivity tests confirmed the contribution of the 30-gene panel to the prediction model (likelihood ratio test; df=31, P=0.03). CONCLUSIONS: Our novel computational program (AtheroSpectrum) identified a specific gene expression profile associated with inflammatory macrophage foam cells. A subset of 30 genes expressed in circulating monocytes jointly contributed to prediction of symptomatic atherosclerotic vascular disease. Incorporating a pathogenic foaming gene set with known risk factors can significantly strengthen the power to predict ASCVD risk. Our programs may facilitate both mechanistic investigations and development of therapeutic and prognostic strategies for ASCVD risk.

COVID-19 Symptoms by Variant Period in the North Carolina COVID-19 Community Research Partnership, North Carolina, USA.

In North Carolina, USA, the SARS-CoV-2 Omicron variant was associated with changing symptomology in daily surveys, including increasing rates of self-reported cough and sore throat and decreased rates of loss of taste and smell. Compared with the pre-Delta period, Delta and Omicron (pre-BA.4/BA.5) variant periods were associated with shorter symptom duration.

swCAM: estimation of subtype-specific expressions in individual samples with unsupervised sample-wise deconvolution.

MOTIVATION: Complex biological tissues are often a heterogeneous mixture of several molecularly distinct cell subtypes. Both subtype compositions and subtype-specific (STS) expressions can vary across biological conditions. Computational deconvolution aims to dissect patterns of bulk tissue data into subtype compositions and STS expressions. Existing deconvolution methods can only estimate averaged STS expressions in a population, while many downstream analyses such as inferring co-expression networks in particular subtypes require subtype expression estimates in individual samples. However, individual-level deconvolution is a mathematically underdetermined problem because there are more variables than observations. RESULTS: We report a sample-wise Convex Analysis of Mixtures (swCAM) method that can estimate subtype proportions and STS expressions in individual samples from bulk tissue transcriptomes. We extend our previous CAM framework to include a new term accounting for between-sample variations and formulate swCAM as a nuclear-norm and ℓ2,1-norm regularized matrix factorization problem. We determine hyperparameter values using cross-validation with random entry exclusion and obtain a swCAM solution using an efficient alternating direction method of multipliers. Experimental results on realistic simulation data show that swCAM can accurately estimate STS expressions in individual samples and successfully extract co-expression networks in particular subtypes that are otherwise unobtainable using bulk data. In two real-world applications, swCAM analysis of bulk RNASeq data from brain tissue of cases and controls with bipolar disorder or Alzheimer's disease identified significant changes in cell proportion, expression pattern and co-expression module in patient neurons. Comparative evaluation of swCAM versus peer methods is also provided. AVAILABILITY AND IMPLEMENTATION: The R Scripts of swCAM are freely available at https://github.com/Lululuella/swCAM. A user's guide and a vignette are provided. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Metabolomic Profile of the Healthy Eating Index-2015 in the Multiethnic Study of Atherosclerosis.

BACKGROUND: Poor diet quality is a risk factor for type 2 diabetes and cardiovascular disease. However, knowledge of metabolites marking adherence to Dietary Guidelines for Americans (2015 version) are limited. OBJECTIVES: The goal was to determine a pattern of metabolites associated with the Healthy Eating Index (HEI)-2015, which measures adherence to the Dietary Guidelines for Americans. METHODS: The analysis examined 3557 adult men and women from the longitudinal cohort Multiethnic Study of Atherosclerosis (MESA), without known cardiovascular disease and with complete dietary data. Fasting serum specimens and diet and demographic questionnaires were assessed at baseline. Untargeted 1H 1-dimensional nuclei magnetic resonance spectroscopy (600 MHz) was used to generate metabolomics and lipidomics. A metabolome-wide association study specified each spectral feature as outcomes, HEI-2015 score as predictor, adjusting for age, sex, race, and study site in linear regression analyses. Subsequently, hierarchical clustering defined the discrete groups of correlated nuclei magnetic resonance features associated with named metabolites, and the linear regression analysis assessed for associations with HEI-2015 total and component scores. RESULTS: The sample included 50% women with an mean age of 63 years, with 40% identifying as White, 23% as Black, 24% as Hispanic, and 13% as Chinese American. The mean HEI-2015 score was 66. The metabolome-wide association study identified 179 spectral features significantly associated with HEI-2015 score. The cluster analysis identified 7 clusters representing 4 metabolites; HEI-2015 score was significantly associated with all. HEI-2015 score was associated with proline betaine [ß = 0.12 (SE = 0.02); P = 4.70 × 10-13] and was inversely related to proline [ß = -0.13 (SE = 0.02); P = 4.45 × 10-14], 1,5 anhydrosorbitol [ß = -0.08 (SE = 0.02); P = 4.37 × 10-7] and unsaturated fatty acyl chains [ß = 0.08 (SE = 0.02); P = 8.98 × 10-7]. Intake of total fruit, whole grains, and seafood and plant proteins was associated with proline betaine. CONCLUSIONS: Diet quality is significantly associated with unsaturated fatty acyl chains, proline betaine, and proline. Further analysis may clarify the link between diet quality, metabolites, and pathogenesis of cardiometabolic disease.

Associations between Metabolomic Biomarkers of Avocado Intake and Glycemia in the Multi-Ethnic Study of Atherosclerosis.

BACKGROUND: Avocado consumption is linked to better glucose homeostasis, but small associations suggest potential population heterogeneity. Metabolomic data capture the effects of food intake after digestion and metabolism, thus accounting for individual differences in these processes. OBJECTIVES: To identify metabolomic biomarkers of avocado intake and to examine their associations with glycemia. METHODS: Baseline data from 6224 multi-ethnic older adults (62% female) included self-reported avocado intake, fasting glucose and insulin, and untargeted plasma proton nuclear magnetic resonance metabolomic features (metabolomic data were available for a randomly selected subset; N = 3438). Subsequently, incident type 2 diabetes (T2D) was assessed over an ∼18 y follow-up period. A metabolome-wide association study of avocado consumption status (consumer compared with nonconsumer) was conducted, and the relationship of these features with glycemia via cross-sectional associations with fasting insulin and glucose and longitudinal associations with incident T2D was examined. RESULTS: Three highly-correlated spectral features were associated with avocado intake at metabolome-wide significance levels (P < 5.3 ∗ 10-7) and combined into a single biomarker. We did not find evidence that these features were additionally associated with overall dietary quality, nor with any of 47 other food groups (all P > 0.001), supporting their suitability as a biomarker of avocado intake. Avocado intake showed a modest association only with lower fasting insulin (ß = -0.07 +/- 0.03, P = 0.03), an association that was attenuated to nonsignificance when additionally controlling for body mass index (kg/m2). However, our biomarker of avocado intake was strongly associated with lower fasting glucose (ß = -0.22 +/- 0.02, P < 2.0 ∗ 10-16), lower fasting insulin (ß = -0.17 +/- 0.02, P < 2.0 ∗ 10-16), and a lower incidence of T2D (hazard ratio: 0.68; 0.63-074, P < 2.0 ∗ 10-16), even when adjusting for BMI. CONCLUSIONS: Highly significant associations between glycemia and avocado-related metabolomic features, which serve as biomarkers of the physiological impact of dietary intake after digestion and absorption, compared to modest relationships between glycemia and avocado consumption, highlights the importance of considering individual differences in metabolism when considering diet-health relationships.

CLINICAL HETEROGENEITY IN THE AGE OF BIG DATA, ADVANCED ANALYTICS, AND COMPLEXITY THEORY.

Clinical heterogeneity remains a challenge in the practice of medicine and is an underlying motivation for much of biomedical research. Unfortunately, despite an abundance of technologies capable of producing millions of discrete data elements with information about a patient's health status or disease prognosis, our ability to translate those data into meaningful improvements in understanding of clinical heterogeneity is limited. To address this gap, we have applied newer approaches to manifold learning and developed additional and complementary techniques to interrogate and interpret complex, high dimensional omics data. The central premise is that there exist manifolds embedded in high dimensional data that represent fundamental biologic processes that may help address the challenges of clinical heterogeneity. Preliminary evidence from several real-world data sets suggests that these techniques can identify coherent and reproducible manifolds embedded in higher dimensional omics data. Work is currently ongoing to determine the clinical informativeness of these novel data structures.

Eicosapentaenoic Acid Ameliorates Cardiac Fibrosis and Tissue Inflammation in Spontaneously Hypertensive Rats.

Hypertension affects 1 in 3 adults in the United States and leads to left ventricular (LV) concentric hypertrophy, interstitial fibrosis, and increased stiffness. The treatment of cardiac fibrosis remains challenging and empiric. Eicosapentaenoic acid (EPA) is an omega-3 polyunsaturated fatty acid that is highly effective in reducing cardiovascular events in patients and cardiac fibrosis and hypertrophy in animals when administered before pressure overload by promoting the increase of anti-inflammatory M1 macrophages. In this study, we investigated whether EPA mitigates the exacerbation of cardiac remodeling and fibrosis induced by established hypertension, a situation that closely recapitulates a clinical scenario. Twelve-week-old spontaneously hypertensive rats were randomized to eat an EPA-enriched or control diet for 20 weeks. We report that rats eating the EPA-enriched diet exhibited a reduction of interstitial cardiac fibrosis and ameliorated LV diastolic dysfunction despite the continuous increase in blood pressure. However, we found that EPA did not have an impact on cardiac hypertrophy. Interestingly, the EPA diet increased mRNA expression of M2 macrophage marker Mrc1 and interleukin-10 in cardiac tissue. These findings indicated that the antifibrotic effects of EPA are mediated in part by phenotypic polarization of macrophages toward anti-inflammatory M2 macrophages and increases of the anti-inflammatory cytokine, interleukin-10. In summary, EPA prevents the exacerbation of cardiac fibrosis and LV diastolic dysfunction during sustained pressure overload. EPA could represent a novel treatment strategy for hypertensive cardiomyopathy.

Automated Annotation of Untargeted All-Ion Fragmentation LC-MS Metabolomics Data with MetaboAnnotatoR.

Untargeted metabolomics and lipidomics LC-MS experiments produce complex datasets, usually containing tens of thousands of features from thousands of metabolites whose annotation requires additional MS/MS experiments and expert knowledge. All-ion fragmentation (AIF) LC-MS/MS acquisition provides fragmentation data at no additional experimental time cost. However, analysis of such datasets requires reconstruction of parent-fragment relationships and annotation of the resulting pseudo-MS/MS spectra. Here, we propose a novel approach for automated annotation of isotopologues, adducts, and in-source fragments from AIF LC-MS datasets by combining correlation-based parent-fragment linking with molecular fragment matching. Our workflow focuses on a subset of features rather than trying to annotate the full dataset, saving time and simplifying the process. We demonstrate the workflow in three human serum datasets containing 599 features manually annotated by experts. Precision and recall values of 82-92% and 82-85%, respectively, were obtained for features found in the highest-rank scores (1-5). These results equal or outperform those obtained using MS-DIAL software, the current state of the art for AIF data annotation. Further validation for other biological matrices and different instrument types showed variable precision (60-89%) and recall (10-88%) particularly for datasets dominated by nonlipid metabolites. The workflow is freely available as an open-source R package, MetaboAnnotatoR, together with the fragment libraries from Github (https://github.com/gggraca/MetaboAnnotatoR).

Finding Correspondence between Metabolomic Features in Untargeted Liquid Chromatography-Mass Spectrometry Metabolomics Datasets.

Integration of multiple datasets can greatly enhance bioanalytical studies, for example, by increasing power to discover and validate biomarkers. In liquid chromatography-mass spectrometry (LC-MS) metabolomics, it is especially hard to combine untargeted datasets since the majority of metabolomic features are not annotated and thus cannot be matched by chemical identity. Typically, the information available for each feature is retention time (RT), mass-to-charge ratio (m/z), and feature intensity (FI). Pairs of features from the same metabolite in separate datasets can exhibit small but significant differences, making matching very challenging. Current methods to address this issue are too simple or rely on assumptions that cannot be met in all cases. We present a method to find feature correspondence between two similar LC-MS metabolomics experiments or batches using only the features' RT, m/z, and FI. We demonstrate the method on both real and synthetic datasets, using six orthogonal validation strategies to gauge the matching quality. In our main example, 4953 features were uniquely matched, of which 585 (96.8%) of 604 manually annotated features were correct. In a second example, 2324 features could be uniquely matched, with 79 (90.8%) out of 87 annotated features correctly matched. Most of the missed annotated matches are between features that behave very differently from modeled inter-dataset shifts of RT, MZ, and FI. In a third example with simulated data with 4755 features per dataset, 99.6% of the matches were correct. Finally, the results of matching three other dataset pairs using our method are compared with a published alternative method, metabCombiner, showing the advantages of our approach. The method can be applied using M2S (Match 2 Sets), a free, open-source MATLAB toolbox, available at https://github.com/rjdossan/M2S.

Diet Patterns Are Associated with Circulating Metabolites and Lipid Profiles of South Asians in the United States.

BACKGROUND: South Asians are at higher risk for cardiometabolic disease than many other racial/ethnic minority groups. Diet patterns in US South Asians have unique components associated with cardiometabolic disease. OBJECTIVES: We aimed to characterize the metabolites associated with 3 representative diet patterns. METHODS: We included 722 participants in the Mediators of Atherosclerosis in South Asians Living in America (MASALA) cohort study aged 40-84 y without known cardiovascular disease. Fasting serum specimens and diet and demographic questionnaires were collected at baseline and diet patterns previously generated through principal components analysis. LC-MS-based untargeted metabolomic and lipidomic analysis was conducted with targeted integration of known metabolite and lipid signals. Linear regression models of diet pattern factor score and log-transformed metabolites adjusted for age, sex, caloric intake, and BMI and adjusted for multiple comparisons were performed, followed by elastic net linear regression of significant metabolites. RESULTS: There were 443 metabolites of known identity extracted from the profiling data. The "animal protein" diet pattern was associated with 61 metabolites and lipids, including glycerophospholipids phosphatidylethanolamine PE(O-16:1/20:4) and/or PE(P-16:0/20:4) (ß: 0.13; 95% CI: 0.11, 0.14) and N-acyl phosphatidylethanolamines (NAPEs) NAPE(O-18:1/20:4/18:0) and/or NAPE(P-18:0/20:4/18:0) (ß: 0.13; 95% CI: 0.11, 0.14), lysophosphatidylinositol (LPI) (22:6/0:0) (ß: 0.14; 95% CI: 0.12, 0.17), and fatty acid (FA) (22:6) (ß: 0.15; 95% CI: 0.13, 0.17). The "fried snacks, sweets, high-fat dairy" pattern was associated with 12 lipids, including PC(16:0/22:6) (ß: -0.08; 95% CI: -0.09, -0.06) and FA (22:6) (ß: 0.14; 95% CI: -0.17, -0.10). The "fruits, vegetables, nuts, and legumes" pattern was associated with 5 metabolites including proline betaine (ß: 0.17; 95% CI: 0.09, 0.25) (P < 0.0002). CONCLUSIONS: Three predominant dietary patterns in US South Asians are associated with circulating metabolites differentiated by lipids including glycerophospholipids and PUFAs and the amino acid proline betaine.

Survey of Adherence with COVID-19 Prevention Behaviors During the 2020 Thanksgiving and Winter Holidays Among Members of the COVID-19 Community Research Partnership.

Prevention behaviors represent important public health tools to limit spread of SARS-CoV-2. Adherence with recommended public health prevention behaviors among 20000 + members of a COVID-19 syndromic surveillance cohort from the mid-Atlantic and southeastern United States was assessed via electronic survey following the 2020 Thanksgiving and winter holiday (WH) seasons. Respondents were predominantly non-Hispanic Whites (90%), female (60%), and ≥ 50 years old (59%). Non-household members (NHM) were present at 47.1% of Thanksgiving gatherings and 69.3% of WH gatherings. Women were more likely than men to gather with NHM (p < 0.0001). Attending gatherings with NHM decreased with older age (Thanksgiving: 60.0% of participants aged < 30 years to 36.3% aged ≥ 70 years [p-trend < 0.0001]; WH: 81.6% of those < 30 years to 61.0% of those ≥ 70 years [p-trend < 0.0001]). Non-Hispanic Whites were more likely to gather with NHM than were Hispanics or non-Hispanic Blacks (p < 0.0001). Mask wearing, reported by 37.3% at Thanksgiving and 41.9% during the WH, was more common among older participants, non-Hispanic Blacks, and Hispanics when gatherings included NHM. In this survey, most people did not fully adhere to recommended public health safety behaviors when attending holiday gatherings. It remains unknown to what extent failure to observe these recommendations may have contributed to the COVID-19 surges observed following Thanksgiving and the winter holidays in the United States.

Monocyte chemoattractant protein-1 is not predictive of cardiac events in patients with non-low-risk chest pain.

BACKGROUND: Prior studies suggest monocyte chemoattractant protein-1 (MCP-1) may be useful for risk stratifying ED patients with chest pain. We hypothesise that MCP-1 will be predictive of 90-day major adverse cardiovascular events (MACEs) in non-low-risk patients. METHODS: A case-control study was nested within a prospective multicentre cohort (STOP-CP), which enrolled adult patients being evaluated for acute coronary syndrome at eight US EDs from 25 January 2017 to 06 September 2018. Patients with a History, ECG, Age, and Risk factor score (HEAR score) ≥4 or coronary artery disease (CAD), a non-ischaemic ECG, and non-elevated contemporary troponins at 0 and 3 hours were included. Cases were patients with 90-day MACE (all-cause death, myocardial infarction or revascularisation). Controls were patients without MACE selected with frequency matching using age, sex, race, and HEAR score or the presence of CAD. Serum MCP-1 was measured. Sensitivity and specificity were determined for cut-off points of 194 pg/mL, 200 pg/mL, 238 pg/mL and 281 pg/mL. Logistic regression adjusting for age, sex, race, and HEAR score/presence of CAD was used to determine the association between MCP-1 and 90-day MACE. A separate logistic model also included high-sensitivity troponin (hs-cTnT). RESULTS: Among 40 cases and 179 controls, there was no difference in age (p=0.90), sex (p=1.00), race (p=0.85), or HEAR score/presence of CAD (p=0.89). MCP-1 was similar in cases (median 191.9 pg/mL, IQR: 161.8-260.1) and controls (median 196.6 pg/mL, IQR: 163.0-261.1) (p=0.48). At a cut-off point of 194 pg/mL, MCP-1 was 50.0% (95% CI 33.8% to 66.2%) sensitive and 46.9% (95% CI 39.4% to 54.5%) specific for 90-day MACE. After adjusting for covariates, MCP-1 was not associated with 90-day MACE at any cut-off point (at 194 pg/mL, OR 0.88 (95% CI 0.43 to 1.79)). When including hs-cTnT in the model, MCP-1 was not associated with 90-day MACE at any cut-off point (at 194 pg/mL, OR 0.85 (95% CI 0.42 to 1.73)). CONCLUSION: MCP-1 is not predictive of 90-day MACE in patients with non-low-risk chest pain.

Sex and race differences in safety and effectiveness of the HEART pathway accelerated diagnostic protocol for acute chest pain.

BACKGROUND: The HEART Pathway is an accelerated diagnostic protocol for Emergency Department patients with possible acute coronary syndrome. The objective was to compare the safety and effectiveness of the HEART Pathway among women vs men and whites vs non-whites. METHODS: A subgroup analysis of the HEART Pathway Implementation Study was conducted. Adults with chest pain were accrued from November 2013 to January 2016 from 3 Emergency Departments in North Carolina. The primary outcomes were death and myocardial infarction (MI) and hospitalization rates at 30 days. Logistic regression evaluated for interactions of accelerated diagnostic protocol implementation with sex or race and changes in outcomes within subgroups. RESULTS: A total of 8,474 patients were accrued, of which 53.6% were female and 34.0% were non-white. The HEART Pathway identified 32.6% of females as low-risk vs 28.5% of males (P = 002) and 35.6% of non-whites as low-risk vs 28.0% of whites (P < .0001). Among low-risk patients, death or MI at 30 days occurred in 0.4% of females vs 0.5% of males (P = .70) and 0.5% of non-whites vs 0.3% of whites (P = .69). Hospitalization at 30 days was reduced by 6.6% in females (aOR: 0.74, 95% CI: 0.64-0.85), 5.1% in males (aOR: 0.87, 95% CI: 0.75-1.02), 8.6% in non-whites (aOR: 0.72, 95% CI: 0.60-0.86), and 4.5% in whites (aOR: 0.83, 95% CI: 0.73-0.94). Interactions were not significant. CONCLUSION: Women and non-whites are more likely to be classified as low-risk by the HEART Pathway. HEART Pathway implementation is associated with decreased hospitalizations and a very low death and MI rate among low-risk patients regardless of sex or race.

Chronic Ca2+/Calmodulin-Dependent Protein Kinase II Inhibition Rescues Advanced Heart Failure.

Ca2+/calmodulin-dependent protein kinase II (CaMKII) is upregulated in congestive heart failure (CHF), contributing to electrical, structural, and functional remodeling. CaMKII inhibition is known to improve CHF, but its direct cardiac effects in CHF remain unclear. We hypothesized that CaMKII inhibition improves cardiomyocyte function, [Ca2+]i regulation, and ß-adrenergic reserve, thus improving advanced CHF. In a 16-week study, we compared plasma neurohormonal levels and left ventricular (LV)- and myocyte-functional and calcium transient ([Ca2+]iT) responses in male Sprague-Dawley rats (10/group) with CHF induced by isoproterenol (170 mg/kg sq for 2 days). In rats with CHF, we studied the effects of the CaMKII inhibitor KN-93 or its inactive analog KN-92 (n = 4) (70 µg/kg per day, mini-pump) for 4 weeks. Compared with controls, isoproterenol-treated rats had severe CHF with 5-fold-increased plasma norepinephrine and about 50% decreases in ejection fraction (EF) and LV contractility [slope of LV end-systolic pressure-LV end-systolic volume relation (EES)] but increased time constant of LV relaxation (τ). They also showed significantly reduced myocyte contraction [maximum rate of myocyte shortening (dL/dtmax)], relaxation (dL/dtmax), and [Ca2+]iT Isoproterenol superfusion caused significantly fewer increases in dL/dtmax and [Ca2+]iT KN-93 treatment prevented plasma norepinephrine elevation, with increased basal and acute isoproterenol-stimulated increases in EF and EES and decreased τ in CHF. KN-93 treatment preserved normal myocyte contraction, relaxation, [Ca2+]iT, and ß-adrenergic reserve, whereas KN-92 treatment failed to improve LV and myocyte function, and plasma norepinephrine remained high in CHF. Thus, chronic CaMKII inhibition prevented CHF-induced activation of the sympathetic nervous system, restoring normal LV and cardiomyocyte basal and ß-adrenergic-stimulated contraction, relaxation, and [Ca2+]iT, thereby playing a rescue role in advanced CHF. SIGNIFICANCE STATEMENT: We investigated the therapeutic efficacy of late initiation of chronic Ca2+/calmodulin-dependent protein kinase II (CaMKII) inhibition on progression of advanced congestive heart failure (CHF). Chronic CaMKII inhibition prevented CHF-induced activation of the sympathetic nervous system and restored normal intrinsic cardiomyocyte basal and ß-adrenergic receptor-stimulated relaxation, contraction, and [Ca2+]i regulation, leading to reversal of CHF progression. These data provide new evidence that CaMKII inhibition is able and sufficient to rescue a failing heart, and thus cardiac CaMKII inhibition is a promising target for improving CHF treatment.

debCAM: a bioconductor R package for fully unsupervised deconvolution of complex tissues.

SUMMARY: We develop a fully unsupervised deconvolution method to dissect complex tissues into molecularly distinctive tissue or cell subtypes based on bulk expression profiles. We implement an R package, deconvolution by Convex Analysis of Mixtures (debCAM) that can automatically detect tissue/cell-specific markers, determine the number of constituent subtypes, calculate subtype proportions in individual samples and estimate tissue/cell-specific expression profiles. We demonstrate the performance and biomedical utility of debCAM on gene expression, methylation, proteomics and imaging data. With enhanced data preprocessing and prior knowledge incorporation, debCAM software tool will allow biologists to perform a more comprehensive and unbiased characterization of tissue remodeling in many biomedical contexts. AVAILABILITY AND IMPLEMENTATION: http://bioconductor.org/packages/debCAM. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Targeted realignment of LC-MS profiles by neighbor-wise compound-specific graphical time warping with misalignment detection.

MOTIVATION: Liquid chromatography-mass spectrometry (LC-MS) is a standard method for proteomics and metabolomics analysis of biological samples. Unfortunately, it suffers from various changes in the retention times (RT) of the same compound in different samples, and these must be subsequently corrected (aligned) during data processing. Classic alignment methods such as in the popular XCMS package often assume a single time-warping function for each sample. Thus, the potentially varying RT drift for compounds with different masses in a sample is neglected in these methods. Moreover, the systematic change in RT drift across run order is often not considered by alignment algorithms. Therefore, these methods cannot effectively correct all misalignments. For a large-scale experiment involving many samples, the existence of misalignment becomes inevitable and concerning. RESULTS: Here, we describe an integrated reference-free profile alignment method, neighbor-wise compound-specific Graphical Time Warping (ncGTW), that can detect misaligned features and align profiles by leveraging expected RT drift structures and compound-specific warping functions. Specifically, ncGTW uses individualized warping functions for different compounds and assigns constraint edges on warping functions of neighboring samples. Validated with both realistic synthetic data and internal quality control samples, ncGTW applied to two large-scale metabolomics LC-MS datasets identifies many misaligned features and successfully realigns them. These features would otherwise be discarded or uncorrected using existing methods. The ncGTW software tool is developed currently as a plug-in to detect and realign misaligned features present in standard XCMS output. AVAILABILITY AND IMPLEMENTATION: An R package of ncGTW is freely available at Bioconductor and https://github.com/ChiungTingWu/ncGTW. A detailed user's manual and a vignette are provided within the package. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.