Evolution of Microsatellite Alleles in Four Species of Mice (Genus Apodemus)

Microsatellite length variation was investigated at a highly variable microsatellite locus in four species of Apodemus. Information obtained from microsatellite allele sequences was contrasted with allele sizes, which included 18 electromorphs. Additional analysis of a 400-bp unique sequence in the flanking region identified 26 different haplotype sequences or ``true' alleles in the sample. Three molecular mechanisms, namely, (1) addition/deletion of repeats, (2) substitutions and indels in the flanking region, and (3) mutations interrupting the repeat, contributed to the generation of allelic variation. Size homoplasy can be inferred for alleles within populations, from different populations of the same species, and from different species. We propose that microsatellite flanking sequences may be informative markers for investigating mutation processes in microsatellite repeats as well as phylogenetic relationships among alleles, populations, and species. Received: 3 November 1999 / Accepted: 2 May 2000     

Variation of Microsatellite Size Homoplasy Across Electromorphs, Loci, and Populations in Three Invertebrate Species

Size homoplasy was analyzed at microsatellite loci by sequencing electromorphs, that is, variants of the same size (base pairs). This study was conducted using five interrupted and/or compound loci in three invertebrate species, the honey bee Apis mellifera, the bumble bee Bombus terrestris, and the freshwater snail Bulinus truncatus. The 15 electromorphs sequenced turned out to hide 31 alleles (i.e., variants identical in sequence). Variation in the amount of size homoplasy was detected among electromorphs and loci. From one to seven alleles were detected per electromorph, and one locus did not show any size homoplasy in both bee species. The amount of size homoplasy was related to the sequencing effort, since the number of alleles was correlated with the number of copies of electromorphs sequenced, but also with the molecular structure of the core sequence at each locus. Size homoplasy within populations was detected only three times, meaning that size homoplasy was detected mostly among populations. We analyzed population structure, estimating F _st and a genetic distance, based on either electromorphs or alleles. Whereas little difference was found in A. mellifera, uncovering size homoplasy led to a more marked population structure in B. terrestris and B. truncatus. We also showed in A. mellifera that the detection of size homoplasy may alter phylogenetic reconstructions. Received: 21 July 1997 / Accepted: 29 January 1998     

Microsatellite Evolution: Testing the Ascertainment Bias Hypothesis

Previous studies suggest the median allele length of microsatellites is longest in the species from which the markers were derived, suggesting that an ascertainment bias was operating. We have examined whether the size distribution of microsatellite alleles between sheep and cattle is source dependent using a set of 472 microsatellites that can be amplified in both species. For those markers that were polymorphic in both species we report a significantly greater number of markers (P < 0.001) with longer median allele sizes in sheep, regardless of microsatellite origin. This finding suggests that any ascertainment bias operating during microsatellite selection is only a minor contributor to the variation observed. Received: 6 January 1997 / Accepted: 19 May 1997     

A Phylogenetic Perspective on Sequence Evolution in Microsatellite Loci

We examined the evolution of the repeat regions of three noncoding microsatellite loci in 58 species of the Polistinae, a subfamily of wasps that diverged over 140 million years ago. A phylogenetic approach allows two new kinds of approaches to studying microsatellite evolution: character mapping and comparative analysis. The basic repeat structure of the loci was highly conserved, but was often punctuated with imperfections that appear to be phylogenetically informative. Repeat numbers evolved more rapidly than other changes in the repeat region. Changes in number of repeats among species seem consistent with the stepwise mutation model, which is based on slippage during replication as the main source of mutations. Changes in repeat numbers can occur even when there are very few tandem repeats but longer repeats, especially perfect repeats led to greater rates of evolutionary change. Species phylogenetically closer to the one from which we identified the loci had longer stretches of uninterrupted repeats and more different motifs, but not longer total repeat regions. The number of perfect repeats increased more often than it decreased. However, there was no evidence that some species have consistently greater numbers of repeats across loci than other species have, once ascertainment bias is eliminated. We also found no evidence for a population size effect posited by one form of the directionality hypothesis. Overall, phylogenetic variation in repeat regions can be explained by adding neutral evolution to what is already known about the mutation process. The life cycle of microsatellites appears to reflect a balance between growth by slippage and degradation by an essentially irreversible accumulation of imperfections. Received: 13 April 1999 / Accepted: 8 September 1999     

Variable Substitution Rates of the 18 Domain Sequences in Artemia Hemoglobin

