Biochemical evolution

Mitochondrial DNA marker: A PCR approach for tracking rat (Rattus rattus and Rattus norvegicus) fecal pollution in surface water systems

Rats act as reservoirs for a wide range of zoonotic pathogens and can negatively impact human health. In this study, we developed a novel dye base mitochondrial DNA (mtDNA) PCR-assay (RatMt) specifically targeting a 180 bp fragment of the NADH dehydrogenase subunit 2 gene for detecting fecal pollution from two species of rats (Rattus rattus and Rattus norvegicus) in environmental samples. Estimation of Escherichia coli concentrations in Rattus norvegicus fecal pellets suggested that there were approximately 2.24 × 10 4 ± 4.86 × 103 MPN/g of fecal pellet. The RatMt PCR assay was robust, had a detection limit of rat feces in water of 0.274 ± 0.14 mg/100 mL and was 100 % specific for detecting Rattus rattus and Rattus norvegicus fecal mtDNA. Fecal Indicator Bacteria (FIB) along an urbanized gradient in Pensacola-Bay was assessed by the IDEXX Colilert™ - 18 and indicated that the majority of the fifteen sampling sites in the Pensacola-Bay area had E. coli concentrations >410 MPN/100 mL. Rattus rattus and Rattus norvegicus mtDNA were detected in all the urban marine sites, three of the urban freshwater sampling areas, and three of the forested sampling sites. The RatMt PCR assay is a useful tool for rapidly detecting Rattus rattus and Rattus norvegicus fecal pollution in environmental samples.

When colors mislead: Genomics and bioacoustics prompt re-classification of Asian flycatcher radiation (Aves: Niltavinae)

Traditional classification of many animals, including birds, has been highly dependent on external morphological characters like plumage coloration. However, both bioacoustics and genetic or genomic data have revolutionized our understanding of the relationships of certain lineages and led to sweeping taxonomic re-organizations. In this study, we present a case of erroneous delimitation of genus boundaries in the species-rich flycatcher subfamily Niltavinae. Genera within this subfamily have historically been delineated based on blue versus brown male body plumage until recent studies based on a few mitochondrial and nuclear loci unearthed several cases of generic misclassification. Here we use extensive bioacoustic data from 43 species and genomic data from 28 species for a fundamental reclassification of species in the Niltavinae. Our study reveals that song is an important trait to classify these birds even at the genus level, whereas plumage traits exhibit ample convergence and have led to numerous historic misattributions. Our taxonomic re-organization leads to new biogeographic limits of major genera, such that the genus Cyornis now only extends as far east as the islands of Sulawesi, Sula, and Banggai, whereas Eumyias is redefined to extend far beyond Wallace's Line to the islands of Seram and Timor. Our conclusions advise against an over-reliance on morphological traits and underscore the importance of integrative datasets.

Moderating the neutralist-selectionist debate: exactly which propositions are we debating, and which arguments are valid?

Half a century after its foundation, the neutral theory of molecular evolution continues to attract controversy. The debate has been hampered by the coexistence of different interpretations of the core proposition of the neutral theory, the 'neutral mutation-random drift' hypothesis. In this review, we trace the origins of these ambiguities and suggest potential solutions. We highlight the difference between the original, the revised and the nearly neutral hypothesis, and re-emphasise that none of them equates to the null hypothesis of strict neutrality. We distinguish the neutral hypothesis of protein evolution, the main focus of the ongoing debate, from the neutral hypotheses of genomic and functional DNA evolution, which for many species are generally accepted. We advocate a further distinction between a narrow and an extended neutral hypothesis (of which the latter posits that random non-conservative amino acid substitutions can cause non-ecological phenotypic divergence), and we discuss the implications for evolutionary biology beyond the domain of molecular evolution. We furthermore point out that the debate has widened from its initial focus on point mutations, and also concerns the fitness effects of large-scale mutations, which can alter the dosage of genes and regulatory sequences. We evaluate the validity of neutralist and selectionist arguments and find that the tested predictions, apart from being sensitive to violation of underlying assumptions, are often derived from the null hypothesis of strict neutrality, or equally consistent with the opposing selectionist hypothesis, except when assuming molecular panselectionism. Our review aims to facilitate a constructive neutralist-selectionist debate, and thereby to contribute to answering a key question of evolutionary biology: what proportions of amino acid and nucleotide substitutions and polymorphisms are adaptive?

Evolutionary genetics of the mitochondrial genome: insights from Drosophila

Mitochondria are key to energy conversion in virtually all eukaryotes. Intriguingly, despite billions of years of evolution inside the eukaryote, mitochondria have retained their own small set of genes involved in the regulation of oxidative phosphorylation (OXPHOS) and protein translation. Although there was a long-standing assumption that the genetic variation found within the mitochondria would be selectively neutral, research over the past 3 decades has challenged this assumption. This research has provided novel insight into the genetic and evolutionary forces that shape mitochondrial evolution and broader implications for evolutionary ecological processes. Many of the seminal studies in this field, from the inception of the research field to current studies, have been conducted using Drosophila flies, thus establishing the species as a model system for studies in mitochondrial evolutionary biology. In this review, we comprehensively review these studies, from those focusing on genetic processes shaping evolution within the mitochondrial genome, to those examining the evolutionary implications of interactions between genes spanning mitochondrial and nuclear genomes, and to those investigating the dynamics of mitochondrial heteroplasmy. We synthesize the contribution of these studies to shaping our understanding of the evolutionary and ecological implications of mitochondrial genetic variation.

Complexity and self-organization in the evolution of cell polarization

