Identifiability and estimability of Bayesian linear and nonlinear crossed random effects models.

Crossed random effects models (CREMs) are particularly useful in longitudinal data applications because they allow researchers to account for the impact of dynamic group membership on individual outcomes. However, no research has determined what data conditions need to be met to sufficiently identify these models, especially the group effects, in a longitudinal context. This is a significant gap in the current literature as future applications to real data may need to consider these conditions to yield accurate and precise model parameter estimates, specifically for the group effects on individual outcomes. Furthermore, there are no existing CREMs that can model intrinsically nonlinear growth. The goals of this study are to develop a Bayesian piecewise CREM to model intrinsically nonlinear growth and evaluate what data conditions are necessary to empirically identify both intrinsically linear and nonlinear longitudinal CREMs. This study includes an applied example that utilizes the piecewise CREM with real data and three simulation studies to assess the data conditions necessary to estimate linear, quadratic, and piecewise CREMs. Results show that the number of repeated measurements collected on groups impacts the ability to recover the group effects. Additionally, functional form complexity impacted data collection requirements for estimating longitudinal CREMs.

Multiple augmented reduced rank regression for pan-cancer analysis.

Statistical approaches that successfully combine multiple datasets are more powerful, efficient, and scientifically informative than separate analyses. To address variation architectures correctly and comprehensively for high-dimensional data across multiple sample sets (ie, cohorts), we propose multiple augmented reduced rank regression (maRRR), a flexible matrix regression and factorization method to concurrently learn both covariate-driven and auxiliary structured variations. We consider a structured nuclear norm objective that is motivated by random matrix theory, in which the regression or factorization terms may be shared or specific to any number of cohorts. Our framework subsumes several existing methods, such as reduced rank regression and unsupervised multimatrix factorization approaches, and includes a promising novel approach to regression and factorization of a single dataset (aRRR) as a special case. Simulations demonstrate substantial gains in power from combining multiple datasets, and from parsimoniously accounting for all structured variations. We apply maRRR to gene expression data from multiple cancer types (ie, pan-cancer) from The Cancer Genome Atlas, with somatic mutations as covariates. The method performs well with respect to prediction and imputation of held-out data, and provides new insights into mutation-driven and auxiliary variations that are shared or specific to certain cancer types.

Prognostic Performance of Hematological and Serum Iron and Metabolite Indices for Detection of Early Iron Deficiency Induced Metabolic Brain Dysfunction in Infant Rhesus Monkeys.

BACKGROUND: The current pediatric practice of monitoring for infantile iron deficiency (ID) via hemoglobin (Hgb) screening at one y of age does not identify preanemic ID nor protect against later neurocognitive deficits. OBJECTIVES: To identify biomarkers of iron-related metabolic alterations in the serum and brain and determine the sensitivity of conventional iron and heme indices for predicting risk of brain metabolic dysfunction using a nonhuman primate model of infantile ID. METHODS: Simultaneous serum iron and RBC indices, and serum and cerebrospinal fluid (CSF) metabolomic profiles were determined in 20 rhesus infants, comparing iron sufficient (IS; N = 10) and ID (N = 10) infants at 2 and 4 mo of age. RESULTS: Reticulocyte hemoglobin (RET-He) was lower at 2 wk in the ID group. Significant IS compared with ID differences in serum iron indices were present at 2 mo, but Hgb and RBC indices differed only at 4 mo (P < 0.05). Serum and CSF metabolomic profiles of the ID and IS groups differed at 2 and 4 mo (P < 0.05). Key metabolites, including homostachydrine and stachydrine (4-5-fold lower at 4 mo in ID group, P < 0.05), were altered in both serum and CSF. Iron indices and RET-He at 2 mo, but not Hgb or other RBC indices, were correlated with altered CSF metabolic profile at 4 mo and had comparable predictive accuracy (area under the receiver operating characteristic curve scores, 0.75-0.80). CONCLUSIONS: Preanemic ID at 2 mo was associated with metabolic alterations in serum and CSF in infant monkeys. Among the RBC indices, only RET-He predicted the future risk of abnormal CSF metabolic profile with a predictive accuracy comparable to serum iron indices. The concordance of homostachydrine and stachydrine changes in serum and CSF indicates their potential use as early biomarkers of brain metabolic dysfunction in infantile ID.

