Though with constraints imposed by endosymbiosis,preferential attachment is still a plausible mechanism responsible for the evolution of the chloroplast metabolic network

Chloroplasts evolved as a result of endosymbiosis, during which sophisticated mechanisms evolved to translocate nucleus‐encoded plastid‐targeted enzymes into the chloroplast to form the chloroplast metabolic network. Given the constraints and complexity of endosymbiosis, will preferential attachment still be a plausible mechanism for chloroplast metabolic network evolution? We answer this question by analysing the metabolic network properties of the chloroplast and a cyanobacterium, Synechococcus sp. WH8102 (syw). First, we found that enzymes related to more ancient pathways are more connected, and synthetases have the highest connectivity. Most of the enzymes shared by the two densest cores between the chloroplast and syw are synthetases. Second, the highly conserved functional modules mainly consist of highly connected enzymes. Finally, isozymes and enzymes from endosymbiotic gene transfer (EGT) were distributed mainly in conserved modules and showed higher connectivity than nonisozymes or non‐EGT enzymes. These results suggest that even with severe evolutionary constraints imposed by endosymbiosis, preferential attachment is still a plausible mechanism responsible for the evolution of the chloroplast metabolic network. However, the current analysis may not completely differentiate whether the chloroplast network properties reflect the evolution of the chloroplast network through preferential attachment or has been inherited from its cyanobacterial ancestor. To fully differentiate these two possibilities, further analyses of the metabolic network structure properties of organisms at various intermediate evolutionary stages between cyanobacteria and the chloroplast are needed.     

Trade-Offs Predicted by Metabolic Network Structure Give Rise to Evolutionary Specialization and Phenotypic Diversification

Mitigating trade-offs between different resource-utilization functions is key to an organism’s ecological and evolutionary success. These trade-offs often reflect metabolic constraints with a complex molecular underpinning; therefore, their consequences for evolutionary processes have remained elusive. Here, we investigate how metabolic architecture induces resource-utilization constraints and how these constraints, in turn, elicit evolutionary specialization and diversification. Guided by the metabolic network structure of the bacterium Lactococcus cremoris, we selected two carbon sources (fructose and galactose) with predicted coutilization constraints. By evolving L. cremoris on either fructose, galactose, or a mix of both sugars, we imposed selection favoring divergent metabolic specializations or coutilization of both resources, respectively. Phenotypic characterization revealed the evolution of either fructose or galactose specialists in the single-sugar treatments. In the mixed-sugar regime, we observed adaptive diversification: both specialists coexisted, and no generalist evolved. Divergence from the ancestral phenotype occurred at key pathway junctions in the central carbon metabolism. Fructose specialists evolved mutations in the fbp and pfk genes that appear to balance anabolic and catabolic carbon fluxes. Galactose specialists evolved increased expression of pgmA (the primary metabolic bottleneck of galactose metabolism) and silencing of ptnABCD (the main glucose transporter) and ldh (regulator/enzyme of downstream carbon metabolism). Overall, our study shows how metabolic network architecture and historical contingency serve to predict targets of selection and inform the functional interpretation of evolved mutations. The elucidation of the relationship between molecular constraints and phenotypic trade-offs contributes to an integrative understanding of evolutionary specialization and diversification.     

METABOLIC FLUX IS A DETERMINANT OF THE EVOLUTIONARY RATES OF ENZYME‐ENCODING GENES

Relationships between evolutionary rates and gene properties on a genomic, functional, pathway, or system level are being explored to unravel the principles of the evolutionary process. In particular, functional network properties have been analyzed to recognize the constraints they may impose on the evolutionary fate of genes. Here we took as a case study the core metabolic network in human erythrocytes and we analyzed the relationship between the evolutionary rates of its genes and the metabolic flux distribution throughout it. We found that metabolic flux correlates with the ratio of nonsynonymous to synonymous substitution rates. Genes encoding enzymes that carry high fluxes have been more constrained in their evolution, while purifying selection is more relaxed in genes encoding enzymes carrying low metabolic fluxes. These results demonstrate the importance of considering the dynamical functioning of gene networks when assessing the action of selection on system‐level properties.     

Genomic signatures of local adaptation in the Drosophila immune response

