Integration of transcriptomes of senescent cell models with multi-tissue patient samples reveals reduced COL6A3 as an inducer of senescence.

Senescent cells are a major contributor to age-dependent cardiovascular tissue dysfunction, but knowledge of their in vivo cell markers and tissue context is lacking. To reveal tissue-relevant senescence biology, we integrate the transcriptomes of 10 experimental senescence cell models with a 224 multi-tissue gene co-expression network based on RNA-seq data of seven tissues biopsies from ∼600 coronary artery disease (CAD) patients. We identify 56 senescence-associated modules, many enriched in CAD GWAS genes and correlated with cardiometabolic traits-which supports universality of senescence gene programs across tissues and in CAD. Cross-tissue network analyses reveal 86 candidate senescence-associated secretory phenotype (SASP) factors, including COL6A3. Experimental knockdown of COL6A3 induces transcriptional changes that overlap the majority of the experimental senescence models, with cell-cycle arrest linked to modulation of DREAM complex-targeted genes. We provide a transcriptomic resource for cellular senescence and identify candidate biomarkers, SASP factors, and potential drivers of senescence in human tissues.

Molecular differences in brain regional vulnerability to aging between males and females.

Background: Aging-related cognitive decline is associated with brain structural changes and synaptic loss. However, the molecular mechanisms of cognitive decline during normal aging remain elusive. Results: Using the GTEx transcriptomic data from 13 brain regions, we identified aging-associated molecular alterations and cell-type compositions in males and females. We further constructed gene co-expression networks and identified aging-associated modules and key regulators shared by both sexes or specific to males or females. A few brain regions such as the hippocampus and the hypothalamus show specific vulnerability in males, while the cerebellar hemisphere and the anterior cingulate cortex regions manifest greater vulnerability in females than in males. Immune response genes are positively correlated with age, whereas those involved in neurogenesis are negatively correlated with age. Aging-associated genes identified in the hippocampus and the frontal cortex are significantly enriched for gene signatures implicated in Alzheimer's disease (AD) pathogenesis. In the hippocampus, a male-specific co-expression module is driven by key synaptic signaling regulators including VSNL1, INA, CHN1 and KCNH1; while in the cortex, a female-specific module is associated with neuron projection morphogenesis, which is driven by key regulators including SRPK2, REPS2 and FXYD1. In the cerebellar hemisphere, a myelination-associated module shared by males and females is driven by key regulators such as MOG, ENPP2, MYRF, ANLN, MAG and PLP1, which have been implicated in the development of AD and other neurodegenerative diseases. Conclusions: This integrative network biology study systematically identifies molecular signatures and networks underlying brain regional vulnerability to aging in males and females. The findings pave the way for understanding the molecular mechanisms of gender differences in developing neurodegenerative diseases such as AD.

Multiscale network analysis identifies potential receptors for SARS-CoV-2 and reveals their tissue-specific and age-dependent expression.

The coronavirus disease 2019 (COVID-19) pandemic has affected tens of millions of individuals and caused hundreds of thousands of deaths worldwide. Here, we present a comprehensive, multiscale network analysis of the transcriptional response to the virus. In particular, we focused on key regulators, cell receptors, and host processes that were hijacked by the virus for its advantage. ACE2-controlled processes involved CD300e (a TYROBP receptor) as a key regulator and the activation of IL-2 pro-inflammatory cytokine signaling. We further investigated the age dependency of such receptors in different tissues. In summary, this study provides novel insights into the gene regulatory organization during the SARS-CoV-2 infection and the tissue-specific, age-dependent expression of the cell receptors involved in COVID-19.

The landscape of human tissue and cell type specific expression and co-regulation of senescence genes.

