Proteomic prediction of incident heart failure and its main subtypes.

AIM: To examine the ability of serum proteins in predicting future heart failure (HF) events, including HF with reduced or preserved ejection fraction (HFrEF or HFpEF), in relation to event time, and with or without considering established HF-associated clinical variables. METHODS AND RESULTS: In the prospective population-based Age, Gene/Environment Susceptibility Reykjavik Study (AGES-RS), 440 individuals developed HF after their first visit with a median follow-up of 5.45 years. Among them, 167 were diagnosed with HFrEF and 188 with HFpEF. A least absolute shrinkage and selection operator regression model with non-parametric bootstrap were used to select predictors from an analysis of 4782 serum proteins, and several pre-established clinical parameters linked to HF. A subset of 8-10 distinct or overlapping serum proteins predicted different future HF outcomes, and C-statistics were used to assess discrimination, revealing proteins combined with a C-index of 0.80 for all incident HF, 0.78 and 0.80 for incident HFpEF or HFrEF, respectively. In the AGES-RS, protein panels alone encompassed the risk contained in the clinical information and improved the performance characteristics of prediction models based on N-terminal pro-B-type natriuretic peptide and clinical risk factors. Finally, the protein predictors performed particularly well close to the time of an HF event, an outcome that was replicated in the Cardiovascular Health Study. CONCLUSION: A small number of circulating proteins accurately predicted future HF in the AGES-RS cohort of older adults, and they alone encompass the risk information found in a collection of clinical data. Incident HF events were predicted up to 8 years, with predictor performance significantly improving for events occurring less than 1 year ahead, a finding replicated in an external cohort study.

A proteogenomic signature of age-related macular degeneration in blood.

Age-related macular degeneration (AMD) is one of the most common causes of visual impairment in the elderly, with a complex and still poorly understood etiology. Whole-genome association studies have discovered 34 genomic regions associated with AMD. However, the genes and cognate proteins that mediate the risk, are largely unknown. In the current study, we integrate levels of 4782 human serum proteins with all genetic risk loci for AMD in a large population-based study of the elderly, revealing many proteins and pathways linked to the disease. Serum proteins are also found to reflect AMD severity independent of genetics and predict progression from early to advanced AMD after five years in this population. A two-sample Mendelian randomization study identifies several proteins that are causally related to the disease and are directionally consistent with the observational estimates. In this work, we present a robust and unique framework for elucidating the pathobiology of AMD.

A genome-wide association study of serum proteins reveals shared loci with common diseases.

With the growing number of genetic association studies, the genotype-phenotype atlas has become increasingly more complex, yet the functional consequences of most disease associated alleles is not understood. The measurement of protein level variation in solid tissues and biofluids integrated with genetic variants offers a path to deeper functional insights. Here we present a large-scale proteogenomic study in 5,368 individuals, revealing 4,035 independent associations between genetic variants and 2,091 serum proteins, of which 36% are previously unreported. The majority of both cis- and trans-acting genetic signals are unique for a single protein, although our results also highlight numerous highly pleiotropic genetic effects on protein levels and demonstrate that a protein's genetic association profile reflects certain characteristics of the protein, including its location in protein networks, tissue specificity and intolerance to loss of function mutations. Integrating protein measurements with deep phenotyping of the cohort, we observe substantial enrichment of phenotype associations for serum proteins regulated by established GWAS loci, and offer new insights into the interplay between genetics, serum protein levels and complex disease.

Coding and regulatory variants are associated with serum protein levels and disease.

Circulating proteins can be used to diagnose and predict disease-related outcomes. A deep serum proteome survey recently revealed close associations between serum protein networks and common disease. In the current study, 54,469 low-frequency and common exome-array variants were compared to 4782 protein measurements in the serum of 5343 individuals from the AGES Reykjavik cohort. This analysis identifies a large number of serum proteins with genetic signatures overlapping those of many diseases. More specifically, using a study-wide significance threshold, we find that 2021 independent exome array variants are associated with serum levels of 1942 proteins. These variants reside in genetic loci shared by hundreds of complex disease traits, highlighting serum proteins' emerging role as biomarkers and potential causative agents of a wide range of diseases.

Serum levels of ACE2 are higher in patients with obesity and diabetes.

