Methylation entropy landscape of Chinese long-lived individuals reveals lower epigenetic noise related to human healthy aging.

The transition from ordered to noisy is a significant epigenetic signature of aging and age-related disease. As a paradigm of healthy human aging and longevity, long-lived individuals (LLI, >90 years old) may possess characteristic strategies in coping with the disordered epigenetic regulation. In this study, we constructed high-resolution blood epigenetic noise landscapes for this cohort by a methylation entropy (ME) method using whole genome bisulfite sequencing (WGBS). Although a universal increase in global ME occurred with chronological age in general control samples, this trend was suppressed in LLIs. Importantly, we identified 38,923 genomic regions with LLI-specific lower ME (LLI-specific lower entropy regions, for short, LLI-specific LERs). These regions were overrepresented in promoters, which likely function in transcriptional noise suppression. Genes associated with LLI-specific LERs have a considerable impact on SNP-based heritability of some aging-related disorders (e.g., asthma and stroke). Furthermore, neutrophil was identified as the primary cell type sustaining LLI-specific LERs. Our results highlight the stability of epigenetic order in promoters of genes involved with aging and age-related disorders within LLI epigenomes. This unique epigenetic feature reveals a previously unknown role of epigenetic order maintenance in specific genomic regions of LLIs, which helps open a new avenue on the epigenetic regulation mechanism in human healthy aging and longevity.

Scrutiny of genome-wide somatic mutation profiles in centenarians identifies the key genomic regions for human longevity.

Somatic mutations accumulate with age and are associated closely with human health, their characterization in longevity cohorts remains largely unknown. Here, by analyzing whole genome somatic mutation profiles in 73 centenarians and 51 younger controls in China, we found that centenarian genomes are characterized by a markedly skewed distribution of somatic mutations, with many genomic regions being specifically conserved but displaying a high function potential. This, together with the observed more efficient DNA repair ability in the long-lived individuals, supports the existence of key genomic regions for human survival during aging, with their integrity being of essential to human longevity.

Longevity-Associated Transcription Factor ATF7 Promotes Healthspan by Suppressing Cellular Senescence and Systematic Inflammation.

Aging is characterized by persistent low-grade systematic inflammation, which is largely responsible for the occurrence of various age-associated diseases. We and others have previously reported that long-lived people (such as centenarians) can delay the onset of or even escape certain major age-related diseases. Here, by screening blood transcriptome and inflammatory profiles, we found that long-lived individuals had a relatively lower inflammation level (IL6, TNFα), accompanied by up-regulation of activating transcription factor 7 (ATF7). Interestingly, ATF7 expression was gradually reduced during cellular senescence. Loss of ATF7 induced cellular senescence, while overexpression delayed senescence progress and senescence-associated secretory phenotype (SASP) secretion. We showed that the anti-senescence effects of ATF7 were achieved by inhibiting nuclear factor kappa B (NF-κB) signaling and increasing histone H3K9 dimethylation (H3K9me2). In Caenorhabditis elegans, ATF7 overexpression significantly suppressed aging biomarkers and extended lifespan. Our findings suggest that ATF7 is a longevity-promoting factor that lowers cellular senescence and inflammation in long-lived individuals.

Mitogenome evidence shows two radiation events and dispersals of matrilineal ancestry from northern coastal China to the Americas and Japan.

Although it is widely recognized that the ancestors of Native Americans (NAs) primarily came from Siberia, the link between mitochondrial DNA (mtDNA) lineage D4h3a (typical of NAs) and D4h3b (found so far only in East China and Thailand) raises the possibility that the ancestral sources for early NAs were more variegated than hypothesized. Here, we analyze 216 contemporary (including 106 newly sequenced) D4h mitogenomes and 39 previously reported ancient D4h data. The results reveal two radiation events of D4h in northern coastal China, one during the Last Glacial Maximum and the other within the last deglaciation, which facilitated the dispersals of D4h sub-branches to different areas including the Americas and the Japanese archipelago. The coastal distributions of the NA (D4h3a) and Japanese lineages (D4h1a and D4h2), in combination with the Paleolithic archaeological similarities among Northern China, the Americas, and Japan, lend support to the coastal dispersal scenario of early NAs.

