Evolutionary diversification of multigene families: allelic selection of toxins in predatory cone snails

Oxidative folding of conotoxins sharing an identical disulfide bridging framework

Conotoxins are short, disulfide-rich peptide neurotoxins produced in the venom of predatory marine cone snails. It is generally accepted that an estimated 100,000 unique conotoxins fall into only a handful of structural groups, based on their disulfide bridging frameworks. This unique molecular diversity poses a protein folding problem of relationships between hypervariability of amino acid sequences and mechanism(s) of oxidative folding. In this study, we present a comparative analysis of the folding properties of four conotoxins sharing an identical pattern of cysteine residues forming three disulfide bridges, but otherwise differing significantly in their primary amino acid sequence. Oxidative folding properties of M-superfamily conotoxins GIIIA, PIIIA, SmIIIA and RIIIK varied with respect to kinetics and thermodynamics. Based on rates for establishing the steady-state distribution of the folding species, two distinct folding mechanisms could be distinguished: first, rapid-collapse folding characterized by very fast, but low-yield accumulation of the correctly folded form; and second, slow-rearrangement folding resulting in higher accumulation of the properly folded form via the reshuffling of disulfide bonds within folding intermediates. Effects of changing the folding conditions indicated that the rapid-collapse and the slow-rearrangement mechanisms were mainly determined by either repulsive electrostatic or productive noncovalent interactions, respectively. The differences in folding kinetics for these two mechanisms were minimized in the presence of protein disulfide isomerase. Taken together, folding properties of conotoxins from the M-superfamily presented in this work and from the O-superfamily published previously suggest that conotoxin sequence diversity is also reflected in their folding properties, and that sequence information rather than a cysteine pattern determines the in vitro folding mechanisms of conotoxins.

Molecular diversity and accelerated evolution of C-type lectin-like proteins from snake venom

A number of C-type lectin-like proteins that affect thrombosis and hemostasis by inhibiting or activating specific platelet membrane receptors or blood coagulation factors have been isolated from the venom of various snake species and characterized and more than 80 have been sequenced. Recent data on the primary sequences and 3D structures of C-type lectins and C-type lectin-like proteins from snake venoms have enabled us to analyze their molecular evolution. Statistical analysis of their cDNA sequences shows that C-type lectin-like proteins, with some exceptions, have evolved in an accelerated manner to acquire their diverse functions. Phylogenetic analysis shows that the A and B chains of C-type lectin-like proteins are clearly separated from C-type lectins and that the A and B chains are further divided into a group of platelet receptor-binding proteins and a group of coagulation factor-binding proteins. Elucidation of the tertiary structures of several C-type lectin-like proteins led to the discovery of a unique domain-swapping interaction between heterodimeric subunits, which creates a concave surface for ligand binding.

Were arachnids the first to use combinatorial peptide libraries?

Spiders, scorpions, and cone snails are remarkable for the extent and diversity of gene-encoded peptide neurotoxins that are expressed in their venom glands. These toxins are produced in the form of structurally constrained combinatorial peptide libraries in which there is hypermutation of essentially all residues in the mature-toxin sequence with the exception of a handful of strictly conserved cysteines that direct the three-dimensional fold of the toxin. This gene-based combinatorial peptide library strategy appears to have been first implemented by arachnids almost 400 million years ago, long before cone snails evolved a similar mechanism for generating peptide diversity.

Patterns of diversifying selection in the phytotoxin-like scr74 gene family of Phytophthora infestans

Phytophthora infestans, the organism responsible for the Irish famine, causes late blight, a re-emerging disease of potato and tomato. Little is known about the molecular evolution of P. infestans genes. To identify candidate effector genes (virulence or avirulence genes) that may have co-evolved with the host, we mined expressed sequence tag (EST) data from infection stages of P. infestans for secreted and potentially polymorphic genes. This led to the identification of scr74, a gene that encodes a predicted 74-amino acid secreted cysteine-rich protein with similarity to the Phytophthora cactorum phytotoxin PcF. The expression of scr74 was upregulated approximately 60-fold 2 to 4 days after inoculation of tomato and was also significantly induced during early stages of colonization of potato. The scr74 gene was found to belong to a highly polymorphic gene family within P. infestans with 21 different sequences identified. Using the approximate and maximum likelihood (ML) methods, we found that diversifying selection likely caused the extensive polymorphism observed within the scr74 gene family. Pairwise comparisons of 17 scr74 sequences revealed elevated ratios of nonsynonymous to synonymous nucleotide-substitution rates, particularly in the mature region of the proteins. Using ML, all 21 polymorphic amino acid sites were identified to be under diversifying selection. Of these 21 amino acids, 19 are located in the mature protein region, suggesting that selection may have acted on the functional portions of the proteins. Further investigation of gene copy number and organization revealed that the scr74 gene family comprises at least three copies located in a region of no more than 300 kb of the P. infestans genome. We found evidence that recombination contributed to sequence divergence within at least one gene locus. These results led us to propose an evolutionary model that involves gene duplication and recombination, followed by functional divergence of scr74 genes. This study provides support for using diversifying selection as a criterion for identifying candidate effector genes from sequence databases.

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.

