The evolution of homing endonuclease genes and group I introns in nuclear rDNA

Group I introns are autonomous genetic elements that can catalyze their own excision from pre-RNA. Understanding how group I introns move in nuclear ribosomal (r)DNA remains an important question in evolutionary biology. Two models are invoked to explain group I intron movement. The first is termed homing and results from the action of an intron-encoded homing endonuclease that recognizes and cleaves an intronless allele at or near the intron insertion site. Alternatively, introns can be inserted into RNA through reverse splicing. Here, we present the sequences of two large group I introns from fungal nuclear rDNA, which both encode putative full-length homing endonuclease genes (HEGs). Five remnant HEGs in different fungal species are also reported. This brings the total number of known nuclear HEGs from 15 to 22. We determined the phylogeny of all known nuclear HEGs and their associated introns. We found evidence for intron-independent HEG invasion into both homologous and heterologous introns in often distantly related lineages, as well as the "switching" of HEGs between different intron peripheral loops and between sense and antisense strands of intron DNA. These results suggest that nuclear HEGs are frequently mobilized. HEG invasion appears, however, to be limited to existing introns in the same or neighboring sites. To study the intron-HEG relationship in more detail, the S943 group I intron in fungal small-subunit rDNA was used as a model system. The S943 HEG is shown to be widely distributed as functional, inactivated, or remnant ORFs in S943 introns.     

Molecular evolution of ubiquitin genes

Ubiquitin is a singular protein with multiple functions. It is probably the most slowly evolving protein known, is encoded by genes with a unique structure, and provides an intriguing case study for various aspects of molecular evolution. In particular, the multiple ubiquitin-coding repeats which have been characterized in man, yeast and a slime mould graphically illustrate the dynamics of concerted evolution, but cast doubts on the effectiveness of this process for unlinked arrays in this repeat family.     

Molecular evolution of a multigene family in group A streptococci

The emm genes are members of a gene family in group A streptococci (GAS) that encode for antiphagocytic cell-surface proteins and/or immunoglobulin-binding proteins. Previously sequenced genes in this family have been named "emm," "fcrA," "enn," "arp," "protH," and "mrp"; herein they will be referred to as the "emm gene family." The genes in the emm family are located in a cluster occupying 3-6 kb between the genes mry and scpA on the chromosome of Streptococcus pyogenes. Most GAS strains contain one to three tandemly arranged copies of emm-family genes in the cluster, but the alleles within the cluster vary among different strains. Phylogenetic analysis of the conserved sequences at the 3' end of these genes differentiates all known members of this family into four evolutionarily distinct emm subfamilies. As a starting point to analyze how the different subfamilies are related evolutionarily, the structure of the emm chromosomal region was mapped in a number of diverse GAS strains by using subfamily-specific primers in the polymerase chain reaction. Nine distinct chromosomal patterns of the genes in the emm gene cluster were found. These nine chromosomal patterns support a model for the evolution of the emm gene family in which gene duplication followed by sequence divergence resulted in the generation of four major-gene subfamilies in this locus.      

Molecular evolution of cycloidea-like genes in Fabaceae

The cycloidea (CYC) gene controls floral symmetry in snapdragon (Antirrhinum majus). We investigated the evolution of CYC-like genes in some species of legumes that have zygomorphic flowers. Two to four CYC-like genes were isolated from a single species. The results of NJ and ML analyses indicate that CYC-like genes in legumes group into two monophyletic clades; one group consists of eight CYC-like genes (Clade 1) and the other contains three CYC-like genes and TB1 of maize (Clade 2). These phylogenetic trees and the Shimodaira–Hasegawa test suggest that Clade 1 is a sister of the original CYC group (Clade 3). Moreover, the result of the GeneTree analysis showed that the CYC-like genes experienced repeated duplication events during the evolution of legumes. We herein speculate as to the role of CYC-like genes in legumes and discuss the evolutionary processes that these genes have undergone. Current address (Jun Yokoyama and Masayuki Maki): Division of Ecology and Evolutionary Biology, Graduate School of Life Sciences, Tohoku University, Aoba, Sendai 980-8578, Japan     

