Aging is associated with progressive tissue dysfunction, leading to frailty and mortality. Characterizing aging features, such as changes in gene expression and dynamics, shared across tissues or specific to each tissue, is crucial for understanding systemic and local factors contributing to the aging process. We performed RNA sequencing on 13 tissues at six different ages in male and female African turquoise killifish, the shortest-lived vertebrate that can be raised in captivity. This comprehensive, sex-balanced 'atlas' dataset revealed varying strength of sex-age interactions across killifish tissues and age-altered genes and biological pathways that are evolutionarily conserved in mice and humans. We discovered a female-biased myeloid shift with age in the killifish hematopoietic organ, developed tissue-specific 'transcriptomic clocks' and identified biomarkers predictive of chronological age. We showed the importance of sex-specific clocks for selected tissues, validated the tissue clocks with an independent transcriptomic dataset and used them to evaluate different lifespan interventions in the killifish. Our work provides a comprehensive resource for studying aging dynamics across tissues in the killifish, a powerful vertebrate aging model.

Engulfment by macrophages is critical for waste clearance in the vertebrate brain. Understanding clearance mechanisms may open new therapeutic possibilities to counter brain aging and neurodegenerative diseases. However, few in vivo models exist to study engulfment in the brain and characterize this process during aging and across species. Here we present a genetic model for secretion of a fluorescent protein by neurons in the brain of the African turquoise killifish, the shortest-lived vertebrate that can be bred in captivity. We use this model to identify a population of brain macrophages in the killifish responsible for engulfment of material from the brain extracellular space. Intriguingly, many of these cells bear similarities to mammalian border-associated and monocyte-derived macrophages, rare subsets of macrophages in mouse and human brains noted for their engulfment capabilities. We also find that in our model, killifish brain macrophages decline in engulfment capacity with age. This work highlights how vertebrate brain macrophages, particularly those at brain border regions, can play a critical role in clearance and provides an opportunity to test interventions that can boost engulfment by these macrophages to promote brain resilience in old age and disease.

The state of genome-wide chromatin accessibility in cells, tissues, or organisms can be investigated with a technique called assay for transposase-accessible chromatin using sequencing (ATAC-seq). ATAC-seq is a powerful approach for profiling the epigenomic landscape of cells using very low input materials. Analysis of chromatin accessibility data allows for prediction of gene expression and identification of regulatory elements such as potential enhancers and specific transcription-factor binding sites. Here, we describe an optimized ATAC-seq protocol for the preparation of isolated nuclei and subsequent next-generation sequencing from whole embryos and tissues of the African turquoise killifish (Nothobranchius furzeri). Importantly, we provide an overview of a pipeline for processing and analyzing ATAC-seq data from killifish.

Suspended animation states allow organisms to survive extreme environments. The African turquoise killifish has evolved diapause as a form of suspended development to survive a complete drought. However, the mechanisms underlying the evolution of extreme survival states are unknown. To understand diapause evolution, we performed integrative multi-omics (gene expression, chromatin accessibility, and lipidomics) in the embryos of multiple killifish species. We find that diapause evolved by a recent remodeling of regulatory elements at very ancient gene duplicates (paralogs) present in all vertebrates. CRISPR-Cas9-based perturbations identify the transcription factors REST/NRSF and FOXOs as critical for the diapause gene expression program, including genes involved in lipid metabolism. Indeed, diapause shows a distinct lipid profile, with an increase in triglycerides with very-long-chain fatty acids. Our work suggests a mechanism for the evolution of complex adaptations and offers strategies to promote long-term survival by activating suspended animation programs in other species.

Protein aggregation is a hallmark of age-related neurodegeneration. Yet, aggregation during normal aging and in tissues other than the brain is poorly understood. Here, we leverage the African turquoise killifish to systematically profile protein aggregates in seven tissues of an aging vertebrate. Age-dependent aggregation is strikingly tissue specific and not simply driven by protein expression differences. Experimental interrogation in killifish and yeast, combined with machine learning, indicates that this specificity is linked to protein-autonomous biophysical features and tissue-selective alterations in protein quality control. Co-aggregation of protein quality control machinery during aging may further reduce proteostasis capacity, exacerbating aggregate burden. A segmental progeria model with accelerated aging in specific tissues exhibits selectively increased aggregation in these same tissues. Intriguingly, many age-related protein aggregates arise in wild-type proteins that, when mutated, drive human diseases. Our data chart a comprehensive landscape of protein aggregation during vertebrate aging and identify strong, tissue-specific associations with dysfunction and disease.