Metabolic Modeling Identifies a Novel Molecular Type of Glioblastoma Associated with Good Prognosis.

Glioblastoma (GBM) is one of the most aggressive forms of cancer. Although IDH1 mutation indicates a good prognosis and a potential target for treatment, most GBMs are IDH1 wild-type. Identifying additional molecular markers would help to generate personalized therapies and improve patient outcomes. Here, we used our recently developed metabolic modeling method (genome-wide precision metabolic modeling, GPMM) to investigate the metabolic profiles of GBM, aiming to identify additional novel molecular markers for this disease. We systematically analyzed the metabolic reaction profiles of 149 GBM samples lacking IDH1 mutation. Forty-eight reactions showing significant association with prognosis were identified. Further analysis indicated that the purine recycling, nucleotide interconversion, and folate metabolism pathways were the most robust modules related to prognosis. Considering the three pathways, we then identified the most significant GBM type for a better prognosis, namely N+P-. This type presented high nucleotide interconversion (N+) and low purine recycling (P-). N+P--type exhibited a significantly better outcome (log-rank p = 4.7 × 10-7) than that of N-P+. GBM patients with the N+P--type had a median survival time of 19.6 months and lived 65% longer than other GBM patients. Our results highlighted a novel molecular type of GBM, which showed relatively high frequency (26%) in GBM patients lacking the IDH1 mutation, and therefore exhibits potential in GBM prognostic assessment and personalized therapy.

The landscape of aging.

Aging is characterized by a progressive deterioration of physiological integrity, leading to impaired functional ability and ultimately increased susceptibility to death. It is a major risk factor for chronic human diseases, including cardiovascular disease, diabetes, neurological degeneration, and cancer. Therefore, the growing emphasis on "healthy aging" raises a series of important questions in life and social sciences. In recent years, there has been unprecedented progress in aging research, particularly the discovery that the rate of aging is at least partly controlled by evolutionarily conserved genetic pathways and biological processes. In an attempt to bring full-fledged understanding to both the aging process and age-associated diseases, we review the descriptive, conceptual, and interventive aspects of the landscape of aging composed of a number of layers at the cellular, tissue, organ, organ system, and organismal levels.

Multiple time-series expression trajectories imply dynamic functional changes during cellular senescence.

Cellular senescence is a dynamic process driven by epigenetic and genetic changes. Although some transcriptomic signatures of senescent cells have been discovered, how these senescence-related signals change over time remains largely unclear. Here, we profiled the transcriptome dynamics of human dermal fibroblast (HDF) cells in successive stages of growth from proliferation to senescence. Based on time-series expression profile analysis, we discovered four trajectories (C1, C2, C3, C4) that are dynamically expressed as senescence progresses. While some genes were continuously up-regulated (C4) or down-regulated (C2) with aging, other genes did not change linearly with cell proliferation, but remained stable until entering the senescent state (C1, C3). Further analysis revealed that the four modes were enriched in different biological pathways, including regulation of cellular senescence. These findings provide a new perspective on understanding the dynamic regulatory mechanism of cellular senescence.

Mechanism of Thymus Rejuvenation in Elderly Macaques.

Senile thymus atrophy is an important factor leading to decreased immune function. Repairing the atrophic thymus tissue structure, rebuilding immune function, and replenishing the number of exogenous stem cells may be ideal methods. In this study, bone marrow mesenchymal stem cells were intravenously infused into elderly macaques. We found that thymus volume was substantially increased, some thymus tissue regeneration was observed, the degree of thymus tissue fibrosis decreased, collagen fiber deposition decreased, cortical and medulla structures emerged gradually, the number of apoptotic cells decreased significantly, and the expression of apoptosis-related proteins decreased. For the effects of stem cell therapy on aging-related genes, we performed transcriptomic analysis of thymus tissue. The results show the expression pattern of the tissue transcriptome tended to be similar to the thymus expression pattern in young macaques compared with the elderly group, reverse aging-related proteins. Based on the results, it is suggested that stem cell therapy is an ideal method to prevent or reverse the aging of the thymus.

