A Multiomics, Molecular Atlas of Breast Cancer Survivors.

Breast cancer imposes a significant burden globally. While the survival rate is steadily improving, much remains to be elucidated. This observational, single time point, multiomic study utilizing genomics, proteomics, targeted and untargeted metabolomics, and metagenomics in a breast cancer survivor (BCS) and age-matched healthy control cohort (N = 100) provides deep molecular phenotyping of breast cancer survivors. In this study, the BCS cohort had significantly higher polygenic risk scores for breast cancer than the control group. Carnitine and hexanoyl carnitine were significantly different. Several bile acid and fatty acid metabolites were significantly dissimilar, most notably the Omega-3 Index (O3I) (significantly lower in BCS). Proteomic and metagenomic analyses identified group and pathway differences, which warrant further investigation. The database built from this study contributes a wealth of data on breast cancer survivorship where there has been a paucity, affording the ability to identify patterns and novel insights that can drive new hypotheses and inform future research. Expansion of this database in the treatment-naïve, newly diagnosed, controlling for treatment confounders, and through the disease progression, can be leveraged to profile and contextualize breast cancer and breast cancer survivorship, potentially leading to the development of new strategies to combat this disease and improve the quality of life for its victims.

A data-driven approach to improve wellness and reduce recurrence in cancer survivors.

For many cancer survivors, toxic side effects of treatment, lingering effects of the aftermath of disease and cancer recurrence adversely affect quality of life (QoL) and reduce healthspan. Data-driven approaches for quantifying and improving wellness in healthy individuals hold great promise for improving the lives of cancer survivors. The data-driven strategy will also guide personalized nutrition and exercise recommendations that may help prevent cancer recurrence and secondary malignancies in survivors.

Integrated unbiased multiomics defines disease-independent placental clusters in common obstetrical syndromes.

BACKGROUND: Placental dysfunction, a root cause of common syndromes affecting human pregnancy, such as preeclampsia (PE), fetal growth restriction (FGR), and spontaneous preterm delivery (sPTD), remains poorly defined. These common, yet clinically disparate obstetrical syndromes share similar placental histopathologic patterns, while individuals within each syndrome present distinct molecular changes, challenging our understanding and hindering our ability to prevent and treat these syndromes. METHODS: Using our extensive biobank, we identified women with severe PE (n = 75), FGR (n = 40), FGR with a hypertensive disorder (FGR + HDP; n = 33), sPTD (n = 72), and two uncomplicated control groups, term (n = 113), and preterm without PE, FGR, or sPTD (n = 16). We used placental biopsies for transcriptomics, proteomics, metabolomics data, and histological evaluation. After conventional pairwise comparison, we deployed an unbiased, AI-based similarity network fusion (SNF) to integrate the datatypes and identify omics-defined placental clusters. We used Bayesian model selection to compare the association between the histopathological features and disease conditions vs SNF clusters. RESULTS: Pairwise, disease-based comparisons exhibited relatively few differences, likely reflecting the heterogeneity of the clinical syndromes. Therefore, we deployed the unbiased, omics-based SNF method. Our analysis resulted in four distinct clusters, which were mostly dominated by a specific syndrome. Notably, the cluster dominated by early-onset PE exhibited strong placental dysfunction patterns, with weaker injury patterns in the cluster dominated by sPTD. The SNF-defined clusters exhibited better correlation with the histopathology than the predefined disease groups. CONCLUSIONS: Our results demonstrate that integrated omics-based SNF distinctively reclassifies placental dysfunction patterns underlying the common obstetrical syndromes, improves our understanding of the pathological processes, and could promote a search for more personalized interventions.

Modulating innate immune activation states impacts the efficacy of specific Aß immunotherapy.

