New genes from non-coding sequence: the role of de novo protein-coding genes in eukaryotic evolutionary innovation

The origin of novel protein-coding genes de novo was once considered so improbable as to be impossible. In less than a decade, and especially in the last five years, this view has been overturned by extensive evidence from diverse eukaryotic lineages. There is now evidence that this mechanism has contributed a significant number of genes to genomes of organisms as diverse as Saccharomyces, Drosophila, Plasmodium, Arabidopisis and human. From simple beginnings, these genes have in some instances acquired complex structure, regulated expression and important functional roles. New genes are often thought of as dispensable late additions; however, some recent de novo genes in human can play a role in disease. Rather than an extremely rare occurrence, it is now evident that there is a relatively constant trickle of proto-genes released into the testing ground of natural selection. It is currently unknown whether de novo genes arise primarily through an ‘RNA-first’ or ‘ORF-first’ pathway. Either way, evolutionary tinkering with this pool of genetic potential may have been a significant player in the origins of lineage-specific traits and adaptations.     

GeneCensus: genome comparisons in terms of metabolic pathway activity and protein family sharing

We present a prototype of a new database tool, GeneCensus, which focuses on comparing genomes globally, in terms of the collective properties of many genes, rather than in terms of the attributes of a single gene (e.g. sequence similarity for a particular ortholog). The comparisons are presented in a visual fashion over the web at GeneCensus.org. The system concentrates on two types of comparisons: (i) trees based on the sharing of generalized protein families between genomes, and (ii) whole pathway analysis in terms of activity levels. For the trees, we have developed a module (TreeViewer) that clusters genomes in terms of the folds, superfamilies or orthologs—all can be considered as generalized ‘families’ or ‘protein parts’—they share, and compares the resulting trees side-by-side with those built from sequence similarity of individual genes (e.g. a traditional tree built on ribosomal similarity). We also include comparisons to trees built on whole-genome dinucleotide or codon composition. For pathway comparisons, we have implemented a module (PathwayPainter) that graphically depicts, in selected metabolic pathways, the fluxes or expression levels of the associated enzymes (i.e. generalized ‘activities’). One can, consequently, compare organisms (and organism states) in terms of representations of these systemic quantities. Develop ment of this module involved compiling, calculating and standardizing flux and expression information from many different sources. We illustrate pathway analysis for enzymes involved in central metabolism. We are able to show that, to some degree, flux and expression fluctuations have characteristic values in different sections of the central metabolism and that control points in this system (e.g. hexokinase, pyruvate kinase, phosphofructokinase, isocitrate dehydrogenase and citric synthase) tend to be especially variable in flux and expression. Both the TreeViewer and PathwayPainter modules connect to other information sources related to individual-gene or organism properties (e.g. a single-gene structural annotation viewer).     

A dynamic model for replication protein A (RPA) function in DNA processing pathways

Processing of DNA in replication, repair and recombination pathways in cells of all organisms requires the participation of at least one major single-stranded DNA (ssDNA)-binding protein. This protein protects ssDNA from nucleolytic damage, prevents hairpin formation and blocks DNA reannealing until the processing pathway is successfully completed. Many ssDNA-binding proteins interact physically and functionally with a variety of other DNA processing proteins. These interactions are thought to temporally order and guide the parade of proteins that ‘trade places’ on the ssDNA, a model known as ‘hand-off’, as the processing pathway progresses. How this hand-off mechanism works remains poorly understood. Recent studies of the conserved eukaryotic ssDNA-binding protein replication protein A (RPA) suggest a novel mechanism by which proteins may trade places on ssDNA by binding to RPA and mediating conformation changes that alter the ssDNA-binding properties of RPA. This article reviews the structure and function of RPA, summarizes recent studies of RPA in DNA replication and other DNA processing pathways, and proposes a general model for the role of RPA in protein-mediated hand-off.     

The Effects of Extra-Somatic Weapons on the Evolution of Human Cooperation towards Non-Kin

