Network analyses unveil ageing-associated pathways evolutionarily conserved from fungi to animals.

The genetic roots of the diverse paces and shapes of ageing and of the large variations in longevity observed across the tree of life are poorly understood. Indeed, pathways associated with ageing/longevity are incompletely known, both in terms of their constitutive genes/proteins and of their molecular interactions. Moreover, there is limited overlap between the genes constituting these pathways across mammals. Yet, dedicated comparative analyses might still unravel evolutionarily conserved, important pathways associated with longevity or ageing. Here, we used an original strategy with a double evolutionary and systemic focus to analyse protein interactions associated with ageing or longevity during the evolution of five species of Opisthokonta. We ranked these proteins and interactions based on their evolutionary conservation and centrality in past and present protein-protein interaction (PPI) networks, providing a big systemic picture of the evolution of ageing and longevity pathways that identified which pathways emerged in which Opisthokonta lineages, were conserved, and/or central. We confirmed that longevity/ageing-associated proteins (LAPs), be they pro- or anti-longevity, are highly central in extant PPI, consistently with the antagonistic pleiotropy theory of ageing, and identified key antagonistic regulators of ageing/longevity, 52 of which with homologues in humans. While some highly central LAPs were evolutionarily conserved for over a billion years, we report a clear transition in the functionally important components of ageing/longevity within bilaterians. We also predicted 487 novel evolutionarily conserved LAPs in humans, 54% of which are more central than mTOR, and 138 of which are druggable, defining new potential targets for anti-ageing treatments in humans.

Aging at Evolutionary Crossroads: Longitudinal Gene Co-expression Network Analyses of Proximal and Ultimate Causes of Aging in Bats.

How, when, and why do organisms, their tissues, and their cells age remain challenging issues, although researchers have identified multiple mechanistic causes of aging, and three major evolutionary theories have been developed to unravel the ultimate causes of organismal aging. A central hypothesis of these theories is that the strength of natural selection decreases with age. However, empirical evidence on when, why, and how organisms age is phylogenetically limited, especially in natural populations. Here, we developed generic comparisons of gene co-expression networks that quantify and dissect the heterogeneity of gene co-expression in conspecific individuals from different age-classes to provide topological evidence about some mechanical and fundamental causes of organismal aging. We applied this approach to investigate the complexity of some proximal and ultimate causes of aging phenotypes in a natural population of the greater mouse-eared bat Myotis myotis, a remarkably long-lived species given its body size and metabolic rate, with available longitudinal blood transcriptomes. M. myotis gene co-expression networks become increasingly fragmented with age, suggesting an erosion of the strength of natural selection and a general dysregulation of gene co-expression in aging bats. However, selective pressures remain sufficiently strong to allow successive emergence of homogeneous age-specific gene co-expression patterns, for at least 7 years. Thus, older individuals from long-lived species appear to sit at an evolutionary crossroad: as they age, they experience both a decrease in the strength of natural selection and a targeted selection for very specific biological processes, further inviting to refine a central hypothesis in evolutionary aging theories.

The Epistemic Revolution Induced by Microbiome Studies: An Interdisciplinary View.

Many separate fields and practices nowadays consider microbes as part of their legitimate focus. Therefore, microbiome studies may act as unexpected unifying forces across very different disciplines. Here, we summarize how microbiomes appear as novel major biological players, offer new artistic frontiers, new uses from medicine to laws, and inspire novel ontologies. We identify several convergent emerging themes across ecosystem studies, microbial and evolutionary ecology, arts, medicine, forensic analyses, law and philosophy of science, as well as some outstanding issues raised by microbiome studies across these disciplines and practices. An 'epistemic revolution induced by microbiome studies' seems to be ongoing, characterized by four features: (i) an ecologization of pre-existing concepts within disciplines, (ii) a growing interest in systemic analyses of the investigated or represented phenomena and a greater focus on interactions as their root causes, (iii) the intent to use openly multi-scalar interaction networks as an explanatory framework to investigate phenomena to acknowledge the causal effects of microbiomes, (iv) a reconceptualization of the usual definitions of which individuals are worth considering as an explanans or as an explanandum by a given field, which result in a fifth strong trend, namely (v) a de-anthropocentrification of our perception of the world.

