The speciation of conger eel galectins by rapid adaptive evolution

Many cases of accelerated evolution driven by positive Darwinian selection are identified in the genes of venomous and reproductive proteins. This evolutional phenomenon might have important consequences in their gene-products' functions, such as multiple specific toxins for quick immobilization of the prey and the establishment of barriers to fertilization that might lead to speciation, and in the molecular evolution of novel genes. Recently, we analyzed the molecular evolution of two galectins isolated from the skin mucus of conger eel (Conger myriaster), named congerins I and II, by cDNA cloning and X-ray structural analysis, and we found that they have evolved in the rapid adaptive manner to emergence of a new structure including strand-swapping and a unique new ligand-binding site. In this review article we summarize and discuss the molecular evolution, especially the rapid adaptive evolution, and the structure-function relationships of conger eel galectins.     

Accelerated evolution in the protein-coding region of galectin cDNAs, congerin I and congerin II, from skin mucus of conger eel (Conger myriaster).

Two cDNAs encoding galectins named congerins I and II from the skin mucus of conger eel (Conger myriaster) were isolated and sequenced. Comparison of the nucleotide sequences of congerins I and II showed that the sequence similarities of the 5' and 3' untranslated regions (86 and 88%, respectively) were much higher than those of the protein-coding region (73%). The numbers of nucleotide substitutions per site (KN) for the untranslated regions are smaller than the numbers of nucleotide substitutions per synonymous site (KS) for the protein coding region. Furthermore, nonsynonymous nucleotide substitutions have accelerated more frequently than synonymous nucleotide substitutions in the protein coding region (KA/KS = 2.57). These results suggest that accelerated substitutions have occurred in the protein-coding regions of galectin genes to generate diverse galectins with different molecular properties. Northern blot analysis showed that both congerins were expressed not only in the skin tissues but also in the stomach of conger eel.     

Allosteric Regulation of the Carbohydrate-binding Ability of a Novel Conger Eel Galectin by D-Mannoside

Conger eel has two galectins, termed congerins I and II (Con I and II), that function in mucus as biodefense molecules. Con I and II have acquired a novel protein fold via domain swapping and a new ligand-binding site by accelerated evolution, which enables recognition of some marine bacteria. In this study, we identified a new congerin isotype, congerin P (Con-P), from the peritoneal cells of conger eel. Although Con-P displayed obvious homology with galectins, we observed substitution of 7 out of 8 amino acid residues in the carbohydrate recognition domain that are conserved in all other known galectins. To understand the structure-function relationships of this unique galectin, recombinant Con-P was successfully expressed in Escherichia coli by using a Con II-tagged fusion protein system and subsequently characterized. In the presence of d-mannose, Con-P displayed 30-fold greater hemagglutinating activity than Con I; however, no activity was observed without mannose, indicating that d-mannoside can act as a modulator of Con-P. Frontal affinity chromatography analysis showed that activated Con-P, allosterically induced by mannose, displayed affinity for oligomannose-type sugars as well as N-acetyllactosamine-type β-galactosides. Thus, Con-P represents a new member of the galectin family with unique properties.     

Crystal structure of a conger eel galectin (congerin II) at 1.45A resolution: implication for the accelerated evolution of a new ligand-binding site following gene duplication

The crystal structure of congerin II, a galectin family lectin from conger eel, was determined at 1.45A resolution. The previously determined structure of its isoform, congerin I, had revealed a fold evolution via strand swap; however, the structure of congerin II described here resembles other prototype galectins. A comparison of the two congerin genes with that of several other galectins suggests acceralated evolution of both congerin genes following gene duplication. The presence of a Mes (2-[N-morpholino]ethanesulfonic acid) molecule near the carbohydrate-binding site in the crystal structure points to the possibility of an additional binding site in congerin II. The binding site consists of a group of residues that had been replaced following gene duplication suggesting that the binding site was built under selective pressure. Congerin II may be a protein specialized for biological defense with an affinity for target carbohydrates on parasites' cell surface.     

Reconstruction of a probable ancestral form of conger eel galectins revealed their rapid adaptive evolution process for specific carbohydrate recognition