Multiple Augmented Reduced Rank Regression for Pan-Cancer Analysis.

Statistical approaches that successfully combine multiple datasets are more powerful, efficient, and scientifically informative than separate analyses. To address variation architectures correctly and comprehensively for high-dimensional data across multiple sample sets (i.e., cohorts), we propose multiple augmented reduced rank regression (maRRR), a flexible matrix regression and factorization method to concurrently learn both covariate-driven and auxiliary structured variation. We consider a structured nuclear norm objective that is motivated by random matrix theory, in which the regression or factorization terms may be shared or specific to any number of cohorts. Our framework subsumes several existing methods, such as reduced rank regression and unsupervised multi-matrix factorization approaches, and includes a promising novel approach to regression and factorization of a single dataset (aRRR) as a special case. Simulations demonstrate substantial gains in power from combining multiple datasets, and from parsimoniously accounting for all structured variation. We apply maRRR to gene expression data from multiple cancer types (i.e., pan-cancer) from TCGA, with somatic mutations as covariates. The method performs well with respect to prediction and imputation of held-out data, and provides new insights into mutation-driven and auxiliary variation that is shared or specific to certain cancer types.

Intranasal insulin treatment partially corrects the altered gene expression profile in the hippocampus of developing rats with perinatal iron deficiency.

Perinatal iron deficiency (FeD) targets the hippocampus and leads to long-term cognitive deficits. Intranasal insulin administration improves cognitive deficits in adult humans with Alzheimer's disease and type 2 diabetes and could provide benefits in FeD-induced hippocampal dysfunction. To objective was to assess the effects of intranasal insulin administration intranasal insulin administration on the hippocampal transcriptome in a developing rat model of perinatal FeD. Perinatal FeD was induced using low-iron diet from gestational day 3 until postnatal day (P) 7, followed by an iron sufficient (FeS) diet through P21. Intranasal insulin was administered at a dose of 0.3 IU twice daily from P8 to P21. Hippocampi were removed on P21 from FeS control, FeD control, FeS insulin, and FeD insulin groups. Total RNA was isolated and profiled using next-generation sequencing. Gene expression profiles were characterized using custom workflows and expression patterns examined using ingenuity pathways analysis (n = 7-9 per group). Select RNAseq results were confirmed via qPCR. Transcriptomic profiling revealed that mitochondrial biogenesis and flux, oxidative phosphorylation, quantity of neurons, CREB signaling in neurons, and RICTOR-based mTOR signaling were disrupted with FeD and positively affected by intranasal insulin treatment with the most benefit observed in the FeD insulin group. Both perinatal FeD and intranasal insulin administration altered gene expression profile in the developing hippocampus. Intranasal insulin treatment reversed the adverse effects of FeD on many molecular pathways and could be explored as an adjunct therapy in perinatal FeD.

A biomarker discovery framework for childhood anxiety.

Introduction: Anxiety is the most common manifestation of psychopathology in youth, negatively affecting academic, social, and adaptive functioning and increasing risk for mental health problems into adulthood. Anxiety disorders are diagnosed only after clinical symptoms emerge, potentially missing opportunities to intervene during critical early prodromal periods. In this study, we used a new empirical approach to extracting nonlinear features of the electroencephalogram (EEG), with the goal of discovering differences in brain electrodynamics that distinguish children with anxiety disorders from healthy children. Additionally, we examined whether this approach could distinguish children with externalizing disorders from healthy children and children with anxiety. Methods: We used a novel supervised tensor factorization method to extract latent factors from repeated multifrequency nonlinear EEG measures in a longitudinal sample of children assessed in infancy and at ages 3, 5, and 7 years of age. We first examined the validity of this method by showing that calendar age is highly correlated with latent EEG complexity factors (r = 0.77). We then computed latent factors separately for distinguishing children with anxiety disorders from healthy controls using a 5-fold cross validation scheme and similarly for distinguishing children with externalizing disorders from healthy controls. Results: We found that latent factors derived from EEG recordings at age 7 years were required to distinguish children with an anxiety disorder from healthy controls; recordings from infancy, 3 years, or 5 years alone were insufficient. However, recordings from two (5, 7 years) or three (3, 5, 7 years) recordings gave much better results than 7 year recordings alone. Externalizing disorders could be detected using 3- and 5 years EEG data, also giving better results with two or three recordings than any single snapshot. Further, sex assigned at birth was an important covariate that improved accuracy for both disorder groups, and birthweight as a covariate modestly improved accuracy for externalizing disorders. Recordings from infant EEG did not contribute to the classification accuracy for either anxiety or externalizing disorders. Conclusion: This study suggests that latent factors extracted from EEG recordings in childhood are promising candidate biomarkers for anxiety and for externalizing disorders if chosen at appropriate ages.

