Adaptive evolution of genes and gene families

Multigenerational Effects of Graphene Oxide Nanoparticles on Acheta domesticus DNA Stability

The use of nanoparticles like graphene oxide (GO) in nanocomposite industries is growing very fast. There is a strong concern that GO can enter the environment and become nanopollutatnt. Environmental pollutants' exposure usually relates to low concentrations but may last for a long time and impact following generations. Attention should be paid to the effects of nanoparticles, especially on the DNA stability passed on to the offspring. We investigated the multigenerational effects on two strains (wild and long-lived) of house cricket intoxicated with low GO concentrations over five generations, followed by one recovery generation. Our investigation focused on oxidative stress parameters, specifically AP sites (apurinic/apyrimidinic sites) and 8-OHdG (8-hydroxy-2'-deoxyguanosine), and examined the global DNA methylation pattern. Five intoxicated generations were able to overcome the oxidative stress, showing that relatively low doses of GO have a moderate effect on the house cricket (8-OHdG and AP sites). The last recovery generation that experienced a transition from contaminated to uncontaminated food presented greater DNA damage. The pattern of DNA methylation was comparable in every generation, suggesting that other epigenetic mechanisms might be involved.

A first annotated genome sequence for Haliotis midae with genomic insights into abalone evolution and traits of economic importance

Haliotis midae or "perlemoen" is one of five abalone species endemic to South Africa, and being palatable, the only commercially important abalone species with a high international demand. The higher demand for this abalone species has resulted in the decrease of natural stocks due to overexploitation by capture fisheries and poaching. Facilitating aquaculture production of H. midae should assist in minimising the pressure on the wild populations. Here, the draft genome of H. midae has been sequenced, assembled, and annotated. The draft assembly resulted in a total length of 1.5 Gb, contig N50 of 0.238 Mb, scaffold N50 of 0. 238 Mb and GC level of 40%. Gene annotation, combining ab initio and evidence-based pipelines identified 52,280 genes with protein coding potential. The genes identified were used to predict orthologous genes shared among the four other abalone species (H. laevigata, H. rubra, H. discus hannai and H. rufescens) and 4702 orthologous genes were shared across the five species. Among the orthologous genes in abalones, single copy genes were further analysed for signatures of selection and several molecular regulatory proteins involved in developmental functions were found to be under positive selection in specific abalone lineages. Furthermore, whole genome SNP-based phylogenomic assessment was performed to confirm the evolutionary relationship among the considered abalone species with draft genomes, reaffirming that H. midae is closely related to the Australian Greenlip (H. laevigata) and Blacklip (H. rubra). The study assists in the understanding of genes related to various biological systems underscoring the evolution and development of abalones, with potential applications for genetic improvement of commercial stocks.

A computational investigation of cis-gene regulation in evolution

In biological processes involving gene networks, genes regulate other genes that determine the phenotypic traits. Gene regulation plays an important role in evolutionary dynamics. In a genetic algorithm, a trans-gene regulatory mechanism was shown to speed up adaptation and evolution. Here, we examine the effect of cis-gene regulation on an adaptive system. The model is haploid. A chromosome is partitioned into regulatory loci and structural loci. The regulatory genes regulate the expression and functioning of structural genes via the cis-elements in a probabilistic manner. In the simulation, the change in the allele frequency, the mean population fitness and the efficiency of phenotypic selection are monitored. Cis-gene regulation increases adaption and accelerates the evolutionary process in comparison with the case involving absence of gene regulation. Some special features of the simulation results are as follows. A low ratio of regulatory loci and structural loci gives higher adaptation for fixed total number of loci. Plasticity is advantageous beyond a threshold value. Adaptation is better for large number of total loci when the ratio of regulatory loci to structural loci is one. However, it reaches a saturation beyond which the increase in the total loci is not advantageous. Efficiency of the phenotypic selection is higher for larger value of the initial plasticity.

Adaptive evolution driving the young duplications in six Rosaceae species

Background

In plant genomes, high proportions of duplicate copies reveals that gene duplications play an important role in the evolutionary processes of plant species. A series of gene families under positive selection after recent duplication events in plant genomes indicated the evolution of duplicates driven by adaptive evolution. However, the genome-wide evolutionary features of young duplicate genes among closely related species are rarely reported.

Results

