Biomarkers of Inflammation and Longitudinal Evaluation of Lung Function, Physical Activity, and Grip Strength: A Secondary Analysis in the CASCADE Study.

Rationale: Physical activity, lung function, and grip strength are associated with exacerbations, hospitalizations, and mortality in people with chronic obstructive pulmonary disease (COPD). We tested whether baseline inflammatory biomarkers were associated with longitudinal outcomes of these physiologic measurements. Methods: The COPD Activity: Serotonin Transporter, Cytokines, and Depression (CASCADE) study was a prospective observational study of individuals with COPD. A total of 14 inflammatory biomarkers were measured at baseline. Participants were followed for 2 years. We analyzed associations between baseline biomarkers and forced expiratory volume in 1 second (FEV1), physical activity, and grip strength. We used a hierarchical hypothesis testing procedure to reduce type I error. We used Pearson correlations to test associations between baseline biomarkers and longitudinal changes in the outcomes of interest. We used Fisher's linear discriminant analysis to test if linear combinations of baseline biomarkers predict rapid FEV1 decline. Finally, we used linear mixed modeling to test associations between baseline biomarkers and outcomes of interest at baseline, year 1, and year 2; models were adjusted for age, smoking status, baseline biomarkers, and FEV1. Results: A total of 302 participants (age 67.5 ± 8.5 years, 19.5% female, 28.5% currently smoking) were included. Baseline biomarkers were not associated with longitudinal changes in grip strength, physical activity, or rapid FEV1 decline. Higher interleukin-6 and C-reactive protein were associated with lower physical activity at baseline and these relationships persisted at year 1 and year 2. Conclusion: Baseline inflammatory biomarkers did not predict changes in lung function or physical activity, but higher inflammatory biomarkers were associated with persistently low levels of physical activity.

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.

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.

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.

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.

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.

Identification of shared and disease-specific host gene-microbiome associations across human diseases using multi-omic integration.

While gut microbiome and host gene regulation independently contribute to gastrointestinal disorders, it is unclear how the two may interact to influence host pathophysiology. Here we developed a machine learning-based framework to jointly analyse paired host transcriptomic (n = 208) and gut microbiome (n = 208) profiles from colonic mucosal samples of patients with colorectal cancer, inflammatory bowel disease and irritable bowel syndrome. We identified associations between gut microbes and host genes that depict shared as well as disease-specific patterns. We found that a common set of host genes and pathways implicated in gastrointestinal inflammation, gut barrier protection and energy metabolism are associated with disease-specific gut microbes. Additionally, we also found that mucosal gut microbes that have been implicated in all three diseases, such as Streptococcus, are associated with different host pathways in each disease, suggesting that similar microbes can affect host pathophysiology in a disease-specific manner through regulation of different host genes. Our framework can be applied to other diseases for the identification of host gene-microbiome associations that may influence disease outcomes.

BIDIMENSIONAL LINKED MATRIX FACTORIZATION FOR PAN-OMICS PAN-CANCER ANALYSIS.

Several modern applications require the integration of multiple large data matrices that have shared rows and/or columns. For example, cancer studies that integrate multiple omics platforms across multiple types of cancer, pan-omics pan-cancer analysis, have extended our knowledge of molecular heterogeneity beyond what was observed in single tumor and single platform studies. However, these studies have been limited by available statistical methodology. We propose a flexible approach to the simultaneous factorization and decomposition of variation across such bidimensionally linked matrices, BIDIFAC+. BIDIFAC+ decomposes variation into a series of low-rank components that may be shared across any number of row sets (e.g., omics platforms) or column sets (e.g., cancer types). This builds on a growing literature for the factorization and decomposition of linked matrices which has primarily focused on multiple matrices that are linked in one dimension (rows or columns) only. Our objective function extends nuclear norm penalization, is motivated by random matrix theory, gives a unique decomposition under relatively mild conditions, and can be shown to give the mode of a Bayesian posterior distribution. We apply BIDIFAC+ to pan-omics pan-cancer data from TCGA, identifying shared and specific modes of variability across four different omics platforms and 29 different cancer types.

