The progressive decline of physiological function and the increased risk of age-related diseases challenge healthy aging. Multiple anti-aging manipulations, such as senolytics, have proven beneficial for health; however, the biomarkers that label in vivo senescence at systemic levels are lacking, thus hindering anti-aging applications. In this study, we generate a Glb1+/m‒Glb1-2A-mCherry (GAC) reporter allele at the Glb1 gene locus, which encodes lysosomal β-galactosidase-an enzyme elevated in tissues of old mice. A linear correlation between GAC signal and chronological age is established in a cohort of middle-aged (9 to 13 months) Glb1+/m mice. The high GAC signal is closely associated with cardiac hypertrophy and a shortened lifespan. Moreover, the GAC signal is exponentially increased in pathological senescence induced by bleomycin in the lung. Senolytic dasatinib and quercetin (D + Q) reduce GAC signal in bleomycin treated mice. Thus, the Glb1-2A-mCherry reporter mice monitors systemic aging and function decline, predicts lifespan, and may facilitate the understanding of aging mechanisms and help in the development of anti-aging interventions.

TGF-β signaling is necessary for CD8+ T cell differentiation into tissue resident memory T cells (TRM). Although higher frequency of CD8+ TRM cells in the tumor microenvironment is associated with better prognosis, TGF-β-blockade typically improves rather than worsens outcomes. Here we show that in a mouse melanoma model, in the tumor-draining lymph nodes (TDLN) rather than in the tumors themselves, stem-like CD8+ T cells differentiate into TRMs in a TGF-β and tumor antigen dependent manner. Following vaccination against a melanoma-specific epitope, most tumour-specific CD8+ T cells are maintained in a stem-like state, but a proportion of cells lost TRM status and differentiate into CX3CR1+ effector CD8+ T cells in the TDLN, which are subsequently migrating into the tumours. Disruption of TGF-β signaling changes the dynamics of these developmental processes, with the net result of improving effector CD8+ T cell migration into the tumours. In summary, TDLN stem-like T cells transiently switch from a TGF-β-dependent TRM differentiation program to an anti-tumor migratory effector development upon vaccination, which transition can be facilitated by targeted TGF-β blockade.

Fasting glucose and insulin are intermediate traits for type 2 diabetes. Here we explore the role of coding variation on these traits by analysis of variants on the HumanExome BeadChip in 60,564 non-diabetic individuals and in 16,491 T2D cases and 81,877 controls. We identify a novel association of a low-frequency nonsynonymous SNV in GLP1R (A316T; rs10305492; MAF=1.4%) with lower FG (β=-0.09±0.01 mmol l(-1), P=3.4 × 10(-12)), T2D risk (OR[95%CI]=0.86[0.76-0.96], P=0.010), early insulin secretion (β=-0.07±0.035 pmolinsulin mmolglucose(-1), P=0.048), but higher 2-h glucose (β=0.16±0.05 mmol l(-1), P=4.3 × 10(-4)). We identify a gene-based association with FG at G6PC2 (pSKAT=6.8 × 10(-6)) driven by four rare protein-coding SNVs (H177Y, Y207S, R283X and S324P). We identify rs651007 (MAF=20%) in the first intron of ABO at the putative promoter of an antisense lncRNA, associating with higher FG (β=0.02±0.004 mmol l(-1), P=1.3 × 10(-8)). Our approach identifies novel coding variant associations and extends the allelic spectrum of variation underlying diabetes-related quantitative traits and T2D susceptibility.

Reduced glomerular filtration rate defines chronic kidney disease and is associated with cardiovascular and all-cause mortality. We conducted a meta-analysis of genome-wide association studies for estimated glomerular filtration rate (eGFR), combining data across 133,413 individuals with replication in up to 42,166 individuals. We identify 24 new and confirm 29 previously identified loci. Of these 53 loci, 19 associate with eGFR among individuals with diabetes. Using bioinformatics, we show that identified genes at eGFR loci are enriched for expression in kidney tissues and in pathways relevant for kidney development and transmembrane transporter activity, kidney structure, and regulation of glucose metabolism. Chromatin state mapping and DNase I hypersensitivity analyses across adult tissues demonstrate preferential mapping of associated variants to regulatory regions in kidney but not extra-renal tissues. These findings suggest that genetic determinants of eGFR are mediated largely through direct effects within the kidney and highlight important cell types and biological pathways.