Human cooperation and altruism towards non-kin is a major evolutionary puzzle, as is ‘strong reciprocity’ where no present or future rewards accrue to the co-operator/altruist. Here, we test the hypothesis that the development of extra-somatic weapons could have influenced the evolution of human cooperative behaviour, thus providing a new explanation for these two puzzles. Widespread weapons use could have made disputes within hominin groups far more lethal and also equalized power between individuals. In such a cultural niche non-cooperators might well have become involved in such lethal disputes at a higher frequency than cooperators, thereby increasing the relative fitness of genes associated with cooperative behaviour. We employ two versions of the evolutionary Iterated Prisoner''s Dilemma (IPD) model – one where weapons use is simulated and one where it is not. We then measured the performance of 25 IPD strategies to evaluate the effects of weapons use on them. We found that cooperative strategies performed significantly better, and non-cooperative strategies significantly worse, under simulated weapons use. Importantly, the performance of an ‘Always Cooperate’ IPD strategy, equivalent to that of ‘strong reciprocity’, improved significantly more than that of all other cooperative strategies. We conclude that the development of extra-somatic weapons throws new light on the evolution of human altruistic and cooperative behaviour, and particularly ‘strong reciprocity’. The notion that distinctively human altruism and cooperation could have been an adaptive trait in a past environment that is no longer evident in the modern world provides a novel addition to theory that seeks to account for this major evolutionary puzzle.     

Cinnamic Aldehyde,the main monomer component of Cinnamon,exhibits anti‐inflammatory property in OA synovial fibroblasts via TLR4/MyD88 pathway

Cinnamon is a wildly used traditional Chinese herbal medicine for osteoarthritis (OA) treatment, but the underlying mechanism remains ambiguous. The purpose of this study is to explore the mechanism of cinnamic aldehyde (CA), a bioactive substance extracted from Cinnamon, on synovial inflammation in OA. A total of 144 CA‐OA co‐targeted genes were identified by detect databases (PubChem, HIT, TCMSP, TTD, DrugBank and GeneCards). The results of GO enrichment analysis indicated that these co‐targeted genes have participated in many biological processes including ‘inflammatory response’, ‘cellular response to lipopolysaccharide’, ‘response to drug’, ‘immune response’, ‘lipopolysaccharide‐mediated signalling pathway’, etc. KEGG pathway analysis showed these co‐targeted genes were mainly enriched in ‘Toll‐like receptor signalling pathway’, ‘TNF signalling pathway’, ‘NF‐kappa B signalling pathway’, etc. Molecular docking demonstrated that CA could successfully bind to TLR2 and TLR4. The results of in vitro experiments showed no potential toxicity of 10, 20 and 50 μM/L CA on human OA FLS, and CA can significantly inhibit the inflammation in LPS‐induced human FLS. Further experimental mechanism evidence confirmed CA can inhibited the inflammation in LPS‐induced human OA FLS via blocking the TLR4/MyD88 signalling pathway. Our results demonstrated that CA exhibited strong anti‐inflammation effect in OA FLS through blocking the activation of TLR4/MyD88 signalling pathway, suggesting its potential as a hopeful candidate for the development of novel agents for the treatment of OA.     

Development of a prognostic model based on an immunogenomic landscape analysis of medulloblastoma

Medulloblastoma (MB) is one of the most common central nervous system tumors in children. At present, the vital role of immune abnormalities has been proved in tumorigenesis and progression. However, the immune mechanism in MB is still poorly understood. In the present study, 51 differentially expressed immune-related genes (DE-IRGs) and 226 survival associated immune-related genes (Sur-IRGs) were screened by an integrated analysis of multi-array. Moreover, the potential pathways were enriched by functional analysis, such as ‘cytokine–cytokine receptor interaction’, ‘Ras signaling pathway’, ‘PI3K-Akt signaling pathway’ and ‘pathways in cancer’. Furthermore, 10 core IRGs were identified from DE-IRGs and Sur-IRGs. And the potential regulatory mechanisms of core IRGs were also explored. Additionally, a new prognostic model, including 7 genes (HDGF, CSK, PNOC, S100A13, RORB, FPR1, and ICAM2) based on IRGs, was established by multivariable COX analysis. In summary, our study revealed the underlying immune mechanism of MB. Moreover, we developed a prognostic model associated with clinical characteristics and could reflect the infiltration of immune cells.     