As environments and pathogen landscapes shift, host defenses must evolve to remain effective. Due to this selection pressure, among-species comparisons of genetic sequence data often find immune genes to be among the fastest evolving genes across the genome. The full extent and nature of these immune adaptations, however, remain largely unexplored. In a recent study, we analyzed patterns of selection within distinct components of the Drosophila melanogaster immune pathway. While we found evidence of positive selection within some immune processes, immune genes were not universally characterized by signatures of strong selection. On the contrary, we even found that some immune functions show greater than expected constraint. Overall these results highlight 2 major factors that appear to play an outsize role in determining a gene's evolutionary rate: the type of pathogen the gene targets and the gene's position within the immune network. These results join a growing body of literature that highlight the complexity of immune adaptation. Rather than there being uniformly strong selection across all immune genes, a combination of pathogen-specificity and host genetic constraints appear to play key roles in determining each immune gene's individual evolutionary trajectory.     

A plant toxin mediated mechanism for the lag in snowshoe hare population recovery following cyclic declines

A necessary condition for a snowshoe hare population to cycle is reduced reproduction after the population declines. But the cause of a cyclic snowshoe hare population's reduced reproduction during the low phase of the cycle, when predator density collapses, is not completely understood. We propose that moderate‐severe browsing by snowshoe hares upon preferred winter‐foods could increase the toxicity of some of the hare's best winter‐foods during the following hare low, with the result being a decline in hare nutrition that could reduce hare reproduction. We used a combination of modeling and experiments to explore this hypothesis. Using the shrub birch Betula glandulosa as the plant of interest, the model predicted that browsing by hares during a hare cycle peak, by increasing the toxicity B. glandulosa twigs during the following hare low, could cause a hare population to cycle. The model's assumptions were verified with assays of dammarane triterpenes in segments of B. glandulosa twigs and captive hare feeding experiments conducted in Alaska during February and March 1986. The model's predictions were tested with estimates of hare density and measurements of B. glandulosa twig growth made at Kluane, Yukon from 1988–2008. The empirical tests supported the model's predictions. Thus, we have concluded that a browsing‐caused increase in twig toxicity that occurs during the hare cycle's low phase could reduce hare reproduction during the low phase of the hare cycle.     

Stronger selective constraint on downstream genes in the oxidative phosphorylation pathway of cetaceans

The oxidative phosphorylation (OXPHOS) pathway is an efficient way to produce energy via adenosine triphosphate (ATP), which is critical for sustaining an energy supply for cetaceans in a hypoxic environment. Several studies have shown that natural selection may shape the evolution of the genes involved in OXPHOS. However, how network architecture drives OXPHOS protein sequence evolution remains poorly explored. Here, we investigated the evolutionary patterns of genes in the OXPHOS pathway across six cetacean genomes within the framework of a functional network. Our results show a negative correlation between the strength of purifying selection and pathway position. This result indicates that downstream genes were subjected to stronger evolutionary constraints than upstream genes, which may be due to the dual function of ATP synthase in the OXPHOS pathway. Additionally, there was a positive correlation between codon usage bias and omega (ω = dN/dS) and a negative correlation with synonymous substitution rate (dS), indicating that the stronger selective constraint on genes (with less biased codon usage) along the OXPHOS pathway is attributable to an increase in the rate of synonymous substitution. Surprisingly, there was no significant correlation between protein–protein interactions and the evolutionary estimates, implying that highly connected enzymes may not always show greater evolutionary constraints. Compared with that observed for terrestrial mammals, we found that the signature of positive selection detected in five genes (ATP5J, LHPP, PPA1, UQCRC1 and UQCRQ) was cetacean‐specific, reflecting the importance of OXPHOS for survival in hypoxic, aquatic environments.     

Beyond topology: coevolution of structure and flux in metabolic networks

Interactions between the structure of a metabolic network and its functional properties underlie its evolutionary diversification, but the mechanism by which such interactions arise remains elusive. Particularly unclear is whether metabolic fluxes that determine the concentrations of compounds produced by a metabolic network, are causally linked to a network's structure or emerge independently of it. A direct empirical study of populations where both structural and functional properties vary among individuals’ metabolic networks is required to establish whether changes in structure affect the distribution of metabolic flux. In a population of house finches (Haemorhous mexicanus), we reconstructed full carotenoid metabolic networks for 442 individuals and uncovered 11 structural variants of this network with different compounds and reactions. We examined the consequences of this structural diversity for the concentrations of plumage‐bound carotenoids produced by flux in these networks. We found that concentrations of metabolically derived, but not dietary carotenoids, depended on network structure. Flux was partitioned similarly among compounds in individuals of the same network structure: within each network, compound concentrations were closely correlated. The highest among‐individual variation in flux occurred in networks with the strongest among‐compound correlations, suggesting that changes in the magnitude, but not the distribution of flux, underlie individual differences in compound concentrations on a static network structure. These findings indicate that the distribution of flux in carotenoid metabolism closely follows network structure. Thus, evolutionary diversification and local adaptations in carotenoid metabolism may depend more on the gain or loss of enzymatic reactions than on changes in flux within a network structure.     