The Artemia hemoglobin is a dimer comprising two nine-domain covalent polymers in quaternary association. Each polymer is encoded by a gene representing nine successive globin domains which have different sequences and are presumed to have been copied originally from a single-domain gene. Two different polymers exist as the result of a complete duplication of the nine-domain gene, allowing the formation of either homodimers or the heterodimer. The total population size of 18 domains comprising nine corresponding pairs, coupled with the probability that they reflect several hundred million years of evolution in the same lineage, provides a unique model in which the process of gene multiplication can be analyzed. The outcome has important implications for the reliability of local molecular clocks. The two polymers differ from each other at 11.7% of amino acid sites; however when corresponding individual domains are compared between polymers, amino acid substitution fluctuates by a factor of 2.7-fold from lowest to highest. This variation is not obvious at the DNA level: Domain pair identity values fluctuate by 1.3-fold. Identity values are, however, uncorrected for multiple substitutions, and both silent and nonsilent changes are pooled. Therefore, to determine the variability in relative substitution rates at the DNA level, we have used the method of Li (1993, J Mol Evol 36:96–99) to determine estimates of nonsynonymous (K _A ) and synonymous (K _S ) substitutions per site for the nine pairs of domains. As expected, the overall level of silent substitutions (K _S of 56.9%) far exceeded nonsilent substitutions (K _A of 6.7%); however, for corresponding domain pairs, K _A fluctuates by 2.3-fold and K _S by 1.7-fold. The large discrepancies reflected in the expressed protein have accrued within a single lineage and the implication is that divergence dates of different genera based on amino acid sequences, even with well-studied proteins of reasonable size, can be wrong by a factor well in excess of 2. Received: 4 June 1997 / Accepted: 17 December 1997     

Evidence of Gene Conversion Events Between Paralogous Sequences Produced by Tetraploidization in Salmoninae Fish

We investigated the occurrence of gene conversions between paralogous sequences of Salmoninae derived from ancestral tetraploidization and their effect on the evolutionary history of DNA sequences. A microsatellite with long flanking regions (750 bp) including both coding and noncoding sequences was analyzed. Microsatellite size polymorphism was used to detect the alleles of both paralogous counterparts and infer linkage arrangement between loci. DNA sequencing of seven Salmoninae species revealed that paralogous sequences were highly differentiated within species, especially for noncoding regions. Ten gene conversion events between paralogous sequences were inferred. While these events appears to have homogenized regions of otherwise highly differential paralogous sequences, they amplified the differentiation among orthologous sequences. Their effects were larger on coding than on noncoding regions. As a consequence, noncoding sequences grouped by orthologous lineages in phylogenetic trees, whereas coding regions grouped by taxa. Based upon these results, we present a model showing how gene conversion events may also result in the PCR amplification of nonorthologous sequences in different taxa, with obvious complications for phylogenetic inferences, comparative mapping, and population genetic studies. Received: 11 October 2000 / Accepted: 18 September 2001     

Identification of Sequences Encoding the Detoxification Metalloisomerase Glyoxalase I in Microbial Genomes from Several Pathogenic Organisms

The ubiquitous glyoxalase system, which is composed of two enzymes, removes cellular cytotoxic methylglyoxal (MG). In an effort to identify critical residues conserved in the evolution of the first enzyme in this system, glyoxalase I (GlxI), as well as the structural implications of sequence alterations in this enzyme, a search of the National Center for Biotechnology Information (NCBI) database of unfinished genomes was undertaken. Eleven putative GlxI sequences from pathogenic organisms were identified and analyses of these sequences in relation to the known and previously identified GlxI enzymes were performed. Several of these sequences show a very high similarity to the Escherichia coli GlxI sequence, most notably the 79% identity of the sequence identified from Yersinia pestis, the causative agent of bubonic plague. In addition to the conservation of residues critical to binding the catalytic metal in all of the proposed GlxI enzymes, four regions in the Homo sapiens GlxI enzyme are absent in all of the bacterial GlxI sequences, with the exception of Pseudomonas putida. Removal of these regions may alter the active-site conformation of the bacterial enzymes in relation to that of the H. sapiens. These differences may be targeted for the development of inhibitors selective to the bacterial enzymes. Received: 13 October 1999 / Accepted: 17 January 2000     

Identification of Glyoxalase I Sequences in Brassica oleracea and Sporobolus stapfianus: Evidence for Gene Duplication Events

Glyoxalase I (GlxI) is the first of two enzymes involved in the cellular detoxification of methylglyoxal. A recent search of the National Center for Biotechnology Information (NCBI) databases with the protein sequence of Salmonella typhimurium GlxI identified two new hypothetical proteins with unassigned function. These two sequences, from Brassica oleracea and Sporobolus stapfianus, have significant sequence similarity to known GlxI sequences, suggesting that these two open reading frames encode for GlxI in these plants. Interestingly, analysis of these two new sequences indicates that they code for a protein composed of two fused monomers, a situation previously found solely in the yeast GlxI enzymes. Received: 10 May 1997 / Accepted: 15 October 1997     