17β-hydroxysteroid dehydrogenase-13 is a hepatocyte-specific, lipid droplet-associated protein. A common loss-of-function variant of HSD17B13 (rs72613567: TA) protects patients against non-alcoholic fatty liver disease with underlying mechanism incompletely understood. In the present study, we identify the serine 33 of 17β-HSD13 as an evolutionally conserved PKA target site and its phosphorylation facilitates lipolysis by promoting its interaction with ATGL on lipid droplets. Targeted mutation of Ser33 to Ala (S33A) decreases ATGL-dependent lipolysis in cultured hepatocytes by reducing CGI-58-mediated ATGL activation. Importantly, a transgenic knock-in mouse strain carrying the HSD17B13 S33A mutation (HSD17B1333A/A) spontaneously develops hepatic steatosis with reduced lipolysis and increased inflammation. Moreover, Hsd17B1333A/A mice are more susceptible to high-fat diet-induced nonalcoholic steatohepatitis. Finally, we find reproterol, a potential 17β-HSD13 modulator and FDA-approved drug, confers a protection against nonalcoholic steatohepatitis via PKA-mediated Ser33 phosphorylation of 17β-HSD13. Therefore, targeting the Ser33 phosphorylation site could represent a potential approach to treat NASH.

X-chromosomal genetic variants are understudied but can yield valuable insights into sexually dimorphic human traits and diseases. We performed a sex-stratified cross-ancestry X-chromosome-wide association meta-analysis of seven kidney-related traits (n = 908,697), identifying 23 loci genome-wide significantly associated with two of the traits: 7 for uric acid and 16 for estimated glomerular filtration rate (eGFR), including four novel eGFR loci containing the functionally plausible prioritized genes ACSL4, CLDN2, TSPAN6 and the female-specific DRP2. Further, we identified five novel sex-interactions, comprising male-specific effects at FAM9B and AR/EDA2R, and three sex-differential findings with larger genetic effect sizes in males at DCAF12L1 and MST4 and larger effect sizes in females at HPRT1. All prioritized genes in loci showing significant sex-interactions were located next to androgen response elements (ARE). Five ARE genes showed sex-differential expressions. This study contributes new insights into sex-dimorphisms of kidney traits along with new prioritized gene targets for further molecular research.

Chronic kidney disease (CKD) is defined by reduced estimated glomerular filtration rate (eGFR). Previous genetic studies have implicated regulatory mechanisms contributing to CKD. Here we present epigenome-wide association studies of eGFR and CKD using whole-blood DNA methylation of 2264 ARIC Study and 2595 Framingham Heart Study participants to identify epigenetic signatures of kidney function. Of 19 CpG sites significantly associated (P < 1e-07) with eGFR/CKD and replicated, five also associate with renal fibrosis in biopsies from CKD patients and show concordant DNA methylation changes in kidney cortex. Lead CpGs at PTPN6/PHB2, ANKRD11, and TNRC18 map to active enhancers in kidney cortex. At PTPN6/PHB2 cg19942083, methylation in kidney cortex associates with lower renal PTPN6 expression, higher eGFR, and less renal fibrosis. The regions containing the 243 eGFR-associated (P < 1e-05) CpGs are significantly enriched for transcription factor binding sites of EBF1, EP300, and CEBPB (P < 5e-6). Our findings highlight kidney function associated epigenetic variation.