Human Cytomegalovirus in breast milk is associated with milk composition, the infant gut microbiome, and infant growth.

Human cytomegalovirus (CMV) is a highly prevalent herpesvirus that is often transmitted to the neonate via breast milk. Postnatal CMV transmission can have negative health consequences for preterm and immunocompromised infants, but any effects on healthy term infants are thought to be benign. Furthermore, the impact of CMV on the composition of the hundreds of bioactive factors in human milk has not been tested. Here, we utilize a cohort of exclusively breastfeeding full term mother-infant pairs to test for differences in the milk transcriptome and metabolome associated with CMV, and the impact of CMV in breast milk on the infant gut microbiome and infant growth. We find upregulation of the indoleamine 2,3- dioxygenase (IDO) tryptophan-to-kynurenine metabolic pathway in CMV+ milk samples, and that CMV+ milk is associated with decreased Bifidobacterium in the infant gut. Our data indicate a complex relationship between milk CMV, milk kynurenine, and infant growth; with kynurenine positively correlated, and CMV viral load negatively correlated, with infant weight-for-length at 1 month of age. These results suggest CMV transmission, CMV-related changes in milk composition, or both may be modulators of full term infant development.

Bayesian predictive modeling of multi-source multi-way data.

A Bayesian approach to predict a continuous or binary outcome from data that are collected from multiple sources with a multi-way (i.e., multidimensional tensor) structure is described. As a motivating example, molecular data from multiple 'omics sources, each measured over multiple developmental time points, as predictors of early-life iron deficiency (ID) in a rhesus monkey model are considered. The method uses a linear model with a low-rank structure on the coefficients to capture multi-way dependence and model the variance of the coefficients separately across each source to infer their relative contributions. Conjugate priors facilitate an efficient Gibbs sampling algorithm for posterior inference, assuming a continuous outcome with normal errors or a binary outcome with a probit link. Simulations demonstrate that the model performs as expected in terms of misclassification rates and correlation of estimated coefficients with true coefficients, with large gains in performance by incorporating multi-way structure and modest gains when accounting for differing signal sizes across the different sources. Moreover, it provides robust classification of ID monkeys for the motivating application.

Multiway sparse distance weighted discrimination.

Modern data often take the form of a multiway array. However, most classification methods are designed for vectors, i.e., 1-way arrays. Distance weighted discrimination (DWD) is a popular high-dimensional classification method that has been extended to the multiway context, with dramatic improvements in performance when data have multiway structure. However, the previous implementation of multiway DWD was restricted to classification of matrices, and did not account for sparsity. In this paper, we develop a general framework for multiway classification which is applicable to any number of dimensions and any degree of sparsity. We conducted extensive simulation studies, showing that our model is robust to the degree of sparsity and improves classification accuracy when the data have multiway structure. For our motivating application, magnetic resonance spectroscopy (MRS) was used to measure the abundance of several metabolites across multiple neurological regions and across multiple time points in a mouse model of Friedreich's ataxia, yielding a four-way data array. Our method reveals a robust and interpretable multi-region metabolomic signal that discriminates the groups of interest. We also successfully apply our method to gene expression time course data for multiple sclerosis treatment. An R implementation is available in the package MultiwayClassification at http://github.com/lockEF/MultiwayClassification.

Care Management Processes Important for High-Quality Diabetes Care.