BACKGROUND: Cellular senescence is a complex stress response that impacts cellular function and organismal health. Multiple developmental and environmental factors, such as intrinsic cellular cues, radiation, oxidative stress, oncogenes, and protein accumulation, activate genes and pathways that can lead to senescence. Enormous efforts have been made to identify and characterize senescence genes (SnGs) in stress and disease systems. However, the prevalence of senescent cells in healthy human tissues and the global SnG expression signature in different cell types are poorly understood. METHODS: This study performed an integrative gene network analysis of bulk and single-cell RNA-seq data in non-diseased human tissues to investigate SnG co-expression signatures and their cell-type specificity. RESULTS: Through a comprehensive transcriptomic network analysis of 50 human tissues in the Genotype-Tissue Expression Project (GTEx) cohort, we identified SnG-enriched gene modules, characterized SnG co-expression patterns, and constructed aggregated SnG networks across primary tissues of the human body. Our network approaches identified 51 SnGs highly conserved across the human tissues, including CDKN1A (p21)-centered regulators that control cell cycle progression and the senescence-associated secretory phenotype (SASP). The SnG-enriched modules showed remarkable cell-type specificity, especially in fibroblasts, endothelial cells, and immune cells. Further analyses of single-cell RNA-seq and spatial transcriptomic data independently validated the cell-type specific SnG signatures predicted by the network analysis. CONCLUSIONS: This study systematically revealed the co-regulated organizations and cell type specificity of SnGs in major human tissues, which can serve as a blueprint for future studies to map senescent cells and their cellular interactions in human tissues.

Transcriptomic Changes Highly Similar to Alzheimer's Disease Are Observed in a Subpopulation of Individuals During Normal Brain Aging.

Aging is a major risk factor for late-onset Alzheimer's disease (LOAD). How aging contributes to the development of LOAD remains elusive. In this study, we examined multiple large-scale transcriptomic datasets from both normal aging and LOAD brains to understand the molecular interconnection between aging and LOAD. We found that shared gene expression changes between aging and LOAD are mostly seen in the hippocampal and several cortical regions. In the hippocampus, the expression of phosphoprotein, alternative splicing and cytoskeleton genes are commonly changed in both aging and AD, while synapse, ion transport, and synaptic vesicle genes are commonly down-regulated. Aging-specific changes are associated with acetylation and methylation, while LOAD-specific changes are more related to glycoprotein (both up- and down-regulations), inflammatory response (up-regulation), myelin sheath and lipoprotein (down-regulation). We also found that normal aging brain transcriptomes from relatively young donors (45-70 years old) clustered into several subgroups and some subgroups showed gene expression changes highly similar to those seen in LOAD brains. Using brain transcriptomic datasets from another cohort of older individuals (>70 years), we found that samples from cognitively normal older individuals clustered with the "healthy aging" subgroup while AD samples mainly clustered with the "AD similar" subgroups. This may imply that individuals in the healthy aging subgroup will likely remain cognitively normal when they become older and vice versa. In summary, our results suggest that on the transcriptome level, aging and LOAD have strong interconnections in some brain regions in a subpopulation of cognitively normal aging individuals. This supports the theory that the initiation of LOAD occurs decades earlier than the manifestation of clinical phenotype and it may be essential to closely study the "normal brain aging" to identify the very early molecular events that may lead to LOAD development.

A comprehensive integrated post-GWAS analysis of Type 1 diabetes reveals enhancer-based immune dysregulation.

Type 1 diabetes (T1D) is an organ-specific autoimmune disease, whereby immune cell-mediated killing leads to loss of the insulin-producing ß cells in the pancreas. Genome-wide association studies (GWAS) have identified over 200 genetic variants associated with risk for T1D. The majority of the GWAS risk variants reside in the non-coding regions of the genome, suggesting that gene regulatory changes substantially contribute to T1D. However, identification of causal regulatory variants associated with T1D risk and their affected genes is challenging due to incomplete knowledge of non-coding regulatory elements and the cellular states and processes in which they function. Here, we performed a comprehensive integrated post-GWAS analysis of T1D to identify functional regulatory variants in enhancers and their cognate target genes. Starting with 1,817 candidate T1D SNPs defined from the GWAS catalog and LDlink databases, we conducted functional annotation analysis using genomic data from various public databases. These include 1) Roadmap Epigenomics, ENCODE, and RegulomeDB for epigenome data; 2) GTEx for tissue-specific gene expression and expression quantitative trait loci data; and 3) lncRNASNP2 for long non-coding RNA data. Our results indicated a prevalent enhancer-based immune dysregulation in T1D pathogenesis. We identified 26 high-probability causal enhancer SNPs associated with T1D, and 64 predicted target genes. The majority of the target genes play major roles in antigen presentation and immune response and are regulated through complex transcriptional regulatory circuits, including those in HLA (6p21) and non-HLA (16p11.2) loci. These candidate causal enhancer SNPs are supported by strong evidence and warrant functional follow-up studies.