Cellular life exhibits order and complexity, which typically increase over the course of evolution. Cell polarization is a well-studied example of an ordering process that breaks the internal symmetry of a cell by establishing a preferential axis. Like many cellular processes, polarization is driven by self-organization, meaning that the macroscopic pattern emerges as a consequence of microscopic molecular interactions at the biophysical level. However, the role of self-organization in the evolution of complex protein networks remains obscure. In this Review, we provide an overview of the evolution of polarization as a self-organizing process, focusing on the model species Saccharomyces cerevisiae and its fungal relatives. Moreover, we use this model system to discuss how self-organization might relate to evolutionary change, offering a shift in perspective on evolution at the microscopic scale.

Twentieth-Century Paleoproteomics: Lessons from Venta Micena Fossils

Simple Summary

Two independent research groups led by Olivares (Spain) and Lowenstein (USA) investigated the immunological reactions of proteins extracted from the controversial Orce skull (VM-0), a 1.3-million-year-old fossil found at the Venta Micena site in Orce, Granada (Spain) and initially believed to come from an unidentified hominin. Work by both groups with polyclonal and monoclonal antibodies showed that proteins from this fossil reacted most strongly to antibodies against modern human proteins. Other hominin and mammal fossils from Venta Micena were also studied.

Abstract

Proteomics methods can identify amino acid sequences in fossil proteins, thus making it possible to determine the ascription or proximity of a fossil to other species. Before mass spectrometry was used to study fossil proteins, earlier studies used antibodies to recognize their sequences. Lowenstein and colleagues, at the University of San Francisco, pioneered the identification of fossil proteins with immunological methods. His group, together with Olivares’s group at the University of Granada, studied the immunological reactions of proteins from the controversial Orce skull fragment (VM-0), a 1.3-million-year-old fossil found at the Venta Micena site in Orce (Granada province, southern Spain) and initially assigned to a hominin. However, discrepancies regarding the morphological features of the internal face of the fossil raised doubts about this ascription. In this article, we review the immunological analysis of the proteins extracted from VM-0 and other Venta Micena fossils assigned to hominins and to other mammals, and explain how these methods helped to determine the species specificity of these fossils and resolve paleontological controversies.

Global genomic and proteomic analysis indicates co-evolution of Neisseria species and with their human host

Neisseria, a genus from the beta-proteobacteria class, is of potential clinical importance. This genus contains both pathogenic and commensal strains. Gonorrhea and meningitis are two major diseases caused by pathogens belonging to this genus. With the increased use of antimicrobial agents against these pathogens they have evolved the antimicrobial resistance capacity making these diseases nearly untreatable. The set of anti-bacterial resistance genes (resistome) and genes associated with signal processing (secretomes) are crucial for the host-microbial interaction. With the virtue of whole-genome sequences and computational biology, it is now possible to study the genomic and proteomic riddles of Neisseria along with their comprehensive evolutionary and metabolic profiling. We have studied relative synonymous codon usage, amino acid usage, reverse ecology, comparative genomics, evolutionary analysis and pathogen-host (Neisseria-human) interaction through bioinformatics analysis. Our analysis revealed the co-evolution of Neisseria genomes with the human host. Moreover, the co-occurrence of Neisseria and humans has been supported through reverse ecology analysis. A differential pattern of the evolutionary rate of resistomes and secretomes was evident among the pathogenic and commensal strains. Comparative genomics supported the presence of virulent genes in both pathogenic and commensal strains of the select genus. Our analysis also indicated a transition from commensal to pathogenic Neisseria strains through the long run of evolution.

Comparative genome analysis among Variovorax species and genome guided aromatic compound degradation analysis emphasizing 4-hydroxybenzoate degradation in Variovorax sp. PAMC26660

Background

While the genus Variovorax is known for its aromatic compound metabolism, no detailed study of the peripheral and central pathways of aromatic compound degradation has yet been reported. Variovorax sp. PAMC26660 is a lichen-associated bacterium isolated from Antarctica. The work presents the genome-based elucidation of peripheral and central catabolic pathways of aromatic compound degradation genes in Variovorax sp. PAMC26660. Additionally, the accessory, core and unique genes were identified among Variovorax species using the pan genome analysis tool. A detailed analysis of the genes related to xenobiotic metabolism revealed the potential roles of Variovorax sp. PAMC26660 and other species in bioremediation.

Results

TYGS analysis, dDDH, phylogenetic placement and average nucleotide identity (ANI) analysis identified the strain as Variovorax sp. Cell morphology was assessed using scanning electron microscopy (SEM). On analysis of the core, accessory, and unique genes, xenobiotic metabolism accounted only for the accessory and unique genes. On detailed analysis of the aromatic compound catabolic genes, peripheral pathway related to 4-hydroxybenzoate (4-HB) degradation was found among all species while phenylacetate and tyrosine degradation pathways were present in most of the species including PAMC26660. Likewise, central catabolic pathways, like protocatechuate, gentisate, homogentisate, and phenylacetyl-CoA, were also present. The peripheral pathway for 4-HB degradation was functionally tested using PAMC26660, which resulted in the growth using it as a sole source of carbon.

Conclusions

Computational tools for genome and pan genome analysis are important to understand the behavior of an organism. Xenobiotic metabolism-related genes, that only account for the accessory and unique genes infer evolution through events like lateral gene transfer, mutation and gene rearrangement. 4-HB, an aromatic compound present among lichen species is utilized by lichen-associated Variovorax sp. PAMC26660 as the sole source of carbon. The strain holds genes and pathways for its utilization. Overall, this study outlines the importance of Variovorax in bioremediation and presents the genomic information of the species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08589-3.

Contrasting gene-level signatures of selection with reproductive fitness