Functional divergence in tandemly duplicated Arabidopsis thaliana trypsin inhibitor genes

In multigene families, variation among loci and alleles can contribute to trait evolution. We explored patterns of functional and genetic variation in six duplicated Arabidopsis thaliana trypsin inhibitor (ATTI) loci. We demonstrate significant variation in constitutive and herbivore-induced transcription among ATTI loci that show, on average, 65% sequence divergence. Significant variation in ATTI expression was also found between two molecularly defined haplotype classes. Population genetic analyses for 17 accessions of A. thaliana showed that six ATTI loci arranged in tandem within 10 kb varied 10-fold in nucleotide diversity, from 0.0009 to 0.0110, and identified a minimum of six recombination events throughout the tandem array. We observed a significant peak in nucleotide and indel polymorphism spanning ATTI loci in the interior of the array, due primarily to divergence between the two haplotype classes. Significant deviation from the neutral equilibrium model for individual genes was interpreted within the context of intergene linkage disequilibrium and correlated patterns of functional differentiation. In contrast to the outcrosser Arabidopsis lyrata for which recombination is observed even within ATTI loci, our data suggest that response to selection was slowed in the inbreeding, annual A. thaliana because of interference among functionally divergent ATTI loci.

Gene expression and feeding ecology: evolution of piscivory in the venomous gastropod genus Conus

Differential expression of gene-family members is typically associated with the specific development of certain tissues and organs, but its importance in the ecological adaptation of organisms has rarely been investigated. Several specialized feeding modes have evolved within the predatory marine gastropod genus Conus, including molluscivory and piscivory. Based on phylogenetic investigations of Conus species, it has been concluded that piscivory arose at least twice in this genus. Moreover, molecular analyses of conotoxin mRNA transcripts reveal that piscivores from independent evolutionary lineages express the same subset of four-loop conotoxins, contrary to phylogenetic expectations. These results demonstrate that differential expression of gene-family members can play a key role in adaptive evolution, particularly during shifts to new ecological niches.

Physico-chemical characterization and synthesis of neuronally active alpha-conotoxins

The high specificity of alpha-conotoxins for different neuronal nicotinic acetylcholine receptors makes them important probes for dissecting receptor subtype selectivity. New sequences continue to expand the diversity and utility of the pool of available alpha-conotoxins. Their identification and characterization depend on a suite of techniques with increasing emphasis on mass spectrometry and microscale chromatography, which have benefited from recent advances in resolution and capability. Rigorous physico-chemical analysis together with synthetic peptide chemistry is a prerequisite for detailed conformational analysis and to provide sufficient quantities of alpha-conotoxins for activity assessment and structure-activity relationship studies.

Structurally conserved α-neurotoxin genes encode functionally diverse proteins in the venom of Naja sputatrix 1

The structure and organization of the genes encoding the long-chain neurotoxins and four other isoforms of weak neurotoxins in the venom of Naja sputatrix are reported. The genes contained three exons interrupted by two introns, a structure similar to other members of the three-finger toxin family. The proteins encoded by these genes, however, show varied affinity towards nicotinic acetylcholine receptors. Phylogenetic analysis of these genes showed that the weak neurotoxin gene is confined to a distinct group. We also observe that specific mutations of the gene provide the diversity in function in these toxins while maintaining a common structural scaffold. This forms the first report where the molecular basis of evolution of postsynaptic neurotoxins from an ancestral gene can be demonstrated using the same species of snake.

cDNA sequence and in vitro folding of GsMTx4, a specific peptide inhibitor of mechanosensitive channels

