Genome-Wide Gene Expression in relation to Age in Large Laboratory Cohorts of Drosophila melanogaster

Programmed ageing: decline of stem cell renewal, immunosenescence, and Alzheimer's disease

The characteristic maximum lifespan varies enormously across animal species from a few hours to hundreds of years. This argues that maximum lifespan, and the ageing process that itself dictates lifespan, are to a large extent genetically determined. Although controversial, this is supported by firm evidence that semelparous species display evolutionarily programmed ageing in response to reproductive and environmental cues. Parabiosis experiments reveal that ageing is orchestrated systemically through the circulation, accompanied by programmed changes in hormone levels across a lifetime. This implies that, like the circadian and circannual clocks, there is a master 'clock of age' (circavital clock) located in the limbic brain of mammals that modulates systemic changes in growth factor and hormone secretion over the lifespan, as well as systemic alterations in gene expression as revealed by genomic methylation analysis. Studies on accelerated ageing in mice, as well as human longevity genes, converge on evolutionarily conserved fibroblast growth factors (FGFs) and their receptors, including KLOTHO, as well as insulin-like growth factors (IGFs) and steroid hormones, as key players mediating the systemic effects of ageing. Age-related changes in these and multiple other factors are inferred to cause a progressive decline in tissue maintenance through failure of stem cell replenishment. This most severely affects the immune system, which requires constant renewal from bone marrow stem cells. Age-related immune decline increases risk of infection whereas lifespan can be extended in germfree animals. This and other evidence suggests that infection is the major cause of death in higher organisms. Immune decline is also associated with age-related diseases. Taking the example of Alzheimer's disease (AD), we assess the evidence that AD is caused by immunosenescence and infection. The signature protein of AD brain, Aβ, is now known to be an antimicrobial peptide, and Aβ deposits in AD brain may be a response to infection rather than a cause of disease. Because some cognitively normal elderly individuals show extensive neuropathology, we argue that the location of the pathology is crucial - specifically, lesions to limbic brain are likely to accentuate immunosenescence, and could thus underlie a vicious cycle of accelerated immune decline and microbial proliferation that culminates in AD. This general model may extend to other age-related diseases, and we propose a general paradigm of organismal senescence in which declining stem cell proliferation leads to programmed immunosenescence and mortality.

Drosophila melanogaster as a model to study age and sex differences in brain injury and neurodegeneration after mild head trauma

Repetitive physical insults to the head, including those that elicit mild traumatic brain injury (mTBI), are a known risk factor for a variety of neurodegenerative conditions including Alzheimer's disease (AD), Parkinson's disease (PD), and chronic traumatic encephalopathy (CTE). Although most individuals who sustain mTBI typically achieve a seemingly full recovery within a few weeks, a subset experience delayed-onset symptoms later in life. As most mTBI research has focused on the acute phase of injury, there is an incomplete understanding of mechanisms related to the late-life emergence of neurodegeneration after early exposure to mild head trauma. The recent adoption of Drosophila-based brain injury models provides several unique advantages over existing preclinical animal models, including a tractable framework amenable to high-throughput assays and short relative lifespan conducive to lifelong mechanistic investigation. The use of flies also provides an opportunity to investigate important risk factors associated with neurodegenerative conditions, specifically age and sex. In this review, we survey current literature that examines age and sex as contributing factors to head trauma-mediated neurodegeneration in humans and preclinical models, including mammalian and Drosophila models. We discuss similarities and disparities between human and fly in aging, sex differences, and pathophysiology. Finally, we highlight Drosophila as an effective tool for investigating mechanisms underlying head trauma-induced neurodegeneration and for identifying therapeutic targets for treatment and recovery.

The genomic basis of copper tolerance in Drosophila is shaped by a complex interplay of regulatory and environmental factors

BACKGROUND

Escalation in industrialization and anthropogenic activity have resulted in an increase of pollutants released into the environment. Of these pollutants, heavy metals such as copper are particularly concerning due to their bio-accumulative nature. Due to its highly heterogeneous distribution and its dual nature as an essential micronutrient and toxic element, the genetic basis of copper tolerance is likely shaped by a complex interplay of genetic and environmental factors.

