Mhc class II B gene evolution in East African cichlid fishes

 A distinctive feature of essential major histocompatibility complex (Mhc) loci is their polymorphism characterized by large genetic distances between alleles and long persistence times of allelic lineages. Since the lineages often span several successive speciations, we investigated the behavior of the Mhc alleles during or close to the speciation phase. We sequenced exon 2 of the class II B locus 4 from 232 East African cichlid fishes representing 32 related species. The divergence times of the (sub)species ranged from 6000 to 8.4 million years. Two types of evolutionary analysis were used to elucidate the pattern of exon 2 sequence divergence. First, phylogenetic methods were applied to reconstruct the most likely evolutionary pathways leading from the last common ancestor of the set to the extant sequences, and to assess the probable mechanisms involved in allelic diversification. Second, pairwise comparisons of sequences were carried out to detect differences seemingly incompatible with origin by nonparallel point mutations. The analysis revealed point mutations to be the most important mechanism behind allelic divergences, with recombination playing only an auxiliary part. Comparison of sequences from related species revealed evidence of random allelic (lineage) losses apparently associated with speciation. Sharing of identical alleles could be demonstrated between species that diverged 2 million years ago. The phylogeny of the exon was incongruent with that of the flanking introns, indicating either a high degree of convergent evolution at the peptide-binding region-encoding sites, or intron homogenization. Received: 6 December 1999 / Revised: 15 February 2000     

Evolution of Duplicated <Emphasis Type="BoldItalic">reggie</Emphasis> Genes in Zebrafish and Goldfish

Invertebrates, tetrapod vertebrates, and fish might be expected to differ in their number of gene copies, possibly due the occurrence of genome duplication events during animal evolution. Reggie (flotillin) genes code for membrane-associated proteins involved in growth signaling in developing and regenerating axons. Until now, there appeared to be only two reggie genes in fruitflies, mammals, and fish. The aim of this research was to search for additional copies of reggie genes in fishes, since a genome duplication might have increased the gene copy number in this group. We report the presence of up to four distinct reggie genes (two reggie-1 and two reggie-2 genes) in the genomes of zebrafish and goldfish. Phylogenetic analyses show that the zebrafish and goldfish sequence pairs are orthologous, and that the additional copies could have arisen through a genome duplication in a common ancestor of bony fish. The presence of novel reggie mRNAs in fish embryos indicates that the newly discovered gene copies are transcribed and possibly expressed in the developing and regenerating nervous system. The intron/exon boundaries of the new fish genes characterized here correspond with those of human genes, both in location and phase. An evolutionary scenario for the evolution of reggie intron-exon structure, where loss of introns appears to be a distinctive trait in invertebrate reggie genes, is presented. Received: 24 January 2001 / Accepted: 27 July 2001     

Evolutionary analysis of the APA genes in the <Emphasis Type="Italic">Phaseolus</Emphasis> genus: wild and cultivated bean species as sources of lectin-related resistance factors?

The APA (Arcelin/Phytohemagglutinin/α-Amylase inhibitor) gene family is composed of various members, present in Phaseolus species and coding for lectin and lectin-related seed proteins having the double role of storage and defense proteins. Here members of the APA family have been identified by immunological, functional, and molecular analyses and representative genes were sequenced in nine wild species of Phaseolus. All taxa possessed at least one member of the true lectin gene. No arcelin type sequences have been isolated from the species examined. Among the wild species studied, only P. costaricensis contained an α-amylase inhibitor (α-AI). In addition P. augusti, P. maculatus, P. microcarpus, and P. oligospermus showed the presence of the lectin-related α-amylase inhibitor-like (AIL) genes and α-AI activity. Data from Southern blot analysis indicated the presence of only one lectin gene in P. parvulus and P. filiformis, while an extensive gene duplication of the APA locus was found in the other Phaseolus species. Phylogenetic analysis carried out on the nucleotide sequences showed the existence of two main clusters and clearly indicated that lectin-related genes originated from a paralogous duplication event preceding the development of the ancestor to the Phaseolus genus. The finding of detectable α-AI activity in species containing AIL genes suggests that exploiting APA genes variability in the Phaseolus genus may represent a valuable tool to find new members that may have acquired insecticidal activities. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Platypus globin genes and flanking loci suggest a new insertional model for beta-globin evolution in birds and mammals

Background  