OBJECTIVE: Identify the improvement in diabetes performance measures and population-based clinical outcomes resulting from changes in care management processes (CMP) in primary care practices over 3 years. RESEARCH DESIGN AND METHODS: This repeated cross-sectional study tracked clinical performance measures for all diabetes patients seen in a cohort of 330 primary care practices in 2017 and 2019. Unit of analysis was patient-year with practice-level CMP exposures. Causal inference is based on dynamic changes in individual CMPs between years by practice. We used the Bayesian method to simultaneously estimate a five-outcome model: A1c, systolic and diastolic blood pressure, guideline-based statin use, and Optimal Diabetes Care (ODC). We control for unobserved time-invariant practice characteristics and secular change. We modeled correlation of errors across outcomes. Statistical significance was identified using 99% Bayesian credible intervals (analogous to P < 0.01). RESULTS: Implementation of 18 of 62 CMPs was associated with statistically significant improvements in patient outcomes. Together, these resulted in 12.1% more patients meeting ODC performance measures. Different CMPs affected different outcomes. Three CMPs accounted for 47% of the total ODC improvement, 68% of A1c decrease, 21% of SBP reduction, and 55% of statin use increase: 1) systems for identifying and reminding patients due for testing, 2) after-visit follow-up by a nonclinician, and 3) guideline-based clinician reminders for preventive services during a clinic visit. CONCLUSIONS: Effective quality improvement in primary care focuses on practice redesign that clearly improves diabetes outcomes. Tailoring CMP adoption in primary care provides effective improvement in ODC performance through focused changes in diabetes outcomes.

Reticulocyte Hemoglobin Equivalent has Comparable Predictive Accuracy as Conventional Serum Iron Indices for Predicting Iron Deficiency and Anemia in a Nonhuman Primate model of Infantile Iron Deficiency.

BACKGROUND: Infantile iron deficiency (ID) causes anemia and compromises neurodevelopment. Current screening relies on hemoglobin (Hgb) determination at 1 year of age, which lacks sensitivity and specificity for timely detection of infantile ID. Low reticulocyte Hgb equivalent (RET-He) indicates ID, but its predictive accuracy relative to conventional serum iron indices is unknown. OBJECTIVES: The objective was to compare diagnostic accuracies of iron indices, red blood cell (RBC) indices, and RET-He for predicting the risk of ID and IDA in a nonhuman primate model of infantile ID. METHODS: Serum iron, total iron binding capacity, unsaturated iron binding capacity, transferrin saturation (TSAT), Hgb, RET-He, and other RBC indices were determined at 2 wk and 2, 4, and 6 mo in breastfed male and female rhesus infants (N = 54). The diagnostic accuracies of RET-He, iron, and RBC indices for predicting the development of ID (TSAT < 20%) and IDA (Hgb < 10 g/dL + TSAT < 20%) were determined using t tests, area under the receiver operating characteristic curve (AUC) analysis, and multiple regression models. RESULTS: Twenty-three (42.6%) infants developed ID and 16 (29.6%) progressed to IDA. All 4 iron indices and RET-He, but not Hgb or RBC indices, predicted future risk of ID and IDA (P < 0.001). The predictive accuracy of RET-He (AUC = 0.78, SE = 0.07; P = 0.003) for IDA was comparable to that of the iron indices (AUC = 0.77-0.83, SE = 0.07; P ≤ 0.002). A RET-He threshold of 25.5 pg strongly correlated with TSAT < 20% and correctly predicted IDA in 10 of 16 infants (sensitivity: 62.5%) and falsely predicted possibility of IDA in only 4 of 38 unaffected infants (specificity: 89.5%). CONCLUSIONS: RET-He is a biomarker of impending ID/IDA in rhesus infants and can be used as a hematological parameter to screen for infantile ID.

Pneumoproteins and biomarkers of inflammation and coagulation do not predict rapid lung function decline in people living with HIV.