Protein aggregation, which can sometimes spread in a prion-like manner, is a hallmark of neurodegenerative diseases. However, whether prion-like aggregates form during normal brain aging remains unknown. Here, we use quantitative proteomics in the African turquoise killifish to identify protein aggregates that accumulate in old vertebrate brains. These aggregates are enriched for prion-like RNA-binding proteins, notably the ATP-dependent RNA helicase DDX5. We validate that DDX5 forms aggregate-like puncta in the brains of old killifish and mice. Interestingly, DDX5's prion-like domain allows these aggregates to propagate across many generations in yeast. In vitro, DDX5 phase separates into condensates. Mutations that abolish DDX5 prion propagation also impair the protein's ability to phase separate. DDX5 condensates exhibit enhanced enzymatic activity, but they can mature into inactive, solid aggregates. Our findings suggest that protein aggregates with prion-like properties form during normal brain aging, which could have implications for the age-dependency of cognitive decline.

Technical advancements over the past two decades have enabled the measurement of the panoply of molecules of cells and tissues including transcriptomes, epigenomes, metabolomes and proteomes at unprecedented resolution. Unbiased profiling of these molecular landscapes in the context of aging can reveal important details about mechanisms underlying age-related functional decline and age-related diseases. However, the high-throughput nature of these experiments creates unique analytical and design demands for robustness and reproducibility. In addition, ‘omic’ experiments are generally onerous, making it crucial to effectively design them to eliminate as many spurious sources of variation as possible as well as account for any biological or technical parameter that may influence such measures. In this Perspective, we provide general guidelines on best practices in the design and analysis of omic experiments in aging research from experimental design to data analysis and considerations for long-term reproducibility and validation of such studies. High-throughput analysis of cellular landscapes is an important tool to decipher the molecular mechanisms driving aging and disease. Here, Singh and Benayoun discuss key considerations in the design and analysis of omic data to gain robust and reproducible insights into the aging process.

The African turquoise killifish is an exciting new vertebrate model for aging studies. A significant challenge for any model organism is the control over its diet in space and time. To address this challenge, we created an automated and networked fish feeding system. Our automated feeder is designed to be open-source, easily transferable, and built from widely available components. Compared to manual feeding, our automated system is highly precise and flexible. As a proof-of-concept for the feeding flexibility of these automated feeders, we define a favorable regimen for growth and fertility for the African killifish and a dietary restriction regimen where both feeding time and quantity are reduced. We show that this dietary restriction regimen extends lifespan in males (but not in females) and impacts the transcriptomes of killifish livers in a sex-specific manner. Moreover, combining our automated feeding system with a video camera, we establish a quantitative associative learning assay to provide an integrative measure of cognitive performance for the killifish. The ability to precisely control food delivery in the killifish opens new areas to assess lifespan and cognitive behavior dynamics and to screen for dietary interventions and drugs in a scalable manner previously impossible with traditional vertebrate model organisms.

Although germline cells are considered to be functionally "immortal," both the germline and supporting somatic cells in the gonad within an organism experience aging. With increased age at parenthood, the age-related decline in reproductive success has become an important biological issue for an aging population. However, molecular mechanisms underlying reproductive aging across sexes in vertebrates remain poorly understood. To decipher molecular drivers of vertebrate gonadal aging across sexes, we perform longitudinal characterization of the gonadal transcriptome throughout the lifespan in the naturally short-lived African turquoise killifish (Nothobranchius furzeri). By combining mRNA-seq and small RNA-seq from 26 individuals, we characterize the aging gonads of young-adult, middle-aged, and old female and male fish. We analyze changes in transcriptional patterns of genes, transposable elements (TEs), and piRNAs. We find that testes seem to undergo only marginal changes during aging. In contrast, in middle-aged ovaries, the time point associated with peak female fertility in this strain, PIWI pathway components are transiently down-regulated, TE transcription is elevated, and piRNA levels generally decrease, suggesting that egg quality may already be declining at middle-age. Furthermore, we show that piRNA ping-pong biogenesis declines steadily with age in ovaries, whereas it is maintained in aging testes. To our knowledge, this data set represents the most comprehensive transcriptomic data set for vertebrate gonadal aging. This resource also highlights important pathways that are regulated during reproductive aging in either ovaries or testes, which could ultimately be leveraged to help restore aspects of youthful reproductive function.