Vertebrate alpha (α)- and beta (β)-globin gene families exemplify the way in which genomes evolve to produce functional complexity. From tandem duplication of a single globin locus, the α- and β-globin clusters expanded, and then were separated onto different chromosomes. The previous finding of a fossil β-globin gene (ω) in the marsupial α-cluster, however, suggested that duplication of the α-β cluster onto two chromosomes, followed by lineage-specific gene loss and duplication, produced paralogous α- and β-globin clusters in birds and mammals. Here we analyse genomic data from an egg-laying monotreme mammal, the platypus (Ornithorhynchus anatinus), to explore haemoglobin evolution at the stem of the mammalian radiation.     

Analysis of Fungal Genomes Reveals Commonalities of Intron Gain or Loss and Functions in Intron-Poor Species

Previous evolutionary reconstructions have concluded that early eukaryotic ancestors including both the last common ancestor of eukaryotes and of all fungi had intron-rich genomes. By contrast, some extant eukaryotes have few introns, underscoring the complex histories of intron–exon structures, and raising the question as to why these few introns are retained. Here, we have used recently available fungal genomes to address a variety of questions related to intron evolution. Evolutionary reconstruction of intron presence and absence using 263 diverse fungal species supports the idea that massive intron reduction through intron loss has occurred in multiple clades. The intron densities estimated in various fungal ancestors differ from zero to 7.6 introns per 1 kb of protein-coding sequence. Massive intron loss has occurred not only in microsporidian parasites and saccharomycetous yeasts, but also in diverse smuts and allies. To investigate the roles of the remaining introns in highly-reduced species, we have searched for their special characteristics in eight intron-poor fungi. Notably, the introns of ribosome-associated genes RPL7 and NOG2 have conserved positions; both intron-containing genes encoding snoRNAs. Furthermore, both the proteins and snoRNAs are involved in ribosome biogenesis, suggesting that the expression of the protein-coding genes and noncoding snoRNAs may be functionally coordinated. Indeed, these introns are also conserved in three-quarters of fungi species. Our study shows that fungal introns have a complex evolutionary history and underappreciated roles in gene expression.     

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.     

Phylogenetic Relationships Within the Genus Equus and the Evolution of α and θ Globin Genes

Sequences of the α1, α2 and θ globin genes from six equid species have been determined to investigate relationships within the genus Equus. Analyses using standard phylogenetic methods, or an approach designed to account for the effects of gene conversion between the α genes, gave broadly similar results and show that the horses diverged from the zebra/ass ancestor ∼2.4 million years ago and that the zebra and ass species arose in a rapid radiation ∼0.9 million years ago. These results from the α genes are corroborated by θ gene data and are in contrast to mitochondrial DNA studies of the phylogeny of this genus, which suggest a more gradual set of speciation events. Received: 22 April 1997 / Accepted: 20 July 1998     

Characterization of Two Alcohol Dehydrogenase (Adh) Loci from the Olive Fruit Fly, Bactrocera (Dacus) oleae and Implications for Adh Duplication in Dipteran Insects

We report the cloning and structural characterization of two Adh loci of the olive fruit fly, Bactrocera oleae. Each of the two genes, named Adh1 and Adh2, consists of three exons and two introns for a total length of 1981 and 988 nucleotides, respectively. Their deduced amino acid sequences of 257 and 258 residues exhibit a 77% identity and display the characteristics of the insect ADH enzymes, which belong to the short-chain dehydrogenases/reductases family. The Adh genes of B. oleae are compared to the two genes of the Mediterranean fly, Ceratitis capitata, the only other species of the Tephritidae family in which the Adh genes have been studied. On the basis of amino acid divergence the four genes form two clusters each containing one gene from each species, as expected if there was one duplication event before speciation. On the basis of nucleotide sequence the four sequences form two clusters each containing the two sequences from the same species, as expected if there was a separate duplication event in each species. To help decide between the two alternatives, we compared at both the amino acid and DNA level the Adh genes of five Drosophila species that are known to carry two such genes and observed that, with only one exception at the amino acid level, conspecific loci cluster together. We conclude that the information we have at present does not allow a firm choice between the hypothesis of a single duplication event that occurred before the split of Bactrocera and Ceratitis from their common ancestor and the hypothesis of two independent duplication events, one in each of the two genera. Received: 30 May 2000 / Accepted: 17 August 2000     

Investigation of the RH locus in gorillas and chimpanzees