Mhc Allelic Diversity and Modern Human Origins

Thirty complete coding sequences of human major histocompatibility complex (Mhc) class II DRB alleles, spanning 237 codons, were analyzed for phylogenetic information using distance, parsimony, and likelihood approaches. Allelic genealogies derived from different parts of the coding sequence (exon 2, the 5′ and 3′ ends of exon 2, respectively, and exons 3–6) were compared. Contrary to prior assertions, a rigorous analysis of allelic genealogies in this gene family cannot be used to justify the claim that the lineage leading to modern humans contained on average at least 100,000 individuals. Phylogenetic inferences based upon the exon 2 region of the DRB loci are complicated by selection and recombination, so this part of the gene does not provide a complete and accurate view of allelic relationships. Attempts to reconstruct human history from genetic data must use realistic models which consider the complicating factors of nonequilibrium populations, recombination, and different patterns of selection. Received: 19 February 1997 / Accepted: 12 June 1997     

Taxonomy of 5 S Ribosomal RNA by the Linguistic Technique: Probing with Mitochondrial and Mammalian Sequences

Linguistic similarities and dissimilarities between 5 S rRNA sequences allowed taxonomical separation of species and classes. Comparisons with the molecule from mammals distinguished fungi and plants from protists and animals. Similarities to mammalians progressively increased from protists to invertebrates and to somatic-type molecules of the vertebrates lineage. In this, deviations were detected in avian, oocyte type, and pseudogene sequences. Among bacteria, actinobacteria were most similar to the mammalians, which could be related to the high frequency of associations among members of these groups. Some archaebacterial species most similar to the mammalians belonged to the Thermoproteales and Halobacteria groups. Comparisons with the soybean mitochondrial molecule revealed high internal homogeneity among plant mitochondria. The eubacterial groups most similar to it were Thermus and Rhodobacteria γ-1 and α-2. Other procedures have already indicated similarities of Rhodobacteria α to mitochondria but the linguistic similarities were on the average higher with the first two groups. Received: 5 August 1996 / Accepted: 9 April 1997     

A Novel Approach to Phylogeny Reconstruction from Protein Sequences

The reliable reconstruction of tree topology from a set of homologous sequences is one of the main goals in the study of molecular evolution. If consistent estimators of distances from a multiple sequence alignment are known, the distance method is attractive because the tree reconstruction is consistent. To obtain a distance estimate d, the observed proportion of differences p (p-distance) is usually ``corrected' for multiple and back substitutions by means of a functional relationship d=f(p). In this paper the conditions under which this correction of p-distances will not alter the selection of the tree topology are specified. When these conditions are not fulfilled the selection of the tree topology may depend on the correction function applied. A novel method which includes estimates of distances not only between sequence pairs, but between triplets, quadruplets, etc., is proposed to strengthen the proper selection of correction function and tree topology. A ``super' tree that includes all tree topologies as special cases is introduced. Received: 17 February 1998 / Accepted: 20 July 1998     

Divergent Human Y-Chromosome Microsatellite Evolution Rates

In this work, we analyze several characteristics influencing the low variability of the microsatellite DYS19 in the major founder Amerindian Y chromosome lineage containing the point mutation DYS199-T. Variation of DYS19 was compared with that of five other Y-linked tetranucleotide repeat loci (DYS389A, DYS389B, DYS390, DYS391, and DYS393) in the DYS199-T lineage. All the other microsatellites showed significantly higher levels of variability than DYS19 as measured by gene diversity and repeat number variance. Moreover, we had previously shown that DYS19 had high diversity in Brazilians and in several other populations worldwide. Thus, the slow DYS19 evolution in the DYS199-T lineage seems to be both locus and allele specific. To understand the slow DYS19 evolutionary rate, the microsatellite loci were compared according to their mapping on the Y chromosome and also on the basis of structural aspects such as the base composition of the repeat motif and flanking regions and the degree of perfection and size (repeat number) of the variable blocks. The only observed difference that might be related to the low DYS19 variability is its small average number of repeats, a value expected to be closer to the founder DYS19 allele in the DYS199-T lineage. These data were also compared to other derived Y lineages. The Tat-C lineage displayed a lower DYS19 variability correlated to a small average repeat number, while in the DYS234-G lineage, a high DYS19 variability was found associated to a larger average repeat number. This approach reveals that evolution of Y microsatellites in lineages defined by slowly evolving markers, such as point mutations, can be greatly influenced by the size (number of repeats of the variable block) of the founder allele in each microsatellite locus. Thus lineage-dating methods using microsatellite variation should be practiced with great care. Received: 7 November 1998 / Accepted: 9 April 1999     