Paleovirology of ‘syncytins’, retroviral env genes exapted for a role in placentation

The development of the emerging field of ‘paleovirology’ allows biologists to reconstruct the evolutionary history of fossil endogenous retroviral sequences integrated within the genome of living organisms and has led to the retrieval of conserved, ancient retroviral genes ‘exapted’ by ancestral hosts to fulfil essential physiological roles, syncytin genes being undoubtedly among the most remarkable examples of such a phenomenon. Indeed, syncytins are ‘new’ genes encoding proteins derived from the envelope protein of endogenous retroviral elements that have been captured and domesticated on multiple occasions and independently in diverse mammalian species, through a process of convergent evolution. Knockout of syncytin genes in mice provided evidence for their absolute requirement for placenta development and embryo survival, via formation by cell–cell fusion of syncytial cell layers at the fetal–maternal interface. These genes of exogenous origin, acquired ‘by chance’ and yet still ‘necessary’ to carry out a basic function in placental mammals, may have been pivotal in the emergence of mammalian ancestors with a placenta from egg-laying animals via the capture of a founding retroviral env gene, subsequently replaced in the diverse mammalian lineages by new env-derived syncytin genes, each providing its host with a positive selective advantage.     

When everything changes at once: finding a new normal after genome duplication

Whole-genome duplication (WGD), which leads to polyploidy, is implicated in adaptation and speciation. But what are the immediate effects of WGD and how do newly polyploid lineages adapt to them? With many studies of new and evolved polyploids now available, along with studies of genes under selection in polyploids, we are in an increasingly good position to understand how polyploidy generates novelty. Here, I will review consistent effects of WGD on the biology of plants, such as an increase in cell size, increased stress tolerance and more. I will discuss how a change in something as fundamental as cell size can challenge the function of some cell types in particular. I will also discuss what we have learned about the short- to medium-term evolutionary response to WGD. It is now clear that some of this evolutionary response may ‘lock in’ traits that happen to be beneficial, while in other cases, it might be more of an ‘emergency response’ to work around physiological changes that are either deleterious, or cannot be undone in the polyploid context. Yet, other traits may return rapidly to a diploid-like state. Polyploids may, by re-jigging many inter-related processes, find a new, conditionally adaptive, normal.     

Phylogenetic Analysis of Glycerol 3-Phosphate Acyltransferases in Opisthokonts Reveals Unexpected Ancestral Complexity and Novel Modern Biosynthetic Components

Glycerolipid synthesis represents a central metabolic process of all forms of life. In the last decade multiple genes coding for enzymes responsible for the first step of the pathway, catalyzed by glycerol 3-phosphate acyltransferase (GPAT), have been described, and characterized primarily in model organisms like Saccharomyces cerevisiae and mice. Notoriously, the fungal enzymes share low sequence identity with their known animal counterparts, and the nature of their homology is unclear. Furthermore, two mitochondrial GPAT isoforms have been described in animal cells, while no such enzymes have been identified in Fungi. In order to determine if the yeast and mammalian GPATs are representative of the set of enzymes present in their respective groups, and to test the hypothesis that metazoan orthologues are indeed absent from the fungal clade, a comparative genomic and phylogenetic analysis was performed including organisms spanning the breadth of the Opisthokonta supergroup. Surprisingly, our study unveiled the presence of ‘fungal’ orthologs in the basal taxa of the holozoa and ‘animal’ orthologues in the basal holomycetes. This includes a novel clade of fungal homologues, with putative peroxisomal targeting signals, of the mitochondrial/peroxisomal acyltransferases in Metazoa, thus potentially representing an undescribed metabolic capacity in the Fungi. The overall distribution of GPAT homologues is suggestive of high relative complexity in the ancestors of the opisthokont clade, followed by loss and sculpting of the complement in the descendent lineages. Divergence from a general versatile metabolic model, present in ancestrally deduced GPAT complements, points to distinctive contributions of each GPAT isoform to lipid metabolism and homeostasis in contemporary organisms like humans and their fungal pathogens.     

Probing small ribosomal subunit RNA helix 45 acetylation across eukaryotic evolution