Recently, many cases of rapid adaptive evolution, which is characterized by the higher substitution rates of nonsynonymous substitutions to synonymous ones, have been identified in the various genes of venomous and biodefense proteins, including the conger eel galectins, congerins I and II (ConI and ConII). To understand the evolutionary process of congerins, we prepared a probable ancestral form, Con-anc, corresponding to the putative amino acid sequence at the divergence of ConI and ConII in phylogenetic tree with 76% and 61% sequence identities to the current proteins, respectively. Con-anc and ConII had comparable thermostability and similar carbohydrate specificities in general, whereas ConI was more thermostable and showed different carbohydrate specificities. Con-anc showed decreased specificity to oligosaccharides with alpha 2,3-sialyl galactose moieties. These suggest that ConI and ConII have evolved via accelerated evolution under significant selective pressure to increase the thermostability and to acquire the activity to bind to alpha2,3-sialyl galactose present in pathogenic bacteria, respectively. Furthermore, comparative mutagenesis analyses of Con-anc and congerins revealed the structural basis for specific recognition of ConII to alpha2,3-sialyl galactose moiety, and strong binding ability of ConI to oligosaccharides including lacto-N-biosyl (Galbeta1-3GlcNAc) or lacto-N-neobiosyl (Galbeta1-4GlcNAc) residues, respectively. Thus, protein engineering using a probable ancestral form presented here is a powerful approach not only to determine the evolutionary process but also to investigate the structure-activity relationships of proteins.     

Quantifying Adaptive Evolution in the Drosophila Immune System

It is estimated that a large proportion of amino acid substitutions in Drosophila have been fixed by natural selection, and as organisms are faced with an ever-changing array of pathogens and parasites to which they must adapt, we have investigated the role of parasite-mediated selection as a likely cause. To quantify the effect, and to identify which genes and pathways are most likely to be involved in the host–parasite arms race, we have re-sequenced population samples of 136 immunity and 287 position-matched non-immunity genes in two species of Drosophila. Using these data, and a new extension of the McDonald-Kreitman approach, we estimate that natural selection fixes advantageous amino acid changes in immunity genes at nearly double the rate of other genes. We find the rate of adaptive evolution in immunity genes is also more variable than other genes, with a small subset of immune genes evolving under intense selection. These genes, which are likely to represent hotspots of host–parasite coevolution, tend to share similar functions or belong to the same pathways, such as the antiviral RNAi pathway and the IMD signalling pathway. These patterns appear to be general features of immune system evolution in both species, as rates of adaptive evolution are correlated between the D. melanogaster and D. simulans lineages. In summary, our data provide quantitative estimates of the elevated rate of adaptive evolution in immune system genes relative to the rest of the genome, and they suggest that adaptation to parasites is an important force driving molecular evolution.     

High-level Expression and Characterization of Fully Active Recombinant Conger Eel Galectins in Eschericia coli

An expression system for recombinant conger eel galectins, congerins I and II, were constructed using the pTV 118N plasmid vector and Escherichia coli. Recombinant congerins I and II could be obtained in the soluble active form with high quantitative yield. Mutation of codons for Val and Leu located in the N-terminal region of Con I increased the expression efficiency. Purification of recombinant proteins were done by only two chromatographical steps from E. coli extract. The purified recombinant congerins were found to be almost the same as the native ones except for the acetyl group at the N-terminus; that is, they showed the same structures and carbohydrate binding activities, suggesting that N-terminal acetyl groups of congerins were not significant for activity.     

Galectin genes: Regulation of expression

In this review we have summarized the more recent studies on the expression of mammalian galectins. One interesting observation that can be made is that in most of microarrays and/or differential display analysis performed in recent years one or more galectins have been picked up. From a critical evaluation of the pertinent studies the main conclusion that can be drawn is that, although it is not yet clear whether the 14 galectins identified so far have functions in common, a striking common feature of all galectins is the strong modulation of their expression during development, differentiation stages and under different physiological or pathological conditions. This suggests that the expression of different galectins is finely tuned and possibly coordinated. In spite of these observations it is rather unexpected that very few studies have been performed on the molecular mechanisms governing the activity of galectin genes. Published in 2004.     

LPS induces gene expression of interleukin-1beta in conger eel (Conger myriaster) macrophages: first cytokine sequence within Anguilliformes

Our aim was to comprehensively analyze the immune response in Anguilliformes macrophages and to survey cytokine genes expressed from them. We therefore used suppression subtractive hybridization (SSH) to randomly clone molecules that are specifically expressed in conger eel (Conger myriaster) macrophages when cells are stimulated by LPS. As a result, we succeeded in identifying a conger eel IL-1beta. This is the first report on cytokines in Anguilliformes, which is the most ancient order in living teleosts.     