Selection leaves signatures in the DNA sequence of genes, with many test statistics devised to detect its action. While these statistics are frequently used to support hypotheses about the adaptive significance of particular genes, the effect these genes have on reproductive fitness is rarely quantified experimentally. Consequently, it is unclear how gene-level signatures of selection are associated with empirical estimates of gene effect on fitness. Eukaryotic datasets that permit this comparison are very limited. Using the model plant Arabidopsis thaliana, for which these resources are available, we calculated seven gene-level substitution and polymorphism-based statistics commonly used to infer selection (dN/dS, NI, DOS, Tajima's D, Fu and Li's D*, Fay and Wu's H, and Zeng's E) and, using knockout lines, compared these to gene-level estimates of effect on fitness. We found that consistent with expectations, essential genes were more likely to be classified as negatively selected. By contrast, using 379 Arabidopsis genes for which data was available, we found no evidence that genes predicted to be positively selected had a significantly different effect on fitness than genes evolving more neutrally. We discuss these results in the context of the analytic challenges posed by Arabidopsis, one of the only systems in which this study could be conducted, and advocate for examination in additional systems. These results are relevant to the evaluation of genome-wide studies across species where experimental fitness data is unavailable, as well as highlighting an increasing need for the latter.

Heterogeneity in conservation of multifunctional partner enzymes with meiotic importance, CDK2 kinase and BRCA1 ubiquitin ligase

The evolution of proteins can be accompanied by changes not only to their amino acid sequences, but also their structural and spatial molecular organization. Comparison of the protein conservation within different taxonomic groups (multifunctional, or highly specific) allows to clarify their specificity and the direction of evolution. Two multifunctional enzymes, cyclin-dependent kinase 2 (CDK2) and BRCA1 ubiquitin ligase, that are partners in some mitotic and meiotic processes were investigated in the present work. Two research methods, bioinformatics and immunocytochemical, were combined to examine the conservation levels of the two enzymes. It has been established that CDK2 is a highly conserved protein in different taxonomic lineages of the eukaryotic tree. Immunocytochemically, a conserved CDK2 pattern was revealed in the meiotic autosomes of five rodent species and partially in domestic turkey and clawed frog. Nevertheless, variable CDK2 distribution was detected at the unsynapsed segments of the rodent X chromosomes. BRCA1 was shown to be highly conserved only within certain mammalian taxa. It was also noted that in those rodent nuclei, where BRCA1 specifically binds to antigens, asynaptic regions of sex chromosomes were positive. BRCA1 staining was not always accompanied by specific binding, and a high nonspecificity in the nucleoplasm was observed. Thus, the studies revealed different conservation of the two enzymes at the level of protein structure as well as at the level of chromosome behavior. This suggests variable rates of evolution due to both size and configuration of the protein molecules and their multifunctionality.

The Missing Expression Level-Evolutionary Rate Anticorrelation in Viruses Does Not Support Protein Function as a Main Constraint on Sequence Evolution

One of the central goals in molecular evolutionary biology is to determine the sources of variation in the rate of sequence evolution among proteins. Gene expression level is widely accepted as the primary determinant of protein evolutionary rate, because it scales with the extent of selective constraints imposed on a protein, leading to the well-known negative correlation between expression level and protein evolutionary rate (the E-R anticorrelation). Selective constraints have been hypothesized to entail the maintenance of protein function, the avoidance of cytotoxicity caused by protein misfolding or nonspecific protein-protein interactions, or both. However, empirical tests evaluating the relative importance of these hypotheses remain scarce, likely due to the nontrivial difficulties in distinguishing the effect of a deleterious mutation on a protein's function versus its cytotoxicity. We realized that examining the sequence evolution of viral proteins could overcome this hurdle. It is because purifying selection against mutations in a viral protein that result in cytotoxicity per se is likely relaxed, whereas purifying selection against mutations that impair viral protein function persists. Multiple analyses of SARS-CoV-2 and nine other virus species revealed a complete absence of any E-R anticorrelation. As a control, the E-R anticorrelation does exist in human endogenous retroviruses where purifying selection against cytotoxicity is present. Taken together, these observations do not support the maintenance of protein function as the main constraint on protein sequence evolution in cellular organisms.

2020 William Allan Award address: genetics as a way of thinking-cultural inheritance from our teachers

This article is based on the address given by the author at the 2020 virtual meeting of The American Society of Human Genetics (ASHG) on October 26, 2020. The video of the original address can be found at the ASHG website.

Nucleotide substitution rates of diatom plastid encoded protein genes are positively correlated with genome architecture

Diatoms are the largest group of heterokont algae with more than 100,000 species. As one of the single-celled photosynthetic organisms that inhabit marine, aquatic and terrestrial ecosystems, diatoms contribute ~ 45% of global primary production. Despite their ubiquity and environmental significance, very few diatom plastid genomes (plastomes) have been sequenced and studied. This study explored patterns of nucleotide substitution rates of diatom plastids across the entire suite of plastome protein-coding genes for 40 taxa representing the major clades. The highest substitution rate was lineage-specific within the araphid 2 taxon Astrosyne radiata and radial 2 taxon Proboscia sp. Rate heterogeneity was also evident in different functional classes and individual genes. Similar to land plants, proteins genes involved in photosynthetic metabolism have lower synonymous and nonsynonymous substitutions rates than those involved in transcription and translation. Significant positive correlations were identified between substitution rates and measures of genomic rearrangements, including indels and inversions, which is a similar result to what was found in legume plants. This work advances the understanding of the molecular evolution of diatom plastomes and provides a foundation for future studies.

Aligning functional network constraint to evolutionary outcomes

BACKGROUND

Functional constraint through genomic architecture is suggested to be an important dimension of genome evolution, but quantitative evidence for this idea is rare. In this contribution, existing evidence and discussions on genomic architecture as constraint for convergent evolution, rapid adaptation, and genic adaptation are summarized into alternative, testable hypotheses. Network architecture statistics from protein-protein interaction networks are then used to calculate differences in evolutionary outcomes on the example of genomic evolution in yeast, and the results are used to evaluate statistical support for these longstanding hypotheses.

RESULTS

