Specialisation of the venom gland proteome in predatory cone snails reveals functional diversification of the conotoxin biosynthetic pathway

Conotoxins, venom peptides from marine cone snails, diversify rapidly as speciation occurs. It has been suggested that each species can synthesize between 1000 and 1900 different toxins with little to no interspecies overlap. Conotoxins exhibit an unprecedented degree of post-translational modifications, the most common one being the formation of disulfide bonds. Despite the great diversity of structurally complex peptides, little is known about the glandular proteins responsible for their biosynthesis and maturation. Here, proteomic interrogations on the Conus venom gland led to the identification of novel glandular proteins of potential importance for toxin synthesis and secretion. A total of 161 and 157 proteins and protein isoforms were identified in the venom glands of Conus novaehollandiae and Conus victoriae, respectively. Interspecies differences in the venom gland proteomes were apparent. A large proportion of the proteins identified function in protein/peptide translation, folding, and protection events. Most intriguingly, however, we demonstrate the presence of a multitude of isoforms of protein disulfide isomerase (PDI), the enzyme catalyzing the formation and isomerization of the native disulfide bond. Investigating whether different PDI isoforms interact with distinct toxin families will greatly advance our knowledge on the generation of cone snail toxins and disulfide-rich peptides in general.     

Conservation and divergence in multigene families: alternatives to selection and drift

It is generally assumed that conservation and divergence of DNA signify function (selection) and no function (drift), respectively. This assumption is based on the view that a mutation is a unique event on a single chromosome, the fate of which depends on selection or drift. Knowledge of the rates, units and biases of widespread mechanisms of non-reciprocal DNA exchange, in particular within multigene families, provides alternative explanations for conservation and divergence, notwithstanding biological function. Such mechanisms of DNA turnover cause continual fluctuations in the copy-number of variant genes in an individual and, hence, promote the gradual and cohesive spread of a variant gene throughout a family (homogenization) and throughout a population (fixation). The dual processes (molecular drive) of homogenization and fixation are inextricably linked. Data are presented of the expected stages of transition in the spread of variant repeats by molecular drive in some non-genic families of DNA, seemingly not under the influence of selection. When a molecularly driven change in a given gene family is accompanied by the coevolution (mediated by selection) of other DNA, RNA or protein molecules that interact with the gene family then biological function is observed to be maintained despite sequence divergence. Conversely, the mechanics of DNA turnover and a turnover bias in favour of ancestral sequences can dramatically retard the rate of sequence change, in the absence of function. Examples of the maintenance of function by molecular coevolution and conservation of sequences in the absence of function, are drawn mainly from the rDNA multigene family.     

Evolutionary diversification in polyamine biosynthesis

Polyamine biosynthesis is an ancient metabolic pathway present in all organisms. Aminopropyltransferases are key enzymes that mediate the synthesis of spermidine, spermine, and thermospermine. The relatively high sequence similarity between aminopropyltransferases and their similarity with putrescine N-methyltransferases (PMT) raises the question of whether they share a common ancestor or have evolved by convergence. Here we show that aminopropyltransferases and PMT are phylogenetically interconnected, and the different activities have been generated by unusually frequent events of diversification of existing functions. Although all spermidine synthases (SPDSs) derive from a common ancestor preceding the separation between prokaryotes and eukaryotes, they have been the origin of a variety of new activities. Among those, spermine synthases (SPMSs) represent a novelty independently arisen at least 3 times, in animals, fungi, and plants. The most parsimonious mechanism would involve the duplication and change of function of preexisting SPDS genes in each phylum. Although spermine is not essential for life, the repeated invention of SPMS and its conservation strongly argues for an evolutionary advantage derived from its presence. Moreover, the appearance of thermospermine synthase (tSPMS) in several genera of Archaea and Bacteria was accompanied by a loss of SPDS, suggesting that the new activity originated as a change of function of this enzyme. Surprisingly, tSPMS was later acquired by plants at an early stage of evolution by horizontal gene transfer and has proven to be essential for vascular development in tracheophytes. Finally, the synthesis of nicotine and tropane alkaloids in Solanales was favored by the origination of a new activity, PMT, as a duplication and change of function from SPDS.     

Malaria multigene families: the price of chronicity