INTRODUCTION: Passive immunotherapies targeting Aß continue to be evaluated as Alzheimer's disease (AD) therapeutics, but there remains debate over the mechanisms by which these immunotherapies work. Besides the amount of preexisting Aß deposition and the type of deposit (compact or diffuse), there is little data concerning what factors, independent of those intrinsic to the antibody, might influence efficacy. Here we (i) explored how constitutive priming of the underlying innate activation states by Il10 and Il6 might influence passive Aß immunotherapy and (ii) evaluated transcriptomic data generated in the AMP-AD initiative to inform how these two cytokines and their receptors' mRNA levels are altered in human AD and an APP mouse model. METHODS: rAAV2/1 encoding EGFP, Il6 or Il10 were delivered by somatic brain transgenesis to neonatal (P0) TgCRND8 APP mice. Then, at 2 months of age, the mice were treated bi-weekly with a high-affinity anti-Aß1-16 mAb5 monoclonal antibody or control mouse IgG until 6 months of age. rAAV mediated transgene expression, amyloid accumulation, Aß levels and gliosis were assessed. Extensive transcriptomic data was used to evaluate the mRNA expression levels of IL10 and IL6 and their receptors in the postmortem human AD temporal cortex and in the brains of TgCRND8 mice, the later at multiple ages. RESULTS: Priming TgCRND8 mice with Il10 increases Aß loads and blocks efficacy of subsequent mAb5 passive immunotherapy, whereas priming with Il6 priming reduces Aß loads by itself and subsequent Aß immunotherapy shows only a slightly additive effect. Transcriptomic data shows that (i) there are significant increases in the mRNA levels of Il6 and Il10 receptors in the TgCRND8 mouse model and temporal cortex of humans with AD and (ii) there is a great deal of variance in individual mouse brain and the human temporal cortex of these interleukins and their receptors. CONCLUSIONS: The underlying immune activation state can markedly affect the efficacy of passive Aß immunotherapy. These results have important implications for ongoing human AD immunotherapy trials, as they indicate that underlying immune activation states within the brain, which may be highly variable, may influence the ability for passive immunotherapy to alter Aß deposition.

Integrative functional genomic analysis of intron retention in human and mouse brain with Alzheimer's disease.

Intron retention (IR) has been implicated in the pathogenesis of complex diseases such as cancers; its association with Alzheimer's disease (AD) remains unexplored. We performed genome-wide analysis of IR through integrating genetic, transcriptomic, and proteomic data of AD subjects and mouse models from the Accelerating Medicines Partnership-Alzheimer's Disease project. We identified 4535 and 4086 IR events in 2173 human and 1736 mouse genes, respectively. Quantitation of IR enabled the identification of differentially expressed genes that conventional exon-level approaches did not reveal. There were significant correlations of intron expression within innate immune genes, like HMBOX1, with AD in humans. Peptides with a high probability of translation from intron-retained mRNAs were identified using mass spectrometry. Further, we established AD-specific intron expression Quantitative Trait Loci, and identified splicing-related genes that may regulate IR. Our analysis provides a novel resource for the search for new AD biomarkers and pathological mechanisms.

Untargeted longitudinal analysis of a wellness cohort identifies markers of metastatic cancer years prior to diagnosis.

We analyzed 1196 proteins in longitudinal plasma samples from participants in a commercial wellness program, including samples collected pre-diagnosis from ten cancer patients and 69 controls. For three individuals ultimately diagnosed with metastatic breast, lung, or pancreatic cancer, CEACAM5 was a persistent longitudinal outlier as early as 26.5 months pre-diagnosis. CALCA, a biomarker for medullary thyroid cancer, was hypersecreted in metastatic pancreatic cancer at least 16.5 months pre-diagnosis. ERBB2 levels spiked in metastatic breast cancer between 10.0 and 4.0 months pre-diagnosis. Our results support the value of deep phenotyping seemingly healthy individuals in prospectively inferring disease transitions.

Multiomic blood correlates of genetic risk identify presymptomatic disease alterations.

Transitions from health to disease are characterized by dysregulation of biological networks under the influence of genetic and environmental factors, often over the course of years to decades before clinical symptoms appear. Understanding these dynamics has important implications for preventive medicine. However, progress has been hindered both by the difficulty of identifying individuals who will eventually go on to develop a particular disease and by the inaccessibility of most disease-relevant tissues in living individuals. Here we developed an alternative approach using polygenic risk scores (PRSs) based on genome-wide association studies (GWAS) for 54 diseases and complex traits coupled with multiomic profiling and found that these PRSs were associated with 766 detectable alterations in proteomic, metabolomic, and standard clinical laboratory measurements (clinical labs) from blood plasma across several thousand mostly healthy individuals. We recapitulated a variety of known relationships (e.g., glutamatergic neurotransmission and inflammation with depression, IL-33 with asthma) and found associations directly suggesting therapeutic strategies (e.g., Ω-6 supplementation and IL-13 inhibition for amyotrophic lateral sclerosis) and influences on longevity (leukemia inhibitory factor, ceramides). Analytes altered in high-genetic-risk individuals showed concordant changes in disease cases, supporting the notion that PRS-associated analytes represent presymptomatic disease alterations. Our results provide insights into the molecular pathophysiology of a range of traits and suggest avenues for the prevention of health-to-disease transitions.