A discriminant function analysis lent statistical support to classifying the yeast interactome into hub, intermediate and peripheral nodes based on network neighborhood connectivity, betweenness centrality, and average shortest path length. Quantitative support for the existence of genomic architecture as a mechanistic basis for evolutionary constraint is then revealed through utilizing these statistical parameters of the protein-protein interaction network in combination with estimators of protein evolution.

CONCLUSIONS

As functional genetic networks are becoming increasingly available, it will now be possible to evaluate functional genetic network constraint against variables describing complex phenotypes and environments, for better understanding of commonly observed deterministic patterns of evolution in non-model organisms. The hypothesis framework and methodological approach outlined herein may help to quantify the extrinsic versus intrinsic dimensions of evolutionary constraint, and result in a better understanding of how fast, effectively, or deterministically organisms adapt.

Collateral fitness effects of mutations

The distribution of fitness effects of mutation plays a central role in constraining protein evolution. The underlying mechanisms by which mutations lead to fitness effects are typically attributed to changes in protein specific activity or abundance. Here, we reveal the importance of a mutation's collateral fitness effects, which we define as effects that do not derive from changes in the protein's ability to perform its physiological function. We comprehensively measured the collateral fitness effects of missense mutations in the Escherichia coli TEM-1 β-lactamase antibiotic resistance gene using growth competition experiments in the absence of antibiotic. At least 42% of missense mutations in TEM-1 were deleterious, indicating that for some proteins collateral fitness effects occur as frequently as effects on protein activity and abundance. Deleterious mutations caused improper posttranslational processing, incorrect disulfide-bond formation, protein aggregation, changes in gene expression, and pleiotropic effects on cell phenotype. Deleterious collateral fitness effects occurred more frequently in TEM-1 than deleterious effects on antibiotic resistance in environments with low concentrations of the antibiotic. The surprising prevalence of deleterious collateral fitness effects suggests they may play a role in constraining protein evolution, particularly for highly expressed proteins, for proteins under intermittent selection for their physiological function, and for proteins whose contribution to fitness is buffered against deleterious effects on protein activity and protein abundance.

Comprehensive study of the genes involved in chlorophyll synthesis and degradation pathways in some monocot and dicot plant species

Chlorophyll (Chl) biosynthesis is one of the most important cellular processes essential for plant photosynthesis. Chl degradation pathway is also important catabolic process occurs during leaf senescence, fruit ripening and under biotic or abiotic stress conditions. Here we have systematically investigated the molecular evolution, gene structure, compositional analysis along with ENc plot, correspondence analysis and codon usage bias of the proteins and encoded genes involved in Chl metabolism from monocots and dicots. The gene and species specific phylogenetic trees using amino acid sequences showed clear clustering formation of the selected species based on monocots and dicots but not supported by 18S rRNA. Nucleotide composition of the encoding genes showed that average GC%, GC1%, GC2% and GC3% were higher in monocots. RSCU analysis depicts that genes from monocots for both pathways and genes for synthesis pathway from dicots only biased to G/C-ending synonymous codons but in degradation pathway most optimal codons (except UUG) in dicots biased to A/U-ending synonymous codons. We found strong evidence of episodic diversifying selection at several amino acid sites in all genes investigated. Conserved domain and gene structures were observed for the genes with varying lengths of introns and exons, involved in Chl metabolism along with some intronless genes within synthesis pathway. ENc and correspondence analyses suggested the mutational or selection constraint on the genes to shape the codon usage. These comprehensive studies may be helpful in further research in molecular phylogenetics and genomics and to better understand the evolutionary dynamics of Chl metabolic pathway.Communicated by Ramaswamy H. Sarma.

Epistatic contributions promote the unification of incompatible models of neutral molecular evolution

Mathematical models of evolution help us understand mechanisms driving protein-sequence change. Previous models recapitulate a disjoint subset of statistical features of natural sequences. We present a neutral evolution model that unifies features including extreme variance of the molecular clock’s tick rate and the observation of an evolutionary Stokes shift, an irreversible effect of mutations in the fitness landscape during sequence evolution. We show that interactions between amino acid sites, which inform our fitness metric, are required to observe these features. These interactions are inferred by using direct coupling analysis, which has been successfully utilized to predict protein structures, dynamics, and complexes from coevolutionary information. We anticipate our model will have applications in phylogenetics, ancestral reconstruction of sequences, and protein design.

We introduce a model of amino acid sequence evolution that accounts for the statistical behavior of real sequences induced by epistatic interactions. We base the model dynamics on parameters derived from multiple sequence alignments analyzed by using direct coupling analysis methodology. Known statistical properties such as overdispersion, heterotachy, and gamma-distributed rate-across-sites are shown to be emergent properties of this model while being consistent with neutral evolution theory, thereby unifying observations from previously disjointed evolutionary models of sequences. The relationship between site restriction and heterotachy is characterized by tracking the effective alphabet dynamics of sites. We also observe an evolutionary Stokes shift in the fitness of sequences that have undergone evolution under our simulation. By analyzing the structural information of some proteins, we corroborate that the strongest Stokes shifts derive from sites that physically interact in networks near biochemically important regions. Perspectives on the implementation of our model in the context of the molecular clock are discussed.

The Phenotype Paradox: Lessons From Natural Transcriptome Evolution on How to Engineer Plants

Plants have evolved genome complexity through iterative rounds of single gene and whole genome duplication. This has led to substantial expansion in transcription factor numbers following preferential retention and subsequent functional divergence of these regulatory genes. Here we review how this simple evolutionary network rewiring process, regulatory gene duplication followed by functional divergence, can be used to inspire synthetic biology approaches that seek to develop novel phenotypic variation for future trait based breeding programs in plants.

The systematics of the Cervidae: a total evidence approach