In this article, Georges Snounou, William Jarra and Peter Preiser discuss the survival strategy of malaria parasites in the light of a novel mechanism of clonal phenotypic variation recently described for a multigene family of Plasmodium yoelii yoelii. The 235 kDa rhoptry proteins (Py235) encoded by these genes may be involved in the selection of red blood cells for invasion by merozoites. The new mechanism may explain the ability of individual parasites to adapt to natural variations in red blood cell subsets, while ensuring that sufficient merozoites escape immune attack, thus maintaining a chronic infection for extended periods. This counterpoints the antigenic variation exemplified by PfEMP1 proteins (a large family of proteins derived from P. falciparum), which operates at the population level. The possibility of manipulating the expression of functionally similar genes in other Plasmodium species could lead to therapies aimed at reducing clinical severity without compromising the acquisition and maintenance of immunity.     

The number of alleles in multigene families

The probability distribution and moments of the number of alleles present in a sample of homologous chromosomes are studied. It is assumed that there are multiple copies of the gene on each chromosome. When there are only two copies per chromosome or when there are only two or three chromosomes, it is possible to use analytic methods to tackle the problem. Otherwise, a simulation method is suggested.     

On the evolution of multigene families

Multigene families are classified into three groups: small families as exemplified by hemoglobin genes of mammals; middlesize multigene families, by genes of mammalian histocompatibility antigens; and large multigene families, by variable region genes of immunoglobulins. Facts and theories on these evolving multigene families are reviewed, with special reference to the population genetics of their concerted evolution. It is shown that multigene families are evolving under continued occurrence of unequal (but homologous) crossing-over and gene conversion, and that mechanisms for maintaining genetic variability are totally different from the conventional models of population genetics. Thus, in view of widespread occurrence of multigene families in genomes of higher organisms, the evolutionary theory based mainly on change of gene frequency at each locus would appear to need considerable revision.     

Inference of selection based on temporal genetic differentiation in the study of highly polymorphic multigene families

