Positive selection is a general phenomenon in the evolution of abalone sperm lysin

Lysin is a 16kDa acrosomal protein used by abalone sperm to create a hole in the egg vitelline envelope (VE). The interaction of lysin with the VE is species-selective and is one step in the multistep fertilization process that restricts heterospecific (cross-species) fertilization. For this reason, the evolution of lysin could play a role in establishing prezygotic reproductive isolation between species. Previously, we sequenced sperm lysin cDNAs from seven California abalone species and showed that positive Darwinian selection promotes their divergence. In this paper an additional 13 lysin sequences are presented representing species from Japan, Taiwan, Australia, New Zealand, South Africa, and Europe. The total of 20 sequences represents the most extensive analysis of a fertilization protein to date. The phylogenetic analysis divides the sequences into two major clades, one composed of species from the northern Pacific (California and Japan) and the other composed of species from other parts of the world. Analysis of nucleotide substitution demonstrates that positive selection is a general process in the evolution of this fertilization protein. Analysis of nucleotide and codon usage bias shows that neither parameter can account for the robust data supporting positive selection. The selection pressure responsible for the positive selection on lysin remains unknown.     

The eosinophil ribonucleases

The eosinophil ribonucleases, eosinophilderived neurotoxin (EDN/RNase 2) and eosinophil cationic protein (ECP/RNase 3) are two closely related proteins with intriguing functional and evolutionary properties. While both EDN and ECP maintain the structural and catalytic residues typical of the RNase A superfamily, the role of ribonuclease activity in the physiologic function of these proteins remains unclear. The biochemistry and physiology of EDN, ECP and the recently discovered ribonuclease k6 (RNase 6) will be reviewed in this chapter.     

Positive Darwinian selection in Vipera palaestinae phospholipase A2 genes is unexpectedly limited to the third exon

The venom of Vipera palaestinae contains a two-component toxin, consisting of an acidic phospholipase A2 (PLA2) and a basic protein. Here we report the cloning and sequence analysis of the complete V. palaestinae PLA2 genes. Since in all Viperidae PLA2 multigene families the 5' and 3' flanking regions are highly conserved, we designed oligonucleotide primers that allow amplification of the whole PLA2 multigene family in a single step. The structural organization of both genes is the same as in the Vipera ammodytes PLA2 multigene family, there being five exons separated by four introns. Comparison of V. palaestinae PLA2 genes with other Viperidae PLA2 genes has shown that the structural organization of the genes and the nucleotide sequence of all introns and flanking regions are highly conserved, whereas the third exon clearly shows a higher number of amino acid replacements, an indication of positive Darwinian selection. The positive Darwinian selection is surprisingly limited to the third exon, in contrast to other Viperidae PLA2 genes, where it is present in all mature protein coding exons.     

Toxicity of an antitumor ribonuclease to Purkinje neurons

Purkinje cell toxicity is one of the characteristic features of the Gordon phenomenon, a syndrome manifested by ataxia, muscular rigidity, paralysis, and tremor that may lead to death (Gordon, 1933). Two members of the RNase superfamily found in humans, EDN (eosinophil-derived neurotoxin) and ECP (eosinophil cationic protein), cause the Gordon phenomenon when injected intraventricularly into guinea pigs or rabbits. We have found that another member of the RNase superfamily, an antitumor protein called onconase, isolated from Rana pipiens oocytes and early embryos, will also cause the Gordon phenomenon when injected into the cerebrospinal fluid of guinea pigs at a dose similar to that of EDN (LD50, 3-4 micrograms). Neurologic abnormalities of onconase-treated animals were indistinguishable from those of EDN-treated animals, and histology showed dramatic Purkinje cell loss in the brains of onconase-treated animals. The neurotoxic activity of onconase correlates with ribonuclease activity. Onconase modified by iodoacetic acid to eliminate 70% and 98% of the ribonuclease activity of the native enzyme displays a similar decrease in ability to cause the Gordon phenomenon. In contrast, the homologous bovine pancreatic RNase A injected intraventricularly at a dose 5000 times greater than the LD50 dose of EDN or onconase is not toxic and does not cause the Gordon phenomenon. A comparison of the RNase activities of EDN, onconase, and bovine pancreatic RNase A using three pancreatic RNA substrates demonstrates that onconase is orders of magnitude less active enzymatically than EDN and RNase A. Thus, another member of the RNase superfamily in addition to EDN and ECP can cause the Gordon phenomenon.(ABSTRACT TRUNCATED AT 250 WORDS)     