Systematic relationships of cervids have been controversial for decades. Despite new input from molecular systematics, consensus could only be partially reached. The initial, gross (sub) classification based on morphology and comparative anatomy was mostly supported by molecular data. The rich fossil record of cervids has never been extensively tested in phylogenetic frameworks concerning potential systematic relationships of fossil cervids to extant cervids. The aim of this work was to investigate the systematic relationships of extant and fossil cervids using molecular and morphological characters and make implications about their evolutionary history based on the phylogenetic reconstructions. To achieve these objectives, molecular data were compiled consisting of five nuclear markers and the complete mitochondrial genome of 50 extant and one fossil cervids. Several analyses using different data partitions, taxon sampling, partitioning schemes, and optimality criteria were undertaken. In addition, the most extensive morphological character matrix for such a broad cervid taxon sampling was compiled including 168 cranial and dental characters of 41 extant and 29 fossil cervids. The morphological and molecular data were analysed in a combined approach and other comprehensive phylogenetic reconstructions. The results showed that most Miocene cervids were more closely related to each other than to any other cervids. They were often positioned between the outgroup and all other cervids or as the sister taxon to Muntiacini. Two Miocene cervids were frequently placed within Muntiacini. Plio- and Pleistocene cervids could often be affiliated to Cervini, Odocoileini or Capreolini. The phylogenetic analyses provide new insights into the evolutionary history of cervids. Several fossil cervids could be successfully related to living representatives, confirming previously assumed affiliations based on comparative morphology and introducing new hypotheses. New systematic relationships were observed, some uncertainties persisted and resolving systematics within certain taxa remained challenging.

The evolutionary plasticity of chromosome metabolism allows adaptation to constitutive DNA replication stress

Many biological features are conserved and thus considered to be resistant to evolutionary change. While rapid genetic adaptation following the removal of conserved genes has been observed, we often lack a mechanistic understanding of how adaptation happens. We used the budding yeast, Saccharomyces cerevisiae, to investigate the evolutionary plasticity of chromosome metabolism, a network of evolutionary conserved modules. We experimentally evolved cells constitutively experiencing DNA replication stress caused by the absence of Ctf4, a protein that coordinates the enzymatic activities at replication forks. Parallel populations adapted to replication stress, over 1000 generations, by acquiring multiple, concerted mutations. These mutations altered conserved features of two chromosome metabolism modules, DNA replication and sister chromatid cohesion, and inactivated a third, the DNA damage checkpoint. The selected mutations define a functionally reproducible evolutionary trajectory. We suggest that the evolutionary plasticity of chromosome metabolism has implications for genome evolution in natural populations and cancer.

White Spot Syndrome Virus Establishes a Novel IE1/JNK/c-Jun Positive Feedback Loop to Drive Replication

Viruses need to hijack and manipulate host proteins to guarantee their replication. Herein, we uncovered that the DNA virus white spot syndrome virus (WSSV) established a novel positive feedback loop by hijacking the host JNK pathway via its immediate-early 1 (IE1) protein to drive replication. Specifically, the WSSV IE1 bound to host JNK, and enhanced JNK autoactivation by autophosphorylation, and in turn, elevated JNK kinase activity to its substrate c-Jun and induced IE1, which resulted in a viral gene-mediated positive feedback loop. Moreover, the activation of this loop is able to induce wsv056, wsv249, and wsv403, in addition to IE1 itself. Disruption of this loop during WSSV infection by knockdown of JNK, c-Jun or IE1 led to an increased survival rate and lower viral burdens in shrimp. Taken together, this loop may provide a potential target for the development of specific antiviral strategies or agents against WSSV infection.

•

Lvc-Jun promotes WSSV IE1 induction via interacting with the promoter of IE1 gene

•

The interaction of IE1-LvJNK enhances the autophosphorylation of LvJNK

•

IE1 hijacks the JNK/c-Jun cascade to create a feedback loop to drive replication

•

wsv056, wsv249, and wsv403 are also benefit from this positive feedback loop

How to resurrect ancestral proteins as proxies for ancient biogeochemistry

Throughout the history of life, enzymes have served as the primary molecular mediators of biogeochemical cycles by catalyzing the metabolic pathways that interact with geochemical substrates. The byproducts of enzymatic activities have been preserved as chemical and isotopic signatures in the geologic record. However, interpretations of these signatures are limited by the assumption that such enzymes have remained functionally conserved over billions of years of molecular evolution. By reconstructing ancient genetic sequences in conjunction with laboratory enzyme resurrection, preserved biogeochemical signatures can instead be related to experimentally constrained, ancestral enzymatic properties. We may thereby investigate instances within molecular evolutionary trajectories potentially tied to significant biogeochemical transitions evidenced in the geologic record. Here, we survey recent enzyme resurrection studies to provide a reasoned assessment of areas of success and common pitfalls relevant to ancient biogeochemical applications. We conclude by considering the Great Oxidation Event, which provides a constructive example of a significant biogeochemical transition that warrants investigation with ancestral enzyme resurrection. This event also serves to highlight the pitfalls of facile interpretation of paleophenotype models and data, as applied to two examples of enzymes that likely both influenced and were influenced by the rise of atmospheric oxygen - RuBisCO and nitrogenase.

Identification and analysis of the human sex-biased genes

Tremendous differences between human sexes are universally observed. Therefore, identifying and analyzing the sex-biased genes are becoming basically important for uncovering the mystery of sex differences and personalized medicine. Here, we presented a computational method to identify sex-biased genes from public gene expression databases. We obtained 1407 female-biased genes (FGs) and 1096 male-biased genes (MGs) across 14 different tissues. Bioinformatics analysis revealed that compared with MGs, FGs have higher evolutionary rate, higher single-nucleotide polymorphism density, less homologous gene numbers and smaller phyletic age. FGs have lower expression level, higher tissue specificity and later expressed stage in body development. Moreover, FGs are highly involved in immune-related functions, whereas MGs are more enriched in metabolic process. In addition, cellular network analysis revealed that MGs have higher degree, more cellular activating signaling and tend to be located in cellular inner space, whereas FGs have lower degree, more cellular repressing signaling and tend to be located in cellular outer space. Finally, the identified sex-biased genes and the discovered biological insights together can be a valuable resource helpful for investigating sex-biased physiology and medicine, for example sex-biased disease diagnosis and therapy, which represents one important aspect of personalized and precision medicine.