Chronic obstructive pulmonary disease (COPD) is among the leading causes of death worldwide and HIV is an independent risk factor for the development of COPD. However, the etiology of this increased risk and means to identify persons with HIV (PWH) at highest risk for COPD have remained elusive. Biomarkers may reveal etiologic pathways and allow better COPD risk stratification. We performed a matched case:control study of PWH in the Strategic Timing of Antiretoviral Treatment (START) pulmonary substudy. Cases had rapid lung function decline (> 40 mL/year FEV1 decline) and controls had stable lung function (+ 20 to - 20 mL/year). The analysis was performed in two distinct groups: (1) those who were virally suppressed for at least 6 months and (2) those with untreated HIV (from the START deferred treatment arm). We used linear mixed effects models to test the relationship between case:control status and blood concentrations of pneumoproteins (surfactant protein-D and club cell secretory protein), and biomarkers of inflammation (IL-6 and hsCRP) and coagulation (d-dimer and fibrinogen); concentrations were measured within ± 6 months of first included spirometry. We included an interaction with treatment group (untreated HIV vs viral suppression) to test if associations varied by treatment group. This analysis included 77 matched case:control pairs in the virally suppressed batch, and 42 matched case:control pairs in the untreated HIV batch (n = 238 total) who were followed for a median of 3 years. Median (IQR) CD4 + count was lowest in the controls with untreated HIV at 674 (580, 838). We found no significant associations between case:control status and pneumoprotein or biomarker concentrations in either virally suppressed or untreated PWH. In this cohort of relatively young, recently diagnosed PWH, concentrations of pneumoproteins and biomarkers of inflammation and coagulation were not associated with subsequent rapid lung function decline.Trial registration: NCT00867048 and NCT01797367.

Lung proteome and metabolome endotype in HIV-associated obstructive lung disease.

Purpose: Obstructive lung disease is increasingly common among persons with HIV, both smokers and nonsmokers. We used aptamer proteomics to identify proteins and associated pathways in HIV-associated obstructive lung disease. Methods: Bronchoalveolar lavage fluid (BALF) samples from 26 persons living with HIV with obstructive lung disease were matched to persons living with HIV without obstructive lung disease based on age, smoking status and antiretroviral treatment. 6414 proteins were measured using SomaScan® aptamer-based assay. We used sparse distance-weighted discrimination (sDWD) to test for a difference in protein expression and permutation tests to identify univariate associations between proteins and forced expiratory volume in 1 s % predicted (FEV1 % pred). Significant proteins were entered into a pathway over-representation analysis. We also constructed protein-driven endotypes using K-means clustering and performed over-representation analysis on the proteins that were significantly different between clusters. We compared protein-associated clusters to those obtained from BALF and plasma metabolomics data on the same patient cohort. Results: After filtering, we retained 3872 proteins for further analysis. Based on sDWD, protein expression was able to separate cases and controls. We found 575 proteins that were significantly correlated with FEV1 % pred after multiple comparisons adjustment. We identified two protein-driven endotypes, one of which was associated with poor lung function, and found that insulin and apoptosis pathways were differentially represented. We found similar clusters driven by metabolomics in BALF but not plasma. Conclusion: Protein expression differs in persons living with HIV with and without obstructive lung disease. We were not able to identify specific pathways differentially expressed among patients based on FEV1 % pred; however, we identified a unique protein endotype associated with insulin and apoptotic pathways.

Human milk variation is shaped by maternal genetics and impacts the infant gut microbiome.

Human milk is a complex mix of nutritional and bioactive components that provide complete nutrition for the infant. However, we lack a systematic knowledge of the factors shaping milk composition and how milk variation influences infant health. Here, we used multi-omic profiling to characterize interactions between maternal genetics, milk gene expression, milk composition, and the infant fecal microbiome in 242 exclusively breastfeeding mother-infant pairs. We identified 487 genetic loci associated with milk gene expression unique to the lactating mammary gland, including loci that impacted breast cancer risk and human milk oligosaccharide concentration. Integrative analyses uncovered connections between milk gene expression and infant gut microbiome, including an association between the expression of inflammation-related genes with IL-6 concentration in milk and the abundance of Bifidobacteria in the infant gut. Our results show how an improved understanding of the genetics and genomics of human milk connects lactation biology with maternal and infant health.

Bayesian Distance Weighted Discrimination.

Distance weighted discrimination (DWD) is a linear discrimination method that is particularly well-suited for classification tasks with high-dimensional data. The DWD coefficients minimize an intuitive objective function, which can solved efficiently using state-of-the-art optimization techniques. However, DWD has not yet been cast into a model-based framework for statistical inference. In this article we show that DWD identifies the mode of a proper Bayesian posterior distribution, that results from a particular link function for the class probabilities and a shrinkage-inducing proper prior distribution on the coefficients. We describe a relatively efficient Markov chain Monte Carlo (MCMC) algorithm to simulate from the true posterior under this Bayesian framework. We show that the posterior is asymptotically normal and derive the mean and covariance matrix of its limiting distribution. Through several simulation studies and an application to breast cancer genomics we demonstrate how the Bayesian approach to DWD can be used to (1) compute well-calibrated posterior class probabilities, (2) assess uncertainty in the DWD coefficients and resulting sample scores, (3) improve power via semi-supervised analysis when not all class labels are available, and (4) automatically determine a penalty tuning parameter within the model-based framework. R code to perform Bayesian DWD is available at https://github.com/lockEF/BayesianDWD.