The human Rh blood-group system is encoded by two homologous genes,RhD andRhCE. TheRH genes in gorillas and chimpanzees were investigated to delineate the phylogeny of the humanRH genes. Southern blot analysis with an exon 7-specific probe suggested that gorillas have more than twoRH genes, as has recently been reported for chimpanzees. Exon 7 was well conserved between humans, gorillas, and chimpanzees, although the exon 7 nucleotide sequences from gorillas were more similar to the humanD gene, whereas the nucleotide sequences of this exon in chimpanzees were more similar to the humanCE gene. The intron between exon 4 and exon 5 is polymorphic and can be used to distinguish the humanD gene from theCE gene. Nucleotide sequencing revealed that the basis for the intron polymorphism is anAlu element inCE which is not present in theD gene. Examination of gorilla and chimpanzee genomic DNA for this intron polymorphism demonstrated that theD intron was present in all the chimpanzees and in all but one gorilla. TheCE intron was found in three of six gorillas, but in none of the seven chimpanzees. Sequence data suggested that theAlu element might have previously been present in the chimpanzeeRH genes but was eliminated by excision or recombination. Conservation of theRhD gene was also apparent from the complete identity between the 3′-noncoding region of the human D cDNA and a gorilla genomic clone, including anAlu element which is present in both species. The data suggest that at least twoRH genes were present in a common ancestor of humans, chimpanzees, and gorillas, and that additionalRH gene duplication has taken place in gorillas and chimpanzees. TheRhCE gene appears to have diverged more thanRhD among primates. In addition, theRhD gene deletion associated with the Rh-negative phenotype in humans seems to have occurred after speciation. Correspondence to: C.M. Westhoff     

Rapid lineage-specific diversification of the mast cell chymase locus during mammalian evolution

Serine proteases constitute the major protein granule content of cells of several hematopoietic cell lineages. A subgroup of these proteases, including the mast cell chymases, neutrophil cathepsin G, and T cell granzymes B to F and N, are in all investigated mammals encoded in one locus, the chymase locus. It is interesting to note that this locus has diversified greatly during the last 95 Myr of mammalian evolution. This divergence is exemplified by the presence of Mcpt8-related genes and multiple β-chymases in the mouse and rat, which lack direct counterparts in primates and in seven functional granzyme genes in the mouse where the human locus has only two. To study the expansion of the locus during rodent evolution and to better understand the evolutionary origin of β-chymases and the Mcpt8-family, we have performed a detailed analysis of the chymase locus of four mammalian species, i.e., human, dog, mouse, and rat. As a result, we report here a second chymase-like gene in dog, Cma2, which clusters with β-chymases in phylogenetic analyses. This finding supports a duplication of the common ancestor for α- and β-chymases before the major radiation of placental mammals, and a loss of the ancestral β-chymase gene sometime during primate evolution. Moreover, we show that in the rat, the Mcpt8-family diversified relatively recently together with sequences related to the β-chymase Mcpt2. Eight novel genes were identified in the duplication region, four of which are predicted to be functional. Duplications of rat granzyme B- and C-like sequences occurred seemingly independently within a similar time frame, but did not give rise to functional genes. Due to the duplications in rat and deletions in the carnivore/primate lineage, the rat chymase locus is approximately 15 and 9 times larger than its counterparts in dog and human, respectively. These findings illustrate the importance of gene duplications in conferring rapid changes in mammalian genomes.     

Origin of new genes: evidence from experimental and computational analyses

Exon shuffling is an essential molecular mechanism for the formation of new genes. Many cases of exon shuffling have been reported in vertebrate genes. These discoveries revealed the importance of exon shuffling in the origin of new genes. However, only a few cases of exon shuffling were reported from plants and invertebrates, which gave rise to the assertion that the intron-mediated recombination mechanism originated very recently. We focused on the origin of new genes by exon shuffling and retroposition. We will first summarize our experimental work, which revealed four new genes in Drosophila, plants, and humans. These genes are 10⁶ to 10⁸ million years old. The recency of these genes allows us to directly examine the origin and evolution of genes in detail. These observations show firstly the importance of exon shuffling and retroposition in the rapid creation of new gene structures. They also show that the resultant chimerical structures appearing as mosaic proteins or as retroposed coding structures with novel regulatory systems, often confer novel functions. Furthermore, these newly created genes appear to have been governed by positive Darwinian selection throughout their history, with rapid changes of amino acid sequence and gene structure in very short periods of evolution. We further analyzed the distribution of intron phases in three non-vertebrate species, Drosophila melanogaster, Caenorhabditis elegans, and Arabidosis thaliana, as inferred from their genome sequences. As in the case of vertebrate genes, we found that intron phases in these species are unevenly distributed with an excess of phase zero introns and a significant excess of symmetric exons. Both findings are consistent with the requirements for the molecular process of exon shuffling. Thus, these non-vertebrate genomes may have also been strongly impacted by exon shuffling in general.     

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.     

