Quantitative trait loci (QTL) analysis of longevity in C57BL/6J by DBA/2J (BXD) recombinant inbred mice

Genetic loci and metabolic states associated with murine epigenetic aging

Changes in DNA methylation (DNAm) are linked to aging. Here, we profile highly conserved CpGs in 339 predominantly female mice belonging to the BXD family for which we have deep longevity and genomic data. We use a 'pan-mammalian' microarray that provides a common platform for assaying the methylome across mammalian clades. We computed epigenetic clocks and tested associations with DNAm entropy, diet, weight, metabolic traits, and genetic variation. We describe the multifactorial variance of methylation at these CpGs and show that high-fat diet augments the age-related changes. Entropy increases with age. The progression to disorder, particularly at CpGs that gain methylation over time, was predictive of genotype-dependent life expectancy. The longer-lived BXD strains had comparatively lower entropy at a given age. We identified two genetic loci that modulate epigenetic age acceleration (EAA): one on chromosome (Chr) 11 that encompasses the Erbb2/Her2 oncogenic region, and the other on Chr19 that contains a cytochrome P450 cluster. Both loci harbor genes associated with EAA in humans, including STXBP4, NKX2-3, and CUTC. Transcriptome and proteome analyses revealed correlations with oxidation-reduction, metabolic, and immune response pathways. Our results highlight concordant loci for EAA in humans and mice, and demonstrate a tight coupling between the metabolic state and epigenetic aging.

Multiomic profiling of the liver across diets and age in a diverse mouse population

We profiled the liver transcriptome, proteome, and metabolome in 347 individuals from 58 isogenic strains of the BXD mouse population across age (7 to 24 months) and diet (low or high fat) to link molecular variations to metabolic traits. Several hundred genes are affected by diet and/or age at the transcript and protein levels. Orthologs of two aging-associated genes, St7 and Ctsd, were knocked down in C. elegans, reducing longevity in wild-type and mutant long-lived strains. The multiomics data were analyzed as segregating gene networks according to each independent variable, providing causal insight into dietary and aging effects. Candidates were cross-examined in an independent diversity outbred mouse liver dataset segregating for similar diets, with ∼80%-90% of diet-related candidate genes found in common across datasets. Together, we have developed a large multiomics resource for multivariate analysis of complex traits and demonstrate a methodology for moving from observational associations to causal connections.

Gene-by-environment modulation of lifespan and weight gain in the murine BXD family

How lifespan and body weight vary as a function of diet and genetic differences is not well understood. Here we quantify the impact of differences in diet on lifespan in a genetically diverse family of female mice, split into matched isogenic cohorts fed a low-fat chow diet (CD, n = 663) or a high-fat diet (HFD, n = 685). We further generate key metabolic data in a parallel cohort euthanized at four time points. HFD feeding shortens lifespan by 12%: equivalent to a decade in humans. Initial body weight and early weight gains account for longevity differences of roughly 4-6 days per gram. At 500 days, animals on a HFD typically gain four times as much weight as control, but variation in weight gain does not correlate with lifespan. Classic serum metabolites, often regarded as health biomarkers, are not necessarily strong predictors of longevity. Our data indicate that responses to a HFD are substantially modulated by gene-by-environment interactions, highlighting the importance of genetic variation in making accurate individualized dietary recommendations.

Lifespan-Associated Gene Expression Signatures of Recombinant BXD Mice Implicates Coro7 and Set in Longevity

Although genetic approaches have identified key genes and pathways that promote longevity, systems-level approaches are less utilized. Here, we took advantage of the wealth of omics data characterizing the BXD family of mice. We associated transcript and peptide levels across five tissues from both female and male BXD isogenic lines with their median lifespan. We identified over 5000 genes that showed a longevity correlation in a given tissue. Surprisingly, we found less than 1% overlap among longevity-correlating genes across tissues and sex. These 1% shared genes consist of 51 genes, of which 13 have been shown to alter lifespan. Only two genes -Coro7 and Set- showed a longevity correlation in all tissues and in both sexes. While differential regulation of aging across tissues and sex has been reported, our systems-level analysis reveals two unique genes that may promote healthy aging in unique sex- and tissue-agnostic manner.

Body weight and high-fat diet are associated with epigenetic aging in female members of the BXD murine family