OBJECTIVE: As severity of outcome in COVID-19 is disproportionately higher among individuals with obesity, smokers, patients with hypertension, kidney disease, chronic pulmonary disease, coronary heart disease (CHD), and/or type 2 diabetes (T2D), serum levels of ACE2, the cellular entry point for the coronavirus SARS-CoV-2, were examined in these high-risk groups. METHODS: Associations of ACE2 levels to smokers and patients with hypertension, T2D, obesity, CHD, or COPD were investigated in a single center population-based study of 5457 Icelanders from the Age, Gene/Environment Susceptibility Reykjavík Study (AGES-RS) of the elderly (mean age 75 ± 6 years), using multiple linear regression analysis. RESULTS: Serum levels of ACE2 were higher in smokers and individuals with T2D and/or obesity while they were unaffected in the other patient groups. CONCLUSION: ACE2 levels are higher in some patient groups with comorbidities linked to COVID-19 including obesity and T2D and as such may have an emerging role as a circulating biomarker for severity of outcome in the disease.

It's in Our Blood: A Glimpse of Personalized Medicine.

Recent advances in protein profiling technology has facilitated simultaneous measurement of thousands of proteins in large population studies, exposing the depth and complexity of the plasma and serum proteomes. This revealed that proteins in circulation were organized into regulatory modules under genetic control and closely associated with current and future common diseases. Unlike networks in solid tissues, serum protein networks comprise members synthesized across different tissues of the body. Genetic analysis reveals that this cross-tissue regulation of the serum proteome participates in systemic homeostasis and mirrors the global disease state of individuals. Here, we discuss how application of this information in routine clinical evaluations may transform the future practice of medicine.

ACE2 levels are altered in comorbidities linked to severe outcome in COVID-19.

AIMS: Severity of outcome in COVID-19 is disproportionately higher among the obese, males, smokers, those suffering from hypertension, kidney disease, coronary heart disease (CHD) and/or type 2 diabetes (T2D). We examined if serum levels of ACE2, the cellular entry point for the coronavirus SARS-CoV-2, were altered in these high-risk groups. METHODS: Associations of serum ACE2 levels to hypertension, T2D, obesity, CHD, smokers and males in a single center population-based study of 5457 Icelanders from the Age, Gene/Environment Susceptibility Reykjavik Study (AGES-RS) of the elderly (mean age 75+/-6 years). RESULTS: Smokers, males, and individuals with T2D or obesity have altered serum levels of ACE2 that may influence productive infection of SARS-CoV-2 in these high-risk groups. CONCLUSION: ACE2 levels are upregulated in some patient groups with comorbidities linked to COVID-19 and as such may have an emerging role as a circulating biomarker for severity of outcome in COVID-19.

Circulating Protein Signatures and Causal Candidates for Type 2 Diabetes.

The increasing prevalence of type 2 diabetes poses a major challenge to societies worldwide. Blood-based factors like serum proteins are in contact with every organ in the body to mediate global homeostasis and may thus directly regulate complex processes such as aging and the development of common chronic diseases. We applied a data-driven proteomics approach, measuring serum levels of 4,137 proteins in 5,438 elderly Icelanders, and identified 536 proteins associated with prevalent and/or incident type 2 diabetes. We validated a subset of the observed associations in an independent case-control study of type 2 diabetes. These protein associations provide novel biological insights into the molecular mechanisms that are dysregulated prior to and following the onset of type 2 diabetes and can be detected in serum. A bidirectional two-sample Mendelian randomization analysis indicated that serum changes of at least 23 proteins are downstream of the disease or its genetic liability, while 15 proteins were supported as having a causal role in type 2 diabetes.

Tau overexpression exacerbates neuropathology after repeated mild head impacts in male mice.

Repeated mild traumatic brain injury (rmTBI) can lead to development of chronic traumatic encephalopathy (CTE), which is characterized by progressive neurodegeneration with presence of white matter damage, gliosis and hyper-phosphorylated tau. While animal models of rmTBI have been documented, few characterize the molecular pathogenesis and expression profiles of relevant injured brain regions. Additionally, while the usage of transgenic tau mice in rmTBI is prevalent, the effects of tau on pathological outcomes has not been well studied. Here we characterized a 42-impact closed-head rmTBI paradigm on 3-4 month old male C57BL/6 (WT) and Tau-overexpressing mice (Tau58.4). This injury paradigm resulted in chronic gliosis, T-cell infiltration, and demyelination of the optic nerve and associated white matter tracts at 1-month post-injury. At 3-months post-injury, Tau58.4 mice showed progressive neuroinflammation and neurodegeneration in multiple brain regions compared to WT mice. Corresponding to histopathology, RNAseq of the optic nerve tract at 1-month post-injury showed significant upregulation of inflammatory pathways and downregulation of myelin synthetic pathways in both genotypes. However, Tau58.4 mice showed additional changes in neurite development, protein processing, and cell stress. Comparisons with published transcriptomes of human Alzheimer's Disease and CTE revealed common signatures including neuroinflammation and downregulation of protein phosphatases. We next investigated the demyelination and T-cell infiltration phenotypes to determine whether these offer potential avenues for therapeutic intervention. Tau58.4 mice were treated with the histamine H3 receptor antagonist GSK239512 for 1-month post-injury to promote remyelination of white matter lesions. This restored myelin gene expression to sham levels but failed to repair the histopathologic lesions. Likewise, injured T-cell-deficient Rag2/Il2rg (R2G2) mice also showed evidence for inflammation and loss of myelin. However, unlike immune-competent mice, R2G2 mice had altered myeloid cell gene expression and fewer demyelinated lesions. Together this data shows that rmTBI leads to chronic white matter inflammatory demyelination and axonal loss exacerbated by human tau overexpression but suggests that immune-suppression and remyelination alone are insufficient to reverse damage.