Unusual and Strongly Structured Sequence Variation in a Complex Satellite DNA Family from the Nematode Meloidogyne chitwoodi

An AluI satellite DNA family has been isolated in the genome of the root-knot nematode Meloidogyne chitwoodi. This repeated sequence was shown to be present at approximately 11,400 copies per haploid genome, and represents about 3.5% of the total genomic DNA. Nineteen monomers were cloned and sequenced. Their length ranged from 142 to 180 bp, and their A + T content was high (from 65.7 to 79.1%), with frequent runs of As and Ts. An unexpected heterogeneity in primary structure was observed between monomers, and multiple alignment analysis showed that the 19 repeats could be unambiguously clustered in six subfamilies. A consensus sequence has been deduced for each subfamily, within which the number of positions conserved is very high, ranging from 86.7% to 98.6%. Even though blocks of conserved regions could be observed, multiple alignment of the six consensus sequences did not enable the establishment of a general unambiguous consensus sequence. Screening of the six consensus sequences for evidence of internal repeated subunits revealed a 6-bp motif (AAATTT), present in both direct and inverted orientation. This motif was found up to nine times in the consensus sequences, also with the occurrence of degenerated subrepeats. Along with the meiotic parthenogenetic mode of reproduction of this nematode, such structural features may argue for the evolution of this satellite DNA family either (1) from a common ancestral sequence by amplification followed by mechanisms of sequence divergence, or (2) through independent mutations of the ancestral sequence in isolated amphimictic nematode populations and subsequent hybridization events. Overall, our results suggest the ancient origin of this satellite DNA family, and may reflect for M. chitwoodi a phylogenetic position close to the ancestral amphimictic forms of root-knot nematodes. Received: 23 April 1997 / Accepted: 9 July 1997     

Compositional Properties of Homologous Coding Sequences from Plants

In this work, we investigated (1) the compositional distributions of all available nuclear coding sequences (and of their three codon positions) of six dicots and four Gramineae; this considerably expanded our knowledge about the differences previously seen between these two groups of plants; (2) the compositional correlations of homologous genes from dicots and from Gramineae, as well as from both groups; all correlations were characterized by very good coefficients, with slopes close to unity in the former two cases and very high in the last; (3) the compositional transition that accompanied the emergence of Gramineae from an ancestral monocot; (4) the compositional correlations between exons and introns, which were very good in Gramineae, but only poor to good in dicots; and (5) the compositional profiles of homologous genes from angiosperms, which were characterized by a series of peaks (exons) and valleys (introns) separated by 15–20% GC. The conservative and transitional modes of compositional evolution in plant genes and their general implications are discussed. Received: 24 June 1997 / Accepted: 20 August 1997     

The Structure and Gene Repertoire of an Ancient Red Algal Plastid Genome

Photosynthetic eukaryotes can, according to features of their chloroplasts, be divided into two major groups: the red and the green lineage of plastid evolution. To extend the knowledge about the evolution of the red lineage we have sequenced and analyzed the chloroplast genome (cp-genome) of Cyanidium caldarium RK1, a unicellular red alga (AF022186). The analysis revealed that this genome shows several unusual structural features, such as a hypothetical hairpin structure in a gene-free region and absence of large repeat units. We provide evidence that this structural organization of the cp-genome of C. caldarium may be that of the most ancient cp-genome so far described. We also compared the cp-genome of C. caldarium to the other known cp-genomes of the red lineage. The cp-genome of C. caldarium cannot be readily aligned with that of Porphyra purpurea, a multicellular red alga, or Guillardia theta due to a displacement of a region of the cp-genome. The phylogenetic tree reveals that the secondary endosymbiosis, through which G. theta evolved, took place after the separation of the ancestors of C. caldarium and P. purpurea. We found several genes unique to the cp-genome of C. caldarium. Five of them seem to be involved in the building of bacterial cell envelopes and may be responsible for the thermotolerance of the chloroplast of this alga. Two additional genes may play a role in stabilizing the photosynthetic machinery against salt stress and detoxification of the chloroplast. Thus, these genes may be unique to the cp-genome of C. caldarium and may be required for the endurance of the extreme living conditions of this alga. Received: 3 June 2000 / Accepted: 18 July 2000     

Origin of the NEFA and Nuc Signal Sequences