NAT10 is an essential enzyme that catalyzes N⁴-acetylcytidine (ac⁴C) in eukaryotic transfer RNA and 18S ribosomal RNA. Recent studies suggested that rRNA acetylation is dependent on SNORD13, a box C/D small nucleolar RNA predicted to base-pair with 18S rRNA via two antisense elements. However, the selectivity of SNORD13-dependent cytidine acetylation and its relationship to NAT10’s essential function remain to be defined. Here, we demonstrate that SNORD13 is required for acetylation of a single cytidine of human and zebrafish 18S rRNA. In-depth characterization revealed that SNORD13-dependent ac⁴C is dispensable for human cell growth, ribosome biogenesis, translation and development. This loss of function analysis inspired a cross-evolutionary survey of the eukaryotic rRNA acetylation ‘machinery’ that led to the characterization of many novel metazoan SNORD13 genes. This includes an atypical SNORD13-like RNA in Drosophila melanogaster which guides ac⁴C to 18S rRNA helix 45 despite lacking one of the two rRNA antisense elements. Finally, we discover that Caenorhabditis elegans 18S rRNA is not acetylated despite the presence of an essential NAT10 homolog. Our findings shed light on the molecular mechanisms underlying SNORD13-mediated rRNA acetylation across eukaryotic evolution and raise new questions regarding the biological and evolutionary relevance of this highly conserved rRNA modification.     

Bridging the membrane lipid divide: bacteria of the FCB group superphylum have the potential to synthesize archaeal ether lipids

Archaea synthesize membranes of isoprenoid lipids that are ether-linked to glycerol-1-phosphate (G1P), while Bacteria/Eukarya produce membranes consisting of fatty acids ester-bound to glycerol-3-phosphate (G3P). This dichotomy in membrane lipid composition (i.e., the ‘lipid divide’) is believed to have arisen after the Last Universal Common Ancestor (LUCA). A leading hypothesis is that LUCA possessed a heterochiral ‘mixed archaeal/bacterial membrane’. However, no natural microbial representatives supporting this scenario have been shown to exist today. Here, we demonstrate that bacteria of the Fibrobacteres–Chlorobi–Bacteroidetes (FCB) group superphylum encode a putative archaeal pathway for ether-bound isoprenoid membrane lipids in addition to the bacterial fatty acid membrane pathway. Key genes were expressed in the environment and their recombinant expression in Escherichia coli resulted in the formation of a ‘mixed archaeal/bacterial membrane’. Genomic evidence and biochemical assays suggest that the archaeal-like lipids of members of the FCB group could possess either a G1P or G3P stereochemistry. Our results support the existence of ‘mixed membranes’ in natural environments and their stability over a long period in evolutionary history, thereby bridging a once-thought fundamental divide in biology.Subject terms: Water microbiology, Cellular microbiology, Molecular evolution, Metagenomics     

Systems-Wide Prediction of Enzyme Promiscuity Reveals a New Underground Alternative Route for Pyridoxal 5’-Phosphate Production in E. coli

Recent insights suggest that non-specific and/or promiscuous enzymes are common and active across life. Understanding the role of such enzymes is an important open question in biology. Here we develop a genome-wide method, PROPER, that uses a permissive PSI-BLAST approach to predict promiscuous activities of metabolic genes. Enzyme promiscuity is typically studied experimentally using multicopy suppression, in which over-expression of a promiscuous ‘replacer’ gene rescues lethality caused by inactivation of a ‘target’ gene. We use PROPER to predict multicopy suppression in Escherichia coli, achieving highly significant overlap with published cases (hypergeometric p = 4.4e-13). We then validate three novel predicted target-replacer gene pairs in new multicopy suppression experiments. We next go beyond PROPER and develop a network-based approach, GEM-PROPER, that integrates PROPER with genome-scale metabolic modeling to predict promiscuous replacements via alternative metabolic pathways. GEM-PROPER predicts a new indirect replacer (thiG) for an essential enzyme (pdxB) in production of pyridoxal 5’-phosphate (the active form of Vitamin B₆), which we validate experimentally via multicopy suppression. We perform a structural analysis of thiG to determine its potential promiscuous active site, which we validate experimentally by inactivating the pertaining residues and showing a loss of replacer activity. Thus, this study is a successful example where a computational investigation leads to a network-based identification of an indirect promiscuous replacement of a key metabolic enzyme, which would have been extremely difficult to identify directly.     