RESULTS

In this study, we utilized the natural variation present in multiple populations of Drosophila melanogaster collected across Europe to screen for variation in copper tolerance. We found that latitude and the degree of urbanization at the collection sites, rather than any other combination of environmental factors, were linked to copper tolerance. While previously identified copper-related genes were not differentially expressed in tolerant vs. sensitive strains, genes involved in metabolism, reproduction, and protease induction contributed to the differential stress response. Additionally, the greatest transcriptomic and physiological responses to copper toxicity were seen in the midgut, where we found that preservation of gut acidity is strongly linked to greater tolerance. Finally, we identified transposable element insertions likely to play a role in copper stress response.

CONCLUSIONS

Overall, by combining genome-wide approaches with environmental association analysis, and functional analysis of candidate genes, our study provides a unique perspective on the genetic and environmental factors that shape copper tolerance in natural D. melanogaster populations and identifies new genes, transposable elements, and physiological traits involved in this complex phenotype.

EPGAT: Gene Essentiality Prediction With Graph Attention Networks

Identifying essential genes and proteins is a critical step towards a better understanding of human biology and pathology. Computational approaches helped to mitigate experimental constraints by exploring machine learning (ML) methods and the correlation of essentiality with biological information, especially protein-protein interaction (PPI) networks, to predict essential genes. Nonetheless, their performance is still limited, as network-based centralities are not exclusive proxies of essentiality, and traditional ML methods are unable to learn from non-euclidean domains such as graphs. Given these limitations, we proposed EPGAT, an approach for Essentiality Prediction based on Graph Attention Networks (GATs), which are attention-based Graph Neural Networks (GNNs), operating on graph-structured data. Our model directly learns gene essentiality patterns from PPI networks, integrating additional evidence from multiomics data encoded as node attributes. We benchmarked EPGAT for four organisms, including humans, accurately predicting gene essentiality with ROC AUC score ranging from 0.78 to 0.97. Our model significantly outperformed network-based and shallow ML-based methods and achieved a very competitive performance against the state-of-the-art node2vec embedding method. Notably, EPGAT was the most robust approach in scenarios with limited and imbalanced training data. Thus, the proposed approach offers a powerful and effective way to identify essential genes and proteins.

Age-Related Changes of Gene Expression Profiles in Drosophila

An individual’s gene expression profile changes throughout their life. This change in gene expression is shaped by differences in physiological needs and functions between the younger and older organism. Despite intensive studies, the aging process is not fully understood, and several genes involved in this process may remain to be identified. Here we report a transcriptomic analysis of Drosophila melanogaster using microarrays. We compared the expression profiles of two-day-old female adult flies with those of 45-day-old flies. We identified 1184 genes with pronounced differences in expression level between young and old age groups. Most genes involved in muscle development/maintenance that display different levels of expression with age were downregulated in older flies. Many of these genes contributed to sarcomere formation and function. Several of these genes were functionally related to direct and indirect flight muscles; some of them were exclusively expressed in these muscles. Conversely, several genes involved in apoptosis processes were upregulated in aging flies. In addition, several genes involved in resistance to toxic chemicals were upregulated in aging flies, which is consistent with a global upregulation of the defense response system in aging flies. Finally, we randomly selected 12 genes among 232 genes with unknown function and generated transgenic flies expressing recombinant proteins fused with GFP protein to determine their subcellular expression. We also found that the knockdown of some of those 12 genes can affect the lifespan of flies.

The Y chromosome may contribute to sex-specific ageing in Drosophila

Heterochromatin suppresses repetitive DNA, and a loss of heterochromatin has been observed in aged cells of several species, including humans and Drosophila. Males often contain substantially more heterochromatic DNA than females, due to the presence of a large, repeat-rich Y chromosome, and male flies generally have shorter average life spans than females. Here we show that repetitive DNA becomes de-repressed more rapidly in old male flies relative to females, and repeats on the Y chromosome are disproportionally mis-expressed during aging. This is associated with a loss of heterochromatin at repetitive elements during aging in male flies, and a general loss of repressive chromatin in aged males away from pericentromeric regions and the Y. By generating flies with different sex chromosome karyotypes (XXY females; X0 and XYY males), we show that repeat de-repression and average lifespan is correlated with the number of Y chromosomes. This suggests that sex-specific chromatin differences may contribute to sex-specific aging in flies.