iREAD: a tool for intron retention detection from RNA-seq data.

BACKGROUND: Intron retention (IR) has been traditionally overlooked as 'noise' and received negligible attention in the field of gene expression analysis. In recent years, IR has become an emerging field for interrogating transcriptomes because it has been recognized to carry out important biological functions such as gene expression regulation and it has been found to be associated with complex diseases such as cancers. However, methods for detecting IR today are limited. Thus, there is a need to develop novel methods to improve IR detection. RESULTS: Here we present iREAD (intron REtention Analysis and Detector), a tool to detect IR events genome-wide from high-throughput RNA-seq data. The command line interface for iREAD is implemented in Python. iREAD takes as input a BAM file, representing the transcriptome, and a text file containing the intron coordinates of a genome. It then 1) counts all reads that overlap intron regions, 2) detects IR events by analyzing the features of reads such as depth and distribution patterns, and 3) outputs a list of retained introns into a tab-delimited text file. iREAD provides significant added value in detecting IR compared with output from IRFinder with a higher AUC on all datasets tested. Both methods showed low false positive rates and high false negative rates in different regimes, indicating that use together is generally beneficial. The output from iREAD can be directly used for further exploratory analysis such as differential intron expression and functional enrichment. The software is freely available at https://github.com/genemine/iread. CONCLUSION: Being complementary to existing tools, iREAD provides a new and generic tool to interrogate poly-A enriched transcriptomic data of intron regions. Intron retention analysis provides a complementary approach for understanding transcriptome.

Identifying Personalized Metabolic Signatures in Breast Cancer.

Cancer cells are adept at reprogramming energy metabolism, and the precise manifestation of this metabolic reprogramming exhibits heterogeneity across individuals (and from cell to cell). In this study, we analyzed the metabolic differences between interpersonal heterogeneous cancer phenotypes. We used divergence analysis on gene expression data of 1156 breast normal and tumor samples from The Cancer Genome Atlas (TCGA) and integrated this information with a genome-scale reconstruction of human metabolism to generate personalized, context-specific metabolic networks. Using this approach, we classified the samples into four distinct groups based on their metabolic profiles. Enrichment analysis of the subsystems indicated that amino acid metabolism, fatty acid oxidation, citric acid cycle, androgen and estrogen metabolism, and reactive oxygen species (ROS) detoxification distinguished these four groups. Additionally, we developed a workflow to identify potential drugs that can selectively target genes associated with the reactions of interest. MG-132 (a proteasome inhibitor) and OSU-03012 (a celecoxib derivative) were the top-ranking drugs identified from our analysis and known to have anti-tumor activity. Our approach has the potential to provide mechanistic insights into cancer-specific metabolic dependencies, ultimately enabling the identification of potential drug targets for each patient independently, contributing to a rational personalized medicine approach.

A systems approach to clinical oncology uses deep phenotyping to deliver personalized care.

Cancer encompasses a complex, heterogeneous and dynamic group of diseases that arise from perturbations to multiple biological networks within the body. A systems biology-based approach would help to decipher this complexity, to deeply characterize the pathophysiology of the disease and to stratify cancers into appropriate molecular subtypes to facilitate the development of personalized therapies. Technological advances made over the past decade have enabled multiscale, longitudinal measurements ('snapshots') of human biology, from single-cell analyses to whole-body monitoring. In this Perspective, we discuss some of these technologies and how they have (and will) contributed to our understanding of cancer biology as well as to the development of early diagnostics and personalized therapies. We argue that the integration of molecular profiling of cancerous tissues with deep, longitudinal profiling of the physiological state of an individual ('deep phenotyping') is key to understanding the prevention, initiation, progression and response to treatment of cancers. Systems biology-based approaches can provide an unprecedented trove of data for early detection of disease transitions, prediction of therapeutic responses and clinical outcomes, and for the design of personalized treatments.