Co-regulatory networks of human serum proteins link genetics to disease.

Proteins circulating in the blood are critical for age-related disease processes; however, the serum proteome has remained largely unexplored. To this end, 4137 proteins covering most predicted extracellular proteins were measured in the serum of 5457 Icelanders over 65 years of age. Pairwise correlation between proteins as they varied across individuals revealed 27 different network modules of serum proteins, many of which were associated with cardiovascular and metabolic disease states, as well as overall survival. The protein modules were controlled by cis- and trans-acting genetic variants, which in many cases were also associated with complex disease. This revealed co-regulated groups of circulating proteins that incorporated regulatory control between tissues and demonstrated close relationships to past, current, and future disease states.

Common dysregulation network in the human prefrontal cortex underlies two neurodegenerative diseases.

Using expression profiles from postmortem prefrontal cortex samples of 624 dementia patients and non-demented controls, we investigated global disruptions in the co-regulation of genes in two neurodegenerative diseases, late-onset Alzheimer's disease (AD) and Huntington's disease (HD). We identified networks of differentially co-expressed (DC) gene pairs that either gained or lost correlation in disease cases relative to the control group, with the former dominant for both AD and HD and both patterns replicating in independent human cohorts of AD and aging. When aligning networks of DC patterns and physical interactions, we identified a 242-gene subnetwork enriched for independent AD/HD signatures. This subnetwork revealed a surprising dichotomy of gained/lost correlations among two inter-connected processes, chromatin organization and neural differentiation, and included DNA methyltransferases, DNMT1 and DNMT3A, of which we predicted the former but not latter as a key regulator. To validate the inter-connection of these two processes and our key regulator prediction, we generated two brain-specific knockout (KO) mice and show that Dnmt1 KO signature significantly overlaps with the subnetwork (P = 3.1 × 10(-12)), while Dnmt3a KO signature does not (P = 0.017).

Integrated systems approach identifies genetic nodes and networks in late-onset Alzheimer's disease.

The genetics of complex disease produce alterations in the molecular interactions of cellular pathways whose collective effect may become clear through the organized structure of molecular networks. To characterize molecular systems associated with late-onset Alzheimer's disease (LOAD), we constructed gene-regulatory networks in 1,647 postmortem brain tissues from LOAD patients and nondemented subjects, and we demonstrate that LOAD reconfigures specific portions of the molecular interaction structure. Through an integrative network-based approach, we rank-ordered these network structures for relevance to LOAD pathology, highlighting an immune- and microglia-specific module that is dominated by genes involved in pathogen phagocytosis, contains TYROBP as a key regulator, and is upregulated in LOAD. Mouse microglia cells overexpressing intact or truncated TYROBP revealed expression changes that significantly overlapped the human brain TYROBP network. Thus the causal network structure is a useful predictor of response to gene perturbations and presents a framework to test models of disease mechanisms underlying LOAD.

A comprehensive characterization of genome-wide copy number aberrations in colorectal cancer reveals novel oncogenes and patterns of alterations.