Transcriptome Analysis of Long-lived Drosophila melanogaster E(z) Mutants Sheds Light on the Molecular Mechanisms of Longevity

The E(z) histone methyltransferase heterozygous mutation in Drosophila is known to increase lifespan and stress resistance. However, the longevity mechanisms of E(z) mutants have not been revealed. Using genome-wide transcriptome analysis, we demonstrated that lifespan extension, increase of resistance to hyperthermia, oxidative stress and endoplasmic reticulum stress, and fecundity enhancement in E(z) heterozygous mutants are accompanied by changes in the expression level of 239 genes (p < 0.05). Our results demonstrated sex-specific effects of E(z) mutation on gene expression, which, however, did not lead to differences in lifespan extension in both sexes. We observed that a mutation in an E(z) gene leads to perturbations in gene expression, most of which participates in metabolism, such as Carbohydrate metabolism, Lipid metabolism, Drug metabolism, Nucleotide metabolism. Age-dependent changes in the expression of genes involved in pathways related to immune response, cell cycle, and ribosome biogenesis were found.

Aging in the Drosophila ovary: contrasting changes in the expression of the piRNA machinery and mitochondria but no global release of transposable elements

Background

Evolutionary theory indicates that the dynamics of aging in the soma and reproductive tissues may be distinct. This difference arises from the fact that only the germline lineage establishes future generations. In the soma, changes in the landscape of heterochromatin have been proposed to have an important role in aging. This is because redistribution of heterochromatin during aging has been linked to the derepression of transposable elements and an overall loss of somatic gene regulation. A role for changes in the chromatin landscape in the aging of reproductive tissues is less well established. Whether or not epigenetic factors, such as heterochromatin marks, are perturbed in aging reproductive tissues is of interest because, in special cases, epigenetic variation may be heritable. Using mRNA sequencing data from late-stage egg chambers in Drosophila melanogaster, we characterized the landscape of altered gene and transposable element expression in aged reproductive tissues. This allowed us to test the hypothesis that reproductive tissues may differ from somatic tissues in their response to aging.

Results

We show that age-related expression changes in late-stage egg chambers tend to occur in genes residing in heterochromatin, particularly on the largely heterochromatic 4th chromosome. However, these expression differences are seen as both decreases and increases during aging, inconsistent with a general loss of heterochromatic silencing. We also identify an increase in expression of the piRNA machinery, suggesting an age-related increased investment in the maintenance of genome stability. We further identify a strong age-related reduction in the expression of mitochondrial transcripts. However, we find no evidence for global TE derepression in reproductive tissues. Rather, the observed effects of aging on TEs are primarily strain and family specific.

Conclusions

These results identify unique responses in somatic versus reproductive tissue with regards to aging. As in somatic tissues, female reproductive tissues show reduced expression of mitochondrial genes. In contrast, the piRNA machinery shows increased expression during aging. Overall, these results also indicate that global loss of TE control observed in other studies may be unique to the soma and sensitive to genetic background and TE family.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5668-3) contains supplementary material, which is available to authorized users.

The Neuronal Overexpression of Gclc in Drosophila melanogaster Induces Life Extension With Longevity-Associated Transcriptomic Changes in the Thorax