DNA methylation (DNAm) is shaped by genetic and environmental factors and modulated by aging. Here, we examine interrelations between epigenetic aging, body weight (BW), and life span in 12 isogenic strains from the BXD family of mice that exhibit over twofold variation in longevity. Genome-wide DNAm was assayed in 70 liver specimens from predominantly female cases, 6-25 months old, that were maintained on normal chow or high-fat diet (HFD). We defined subsets of CpG regions associated with age, BW at young adulthood, and strain-by-diet-dependent life span. These age-associated differentially methylated CpG regions (age-DMRs) featured distinct genomic characteristics, with DNAm gains over time occurring in sites such as promoters and exons that have high CpG density and low average methylation. CpG regions associated with BW were enriched in introns, tended to have lower methylation in mice with higher BW, and were inversely correlated with gene expression (i.e., higher mRNA levels in mice with higher BW). CpG regions associated with life span were linked to genes involved in life span modulation, including the telomerase reverse transcriptase gene, Tert, which had both lower methylation and higher expression in long-lived strains. An epigenetic clock defined from age-DMRs revealed accelerated aging in mice belonging to strains with shorter life spans. Both higher BW and the HFD were associated with accelerated epigenetic aging. Our results highlight the age-accelerating effect of heavier BW. Furthermore, we demonstrate that the measure of epigenetic aging derived from age-DMRs can predict genotype and diet-induced differences in life span among female BXD members.

Genetic Basis of Increased Lifespan and Postponed Senescence in Drosophila melanogaster

Limited lifespan and senescence are near-universal phenomena. These quantitative traits exhibit variation in natural populations due to the segregation of many interacting loci and from environmental effects. Due to the complexity of the genetic control of lifespan and senescence, our understanding of the genetic basis of variation in these traits is incomplete. Here, we analyzed the pattern of genetic divergence between long-lived (O) Drosophila melanogaster lines selected for postponed reproductive senescence and unselected control (B) lines. We quantified the productivity of the O and B lines and found that reproductive senescence is maternally controlled. We therefore chose 57 candidate genes that are expressed in ovaries, 49 of which have human orthologs, and assessed the effects of RNA interference in ovaries and accessary glands on lifespan and reproduction. All but one candidate gene affected at least one life history trait in one sex or productivity week. In addition, 23 genes had antagonistic pleiotropic effects on lifespan and productivity. Identifying evolutionarily conserved genes affecting increased lifespan and delayed reproductive senescence is the first step toward understanding the evolutionary forces that maintain segregating variation at these loci in nature and may provide potential targets for therapeutic intervention to delay senescence while increasing lifespan.

The Genetics of Aging: A Vertebrate Perspective

Aging negatively impacts vitality and health. Many genetic pathways that regulate aging were discovered in invertebrates. However, the genetics of aging is more complex in vertebrates because of their specialized systems. This Review discusses advances in the genetic regulation of aging in vertebrates from work in mice, humans, and organisms with exceptional lifespans. We highlight challenges for the future, including sex-dependent differences in lifespan and the interplay between genes and environment. We also discuss how the identification of reliable biomarkers of age and development of new vertebrate models can be leveraged for personalized interventions to counter aging and age-related diseases.

Genetic cartography of longevity in humans and mice: Current landscape and horizons

Aging is a complex and highly variable process. Heritability of longevity among humans and other species is low, and this finding has given rise to the idea that it may be futile to search for DNA variants that modulate aging. We argue that the problem in mapping longevity genes is mainly one of low power and the genetic and environmental complexity of aging. In this review we highlight progress made in mapping genes and molecular networks associated with longevity, paying special attention to work in mice and humans. We summarize 40 years of linkage studies using murine cohorts and 15 years of studies in human populations that have exploited candidate gene and genome-wide association methods. A small but growing number of gene variants contribute to known longevity mechanisms, but a much larger set have unknown functions. We outline these and other challenges and suggest some possible solutions, including more intense collaboration between research communities that use model organisms and human cohorts. Once hundreds of gene variants have been linked to differences in longevity in mammals, it will become feasible to systematically explore gene-by-environmental interactions, dissect mechanisms with more assurance, and evaluate the roles of epistasis and epigenetics in aging. A deeper understanding of complex networks-genetic, cellular, physiological, and social-should position us well to improve healthspan.

Profiling DNA methylation differences between inbred mouse strains on the Illumina Human Infinium MethylationEPIC microarray

