The mate recognition protein gene mediates reproductive isolation and speciation in the Brachionus plicatilis cryptic species complex

ABSTRACT: BACKGROUND: Chemically mediated prezygotic barriers to reproduction likely play an important role in speciation. In facultatively sexual monogonont rotifers from the Brachionus plicatilis cryptic species complex, mate recognition of females by males is mediated by the Mate Recognition Protein (MRP), a globular glycoprotein on the surface of females, encoded by the mmr-b gene family. In this study, we sequenced mmr-b copies from 27 isolates representing 11 phylotypes of the B. plicatilis species complex, examined the mode of evolution and selection of mmr-b, and determined the relationship between mmr-b genetic distance and mate recognition among isolates. RESULTS: Isolates of the B. plicatilis species complex have 1-4 copies of mmr-b, each composed of 2-9 nearly identical tandem repeats. The repeats within a gene copy are generally more similar than are gene copies among phylotypes, suggesting concerted evolution. Compared to housekeeping genes from the same isolates, mmr-b has accumulated only half as many synonymous differences but twice as many non-synonymous differences. Most of the amino acid differences between repeats appear to occur on the outer face of the protein, and these often result in changes in predicted patterns of phosphorylation. However, we found no evidence of positive selection driving these differences. Isolates with the most divergent copies were unable to mate with other isolates and rarely self-crossed. Overall the degree of mate recognition was significantly correlated with the genetic distance of mmr-b. CONCLUSIONS: Discrimination of compatible mates in the B. plicatilis species complex is determined by proteins encoded by closely related copies of a single gene, mmr-b. While concerted evolution of the tandem repeats in mmr-b may function to maintain identity, it can also lead to the rapid spread of a mutation through all copies in the genome and thus to reproductive isolation. The mmr-b gene is evolving rapidly, and novel alleles may be maintained and increase in frequency via asexual reproduction. Our analyses indicate that mate recognition, controlled by MMR-B, may drive reproductive isolation and allow saltational sympatric speciation within the B. plicatilis cryptic species complex, and that this process may be largely neutral.     

Special evolutionary properties of genes encoding a protein with a simple amino acid repeat

We have examined the evolution of a gene, SM50, encoding a component of the spicule matrix, which plays an integral role in the formation of the echinoderm skeleton. This gene was originally characterized in Strongylocentrotus purpuratus and encodes an imperfect tandem repeat of six or seven amino acids. We have analyzed the sequence of this repeat in a number of sea urchin species and have determined that the repeat regions have undergone concerted evolution. There are differences in the repeat region between species, but the overall repeat structure is conserved, suggesting the repeat forms a structural domain important in biomineralization. The inherent conserved amino acid repeat structure promotes concerted evolution due to the high probability of misreplication and unequal crossing-over in the repeated segment of the gene. While there are constraints on the amino acids allowed in the repeat region, there are also variations, so that the sequences observed illustrate the balance between amino acid substitutions and concerted evolution. We have evidence that substitutions can alter the mechanisms of unequal crossing-over, altering the way concerted evolution occurs. The way in which concerted evolution occurred appears to be determined by the degree of sequence similarity between the repeats in a given gene, which influences how unequal crossing over may occur. We have mapped the differences in repeat regions on existing phylogenetic trees and indicate where concerted evolution has taken place. We also confirm an earlier report that Hemicentrotus pulcherrimus fits into the Strongylocentrotus genus and examine the evolution of the H. pulcherrimus SM50 repeat relative to other members of this genus. Received: 31 October 2000 / Accepted: 20 March 2001     

Freshwater sponge silicateins: Comparison of gene sequences and exon-intron structure