Introducing Trait Networks to Elucidate the Fluidity of Organismal Evolution Using Palaeontological Data.

Explaining the evolution of animals requires ecological, developmental, paleontological, and phylogenetic considerations because organismal traits are affected by complex evolutionary processes. Modeling a plurality of processes, operating at distinct time-scales on potentially interdependent traits, can benefit from approaches that are complementary treatments to phylogenetics. Here, we developed an inclusive network approach, implemented in the command line software ComponentGrapher, and analyzed trait co-occurrence of rhinocerotoid mammals. We identified stable, unstable, and pivotal traits, as well as traits contributing to complexes, that may follow to a common developmental regulation, that point to an early implementation of the postcranial Bauplan among rhinocerotoids. Strikingly, most identified traits are highly dissociable, used repeatedly in distinct combinations and in different taxa, which usually do not form clades. Therefore, the genes encoding these traits are likely recruited into novel gene regulation networks during the course of evolution. Our evo-systemic framework, generalizable to other evolved organizations, supports a pluralistic modeling of organismal evolution, including trees and networks.

The Methodology Behind Network Thinking: Graphs to Analyze Microbial Complexity and Evolution.

In the post genomic era, large and complex molecular datasets from genome and metagenome sequencing projects expand the limits of what is possible for bioinformatic analyses. Network-based methods are increasingly used to complement phylogenetic analysis in studies in molecular evolution, including comparative genomics, classification, and ecological studies. Using network methods, the vertical and horizontal relationships between all genes or genomes, whether they are from cellular chromosomes or mobile genetic elements, can be explored in a single expandable graph. In recent years, development of new methods for the construction and analysis of networks has helped to broaden the availability of these approaches from programmers to a diversity of users. This chapter introduces the different kinds of networks based on sequence similarity that are already available to tackle a wide range of biological questions, including sequence similarity networks, gene-sharing networks and bipartite graphs, and a guide for their construction and analyses.

Comparative Transcriptomics Highlights New Features of the Iron Starvation Response in the Human Pathogen Candida glabrata.

In this work, we used comparative transcriptomics to identify regulatory outliers (ROs) in the human pathogen Candida glabrata. ROs are genes that have very different expression patterns compared to their orthologs in other species. From comparative transcriptome analyses of the response of eight yeast species to toxic doses of selenite, a pleiotropic stress inducer, we identified 38 ROs in C. glabrata. Using transcriptome analyses of C. glabrata response to five different stresses, we pointed out five ROs which were more particularly responsive to iron starvation, a process which is very important for C. glabrata virulence. Global chromatin Immunoprecipitation and gene profiling analyses showed that four of these genes are actually new targets of the iron starvation responsive Aft2 transcription factor in C. glabrata. Two of them (HBS1 and DOM34b) are required for C. glabrata optimal growth in iron limited conditions. In S. cerevisiae, the orthologs of these two genes are involved in ribosome rescue by the NO GO decay (NGD) pathway. Hence, our results suggest a specific contribution of NGD co-factors to the C. glabrata adaptation to iron starvation.

MultiTwin: A Software Suite to Analyze Evolution at Multiple Levels of Organization Using Multipartite Graphs.

The inclusion of introgressive processes in evolutionary studies induces a less constrained view of evolution. Network-based methods (like large-scale similarity networks) allow to include in comparative genomics all extrachromosomic carriers (like viruses, the most abundant biological entities on the planet) with their cellular hosts. The integration of several levels of biological organization (genes, genomes, communities, environments) enables more comprehensive analyses of gene sharing and improved sequence-based classifications. However, the algorithmic tools for the analysis of such networks are usually restricted to people with high programming skills. We present an integrated suite of software tools named MultiTwin, aimed at the construction, structuring, and analysis of multipartite graphs for evolutionary biology. Typically, this kind of graph is useful for the comparative analysis of the gene content of genomes in microbial communities from the environment and for exploring patterns of gene sharing, for example between distantly related cellular genomes, pangenomes, or between cellular genomes and their mobile genetic elements. We illustrate the use of this tool with an application of the bipartite approach (using gene family-genome graphs) for the analysis of pathogenicity traits in prokaryotes.