In this study, we conducted a systematic survey of young duplicate genes at genome-wide levels among six Rosaceae species, whose whole-genome sequencing data were successively released in recent years. A total of 35,936 gene families were detected among the six species, in which 60.25% were generated by young duplications. The 21,650 young duplicate gene families could be divided into two expansion types based on their duplication patterns, species-specific and lineage-specific expansions. Our results showed the species-specific expansions advantaging over the lineage-specific expansions. In the two types of expansions, high-frequency duplicate domains exhibited functional preference in response to environmental stresses.

Conclusions

The functional preference of the young duplicate genes in both the expansion types showed that they were inclined to respond to abiotic or biotic stimuli. Moreover, young duplicate genes under positive selection in both species-specific and lineage-specific expansions suggested that they were generated to adapt to the environmental factors in Rosaceae species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07422-7.

LMAP: Lightweight Multigene Analyses in PAML

Background

Uncovering how phenotypic diversity arises and is maintained in nature has long been a major interest of evolutionary biologists. Recent advances in genome sequencing technologies have remarkably increased the efficiency to pinpoint genes involved in the adaptive evolution of phenotypes. Reliability of such findings is most often examined with statistical and computational methods using Maximum Likelihood codon-based models (i.e., site, branch, branch-site and clade models), such as those available in codeml from the Phylogenetic Analysis by Maximum Likelihood (PAML) package. While these models represent a well-defined workflow for documenting adaptive evolution, in practice they can be challenging for researchers having a vast amount of data, as multiple types of relevant codon-based datasets are generated, making the overall process hard and tedious to handle, error-prone and time-consuming.

Results

We introduce LMAP (Lightweight Multigene Analyses in PAML), a user-friendly command-line and interactive package, designed to handle the codeml workflow, namely: directory organization, execution, results gathering and organization for Likelihood Ratio Test estimations with minimal manual user intervention. LMAP was developed for the workstation multi-core environment and provides a unique advantage for processing one, or more, if not all codeml codon-based models for multiple datasets at a time. Our software, proved efficiency throughout the codeml workflow, including, but not limited, to simultaneously handling more than 20 datasets.

Conclusions

We have developed a simple and versatile LMAP package, with outstanding performance, enabling researchers to analyze multiple different codon-based datasets in a high-throughput fashion. At minimum, two file types are required within a single input directory: one for the multiple sequence alignment and another for the phylogenetic tree. To our knowledge, no other software combines all codeml codon substitution models of adaptive evolution. LMAP has been developed as an open-source package, allowing its integration into more complex open-source bioinformatics pipelines. LMAP package is released under GPLv3 license and is freely available at http://lmapaml.sourceforge.net/.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1204-5) contains supplementary material, which is available to authorized users.

Differential Evolutionary Selection and Natural Evolvability Observed in ALT Proteins of Human Filarial Parasites

The abundant larval transcript (ALT-2) protein is present in all members of the Filarioidea, and has been reported as a potential candidate antigen for a subunit vaccine against lymphatic filariasis. To assess the potential for vaccine escape or heterologous protection, we examined the evolutionary selection acting on ALT-2. The ratios of nonsynonymous (K(a)) to synonymous (K(s)) mutation frequencies (ω) were calculated for the alt-2 genes of the lymphatic filariasis agents Brugia malayi and Wuchereria bancrofti and the agents of river blindness and African eyeworm disease Onchocerca volvulus and Loa loa. Two distinct Bayesian models of sequence evolution showed that ALT-2 of W. bancrofti and L. loa were under significant (P<0.05; P < 0.001) diversifying selection, while ALT-2 of B. malayi and O. volvulus were under neutral to stabilizing selection. Diversifying selection as measured by ω values was notably strongest on the region of ALT-2 encoding the signal peptide of L. loa and was elevated in the variable acidic domain of L. loa and W. bancrofti. Phylogenetic analysis indicated that the ALT-2 consensus sequences formed three clades: the first consisting of B. malayi, the second consisting of W. bancrofti, and the third containing both O. volvulus and L. loa. ALT-2 selection was therefore not predictable by phylogeny or pathology, as the two species parasitizing the eye were selected differently, as were the two species parasitizing the lymphatic system. The most immunogenic regions of L. loa and W. bancrofti ALT-2 sequence as modeled by antigenicity prediction analysis did not correspond with elevated levels of diversifying selection, and were not selected differently than predicted antigenic epitopes in B. malayi and O. volvulus. Measurements of ALT-2 evolvability made by χ2 analysis between alleles that were stable (O. volvulus and B. malayi) and those that were under diversifying selection (W. bancrofti and L. loa) indicated significant (P<0.01) deviations from a normal distribution for both W. bancrofti and L. loa. The relationship between evolvability and selection in L. loa followed a second order polynomial distribution (R2 = 0.89), indicating that the two factors relate to one another in accordance with an additional unknown factor. Taken together, these findings indicate discrete evolutionary drivers acting on ALT-2 of the four organisms examined, and the described variation has implications for design of novel vaccines and diagnostic reagents. Additionally, this represents the first mathematical description of evolvability in a naturally occurring setting.