Genomewide landscape of gene–metabolome associations in Escherichia coli

Metabolism is one of the best‐understood cellular processes whose network topology of enzymatic reactions is determined by an organism's genome. The influence of genes on metabolite levels, however, remains largely unknown, particularly for the many genes encoding non‐enzymatic proteins. Serendipitously, genomewide association studies explore the relationship between genetic variants and metabolite levels, but a comprehensive interaction network has remained elusive even for the simplest single‐celled organisms. Here, we systematically mapped the association between > 3,800 single‐gene deletions in the bacterium Escherichia coli and relative concentrations of > 7,000 intracellular metabolite ions. Beyond expected metabolic changes in the proximity to abolished enzyme activities, the association map reveals a largely unknown landscape of gene–metabolite interactions that are not represented in metabolic models. Therefore, the map provides a unique resource for assessing the genetic basis of metabolic changes and conversely hypothesizing metabolic consequences of genetic alterations. We illustrate this by predicting metabolism‐related functions of 72 so far not annotated genes and by identifying key genes mediating the cellular response to environmental perturbations.     

Adenosine accumulation causes metabolic disorders in testes and associates with lower testosterone level in obese mice

Overweight and obese men face numerous health problems, including type 2 diabetes, subfertility, and even infertility. However, few studies have focused on the effects of nutritional status and obesity‐related regulatory signals on fertility deficiency. Our previous observations have shown that the elevation of plasma 5'‐adenosine monophosphate (5'‐AMP) and the accumulation of adenosine in liver and muscle of obese diabetic db/db mice are related to insulin resistance. Here, we found that adenosine accumulation in testis is a common marker of both genetic obesity and high‐fat‐diet induced obese mice. An messenger RNA sequencing analysis indicated that 78 upregulated genes and 155 downregulated genes in the testis of 5'‐AMP‐treated mice overlapped with the same genes in the testis of ob/ob mice, and these genes belonged to the clusters of steroid metabolic process and regulation of hormone levels, respectively. Serum testosterone was reduced in ob/ob and 5'‐AMP‐treated mice. Metabolomic analysis based on ¹H nuclear magnetic resonance showed that the testicular metabolic profiles of ob/ob mice were similar to those of 5'‐AMP treated mice. Exogenous 5'‐AMP inhibited the phosphorylation of AKT and mammalian target of rapamycin signal transduction and reduced the proliferating cell nuclear antigen expressions in testes. Our results suggest that the accumulation of adenosine causes metabolic disorders in testes and associates lower testosterone level in obese mice.     

The evolution of queen pheromones in the ant genus Lasius

Queen pheromones are among the most important chemical messages regulating insect societies yet they remain largely undiscovered, hindering research into interesting proximate and ultimate questions. Identifying queen pheromones in multiple species would give new insight into the selective pressures and evolutionary constraints acting on these ubiquitous signals. Here, we present experimental and comparative evidence that 3‐methylalkanes, hydrocarbons present on the queen's cuticle, are a queen pheromone throughout the ant genus Lasius. Interspecific variation in the chemical profile is consistent with 3‐methylalkanes evolving more slowly than other types of hydrocarbons, perhaps due to differential selection or evolutionary constraints. We argue that the sensory ecology of the worker response imposes strong stabilizing selection on queen pheromones relative to other hydrocarbons. 3‐Methylalkanes are also strongly physiologically and genetically coupled with fecundity in at least one Lasius species, which may translate into evolutionary constraints. Our results highlight how honest signalling could minimize evolutionary conflict over reproduction, promoting the evolution and maintenance of eusociality.     

Sex‐biased dispersal,kin selection and the evolution of sexual conflict

There is growing interest in resolving the curious disconnect between the fields of kin selection and sexual selection. Rankin's (2011, J. Evol. Biol. 24 , 71–81) theoretical study of the impact of kin selection on the evolution of sexual conflict in viscous populations has been particularly valuable in stimulating empirical research in this area. An important goal of that study was to understand the impact of sex‐specific rates of dispersal upon the coevolution of male‐harm and female‐resistance behaviours. But the fitness functions derived in Rankin's study do not flow from his model's assumptions and, in particular, are not consistent with sex‐biased dispersal. Here, we develop new fitness functions that do logically flow from the model's assumptions, to determine the impact of sex‐specific patterns of dispersal on the evolution of sexual conflict. Although Rankin's study suggested that increasing male dispersal always promotes the evolution of male harm and that increasing female dispersal always inhibits the evolution of male harm, we find that the opposite can also be true, depending upon parameter values.     