Metabolic Determinants of Enzyme Evolution in a Genome-Scale Bacterial Metabolic Network

Different genes and proteins evolve at very different rates. To identify the factors that explain these differences is an important aspect of research in molecular evolution. One such factor is the role a protein plays in a large molecular network. Here, we analyze the evolutionary rates of enzyme-coding genes in the genome-scale metabolic network of Escherichia coli to find the evolutionary constraints imposed by the structure and function of this complex metabolic system. Central and highly connected enzymes appear to evolve more slowly than less connected enzymes, but we find that they do so as a by-product of their high abundance, and not because of their position in the metabolic network. In contrast, enzymes catalyzing reactions with high metabolic flux-high substrate to product conversion rates-evolve slowly even after we account for their abundance. Moreover, enzymes catalyzing reactions that are difficult to by-pass through alternative pathways, such that they are essential in many different genetic backgrounds, also evolve more slowly. Our analyses show that an enzyme's role in the function of a metabolic network affects its evolution more than its place in the network's structure. They highlight the value of a system-level perspective for studies of molecular evolution.

[Searching for essential genes in cancer genomes]

The concept of essential genes, whose loss of functionality leads to cell death, is one of the fundamental concepts of genetics and is important for fundamental and applied research. This field is particularly promising in relation to oncology, since the search for genetic vulnerabilities of cancer cells allows us to identify new potential targets for antitumor therapy. The modern biotechnology capacities allow carrying out large-scale projects for sequencing somatic mutations in tumors, as well as directly interfering the genetic apparatus of cancer cells. They provided accumulation of a considerable body of knowledge about genetic variants and corresponding phenotypic manifestations in tumors. In the near future this knowledge will find application in clinical practice. This review describes the main experimental and computational approaches to the search for essential genes, concentrating on the application of these methods in the field of molecular oncology.

Position Matters: Network Centrality Considerably Impacts Rates of Protein Evolution in the Human Protein-Protein Interaction Network

The proteins of any organism evolve at disparate rates. A long list of factors affecting rates of protein evolution have been identified. However, the relative importance of each factor in determining rates of protein evolution remains unresolved. The prevailing view is that evolutionary rates are dominantly determined by gene expression, and that other factors such as network centrality have only a marginal effect, if any. However, this view is largely based on analyses in yeasts, and accurately measuring the importance of the determinants of rates of protein evolution is complicated by the fact that the different factors are often correlated with each other, and by the relatively poor quality of available functional genomics data sets. Here, we use correlation, partial correlation and principal component regression analyses to measure the contributions of several factors to the variability of the rates of evolution of human proteins. For this purpose, we analyzed the entire human protein–protein interaction data set and the human signal transduction network—a network data set of exceptionally high quality, obtained by manual curation, which is expected to be virtually free from false positives. In contrast with the prevailing view, we observe that network centrality (measured as the number of physical and nonphysical interactions, betweenness, and closeness) has a considerable impact on rates of protein evolution. Surprisingly, the impact of centrality on rates of protein evolution seems to be comparable, or even superior according to some analyses, to that of gene expression. Our observations seem to be independent of potentially confounding factors and from the limitations (biases and errors) of interactomic data sets.

Correlates of evolutionary rates in the murine sperm proteome

Background

Protein-coding genes expressed in sperm evolve at different rates. To gain deeper insight into the factors underlying this heterogeneity we examined the relative importance of a diverse set of previously described rate correlates in determining the evolution of murine sperm proteins.

Results

Using partial rank correlations we detected several major rate indicators: Phyletic gene age, numbers of protein-protein interactions, and survival essentiality emerged as particularly important rate correlates in murine sperm proteins. Tissue specificity, numbers of paralogs, and untranslated region lengths also correlate significantly with sperm genes’ evolutionary rates, albeit to a lesser extent. Multifunctionality, coding sequence or average intron lengths, and mean expression level have insignificant or virtually no independent effects on evolutionary rates in murine sperm genes. Gene ontology enrichment analyses of three equally sized murine sperm protein groups classified based on their evolutionary rates indicate strongest sperm-specific functional specialization in the most quickly evolving gene class.

Conclusions

We propose a model according to which slowly evolving murine sperm proteins tend to be constrained by factors such as survival essentiality, network connectivity, and/or broad expression. In contrast, evolutionary change may arise especially in less constrained sperm proteins, which might, moreover, be prone to specialize to reproduction-related functions. Our results should be taken into account in future studies on rate variations of reproductive genes.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1157-6) contains supplementary material, which is available to authorized users.

Bayesian Phylogenetic Estimation of Clade Ages Supports Trans-Atlantic Dispersal of Cichlid Fishes

Divergence-time estimation based on molecular phylogenies and the fossil record has provided insights into fundamental questions of evolutionary biology. In Bayesian node dating, phylogenies are commonly time calibrated through the specification of calibration densities on nodes representing clades with known fossil occurrences. Unfortunately, the optimal shape of these calibration densities is usually unknown and they are therefore often chosen arbitrarily, which directly impacts the reliability of the resulting age estimates. As possible solutions to this problem, two nonexclusive alternative approaches have recently been developed, the “fossilized birth–death” (FBD) model and “total-evidence dating.” While these approaches have been shown to perform well under certain conditions, they require including all (or a random subset) of the fossils of each clade in the analysis, rather than just relying on the oldest fossils of clades. In addition, both approaches assume that fossil records of different clades in the phylogeny are all the product of the same underlying fossil sampling rate, even though this rate has been shown to differ strongly between higher level taxa. We here develop a flexible new approach to Bayesian age estimation that combines advantages of node dating and the FBD model. In our new approach, calibration densities are defined on the basis of first fossil occurrences and sampling rate estimates that can be specified separately for all clades. We verify our approach with a large number of simulated data sets, and compare its performance to that of the FBD model. We find that our approach produces reliable age estimates that are robust to model violation, on par with the FBD model. By applying our approach to a large data set including sequence data from over 1000 species of teleost fishes as well as 147 carefully selected fossil constraints, we recover a timeline of teleost diversification that is incompatible with previously assumed vicariant divergences of freshwater fishes. Our results instead provide strong evidence for transoceanic dispersal of cichlids and other groups of teleost fishes.