Genetic Predisposition Impacts Clinical Changes in a Lifestyle Coaching Program.

Both genetic and lifestyle factors contribute to an individual's disease risk, suggesting a multi-omic approach is essential for personalized prevention. Studies have examined the effectiveness of lifestyle coaching on clinical outcomes, however, little is known about the impact of genetic predisposition on the response to lifestyle coaching. Here we report on the results of a real-world observational study in 2531 participants enrolled in a commercial "Scientific Wellness" program, which combines multi-omic data with personalized, telephonic lifestyle coaching. Specifically, we examined: 1) the impact of this program on 55 clinical markers and 2) the effect of genetic predisposition on these clinical changes. We identified sustained improvements in clinical markers related to cardiometabolic risk, inflammation, nutrition, and anthropometrics. Notably, improvements in HbA1c were akin to those observed in landmark trials. Furthermore, genetic markers were associated with longitudinal changes in clinical markers. For example, individuals with genetic predisposition for higher LDL-C had a lesser decrease in LDL-C on average than those with genetic predisposition for average LDL-C. Overall, these results suggest that a program combining multi-omic data with lifestyle coaching produces clinically meaningful improvements, and that genetic predisposition impacts clinical responses to lifestyle change.

Single-Cell Proteomics Reveal that Quantitative Changes in Co-expressed Lineage-Specific Transcription Factors Determine Cell Fate.

Hematopoiesis provides an accessible system for studying the principles underlying cell-fate decisions in stem cells. Proposed models of hematopoiesis suggest that quantitative changes in lineage-specific transcription factors (LS-TFs) underlie cell-fate decisions. However, evidence for such models is lacking as TF levels are typically measured via RNA expression rather than by analyzing temporal changes in protein abundance. Here, we used single-cell mass cytometry and absolute quantification by mass spectrometry to capture the temporal dynamics of TF protein expression in individual cells during human erythropoiesis. We found that LS-TFs from alternate lineages are co-expressed, as proteins, in individual early progenitor cells and quantitative changes of LS-TFs occur gradually rather than abruptly to direct cell-fate decisions. Importantly, upregulation of a megakaryocytic TF in early progenitors is sufficient to deviate cells from an erythroid to a megakaryocyte trajectory, showing that quantitative changes in protein abundance of LS-TFs in progenitors can determine alternate cell fates.

Multiscale Analysis of Independent Alzheimer's Cohorts Finds Disruption of Molecular, Genetic, and Clinical Networks by Human Herpesvirus.

Investigators have long suspected that pathogenic microbes might contribute to the onset and progression of Alzheimer's disease (AD) although definitive evidence has not been presented. Whether such findings represent a causal contribution, or reflect opportunistic passengers of neurodegeneration, is also difficult to resolve. We constructed multiscale networks of the late-onset AD-associated virome, integrating genomic, transcriptomic, proteomic, and histopathological data across four brain regions from human post-mortem tissue. We observed increased human herpesvirus 6A (HHV-6A) and human herpesvirus 7 (HHV-7) from subjects with AD compared with controls. These results were replicated in two additional, independent and geographically dispersed cohorts. We observed regulatory relationships linking viral abundance and modulators of APP metabolism, including induction of APBB2, APPBP2, BIN1, BACE1, CLU, PICALM, and PSEN1 by HHV-6A. This study elucidates networks linking molecular, clinical, and neuropathological features with viral activity and is consistent with viral activity constituting a general feature of AD.

Targeting TWIST1 through loss of function inhibits tumorigenicity of human glioblastoma.