Distant metastasis is the main cause of breast cancer-related death; however, effective therapeutic strategies targeting metastasis are still scarce. This is largely attributable to the spatiotemporal intratumor heterogeneity during metastasis. Here we show that protein deacetylase SIRT7 is significantly downregulated in breast cancer lung metastases in human and mice, and predicts metastasis-free survival. SIRT7 deficiency promotes breast cancer cell metastasis, while temporal expression of Sirt7 inhibits metastasis in polyomavirus middle T antigen breast cancer model. Mechanistically, SIRT7 deacetylates and promotes SMAD4 degradation mediated by β-TrCP1, and SIRT7 deficiency activates transforming growth factor-β signaling and enhances epithelial-to-mesenchymal transition. Significantly, resveratrol activates SIRT7 deacetylase activity, inhibits breast cancer lung metastases, and increases survival. Our data highlight SIRT7 as a modulator of transforming growth factor-β signaling and suppressor of breast cancer metastasis, meanwhile providing an effective anti-metastatic therapeutic strategy.Metastatic disease is the major reason for breast cancer-related deaths; therefore, a better understanding of this process and its players is needed. Here the authors report the role of SIRT7 in inhibiting SMAD4-mediated breast cancer metastasis providing a possible therapeutic avenue.

Large consortia have revealed hundreds of genetic loci associated with anthropometric traits, one trait at a time. We examined whether genetic variants affect body shape as a composite phenotype that is represented by a combination of anthropometric traits. We developed an approach that calculates averaged PCs (AvPCs) representing body shape derived from six anthropometric traits (body mass index, height, weight, waist and hip circumference, waist-to-hip ratio). The first four AvPCs explain >99% of the variability, are heritable, and associate with cardiometabolic outcomes. We performed genome-wide association analyses for each body shape composite phenotype across 65 studies and meta-analysed summary statistics. We identify six novel loci: LEMD2 and CD47 for AvPC1, RPS6KA5/C14orf159 and GANAB for AvPC3, and ARL15 and ANP32 for AvPC4. Our findings highlight the value of using multiple traits to define complex phenotypes for discovery, which are not captured by single-trait analyses, and may shed light onto new pathways.

Few genome-wide association studies (GWAS) account for environmental exposures, like smoking, potentially impacting the overall trait variance when investigating the genetic contribution to obesity-related traits. Here, we use GWAS data from 51,080 current smokers and 190,178 nonsmokers (87% European descent) to identify loci influencing BMI and central adiposity, measured as waist circumference and waist-to-hip ratio both adjusted for BMI. We identify 23 novel genetic loci, and 9 loci with convincing evidence of gene-smoking interaction (GxSMK) on obesity-related traits. We show consistent direction of effect for all identified loci and significance for 18 novel and for 5 interaction loci in an independent study sample. These loci highlight novel biological functions, including response to oxidative stress, addictive behaviour, and regulatory functions emphasizing the importance of accounting for environment in genetic analyses. Our results suggest that tobacco smoking may alter the genetic susceptibility to overall adiposity and body fat distribution.

Increased levels of the urinary albumin-to-creatinine ratio (UACR) are associated with higher risk of kidney disease progression and cardiovascular events, but underlying mechanisms are incompletely understood. Here, we conduct trans-ethnic (n = 564,257) and European-ancestry specific meta-analyses of genome-wide association studies of UACR, including ancestry- and diabetes-specific analyses, and identify 68 UACR-associated loci. Genetic correlation analyses and risk score associations in an independent electronic medical records database (n = 192,868) reveal connections with proteinuria, hyperlipidemia, gout, and hypertension. Fine-mapping and trans-Omics analyses with gene expression in 47 tissues and plasma protein levels implicate genes potentially operating through differential expression in kidney (including TGFB1, MUC1, PRKCI, and OAF), and allow coupling of UACR associations to altered plasma OAF concentrations. Knockdown of OAF and PRKCI orthologs in Drosophila nephrocytes reduces albumin endocytosis. Silencing fly PRKCI further impairs slit diaphragm formation. These results generate a priority list of genes and pathways for translational research to reduce albuminuria.

Dietary interventions such as intermittent fasting (IF) have emerged as an attractive strategy for cancer therapies; therefore, understanding the underlying molecular mechanisms is pivotal. Here, we find SIRT7 decline markedly attenuates the anti-tumor effect of IF. Mechanistically, AMP-activated protein kinase (AMPK) phosphorylating SIRT7 at T263 triggers further phosphorylation at T255/S259 by glycogen synthase kinase 3β (GSK3β), which stabilizes SIRT7 by decoupling E3 ligase UBR5. SIRT7 hyperphosphorylation achieves anti-tumor activity by disrupting the SKP2-SCF E3 ligase, thus preventing SKP2-mediated K63-linked AKT polyubiquitination and subsequent activation. In contrast, GSK3β-SIRT7 axis is inhibited by EGF/ERK2 signaling, with ERK2 inactivating GSK3β, thus accelerating SIRT7 degradation. Unfavorably, glucose deprivation or chemotherapy hijacks the GSK3β-SIRT7 axis via ERK2, thus activating AKT and ensuring survival. Notably, Trametinib, an FDA-approved MEK inhibitor, enhances the efficacy of combination therapy with doxorubicin and IF. Overall, we have revealed the GSK3β-SIRT7 axis that must be fine-tuned in the face of the energetic and oncogenic stresses in malignancy.