Serum levels of eosinophil cationic protein and eosinophil protein x in pollen atopic patients with stable asthma and its relation with bronchial hyperresponsiveness

Eosinophils are important effector cells in allergic inflammation described in allergic rhinitis (AR) and allergic bronchial asthma (BA). During the pollen season serum levels of eosinophil cationic protein (ECP) and eosinophil X protein/eosinophil-derived neurotoxin (EPX/EDN) are increased in BA. The aim of the present study was to evaluate the serum levels of ECP and EPC in pollen atopic patients with AR and BA during the winter. 92 patients were studied. They were divided into three groups: I 29 patients with AR, II 51 patients with BA and III 12 healthy subjects. Allergic rhinitis and bronchial asthma were diagnosed by routine clinical tests: clinical history, skin tests, total IgE and specific IgE. In addition ECP and EPX were determined in serum. All patients were asymptomatic, stable and without medical treatment. Methacholine challenge test (MCT) was performed in all patients. MCT were positive in 4 patients of group I and 45 patients of group II. ECP levels (ug/l) were: 21 (I), 24 (II) and 7 (III). EPX levels (ug/l) were 35 (I), 45 (II) and 21 (III). Statistical differences (p < 0.01) were observed both in ECP and EPX levels in patients with MCT positive in relation to patients with MCT negative, and in allergic patients (I and II) in comparison with the healthy subjects (III) (p < 0.01). ECP and EPX serum levels are increased in patients with a positive MCT in the winter, out of the pollen season, when patients are asymptomatic, stable and without treatment. This fact suggests that eosinophils play an important role in the pathogenesis of bronchial asthma.     

Concerted evolution in an egg receptor for a rapidly evolving abalone sperm protein

Gamete interactions during fertilization exhibit species specificity. In abalone, the sperm protein lysin species-specifically creates a hole in the egg envelope. Lysin evolves rapidly by positive Darwinian selection. Evolution of the egg receptor for lysin provides the selective pressure for lysin's divergence. The egg receptor for lysin is a tandemly repeated sequence that evolves by concerted evolution. Concerted evolution in the egg receptor could explain the rapid, adaptive evolution in sperm lysin and may provide an underlying molecular mechanism that gives rise to species-specific fertilization.     

Molecular evolution of a duplication: the sex-peptide (Acp70A) gene region of Drosophila subobscura and Drosophila madeirensis

In Drosophila melanogaster, the Acp70A gene, which is involved in the postmating reactions of the female, is a single-copy gene. However, in Drosophila subobscura, the gene is duplicated and both copies are transcribed. To study the molecular evolution of the duplication, a 2.1-kb fragment encompassing both copies of the duplication was sequenced for 10 lines of D. subobscura and one line of Drosophila madeirensis. Estimates of the divergence between the two copies of the duplicated region and between the two species studied, D. subobscura and D. madeirensis, revealed that both copies of the Acp70a gene had evolved independently since their duplication. The ratio of nonsynonymous to silent divergence between copies was generally higher than one. The McDonald and Kreitman test revealed an excess of nonsynonymous changes fixed since the duplication and before the split of the D. subobscura and D. madeirensis lineages. These results point to natural selection driving protein evolution after the duplication. Specifically, adaptive evolution appears to have caused the initial differentiation between copies of the N-terminal parts of the proteins, while purifying selection could be responsible for the high conservation of the C-terminal parts.     

Reciprocal relationship between stem-loop potential and substitution density in retroviral quasispecies under positive Darwinian selection