Genome-wide prediction of the polymorphic Ser gene family in Tetrahymena thermophila based on motif analysis

Even though antigenic variation is employed among parasitic protozoa for host immune evasion, Tetrahymena thermophila, a free-living ciliate, can also change its surface protein antigens. These cysteine-rich glycosylphosphatidylinositol (GPI)-linked surface proteins are encoded by a family of polymorphic Ser genes. Despite the availability of T. thermophila genome, a comprehensive analysis of the Ser family is limited by its high degree of polymorphism. In order to overcome this problem, a new approach was adopted by searching for Ser candidates with common motif sequences, namely length-specific repetitive cysteine pattern and GPI anchor site. The candidate genes were phylogenetically compared with the previously identified Ser genes and classified into subtypes. Ser candidates were often found to be located as tandem arrays of the same subtypes on several chromosomal scaffolds. Certain Ser candidates located in the same chromosomal arrays were transcriptionally expressed at specific T. thermophila developmental stages. These Ser candidates selected by the motif analysis approach can form the foundation for a systematic identification of the entire Ser gene family, which will contribute to the understanding of their function and the basis of T. thermophila antigenic variation.

Positive selection drives accelerated evolution of mosquito salivary genes associated with blood-feeding

The saliva of bloodsucking animals contains dozens to hundreds of proteins that counteract their hosts' haemostasis, inflammation and immunity. It was previously observed that salivary proteins involved in haematophagy are much more divergent in their primary sequence than those of housekeeping function, when comparisons were made between closely related organisms. While this pattern of evolution could result from relaxed selection or drift, it could alternatively be the result of positive selection driven by the intense pressure of the host immune system. We investigated the polymorphism of five different genes associated with blood-feeding in the mosquito Anopheles gambiae and obtained evidence in four genes for sites with signatures of positive selection. These results add salivary gland genes from bloodsucking arthropods to the small list of genes driven by positive selection.

Evolutionary perspectives into placental biology and disease

In all mammals including humans, development takes place within the protective environment of the maternal womb. Throughout gestation, nutrients and waste products are continuously exchanged between mother and fetus through the placenta. Despite the clear importance of the placenta to successful pregnancy and the health of both mother and offspring, relatively little is understood about the biology of the placenta and its role in pregnancy-related diseases. Given that pre- and peri-natal diseases involving the placenta affect millions of women and their newborns worldwide, there is an urgent need to understand placenta biology and development. Here, we suggest that the placenta is an organ under unique selective pressures that have driven its rapid diversification throughout mammalian evolution. The high divergence of the placenta complicates the use of non-human animal models and necessitates an evolutionary perspective when studying its biology and role in disease. We suggest that diversifying evolution of the placenta is primarily driven by intraspecies evolutionary conflict between mother and fetus, and that many pregnancy diseases are a consequence of this evolutionary force. Understanding how maternal-fetal conflict shapes both basic placental and reproductive biology - in all species - will provide key insights into diseases of pregnancy.

Theory of the origin, evolution, and nature of life

Life is an inordinately complex unsolved puzzle. Despite significant theoretical progress, experimental anomalies, paradoxes, and enigmas have revealed paradigmatic limitations. Thus, the advancement of scientific understanding requires new models that resolve fundamental problems. Here, I present a theoretical framework that economically fits evidence accumulated from examinations of life. This theory is based upon a straightforward and non-mathematical core model and proposes unique yet empirically consistent explanations for major phenomena including, but not limited to, quantum gravity, phase transitions of water, why living systems are predominantly CHNOPS (carbon, hydrogen, nitrogen, oxygen, phosphorus, and sulfur), homochirality of sugars and amino acids, homeoviscous adaptation, triplet code, and DNA mutations. The theoretical framework unifies the macrocosmic and microcosmic realms, validates predicted laws of nature, and solves the puzzle of the origin and evolution of cellular life in the universe.