Specific Gain and Loss of Co-Expression Modules in Long-Lived Individuals Indicate a Role of circRNAs in Human Longevity.

Deep RNA sequencing of 164 blood samples collected from long-lived families was performed to investigate the expression patterns of circular RNAs (circRNAs). Unlike that observed in previous studies, circRNA expression in long-lived elderly individuals (98.3 ± 3.4 year) did not exhibit an age-accumulating pattern. Based on weighted circRNA co-expression network analysis, we found that longevous elders specifically gained eight but lost seven conserved circRNA-circRNA co-expression modules (c-CCMs) compared with normal elder controls (spouses of offspring of long-lived individuals, age = 59.3 ± 5.8 year). Further analysis showed that these modules were associated with healthy aging-related pathways. These results together suggest an important role of circRNAs in regulating human lifespan extension.

ETS1 acts as a regulator of human healthy aging via decreasing ribosomal activity.

Adaptation to reduced energy production during aging is a fundamental issue for maintaining healthspan or prolonging life span. Currently, however, the underlying mechanism in long-lived people remains poorly understood. Here, we analyzed transcriptomes of 185 long-lived individuals (LLIs) and 86 spouses of their children from two independent Chinese longevity cohorts and found that the ribosome pathway was significantly down-regulated in LLIs. We found that the down-regulation is likely controlled by ETS1 (ETS proto-oncogene 1), a transcription factor down-regulated in LLIs and positively coexpressed with most ribosomal protein genes (RPGs). Functional assays showed that ETS1 can bind to RPG promoters, while ETS1 knockdown reduces RPG expression and alleviates cellular senescence in human dermal fibroblast (HDF) and embryonic lung fibroblast (IMR-90) cells. As protein synthesis/turnover in ribosomes is an energy-intensive cellular process, the decline in ribosomal biogenesis governed by ETS1 in certain female LLIs may serve as an alternative mechanism to achieve energy-saving and healthy aging.

Why Senescent Cells Are Resistant to Apoptosis: An Insight for Senolytic Development.

Cellular senescence is a process that leads to a state of irreversible cell growth arrest induced by a variety of intrinsic and extrinsic stresses. Senescent cells (SnCs) accumulate with age and have been implicated in various age-related diseases in part via expressing the senescence-associated secretory phenotype. Elimination of SnCs has the potential to delay aging, treat age-related diseases and extend healthspan. However, once cells becoming senescent, they are more resistant to apoptotic stimuli. Senolytics can selectively eliminate SnCs by targeting the SnC anti-apoptotic pathways (SCAPs). They have been developed as a novel pharmacological strategy to treat various age-related diseases. However, the heterogeneity of the SnCs indicates that SnCs depend on different proteins or pathways for their survival. Thus, a better understanding of the underlying mechanisms for apoptotic resistance of SnCs will provide new molecular targets for the development of cell-specific or broad-spectrum therapeutics to clear SnCs. In this review, we discussed the latest research progresses and challenge in senolytic development, described the significance of regulation of senescence and apoptosis in aging, and systematically summarized the SCAPs involved in the apoptotic resistance in SnCs.

System-level metabolic modeling facilitates unveiling metabolic signature in exceptional longevity.

Although it is well known that metabolic control plays a crucial role in regulating the health span and life span of various organisms, little is known for the systems metabolic profile of centenarians, the paradigm of human healthy aging and longevity. Meanwhile, how to well characterize the system-level metabolic states in an organism of interest remains to be a major challenge in systems metabolism research. To address this challenge and better understand the metabolic mechanisms of healthy aging, we developed a method of genome-wide precision metabolic modeling (GPMM) which is able to quantitatively integrate transcriptome, proteome and kinetome data in predictive modeling of metabolic networks. Benchmarking analysis showed that GPMM successfully characterized metabolic reprogramming in the NCI-60 cancer cell lines; it dramatically improved the performance of the modeling with an R2 of 0.86 between the predicted and experimental measurements over the performance of existing methods. Using this approach, we examined the metabolic networks of a Chinese centenarian cohort and identified the elevated fatty acid oxidation (FAO) as the most significant metabolic feature in these long-lived individuals. Evidence from serum metabolomics supports this observation. Given that FAO declines with normal aging and is impaired in many age-related diseases, our study suggests that the elevated FAO has potential to be a novel signature of healthy aging of humans.