Silicateins found in spicules of siliceous sponges are proteins that take part in biogenic silica precipitation and determine the morphological features of spicules. The exon-intron structure of the genes encoding four silicatein-α isoforms (−α1, −α2, −α3, and −α4) from an endemic Baikalian sponge Lubomirskia baicalensis was studied. For eight sponge species, including both cosmopolitan (Spongilla lacustris, Ephydatia muelleri, E. fluviatilis) and endemic Baikalian (L. baicalensis, L. incrustans, Baikalospongia intermedia, B. fungiformis, Sw. papyracea) species, seventeen partial sequences of different silicatein isoform genes were determined. It was shown that cosmopolitan and endemic Baikalian sponges differ from each other in gene structure, in particular, in intron length. Among Baikalian sponges, silicatein-α1 genes had the highest variation of intron length, and silicatein-α4 genes were the most conservative. A phylogenetic analysis based on amino acid sequences of different silicatein isoforms identified four distinct clusters within the freshwater sponge clade. An analysis based on exon-intron gene sequences enables discrimination between different sponge species within the clusters.     

Axenic culture of <Emphasis Type="Italic">Brachionus plicatilis</Emphasis> using antibiotics

The rotifer Brachionus plicatilis culture is composed of complex microcosms including bacteria, protozoans, algae, and fungi. Previous studies reported methods to establish axenic rotifer cultures, but further refinement of these techniques is needed, for molecular biological research which requires pure culture to isolate nucleic acids from rotifers only. In order to render rotifer culture axenic, we tested five antibiotics: ampicillin (Amp), chloramphenicol (Cp), kanamycin (Km), nalidixic acid (Na), and streptomycin (Sm) at 30–100 μg/ml. Except for Cp, which reduces rotifer reproduction, all other antibiotics at the tested concentrations did not affect rotifer reproduction or show any toxic effects. A rotifer disinfection method was finally established by treating the resting eggs with 0.25% (w/v) sodium hypochlorite (NaOCl) for 3 min, washing with sterilized sea water, and then exposing the neonates to an Amp, Km, Na, and Sm mixture. Using four nutrient media, we confirmed that this protocol renders the rotifer culture bacterial and fungus free. The axenic rotifer culture generated here is useful not only for genetic analysis of Brachionus plicatilis, but for studying the rotifer life cycle without bacterial influence.     

A novel growth-promoting protein in the conditioned media from the rotifer <Emphasis Type="Italic">Brachionus plicatilis</Emphasis> at an early exponential growth phase

We confirmed the existence of growth-promoting substances in the conditioned media (CM) from the rotifer Brachionus plicatilis at an early exponential growth phase and isolated a novel protein with a growth-promoting activity from the crude extract (CE) of rotifer cells. CM was prepared from the culture media where rotifers had been cultured at an early exponential growth phase and filtered through a 0.22-μm filter membrane. The growth-promoting activity was determined using rotifers in CM for 5 days. As a result, the increase of rotifers added with CM was significantly higher than that of the control in artificial seawater (P < 0.001). Moreover, the growth-promoting activity of CM was dose-dependent and inactivated by heat treatment at 80°C for 60 min. Meanwhile, CM filtered through a <10 kDa ultrafiltration membrane showed a low activity, whereas proteinase K treatment resulted in a complete inactivation. These results suggest that the rotifer secrets growth-promoting proteins into CM. CE also contained a protein with the activity and properties similar to those found in CM. Then, CE was subjected to purification of a growth-promoting protein for convenience using various types of chromatography after fractionation with 30–80% saturated ammonium sulfate. Subsequently, a protein with an approximate molecular weight of 25000 was isolated, and its N-terminal amino acid sequence was determined to be PAVVDFTAVWFGPLQMIKP. An orthologue was found in the EST database of B. plicatilis, the full sequence of which showed about 50% identity to the corresponding regions of thioredoxins from other organisms.     

Widespread Gene Conversion of Alpha-2-Fucosyltransferase Genes in Mammals