A Compendium of Age-Related PheWAS and GWAS Traits for Human Genetic Association Studies, Their Networks and Genetic Correlations.

The rich data from the genome-wide association studies (GWAS) and phenome-wide association studies (PheWAS) offer an unprecedented opportunity to identify the biological underpinnings of age-related disease (ARD) risk and multimorbidity. Surprisingly, however, a comprehensive list of ARDs remains unavailable due to the lack of a clear definition and selection criteria. We developed a method to identify ARDs and to provide a compendium of ARDs for genetic association studies. Querying 1,358 electronic medical record-derived traits, we first defined ARDs and age-related traits (ARTs) based on their prevalence profiles, requiring a unimodal distribution that shows an increasing prevalence after the age of 40 years, and which reaches a maximum peak at 60 years of age or later. As a result, we identified a list of 463 ARDs and ARTs in the GWAS and PheWAS catalogs. We next translated the ARDs and ARTs to their respective 276 Medical Subject Headings diseases and 45 anatomy terms. The most abundant disease categories are neoplasms (48 terms), cardiovascular diseases (44 terms), and nervous system diseases (27 terms). Employing data from a human symptoms-disease network, we found 6 symptom-shared disease groups, representing cancers, heart diseases, brain diseases, joint diseases, eye diseases, and mixed diseases. Lastly, by overlaying our ARD and ART list with genetic correlation data from the UK Biobank, we found 54 phenotypes in 2 clusters with high genetic correlations. Our compendium of ARD and ART is a highly useful resource, with broad applicability for studies of the genetics of aging, ARD, and multimorbidity.

Age-related telomere attrition causes aberrant gene expression in sub-telomeric regions.

Telomere attrition has been proposed as a biomarker and causal factor in aging. In addition to causing cellular senescence and apoptosis, telomere shortening has been found to affect gene expression in subtelomeric regions. Here, we analyzed the distribution of age-related differentially expressed genes from the GTEx RNA sequencing database of 54 tissue types from 979 human subjects and found significantly more upregulated than downregulated genes in subtelomeric regions as compared to the genome-wide average. Our data demonstrate spatial relationships between telomeres and gene expression in aging.

An integrative multiomic network model links lipid metabolism to glucose regulation in coronary artery disease.

Elevated plasma cholesterol and type 2 diabetes (T2D) are associated with coronary artery disease (CAD). Individuals treated with cholesterol-lowering statins have increased T2D risk, while individuals with hypercholesterolemia have reduced T2D risk. We explore the relationship between lipid and glucose control by constructing network models from the STARNET study with sequencing data from seven cardiometabolic tissues obtained from CAD patients during coronary artery by-pass grafting surgery. By integrating gene expression, genotype, metabolomic, and clinical data, we identify a glucose and lipid determining (GLD) regulatory network showing inverse relationships with lipid and glucose traits. Master regulators of the GLD network also impact lipid and glucose levels in inverse directions. Experimental inhibition of one of the GLD network master regulators, lanosterol synthase (LSS), in mice confirms the inverse relationships to glucose and lipid levels as predicted by our model and provides mechanistic insights.

Improved Human Age Prediction by Using Gene Expression Profiles From Multiple Tissues.

Studying transcriptome chronological change from tissues across the whole body can provide valuable information for understanding aging and longevity. Although there has been research on the effect of single-tissue transcriptomes on human aging or aging in mice across multiple tissues, the study of human body-wide multi-tissue transcriptomes on aging is not yet available. In this study, we propose a quantitative model to predict human age by using gene expression data from 46 tissues generated by the Genotype-Tissue Expression (GTEx) project. Specifically, the biological age of a person is first predicted via the gene expression profile of a single tissue. Then, we combine the gene expression profiles from two tissues and compare the predictive accuracy between single and two tissues. The best performance as measured by the root-mean-square error is 3.92 years for single tissue (pituitary), which deceased to 3.6 years when we combined two tissues (pituitary and muscle) together. Different tissues have different potential in predicting chronological age. The prediction accuracy is improved by combining multiple tissues, supporting that aging is a systemic process involving multiple tissues across the human body.