Improved network performance via antagonism: From synthetic rescues to multi‐drug combinations

Recent research shows that a faulty or sub‐optimally operating metabolic network can often be rescued by the targeted removal of enzyme‐coding genes – the exact opposite of what traditional gene therapy would suggest. Predictions go as far as to assert that certain gene knockouts can restore the growth of otherwise nonviable gene‐deficient cells. Many questions follow from this discovery: What are the underlying mechanisms? How generalizable is this effect? What are the potential applications? Here, I approach these questions from the perspective of compensatory perturbations on networks. Relations are drawn between such synthetic rescues and naturally occurring cascades of reaction inactivation, as well as their analogs in physical and other biological networks. I specially discuss how rescue interactions can lead to the rational design of antagonistic drug combinations that select against resistance and how they can illuminate medical research on cancer, antibiotics, and metabolic diseases. Editor's suggested further reading in BioEssays The evolutionary context of robust and redundant cell biological mechanisms Abstract Reprogramming cell fates: reconciling rarity with robustness Abstract     

Evolutionary pattern of the regulatory network for flower development: Insights gained from a comparison of two Arabidopsis species

Abstract Previous studies on Arabidopsis thaliana and other model plants have indicated that the development of a flower is controlled by a regulatory network composed of genes and the interactions among them. Studies on the evolution of this network will therefore help understand the genetic basis that underlies flower evolution. In this study, by reviewing the most recent published work, we added 31 genes into the previously proposed regulatory network for flower development. Thus, the number of genes reached 60. We then compared the composition, structure, and evolutionary rate of these genes between A. thaliana and one of its allies, A. lyrata. We found that two genes (FLC and MAF2) show 1: 2 and 2: 2 relationships between the two species, suggesting that they have experienced independent, post‐speciation duplications. Of the remaining 58 genes, 35 (60.3%) have diverged in exon–intron structure and, consequently, code for proteins with different sequence features and functions. Molecular evolutionary analyses further revealed that, although most floral genes have evolved under strong purifying selection, some have evolved under relaxed or changed constraints, as evidenced by the elevation of nonsynonymous substitution rates and/or the presence of positively selected sites. Taken together, these results suggest that the regulatory network for flower development has evolved rather rapidly, with changes in the composition, structure, and functional constraint of genes, as well as the interactions among them, being the most important contributors.     

Signatures of competition and strain structure within the major blood‐stage antigen of Plasmodium falciparum in a local community in Ghana

The concept of niche partitioning has received considerable theoretical attention at the interface of ecology and evolution of infectious diseases. Strain theory postulates that pathogen populations can be structured into distinct nonoverlapping strains by frequency‐dependent selection in response to intraspecific competition for host immune space. The malaria parasite Plasmodium falciparum presents an opportunity to investigate this phenomenon in nature, under conditions of high recombination rate and extensive antigenic diversity. The parasite's major blood‐stage antigen, PfEMP1, is encoded by the hyperdiverse var genes. With a dataset that includes thousands of var DBLα sequence types sampled from asymptomatic cases within an area of high endemicity in Ghana, we address how var diversity is distributed within isolates and compare this to the distribution of microsatellite allelic diversity within isolates to test whether antigenic and neutral regions of the genome are structured differently. With respect to var DBLα sequence types, we find that on average isolates exhibit significantly lower overlap than expected randomly, but that there also exists frequent pairs of isolates that are highly related. Furthermore, the linkage network of var DBLα sequence types reveals a pattern of nonrandom modularity unique to these antigenic genes, and we find that modules of highly linked DBLα types are not explainable by neutral forces related to var recombination constraints, microsatellite diversity, sampling location, host age, or multiplicity of infection. These findings of reduced overlap and modularity among the var antigenic genes are consistent with a role for immune selection as proposed by strain theory. Identifying the evolutionary and ecological dynamics that are responsible for the nonrandom structure in P. falciparum antigenic diversity is important for designing effective intervention in endemic areas.     

DOLLO'S LAW AND THE IRREVERSIBILITY OF DIGIT LOSS IN BACHIA