Diapause is a state of suspended development that helps organisms survive extreme environments. How diapause protects living organisms is largely unknown. Using the African turquoise killifish (Nothobranchius furzeri), we show that diapause preserves complex organisms for extremely long periods of time without trade-offs for subsequent adult growth, fertility, and life span. Transcriptome analyses indicate that diapause is an active state, with dynamic regulation of metabolism and organ development genes. The most up-regulated genes in diapause include Polycomb complex members. The chromatin mark regulated by Polycomb, H3K27me3, is maintained at key developmental genes in diapause, and the Polycomb member CBX7 mediates repression of metabolism and muscle genes in diapause. CBX7 is functionally required for muscle preservation and diapause maintenance. Thus, vertebrate diapause is a state of suspended life that is actively maintained by specific chromatin regulators, and this has implications for long-term organism preservation.

Transient receptor potential (TRP) channels form a family of polymodal cation channels gated by thermal, mechanical, or chemical stimuli, with many of them involved in the control of proliferation, apoptosis, or cell cycle. From an evolutionary point of view, TRP family is characterized by high conservation of duplicated genes originating from whole-genome duplication at the onset of vertebrates. The conservation of such "ohnolog" genes is theoretically linked to an increased probability of generating phenotypes deleterious for the organism upon gene mutation. We aimed to test experimentally the hypothesis that TRP mutations, in particular gain-of-function, could be involved in the generation of deleterious phenotypes involved in cancer, such as gain of invasiveness. Indeed, a number of TRP channels have been linked to cancer progression, and exhibit changes in expression levels in various types of cancers. However, TRP mutations in cancer have been poorly documented. We focused on 2 TRPV family members, TRPV4 and TRPV6, and studied the effect of putative gain-of-function mutations on invasiveness properties. TRPV channels have a C-terminal calmodulin-binding domain (CaMBD) that has important functions for regulating protein function, through different mechanisms depending on the channel (channel inactivation/potentiation, cytoskeleton regulation). We studied the effect of mutations mimicking constitutive phosphorylation in TRPV4 and TRPV6 CaMBDs: TRPV4 S823D, S824D and T813D, TRPV6 S691D, S692D and T702. We found that most of these mutants induced a strong gain of invasiveness of colon adenocarcinoma SW480 cells, both for TRPV4 and TRPV6. While increased invasion with TRPV6 S692D and T702D mutants was correlated to increased mutant channel activity, it was not the case for TRPV4 mutants, suggesting different mechanisms with the same global effect of gain in deleterious phenotype. This highlights the potential importance to search for TRP mutations involved in cancer.

All vertebrates including human have evolved from an ancestor that underwent two rounds of whole genome duplication (2R-WGD). In addition, teleost fish underwent an additional third round of genome duplication (3R-WGD). The genes retained from these genome duplications, so-called ohnologs, have been instrumental in the evolution of vertebrate complexity, development and susceptibility to genetic diseases. However, the identification of vertebrate ohnologs has been challenging, due to lineage specific genome rearrangements since 2R- and 3R-WGD. We previously identified vertebrate ohnologs using a novel synteny comparison across multiple genomes. Here, we refine and apply this approach on 27 vertebrate genomes to identify ohnologs from both 2R- and 3R-WGD, while taking into account the phylogenetically biased sampling of available species. We assemble vertebrate ohnolog pairs and families in an expanded OHNOLOGS v2 database. We find that teleost fish have retained more 2R-WGD ohnologs than mammals and sauropsids, and that these 2R-ohnologs have retained significantly more ohnologs from the subsequent 3R-WGD than genes without 2R-ohnologs. Interestingly, species with fewer extant genes, such as sauropsids, have retained similar or higher proportions of ohnologs. OHNOLOGS v2 should allow deeper evolutionary genomic analysis of the impact of WGD on vertebrates and can be freely accessed at http://ohnologs.curie.fr.

Aging is accompanied by the functional decline of tissues. However, a systematic study of epigenomic and transcriptomic changes across tissues during aging is missing. Here, we generated chromatin maps and transcriptomes from four tissues and one cell type from young, middle-aged, and old mice-yielding 143 high-quality data sets. We focused on chromatin marks linked to gene expression regulation and cell identity: histone H3 trimethylation at lysine 4 (H3K4me3), a mark enriched at promoters, and histone H3 acetylation at lysine 27 (H3K27ac), a mark enriched at active enhancers. Epigenomic and transcriptomic landscapes could easily distinguish between ages, and machine-learning analysis showed that specific epigenomic states could predict transcriptional changes during aging. Analysis of data sets from all tissues identified recurrent age-related chromatin and transcriptional changes in key processes, including the up-regulation of immune system response pathways such as the interferon response. The up-regulation of the interferon response pathway with age was accompanied by increased transcription and chromatin remodeling at specific endogenous retroviral sequences. Pathways misregulated during mouse aging across tissues, notably innate immune pathways, were also misregulated with aging in other vertebrate species-African turquoise killifish, rat, and humans-indicating common signatures of age across species. To date, our data set represents the largest multitissue epigenomic and transcriptomic data set for vertebrate aging. This resource identifies chromatin and transcriptional states that are characteristic of young tissues, which could be leveraged to restore aspects of youthful functionality to old tissues.

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.