Multiscale causal networks identify VGF as a key regulator of Alzheimer's disease.

Though discovered over 100 years ago, the molecular foundation of sporadic Alzheimer's disease (AD) remains elusive. To better characterize the complex nature of AD, we constructed multiscale causal networks on a large human AD multi-omics dataset, integrating clinical features of AD, DNA variation, and gene- and protein-expression. These probabilistic causal models enabled detection, prioritization and replication of high-confidence master regulators of AD-associated networks, including the top predicted regulator, VGF. Overexpression of neuropeptide precursor VGF in 5xFAD mice partially rescued beta-amyloid-mediated memory impairment and neuropathology. Molecular validation of network predictions downstream of VGF was also achieved in this AD model, with significant enrichment for homologous genes identified as differentially expressed in 5xFAD brains overexpressing VGF. Our findings support a causal role for VGF in protecting against AD pathogenesis and progression.

Transcriptome analysis reveals the difference between "healthy" and "common" aging and their connection with age-related diseases.

A key goal of aging research was to understand mechanisms underlying healthy aging and develop methods to promote the human healthspan. One approach is to identify gene regulations unique to healthy aging compared with aging in the general population (i.e., "common" aging). Here, we leveraged Genotype-Tissue Expression (GTEx) project data to investigate "healthy" and "common" aging gene expression regulations at a tissue level in humans and their interconnection with diseases. Using GTEx donors' disease annotations, we defined a "healthy" aging cohort for each tissue. We then compared the age-associated genes derived from this cohort with age-associated genes from the "common" aging cohort which included all GTEx donors; we also compared the "healthy" and "common" aging gene expressions with various disease-associated gene expressions to elucidate the relationships among "healthy," "common" aging and disease. Our analyses showed that 1. GTEx "healthy" and "common" aging shared a large number of gene regulations; 2. Despite the substantial commonality, "healthy" and "common" aging genes also showed distinct function enrichment, and "common" aging genes had a higher enrichment for disease genes; 3. Disease-associated gene regulations were overall different from aging gene regulations. However, for genes regulated by both, their regulation directions were largely consistent, implying some aging processes could increase the susceptibility to disease development; and 4. Possible protective mechanisms were associated with some "healthy" aging gene regulations. In summary, our work highlights several unique features of GTEx "healthy" aging program. This new knowledge could potentially be used to develop interventions to promote the human healthspan.

Human geroprotector discovery by targeting the converging subnetworks of aging and age-related diseases.

A key goal of geroscience research is to identify effective interventions to extend human healthspan, the years of healthy life. Currently, majority of the geroprotectors are found by screening compounds in model organisms; whether they will be effective in humans is largely unknown. Here we present a new strategy called ANDRU (aging network based drug discovery) to help the discovery of human geroprotectors. It first identifies human aging subnetworks that putatively function at the interface between aging and age-related diseases; it then screens for pharmacological interventions that may "reverse" the age-associated transcriptional changes occurred in these subnetworks. We applied ANDRU to human adipose gene expression data from the Genotype Tissue Expression (GTEx) project. For the top 31 identified compounds, 19 of them showed at least some evidence supporting their function in improving metabolic traits or lifespan, which include type 2 diabetes drugs such as pioglitazone. As the query aging genes were refined to the ones with more intimate links to diseases, ANDRU identified more meaningful drug hits than the general approach without considering the underlying network structures. In summary, ANDRU represents a promising human data-driven strategy that may speed up the discovery of interventions to extend human healthspan.

A probabilistic multi-omics data matching method for detecting sample errors in integrative analysis.