To develop a comprehensive overview of copy number aberrations (CNAs) in stage-II/III colorectal cancer (CRC), we characterized 302 tumors from the PETACC-3 clinical trial. Microsatellite-stable (MSS) samples (n = 269) had 66 minimal common CNA regions, with frequent gains on 20 q (72.5%), 7 (41.8%), 8 q (33.1%) and 13 q (51.0%) and losses on 18 (58.6%), 4 q (26%) and 21 q (21.6%). MSS tumors have significantly more CNAs than microsatellite-instable (MSI) tumors: within the MSI tumors a novel deletion of the tumor suppressor WWOX at 16 q23.1 was identified (p<0.01). Focal aberrations identified by the GISTIC method confirmed amplifications of oncogenes including EGFR, ERBB2, CCND1, MET, and MYC, and deletions of tumor suppressors including TP53, APC, and SMAD4, and gene expression was highly concordant with copy number aberration for these genes. Novel amplicons included putative oncogenes such as WNK1 and HNF4A, which also showed high concordance between copy number and expression. Survival analysis associated a specific patient segment featured by chromosome 20 q gains to an improved overall survival, which might be due to higher expression of genes such as EEF1B2 and PTK6. The CNA clustering also grouped tumors characterized by a poor prognosis BRAF-mutant-like signature derived from mRNA data from this cohort. We further revealed non-random correlation between CNAs among unlinked loci, including positive correlation between 20 q gain and 8 q gain, and 20 q gain and chromosome 18 loss, consistent with co-selection of these CNAs. These results reinforce the non-random nature of somatic CNAs in stage-II/III CRC and highlight loci and genes that may play an important role in driving the development and outcome of this disease.

Inferring causal genomic alterations in breast cancer using gene expression data.

BACKGROUND: One of the primary objectives in cancer research is to identify causal genomic alterations, such as somatic copy number variation (CNV) and somatic mutations, during tumor development. Many valuable studies lack genomic data to detect CNV; therefore, methods that are able to infer CNVs from gene expression data would help maximize the value of these studies. RESULTS: We developed a framework for identifying recurrent regions of CNV and distinguishing the cancer driver genes from the passenger genes in the regions. By inferring CNV regions across many datasets we were able to identify 109 recurrent amplified/deleted CNV regions. Many of these regions are enriched for genes involved in many important processes associated with tumorigenesis and cancer progression. Genes in these recurrent CNV regions were then examined in the context of gene regulatory networks to prioritize putative cancer driver genes. The cancer driver genes uncovered by the framework include not only well-known oncogenes but also a number of novel cancer susceptibility genes validated via siRNA experiments. CONCLUSIONS: To our knowledge, this is the first effort to systematically identify and validate drivers for expression based CNV regions in breast cancer. The framework where the wavelet analysis of copy number alteration based on expression coupled with the gene regulatory network analysis, provides a blueprint for leveraging genomic data to identify key regulatory components and gene targets. This integrative approach can be applied to many other large-scale gene expression studies and other novel types of cancer data such as next-generation sequencing based expression (RNA-Seq) as well as CNV data.

Predictive genes in adjacent normal tissue are preferentially altered by sCNV during tumorigenesis in liver cancer and may rate limiting.

BACKGROUND: In hepatocellular carcinoma (HCC) genes predictive of survival have been found in both adjacent normal (AN) and tumor (TU) tissues. The relationships between these two sets of predictive genes and the general process of tumorigenesis and disease progression remains unclear. METHODOLOGY/PRINCIPAL FINDINGS: Here we have investigated HCC tumorigenesis by comparing gene expression, DNA copy number variation and survival using ∼250 AN and TU samples representing, respectively, the pre-cancer state, and the result of tumorigenesis. Genes that participate in tumorigenesis were defined using a gene-gene correlation meta-analysis procedure that compared AN versus TU tissues. Genes predictive of survival in AN (AN-survival genes) were found to be enriched in the differential gene-gene correlation gene set indicating that they directly participate in the process of tumorigenesis. Additionally the AN-survival genes were mostly not predictive after tumorigenesis in TU tissue and this transition was associated with and could largely be explained by the effect of somatic DNA copy number variation (sCNV) in cis and in trans. The data was consistent with the variance of AN-survival genes being rate-limiting steps in tumorigenesis and this was confirmed using a treatment that promotes HCC tumorigenesis that selectively altered AN-survival genes and genes differentially correlated between AN and TU. CONCLUSIONS/SIGNIFICANCE: This suggests that the process of tumor evolution involves rate-limiting steps related to the background from which the tumor evolved where these were frequently predictive of clinical outcome. Additionally treatments that alter the likelihood of tumorigenesis occurring may act by altering AN-survival genes, suggesting that the process can be manipulated. Further sCNV explains a substantial fraction of tumor specific expression and may therefore be a causal driver of tumor evolution in HCC and perhaps many solid tumor types.

microRNA-122 as a regulator of mitochondrial metabolic gene network in hepatocellular carcinoma.