The alpha-2-fucosyltransferases (α2FTs) are enzymes involved in the biosynthesis of α2fucosylated glycan structures. In mammalian genomes, there are three α2FT genes located in tandem—FUT1, FUT2, and Sec1—each contained within a single exon. It has been suggested that these genes originated from two successive duplications, with FUT1 being generated first and FUT2 and Sec1 second. Despite gene conversion being considered the main mechanism of concerted evolution in gene families, previous studies of primates α2FTs failed to detect it, although the occurrence of gene conversion between FUT2 and Sec1 was recently reported in a human allele. The primary aim of our work was to initiate a broader study on the molecular evolution of mammalian α2FTs. Sequence comparison leads us to confirm that the three genes appeared by two rounds of duplication. In addition, we were able to detect multiple gene-conversion events at the base of primates and within several nonprimate species involving FUT2 and Sec1. Gene conversion involving FUT1 and either FUT2 or Sec1 was also detected in rabbit. The extent of gene conversion between the α2FTs genes appears to be species-specific, possibly related to functional differentiation of these genes. With the exception of rabbits, gene conversion was not observed in the region coding the C-terminal part of the catalytic domain. In this region, the number of amino acids that are identical between FUT1 and FUT2, but different in Sec1, is higher than in other parts of the protein. The biologic meaning of this observation may be related to functional constraints. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Development of genomic resources for the phylogenetic analysis of the <Emphasis Type="Italic">Brachionus plicatilis</Emphasis> species complex (Rotifera: Monogononta)

There has been a substantial leap forward in the quantity and quality of genomic resources available for research on rotifers in recent years. We used available genomic and bioinformatics resources to identify variable regions of the genome to design PCR and sequencing primers for the Brachionus plicatilis complex. We then tested their suitability for the study of systematics of this group. Eight amplified successfully for members of the B. plicatilis complex. We sequenced the amplified products, constructed a concatenated alignment (5,511 bp) and carried out phylogenetic analyses. The resulting tree, based on mitochondrial and nuclear genes for 11 clones of six species, was well supported even at relatively deep nodes, contrasting with results of previous studies using only mitochondrial genes, which provide little phylogenetic information for the deepest nodes. The same procedure could be used to design primers for more conserved regions to be used in a wider range of rotifer groups.     

Molecular evolution of the membrane associated progesterone receptor in the <Emphasis Type="Italic">Brachionus plicatilis</Emphasis> (Rotifera,Monogononta) species complex

Many studies have investigated physiological roles of the membrane associated progesterone receptor (MAPR), but little is known of its evolution. Marked variations in response to exogenous progesterone have been reported for four brachionid rotifer species, suggesting differences in progesterone signaling and reception. Here we report sequence variation for the MAPR gene in the Brachionus plicatilis species complex. Phylogenetic analysis of this receptor is compared with relatedness based on cytochrome c oxidase subunit 1 sequences. Nonsynonymous to synonymous site substitution rate ratios, amino acid divergence, and variations in predicted phosphorylation sites are examined to assess evolution of the MAPR among brachionid clades.     

Multisequence comparisons in protein coding genes

A very powerful method for detecting functional constraints operative in biological macromolecules is presented. This method entails performing a base permanence analysis of protein coding genes at each codon position simultaneously in different species. It calculates the degree of permanence of subregions of the gene by dividing it into segments,c codons long, counting how many sites remain unchanged in each segment among all species compared. By comparing the base permanence among several sequences with the expectations based on a stochastic evolutionary process, gene regions showing different degrees of conservation can be selected. This means that wherever the permanence deviates significantly from the expected value generated by the simulation, the corresponding regions are considered “constrained” or “hypervariable”. The constrained regions are of two types: α and β. The α regions result from constraints at the amino acid level, whereas the β regions are those probably involved in “control” processing. The method has been applied to mitochondrial genes coding for subunit 6 of the ATPase and subunit 1 of the cytochrome oxidase in four mammalian species: human, rat, mouse, and cow. In the two mitochondrial genes a few regions that are highly conserved in all codon positions have been identified. Among these regions a sequence, common to both genes, that is complementary to a strongly conserved region of 12S rRNA has been found. This method can also be of great help in studying molecular evolution mechanisms.     