BACKGROUND:The annual killifish Austrofundulus limnaeus inhabits ephemeral ponds in northern Venezuela, South America, and is an emerging extremophile model for vertebrate diapause, stress tolerance, and evolution. Embryos of A. limnaeus regularly experience extended periods of desiccation and anoxia as a part of their natural history and have unique metabolic and developmental adaptations. Currently, there are limited genomic resources available for gene expression and evolutionary studies that can take advantage of A. limnaeus as a unique model system.RESULTS:We describe the first draft genome sequence of A. limnaeus. The genome was assembled de novo using a merged assembly strategy and was annotated using the NCBI Eukaryotic Annotation Pipeline. We show that the assembled genome has a high degree of completeness in genic regions that is on par with several other teleost genomes. Using RNA-seq and phylogenetic-based approaches, we identify several candidate genes that may be important for embryonic stress tolerance and post-diapause development in A. limnaeus. Several of these genes include heat shock proteins that have unique expression patterns in A. limnaeus embryos and at least one of these may be under positive selection.CONCLUSION:The A. limnaeus genome is the first South American annual killifish genome made publicly available. This genome will be a valuable resource for comparative genomics to determine the genetic and evolutionary mechanisms that support the unique biology of annual killifishes. In a broader context, this genome will be a valuable tool for exploring genome-environment interactions and their impacts on vertebrate physiology and evolution.

Learning causal networks from large-scale genomic data remains challenging in absence of time series or controlled perturbation experiments. We report an information- theoretic method which learns a large class of causal or non-causal graphical models from purely observational data, while including the effects of unobserved latent variables, commonly found in many genomic datasets. Starting from a complete graph, the method iteratively removes dispensable edges, by uncovering significant information contributions from indirect paths, and assesses edge-specific confidences from randomization of available data. The remaining edges are then oriented based on the signature of causality in observational data. The approach and associated algorithm, miic, outperform earlier methods on a broad range of benchmark networks. Causal network reconstructions are presented at different biological size and time scales, from gene regulation in single cells to whole genome duplication in tumor development as well as long term evolution of vertebrates. Miic is publicly available at https://github.com/miicTeam/MIIC.

The discovery that heterozygous and homozygous mutations in the gene encoding progranulin are causally linked to frontotemporal dementia and lysosomal storage disease, respectively, reveals previously unrecognized roles of the progranulin protein in regulating lysosome biogenesis and function. Given the importance of lysosomes in cellular homeostasis, it is not surprising that progranulin deficiency has pleiotropic effects on neural circuit development and maintenance, stress response, innate immunity and ageing. This Progress article reviews recent advances in progranulin biology emphasizing its roles in lysosomal function and brain innate immunity, and outlines future avenues of investigation that may lead to new therapeutic approaches for neurodegeneration.

Whole genome duplications (WGD) have now been firmly established in all major eukaryotic kingdoms. In particular, all vertebrates descend from two rounds of WGDs, that occurred in their jawless ancestor some 500 MY ago. Paralogs retained from WGD, also coined 'ohnologs' after Susumu Ohno, have been shown to be typically associated with development, signaling and gene regulation. Ohnologs, which amount to about 20 to 35% of genes in the human genome, have also been shown to be prone to dominant deleterious mutations and frequently implicated in cancer and genetic diseases. Hence, identifying ohnologs is central to better understand the evolution of vertebrates and their susceptibility to genetic diseases. Early computational analyses to identify vertebrate ohnologs relied on content-based synteny comparisons between the human genome and a single invertebrate outgroup genome or within the human genome itself. These approaches are thus limited by lineage specific rearrangements in individual genomes. We report, in this study, the identification of vertebrate ohnologs based on the quantitative assessment and integration of synteny conservation between six amniote vertebrates and six invertebrate outgroups. Such a synteny comparison across multiple genomes is shown to enhance the statistical power of ohnolog identification in vertebrates compared to earlier approaches, by overcoming lineage specific genome rearrangements. Ohnolog gene families can be browsed and downloaded for three statistical confidence levels or recompiled for specific, user-defined, significance criteria at http://ohnologs.curie.fr/. In the light of the importance of WGD on the genetic makeup of vertebrates, our analysis provides a useful resource for researchers interested in gaining further insights on vertebrate evolution and genetic diseases.