Tracking the Rules of Transmission and Introgression with Networks.

Understanding how an animal organism and its gut microbes form an integrated biological organization, known as a holobiont, is becoming a central issue in biological studies. Such an organization inevitably involves a complex web of transmission processes that occur on different scales in time and space, across microbes and hosts. Network-based models are introduced in this chapter to tackle aspects of this complexity and to better take into account vertical and horizontal dimensions of transmission. Two types of network-based models are presented, sequence similarity networks and bipartite graphs. One interest of these networks is that they can consider a rich diversity of important players in microbial evolution that are usually excluded from evolutionary studies, like plasmids and viruses. These methods bring forward the notion of "gene externalization," which is defined as the presence of redundant copies of prokaryotic genes on mobile genetic elements (MGEs), and therefore emphasizes a related although distinct process from lateral gene transfer between microbial cells. This chapter introduces guidelines to the construction of these networks, reviews their analysis, and illustrates their possible biological interpretations and uses. The application to human gut microbiomes shows that sequences present in a higher diversity of MGEs have both biased functions and a broader microbial and human host range. These results suggest that an "externalized gut metagenome" is partly common to humans and benefits the gut microbial community. We conclude that testing relationships between microbial genes, microbes, and their animal hosts, using network-based methods, could help to unravel additional mechanisms of transmission in holobionts.

Bipartite Network Analysis of Gene Sharings in the Microbial World.

Extensive microbial gene flows affect how we understand virology, microbiology, medical sciences, genetic modification, and evolutionary biology. Phylogenies only provide a narrow view of these gene flows: plasmids and viruses, lacking core genes, cannot be attached to cellular life on phylogenetic trees. Yet viruses and plasmids have a major impact on cellular evolution, affecting both the gene content and the dynamics of microbial communities. Using bipartite graphs that connect up to 149,000 clusters of homologous genes with 8,217 related and unrelated genomes, we can in particular show patterns of gene sharing that do not map neatly with the organismal phylogeny. Homologous genes are recycled by lateral gene transfer, and multiple copies of homologous genes are carried by otherwise completely unrelated (and possibly nested) genomes, that is, viruses, plasmids and prokaryotes. When a homologous gene is present on at least one plasmid or virus and at least one chromosome, a process of "gene externalization," affected by a postprocessed selected functional bias, takes place, especially in Bacteria. Bipartite graphs give us a view of vertical and horizontal gene flow beyond classic taxonomy on a single very large, analytically tractable, graph that goes beyond the cellular Web of Life.

Bipartite graph analyses reveal interdomain LGT involving ultrasmall prokaryotes and their divergent, membrane-related proteins.

Based on their small size and genomic properties, ultrasmall prokaryotic groups like the Candidate Phyla Radiation have been proposed as possible symbionts dependent on other bacteria or archaea. In this study, we use a bipartite graph analysis to examine patterns of sequence similarity between draft and complete genomes from ultrasmall bacteria and other complete prokaryotic genomes, assessing whether the former group might engage in significant gene transfer (or even endosymbioses) with other community members. Our results provide preliminary evidence for many lateral gene transfers with other prokaryotes, including members of the archaea, and report the presence of divergent, membrane-associated proteins among these ultrasmall taxa. In particular, these divergent genes were found in TM6 relatives of the intracellular parasite Babela massiliensis.

Network-Thinking: Graphs to Analyze Microbial Complexity and Evolution.

The tree model and tree-based methods have played a major, fruitful role in evolutionary studies. However, with the increasing realization of the quantitative and qualitative importance of reticulate evolutionary processes, affecting all levels of biological organization, complementary network-based models and methods are now flourishing, inviting evolutionary biology to experience a network-thinking era. We show how relatively recent comers in this field of study, that is, sequence-similarity networks, genome networks, and gene families-genomes bipartite graphs, already allow for a significantly enhanced usage of molecular datasets in comparative studies. Analyses of these networks provide tools for tackling a multitude of complex phenomena, including the evolution of gene transfer, composite genes and genomes, evolutionary transitions, and holobionts.

Automatic detection of anchor points for multiple sequence alignment.