Multiomic profiling of iron-deficient infant monkeys reveals alterations in neurologically important biochemicals in serum and cerebrospinal fluid before the onset of anemia.

The effects of iron deficiency (ID) during infancy extend beyond the hematologic compartment and include short- and long-term adverse effects on many tissues including the brain. However, sensitive biomarkers of iron-dependent brain health are lacking in humans. To determine whether serum and cerebrospinal fluid (CSF) biomarkers of ID-induced metabolic dysfunction are concordant in the pre/early anemic stage of ID before anemia in a nonhuman primate model of infantile iron deficiency anemia (IDA). ID (n = 7), rhesus infants at 4 mo (pre-anemic period) and 6 mo of age (anemic) were examined. Hematological, metabolomic, and proteomic profiles were generated via HPLC/MS at both time points to discriminate serum biomarkers of ID-induced brain metabolic dysfunction. We identified 227 metabolites and 205 proteins in serum. Abnormalities indicating altered liver function, lipid dysregulation, and increased acute phase reactants were present in ID. In CSF, we measured 210 metabolites and 1,560 proteins with changes in ID infants indicative of metabolomic and proteomic differences indexing disrupted synaptogenesis. Systemic and CSF proteomic and metabolomic changes were present and concurrent in the pre-anemic and anemic periods. Multiomic serum and CSF profiling uncovered pathways disrupted by ID in both the pre-anemic and anemic stages of infantile IDA, including evidence for hepatic dysfunction and activation of acute phase response. Parallel changes observed in serum and CSF potentially provide measurable serum biomarkers of ID that reflect at-risk brain processes prior to progression to clinical anemia.

Correcting iron deficiency anemia with iron dextran alters the serum metabolomic profile of the infant Rhesus Monkey.

BACKGROUND: The effects of infantile iron deficiency anemia (IDA) extend beyond hematological indices and include short- and long-term adverse effects on multiple cells and tissues. IDA is associated with an abnormal serum metabolomic profile, characterized by altered hepatic metabolism, lowered NAD flux, increased nucleoside levels, and a reduction in circulating dopamine levels. OBJECTIVES: The objective of this study was to determine whether the serum metabolomic profile is normalized after rapid correction of IDA using iron dextran injections. METHODS: Blood was collected from iron-sufficient (IS; n = 10) and IDA (n = 12) rhesus infants at 6 months of age. IDA infants were then administered iron dextran and vitamin B via intramuscular injections at weekly intervals for 2 to 8 weeks. Their hematological and metabolomic statuses were evaluated following treatment and compared with baseline and a separate group of age-matched IS infants (n = 5). RESULTS: Serum metabolomic profiles assessed at baseline and after treatment via HPLC/MS using isobaric standards identified 654 quantifiable metabolites. At baseline, 53 metabolites differed between IS and IDA infants. Iron treatment restored traditional hematological indices, including hemoglobin and mean corpuscular volume, into the normal range, but the metabolite profile in the IDA group after iron treatment was markedly altered, with 323 metabolites differentially expressed when compared with an infant's own baseline profile. CONCLUSIONS: Rapid correction of IDA with iron dextran resulted in extensive metabolic changes across biochemical pathways indexing the liver function, bile acid release, essential fatty acid production, nucleoside release, and several neurologically important metabolites. The results highlight the importance of a cautious approach when developing a route and regimen of iron repletion to treat infantile IDA.

A Pan-Cancer and Polygenic Bayesian Hierarchical Model for the Effect of Somatic Mutations on Survival.