Molecular evolution of genes in avian genomes

Background

Obtaining a draft genome sequence of the zebra finch (Taeniopygia guttata), the second bird genome to be sequenced, provides the necessary resource for whole-genome comparative analysis of gene sequence evolution in a non-mammalian vertebrate lineage. To analyze basic molecular evolutionary processes during avian evolution, and to contrast these with the situation in mammals, we aligned the protein-coding sequences of 8,384 1:1 orthologs of chicken, zebra finch, a lizard and three mammalian species.

Results

We found clear differences in the substitution rate at fourfold degenerate sites, being lowest in the ancestral bird lineage, intermediate in the chicken lineage and highest in the zebra finch lineage, possibly reflecting differences in generation time. We identified positively selected and/or rapidly evolving genes in avian lineages and found an over-representation of several functional classes, including anion transporter activity, calcium ion binding, cell adhesion and microtubule cytoskeleton.

Conclusions

Focusing specifically on genes of neurological interest and genes differentially expressed in the unique vocal control nuclei of the songbird brain, we find a number of positively selected genes, including synaptic receptors. We found no evidence that selection for beneficial alleles is more efficient in regions of high recombination; in fact, there was a weak yet significant negative correlation between ω and recombination rate, which is in the direction predicted by the Hill-Robertson effect if slightly deleterious mutations contribute to protein evolution. These findings set the stage for studies of functional genetics of avian genes.     

Molecular evolution of sex-biased genes in Drosophila

Studies of morphology, interspecific hybridization, protein/DNA sequences, and levels of gene expression have suggested that sex-related characters (particularly those involved in male reproduction) evolve rapidly relative to non-sex-related characters. Here we report a general comparison of evolutionary rates of sex-biased genes using data from cDNA microarray experiments and comparative genomic studies of Drosophila. Comparisons of nonsynonymous/synonymous substitution rates (d(N)/d(S)) between species of the D. melanogaster subgroup revealed that genes with male-biased expression had significantly faster rates of evolution than genes with female-biased or unbiased expression. The difference was caused primarily by a higher d(N) in the male-biased genes. The same pattern was observed for comparisons among more distantly related species. In comparisons between D. melanogaster and D. pseudoobscura, genes with highly biased male expression were significantly more divergent than genes with highly biased female expression. In many cases, orthologs of D. melanogaster male-biased genes could not be identified in D. pseudoobscura through a Blast search. In contrast to the male-biased genes, there was no clear evidence for accelerated rates of evolution in female-biased genes, and most comparisons indicated a reduced rate of evolution in female-biased genes relative to unbiased genes. Male-biased genes did not show an increased ratio of nonsynonymous/synonymous polymorphism within D. melanogaster, and comparisons of polymorphism/divergence ratios suggest that the rapid evolution of male-biased genes is caused by positive selection.     

Molecular evolution of nitrate reductase genes

To understand the evolutionary mechanisms and relationships of nitrate reductases (NRs), the nucleotide sequences encoding 19 nitrate reductase (NR) genes from 16 species of fungi, algae, and higher plants were analyzed. The NR genes examined show substantial sequence similarity, particularly within functional domains, and large variations in GC content at the third codon position and intron number. The intron positions were different between the fungi and plants, but conserved within these groups. The overall and nonsynonymous substitution rates among fungi, algae, and higher plants were estimated to be 4.33 × 10⁻¹⁰ and 3.29 × 10⁻¹⁰ substitutions per site per year. The three functional domains of NR genes evolved at about one-third of the rate of the N-terminal and the two hinge regions connecting the functional domains. Relative rate tests suggested that the nonsynonymous substitution rates were constant among different lineages, while the overall nucleotide substitution rates varied between some lineages. The phylogenetic trees based on NR genes correspond well with the phylogeny of the organisms determined from systematics and other molecular studies. Based on the nonsynonymous substitution rate, the divergence time of monocots and dicots was estimated to be about 340 Myr when the fungi–plant or algae–higher plant divergence times were used as reference points and 191 Myr when the rice–barley divergence time was used as a reference point. These two estimates are consistent with other estimates of divergence times based on these reference points. The lack of consistency between these two values appears to be due to the uncertainty of the reference times. Received: 10 April 1995 / Accepted: 10 September 1995     