Nucleic acids have the potential to form intrastrand stem-loops if complementary bases are suitably located. Computer analyses of poliovirus and retroviral RNAs have revealed a reciprocal relationship between "statistically significant" stem-loop potential and "sequence variability." The statistically significant stem-loop potential of a nucleic acid segment has been defined as a function of the difference between the folding energy of the natural segment (FONS) and the mean folding energy of a set of randomized (shuffled) versions of the natural segment (FORS-M). Since FONS is dependent on both base composition and base order, whereas FORS-M is solely dependent on base composition (a genomic characteristic), it follows that statistically significant stem-loop potential (FORS-D) is a function of base order (a local characteristic). In retroviral genomes, as in all DNA genomes studied, positive FORS-D values are widely distributed. Thus there have been pressures on base order both to encode specific functions and to encode stem-loops. As in the case of DNA genomes under positive Darwinian selection pressure, in HIV-1 specific function appears to dominate in rapidly evolving regions. Here high sequence variability, expressed as substitution density (not indel density), is associated with negative FORS-D values (impaired base-order-dependent stem-loop potential). This suggests that in these regions HIV-1 genomes are under positive selection pressure by host defenses. The general function of stem-loops is recombination. This is a vital process if, from among members of viral "quasispecies," functional genomes are to be salvaged. Thus, for rapidly evolving RNA genomes, it is as important to conserve base-order-dependent stem-loop potential as to conserve other functions.     

Population dynamics of sexual and resource competition

The organization and sequence of genes encoding the alpha 1-proteinase inhibitor (alpha 1PI), a major serine proteinase inhibitor of the mammalian bloodstream, have been compared in several species, including murine rodents (genus Mus). Analysis of gene copy number indicates that amplification of alpha 1PI genes occurred at some time during evolution of the Mus genus, leading to fixation of a family of about three to five genes in several existing species (e.g., M. domesticus and M. saxicola), and only a single gene in others (e.g., M. caroli). A phylogeny for the various mammalian alpha 1PI mRNAs was constructed based upon synonymous substitutions within coding regions. The mRNAs in different murine species diverged from a common ancestor before the formation of the first species lineages of the Mus genus, i.e., about 10-13 million years ago. Thus, alpha 1PI gene amplification must have occurred prior to Mus speciation; gene families were retained in some, but not all, murine species. The reactive center region of the alpha 1PI polypeptide, which determines target protease specificity, has diverged rapidly during evolution of the Mus species, but not during evolution of other mammalian species included in the analysis. It is likely that this accelerated evolution of the reactive center, which has been noted previously for serine proteinase inhibitors, was driven by some sort of a positive Darwinian selection that was exerted in a taxon-specific manner. We suggest that evolution of alpha 1PI genes of murine rodents has been characterized by both modification of gene copy number and rapid reactive center divergence. These processes may have resulted in a broadened repertoire of proteinase inhibitors that was evolutionarily advantageous during Mus speciation.     

Evolution of mammalian X-linked and autosomal Pgk and Pdh E1 alpha subunit genes

The phylogeny and substitution rates of the mammalian X chromosome-located and autosomal phosphoglycerate kinase and pyruvate dehydrogenase genes were investigated. Compatibility analysis was used to show reticulate evolution in these genes. Analysis of the marsupial, mouse, and human phosphoglycerate kinase genes suggests that at least two recombination events have taken place, one occurring about the time of the placental-marsupial split involving exons 1-5 and the other before the primate-rodent split involving exons 9-10. Similar analysis of the pyruvate dehydrogenase genes indicates a recombination event involving exons 2-3 at a time before the primate-rodent split and a gene conversion between exons 3-4 in the human somatic and testis-specific pyruvate dehydrogenase genes after the primate-rodent split. This demonstrates that genetic exchange can occur between paralogous genes at widely separated chromosomal locations. Estimation of nucleotide substitution rates in these genes confirmed a higher substitution rate in the pyruvate dehydrogenase genes. In the phosphoglycerate kinase genes, there is no difference between the substitution rates in mice and humans and between the X chromosome- and autosome-located genes. A greater substitution rate was noted in the mouse autosomal pyruvate dehydrogenase gene when compared with the other mouse and human genes. This may be a result of either directional natural selection or a relaxation of functional constraint at this specific gene.     