A natural history study to track brain and spinal cord changes in individuals with Friedreich's ataxia: TRACK-FA study protocol.

INTRODUCTION: Drug development for neurodegenerative diseases such as Friedreich's ataxia (FRDA) is limited by a lack of validated, sensitive biomarkers of pharmacodynamic response in affected tissue and disease progression. Studies employing neuroimaging measures to track FRDA have thus far been limited by their small sample sizes and limited follow up. TRACK-FA, a longitudinal, multi-site, and multi-modal neuroimaging natural history study, aims to address these shortcomings by enabling better understanding of underlying pathology and identifying sensitive, clinical trial ready, neuroimaging biomarkers for FRDA. METHODS: 200 individuals with FRDA and 104 control participants will be recruited across seven international study sites. Inclusion criteria for participants with genetically confirmed FRDA involves, age of disease onset ≤ 25 years, Friedreich's Ataxia Rating Scale (FARS) functional staging score of ≤ 5, and a total modified FARS (mFARS) score of ≤ 65 upon enrolment. The control cohort is matched to the FRDA cohort for age, sex, handedness, and years of education. Participants will be evaluated at three study visits over two years. Each visit comprises of a harmonized multimodal Magnetic Resonance Imaging (MRI) and Spectroscopy (MRS) scan of the brain and spinal cord; clinical, cognitive, mood and speech assessments and collection of a blood sample. Primary outcome measures, informed by previous neuroimaging studies, include measures of: spinal cord and brain morphometry, spinal cord and brain microstructure (measured using diffusion MRI), brain iron accumulation (using Quantitative Susceptibility Mapping) and spinal cord biochemistry (using MRS). Secondary and exploratory outcome measures include clinical, cognitive assessments and blood biomarkers. DISCUSSION: Prioritising immediate areas of need, TRACK-FA aims to deliver a set of sensitive, clinical trial-ready neuroimaging biomarkers to accelerate drug discovery efforts and better understand disease trajectory. Once validated, these potential pharmacodynamic biomarkers can be used to measure the efficacy of new therapeutics in forestalling disease progression. CLINICAL TRIAL REGISTRATION: ClinicalTrails.gov Identifier: NCT04349514.

sJIVE: Supervised Joint and Individual Variation Explained.

Analyzing multi-source data, which are multiple views of data on the same subjects, has become increasingly common in molecular biomedical research. Recent methods have sought to uncover underlying structure and relationships within and/or between the data sources, and other methods have sought to build a predictive model for an outcome using all sources. However, existing methods that do both are presently limited because they either (1) only consider data structure shared by all datasets while ignoring structures unique to each source, or (2) they extract underlying structures first without consideration to the outcome. The proposed method, supervised joint and individual variation explained (sJIVE), can simultaneously (1) identify shared (joint) and source-specific (individual) underlying structure and (2) build a linear prediction model for an outcome using these structures. These two components are weighted to compromise between explaining variation in the multi-source data and in the outcome. Simulations show sJIVE to outperform existing methods when large amounts of noise are present in the multi-source data. An application to data from the COPDGene study explores gene expression and proteomic patterns associated with lung function.

Tandem mass tag proteomic and untargeted metabolomic profiling reveals altered serum and CSF biochemical datasets in iron deficient monkeys.