We built a novel Bayesian hierarchical survival model based on the somatic mutation profile of patients across 50 genes and 27 cancer types. The pan-cancer quality allows for the model to "borrow" information across cancer types, motivated by the assumption that similar mutation profiles may have similar (but not necessarily identical) effects on survival across different tissues of origin or tumor types. The effect of a mutation at each gene was allowed to vary by cancer type, whereas the mean effect of each gene was shared across cancers. Within this framework, we considered 4 parametric survival models (normal, log-normal, exponential, and Weibull), and we compared their performance via a cross-validation approach in which we fit each model on training data and estimate the log-posterior predictive likelihood on test data. The log-normal model gave the best fit, and we investigated the partial effect of each gene on survival via a forward selection procedure. Through this we determined that mutations at TP53 and FAT4 were together the most useful for predicting patient survival. We validated the model via simulation to ensure that our algorithm for posterior computation gave nominal coverage rates. The code used for this analysis can be found at https://github.com/sarahsamorodnitsky/Pan-Cancer-Survival-Modeling.git, and the results are summarized at http://ericfrazerlock.com/surv_figs/SurvivalDisplay.html.

Integrative factorization of bidimensionally linked matrices.

Advances in molecular "omics" technologies have motivated new methodologies for the integration of multiple sources of high-content biomedical data. However, most statistical methods for integrating multiple data matrices only consider data shared vertically (one cohort on multiple platforms) or horizontally (different cohorts on a single platform). This is limiting for data that take the form of bidimensionally linked matrices (eg, multiple cohorts measured on multiple platforms), which are increasingly common in large-scale biomedical studies. In this paper, we propose bidimensional integrative factorization (BIDIFAC) for integrative dimension reduction and signal approximation of bidimensionally linked data matrices. Our method factorizes data into (a) globally shared, (b) row-shared, (c) column-shared, and (d) single-matrix structural components, facilitating the investigation of shared and unique patterns of variability. For estimation, we use a penalized objective function that extends the nuclear norm penalization for a single matrix. As an alternative to the complicated rank selection problem, we use results from the random matrix theory to choose tuning parameters. We apply our method to integrate two genomics platforms (messenger RNA and microRNA expression) across two sample cohorts (tumor samples and normal tissue samples) using the breast cancer data from the Cancer Genome Atlas. We provide R code for fitting BIDIFAC, imputing missing values, and generating simulated data.

Early-Life Iron Deficiency and Its Natural Resolution Are Associated with Altered Serum Metabolomic Profiles in Infant Rhesus Monkeys.

BACKGROUND: Iron deficiency is the most common nutrient deficiency in human infants aged 6 to 24 mo, and negatively affects many cellular metabolic processes, including energy production, electron transport, and oxidative degradation of toxins. There can be persistent influences on long-term metabolic health beyond its acute effects. OBJECTIVES: The objective was to determine how iron deficiency in infancy alters the serum metabolomic profile and to test whether these effects persist after the resolution of iron deficiency in a nonhuman primate model of spontaneous iron deficiency. METHODS: Blood was collected from naturally iron-sufficient (IS; n = 10) and iron-deficient (ID; n = 10) male and female infant rhesus monkeys (Macaca mulatta) at 6 mo of age. Iron deficiency resolved without intervention upon feeding of solid foods, and iron status was re-evaluated at 12 mo of age from the IS and formerly ID monkeys using hematological and other indices; sera were metabolically profiled using HPLC/MS and GC/MS with isobaric standards for identification and quantification at both time points. RESULTS: A total of 413 metabolites were measured, with differences in 40 metabolites identified between IS and ID monkeys at 6 mo (P$\le $ 0.05). At 12 mo, iron-related hematological parameters had returned to normal, but the formerly ID infants remained metabolically distinct from the age-matched IS infants, with 48 metabolites differentially expressed between the groups. Metabolomic profiling indicated altered liver metabolites, differential fatty acid production, increased serum uridine release, and atypical bile acid production in the ID monkeys. CONCLUSIONS: Pathway analyses of serum metabolites provided evidence of a hypometabolic state, altered liver function, differential essential fatty acid production, irregular uracil metabolism, and atypical bile acid production in ID infants. Many metabolites remained altered after the resolution of ID, suggesting long-term effects on metabolic health.