Molecular evolution: Codes,clocks, genes and genomes

The discoveries, advancements and continuing controversies in the field of molecular evolution are reviewed. Topics summarized are (1) the evolution of the genetic code, (2) gene evolution including the demonstration of homology, estimation of sequence divergence, phylogenetic trees, the molecular clock and the origin of genes and gene families by various genetic mechanisms, and (3) eukaryotic genome evolution, including the highly repeated satellite sequences, the interspersed and potentially mobile repeated sequences and the unique sequence fraction of the genome.     

Molecular evolution of immunoglobulin superfamily genes in primates

Genes of the immunoglobulin superfamily (IgSF) have a wide variety of cellular activities. In this study, we investigated molecular evolution of IgSF genes in primates by comparing orthologous sequences of 249 IgSF genes among human, chimpanzee, orangutan, rhesus macaque, and common marmoset. To evaluate the non-synonymous/synonymous substitution ratio (ω), we applied Bn-Bs program and PAML program. IgSF genes were classified into 11 functional categories based on the Gene Ontology (GO) database. Among them, IgSF genes in three functional categories, immune system process (GO:0002376), defense response (GO:0006952), and multi-organism process (GO:0051704), which are tightly linked to the regulation of immune system had much higher values of ω than genes in the other GO categories. In addition, we estimated the average values of ω for each primate lineage. Although each primate lineage had comparable average values of ω, the human lineage showed the lowest ω value for the immune-related genes. Furthermore, 11 IgSF genes, SIGLEC5, SLAMF6, CD33, CD3E, CEACAM8, CD3G, FCER1A, CD48, CD4, TIM4, and FCGR2A, were implied to have been under positive selective pressure during the course of primate evolution. Further sequence analyses of CD3E and CD3G from 23 primate species suggested that the Ig domains of CD3E and CD3G underwent the positive Darwinian selection.     

Molecular evolution of pathogenicity-island genes in Pseudomonas viridiflava

The bacterial pathogen Pseudomonas viridiflava possesses two pathogenicity islands (PAIs) that share many gene homologs, but are structurally and phenotypically differentiated (T-PAI and S-PAI). These PAIs are paralogous, but only one is present in each isolate. While this dual presence/absence polymorphism has been shown to be maintained by balancing selection, little is known about the molecular evolution of individual genes on the PAIs. Here we investigate genetic variation of 12 PAI gene loci (7 on T-PAI and 5 on S-PAI) in 96 worldwide isolates of P. viridiflava. These genes include avirulence genes (hopPsyA and avrE), their putative chaperones (shcA and avrF), and genes encoding the type III outer proteins (hrpA, hrpZ, and hrpW). Average nucleotide diversities in these genes (pi = 0.004-0.020) were close to those in the genetic background. Large numbers of recombination events were found within PAIs and a sign of positive selection was detected in avrE. These results suggest that the PAI genes are evolving relatively freely from each other on the PAIs, rather than as a single unit under balancing selection. Evolutionarily stable PAIs may be preferable in this species because preexisting genetic variation enables P. viridiflava to respond rapidly to natural selection.     

Molecular considerations in the evolution of bacterial genes

Summary Synonymous and nonsynonymous substitution rates at the loci encoding glyceraldehyde-3-phosphate dehydrogenase (gap) and outer membrane protein 3A (ompA) were examined in 12 species of enteric bacteria. By examining homologous sequences in species of varying degrees of relatedness and of known phylogenetic relationships, we analyzed the patterns of synonymous and nonsynonymous substitutions within and among these genes. Although both loci accumulate synonymous substitutions at reduced rates due to codon usage bias, portions of thegap andompA reading frames show significant deviation in synonymous substitution rates not attributable to local codon bias. A paucity of synonymous substitutions in portions of theompA gene may reflect selection for a novel mRNA secondary structure. In addition, these studies allow comparisons of homologous protein-coding sequences (gap) in plants, animals, and bacteria, revealing differences in evolutionary constraints on this glycolytic enzyme in these lineages.     

Molecular evolution of Drosophila odorant receptor genes