BACKGROUND: determining beforehand specific positions to align (anchor points) has proved valuable for the accuracy of automated multiple sequence alignment (MSA) software. This feature can be used manually to include biological expertise, or automatically, usually by pairwise similarity searches. Multiple local similarities are be expected to be more adequate, as more biologically relevant. However, even good multiple local similarities can prove incompatible with the ordering of an alignment. RESULTS: we use a recently developed algorithm to detect multiple local similarities, which returns subsets of positions in the sequences sharing similar contexts of appearence. In this paper, we describe first how to get, with the help of this method, subsets of positions that could form partial columns in an alignment. We introduce next a graph-theoretic algorithm to detect (and remove) positions in the partial columns that are inconsistent with a multiple alignment. Partial columns can be used, for the time being, as guide only by a few MSA programs: ClustalW 2.0, DIALIGN 2 and T-Coffee. We perform tests on the effect of introducing these columns on the popular benchmark BAliBASE 3. CONCLUSIONS: we show that the inclusion of our partial alignment columns, as anchor points, improve on the whole the accuracy of the aligner ClustalW on the benchmark BAliBASE 3.

MS4--Multi-Scale Selector of Sequence Signatures: an alignment-free method for classification of biological sequences.

BACKGROUND: While multiple alignment is the first step of usual classification schemes for biological sequences, alignment-free methods are being increasingly used as alternatives when multiple alignments fail. Subword-based combinatorial methods are popular for their low algorithmic complexity (suffix trees ...) or exhaustivity (motif search), in general with fixed length word and/or number of mismatches. We developed previously a method to detect local similarities (the N-local decoding) based on the occurrences of repeated subwords of fixed length, which does not impose a fixed number of mismatches. The resulting similarities are, for some "good" values of N, sufficiently relevant to form the basis of a reliable alignment-free classification. The aim of this paper is to develop a method that uses the similarities detected by N-local decoding while not imposing a fixed value of N. We present a procedure that selects for every position in the sequences an adaptive value of N, and we implement it as the MS4 classification tool. RESULTS: Among the equivalence classes produced by the N-local decodings for all N, we select a (relatively) small number of "relevant" classes corresponding to variable length subwords that carry enough information to perform the classification. The parameter N, for which correct values are data-dependent and thus hard to guess, is here replaced by the average repetitivity kappa of the sequences. We show that our approach yields classifications of several sets of HIV/SIV sequences that agree with the accepted taxonomy, even on usually discarded repetitive regions (like the non-coding part of LTR). CONCLUSIONS: The method MS4 satisfactorily classifies a set of sequences that are notoriously hard to align. This suggests that our approach forms the basis of a reliable alignment-free classification tool. The only parameter kappa of MS4 seems to give reasonable results even for its default value, which can be a great advantage for sequence sets for which little information is available.

DIALIGN-TX and multiple protein alignment using secondary structure information at GOBICS.

We introduce web interfaces for two recent extensions of the multiple-alignment program DIALIGN. DIALIGN-TX combines the greedy heuristic previously used in DIALIGN with a more traditional 'progressive' approach for improved performance on locally and globally related sequence sets. In addition, we offer a version of DIALIGN that uses predicted protein secondary structures together with primary sequence information to construct multiple protein alignments. Both programs are available through 'Göttingen Bioinformatics Compute Server' (GOBICS).

A min-cut algorithm for the consistency problem in multiple sequence alignment.

MOTIVATION: Multiple sequence alignments can be constructed on the basis of pairwise local sequence similarities. This approach is rather flexible and can combine the advantages of global and local alignment methods. The restriction to pairwise alignments as building blocks, however, can lead to misalignments since weak homologies may be missed if only pairs of sequences are compared. RESULTS: Herein, we propose a graph-theoretical approach to find local multiple sequence similarities. Starting with pairwise alignments produced by DIALIGN, we use a min-cut algorithm to find potential (partial) alignment columns that we use to construct a final multiple alignment. On real and simulated benchmark data, our approach consistently outperforms the standard version of DIALIGN where local pairwise alignments are greedily incorporated into a multiple alignment. AVAILABILITY: The prototype is freely available under GNU Public Licence from E.C.