A divergent non-classical class I gene conserved in salmonids

 Complementary DNA for two class I genes of the rainbow trout, Oncorhynchus mykiss, were characterized. MhcOnmy-UBA*01 is similar to Onmy-UA-C32 and the classical major histocompatibility complex class I genes of other fish species, whereas Onmy-UAA*01 is divergent from all class I genes so far characterized. Onmy-UAA*01 is expressed at lower levels than Onmy-UBA*01. Although Onmy-UAA*01 exhibits restriction fragment length polymorphism on Southern blotting, the encoded protein is highly conserved. Two allotypes, which differ only by substitution at amino acid position 223 of the α₃ domain, have been defined. Onmy-UAA*01 has an exon-intron organization like other class I genes and contains a Tc1-like transposon element in intron III. Orthologues of Onmy-UAA*01 have been characterized in four other species of salmonid. Between four species of Oncorhynchus, UAA*01 proteins differ by only 2–6 amino acids, whereas comparison of Oncorhynchus with Salmo trutta (brown trout) reveals 14–16 amino acid differences. The Onmy-UAA*01 gene has properties indicative of a particularly divergent non-classical class I gene. Received: 22 September 1998 / Revised: 24 November 1998     

Many gene-regulatory proteins appear to have a similar α-helical fold that binds DNA and evolved from a common precursor

Summary Amino acid and DNA sequence comparisons suggest that many sequence-specific DNA-binding proteins have in common and homologous region of about 22 amino acids. This region corresponds to two consecutive α-helices that occur in bot Cro and cI repressor proteins of bacteriophage λ and in catabolite gene activator protein ofEscherichia coli and are presumed to interact with DNA. The results obtained here suggest that this α-helical DNA-binding fold occurs in many proteins that regulate gene expression. It also appears that this DNA-binding unit evolved from a common evolutionary precursor.     

Characterization of Maltase Clusters in the Genus <Emphasis Type="Italic">Drosophila</Emphasis>

To reveal evolutionary history of maltase gene family in the genus Drosophila, we undertook a bioinformatics study of maltase genes from available genomes of 12 Drosophila species. Molecular evolution of a closely related glycoside hydrolase, the α-amylase, in Drosophila has been extensively studied for a long time. The α-amylases were even used as a model of evolution of multigene families. On the other hand, maltase, i.e., the α-glucosidase, got only scarce attention. In this study, we, therefore, investigated spatial organization of the maltase genes in Drosophila genomes, compared the amino acid sequences of the encoded enzymes and analyzed the intron/exon composition of orthologous genes. We found that the Drosophila maltases are more numerous than previously thought (ten instead of three genes) and are localized in two clusters on two chromosomes (2L and 2R). To elucidate the approximate time line of evolution of the clusters, we estimated the order and dated duplication of all the 10 genes. Both clusters are the result of ancient series of subsequent duplication events, which took place from 352 to 61 million years ago, i.e., well before speciation to extant Drosophila species. Also observed was a remarkable intron/exon composition diversity of particular maltase genes of these clusters, probably a result of independent intron loss after duplication of intron-rich gene ancestor, which emerged well before speciation in a common ancestor of all extant Drosophila species.     

Morphological Stasis of Two Species Belonging to the L-morphotype in the <Emphasis Type="Italic">Brachionus plicatilis</Emphasis> Species Complex

Detection and characterization of sibling species complexes in zooplankton are critical to understanding their ecological responses and patterns of evolution. The taxon Brachionus plicatilis is a complex of at least 14 species with three major, deeply diverged clades, which are morphologically distinct. We studied morphometric differences between two species – B. plicatilis sensu stricto and B. ‘Manjavacas’– which belong to the L-(large) morphotype and often co-occur in ponds or lakes. B. plicatilis s.s. was on average 6% longer than B. ‘Manjavacas’. They differed significantly in the measurements related to lorica spines. A significant discriminant function relating spine measurements was found, however, individuals from each species showed extensive overlap. Our morphometric data provide additional evidence for the species status of B. plicatilis s.s. and B. ‘Manjavacas’. Since these are ancient species, our results support that a morphological stasis occurs in these taxa. We identified COI restriction sites for PvuII and KpnI which are diagnostic for B. ‘Manjavacas’ and B. plicatilis s.s., respectively. We conclude that morphometry is not useful in classifying the two species. At present, this can only be done reliably using molecular methods.     