Molecular basis for the dominant white phenotype in the domestic pig

The change of phenotypic traits in domestic animals and crops as a response to selective breeding mimics the much slower evolutionary change in natural populations. Here, we describe that the dominant white phenotype in domestic pigs is caused by two mutations in the KIT gene encoding the mast/stem cell growth factor receptor (MGF), one gene duplication associated with a partially dominant phenotype and a splice mutation in one of the copies leading to the fully dominant allele. The splice mutation is a G to A substitution in the first nucleotide of intron 17 and leads to skipping of exon 17. The duplication is most likely a regulatory mutation affecting KIT expression, whereas the splice mutation is expected to cause a receptor with impaired or absent tyrosine kinase activity. Immunocytochemistry showed that this variant form is expressed in 17- to 19-day-old pig embryos. Hundreds of millions of white pigs around the world are assumed to be heterozygous or homozygous for the two mutations. [The EMBL accession numbers for porcine KIT1*0101, KIT1*0202, KIT2*0202, and KIT2*0101 are AJ223228-AJ223231, respectively.]     

Mycobacteriosis caused by Mycobacterium genavense in birds kept in a zoo: 11-year survey

Nucleotide sequences encoding GroEL homologs of intracellular symbionts in three closely interrelated aphids were compared with one another and that for GroEL. It was suggested that in these proteins a particular position is highly susceptible to amino acid substitution, through which the GroEL homologs of symbionts seemed to have acquired a unique function on top of the activity as molecular chaperone. This may represent a rare example of non-neutral evolution of molecule under the positive selection pressure.     

Methods for the detection of non-random base substitution in virus genes: models of synonymous nucleotide substitution in picornavirus genes

A substantial fraction of phylogenetic divergence between closely related RNA virus genes is generally accounted for by synonymous (non-amino acid changing) point mutation. Viral evolution may be a complicated phenomena, governed by many different processes. However in this study we ask whether there are any properties in the patterns of synonymous nucleotide substitutions in three different Picornavirus genes that permit the process of accumulation of synonymous point mutation in these genes to be distinguished from some of the simplest most basic evolutionary models. We conclude that while the observed patterns in the occurrence of synonymous point substitution are consistent with those predicted by a model in which base mutation is equi-probable along a gene, and the probability of synonymous substitution determined only by local codon usage, some patterns in the actual nucleotides exchanged remain to be explained.     

Evolution of indirect reciprocity by image scoring

Darwinian evolution has to provide an explanation for cooperative behaviour. Theories of cooperation are based on kin selection (dependent on genetic relatedness), group selection and reciprocal altruism. The idea of reciprocal altruism usually involves direct reciprocity: repeated encounters between the same individuals allow for the return of an altruistic act by the recipient. Here we present a new theoretical framework, which is based on indirect reciprocity and does not require the same two individuals ever to meet again. Individual selection can nevertheless favour cooperative strategies directed towards recipients that have helped others in the past. Cooperation pays because it confers the image of a valuable community member to the cooperating individual. We present computer simulations and analytic models that specify the conditions required for evolutionary stability of indirect reciprocity. We show that the probability of knowing the 'image' of the recipient must exceed the cost-to-benefit ratio of the altruistic act. We propose that the emergence of indirect reciprocity was a decisive step for the evolution of human societies.     

Identification and functional analysis of sulfonylurea receptor 1 variants in Japanese patients with NIDDM