A total of 752 odorant receptor (Or) genes, including pseudogenes, were identified in 11 Drosophila species and named after their orthologs in Drosophila melanogaster. The 813 Or genes, including 61 from D. melanogaster, were classified into 59 orthologous groups that are well supported by gene phylogeny. By reconciling with the gene family phylogeny, we estimated the number of gene duplication/loss events and intron gain/loss events in the species phylogeny. We found that these events are particularly frequent in Drosophila grimshawi, Drosophila willistoni, and obscura group. More than half of the duplicated genes stay as tandem arrays, whose size range from 2 to 8. These genes vary in sequence and some likely underwent positive selection, indicating that the gene duplication was important for flies to acquire new olfactory functions. We hypothesize that Or genes conferred the basic olfactory repertoire to ancestral flies before the speciation of the Drosophila and Sophophora subgenera about 40 Mya. This repertoire has been largely maintained in the current species, whereas lineage-specific gene duplication seems to have led to additional specialization in some species in response to specific ecological conditions.     

Molecular evolution of bat color vision genes

The two suborders of bats, Megachiroptera (megabats) and Microchiroptera(microbats), use different sensory modalities for perceivingtheir environment. Megabats are crepuscular and rely on a well-developedeyes and visual pathway, whereas microbats occupy a nocturnalniche and use acoustic orientation or echolocation more thanvision as the major means of perceiving their environment. Inview of the differences associated with their sensory systems,we decided to investigate the function and evolution of colorvision (opsin genes) in these two suborders of bats. The middle/longwavelength (M/L) and short wavelength (S) opsin genes were sequencedfrom two frugivorous species of megabats, Haplonycteris fischeriand Pteropus dasymallus formosus, and one insectivorous speciesof microbat, Myotis velifer. Contrary to the situation in primates,where many nocturnal species have lost the functional S opsingene, both crepuscular and strictly nocturnal species of batsthat we examined have functional M/L and S opsin genes. Surprisingly,the S opsin in these bats may be sensitive to UV light, whichis relatively more abundant at dawn and at dusk. The M/L opsinin these bats appears to be the L type, which is sensitive tored and may be helpful for identifying fruits among leaves orfor other purposes. Most interestingly, H. fischeri has a recentduplication of the M/L opsin gene, representing to date theonly known case of opsin gene duplication in non-primate mammals.Some of these observations are unexpected and may provide insightsinto the effect of nocturnal life on the evolution of opsingenes in mammals and the evolution of the life history traitsof bats in general.     

Molecular evolution of human visual pigment genes

By comparing the published DNA sequences for (a) the genes encoding the human visual color pigments (red, green, and blue) with (b) the genes encoding human, bovine, and Drosophila rhodopsins, a phylogenetic tree for the mammalian pigment genes has been constructed. This evolutionary tree shows that the common ancestor of the visual color pigment genes diverged first from that of the rhodopsin genes; then the common ancestor of the red and green pigment genes and the ancestor of the blue pigment gene diverged; and finally the red and green pigment genes diverged from each other much more recently. Nucleotide substitutions in the rhodopsin genes are best explained by the neutral theory of molecular evolution. However, important functional adaptations seem to have occurred twice during the evolution of the color pigment genes in humans: first, to the common ancestor of the three color pigment genes after its divergence from the ancestor of the rhodopsin gene and, second, to the ancestor of the red pigment gene after its divergence from that of the green pigment gene.     

Molecular evolution of the murine tspy genes

Molecular aspects of murine evolution were studied by sequencing, and subsequently comparing, introns of the Y-chromosomal tspy genes from Apodemus agrarius, A. sylvaticus, A. flavicollis, Mus platythrix (subgenus Pyromys), M. booduga (subgenus Leggada), and from species of the subgenus Mus, including M. cervicolor, M. macedonicus and M. spretus. Estimates of nucleotide substitution rates in these lineages were in perfect agreement with phylogenetic data previously published by She et al. (1990), Catzeflis et al. (1992; 1993), and Lyon et al. (1996). The only exception was provided by a comparatively late divergence of M. spretus and M. macedonicus. Our data also suggest that M. booduga diverged from the subgenus Mus about 3 Myr ago.     