Comparative analysis of evolutionary dynamics of genes encoding leucine-rich repeat receptor-like kinase between rice and Arabidopsis

The leucine-rich repeat (LRR) receptor kinase (RLK) proteins constitute a large superfamily in the plant genome, and carry out key functions in a variety of biological pathways. In an effort to determine the evolutionary fate of members of a large gene family such as plant LRR RLK proteins we conducted in silico analysis using complete genome sequencing datasets, genome-wide transcriptome databases, and bioinformatics tools. A total of 292 and 165 LRR RLK genes were retrieved from the rice and Arabidopsis genomes, respectively, formed by diverse duplication events for gene expansion. The phylogenic analyses of the LRR RLK genes suggested combinations of LRR domains and RLK domains in the ancient plant genome prior to the divergence of rice and Arabidopsis, followed by massive independent expansions during speciation. The somewhat high frequencies (50-73%) of expressional divergence of members of duplicate gene pairs formed by whole/segmental genome duplication (W/SGD) and tandem duplication (TD) events of Arabidopsis and TD events of rice support the idea of their functional diversity for gene retention. By contrast, a relatively low degree (at least 20%) of members of rice LRR RLK gene pairs formed by W/SGD appear to be divergent in expression following the duplication event. At least 7 pairs of co-expressed gene clusters, including each of the tentative orthologous LRR RLK genes between rice and Arabidopsis, were enriched to an orthologous set between members of each of the pairs as compared to those of the random pairs, suggesting some degree of functional conservation of individual genes. These results may shed some light on the crucial functions of the plant LRR RLK genes with regard to a variety of biological processes.

A candidate subspecies discrimination system involving a vomeronasal receptor gene with different alleles fixed in M. m. domesticus and M. m. musculus

Assortative mating, a potentially efficient prezygotic reproductive barrier, may prevent loss of genetic potential by avoiding the production of unfit hybrids (i.e., because of hybrid infertility or hybrid breakdown) that occur at regions of secondary contact between incipient species. In the case of the mouse hybrid zone, where two subspecies of Mus musculus (M. m. domesticus and M. m. musculus) meet and exchange genes to a limited extent, assortative mating requires a means of subspecies recognition. We based the work reported here on the hypothesis that, if there is a pheromone sufficiently diverged between M. m. domesticus and M. m. musculus to mediate subspecies recognition, then that process must also require a specific receptor(s), also sufficiently diverged between the subspecies, to receive the signal and elicit an assortative mating response. We studied the mouse V1R genes, which encode a large family of receptors in the vomeronasal organ (VNO), by screening Perlegen SNP data and identified one, Vmn1r67, with 24 fixed SNP differences most of which (15/24) are nonsynonymous nucleotide substitutions between M. m. domesticus and M. m. musculus. We observed substantial linkage disequilibrium (LD) between Vmn1r67 and Abpa27, a mouse salivary androgen-binding protein gene that encodes a proteinaceous pheromone (ABP) capable of mediating assortative mating, perhaps in conjunction with its bound small lipophilic ligand. The LD we observed is likely a case of association rather than residual physical linkage from a very recent selective sweep, because an intervening gene, Vmn1r71, shows significant intra(sub)specific polymorphism but no inter(sub)specific divergence in its nucleotide sequence. We discuss alternative explanations of these observations, for example that Abpa27 and Vmn1r67 are coevolving as signal and receptor to reinforce subspecies hybridization barriers or that the unusually divergent Vmn1r67 allele was not a product of fast positive selection, but was derived from an introgressed allele, possibly from Mus spretus.

Molecular evolutionary analysis of seminal receptacle sperm storage organ genes of Drosophila melanogaster

Sperm storage organs are common and broadly distributed among animal taxa. However, little is known about how these organs function at the molecular level. Additionally, there is a paucity of knowledge about the evolution of genes expressed in these organs. This investigation is an evolutionary expressed sequence tag (EST) study of genes expressed in the seminal receptacle, one of the sperm storage organs in Drosophila. The incidence of positive selection is higher for the seminal receptacle genes than Drosophila reproductive genes as a whole, but lower than genes associated with the spermatheca, a second type of Drosophila sperm storage organ. By identifying overrepresented classes of proteins and classes for which sperm storage function is suggested by the nature of the proteins, candidate genes were discovered. These candidates belong to protein classes such as muscle contraction, odorant binding and odorant receptor, protease inhibitor and immunity.

Evolutionary fingerprinting of genes