The human protein NEFA binds calcium, contains a leucine zipper repeat that does not form a homodimer, and is proposed (along with the homologous Nuc protein) to have a common evolutionary history with an EF-hand ancestor. We have isolated and characterized the N-terminal domain of NEFA that contains a signal sequence inferred from both endoproteinase Asp-N (Asp-N) and tryptic digests. Analysis of this N-terminal sequence shows significant similarity to the conserved multiple domains of the mitochondrial carrier family (MCF) proteins. The leader sequence of Nuc is, however, most similar to the signal sequences of membrane and/or secreted proteins (e.g., mouse insulin-like growth factor receptor). We suggest that the divergent NEFA and Nuc N-terminal sequences may have independent origins and that the common high hydrophobicity governs their targeting to the ER. These results provide insights into signal sequence evolution and the multiple origins of protein targeting. Received: 20 February 1997 / Accepted: 28 July 1997     

Allelic variation in the squirrel monkey x-linked color vision gene: biogeographical and behavioral correlates

Most Neotropical primate species possess a polymorphic X-linked and a monomorphic autosomal color vision gene. Consequently, populations are composed of both dichromatics and trichromatics. Most theories on the maintenance of this genetic system revolve around possible advantages for foraging ecology. To examine the issue from a different angle, we compared the numbers and relative frequencies of alleles at the X-linked locus among three species of Saimiri representing a wide range of geographical and behavioral variation in the genus. Exons 3, 4, and 5 of the X-linked opsin gene were sequenced for a large number of X chromosomes for all three species. Several synonymous mutations were detected in exons 4 and 5 for the originally reported alleles but only a single nonsynonymous change was detected. Two alleles were found that appeared to be the result of recombination events. The low occurrence of recombinant alleles and absence of mutations in the amino acids critical for spectral tuning indicates that stabilizing selection acts to maintain the combinations of critical sites specific to each allele. Allele frequencies were approximately the same for all Saimiri species, with a slight but significant difference between S. boliviensis and S. oerstedii. No apparent correlation exists between allele frequencies and behavioral or biogeographical differences between species, casting doubt on the speculation that the spectral sensitivities of the alleles have been maintained because they are specifically well-tuned to Saimiri visual ecology. Rather, the spectral tuning peaks might have been maintained because they are as widely spaced as possible within the limited range of middlewave to longwave spectra useful to all primates. This arrangement creates a balance between maximizing the distance between spectral tuning peaks (allowing the color opponency of the visual system to distinguish between peaks) and maximizing the number of alleles within a limited range (yielding the greatest possible frequency of heterozygotes).     

rBAT is an Amino Acid Exchanger with Variable Stoichiometry

The rBAT protein, when expressed in Xenopus oocytes, was previously shown to reproduce the selectivity of the Na⁺-independent neutral and basic amino acid transport system called b^o,+. More recently, the capacity of rBAT to generate a transmembrane current was demonstrated when addition of neutral amino acids stimulated the efflux of cations (presumably basic amino acids) in rBAT-injected oocytes. In the present paper, aminoisobutyric acid (AIB), a neutral amino acid analogue, was shown to induce outward currents (efflux of basic amino acids) through rBAT similar to those caused by alanine in terms of affinity, maximal currents and I-V curves. Despite generating similar currents, the AIB transport rate was more than 30 times lower than that of alanine, thus challenging the assumption that rBAT functions as a classical exchanger. Experiments using a cut-open oocyte voltage clamp demonstrated that AIB was capable of stimulating rBAT-mediated currents from either side of the membrane. AIB, like alanine, was able to stimulate the efflux of radiolabeled alanine and arginine while no rBAT-mediated efflux was measurable in the absence of external rBAT substrates. These results demonstrate that (i) the presence of amino acids is required on both sides of the membrane for rBAT to mediate amino acid flux and thus rBAT must be some type of exchanger but (ii) rBAT-mediated amino acid influx is not stoichiometrically related to the efflux. A model of a ``double gated pore' is proposed to account for these properties of rBAT, which contravene standard models of exchangers and other transporters. Received: 15 June 1995/Revised: 21 September 1995     

Codon Bias and Mutability in HIV Sequences

A survey of the patterns of synonymous codon preference in the HIV env gene reveals a correlation between the codon bias and the mutability requirements of different regions of the protein. At hypervariable regions in gp120 one finds a greater proportion of codons that tend to mutate nonsynonymously, but to a target that is similar in hydrophobicity and volume. We argue that this strategy results from a compromise between the selective pressure placed on the virus by the induced immune response, which favors amino acid substitutions in the complementarity determining regions, and the negative selection against missense mutations that violate structural constraints of the env protein. Received: 9 June 1997 / Accepted: 25 May 1998