Some effects of aging in animals are tissue-specific. In D. melanogaster neuronal overexpression of Gclc increases lifespan and improves certain physiological parameters associated with health benefits such as locomotor activity, circadian rhythmicity, and stress resistance. Our previous transcriptomic analyses of Drosophila heads, primarily composed of neuronal tissue, revealed significant changes in expression levels of genes involved in aging-related signaling pathways (Jak-STAT, MAPK, FOXO, Notch, mTOR, TGF-beta), translation, protein processing in endoplasmic reticulum, proteasomal degradation, glycolysis, oxidative phosphorylation, apoptosis, regulation of circadian rhythms, differentiation of neurons, synaptic plasticity, and transmission. Considering that various tissues age differently and age-related gene expression changes are tissue-specific, we investigated the effects of neuronal Gclc overexpression on gene expression levels in the imago thorax, which is primarily composed of muscles. A total of 58 genes were found to be differentially expressed between thoraces of control and Gclc overexpressing flies. The Gclc level demonstrated associations with expression of genes involved in the circadian rhythmicity, the genes in categories related to the muscle system process and the downregulation of genes involved in proteolysis. Most of the functional categories altered by Gclc overexpression related to metabolism including Drug metabolism, Metabolism of xenobiotics by cytochrome P450, Glutathione metabolism, Starch and sucrose metabolism, Citrate cycle (TCA cycle), One carbon pool by folate. Thus, the transcriptomic changes caused by neuron-specific Gclc overexpression in the thorax were less pronounced than in the head and affected pathways also differed from previous results. Although these pathways don't belong to the canonical longevity pathways, we suggest that they could participate in the delay of aging of Gclc overexpressing flies.

Ribosomal DNA harbors an evolutionarily conserved clock of biological aging

The ribosomal DNA (rDNA) is the most evolutionarily conserved segment of the genome and gives origin to the nucleolus, an energy intensive nuclear organelle and major hub influencing myriad molecular processes from cellular metabolism to epigenetic states of the genome. The rDNA/nucleolus has been directly and mechanistically implicated in aging and longevity in organisms as diverse as yeasts, Drosophila, and humans. The rDNA is also a significant target of DNA methylation that silences supernumerary rDNA units and regulates nucleolar activity. Here, we introduce an age clock built exclusively with CpG methylation within the rDNA. The ribosomal clock is sufficient to accurately estimate individual age within species, is responsive to genetic and environmental interventions that modulate life-span, and operates across species as distant as humans, mice, and dogs. Further analyses revealed a significant excess of age-associated hypermethylation in the rDNA relative to other segments of the genome, and which forms the basis of the rDNA clock. Our observations identified an evolutionarily conserved marker of aging that is easily ascertained, grounded on nucleolar biology, and could serve as a universal marker to gauge individual age and response to interventions in humans as well as laboratory and wild organisms across a wide diversity of species.

A 24 h Age Difference Causes Twice as Much Gene Expression Divergence as 100 Generations of Adaptation to a Novel Environment

Gene expression profiling is one of the most reliable high-throughput phenotyping methods, allowing researchers to quantify the transcript abundance of expressed genes. Because many biotic and abiotic factors influence gene expression, it is recommended to control them as tightly as possible. Here, we show that a 24 h age difference of Drosophilasimulans females that were subjected to RNA sequencing (RNA-Seq) five and six days after eclosure resulted in more than 2000 differentially expressed genes. This is twice the number of genes that changed expression during 100 generations of evolution in a novel hot laboratory environment. Importantly, most of the genes differing in expression due to age introduce false positives or negatives if an adaptive gene expression analysis is not controlled for age. Our results indicate that tightly controlled experimental conditions, including precise developmental staging, are needed for reliable gene expression analyses, in particular in an evolutionary framework.

Analysis of natural female post-mating responses of Anopheles gambiae and Anopheles coluzzii unravels similarities and differences in their reproductive ecology

Anopheles gambiae and An. coluzzii, the two most important malaria vectors in sub-Saharan Africa, are recently radiated sibling species that are reproductively isolated even in areas of sympatry. In females from these species, sexual transfer of male accessory gland products, including the steroid hormone 20-hydroxyecdysone (20E), induces vast behavioral, physiological, and transcriptional changes that profoundly shape their post-mating ecology, and that may have contributed to the insurgence of post-mating, prezygotic reproductive barriers. As these barriers can be detected by studying transcriptional changes induced by mating, we set out to analyze the post-mating response of An. gambiae and An. coluzzii females captured in natural mating swarms in Burkina Faso. While the molecular pathways shaping short- and long-term mating-induced changes are largely conserved in females from the two species, we unravel significant inter-specific differences that suggest divergent regulation of key reproductive processes such as egg development, processing of seminal secretion, and mating behavior, that may have played a role in reproductive isolation. Interestingly, a number of these changes occur in genes previously shown to be regulated by the sexual transfer of 20E and may be due to divergent utilization of this steroid hormone in the two species.