Over time, natural selection molds every gene into a unique mosaic of sites evolving rapidly or resisting change-an "evolutionary fingerprint" of the gene. Aspects of this evolutionary fingerprint, such as the site-specific ratio of nonsynonymous to synonymous substitution rates (dN/dS), are commonly used to identify genetic features of potential biological interest; however, no framework exists for comparing evolutionary fingerprints between genes. We hypothesize that protein-coding genes with similar protein structure and/or function tend to have similar evolutionary fingerprints and that comparing evolutionary fingerprints can be useful for discovering similarities between genes in a way that is analogous to, but independent of, discovery of similarity via sequence-based comparison tools such as Blast. To test this hypothesis, we develop a novel model of coding sequence evolution that uses a general bivariate discrete parameterization of the evolutionary rates. We show that this approach provides a better fit to the data using a smaller number of parameters than existing models. Next, we use the model to represent evolutionary fingerprints as probability distributions and present a methodology for comparing these distributions in a way that is robust against variations in data set size and divergence. Finally, using sequences of three rapidly evolving RNA viruses (HIV-1, hepatitis C virus, and influenza A virus), we demonstrate that genes within the same functional group tend to have similar evolutionary fingerprints. Our framework provides a sound statistical foundation for efficient inference and comparison of evolutionary rate patterns in arbitrary collections of gene alignments, clustering homologous and nonhomologous genes, and investigation of biological and functional correlates of evolutionary rates.

Evolutionary aspects of functional and pseudogene members of the phytochrome gene family in Scots pine

According to the neutral theory of evolution, mutation and genetic drift are the only forces that shape unconstrained, neutral, gene evolution. Thus, pseudogenes (which often evolve neutrally) provide opportunities to obtain direct estimates of mutation rates that are not biased by selection, and gene families comprising functional and pseudogene members provide useful material for both estimating neutral mutation rates and identifying sites that appear to be under positive or negative selection pressures. Conifers could be very useful for such analyses since they have large and complex genomes. There is evidence that pseudogenes make significant contributions to the size and complexity of gene families in pines, although few studies have examined the composition and evolution of gene families in conifers. In this work, I examine the complexity and rates of mutation of the phytochrome gene family in Pinus sylvestris and show that it includes not only functional genes but also pseudogenes. As expected, the functional PHYO does not appear to have evolved neutrally, while phytochrome pseudogenes show signs of unconstrained evolution.

Mammalian evolution and biomedicine: new views from phylogeny

Recent progress resolving the phylogenetic relationships of the major lineages of mammals has had a broad impact in evolutionary biology, comparative genomics and the biomedical sciences. Novel insights into the timing and historical biogeography of early mammalian diversification have resulted from a new molecular tree for placental mammals coupled with dating approaches that relax the assumption of the molecular clock. We highlight the numerous applications to come from a well-resolved phylogeny and genomic prospecting in multiple lineages of mammals, from identifying regulatory elements in mammalian genomes to assessing the functional consequences of mutations in human disease loci and those driving adaptive evolution.

Diversification and specialization of HIV protease function during in vitro evolution

Our goal is to understand how enzymes adapt to utilize novel substrates. We and others have shown that directed evolution tends to generate enzyme variants with broadened substrate specificity. Broad-specificity enzymes are generally deleterious to living cells, so this observed trend might be an artifact of the most commonly employed high throughput screens. Here, we demonstrate a more natural and effective screening strategy for directed evolution. The gene encoding model enzyme HIV protease was randomly mutated, and the resulting library was expressed in Escherichia coli cells to eliminate cytotoxic broad-specificity variants. The surviving variants were screened for clones with activity against a reporter enzyme. The wild-type human immunodeficiency virus type I protease (HIV PR) is cytotoxic and exhibits no detectable activity in reactions with beta-galactosidase (BGAL). In contrast, the selected variants were nontoxic and exhibited greater activity and specificity against BGAL than did the wild-type HIV PR in reactions with any substrate. A single round of whole gene random mutagenesis and conventional high-throughput screening does not usually effect complete inversions of substrate specificity. This suggests that a combination of positive and purifying selection engenders more rapid adaptation than positive selection alone.

Highly polymorphic vitelline-coat protein HaVC80 from the ascidian, Halocynthia aurantium: structural analysis and involvement in self/nonself recognition during fertilization