Evolution of the scrambled germline gene encoding α-telomere binding protein in three hypotrichous ciliates

The micronuclear genes encoding α-telomere-binding protein (αTP) in Oxytricha trifallax and Stylonychia mytilus contain multiple internal eliminated segments, or IESs, that divide the gene into multiple parts called macronuclear destined segments, or MDSs. The MDSs have become disordered, or scrambled, during evolution. The scrambled structures of the αTP genes in Oxytricha trifallax and S. mytilus have been compared with the previously published scrambled structure of the αTP gene in O. nova. The scrambled patterns of the αTP gene in the three species are similar but show significant differences. The micronuclear genes in O. nova and S. mytilus consist of 13 IESs and 14 MDSs, but the gene in O. trifallax is divided into three additional MDSs by the presence of three additional IESs, believed to have been inserted into the O. trifallaxαTP gene after divergence of O. trifallax from the other two species. Corresponding IESs among the three species have shifted along the DNA during evolution, presumably by a mutational mechanism that changes the short repeat sequences that flank IESs. The IESs also have changed markedly in length by insertion and/or deletion of nucleotides. Comparison of the putative αTP amino acid sequences in the three species reveals three conserved and three nonconserved domains. The 5′ nontranslated regions of the gene-sized molecules encoding αTP contain several conserved segments, and the 3′ nontranscribed trailer contains one conserved segment. Received: 29 May 1998; in revised form: 3 August 1998 / Accepted: 18 August 1998     

Isolation of a gene from Leuconostoc citreum B/110-1-2 encoding a novel dextransucrase enzyme

The amplicon encoding dextransucrase DSR-F from Leuconostoc citreum B/110-1-2, a novel sucrose glucosyltransferase (GTF)-specific for α-1,6 and α-1,3 glucosidic bond synthesis, with α-1,4 branching was cloned, sequenced, and expressed into Escherichia coli JM109. Recombinant enzyme catalyzed oligosaccharides synthesis from sucrose as donor and maltose acceptor. The dsrF gene encodes for a protein (DSR-F) of 1,528 amino acids, with a theoretical molecular mass of 170447.72 Da (~170 kDa). From amino acid sequence comparison, it appears that DSR-F possesses the same domains as those described for GTFs. However, the variable region is longer than in other GTFs (by 100 amino acids) and two APY repeats (a 79 residue long motif with a high number of conserved glycine and aromatic residues, characterized by the presence of the three consecutive residues Ala, Pro, and Tyr) were identified in the glucan binding domain. The DSR-F catalytic domain possesses the catalytic triad involved in the glucosyl enzyme formation. The amino acid sequence of this domain shares a 56% identity with catalytic domain of the alternansucrase ASR from L. citreum NRRL B-1355 and with the catalytic domain of a putative alternansucrase sequence found in the genome of L. citreum KM20. A truncated active variant DSR-F-∆SP-∆GBD of 1,251 amino acids, with a molecular mass of 145 544 Da (~145 kDa), was obtained.     

Origins of gene, genetic code, protein and life: comprehensive view of life systems from a GNC-SNS primitive genetic code hypothesis

We have investigated the origin of genes, the genetic code, proteins and life using six indices (hydropathy, α-helix, β-sheet and β-turn formabilities, acidic amino acid content and basic amino acid content) necessary for appropriate three-dimensional structure formation of globular proteins. From the analysis of microbial genes, we have concluded that newly-born genes are products of nonstop frames (NSF) on antisense strands of microbial GC-rich genes [GC-NSF(a)] and from SNS repeating sequences [(SNS)_n] similar to the GC-NSF(a) (S and N mean G or C and either of four bases, respectively). We have also proposed that the universal genetic code used by most organisms on the earth presently could be derived from a GNC-SNS primitive genetic code. We have further presented the [GADV]-protein world hypothesis of the origin of life as well as a hypothesis of protein production, suggesting that proteins were originally produced by random peptide formation of amino acids restricted in specific amino acid compositions termed as GNC-, SNS and GC-NSF(a)-0th order structures of proteins. The [GADV]-protein world hypothesis is primarily derived from the GNC-primitive genetic code hypothesis. It is also expected that basic properties of extant genes and proteins could be revealed by considerations based on the scenario with four stages This review is a modified English version of the paper, which was written in Japanese and published inViva Origino 2001 29 66–85.     