Molecular evolution of plant immune system genes

Molecular population genetic studies are providing new perspectives on the evolution of genes that confer resistance to pathogens and herbivores. Here, we compare the evolutionary history of different components of the defense response (detection, signaling and response) and of genes with parallel function in plants and Drosophila. A review of the literature indicates that the dominant form of selection acting on defense genes (balancing, positive and purifying) differs among components of defense. Sampling of particular classes of genes and genes from non-model organisms, however, remains limited. Future studies combining molecular evolutionary analyses with ecological genetic and functional analyses should better reveal how natural selection has shaped defense gene evolution.     

Molecular phylogenies and evolution of crt genes in algae

The carotenoids constitute the most widespread class of pigments in nature. Most previous work has concentrated on the identification and characterization of their chemical physical properties and bioavailability. In recent years, significant amounts of research have been conducted in an attempt to analyze the genes and the molecular regulation of the genes involved in the biosynthesis of carotenoids. However, it is important not to lose sight of the early evolution of carotenoid biosynthesis. One of the major obstacles in understanding the evolution of the respective enzymes and their patterns of selection is a lack of a well-supported phylogenic analysis. In the present research, a major long-term objective was to provide a clearer picture of the evolutionary history of genes, together with an evaluation of the patterns of selection in algae. These phylogenies will be important in studies characterizing the evolution of algae. The gene sequences of the enzymes involved in the major steps of the carotenoid biosynthetic pathway in algae (cyanobacteria, rhofophyta, chlorophyta) have been analyzed. Phylogenetic relationships among protein-coding DNA sequences were reconstructed by neighbor-joining (NJ) analysis for the respective carotenoid biosynthetic pathway genes (crt) in algae. The analysis also contains an estimation of the rate of nonsynonymous nucleotide substitutions per nonsynonymous site (d(N)), synonymous nucleotide substitution per synonymous site (d(S)), and the ratio of nonsynonmous (d(N)/d(S)) for the test of selection patterns. The phylogenetic trees show that the taxa of some genera have a closer evolutionary relationship with other genera in some gene sequences, which suggests a common ancient origin and that lateral gene transfer has occurred among unrelated genera. The d(N) values of crt genes in the early pathway are relatively low, while those of the following steps are slightly higher, while the d(N) values of crt genes in chlorophyta are higher than those in cyanobacteria. Most of the d(N)/d(S) values exceed 1. The phylogenetic analysis revealed that lateral gene transfer may have taken place across algal genomes and the d(N) values suggest that most of the early crt genes are well conserved compared to the later crt genes. Furthermore, d(N) values also revealed that the crt genes of chlorophyta are more evolutionary than cyanobacteria. The amino acids' changes are mostly adaptive evolution under the influence of positive diversity selection.     

Molecular evolution of group II phospholipases A2

The nucleotide sequences of 13 cDNAs encoding group II phospholipases A₂ (PLA₂ ^S), which are from viperidae snake venoms and from mammalian sources, were aligned and analyzed by phylogenetic trees constructed using various components of the sequences. The evolutionary trees derived from the combined sequences of the untranslated (5 and 3) region and the signal peptide region of cDNAs were in accord with the consequences from taxonomy. In contrast, the evolutionary trees from the mature protein-coding region sequences of cDNAs and from the amino acid sequences showed random patterns. These observations indicated that the mature protein-coding region has evolved through a process differently from the untranslated and signal peptide regions. The trees built from the nucleotide differences at each of three positions of codons in the mature protein-coding region suggested that snakevenom-gland PLA₂ genes have evolved via a process different from mammalian PLA₂ genes. The occurrence of accelerated evolution has been recently discovered in Trimeresurus flavoviridis venom-gland group II PLA₂ isozyme genes (Nakashima et al. 1993, Proc Natl Acad Sci USA 90:5964–5968), so the present phylogenetic analysis together with the estimation of nucleotide divergence of cDNAs provides further evidence that snakevenom-group II PLA₂ isozyme genes have evolved by accelerated evolution to gain diverse physiological activities. Correspondence to: M. Ohno