VarSAn: associating pathways with a set of genomic variants using network analysis

There is a pressing need today to mechanistically interpret sets of genomic variants associated with diseases. Here we present a tool called ‘VarSAn’ that uses a network analysis algorithm to identify pathways relevant to a given set of variants. VarSAn analyzes a configurable network whose nodes represent variants, genes and pathways, using a Random Walk with Restarts algorithm to rank pathways for relevance to the given variants, and reports P-values for pathway relevance. It treats non-coding and coding variants differently, properly accounts for the number of pathways impacted by each variant and identifies relevant pathways even if many variants do not directly impact genes of the pathway. We use VarSAn to identify pathways relevant to variants related to cancer and several other diseases, as well as drug response variation. We find VarSAn''s pathway ranking to be complementary to the standard approach of enrichment tests on genes related to the query set. We adopt a novel benchmarking strategy to quantify its advantage over this baseline approach. Finally, we use VarSAn to discover key pathways, including the VEGFA-VEGFR2 pathway, related to de novo variants in patients of Hypoplastic Left Heart Syndrome, a rare and severe congenital heart defect.     

Diversification of Two Lineages of Symbiotic Photobacterium

Understanding of processes driving bacterial speciation requires examination of closely related, recently diversified lineages. To gain an insight into diversification of bacteria, we conducted comparative genomic analysis of two lineages of bioluminescent symbionts, Photobacterium leiognathi and ‘P. mandapamensis’. The two lineages are evolutionary and ecologically closely related. Based on the methods used in bacterial taxonomy for classification of new species (DNA-DNA hybridization and ANI), genetic relatedness of the two lineages is at a cut-off point for species delineation. In this study, we obtained the whole genome sequence of a representative P. leiognathi strain lrivu.4.1, and compared it to the whole genome sequence of ‘P. mandapamensis’ svers.1.1. Results of the comparative genomic analysis suggest that P. leiognathi has a more plastic genome and acquired genes horizontally more frequently than ‘P. mandapamensis’. We predict that different rates of recombination and gene acquisition contributed to diversification of the two lineages. Analysis of lineage-specific sequences in 25 strains of P. leiognathi and ‘P. mandapamensis’ found no evidence that bioluminescent symbioses with specific host animals have played a role in diversification of the two lineages.     

Eco-evolutionary feedbacks,adaptive dynamics and evolutionary rescue theory

Adaptive dynamics theory has been devised to account for feedbacks between ecological and evolutionary processes. Doing so opens new dimensions to and raises new challenges about evolutionary rescue. Adaptive dynamics theory predicts that successive trait substitutions driven by eco-evolutionary feedbacks can gradually erode population size or growth rate, thus potentially raising the extinction risk. Even a single trait substitution can suffice to degrade population viability drastically at once and cause ‘evolutionary suicide’. In a changing environment, a population may track a viable evolutionary attractor that leads to evolutionary suicide, a phenomenon called ‘evolutionary trapping’. Evolutionary trapping and suicide are commonly observed in adaptive dynamics models in which the smooth variation of traits causes catastrophic changes in ecological state. In the face of trapping and suicide, evolutionary rescue requires that the population overcome evolutionary threats generated by the adaptive process itself. Evolutionary repellors play an important role in determining how variation in environmental conditions correlates with the occurrence of evolutionary trapping and suicide, and what evolutionary pathways rescue may follow. In contrast with standard predictions of evolutionary rescue theory, low genetic variation may attenuate the threat of evolutionary suicide and small population sizes may facilitate escape from evolutionary traps.     

Drug Resistance Missense Mutations in Cancer Are Subject to Evolutionary Constraints