BACKGROUND: Data errors, including sample swapping and mis-labeling, are inevitable in the process of large-scale omics data generation. Data errors need to be identified and corrected before integrative data analyses where different types of data are merged on the basis of the annotated labels. Data with labeling errors dampen true biological signals. More importantly, data analysis with sample errors could lead to wrong scientific conclusions. We developed a robust probabilistic multi-omics data matching procedure, proMODMatcher, to curate data and identify and correct data annotation and errors in large databases. RESULTS: Application to simulated datasets suggests that proMODMatcher achieved robust statistical power even when the number of cis-associations was small and/or the number of samples was large. Application of our proMODMatcher to multi-omics datasets in The Cancer Genome Atlas and International Cancer Genome Consortium identified sample errors in multiple cancer datasets. Our procedure was not only able to identify sample-labeling errors but also to unambiguously identify the source of the errors. Our results demonstrate that these errors should be identified and corrected before integrative analysis. CONCLUSIONS: Our results indicate that sample-labeling errors were common in large multi-omics datasets. These errors should be corrected before integrative analysis.

Temporal genetic association and temporal genetic causality methods for dissecting complex networks.

A large amount of panomic data has been generated in populations for understanding causal relationships in complex biological systems. Both genetic and temporal models can be used to establish causal relationships among molecular, cellular, or phenotypical traits, but with limitations. To fully utilize high-dimension temporal and genetic data, we develop a multivariate polynomial temporal genetic association (MPTGA) approach for detecting temporal genetic loci (teQTLs) of quantitative traits monitored over time in a population and a temporal genetic causality test (TGCT) for inferring causal relationships between traits linked to the locus. We apply MPTGA and TGCT to simulated data sets and a yeast F2 population in response to rapamycin, and demonstrate increased power to detect teQTLs. We identify a teQTL hotspot locus interacting with rapamycin treatment, infer putative causal regulators of the teQTL hotspot, and experimentally validate RRD1 as the causal regulator for this teQTL hotspot.

A next generation sequencing based approach to identify extracellular vesicle mediated mRNA transfers between cells.

BACKGROUND: Exosomes and other extracellular vesicles (EVs) have emerged as an important mechanism of cell-to-cell communication. However, previous studies either did not fully resolve what genetic materials were shuttled by exosomes or only focused on a specific set of miRNAs and mRNAs. A more systematic method is required to identify the genetic materials that are potentially transferred during cell-to-cell communication through EVs in an unbiased manner. RESULTS: In this work, we present a novel next generation of sequencing (NGS) based approach to identify EV mediated mRNA exchanges between co-cultured adipocyte and macrophage cells. We performed molecular and genomic profiling and jointly considered data from RNA sequencing (RNA-seq) and genotyping to track the "sequence varying mRNAs" transferred between cells. We identified 8 mRNAs being transferred from macrophages to adipocytes and 21 mRNAs being transferred in the opposite direction. These mRNAs represented biological functions including extracellular matrix, cell adhesion, glycoprotein, and signal peptides. CONCLUSIONS: Our study sheds new light on EV mediated RNA communications between adipocyte and macrophage cells, which may play a significant role in developing insulin resistance in diabetic patients. This work establishes a new method that is applicable to examining genetic material exchanges in many cellular systems and has the potential to be extended to in vivo studies as well.

Human Pancreatic ß Cell lncRNAs Control Cell-Specific Regulatory Networks.

Recent studies have uncovered thousands of long non-coding RNAs (lncRNAs) in human pancreatic ß cells. ß cell lncRNAs are often cell type specific and exhibit dynamic regulation during differentiation or upon changing glucose concentrations. Although these features hint at a role of lncRNAs in ß cell gene regulation and diabetes, the function of ß cell lncRNAs remains largely unknown. In this study, we investigated the function of ß cell-specific lncRNAs and transcription factors using transcript knockdowns and co-expression network analysis. This revealed lncRNAs that function in concert with transcription factors to regulate ß cell-specific transcriptional networks. We further demonstrate that the lncRNA PLUTO affects local 3D chromatin structure and transcription of PDX1, encoding a key ß cell transcription factor, and that both PLUTO and PDX1 are downregulated in islets from donors with type 2 diabetes or impaired glucose tolerance. These results implicate lncRNAs in the regulation of ß cell-specific transcription factor networks.