Differences between sexes contribute to variation in the levels of fasting glucose and insulin. Epidemiological studies established a higher prevalence of impaired fasting glucose in men and impaired glucose tolerance in women, however, the genetic component underlying this phenomenon is not established. We assess sex-dimorphic (73,089/50,404 women and 67,506/47,806 men) and sex-combined (151,188/105,056 individuals) fasting glucose/fasting insulin genetic effects via genome-wide association study meta-analyses in individuals of European descent without diabetes. Here we report sex dimorphism in allelic effects on fasting insulin at IRS1 and ZNF12 loci, the latter showing higher RNA expression in whole blood in women compared to men. We also observe sex-homogeneous effects on fasting glucose at seven novel loci. Fasting insulin in women shows stronger genetic correlations than in men with waist-to-hip ratio and anorexia nervosa. Furthermore, waist-to-hip ratio is causally related to insulin resistance in women, but not in men. These results position dissection of metabolic and glycemic health sex dimorphism as a steppingstone for understanding differences in genetic effects between women and men in related phenotypes.

Statins effectively lower LDL cholesterol levels in large studies and the observed interindividual response variability may be partially explained by genetic variation. Here we perform a pharmacogenetic meta-analysis of genome-wide association studies (GWAS) in studies addressing the LDL cholesterol response to statins, including up to 18,596 statin-treated subjects. We validate the most promising signals in a further 22,318 statin recipients and identify two loci, SORT1/CELSR2/PSRC1 and SLCO1B1, not previously identified in GWAS. Moreover, we confirm the previously described associations with APOE and LPA. Our findings advance the understanding of the pharmacogenetic architecture of statin response.

Despite existing reports on differential DNA methylation in type 2 diabetes (T2D) and obesity, our understanding of its functional relevance remains limited. Here we show the effect of differential methylation in the early phases of T2D pathology by a blood-based epigenome-wide association study of 4808 non-diabetic Europeans in the discovery phase and 11,750 individuals in the replication. We identify CpGs in LETM1, RBM20, IRS2, MAN2A2 and the 1q25.3 region associated with fasting insulin, and in FCRL6, SLAMF1, APOBEC3H and the 15q26.1 region with fasting glucose. In silico cross-omics analyses highlight the role of differential methylation in the crosstalk between the adaptive immune system and glucose homeostasis. The differential methylation explains at least 16.9% of the association between obesity and insulin. Our study sheds light on the biological interactions between genetic variants driving differential methylation and gene expression in the early pathogenesis of T2D.

Elevated serum urate levels can cause gout, an excruciating disease with suboptimal treatment. Previous GWAS identified common variants with modest effects on serum urate. Here we report large-scale whole-exome sequencing association studies of serum urate and kidney function among ≤19,517 European ancestry and African-American individuals. We identify aggregate associations of low-frequency damaging variants in the urate transporters SLC22A12 (URAT1; p = 1.3 × 10-56) and SLC2A9 (p = 4.5 × 10-7). Gout risk in rare SLC22A12 variant carriers is halved (OR = 0.5, p = 4.9 × 10-3). Selected rare variants in SLC22A12 are validated in transport studies, confirming three as loss-of-function (R325W, R405C, and T467M) and illustrating the therapeutic potential of the new URAT1-blocker lesinurad. In SLC2A9, mapping of rare variants of large effects onto the predicted protein structure reveals new residues that may affect urate binding. These findings provide new insights into the genetic architecture of serum urate, and highlight molecular targets in SLC22A12 and SLC2A9 for lowering serum urate and preventing gout.