Genome-wide identification of genes essential for podocyte cytoskeletons based on single-cell RNA sequencing

Gene expression differs substantially among individual cells of the same type. We speculate that genes that are expressed in all but a portion of cells of a given cell type would be likely essential and required for either the cell survival (housekeeping) or for the cell type's unique structure and function, enabling the organism to survive. Here, we performed RNA-seq of 20 mouse podocytes using the Fluidigm C1 system and identified 335 genes that were expressed in all of them. Among them, 239 genes were also expressed in mesangial and endothelial cells and were involved in energy metabolism, protein synthesis, etc., as housekeeping genes. In contrast, 92 genes were preferentially expressed in podocytes (over five-fold versus expression in mesangial and endothelial cells) and are, therefore, the essential candidate genes specific for podocytes. Assessments by bioinformatics, conserved expression in human podocytes, and association with injury/disease all support the essentiality of these genes for podocytes. Factually, 27 of the 92 genes are already known to be essential for podocyte structure and function. Thirty-seven novel genes were functionally analyzed by siRNA silencing, and we found that a deficiency of 30 genes led to either cytoskeletal injury (FGFR1, AOX1, AIF1L, HAUS8, RAB3B, LPIN2, GOLIM4, CERS6, ARHGEF18, ARPC1A, SRGAP1, ITGB5, ILDR2, MPP5, TSC22D1, DNAJC11, SEPT10, MOCS2, FNBP1L, and TMOD3) or significant downregulation of CD2AP and synaptopodin (IFT80, MYOM2, ANXA4, CYB5R4, GPC1, ZNF277, NSF, ITGAV, CRYAB, and MTSS1). Thus, the list of genes essential for podocyte cytoskeletons is expanded by single-cell RNA sequencing. It appears that podocyte-specific essential genes are mainly associated with podocyte cytoskeletons.

The ancestral levels of transcription and the evolution of sexual phenotypes in filamentous fungi

Changes in gene expression have been hypothesized to play an important role in the evolution of divergent morphologies. To test this hypothesis in a model system, we examined differences in fruiting body morphology of five filamentous fungi in the Sordariomycetes, culturing them in a common garden environment and profiling genome-wide gene expression at five developmental stages. We reconstructed ancestral gene expression phenotypes, identifying genes with the largest evolved increases in gene expression across development. Conducting knockouts and performing phenotypic analysis in two divergent species typically demonstrated altered fruiting body development in the species that had evolved increased expression. Our evolutionary approach to finding relevant genes proved far more efficient than other gene deletion studies targeting whole genomes or gene families. Combining gene expression measurements with knockout phenotypes facilitated the refinement of Bayesian networks of the genes underlying fruiting body development, regulation of which is one of the least understood processes of multicellular development.

MOLECULAR GENETIC-DISTANCE ESTIMATES AMONG THE URSIDAE AS INDICATED BY ONE- AND TWO-DIMENSIONAL PROTEIN ELECTROPHORESIS

Evolutionary relationships among eight species of Ursidae (including the giant panda) relative to two Procyonidae species (raccoon and red panda) were estimated based on the extent of electrophoretic variation of 289 radiolabelled fibroblast proteins resolved by two-dimensional gel electrophoresis and among 44 isozyme loci resolved by one-dimensional electrophoresis. Allelic differences among these species were converted to genetic distances, and phenetic trees were constructed. In addition, the electrophoretic data were coded as unit characters, and minimum-length trees were derived based on the Wagner method using maximum parsimony. Regardless of the tree-building method employed, the data sets agreed on the following branching sequence: between 22.4 and 32.3 million years (MY) ago, the ancestors of the procyonids and the ursids split into two lineages. Within 10 MY, the red panda split from the line that led to the raccoon. An ancestor of the giant panda split from the ursid line 18-22 MY ago, and the South American spectacled bear split from the line leading to ursine bears 10.5-15.0 MY B.P. A group of six closely related ursine bears (brown bear, polar bear, Asiatic black bear, Malayan sun bear, American black bear, and sloth bear) diverged from a common ancestor during the past 4-8 MY. Much of this ursine radiation was not resolved by our results, with the exception of a recent (2-3 MY B.P.) divergence of brown bear and polar bear. The topological concordance of the data sets from one- and two-dimensional electrophoresis supports the usefulness of these procedures for evolutionary inference and provides additional precision to the reconstruction of divergence nodes of this carnivore group.

EVOLUTION OF ATLANTIC AND PACIFIC COD: LOSS OF GENETIC VARIATION AND GENE EXPRESSION IN PACIFIC COD

An allozyme investigation of 41 protein-coding loci in two morphologically similar fishes, Atlantic and Pacific cod, indicates that Pacific cod experienced a severe population bottleneck that led to the loss of gene diversity and gene expression. Pacific cod possesses a significantly lesser amount of gene diversity (H = 0.032) than Atlantic cod (H = 0.125) and lacks gene expression for Me-3. Allele-frequency distributions differ between species as predicted by neutral theory: Atlantic cod has a U-shaped distribution, which is expected for populations in drift-mutation equilibrium, whereas Pacific cod has a J-shaped distribution with an excess of low-frequency alleles. This excess may be explained by the appearance of new alleles through mutation which have not yet reached intermediate frequencies through drift. The population bottleneck in Pacific cod was most likely associated with founder populations that dispersed into the Pacific Ocean after the Bering Strait opened. Under the molecular-clock hypothesis a Nei genetic distance of 0.415 (based on 41 loci) suggests that Pacific cod dispersed into the Pacific Ocean soon after the Bering Strait opened in the mid-Pliocene, 3.0 to 3.5 million years ago.