Protein engineering of conger eel galectins by tracing of molecular evolution using probable ancestral mutants

Background  

Conger eel galectins, congerin I (ConI) and congerin II (ConII), show the different molecular characteristics resulting from accelerating evolution. We recently reconstructed a probable ancestral form of congerins, Con-anc. It showed properties similar to those of ConII in terms of thermostability and carbohydrate recognition specificity, although it shares a higher sequence similarity with ConI than ConII.     

Molecular Specificity,Convergence and Constraint Shape Adaptive Evolution in Nutrient-Poor Environments

One of the central goals of evolutionary biology is to explain and predict the molecular basis of adaptive evolution. We studied the evolution of genetic networks in Saccharomyces cerevisiae (budding yeast) populations propagated for more than 200 generations in different nitrogen-limiting conditions. We find that rapid adaptive evolution in nitrogen-poor environments is dominated by the de novo generation and selection of copy number variants (CNVs), a large fraction of which contain genes encoding specific nitrogen transporters including PUT4, DUR3 and DAL4. The large fitness increases associated with these alleles limits the genetic heterogeneity of adapting populations even in environments with multiple nitrogen sources. Complete identification of acquired point mutations, in individual lineages and entire populations, identified heterogeneity at the level of genetic loci but common themes at the level of functional modules, including genes controlling phosphatidylinositol-3-phosphate metabolism and vacuole biogenesis. Adaptive strategies shared with other nutrient-limited environments point to selection of genetic variation in the TORC1 and Ras/PKA signaling pathways as a general mechanism underlying improved growth in nutrient-limited environments. Within a single population we observed the repeated independent selection of a multi-locus genotype, comprised of the functionally related genes GAT1, MEP2 and LST4. By studying the fitness of individual alleles, and their combination, as well as the evolutionary history of the evolving population, we find that the order in which these mutations are acquired is constrained by epistasis. The identification of repeatedly selected variation at functionally related loci that interact epistatically suggests that gene network polymorphisms (GNPs) may be a frequent outcome of adaptive evolution. Our results provide insight into the mechanistic basis by which cells adapt to nutrient-limited environments and suggest that knowledge of the selective environment and the regulatory mechanisms important for growth and survival in that environment greatly increase the predictability of adaptive evolution.     

Effects of the piezo-tolerance of cultured deep-sea eel cells on survival rates, cell proliferation, and cytoskeletal structures

We investigated the pressure tolerance of deep-sea eel (Simenchelys parasiticus; habitat depth, 366–2,630 m) cells, conger eel (Conger myriaster) cells, and mouse 3T3-L1 cells. Although there were no living mouse 3T3-L1 and conger eel cells after 130 MPa (0.1 MPa = 1 bar) hydrostatic pressurization for 20 min, all deep-sea eel cells remained alive after being subjected to pressures up to 150 MPa for 20 min. Pressurization at 40 MPa for 20 min induced disruption of actin and tubulin filaments with profound cell-shape changes in the mouse and conger eel cells. In the deep-sea eel cells, microtubules and some actin filaments were disrupted after being subjected to hydrostatic pressure of 100 MPa and greater for 20 min. Conger eel cells were sensitive to pressure and did not grow at 10 MPa. Mouse 3T3-L1 cells grew faster under pressure of 5 MPa than at atmospheric pressure and stopped growing at 18 MPa. Deep-sea eel cells were capable of growth in pressures up to 25 MPa and stopped growing at 30 MPa. Deep-sea eel cells required 4 h at 20 MPa to finish the M phase, which was approximately fourfold the time required under atmospheric conditions.     

High-level expression and characterization of fully active recombinant conger eel galectins in Eschericia coli

An expression system for recombinant conger eel galectins, congerins I and II, were constructed using the pTV 118N plasmid vector and Escherichia coli. Recombinant congerins I and II could be obtained in the soluble active form with high quantitative yield. Mutation of codons for Val and Leu located in the N-terminal region of Con I increased the expression efficiency. Purification of recombinant proteins were done by only two chromatographical steps from E. coli extract. The purified recombinant congerins were found to be almost the same as the native ones except for the acetyl group at the N-terminus; that is, they showed the same structures and carbohydrate binding activities, suggesting that N-terminal acetyl groups of congerins were not significant for activity.     

Perspectives of genomic diversification and molecular recombination towards R-gene evolution in plants