h-TBP: an approach based on intron-length polymorphism for the rapid isolation and characterization of the multiple members of the β-tubulin gene family in <Emphasis Type="Italic">Camelina sativa</Emphasis> (L.) Crantz

We have developed a new version of the cTBP (combinatorial tubulin-based polymorphism) method, a previously described approach based on intron-length polymorphism (ILP), to rapidly characterize the β-tubulin gene family of Camelina sativa (L.) Crantz, a plant species of importance for oil production but still largely unexplored at genomic level. The method, named h-TBP, allows the rapid cloning of the β-tubulin genomic sequences that encompass the two introns, invariantly present at fixed positions within the coding region of the vast majority of the plant species. The β-tubulin sequences cloned by h-TBP also comprise part of exon1 and exon3 and the whole sequence of exon2. The h-TBP method has then been used to isolate, clone and characterize the β-tubulin gene family of C. sativa, composed of at least 20 different β-tubulin isotypes, named CsTUB1 through CsTUB20. The relatively high number of β-tubulin genes has been further substantiated by Southern-blot analysis. Comparison of the β-tubulin exon sequences of C. sativa with those of Arabidopsis thaliana, the closest relative among crucifers, defines distinct groups of putative orthologous genes, identified by a UPMGA cluster analysis. Analysis of the C. sativa β-tubulin intron sequences reveals some molecular features that can provide the first hints for the understanding of intron plasticity and evolution. From a more immediate perspective, these data provide the first substantial contribution to the characterization of the largely unexplored genome of C. sativa, and the tools for assisting programmes of breeding and selection of the most productive plants.     

OrthoParaMap: Distinguishing orthologs from paralogs by integrating comparative genome data and gene phylogenies

Background  

In eukaryotic genomes, most genes are members of gene families. When comparing genes from two species, therefore, most genes in one species will be homologous to multiple genes in the second. This often makes it difficult to distinguish orthologs (separated through speciation) from paralogs (separated by other types of gene duplication). Combining phylogenetic relationships and genomic position in both genomes helps to distinguish between these scenarios. This kind of comparison can also help to describe how gene families have evolved within a single genome that has undergone polyploidy or other large-scale duplications, as in the case of Arabidopsis thaliana – and probably most plant genomes.     

Evolutionary History and Mode of the <Emphasis Type="Italic">amylase</Emphasis> Multigene Family in <Emphasis Type="Italic">Drosophila</Emphasis>

Previous studies indicate that the tandemly repeated members of the amylase (Amy) gene family evolved in a concerted manner in the melanogaster subgroup and in some other species. In this paper, we analyzed all of the 49 active and complete Amy gene sequences in Drosophila, mostly from subgenus Sophophora. Phylogenetic analysis indicated that the two types of diverged Amy genes in the Drosophila montium subgroup and Drosophila ananassae, which are located in distant chromosomal regions from each other, originated independently in different evolutionary lineages of the melanogaster group after the split of the obscura and melanogaster groups. One of the two clusters was lost after duplication in the melanogaster subgroup. Given the time, 24.9 mya, of divergence between the obscura and the melanogaster groups (Russo et al. 1995), the two duplication events were estimated to occur at about 13.96 ± 1.93 and 12.38 ± 1.76 mya in the montium subgroup and D. ananassae, respectively. An accelerated rate of amino acid changes was not observed in either lineage after these gene duplications. However, the G+C contents at the third codon positions (GC3) decreased significantly along one of the two Amy clusters both in the montium subgroup and in D. ananassae right after gene duplication. Furthermore, one of the two types of the Amy genes with a lower GC3 content has lost a specific regulatory element within the montium subgroup species and D. ananassae. While the tandemly repeated members evolved in a concerted manner, the two types of diverged Amy genes in Drosophila experienced frequent gene duplication, gene loss, and divergent evolution following the model of a birth-and-death process.     

Adaptive Evolution and Frequent Gene Conversion in the Brain Expressed X-Linked Gene Family in Mammals