Tumorigenesis involves multistep genetic alterations. To elucidate the microRNA (miRNA)-gene interaction network in carcinogenesis, we examined their genome-wide expression profiles in 96 pairs of tumor/non-tumor tissues from hepatocellular carcinoma (HCC). Comprehensive analysis of the coordinate expression of miRNAs and mRNAs reveals that miR-122 is under-expressed in HCC and that increased expression of miR-122 seed-matched genes leads to a loss of mitochondrial metabolic function. Furthermore, the miR-122 secondary targets, which decrease in expression, are good prognostic markers for HCC. Transcriptome profiling data from additional 180 HCC and 40 liver cirrhotic patients in the same cohort were used to confirm the anti-correlation of miR-122 primary and secondary target gene sets. The HCC findings can be recapitulated in mouse liver by silencing miR-122 with antagomir treatment followed by gene-expression microarray analysis. In vitro miR-122 data further provided a direct link between induction of miR-122-controlled genes and impairment of mitochondrial metabolism. In conclusion, miR-122 regulates mitochondrial metabolism and its loss may be detrimental to sustaining critical liver function and contribute to morbidity and mortality of liver cancer patients.

Complexity in common diseases: big biology for all.

The Commentary of Mills and Sykes in Science Translational Medicine presented their thesis on the advent of high-throughput technologies and the dangers they may represent for the future of biomedical research. In response, we argue that true progress on the diagnosis and treatment of common human diseases will require the advent of big biology and its deep integration with focused research as practiced in both academic and industrial institutions.

Characterizing dynamic changes in the human blood transcriptional network.

Gene expression data generated systematically in a given system over multiple time points provides a source of perturbation that can be leveraged to infer causal relationships among genes explaining network changes. Previously, we showed that food intake has a large impact on blood gene expression patterns and that these responses, either in terms of gene expression level or gene-gene connectivity, are strongly associated with metabolic diseases. In this study, we explored which genes drive the changes of gene expression patterns in response to time and food intake. We applied the Granger causality test and the dynamic Bayesian network to gene expression data generated from blood samples collected at multiple time points during the course of a day. The simulation result shows that combining many short time series together is as powerful to infer Granger causality as using a single long time series. Using the Granger causality test, we identified genes that were supported as the most likely causal candidates for the coordinated temporal changes in the network. These results show that PER1 is a key regulator of the blood transcriptional network, in which multiple biological processes are under circadian rhythm regulation. The fasted and fed dynamic Bayesian networks showed that over 72% of dynamic connections are self links. Finally, we show that different processes such as inflammation and lipid metabolism, which are disconnected in the static network, become dynamically linked in response to food intake, which would suggest that increasing nutritional load leads to coordinate regulation of these biological processes. In conclusion, our results suggest that food intake has a profound impact on the dynamic co-regulation of multiple biological processes, such as metabolism, immune response, apoptosis and circadian rhythm. The results could have broader implications for the design of studies of disease association and drug response in clinical trials.

The effect of food intake on gene expression in human peripheral blood.

Human gene expression traits have been shown to be dependent on gender, age and time of day in blood and other tissues. However, other factors that may impact gene expression have not been systematically explored. For example, in studies linking blood gene expression to obesity related traits, whether the fasted or fed state will be the most informative is an open question. Here, we employed a two-arm cross-over design to perform a genome-wide survey of gene expression in human peripheral blood to address explicitly this type of question. We were able to distinguish expression changes due to individual and time-specific effects from those due to food intake. We demonstrate that the transcriptional response to food intake is robust by constructing a classifier from the gene expression traits with >90% accuracy classifying individuals as being in the fasted or fed state. Gene expression traits that were best able to discriminate the fasted and fed states were more heritable and achieved greater coherence with respect to pathways associated with metabolic traits. The connectivity structure among gene expression traits was explored in the context of coexpression networks. Changes in the connectivity structure were observed between the fasted and fed states. We demonstrate that differential expression and differential connectivity are two complementary ways to characterize changes between fasted and fed states. Both gene sets were significantly enriched for genes associated with obesity related traits. Our results suggest that the pair of fasted/fed blood expression profiles provide more comprehensive information about an individual's metabolic states.

Validation of candidate causal genes for obesity that affect shared metabolic pathways and networks.

A principal task in dissecting the genetics of complex traits is to identify causal genes for disease phenotypes. We previously developed a method to infer causal relationships among genes through the integration of DNA variation, gene transcription and phenotypic information. Here we have validated our method through the characterization of transgenic and knockout mouse models of genes predicted to be causal for abdominal obesity. Perturbation of eight out of the nine genes, with Gas7, Me1 and Gpx3 being newly confirmed, resulted in significant changes in obesity-related traits. Liver expression signatures revealed alterations in common metabolic pathways and networks contributing to abdominal obesity and overlapped with a macrophage-enriched metabolic network module that is highly associated with metabolic traits in mice and humans. Integration of gene expression in the design and analysis of traditional F(2) intercross studies allows high-confidence prediction of causal genes and identification of pathways and networks involved.