The co-evolutionary arms race between host immune genes and parasite virulence genes is known as Red Queen dynamics. Temporal fluctuations in allele frequencies, or the 'turnover' of alleles at immune genes, are concordant with predictions of the Red Queen hypothesis. Such observations are often taken as evidence of host-parasite co-evolution. Here, we use computer simulations of the Major Histocompatibility Complex (MHC) of guppies (Poecilia reticulata) to study the turnover rate of alleles (temporal genetic differentiation, G'(ST)). Temporal fluctuations in MHC allele frequencies can be ≥≤order of magnitude larger than changes observed at neutral loci. Although such large fluctuations in the MHC are consistent with Red Queen dynamics, simulations show that other demographic and population genetic processes can account for this observation, these include: (1) overdominant selection, (2) fluctuating population size within a metapopulation, and (3) the number of novel MHC alleles introduced by immigrants when there are multiple duplicated genes. Synergy between these forces combined with migration rate and the effective population size can drive the rapid turnover in MHC alleles. We posit that rapid allelic turnover is an inherent property of highly polymorphic multigene families and that it cannot be taken as evidence of Red Queen dynamics. Furthermore, combining temporal samples in spatial F(ST) outlier analysis may obscure the signal of selection.     

Linkage and evolutionary diversification of developmentally regulated multigene families: tandem arrays of the 401/18 chorion gene pair in silkmoths.

The coordinately expressed silkmoth chorion genes, 401 and 18, are closely linked as a pair, in divergent orientation. Analysis of overlapping clones (chromosomal "walk") demonstrated that each of the multiple copies of this gene pair is embedded within a larger deoxyribonucleic acid unit, which is tandemly repeated in a few arrays or possibly a single array. Southern analysis and examination of clones from a single individual moth demonstrated that the repeat units are extensively polymorphic in restriction sites, length, and possibly number, no differential amplification was evident during choriogenesis. Intron and 5'-flanking sequences were shown to be specific for the 401/18 gene pair and not to be present elsewhere in the genome. The spatial distribution of variations in the genes and their flanking sequences were examined.     

Evolutionary diversification of the flowers in angiosperms

Angiosperms and their flowers have greatly diversified into an overwhelming array of forms in the past 135 million years. Diversification was shaped by changes in climate and the biological environment (vegetation, interaction with other organisms) and by internal structural constraints and potentials. This review focuses on the development and structural diversity of flowers and structural constraints. It traces floral diversification in the different organs and organ complexes (perianth, androecium, gynoecium) through the major clades of extant angiosperms. The continuously improved results of molecular phylogenetics provide the framework for this endeavor, which is necessary for the understanding of the biology of the angiosperms and their flowers. Diversification appears to work with innovations and modifications of form. Many structural innovations originated in several clades and in special cases could become key innovations, which likely were hot spots of diversification. Synorganization between organs was an important process to reach new structural levels, from which new diversifications originated. Complexity of synorganization reached peaks in Orchidaceae and Apocynaceae with the independent evolution of pollinaria. Such a review throughout the major clades of angiosperms also shows how superficial and fragmentary our knowledge on floral structure in many clades is. Fresh studies and a multidisciplinary approach are needed.     

Putative gamma-conotoxins in vermivorous cone snails: the case of Conus delessertii

Peptide de7a was purified from the venom of Conus delessertii, a vermivorous cone snail collected in the Yucatan Channel, Mexico. Its amino acid sequence was determined by automatic Edman degradation after reduction and alkylation. The sequence shows six Cys residues arranged in the pattern that defines the O-superfamily of conotoxins, and several post-translationally modified residues. The determination of its molecular mass by means of laser desorption ionization time-of-flight mass spectrometry (average mass, 3170.0 Da) confirmed the chemical data and suggested amidation of the C-terminus. The primary structure (ACKOKNNLCAITgammaMAgammaCCSGFCLIYRCS*; O, hydroxyproline; gamma, gamma-carboxyglutamate; *, amidated C-terminus; calculated average mass, 3169.66 Da) of de7a contains a motif (gammaCCS) that has previously only been found in two other toxins, both from molluscivorous cone snails: TxVIIA from Conus textile and gamma-PnVIIA from Conus pennaceus. These toxins cause depolarization and increased firing of action potentials in molluscan neuronal systems, and toxin gamma-PnVIIA has been shown to act as an agonist of neuronal pacemaker cation currents. The similarities to toxins TxVIIA and gamma-PnVIIA suggest that peptide de7a might also affect voltage-gated nonspecific cation pacemaker channels.     

On the divergence of genes in multigene families

Statistical properties of the amount of divergence of genes in multigene families are studied. The model considered is an infinite-site neutral model with unbiased intrachromosomal conversion, unbiased interchromosomal conversion, and recombination. By considering the time back to the most recent common ancestor of two genes, both the probability of identity and the moments of S, the number of sites that differ between two sampled genes, are obtained. We find that if recombination rates are large or conversion is always interchromosomal, then the expectation of S is 4N mu n where N is the population size, mu is the rate of mutation per generation per gene and n is the number of genes in the gene family, as the conversion rates approach zero, the moments of divergence do not approach the moments of divergence with conversion rates equal to zero, and it is possible for a decrease in the rate of intrachromosomal conversion to result in a higher probability of identity, but a greater mean divergence of the two genes.     

Choosing among alternative trees of multigene families

Estimation of gene trees is the first step in testing alternative hypotheses about the evolution of multigene families. The standard practice for inferring gene family history is to construct trees that meet some objective criteria based on the fit of the character state changes (nucleotide or amino acid changes) to the gene tree. Unfortunately, analysis of character state data can be misleading. In addition, this approach ignores information about the relationships of the species from which the genes have been sampled. In this paper I explore using statistics of fit between the character data and gene trees and the reconciliation of the gene and species trees for choosing among alternative evolutionary hypotheses of gene families. In particular, I advocate a two-pronged strategy for choosing among alternative gene trees. First, the character data are used to define a set of acceptable gene trees (i.e., trees that are not significantly different from the minimum length tree). Next, the set of acceptable gene trees is reconciled with a known species tree, and the gene tree requiring the fewest number of gene duplications and losses is adopted as the best estimate of evolutionary history. The approach is illustrated using three gene families: BMP, EGR, and LDH.     

Adaptive evolution of animal toxin multigene families

Animal toxins comprise a diverse array of proteins that have a variety of biochemical and pharmacological functions. A large number of animal toxins are encoded by multigene families. From studies of several toxin multigene families at the gene level the picture is emerging that most have been functionally diversified by gene duplication and adaptive evolution. The number of pharmacological activities in most toxin multigene families results from their adaptive evolution. The molecular evolution of animal toxins has been analysed in some multigene families, at both the intraspecies and interspecies levels. In most toxin multigene families, the rate of non-synonymous to synonymous substitutions (dN/dS) is higher than one. Thus natural selection has acted to diversify coding sequences and consequently the toxin functions. The selection pressure for the rapid adaptive evolution of animal toxins is the need for quick immobilization of the prey in classical predator and prey interactions. Currently available evidence for adaptive evolution in animal toxin multigene families will be considered in this review.     

Evolutionary diversification of DNA methyltransferases in eukaryotic genomes

In eukaryotes, C5-cytosine methylation is a common mechanism associated with a variety of functions such as gene regulation or control of genomic stability. Different subfamilies of eukaryotic methyltransferases (MTases) have been identified, mainly in metazoa, plants, and fungi. In this paper, we used hidden Markov models to detect MTases in completed or almost completed eukaryotic genomes, including different species of Protozoa. A phylogenetic analysis of MTases enabled us to define six subfamilies of MTases, including two new subfamilies. The dnmt1 subfamily that includes all the known MTases with a maintenance activity seems to be absent in the Protozoa. The dnmt2 subfamily seems to be the most widespread, being present even in the nonmethylated Dictyostelium discoideum. We also found two dnmt2 members in the bacterial genus Geobacter, suggesting that horizontal transfers of MTases occurred between eukaryotes and prokaryotes. Even if the direction of transfer cannot be determined, this relationship might be useful for understanding the function of this enigmatic subfamily of MTases. Globally, our analysis reveals a great diversity of MTases in eukaryotes, suggesting the existence of different methylation systems. Our results also suggest acquisitions and losses of different MTases in every eukaryotic lineage studied and that some eukaryotes appear to be devoid of methylation.     

Evolutionary diversification of protein-coding genes of hantaviruses

Phylogenetic analyses of the S:, M, and L: genes of the hantaviruses (Bunyaviridae: Hantavirus) revealed three well-differentiated clades corresponding to viruses parasitic on three subfamilies (Murinae, Arvicolinae, and Sigmodontinae) of the rodent family Muridae. In rooted trees of M: and L: genes, the viruses with hosts belonging to Murinae formed an outgroup to those with hosts in Arvicolinae and Sigmodontinae. This phylogeny corresponded with a phylogeny of the murid subfamilies based on mitochondrial cytochrome b sequences, supporting the hypothesis that hantaviruses have coevolved with their mammalian hosts at least since the common ancestor of these three subfamilies, which probably occurred about 50 MYA. The nucleocapsid protein (encoded by the S: gene) differentiated among the viruses parasitic on the three subfamilies in such a way that a high frequency of amino acid residue charge changes occurred in a hypervariable (HV) portion of the molecule, and nonsynonymous nucleotide differences causing amino acid charge changes in the HV region occurred significantly more frequently than expected under random substitution. Along with evidence that at least in some hantaviruses the HV region is a target for host antibodies and the known importance of charged residues in determining antibody epitopes, these results suggest that changes in the HV region may represent adaptation to host-specific characteristics of the immune response.     

Toxins from cone snails: properties,applications and biotechnological production

Cone snails are marine predators that use venoms to immobilize their prey. The venoms of these mollusks contain a cocktail of peptides that mainly target different voltage- and ligand-gated ion channels. Typically, conopeptides consist of ten to 30 amino acids but conopeptides with more than 60 amino acids have also been described. Due to their extraordinary pharmacological properties, conopeptides gained increasing interest in recent years. There are several conopeptides used in clinical trials and one peptide has received approval for the treatment of pain. Accordingly, there is an increasing need for the production of these peptides. So far, most individual conopeptides are synthesized using solid phase peptide synthesis. Here, we describe that at least some of these peptides can be obtained using prokaryotic or eukaryotic expression systems. This opens the possibility for biotechnological production of also larger amounts of long chain conopeptides for the use of these peptides in research and medical applications.