The Illumina Infinium MethylationEPIC provides an efficient platform for profiling DNA methylation in humans at over 850,000 CpGs. Model organisms such as mice do not currently benefit from an equivalent array. Here we used this array to measure DNA methylation in mice. We defined probes targeting conserved regions and performed differential methylation analysis and compared between the array-based assay and affinity-based DNA sequencing of methyl-CpGs (MBD-seq) and reduced representation bisulfite sequencing. Mouse samples consisted of 11 liver DNA from two strains, C57BL/6J (B6) and DBA/2J (D2), that varied widely in age. Linear regression was applied to detect differential methylation. In total, 13,665 probes (1.6% of total probes) aligned to conserved CpGs. Beta-values (β-value) for these probes showed a distribution similar to that in humans. Overall, there was high concordance in methylation signal between the EPIC array and MBD-seq (Pearson correlation r = 0.70, p-value < 0.0001). However, the EPIC probes had higher quantitative sensitivity at CpGs that are hypo- (β-value < 0.3) or hypermethylated (β-value > 0.7). In terms of differential methylation, no EPIC probe detected a significant difference between age groups at a Benjamini-Hochberg threshold of 10%, and the MBD-seq performed better at detecting age-dependent change in methylation. However, the top most significant probe for age (cg13269407; uncorrected p-value = 1.8 x 10-5) is part of the clock CpGs used to estimate the human epigenetic age. For strain, 219 EPIC probes detected significant differential methylation (FDR cutoff 10%) with ~80% CpGs associated with higher methylation in D2. This higher methylation profile in D2 compared to B6 was also replicated by the MBD-seq data. To summarize, we found only a small subset of EPIC probes that target conserved sites. However, for this small subset the array provides a reliable assay of DNA methylation and can be effectively used to measure differential methylation in mice.

Complex Genetics of Behavior: BXDs in the Automated Home-Cage

This chapter describes a use case for the genetic dissection and automated analysis of complex behavioral traits using the genetically diverse panel of BXD mouse recombinant inbred strains. Strains of the BXD resource differ widely in terms of gene and protein expression in the brain, as well as in their behavioral repertoire. A large mouse resource opens the possibility for gene finding studies underlying distinct behavioral phenotypes, however, such a resource poses a challenge in behavioral phenotyping. To address the specifics of large-scale screening we describe how to investigate: (1) how to assess mouse behavior systematically in addressing a large genetic cohort, (2) how to dissect automation-derived longitudinal mouse behavior into quantitative parameters, and (3) how to map these quantitative traits to the genome, deriving loci underlying aspects of behavior.

A Chromosome 13 locus is associated with male-specific mortality in mice

BACKGROUND AND AIM

Mortality is a highly complex trait influenced by a wide array of genetic factors.

METHODS

We examined a population of 1200 mice that were F2 generation offspring of a 4-way reciprocal cross between C57BL6/J and DBA2/J strains. Animals were sacrificed at age 200, 500, or 800 days and genotyped at 96 markers. The 800 days old cohort, which were the survivors of a much larger breeding group, were examined for enriched frequency of alleles that benefit survival and depletion of alleles that reduce survival.

RESULTS

Loci on Chr 13 in males and on Chr X in females were significantly distorted from Mendelian expectations, even after conservative correction for multiple testing. DBA2/J alleles between 35 and 80 Mb on Chr 13 were underrepresented in the age 800 male animals. D2 genotypes in this region were also associated with premature death during behavioral testing. Furthermore, confirmatory analysis showed BXD recombinant inbred strains carrying the D2 alleles in this region had shorter median survival. Exploration of available pathology data indicated that a syndrome involving dental malocclusions, pancreatic islet hypertrophy, and kidney lipidosis may have mediated the effects of DBA alleles on mortality specifically in male mice. The heterozygote advantage locus on the X Chr was not found to be associated with any pathology.

CONCLUSIONS

These results suggest a novel locus influencing survival in the B6/D2 genetic background, perhaps via a metabolic disorder that emerges by 200 days of age in male animals.

The Genomic Basis of Postponed Senescence in Drosophila melanogaster

Natural populations harbor considerable genetic variation for lifespan. While evolutionary theory provides general explanations for the existence of this variation, our knowledge of the genes harboring naturally occurring polymorphisms affecting lifespan is limited. Here, we assessed the genetic divergence between five Drosophila melanogaster lines selected for postponed senescence for over 170 generations (O lines) and five lines from the same base population maintained at a two week generation interval for over 850 generations (B lines). On average, O lines live 70% longer than B lines, are more productive at all ages, and have delayed senescence for other traits than reproduction. We performed population sequencing of pools of individuals from all B and O lines and identified 6,394 genetically divergent variants in or near 1,928 genes at a false discovery rate of 0.068. A 2.6 Mb region at the tip of the X chromosome contained many variants fixed for alternative alleles in the two populations, suggestive of a hard selective sweep. We also assessed genome wide gene expression of O and B lines at one and five weeks of age using RNA sequencing and identified genes with significant (false discovery rate < 0.05) effects on gene expression with age, population and the age by population interaction, separately for each sex. We identified transcripts that exhibited the transcriptional signature of postponed senescence and integrated the gene expression and genetic divergence data to identify 98 (175) top candidate genes in females (males) affecting postponed senescence and increased lifespan. While several of these genes have been previously associated with Drosophila lifespan, most are novel and constitute a rich resource for future functional validation.