Molecular Characterization of Mn-superoxide Dismutase and Gene Expression Studies in Dietary Restricted <Emphasis Type="Italic">Brachionus plicatilis</Emphasis> Rotifers

Superoxide dismutases (SODs) promote a conversion of harmful reactive oxygen species (ROS) to relatively moderate forms, resulting in the extension of lifespan in the nematode Caenorhabditis elegans under caloric restriction. The lifespan of the rotifer Brachionus plicatilis is also markedly extended by caloric restriction. We, therefore, cloned cDNA encoding SOD activated with Mn (Mn SOD) from B. plicatilis and examined its expression pattern in rotifers raised with energy restricted diet. The full length deduced amino acid sequence of the rotifer Mn SOD showed 61% identity with the C. elegans ortholog. Four amino acid residues that are essential to the binding of this enzyme to Mn were conserved in the rotifer Mn SOD. Subsequently we examined the mRNA expression patterns of Mn SOD using highly sensitive quantitative real-time PCR for various rotifer populations that are likely to differ in their lifespans in experiments on calorie restricted diets. The accumulated mRNA levels of Mn SOD were found to increase in supposedly long-lived rotifers. These results suggest that Mn SOD is possibly related to the aging of B. plicatilis.     

Concerted evolution of two novel protein families in Caenorhabditis species

Among a large number of homologous gene clusters in C. elegans, two gene families that appear to undergo concerted evolution were studied in detail. Both gene families are nematode specific and encode small secreted proteins of unknown function. For both families in three Caenorhabditis species, concerted groups of genes are characterized by close genomic proximity and by genes in inverted orientation. The rate of protein evolution in one of the two families could be calibrated by comparison with a closely related nonconcerted singleton gene with one-to-one orthologs in all three species. This comparison suggests that protein evolution in concerted gene clusters is two- to sevenfold accelerated. A broader survey of clustered gene families, focused on adjacent inverted gene pairs, identified an additional seven families in which concerted evolution probably occurs. All nine identified families encode relatively small proteins, eight of them encode putative secreted proteins, and most of these have very unusual amino acid composition or sequence. I speculate that these genes encode rapidly evolving antimicrobial peptides.     

Molecular cloning and characterization of five genes encoding pentatricopeptide repeat proteins from Upland cotton (<Emphasis Type="Italic">Gossypium hirsutum</Emphasis> L.)

The pentatricopeptide repeat (PPR) protein family is one of the largest and most complex families in plants. These proteins contain multiple 35-amino acid repeats that are proposed to form a super helix capable of binding RNA. PPR proteins have been implicated in many crucial functions broadly involving organelle biogenesis and plant development. In this study, we identified many genes encoding PPR protein in Upland cotton through an extensive survey of the database of Gossypium hirsutum. Furthermore, we isolated five full-length cDNA of PPR genes from G. hirsutum 0-613-2R which were named GhPPR1–GhPPR5. Domain analysis revealed that the deduced amino acid sequences of GhPPR1–5 contained from 5 to 10 PPR motifs and those PPR proteins were divided into two different PPR subfamilies. GhPPR1–2 belonged to the PLS subfamily and GhPPR3–5 belonged to the P subfamily. Phylogenetic analysis of the five GhPPR proteins and 18 other plant PPR proteins also revealed that the same subfamily clustered together. All five GhPPR genes were differentially but constitutively expressed in roots, stems, leaves, pollens, and fibers based on the gene expression analysis by real-time quantitative RT-PCR. This study is the first report and analysis of genes encoding PPR proteins in cotton.