A pair of long intergenic non-coding RNA LINC00887 variants act antagonistically to control Carbonic Anhydrase IX transcription upon hypoxia in tongue squamous carcinoma progression.

BACKGROUND: Long noncoding RNAs (lncRNAs) are important regulators in tumor progression. However, their biological functions and underlying mechanisms in hypoxia adaptation remain largely unclear. RESULTS: Here, we established a correlation between a Chr3q29-derived lncRNA gene and tongue squamous carcinoma (TSCC) by genome-wide analyses. Using RACE, we determined that two novel variants of this lncRNA gene are generated in TSCC, namely LINC00887_TSCC_short (887S) and LINC00887_TSCC_long (887L). RNA-sequencing in 887S or 887L loss-of-function cells identified their common downstream target as Carbonic Anhydrase IX (CA9), a gene known to be upregulated by hypoxia during tumor progression. Mechanistically, our results showed that the hypoxia-augmented 887S and constitutively expressed 887L functioned in opposite directions on tumor progression through the common target CA9. Upon normoxia, 887S and 887L interacted. Upon hypoxia, the two variants were separated. Each RNA recognized and bound to their responsive DNA cis-acting elements on CA9 promoter: 887L activated CA9's transcription through recruiting HIF1α, while 887S suppressed CA9 through DNMT1-mediated DNA methylation. CONCLUSIONS: We provided hypoxia-permitted functions of two antagonistic lncRNA variants to fine control the hypoxia adaptation through CA9.

Complete mitogenomes document substantial genetic contribution from the Eurasian Steppe into northern Pakistani Indo-Iranian speakers.

To elucidate whether Bronze Age population dispersals from the Eurasian Steppe to South Asia contributed to the gene pool of Indo-Iranian-speaking groups, we analyzed 19,568 mitochondrial DNA (mtDNA) sequences from northern Pakistani and surrounding populations, including 213 newly generated mitochondrial genomes (mitogenomes) from Iranian and Dardic groups, both speakers from the ancient Indo-Iranian branch in northern Pakistan. Our results showed that 23% of mtDNA lineages with west Eurasian origin arose in situ in northern Pakistan since ~5000 years ago (kya), a time depth very close to the documented Indo-European dispersals into South Asia during the Bronze Age. Together with ancient mitogenomes from western Eurasia since the Neolithic, we identified five haplogroups (~8.4% of maternal gene pool) with roots in the Steppe region and subbranches arising (age ~5-2 kya old) in northern Pakistan as genetic legacies of Indo-Iranian speakers. Some of these haplogroups, such as W3a1b that have been found in the ancient samples from the late Bronze Age to the Iron Age period individuals of Swat Valley northern Pakistan, even have sub-lineages (age ~4 kya old) in the southern subcontinent, consistent with the southward spread of Indo-Iranian languages. By showing that substantial genetic components of Indo-Iranian speakers in northern Pakistan can be traced to Bronze Age in the Steppe region, our study suggests a demographic link with the spread of Indo-Iranian languages, and further highlights the corridor role of northern Pakistan in the southward dispersal of Indo-Iranian-speaking groups.

Exploring European ancestry among the Kalash population: a mitogenomic perspective.