TWIST1 (TW) is a bHLH transcription factor (TF) and master regulator of the epithelial-to-mesenchymal transition (EMT). In vitro, TW promotes mesenchymal change, invasion, and self-renewal in glioblastoma (GBM) cells. However, the potential therapeutic relevance of TW has not been established through loss-of-function studies in human GBM cell xenograft models. The effects of TW loss of function (gene editing and knockdown) on inhibition of tumorigenicity of U87MG and GBM4 glioma stem cells were tested in orthotopic xenograft models and conditional knockdown in established flank xenograft tumors. RNAseq and the analysis of tumors investigated putative TW-associated mechanisms. Multiple bioinformatic tools revealed significant alteration of ECM, membrane receptors, signaling transduction kinases, and cytoskeleton dynamics leading to identification of PI3K/AKT signaling. We experimentally show alteration of AKT activity and periostin (POSTN) expression in vivo and/or in vitro. For the first time, we show that effect of TW knockout inhibits AKT activity in U87MG cells in vivo independent of PTEN mutation. The clinical relevance of TW and candidate mechanisms was established by analysis of the TCGA and ENCODE databases. TW expression was associated with decreased patient survival and LASSO regression analysis identified POSTN as one of top targets of TW in human GBM. While we previously demonstrated the role of TW in promoting EMT and invasion of glioma cells, these studies provide direct experimental evidence supporting protumorigenic role of TW independent of invasion in vivo and the therapeutic relevance of targeting TW in human GBM. Further, the role of TW driving POSTN expression and AKT signaling suggests actionable targets, which could be leveraged to mitigate the oncogenic effects of TW in GBM.

Digitizing omics profiles by divergence from a baseline.

Data collected from omics technologies have revealed pervasive heterogeneity and stochasticity of molecular states within and between phenotypes. A prominent example of such heterogeneity occurs between genome-wide mRNA, microRNA, and methylation profiles from one individual tumor to another, even within a cancer subtype. However, current methods in bioinformatics, such as detecting differentially expressed genes or CpG sites, are population-based and therefore do not effectively model intersample diversity. Here we introduce a unified theory to quantify sample-level heterogeneity that is applicable to a single omics profile. Specifically, we simplify an omics profile to a digital representation based on the omics profiles from a set of samples from a reference or baseline population (e.g., normal tissues). The state of any subprofile (e.g., expression vector for a subset of genes) is said to be "divergent" if it lies outside the estimated support of the baseline distribution and is consequently interpreted as "dysregulated" relative to that baseline. We focus on two cases: single features (e.g., individual genes) and distinguished subsets (e.g., regulatory pathways). Notably, since the divergence analysis is at the individual sample level, dysregulation can be analyzed probabilistically; for example, one can estimate the probability that a gene or pathway is divergent in some population. Finally, the reduction in complexity facilitates a more "personalized" and biologically interpretable analysis of variation, as illustrated by experiments involving tissue characterization, disease detection and progression, and disease-pathway associations.

Genome-wide pleiotropy analysis of neuropathological traits related to Alzheimer's disease.

BACKGROUND: Simultaneous consideration of two neuropathological traits related to Alzheimer's disease (AD) has not been attempted in a genome-wide association study. METHODS: We conducted genome-wide pleiotropy analyses using association summary statistics from the Beecham et al. study (PLoS Genet 10:e1004606, 2014) for AD-related neuropathological traits, including neuritic plaque (NP), neurofibrillary tangle (NFT), and cerebral amyloid angiopathy (CAA). Significant findings were further examined by expression quantitative trait locus and differentially expressed gene analyses in AD vs. control brains using gene expression data. RESULTS: Genome-wide significant pleiotropic associations were observed for the joint model of NP and NFT (NP + NFT) with the single-nucleotide polymorphism (SNP) rs34487851 upstream of C2orf40 (alias ECRG4, P = 2.4 × 10-8) and for the joint model of NFT and CAA (NFT + CAA) with the HDAC9 SNP rs79524815 (P = 1.1 × 10-8). Gene-based testing revealed study-wide significant associations (P ≤ 2.0 × 10-6) for the NFT + CAA outcome with adjacent genes TRAPPC12, TRAPPC12-AS1, and ADI1. Risk alleles of proxy SNPs for rs79524815 were associated with significantly lower expression of HDAC9 in the brain (P = 3.0 × 10-3), and HDAC9 was significantly downregulated in subjects with AD compared with control subjects in the prefrontal (P = 7.9 × 10-3) and visual (P = 5.6 × 10-4) cortices. CONCLUSIONS: Our findings suggest that pleiotropy analysis is a useful approach to identifying novel genetic associations with complex diseases and their endophenotypes. Functional studies are needed to determine whether ECRG4 or HDAC9 is plausible as a therapeutic target.

The building blocks of successful translation of proteomics to the clinic.