Several tumour types are sensitive to deactivation of just one or very few genes that are constantly active in the cancer cells, a phenomenon that is termed ‘oncogene addiction’. Drugs that target the products of those oncogenes can yield a temporary relief, and even complete remission. Unfortunately, many patients receiving oncogene-targeted therapies relapse on treatment. This often happens due to somatic mutations in the oncogene (‘resistance mutations’). ‘Compound mutations’, which in the context of cancer drug resistance are defined as two or more mutations of the drug target in the same clone may lead to enhanced resistance against the most selective inhibitors. Here, it is shown that the vast majority of the resistance mutations occurring in cancer patients treated with tyrosin kinase inhibitors aimed at three different proteins follow an evolutionary pathway. Using bioinformatic analysis tools, it is found that the drug-resistance mutations in the tyrosine kinase domains of Abl1, ALK and exons 20 and 21 of EGFR favour transformations to residues that can be identified in similar positions in evolutionary related proteins. The results demonstrate that evolutionary pressure shapes the mutational landscape in the case of drug-resistance somatic mutations. The constraints on the mutational landscape suggest that it may be possible to counter single drug-resistance point mutations. The observation of relatively many resistance mutations in Abl1, but not in the other genes, is explained by the fact that mutations in Abl1 tend to be biochemically conservative, whereas mutations in EGFR and ALK tend to be radical. Analysis of Abl1 compound mutations suggests that such mutations are more prevalent than hitherto reported and may be more difficult to counter. This supports the notion that such mutations may provide an escape route for targeted cancer drug resistance.     

Compartmentation in plant metabolism

Cell fractionation and immunohistochemical studies in the last 40 years have revealed the extensive compartmentation of plant metabolism. In recent years, new protein mass spectrometry and fluorescent-protein tagging technologies have accelerated the flow of information, especially for Arabidopsis thaliana, but the intracellular locations of the majority of proteins in the plant proteome are still not known. Prediction programs that search for targeting information within protein sequences can be applied to whole proteomes, but predictions from different programs often do not agree with each other or, indeed, with experimentally determined results. The compartmentation of most pathways of primary metabolism is generally covered in plant physiology textbooks, so the focus here is mainly on newly discovered metabolic pathways in plants or pathways that have recently been revised. Ultimately, all of the pathways of plant metabolism are interconnected, and a major challenge facing plant biochemists is to understand the regulation and control of metabolic networks. One of the best-characterized networks links sucrose synthesis in the cytosol with photosynthetic CO(2) fixation and starch synthesis in the chloroplasts. One of the key features of this network is how the transport of pathway intermediates and signal metabolites across the chloroplast envelope conveys information between the two compartments, influencing the regulation of several enzymes to co-ordinate fluxes through the different pathways. It is widely accepted that chloroplasts and mitochondria originated from prokaryotic endosymbionts, and that new transporters and regulatory networks evolved to integrate metabolism in these organelles with the rest of the cell. Curiously, the present-day locations of many metabolic pathways within the cell often do not reflect their evolutionary origin, and there is evidence of extensive shuffling of enzymes and whole pathways between compartments during the evolution of plants.     

Network-Based Data Integration for Selecting Candidate Virulence Associated Proteins in the Cereal Infecting Fungus Fusarium graminearum

The identification of virulence genes in plant pathogenic fungi is important for understanding the infection process, host range and for developing control strategies. The analysis of already verified virulence genes in phytopathogenic fungi in the context of integrated functional networks can give clues about the underlying mechanisms and pathways directly or indirectly linked to fungal pathogenicity and can suggest new candidates for further experimental investigation, using a ‘guilt by association’ approach. Here we study 133 genes in the globally important Ascomycete fungus Fusarium graminearum that have been experimentally tested for their involvement in virulence. An integrated network that combines information from gene co-expression, predicted protein-protein interactions and sequence similarity was employed and, using 100 genes known to be required for virulence, we found a total of 215 new proteins potentially associated with virulence of which 29 are annotated as hypothetical proteins. The majority of these potential virulence genes are located in chromosomal regions known to have a low recombination frequency. We have also explored the taxonomic diversity of these candidates and found 25 sequences, which are likely to be fungal specific. We discuss the biological relevance of a few of the potentially novel virulence associated genes in detail. The analysis of already verified virulence genes in phytopathogenic fungi in the context of integrated functional networks can give clues about the underlying mechanisms and pathways directly or indirectly linked to fungal pathogenicity and can suggest new candidates for further experimental investigation, using a ‘guilt by association’ approach.