Several recent studies conclude that exceptions to Dollo's law are more common than used to be thought. If the claims are true this would change our view on the role of developmental constraints in the evolution of body plans. One study claims the reevolution of lost digits in the lizard genus Bachia ( Kohlsdorf and Wagner 2006 ). We evaluate this claim. We conclude that the proposed molecular phylogenetic tree is in conflict with evolutionary mechanisms concerning the biogeography of lizards and with morphology‐based phylogenies. A reanalysis of the molecular data does not support the topology of the published tree. Furthermore, two implicit assumptions, that digit numbers are fixed and that polydactyly evolves independently from other characters, are incorrect. We conclude that there is no convincing support for reevolution of digits in Bachia. We discuss our findings in the light of the current evidence for the reversal of losses of complex traits. We conclude that in metazoans, exceptions to Dollo's law are mainly found among meristic traits that originate relatively late during embryogenesis, when developmental systems are more compartmentalized. Finally, our study shows that phylogenetic analyses should incorporate evolutionary mechanisms including constraints, variation, and selection, not only for correct phylogenetic reconstruction, but also for correct evolutionary inference.     

The miR‐379/miR‐410 cluster at the imprinted Dlk1‐Dio3 domain controls neonatal metabolic adaptation

In mammals, birth entails complex metabolic adjustments essential for neonatal survival. Using a mouse knockout model, we identify crucial biological roles for the miR‐379/miR‐410 cluster within the imprinted Dlk1‐Dio3 region during this metabolic transition. The miR‐379/miR‐410 locus, also named C14MC in humans, is the largest known placental mammal‐specific miRNA cluster, whose 39 miRNA genes are expressed only from the maternal allele. We found that heterozygote pups with a maternal—but not paternal—deletion of the miRNA cluster display partially penetrant neonatal lethality with defects in the maintenance of energy homeostasis. This maladaptive metabolic response is caused, at least in part, by profound changes in the activation of the neonatal hepatic gene expression program, pointing to as yet unidentified regulatory pathways that govern this crucial metabolic transition in the newborn's liver. Not only does our study highlight the physiological importance of miRNA genes that recently evolved in placental mammal lineages but it also unveils additional layers of RNA‐mediated gene regulation at the Dlk1‐Dio3 domain that impose parent‐of‐origin effects on metabolic control at birth and have likely contributed to mammal evolution.     

Why only some plants emit isoprene

Isoprene (2‐methyl‐1,3‐butadiene) is emitted from many plants and it appears to have an adaptive role in protecting leaves from abiotic stress. However, only some species emit isoprene. Isoprene emission has appeared and been lost many times independently during the evolution of plants. As an example, our phylogenetic analysis shows that isoprene emission is likely ancestral within the family Fabaceae (= Leguminosae), but that it has been lost at least 16 times and secondarily gained at least 10 times through independent evolutionary events. Within the division Pteridophyta (ferns), we conservatively estimate that isoprene emissions have been gained five times and lost two times through independent evolutionary events. Within the genus Quercus (oaks), isoprene emissions have been lost from one clade, but replaced by a novel type of light‐dependent monoterpene emissions that uses the same metabolic pathways and substrates as isoprene emissions. This novel type of monoterpene emissions has appeared at least twice independently within Quercus, and has been lost from 9% of the individuals within a single population of Quercus suber. Gain and loss of gene function for isoprene synthase is possible through relatively few mutations. Thus, this trait appears frequently in lineages; but, once it appears, the time available for evolutionary radiation into environments that select for the trait is short relative to the time required for mutations capable of producing a non‐functional isoprene synthase gene. The high frequency of gains and losses of the trait and its heterogeneous taxonomic distribution in plants may be explained by the relatively few mutations necessary to produce or lose the isoprene synthase gene combined with the assumption that isoprene emission is advantageous in a narrow range of environments and phenotypes.     

Basic and applied uses of genome‐scale metabolic network reconstructions of Escherichia coli

The genome‐scale model (GEM) of metabolism in the bacterium Escherichia coli K‐12 has been in development for over a decade and is now in wide use. GEM‐enabled studies of E. coli have been primarily focused on six applications: (1) metabolic engineering, (2) model‐driven discovery, (3) prediction of cellular phenotypes, (4) analysis of biological network properties, (5) studies of evolutionary processes, and (6) models of interspecies interactions. In this review, we provide an overview of these applications along with a critical assessment of their successes and limitations, and a perspective on likely future developments in the field. Taken together, the studies performed over the past decade have established a genome‐scale mechanistic understanding of genotype–phenotype relationships in E. coli metabolism that forms the basis for similar efforts for other microbial species. Future challenges include the expansion of GEMs by integrating additional cellular processes beyond metabolism, the identification of key constraints based on emerging data types, and the development of computational methods able to handle such large‐scale network models with sufficient accuracy.