The peptide GsMTx4 from the tarantula venom (Grammostola spatulata) inhibits mechanosensitive ion channels. In this work, we report the cDNA sequence encoding GsMTx4. The gene is translated as a precursor protein of 80 amino acids. The first 21 amino acids are a predicted signal sequence and the C-terminal residues are a signal for amidation. An arginine residue adjacent to the N-terminal glycine of GsMTx4 is the cleavage site for release. The resulting peptide is 34 amino acids in length with a C-terminal phenylalanine and not a serine-alanine previously identified [J. Gen. Physiol. 115 (2000) 583]. We chemically synthesized this peptide and folded it in 0.1 M Tris, pH 7.9 with oxidized/reduced glutathione (1/10). Properties of the synthetic peptide were identical to the wild type for high performance liquid chromatography (HPLC), mass spectrometry, CD, and NMR. We also cloned GsMTx4 in a thioredoxin fusion protein system containing six histidines. Nickel affinity columns allowed rapid purification and folding occurred in conditions described above with 0.5 M guanidiniumHCl present. Thrombin cleavage liberated GsMTx4 with three extra amino acids at the N-terminus. The retention time in HPLC analysis and the CD spectrum was similar to wild type. Both the synthetic and cloned peptides were active in the patch clamp assay.

Population genetics of tandem trypsin inhibitor genes in Arabidopsis species with contrasting ecology and life history

Duplicated genes are important in the evolution and ecology of plant-defences because herbivore and pathogen attack can be countered via functional diversification at two levels: among duplicated loci and within loci. We explore molecular sequence variation for two members of a defence-related gene family, Arabidopsis thaliana trypsin inhibitors (ATTI), in A. thaliana and a closely related species, A. lyrata subspp. petraea. A worldwide sample of the inbreeding annual A. thaliana had less genetic variation at two ATTI loci (piTOTAL <or= 0.0006) than observed previously at other functional loci. A significant excess of high frequency derived alleles in the signal sequence and 5' UTR of ATTI2 was consistent with a model of positive selection. However, demographic processes such as population subdivision and expansion, both likely to have occurred in A. thaliana during the last 10 000 years, can also give rise to similar deviations from neutrality. A single population of A. lyrata subspp. petraea in Germany had up to an order of magnitude more standing genetic variation at ATTI loci than the species-wide sample of A. thaliana. Although the level of variability for ATTI1 and ATTI2 within this single population was similar to, or even greater than, observed species-wide diversity for other loci in A. lyrata, there was little evidence to reject an equilibrium neutral model. A spatially explicit sample of 87 A. lyrata subspp. petraea individuals detected outbreeding (FIS = -0.16; FIT = -0.15) but little population subdivision (FST = 0.006) in this self-incompatible perennial herb. Genetic differences between Arabidopsis species were consistent with, but not fully explained by, divergence in ecology and life history. Diversification appears to have occurred in different functional domains for the tandemly duplicated ATTI1 and ATTI2 genes; the majority of fixed replacements in ATTI1 surround the enzyme binding site of the mature protein, whereas in ATTI2 most functional evolutionary change is located in the signal peptide. This pattern is consistent with a hypothesis of subfunctionalization in trypsin inhibitory function.

Purifying selection and birth-and-death evolution in the histone H4 gene family

Histones are small basic proteins encoded by a multigene family and are responsible for the nucleosomal organization of chromatin in eukaryotes. Because of the high degree of protein sequence conservation, it is generally believed that histone genes are subject to concerted evolution. However, purifying selection can also generate a high degree of sequence homogeneity. In this study, we examined the long-term evolution of histone H4 genes to determine whether concerted evolution or purifying selection was the major factor for maintaining sequence homogeneity. We analyzed the proportion (p(S)) of synonymous nucleotide differences between the H4 genes from 59 species of fungi, plants, animals, and protists and found that p(S) is generally very high and often close to the saturation level (p(S) ranging from 0.3 to 0.6) even though protein sequences are virtually identical for all H4 genes. A small proportion of genes showed a low level of p(S) values, but this appeared to be caused by recent gene duplication. Our findings suggest that the members of this gene family evolve according to the birth-and-death model of evolution under strong purifying selection. Using histone-like genes in archaebacteria as outgroups, we also showed that H1, H2A, H2B, H3, and H4 histone genes in eukaryotes form separate clusters and that these classes of genes diverged nearly at the same time, before the eukaryotic kingdoms diverged.

Cone venom--from accidental stings to deliberate injection

Cone snails have long been of note due to their colorful shells and deadly venom. Over the years, a number of people who have encountered these molluscs have been injured or killed by their sting. Biochemical analysis of the venom components has revealed a plethora of peptides and proteins that target a variety of receptors and ion channels. Pharmaceutical companies are now utilizing the selectivity and potency of Conus-derived peptides to develop novel medications for pain, epilepsy and other disorders.