EVOLUTIONARY ASSOCIATIONS OF BROOD PARASITIC FINCHES (VIDUA) AND THEIR HOST SPECIES: ANALYSES OF MITOCHONDRIAL DNA RESTRICTION SITES

The species-specific associations of the African brood parasitic finches Vidua with their estrildid finch host species may have originated by cospeciation with the host species or by later colonizations of new hosts. Predictions of these alternative models were tested in two species groups of brood parasites (indigobirds, paradise whydahs) and their hosts. Phylogenetic analyses suggested that the brood parasites and their hosts did not speciate in parallel. The parasitic indigobirds share mitochondrial haplotypes with each other, and species limits in both indigobirds and paradise whydahs do not correspond with their gene trees. Different parasite species within a region are more closely related to each other than any is to parasites that are associated with its same host species in other regions of Africa. There is little genetic difference between parasite species D̂_i,j < 0.001 in the indigobirds, D̂_i,j = 0.01 in the whydahs). Genetic distances D̂_i,j between the parasite species are less than the genetic distances between their corresponding host species in all parasite-host comparisons, and average only 7.2% as large in the indigobirds as in their hosts and 42% as large in the paradise whydahs as in their hosts. A phylogenetic model that allows ancestral haplotype polymorphisms to be retained in descendant species was compared to a constraint model of species monophyly requiring all but the one ancestral haplotype to be independently derived within each species. The constraint model increases the length of the indigobird tree by 50% over that of the model of retained ancestral polymorphisms; the difference is statistically significant. Both phylogenetic and distance analyses indicate that the brood parasites have become associated with their host species through host switches and independent colonizations of the hosts, rather than through parallel cospeciation with them. The molecular genetic results are supported by recent discoveries of additional host species that are associated with the indigobirds in the field and by variation in the species-specific song behaviors of the brood parasites.

PHYLOGENETIC ESTIMATION OF PLASMID EXCHANGE IN BACTERIA

The existence of differential horizontal gene transfer may be assessed by comparing the phylogenetic trees derived from two different genes. We use this concept to estimate quantitatively the amount of plasmid exchange that has occurred in a bacterial population. By means of computer simulations we studied the effect of gene transfer on the topological distortion between two phylogenetic trees: one obtained from an euchromosomal gene and another from a plasmid-borne sequence, which may be subjected to horizontal transfer. The basic assumptions of our simulations were (a) that plasmid exchange had occurred recently (after the last population split); and (b) that either the amount of chromosomal horizontal exchange was negligible or that it was only a fraction of the amount of plasmid exchange in which case we will be estimating relative amounts of plasmid transfer. We found that the topological difference between two such trees is a function of the number of plasmid exchange events that have occurred. It can be explained by a logistic model that relates the average distortion index between two trees (d_T ) to the number of transfer events (x). The behavior remains the same under different conditions that were tested (symmetry of the topology, number of taxa in the tree, effect of reconstruction errors, mutation after plasmid transfer). We have also tried our method on empirical data from the literature and estimated the amount of gene transfer that may have occurred among Sym plasmids in agricultural field populations of Rhizobium leguminosarum biovar phaseoli. We found that between 15.77 to 29.98% of all genetic types in these populations have been either the source or the target of a plasmid transfer event. When the comparisons were made among trees derived exclusively from plasmid probes this value dropped to 2.00%. Phylogenetic trees derived from symbiotic and nonsymbiotic sequences were also used to infer the number of gene transfer events among 11 isolates from R. galegae. The estimated number of transfer events of symbiotic sequences was 10.515 (although we do not know out of how many genetic types). We concluded that intraspecific transfer of symbiotic sequences is widespread in these two species of the genus Rhizobium.

EVOLUTION OF β-GLUCURONIDASE REGULATION IN THE GENUS MUS

Despite the central role suggested for regulatory mutations in many evolutionary scenarios, there is relatively little information available about the type and extent of regulatory differences between species, or to what extent differences between species are independent of variation within species. To address this issue we have studied the regulatory system of β-glucuronidase, a gene implicated in a murine androgen-inducible pheromone-signalling system. We examined the changes in β-glucuronidase hormonal regulation which have occurred during the radiation of a group of 12 closely related species of mice by assaying β-glucuronidase activity in six different tissues after treatment with estrogen, and with androgen alone and in combination with either estrogen or growth hormone. We also examined in some detail the extent of variation in regulatory responses within species. We found extensive variation in regulatory phenotypes both within and among the species surveyed, suggesting that many of the species examined are currently polymorphic for various regulatory factors that affect inducibility of β-glucuronidase. The variation we observed reflects changes in the ability of the β-glucuronidase gene to respond to hormonal influences, rather than changes in aspects of the hormonal signalling system exterior to the gene. The marked differences among species in the renal and uterine responses to hormonal induction of β-glucuronidase are not easily related to the phylogeny of the genus Mus. If hormonal induction of the gene for β-glucuronidase is subject to natural selection, it appears to be subject to widely fluctuating selective forces. We review evidence that the apparently disorderly evolution of the hormonal responsiveness of β-glucuronidase does not appear to be a unique property of this regulatory system. In contrast to the evolution of many protein sequences, which are tightly correlated with phylogeny and proceed at a relatively constant rate, some, perhaps many, regulatory phenotypes are in rapid evolutionary flux, providing an extensive range of phenotypes upon which selection can act.