With a population of around 4 000 individuals, the Kalash people have been living in the Hindu-Kush mountain valleys of present-day northern Pakistan for centuries. Due to their mysterious origin and fairer European complexion, the genetic history of this ethnic group has been investigated previously using different markers. To date, however, the maternal genetic architecture has not been systematically dissected based on high-resolution complete mitochondrial genomes (mitogenomes), making their maternal genetic history, especially their genetic connection with Europeans from a matrilineal perspective, unclear. To unravel this issue, we analyzed mitogenome data of 34 Kalash samples together with 6 075 individuals from across Eurasia. Our results indicated exclusive western Eurasian origin of the Kalash people, represented by eight haplogroups. Among these haplogroups, J2b1a7a and R0a5a (accounting for ~50% of the Kalash gene pool) displayed in situ differentiations in the Kalash and could be traced to the Mediterranean region. Age estimations suggested these haplogroups arose in the Kalash population ~2.26 and 3.01 thousand years ago (kya), a time frame consistent with the invasion of Alexander III of Macedon to the region. One possible explanation for the maternal genetic contribution from Europeans to the Kalash people would be the involvement of women in foreign campaigns of ancient Greek warfare, followed by a founder effect. Our study thus sheds important light on the genetic origin of the Kalash community of Pakistan.

Comparative analysis of long noncoding RNAs in long-lived mammals provides insights into natural cancer-resistance.

Mouse and rats are staple model organisms that have been traditionally used for oncological studies; however, their short lifespan and highly prone to cancers limit their utilizationsin understanding the mechanisms of cancer resistance. In recent years, several studies of the non-standard long-lived mammalian species like naked mole rat (NMR) have provided new insights of mechanisms in natural anti-cancer. How long-lived species genetically maintain longevity and cancer-resistance remains largely elusive. To better understand the underlying anti-cancer mechanisms in long-lived mammals, we genome widely identified long noncoding RNA (lncRNA) transcripts of two longevous mammals, bowhead whale (BW, Balaena mysticetus) and Brandt's bat (BB, Myotis brandtii) and featured their sequence traits, expression patterns, and their correlations with cancer-resistance. Similar with naked mole rat (NMR, Heterocephalus glaber), the most long-lived rodent, BW and BB lncRNAs show low sequence conservation and dynamic expressions among tissues and physiological stages. By utilizing k-mers clustering, 75-136 of BW, BB and NMR lncRNAs were found in close relation (Pearson's r ≥0.9, p < 0.01) with human ageing diseases related lncRNAs (HAR-Lncs). In addition, we observed thousands of BB and BW lncRNAs strongly co-expressed (r > 0.8 or r <-0.8, p < 0.01) with potential tumour suppressors, indicating that lncRNAs are potentially involved in anti-cancer regulation in long-lived mammals. Our study provides the basis for lncRNA researches in perspectives of evolution and anti-cancer studies. Abbreviations: BW: bowhead whale; BB: Brandt's bat; NMR: naked mole rat; LLM: long-lived mammal; HTS: human tumour-suppressors; PTS: potential tumour suppressor; ARD: ageing related diseases; HAR-Lncs: lncRNAs that related with human ageing diseases; Kmer-lncs: lncRNAs in long-lived mammal species that corelated (Pearson'sr ≥0.9, p < 0.01) with the 10 HAR-Lncs by k-mers clustering; All-lncs: all the lncRNAs in long-lived mammal species; SDE-lncs: significant differentially expressed lncRNAs.

Identification of DNA N6-methyladenine sites by integration of sequence features.

BACKGROUND: An increasing number of nucleic acid modifications have been profiled with the development of sequencing technologies. DNA N6-methyladenine (6mA), which is a prevalent epigenetic modification, plays important roles in a series of biological processes. So far, identification of DNA 6mA relies primarily on time-consuming and expensive experimental approaches. However, in silico methods can be implemented to conduct preliminary screening to save experimental resources and time, especially given the rapid accumulation of sequencing data. RESULTS: In this study, we constructed a 6mA predictor, p6mA, from a series of sequence-based features, including physicochemical properties, position-specific triple-nucleotide propensity (PSTNP), and electron-ion interaction pseudopotential (EIIP). We performed maximum relevance maximum distance (MRMD) analysis to select key features and used the Extreme Gradient Boosting (XGBoost) algorithm to build our predictor. Results demonstrated that p6mA outperformed other existing predictors using different datasets. CONCLUSIONS: p6mA can predict the methylation status of DNA adenines, using only sequence files. It may be used as a tool to help the study of 6mA distribution pattern. Users can download it from https://github.com/Konglab404/p6mA.