The sulfonylurea receptor 1 (SUR1) is an essential regulatory subunit of the beta-cell ATP-sensitive K+ channel (K[ATP]). The possible role of SUR1 gene mutation(s) in the development of NIDDM remains controversial as both a positive association and negative linkage results have been reported. Therefore, we examined the SUR1 gene at the single nucleotide level with single strand conformation polymorphism analysis in 100 Japanese NIDDM patients. We identified a total of five amino acid substitutions and 17 silent mutations by examining all 39 exons of this gene. Two rare novel mutations, D811N in exon 20 and R835C in exon 21, were identified in the first nucleotide-binding fold (NBF), a functionally important region of SUR1, in one patient each, both heterozygotes. To analyze possible functional alterations, we reconstituted the mutant K(ATP) by coexpressing beta-cell inward rectifier (BIR) (Kir 6.2), a channel subunit of K(ATP), and mutant SUR1 in HEK293T and COS-7 cells. As demonstrated by the patch clamp technique and rubidium (Rb+) efflux studies, neither mutation alters the properties of channel activities. Two other rare missense mutations, R275Q in exon 6 and V560M in exon 12, were also identified. The R275Q substitution was not found in 67 control subjects, and V560M was present in three control subjects. Neither of these substitutions appeared to cosegregate with NIDDM in the probands' families. A previously reported S1370A substitution located in the second NBF was also common in the Japanese subjects (allelic frequency 0.37), and was found at an equal frequency in nondiabetic control subjects. In conclusion, SUR1 mutations impairing K(ATP) function do not appear to be major determinants of NIDDM susceptibility in Japanese.     

Correlates of recency of eye examination among elderly African-Americans

The E2 gene of the branched-chain alpha-keto acid dehydrogenase (BCKDH) complex was studied at the molecular level in three patients with intermittent maple syrup urine disease (MSUD). All three patients had higher BCKDH activity than did those with the classical phenotype. In the first patient, a single base substitution from A to G in intron 8 created a new 5' splice site and caused an insertion of 126 nucleotides between exons 8 and 9 by activating an upstream cryptic 3' splice site in the same intron. The predicted mRNA encoded a truncated protein with 282 amino acids including 4 novel ones at the carboxyl terminus, compared with the normal protein with 421 amino acids. In vitro, the region from the patient but not from a normal control was recognized and was recovered as a novel exon, indicating that the single substitution was responsible for incorporation of the region into mRNA. This mutation probably supports an exon definition model in which the spliceosome recognizes a 3' splice site and then scans downstream for an acceptable 5' splice site, thereby defining an exon. The second patient was homozygous for a G to T transversion at nucleotide 1463 in exon 11, which predicted a substitution of the termination codon by a leucine residue and the addition of 7 extra amino acids at the carboxyl terminus. For each mutation, these two patients were homozygous and their parents were heterozygous. The third patient was a compound heterozygote for a C to G transversion at nucleotide 309 in exon 4 and a G to A transition at nucleotide 1165 in exon 9, causing an Ile-to-Met substitution at amino acid 37 and a Gly-to-Ser substitution at amino acid 323, respectively. Taken together, these results indicate that the molecular basis of intermittent phenotype MSUD in some patients can be due to mutations in the E2 gene, giving rise to a low but significant residual activity of the BCKDH complex.     

Molecular phylogeny of Fv1

Alleles at the Fv1 gene of inbred mice confer resistance to infection and spread of vertically or horizontally transmitted murine leukemia viruses (MuLV). The nucleotide sequence of Fv1 bears similarity to the gag of a human endogenous retrovirus, HERV-L, but is more closely related to the gag-coding sequence of a newly described class of HERV-L-related mouse endogenous retroviruses designated MuERV-L. Both observations suggest an origin of Fv1 from endogenous gag sequences. The molecular definition of Fv1 provided an opportunity to determine the phylogeny of the gene among wild mice and its relation to MuERV-L. PCR primers, chosen to include most of the coding region of Fv1 for both the n and b alleles, were used to amplify sequences from animals of the genus Mus, which were then sequenced. Closely related products were obtained from almost all animals examined that evolved after the separation from Rattus, in which the homologous gene was shown to be absent. A phylogenetic tree generated with Fv1 sequence data differs noticeably from that developed with sequence data from other genes. In addition, non-synonymous changes were found to be present twice as frequently as synonymous changes, a fact that departs from the standard behavior of a structural gene. These observations suggest that the Fv1 gene may have been subjected to possible horizontal transfers as well as to positive Darwinian selection.