Genetic Regulation of Female Sexual Maturation and Longevity Through Circulating IGF1

We previously reported that insulin-like growth factor 1 (IGF1) was involved in coregulating female sexual maturation and longevity. To understand the underlying genetic mechanisms, based on the strain survey assays of development and aging traits, we crossed two mouse strains, KK/HIJ and PL/J, and produced 307 female F2 mice. We observed the age of vaginal patency (AVP) and the life span of these females. We also measured circulating IGF1 level at 7, 16, 24, 52, and 76 weeks. IGF1 level at 7 weeks significantly correlated with AVP. IGF1 levels at ages of 52 and 76 weeks negatively correlated with longevity (p ≤ .05). A gene mapping study found 22, 4 ,and 3 quantitative trait loci for IGF1, AVP, and life span, respectively. Importantly, the colocalization of IGF1, AVP, and life span quantitative trait loci in the distal region of chromosome 2 suggests this locus carries gene(s) that could regulate IGF1, AVP, and life span. In this region, proprotein convertase subtilisin/kexin type 2 has been found to be associated with female sexual maturation in a human genome-wide association study. We verified the roles of proprotein convertase subtilisin/kexin type 2 in regulating IGF1 and AVP by showing that depletion of proprotein convertase subtilisin/kexin type 2 significantly reduced IGF1 and delayed AVP in mice, suggesting that it also might be involved in the regulation of aging.

Microarray-assisted fine mapping of quantitative trait loci on chromosome 15 for susceptibility to seizure-induced cell death in mice

Prior studies with crosses of the FVB/NJ (FVB; seizure-induced cell death-susceptible) mouse and the C57BL/6J (B6; seizure-induced cell death-resistant) mouse revealed the presence of a quantitative trait locus (QTL) on chromosome 15 that influenced susceptibility to kainic acid-induced cell death (Sicd2). In an earlier study, we confirmed that the Sicd2 interval harbors gene(s) conferring strong protection against seizure-induced cell death through the creation of the FVB.B6-Sicd2 congenic strain, and created three interval-specific congenic lines (ISCLs) that encompass Sicd2 on chromosome 15 to fine-map this locus. To further localise this Sicd2 QTL, an additional congenic line carrying overlapping intervals of the B6 segment was created (ISCL-4), and compared with the previously created ISCL-1-ISCL-3 and assessed for seizure-induced cell death phenotype. Whereas all of the ISCLs showed reduced cell death associated with the B6 phenotype, ISCL-4, showed the most extensive reduction in seizure-induced cell death throughout all hippocampal subfields. In order to characterise the susceptibility loci on Sicd2 by use of this ISCL and identify compelling candidate genes, we undertook an integrative genomic strategy of comparing exon transcript abundance in the hippocampus of this newly developed chromosome 15 subcongenic line (ISCL-4) and FVB-like littermates. We identified 10 putative candidate genes that are alternatively spliced between the strains and may govern strain-dependent differences in susceptibility to seizure-induced excitotoxic cell death. These results illustrate the importance of identifying transcriptomics variants in expression studies, and implicate novel candidate genes conferring susceptibility to seizure-induced cell death.

Genetic variation in responses to dietary restriction--an unbiased tool for hypothesis testing

Dietary restriction (DR) extends lifespan in a wide range of animal models. A major obstacle to understanding how DR modulates lifespan and aging-related dysfunction is the multiplicity of physiological and molecular changes associated with DR. Unraveling their importance to the longevity effect of DR remains a major challenge. In this perspective, we review the marked genetic variation in the response to DR of multiple recombinant inbred (RI) mouse strains. We illustrate how this genetic variation can be exploited to probe the mechanisms mediating lifespan extension by DR, as well as uncover its limits as an intervention. RI strains exhibit marked variation in their lifespan as well as physiological responses to DR. Quantitative genetic and statistical tools can use this phenotypic variation to probe the importance of physiological and molecular changes that have been hypothesized to play roles in DR-mediated lifespan extension.

Concise review: hematopoietic stem cell aging, life span, and transplantation