This work examines the molecular evolution of a brain-expressed X-linked gene family in the mammalian genomes of human, chimp, macaque, mouse, rat, dog, and cow. The gene structures are well conserved across family members and among the mammals in that all five members have three exons with the first two exons untranslated. Furthermore, the five members are arranged tandemly on chromosome X with Bex5, Bex1, Bex2 on the negative strand and Bex4, Bex3 on the positive strand, and this physical arrangement remains conserved among species. Sequence analyses indicate that gene conversion has been frequent and ongoing among Bex1-4, occurring in multiple species independently. All gene conversions in different species between Bex1 and Bex4, and between Bex2 and Bex3, appear to be limited to the upstream regions of the third exon, whereas the gene conversions occurred independently in different species between Bex1 and Bex2 and cover only the third exon. Bex5 appears to have little exchange of genetic information with other members, possibly due to its distance from other members. The GC content decreases from 5′-UTR, intron 1, intron 2, coding region, to 3′-UTR, reflecting faithfully the frequency of gene conversion in different regions of the Bex genes. Sequence analyses also suggest that both relaxed selective constraint and positive selection have acted on the Bex members after duplication. In particular, Bex3 shows strong evidence of positive selection and seems to have evolved a new gene function after gene duplication.     

The Mutually Exclusive Flip and Flop Exons of AMPA Receptor Genes Were Derived from an Intragenic Duplication in the Vertebrate Lineage

The incomplete correlation between the organismal complexities and the number of genes among eukaryotic organisms can be partially explained by multiple protein products of a gene created by alternative splicing. One type of alternative splicing involves alternative selection of mutually exclusive exons and creates protein products with substitution of one segment of the amino acid sequence for another. To elucidate the evolution of the mutually exclusive 115-bp exons, designated flip and flop, of vertebrate AMPA receptor genes, the gene structures of chordate (tunicate, cephalochordate, and vertebrate) and protostome (Drosophila and Caenorhabditis elegans) AMPA receptor subunits were compared. Phylogenetic analysis supports that the vertebrate flip and flop exons evolved from a common sequence. Flip and flop exons exist in all vertebrate AMPA receptor genes but only one 115-bp exon is present in the genes of tunicates and cephalochordates, suggesting that the exon duplication event occurred at the ancestral vertebrate AMPA receptor gene after the separation of vertebrates from primitive chordates. The structures of animal AMPA receptor genes also suggest that an intron insertion to separate the primordial flip/flop exon from the M4-coding exon occurred before the exon duplication event and probably at the chordate lineage. [Reviewing Editor: Dr. Manyuan Long]     

A Nested Alpha-Amylase Gene in <Emphasis Type="Italic">Drosophila ananassae</Emphasis>

The amylase gene family of Drosophila ananassae consists in seven copies, scattered on several chromosomal arms. We have evidenced that a member of the family, Amy35, lies within an intron of a gene homologous to the CG14696 gene of D. melanogaster. This nested arrangement seems restricted to the D. ananassae subgroup. The nested and the nest genes are encoded on opposite strands. Both are actively transcribed in the midgut at the same time, raising the possibility of interference between their mRNAs. Our data also help to elucidate the history of the Amy family, suggesting that Amy35 arose by duplication and translocation from another ancestral locus, into a formerly short intron, in an ancestor of the subgroup.     

Molecular Aspects of Intron Evolution in Dynein Encoding Mega-Genes on The Heterochromatic Y Chromosome of Drosophila sp.

Fertility genes on the heterochromatic Y chromosome of various Drosophilaspecies are unique for several reasons. Most of them are megabase-sized. Their expression is restricted to premeiotic spermatocytes and often associated with unfolding of huge species-specific lampbrush loops. Molecular analysis of the orthologous dynein genes Dhc-Yh3, DhDhc7(Y)and DeDhc7(Y)on the Y chromosome of the three species D. melanogaster, D. hydeiand D. eohydei, respectively, revealed that the megabase gene size as well as the species-specific morphology of the corresponding lampbrush loops kl-5, Threadsand diffuse loopsresult from huge introns and their specific sequence composition, whereas the majority of all 20 introns in each of the three genes is in a size of 45–72 bp. The loop-specifying introns are extreme exceptions due to extended assemblies of degenerated transposable elements and/or large clusters of satellite DNAs. Here we use sequence information from the complete intron sets of three orthologous Y chromosomal dynein genes to deduce a scenario for an evolutionary pathway leading to the megabase-sized genes on the heterochromatic Y chromosome of Drosophila. The obvious bias between very small and species-specific mega introns is explained as the result of an autocatalytic mode of intron growth. An initial coincidental hit by a single transposable element extends the size of a 50 bp intron for about two orders of magnitude and determines it for preferential extension by similar insertion events. This phase of continuous moderate growth is followed by rapid size enlargements by repeating amplifications generating extended clusters of satellite DNA. Size control by recombination, on the other hand, is suppressed in Drosophilamales by achiasmatic meiosis. This revised version was published online in July 2006 with corrections to the Cover Date.