Plants are under strong evolutionary pressure in developing new and noble R genes to recognize pathogen avirulence (avr) determinants and bring about stable defense for generation after generations. Duplication, sequence variation by mutation, disparity in the length and structure of leucine rich repeats etc., causes tremendous variations within and among R genes in a plant thereby developing diverse recognitional specificity suitable enough for defense against new pathogens. Recent studies on genome sequencing, diversity and population genetics in different plants have thrown new insights on the molecular evolution of these genes. Tandem and segmental duplication are important factors in R gene abundance as inferred from the distribution of major nucleotide binding site-leucine rich repeats (NBS-LRRs) type R-genes in plant genomes. Likewise, R-gene evolution is also thought to be facilitated by cluster formation thereby causing recombination and sequence exchange and resulting in haplotypic diversity. Population studies have further proven that balancing selection is responsible for the maintenance of allelic diversity in R genes. In this review, we emphasize and discuss on improved perspectives towards the molecular mechanisms and selection pressure responsible for the evolution of NBS-LRR class resistance genes in plants.     

Inducible galectins are expressed in the inflamed pharynx of the ascidian Ciona intestinalis

Although ascidians belong to a key group in chordate phylogenesis, amino acid sequences of Ciona intestinalis galectin-CRDs (CiLgals-a and -b) have been retained too divergent from vertebrate galectins. In the present paper, to contribute in disclosing Bi-CRD galectin evolution a novel attempt was carried out on CiLgals-a and -b CRDs phylogenetic analysis, and their involvement in ascidian inflammatory responses was shown. CiLgals resulted aligned with Bi-CRD galectins from vertebrates (Xenopus tropicalis, Gallus gallus, Mus musculus, Homo sapiens), cephalochordates (Branchiostoma floridae), echinoderms (Strongylocentrotus purpuratus) and a mono-CRD galectin from the ascidian Clavelina picta. The CiLgals-a N-terminal and C-terminal CRDs contain the signature sequence involved in carbohydrate binding, whereas the CiLgals-b C-CRD presents only three out of seven key aminoacids and it could not be suitable as sugar binding motif. Sequence similarity between clusters suggests an evolutionary model based on CRD domain gene duplication and sequence diversification. In particular CiLgals-b N-CRD and C-CRD were similar to each other and both grouped with the ascidian C. picta mono-CRD. Homology modeling process shows a CiLgals molecular structure superimposed to chicken and mouse galectins. The CiLgals-a and CiLgals-b genes were upregulated by LPS inoculation suggesting that they are inducible and expressed in the inflamed pharynx as revealed by real-time PCR analysis. Finally, in situ hybridization and immunohistochemical assays showed their localization in the inflamed tissues, while immunoblotting analysis indicated that CiLgals can form oligomers.     

Isolation and characterization of microsatellite loci from the Japanese conger eel Conger myriaster

Five polymorphic microsatellite markers were developed in the Japanese conger eel (Conger myriaster) to understand the recruitment mechanisms and population structure of this species, the spawning sites and migration routes of which are unknown. Three of the loci were highly polymorphic and will serve as powerful tools for detailed population studies of the Japanese conger eel. One of the other two loci showed significant departure from Hardy–Weinberg equilibrium and the other showed low allelic diversity.     

Adaptive gene expression divergence inferred from population genomics

Detailed studies of individual genes have shown that gene expression divergence often results from adaptive evolution of regulatory sequence. Genome-wide analyses, however, have yet to unite patterns of gene expression with polymorphism and divergence to infer population genetic mechanisms underlying expression evolution. Here, we combined genomic expression data—analyzed in a phylogenetic context—with whole genome light-shotgun sequence data from six Drosophila simulans lines and reference sequences from D. melanogaster and D. yakuba. These data allowed us to use molecular population genetics to test for neutral versus adaptive gene expression divergence on a genomic scale. We identified recent and recurrent adaptive evolution along the D. simulans lineage by contrasting sequence polymorphism within D. simulans to divergence from D. melanogaster and D. yakuba. Genes that evolved higher levels of expression in D. simulans have experienced adaptive evolution of the associated 3′ flanking and amino acid sequence. Concomitantly, these genes are also decelerating in their rates of protein evolution, which is in agreement with the finding that highly expressed genes evolve slowly. Interestingly, adaptive evolution in 5′ cis-regulatory regions did not correspond strongly with expression evolution. Our results provide a genomic view of the intimate link between selection acting on a phenotype and associated genic evolution.