Self-renewal and multilineage differentiation of stem cells are keys to the lifelong homeostatic maintenance of tissues and organs. Hematopoietic aging, characterized by immunosenescence, proinflammation, and anemia, is attributed to age-associated changes in the number and function of hematopoietic stem cells (HSCs) and their microenvironmental niche. Genetic variants and factors regulating stem cell aging are correlatively or causatively associated with overall organismal aging and longevity. Translational use of HSCs for transplantation and gene therapy demands effective methods for stem cell expansion. Targeting the molecular pathways involved in HSC self-renewal, proliferation, and homing has led to enhanced expansion and engraftment of stem cells upon transplantation. HSC transplantation is less effective in elderly people, even though this is the demographic with the greatest need for this form of treatment. Thus, understanding the biological changes in the aging of stem cells as well as local and systematic environments will improve the efficacy of aged stem cells for regenerative medicine and ultimately facilitate improved health and life spans.

Pleiotropy of segregating genetic variants that affect honey bee worker life expectancy

In contrast to many other complex traits, the natural genetic architecture of life expectancy has not been intensely studied, particularly in non-model organisms, such as the honey bee (Apis mellifera L.). Multiple factors that determine honey bee worker lifespan have been identified and genetic analyses have been performed on some of those traits. Several of the traits are included in a suite of correlated traits that form the pollen hoarding syndrome, which was named after the behavior to store surplus pollen in the nest and is tied to social evolution. Here, seven quantitative trait loci that had previously been identified for their effects on different aspects of the pollen hoarding syndrome were studied for their genetic influence on the survival of adult honey bee workers. To gain a more comprehensive understanding of the genetic architecture of worker longevity, a panel of 280 additional SNP markers distributed across the genome was also tested. Allelic distributions were compared between young and old bees in two backcross populations of the bi-directionally selected high- and low-pollen hoarding strain. Our results suggest a pleiotropic effect of at least one of the behavioral quantitative trait loci on worker longevity and one significant and several other putative genetic effects in other genomic regions. At least one locus showed evidence for strong antagonistic pleiotropy and several others suggested genetic factors that influence pre-emergence survival of worker honey bees. Thus, the predicted association between worker lifespan and the pollen hoarding syndrome was supported at the genetic level and the magnitude of the identified effects also strengthened the view that naturally segregating genetic variation can have major effects on age-specific survival probability in the wild.

Mapping of quantitative trait loci controlling lifespan in the short-lived fish Nothobranchius furzeri--a new vertebrate model for age research

The African annual fish Nothobranchius furzeri emerged as a new model for age research over recent years. Nothobranchius furzeri show an exceptionally short lifespan, age-dependent cognitive/behavioral decline, expression of age-related biomarkers, and susceptibility to lifespan manipulation. In addition, laboratory strains differ largely in lifespan. Here, we set out to study the genetics of lifespan determination. We crossed a short- to a long-lived strain, recorded lifespan, and established polymorphic markers. On the basis of genotypes of 411 marker loci in 404 F(2) progeny, we built a genetic map comprising 355 markers at an average spacing of 5.5 cM, 22 linkage groups (LGs) and 1965 cM. By combining marker data with lifespan values, we identified one genome-wide highly significant quantitative trait locus (QTL) on LG 9 (P < 0.01), which explained 11.3% of the F(2) lifespan variance, and three suggestive QTLs on LG 11, 14, and 17. We characterized the highly significant QTL by synteny analysis, because a genome sequence of N. furzeri was not available. We located the syntenic region on medaka chromosome 5, identified candidate genes, and performed fine mapping, resulting in a c. 40% reduction of the initial 95% confidence interval. We show both that lifespan determination in N. furzeri is polygenic, and that candidate gene detection is easily feasible by cross-species analysis. Our work provides first results on the way to identify loci controlling lifespan in N. furzeri and illustrates the potential of this vertebrate species as a genetic model for age research.

Mice as a mammalian model for research on the genetics of aging

Mice are an ideal mammalian model for studying the genetics of aging: considerable resources are available, the generation time is short, and the environment can be easily controlled, an important consideration when performing mapping studies to identify genes that influence lifespan and age-related diseases. In this review we highlight some salient contributions of the mouse in aging research: lifespan intervention studies in the Interventions Testing Program of the National Institute on Aging; identification of the genetic underpinnings of the effects of calorie restriction on lifespan; the Aging Phenome Project at the Jackson Laboratory, which has submitted multiple large, freely available phenotyping datasets to the Mouse Phenome Database; insights from spontaneous and engineered mouse mutants; and complex traits analyses identifying quantitative trait loci that affect lifespan. We also show that genomewide association peaks for lifespan in humans and lifespan quantitative loci for mice map to homologous locations in the genome. Thus, the vast bioinformatic and genetic resources of the mouse can be used to screen candidate genes identified in both mouse and human mapping studies, followed by functional testing, often not possible in humans, to determine their influence on aging.