Ascidians release sperm and eggs simultaneously, but self-fertilization is effectively blocked by unknown mechanisms. We previously reported that a 70-kDa sperm receptor HrVC70 on the egg vitelline coat (VC) consisting of 12 EGF-like repeats is a candidate self/nonself recognition molecule during fertilization of the ascidian, Halocynthia roretzi. Here, we report that Halocynthia aurantium also utilizes a homolog (HaVC80) of HrVC70 as an allorecognizable sperm receptor. HaVC80 is attached to the VC during the acquisition of self-sterility and is detached from the VC by acid treatment, allowing self-fertilization. A cDNA clone of the HaVC80 precursor, HaVC130, consists of 3726 nucleotides and encodes an open reading frame of 1208 amino acids. The structure of HaVC130 is very similar to the HrVC70 precursor HrVC120, but the number of EGF-like repeats of HaVC130/VC80 is one repeat larger than that of HrVC120/VC70. There are several amino acid substitutions between different individuals, and two alleles of the HaVC80 sequence were detected in each individual. Genomic DNA sequence analysis reveals that each EGF-like domain corresponds to a specific exon, and HaVC130 may have been evolutionarily generated from HrVC120 by duplication of the 8th EGF-like repeat. The data support the hypothesis that HaVC80 is a highly polymorphic protein responsible for self-sterility in H. aurantium.

Comparative primate genomics: the year of the chimpanzee

This is the year of the chimpanzee genome. Chimpanzee chromosome 22 has been sequenced and soon will be followed by the whole genome, and thousands of chimpanzee cDNA sequences are available for comparative analysis. Not only does this genomic information allow us to identify human-specific changes in particular genes that are potentially under selection, but also to understand molecular evolutionary dynamics characterizing the two most closely related mammalian genomes sequenced so far. Studies comparing gene expression in chimpanzees and other closely related primates reveal significant species differences in brain, liver and fibroblasts. New empirical data, in combination with models of speciation, are giving insight into how humans and chimpanzees speciated.

The genetic architecture of selection at the human dopamine receptor D4 (DRD4) gene locus

Associations of the seven-repeat (7R) allele of the human dopamine receptor D4 (DRD4) gene with both the personality trait of novelty seeking and attention deficit/hyperactivity disorder have been reported. Recently, on the basis of the unusual DNA sequence organization of the DRD4 7R 48-bp tandem repeat (VNTR), we proposed that the 7R allele originated as a rare mutational event that increased to high frequency by positive selection. We now have resequenced the entire DRD4 locus from 103 individuals homozygous for 2R, 4R, or 7R variants of the VNTR, a method developed to directly estimate haplotype diversity. DNA from individuals of African, European, Asian, North and South American, and Pacific Island ancestry were used. 4R/4R homozygotes exhibit little linkage disequilibrium (LD) over the region examined, with more polymorphisms observed in DNA samples from African individuals. In contrast, the evidence for strong LD surrounding the 7R allele is dramatic, with all 7R/7R individuals (including those from Africa) exhibiting the same alleles at most polymorphic sites. By intra-allelic comparison at 18 high-heterozygosity sites spanning the locus, we estimate that the 7R allele arose prior to the upper Paleolithic era (approximately 40000-50000 years ago). Further, the pattern of recombination at these polymorphic sites is the pattern expected for selection acting at the 7R VNTR itself, rather than at an adjacent site. We propose a model for selection at the DRD4 locus consistent with these observed LD patterns and with the known biochemical and physiological differences between receptor variants.

A perforin-like protein from a marine mollusk

Abalone (gastropod mollusks) express a protein, abMpeg1, which is a homolog of two mammalian proteins that share homology with mammalian perforin, a cytolytic and immune-regulatory protein of lymphocytes. One of the mammalian proteins, Mpeg1, is expressed in mature macrophage and prion-infected mouse brains, while the other, Epcs50, is expressed in ectoplacental cone cells of the invading placenta. Although the functions of these three proteins remain unknown, their structural similarity to mammalian perforin suggests that they may be involved in cell killing, the inflammatory response or tissue invasion. Consistent with these proposed functions, the Mpeg1 gene family shows the signature of positive Darwinian selection (adaptive evolution). The perforin-homology domain of abMpeg1 contains the cytolytic "helix-turn-helix" domain of perforin, supporting the idea that abMpeg1 is a cytolytic protein of the abalone innate immune system. The alpha-helices of abMpeg1 are amphipathic as are those of perforin. The conservation among abMpeg1, mammalian Mpeg1, and Epcs50 shows that Mpeg1 proteins represent a novel, ancient protein family of probable immunological function.