Recently, the first two multiplexed tests using selective reaction monitoring (SRM-MS) mass spectrometry have entered clinical practice. Despite different areas of indication, risk stratification in lung cancer and preterm birth, they share multiple steps in their development strategies. Here we review these strategies and their implications for successful translation of biomarkers to clinical practice. We believe that the identification of blood protein panels for the identification of disease phenotypes is now a reproducible and standard (albeit complex) process.

Comparative transcriptomic analysis of human placentae at term and preterm delivery.

Preterm birth affects 1 out of every 10 infants in the United States, resulting in substantial neonatal morbidity and mortality. Currently, there are few predictive markers and few treatment options to prevent preterm birth. A healthy, functioning placenta is essential to positive pregnancy outcomes. Previous studies have suggested that placental pathology may play a role in preterm birth etiology. Therefore, we tested the hypothesis that preterm placentae may exhibit unique transcriptomic signatures compared to term samples reflective of their abnormal biology leading to this adverse outcome. We aggregated publicly available placental villous microarray data to generate a preterm and term sample dataset (n = 133, 55 preterm placentae and 78 normal term placentae). We identified differentially expressed genes using the linear regression for microarray (LIMMA) package and identified perturbations in known biological networks using Differential Rank Conservation (DIRAC). We identified 129 significantly differentially expressed genes between term and preterm placenta with 96 genes upregulated and 33 genes downregulated (P-value <0.05). Significant changes in gene expression in molecular networks related to Tumor Protein 53 and phosphatidylinositol signaling were identified using DIRAC. We have aggregated a uniformly normalized transcriptomic dataset and have identified novel and established genes and pathways associated with developmental regulation of the placenta and potential preterm birth pathology. These analyses provide a community resource to integrate with other high-dimensional datasets for additional insights in normal placental development and its disruption.

Linkage, whole genome sequence, and biological data implicate variants in RAB10 in Alzheimer's disease resilience.

BACKGROUND: While age and the APOE ε4 allele are major risk factors for Alzheimer's disease (AD), a small percentage of individuals with these risk factors exhibit AD resilience by living well beyond 75 years of age without any clinical symptoms of cognitive decline. METHODS: We used over 200 "AD resilient" individuals and an innovative, pedigree-based approach to identify genetic variants that segregate with AD resilience. First, we performed linkage analyses in pedigrees with resilient individuals and a statistical excess of AD deaths. Second, we used whole genome sequences to identify candidate SNPs in significant linkage regions. Third, we replicated SNPs from the linkage peaks that reduced risk for AD in an independent dataset and in a gene-based test. Finally, we experimentally characterized replicated SNPs. RESULTS: Rs142787485 in RAB10 confers significant protection against AD (p value = 0.0184, odds ratio = 0.5853). Moreover, we replicated this association in an independent series of unrelated individuals (p value = 0.028, odds ratio = 0.69) and used a gene-based test to confirm a role for RAB10 variants in modifying AD risk (p value = 0.002). Experimentally, we demonstrated that knockdown of RAB10 resulted in a significant decrease in Aß42 (p value = 0.0003) and in the Aß42/Aß40 ratio (p value = 0.0001) in neuroblastoma cells. We also found that RAB10 expression is significantly elevated in human AD brains (p value = 0.04). CONCLUSIONS: Our results suggest that RAB10 could be a promising therapeutic target for AD prevention. In addition, our gene discovery approach can be expanded and adapted to other phenotypes, thus serving as a model for future efforts to identify rare variants for AD and other complex human diseases.

A wellness study of 108 individuals using personal, dense, dynamic data clouds.

Personal data for 108 individuals were collected during a 9-month period, including whole genome sequences; clinical tests, metabolomes, proteomes, and microbiomes at three time points; and daily activity tracking. Using all of these data, we generated a correlation network that revealed communities of related analytes associated with physiology and disease. Connectivity within analyte communities enabled the identification of known and candidate biomarkers (e.g., gamma-glutamyltyrosine was densely interconnected with clinical analytes for cardiometabolic disease). We calculated polygenic scores from genome-wide association studies (GWAS) for 127 traits and diseases, and used these to discover molecular correlates of polygenic risk (e.g., genetic risk for inflammatory bowel disease was negatively correlated with plasma cystine). Finally, behavioral coaching informed by personal data helped participants to improve clinical biomarkers. Our results show that measurement of personal data clouds over time can improve our understanding of health and disease, including early transitions to disease states.