Gene expression hallmarks of cellular ageing

Ageing leads to dramatic changes in the physiology of many different tissues resulting in a spectrum of pathology. Nonetheless, many lines of evidence suggest that ageing is driven by highly conserved cell intrinsic processes, and a set of unifying hallmarks of ageing has been defined. Here, we survey reports of age-linked changes in basal gene expression across eukaryotes from yeast to human and identify six gene expression hallmarks of cellular ageing: downregulation of genes encoding mitochondrial proteins; downregulation of the protein synthesis machinery; dysregulation of immune system genes; reduced growth factor signalling; constitutive responses to stress and DNA damage; dysregulation of gene expression and mRNA processing. These encompass widely reported features of ageing such as increased senescence and inflammation, reduced electron transport chain activity and reduced ribosome synthesis, but also reveal a surprising lack of gene expression responses to known age-linked cellular stresses. We discuss how the existence of conserved transcriptomic hallmarks relates to genome-wide epigenetic differences underlying ageing clocks, and how the changing transcriptome results in proteomic alterations where data is available and to variations in cell physiology characteristic of ageing. Identification of gene expression events that occur during ageing across distant organisms should be informative as to conserved underlying mechanisms of ageing, and provide additional biomarkers to assess the effects of diet and other environmental factors on the rate of ageing.

The interplay between immunity and aging in Drosophila

Here, we provide a brief review of the mechanistic connections between immunity and aging—a fundamental biological relationship that remains poorly understood—by considering two intertwined questions: how does aging affect immunity, and how does immunity affect aging? On the one hand, aging contributes to the deterioration of immune function and predisposes the organism to infections (“immuno-senescence”). On the other hand, excessive activation of the immune system can accelerate degenerative processes, cause inflammation and immunopathology, and thus promote aging (“inflammaging”). Interestingly, several recent lines of evidence support the hypothesis that restrained or curbed immune activity at old age (that is, optimized age-dependent immune homeostasis) might actually improve realized immune function and thereby promote longevity. We focus mainly on insights from Drosophila, a powerful genetic model system in which both immunity and aging have been extensively studied, and conclude by outlining several unresolved questions in the field.

Genetic analysis of variation in lifespan using a multiparental advanced intercross Drosophila mapping population

BACKGROUND

Considerable natural variation for lifespan exists within human and animal populations. Genetically dissecting this variation can elucidate the pathways and genes involved in aging, and help uncover the genetic mechanisms underlying risk for age-related diseases. Studying aging in model systems is attractive due to their relatively short lifespan, and the ability to carry out programmed crosses under environmentally-controlled conditions. Here we investigate the genetic architecture of lifespan using the Drosophila Synthetic Population Resource (DSPR), a multiparental advanced intercross mapping population.

RESULTS

We measured lifespan in females from 805 DSPR lines, mapping five QTL (Quantitative Trait Loci) that each contribute 4-5 % to among-line lifespan variation in the DSPR. Each of these QTL co-localizes with the position of at least one QTL mapped in 13 previous studies of lifespan variation in flies. However, given that these studies implicate >90 % of the genome in the control of lifespan, this level of overlap is unsurprising. DSPR QTL intervals harbor 11-155 protein-coding genes, and we used RNAseq on samples of young and old flies to help resolve pathways affecting lifespan, and identify potentially causative loci present within mapped QTL intervals. Broad age-related patterns of expression revealed by these data recapitulate results from previous work. For example, we see an increase in antimicrobial defense gene expression with age, and a decrease in expression of genes involved in the electron transport chain. Several genes within QTL intervals are highlighted by our RNAseq data, such as Relish, a critical immune response gene, that shows increased expression with age, and UQCR-14, a gene involved in mitochondrial electron transport, that has reduced expression in older flies.