The effects of early-life iron deficiency anemia (IDA) extend past the blood and include both short- and long-term adverse effects on many tissues including the brain. Prior to IDA, iron deficiency (ID) can cause similar tissue effects, but a sensitive biomarker of iron-dependent brain health is lacking. To determine serum and CSF biomarkers of ID-induced metabolic dysfunction we performed proteomic and metabolomic analysis of serum and CSF at 4- and 6- months from a nonhuman primate model of infantile IDA. LC/MS/MS analyses identified a total of 227 metabolites and 205 proteins in serum. In CSF, we measured 210 metabolites and 1,560 proteins. Data were either processed from a Q-Exactive (Thermo Scientific, Waltham, MA) through Progenesis QI with accurate mass and retention time comparisons to a proprietary small molecule database and Metlin or with raw files imported directly from a Fusion Orbitrap (Thermo Scientific, Waltham, MA) through Sequest in Proteome Discoverer 2.4.0.305 (Thermo Scientific, Waltham, MA) with peptide matches through the latest Rhesus Macaque HMDB database. Metabolite and protein identifiers, p-values, and q-values were utilized for molecular pathway analysis with Ingenuity Pathways Analysis (IPA). We applied multiway distance weighted discrimination (DWD) to identify a weighted sum of the features (proteins or metabolites) that distinguish ID from IS at 4-months (pre-anemic period) and 6-months of age (anemic).

Lung and Plasma Metabolome in HIV-Associated Obstructive Lung Disease.

BACKGROUND: HIV is a risk factor for obstructive lung disease (OLD), independent of smoking. We used mass spectrometry (MS) approaches to identify metabolomic biomarkers that inform mechanistic pathogenesis of OLD in persons with HIV (PWH). METHODS: We obtained bronchoalveolar lavage fluid (BALF) samples from 52 PWH, in case:control (+OLD/-OLD) pairs matched on age, smoking status, and antiretroviral treatment. Four hundred nine metabolites from 8 families were measured on BALF and plasma samples using a MS-based Biocrates platform. After filtering metabolites with a high proportion of missing values and values below the level of detection, we performed univariate testing using paired t tests followed by false discovery rate corrections. We used distance-weighted discrimination (DWD) to test for an overall difference in the metabolite profile between cases and controls. RESULTS: After filtering, there were 252 BALF metabolites for analysis from 8 metabolite families. DWD testing found that collectively, BALF metabolites differentiated cases from controls, whereas plasma metabolites did not. In BALF samples, we identified 3 metabolites that correlated with OLD at the false discovery rate of 10%; all were in the phosphatidylcholine family. We identified additional BALF metabolites when analyzing lung function as a continuous variable, and these included acylcarnitines, triglycerides, and a cholesterol ester. CONCLUSIONS: Collectively, BALF metabolites differentiate PWH with and without OLD. These included several BALF lipid metabolites. These findings were limited to BALF and were not found in plasma from the same individuals. Phosphatidylcholine, the most common lipid component of surfactant, was the predominant lipid metabolite differentially expressed.

A hierarchical spike-and-slab model for pan-cancer survival using pan-omic data.

BACKGROUND: Pan-omics, pan-cancer analysis has advanced our understanding of the molecular heterogeneity of cancer. However, such analyses have been limited in their ability to use information from multiple sources of data (e.g., omics platforms) and multiple sample sets (e.g., cancer types) to predict clinical outcomes. We address the issue of prediction across multiple high-dimensional sources of data and sample sets by using molecular patterns identified by BIDIFAC+, a method for integrative dimension reduction of bidimensionally-linked matrices, in a Bayesian hierarchical model. Our model performs variable selection through spike-and-slab priors that borrow information across clustered data. We use this model to predict overall patient survival from the Cancer Genome Atlas with data from 29 cancer types and 4 omics sources and use simulations to characterize the performance of the hierarchical spike-and-slab prior. RESULTS: We found that molecular patterns shared across all or most cancers were largely not predictive of survival. However, our model selected patterns unique to subsets of cancers that differentiate clinical tumor subtypes with markedly different survival outcomes. Some of these subtypes were previously established, such as subtypes of uterine corpus endometrial carcinoma, while others may be novel, such as subtypes within a set of kidney carcinomas. Through simulations, we found that the hierarchical spike-and-slab prior performs best in terms of variable selection accuracy and predictive power when borrowing information is advantageous, but also offers competitive performance when it is not. CONCLUSIONS: We address the issue of prediction across multiple sources of data by using results from BIDIFAC+ in a Bayesian hierarchical model for overall patient survival. By incorporating spike-and-slab priors that borrow information across cancers, we identified molecular patterns that distinguish clinical tumor subtypes within a single cancer and within a group of cancers. We also corroborate the flexibility and performance of using spike-and-slab priors as a Bayesian variable selection approach.