Discover the network underlying the connections between aging and age-related diseases.

Although our knowledge of aging has greatly expanded in the past decades, it remains elusive why and how aging contributes to the development of age-related diseases (ARDs). In particular, a global mechanistic understanding of the connections between aging and ARDs is yet to be established. We rely on a network modelling named "GeroNet" to study the connections between aging and more than a hundred diseases. By evaluating topological connections between aging genes and disease genes in over three thousand subnetworks corresponding to various biological processes, we show that aging has stronger connections with ARD genes compared to non-ARD genes in subnetworks corresponding to "response to decreased oxygen levels", "insulin signalling pathway", "cell cycle", etc. Based on subnetwork connectivity, we can correctly "predict" if a disease is age-related and prioritize the biological processes that are involved in connecting to multiple ARDs. Using Alzheimer's disease (AD) as an example, GeroNet identifies meaningful genes that may play key roles in connecting aging and ARDs. The top modules identified by GeroNet in AD significantly overlap with modules identified from a large scale AD brain gene expression experiment, supporting that GeroNet indeed reveals the underlying biological processes involved in the disease.

Proteogenomics connects somatic mutations to signalling in breast cancer.

Somatic mutations have been extensively characterized in breast cancer, but the effects of these genetic alterations on the proteomic landscape remain poorly understood. Here we describe quantitative mass-spectrometry-based proteomic and phosphoproteomic analyses of 105 genomically annotated breast cancers, of which 77 provided high-quality data. Integrated analyses provided insights into the somatic cancer genome including the consequences of chromosomal loss, such as the 5q deletion characteristic of basal-like breast cancer. Interrogation of the 5q trans-effects against the Library of Integrated Network-based Cellular Signatures, connected loss of CETN3 and SKP1 to elevated expression of epidermal growth factor receptor (EGFR), and SKP1 loss also to increased SRC tyrosine kinase. Global proteomic data confirmed a stromal-enriched group of proteins in addition to basal and luminal clusters, and pathway analysis of the phosphoproteome identified a G-protein-coupled receptor cluster that was not readily identified at the mRNA level. In addition to ERBB2, other amplicon-associated highly phosphorylated kinases were identified, including CDK12, PAK1, PTK2, RIPK2 and TLK2. We demonstrate that proteogenomic analysis of breast cancer elucidates the functional consequences of somatic mutations, narrows candidate nominations for driver genes within large deletions and amplified regions, and identifies therapeutic targets.

Inter-tissue coexpression network analysis reveals DPP4 as an important gene in heart to blood communication.

BACKGROUND: Inter-tissue molecular interactions are critical to the function and behavior of biological systems in multicellular organisms, but systematic studies of interactions between tissues are lacking. Also, existing studies of inter-tissue interactions are based on direct gene expression correlations, which can't distinguish correlations due to common genetic architectures versus biochemical or molecular signal exchange between tissues. METHODS: We developed a novel strategy to study inter-tissue interaction by removing effects of genetic regulation of gene expression (genetic decorrelation). We applied our method to the comprehensive atlas of gene expression across nine human tissues in the Genotype-Tissue Expression (GTEx) project to generate novel genetically decorrelated inter-tissue networks. From this we derived modules of genes important in inter-tissue interactions that are likely driven by biological signal exchange instead of their common genetic basis. Importantly we highlighted communication between tissues and elucidated gene activities in one tissue inducing gene expression changes in others. RESULTS: We reveal global unidirectional inter-tissue coordination of specific biological pathways such as protein synthesis. Using our data, we highlighted a clinically relevant example whereby heart expression of DPP4 was coordinated with a gene expression signature characteristic for whole blood proliferation, potentially impacting peripheral stem cell mobilization. We also showed that expression of the poorly characterized FOCAD in heart correlated with protein biosynthetic processes in the lung. CONCLUSIONS: In summary, this is the first resource of human multi-tissue networks enabling the investigation of molecular inter-tissue interactions. With the networks in hand, we may systematically design combination therapies that simultaneously target multiple tissues or pinpoint potential side effects of a drug in other tissues.