CONCLUSIONS

The five QTL we isolate collectively explain a considerable fraction of the genetic variation for female lifespan in the DSPR, and implicate modest numbers of genes. In several cases the candidate loci we highlight reside in biological pathways already implicated in the control of lifespan variation. Thus, our results provide further evidence that functional genetics tests targeting these genes will be fruitful, lead to the identification of natural sequence variants contributing to lifespan variation, and help uncover the mechanisms of aging.

A dataset for assessing temporal changes in gene expression during the aging process of adult Drosophila melanogaster

A Drosophila melanogaster genome-wide transcriptome dataset is available for studies on temporal patterns of gene expression. Gene expression was measured using two-dye color oligonucleotide arrays derived from Version 2 of the Drosophila Genomics Resource Center. A total of 15,158 oligonucleotide probes corresponded to a high proportion of the coding genes in the genome. The source of the flies was a highly genetically heterogeneous population maintained in an overlapping generation population regime. This regime was designed to maintain life history traits so that they were similar to those found in natural populations. Flies collected for the cohorts were obtained in a short period of time in a carefully controlled manner before virgin females and males were allowed to mate. Mated females were introduced into two large population cages in unusually high numbers (approximately 12,000 per cage) for a Drosophila laboratory longevity study. Samples were taken weekly from each cohort for 11 weeks; only a small proportion of surviving flies were present at the last two collection time points and thus they were exceptionally old compared to those collected in early-to-midlife samples. The data set is useful for studies of temporal patterns of gene expression as flies age. The very large size of each cohort, and relatively frequent incidence of temporal samples, allows for a fine-scale study of gene expression from young to very old flies.

Slowed aging during reproductive dormancy is reflected in genome-wide transcriptome changes in Drosophila melanogaster

Background

In models extensively used in studies of aging and extended lifespan, such as C. elegans and Drosophila, adult senescence is regulated by gene networks that are likely to be similar to ones that underlie lifespan extension during dormancy. These include the evolutionarily conserved insulin/IGF, TOR and germ line-signaling pathways. Dormancy, also known as dauer stage in the larval worm or adult diapause in the fly, is triggered by adverse environmental conditions, and results in drastically extended lifespan with negligible senescence. It is furthermore characterized by increased stress resistance and somatic maintenance, developmental arrest and reallocated energy resources. In the fly Drosophila melanogaster adult reproductive diapause is additionally manifested in arrested ovary development, improved immune defense and altered metabolism. However, the molecular mechanisms behind this adaptive lifespan extension are not well understood.

Results

A genome wide analysis of transcript changes in diapausing D. melanogaster revealed a differential regulation of more than 4600 genes. Gene ontology (GO) and KEGG pathway analysis reveal that many of these genes are part of signaling pathways that regulate metabolism, stress responses, detoxification, immunity, protein synthesis and processes during aging. More specifically, gene readouts and detailed mapping of the pathways indicate downregulation of insulin-IGF (IIS), target of rapamycin (TOR) and MAP kinase signaling, whereas Toll-dependent immune signaling, Jun-N-terminal kinase (JNK) and Janus kinase/signal transducer and activator of transcription (JAK/STAT) pathways are upregulated during diapause. Furthermore, we detected transcriptional regulation of a large number of genes specifically associated with aging and longevity.

Conclusions

We find that many affected genes and signal pathways are shared between dormancy, aging and lifespan extension, including IIS, TOR, JAK/STAT and JNK. A substantial fraction of the genes affected by diapause have also been found to alter their expression in response to starvation and cold exposure in D. melanogaster, and the pathways overlap those reported in GO analysis of other invertebrates in dormancy or even hibernating mammals. Our study, thus, shows that D. melanogaster is a genetically tractable model for dormancy in other organisms and effects of dormancy on aging and lifespan.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2383-1) contains supplementary material, which is available to authorized users.