Dynamic evolution of the CYP2ABFGST gene cluster in primates

We have annotated and analyzed sequences from the CYP2ABFGST gene cluster for the chimpanzee and the Rhesus macaque, and compared them in detail to the complete map of the syntenic human cluster. While the macaque cluster has accumulated substantial differences in gene number and organization relative to the other two species, the chimpanzee cluster has remained surprisingly similar to the human cluster, given the evolutionary volatility of this gene family. The comparison will allow us to further refine our model of the evolution of this gene family in mammals.     

The amplification and evolution of orthologous 22-kDa α-prolamin tandemly arrayed genes in <Emphasis Type="Italic">coix</Emphasis>, sorghum and maize genomes

Tandemly arrayed genes (TAGs) account for about one-third of the duplicated genes in eukaryotic genomes. They provide raw genetic material for biological evolution, and play important roles in genome evolution. The 22-kDa prolamin genes in cereal genomes represent typical TAG organization, and provide the good material to investigate gene amplification of TAGs in closely related grass genomes. Here, we isolated and sequenced the Coix 22-kDa prolamin (coixin) gene cluster (283 kb), and carried out a comparative analysis with orthologous 22-kDa prolamin gene clusters from maize and sorghum. The 22-kDa prolamin gene clusters descended from orthologous ancestor genes, but underwent independent gene amplification paths after the separation of these species, therefore varied dramatically in sequence and organization. Our analysis indicated that the gene amplification model of 22-kDa prolamin gene clusters can be divided into three major stages. In the first stage, rare gene duplications occurred from the ancestor gene copy accidentally. In the second stage, rounds of gene amplification occurred by unequal crossing over to form tandem gene array(s). In the third stage, gene array was further diverged by other genomic activities, such as transposon insertions, segmental rearrangements, etc. Unlike their highly conserved sequences, the amplified 22-kDa prolamin genes diverged rapidly at their expression capacities and expression levels. Such processes had no apparent correlation to age or order of amplified genes within TAG cluster, suggesting a fast evolving nature of TAGs after gene amplification. These results provided insights into the amplification and evolution of TAG families in grasses.     

Origin and Evolution of the Gene Family of Proteinaceous Pheromones,the Exocrine Gland-Secreting Peptides,in Rodents

The exocrine-gland secreting peptide (ESP)gene family encodes proteinaceous pheromones that are recognized by the vomeronasal organ in mice. For example, ESP1 is a male pheromone secreted in tear fluid that regulates socio-sexual behavior, and ESP22 is a juvenile pheromone that suppresses adult sexual behavior. The family consists of multiple genes and has been identified only in mouse and rat genomes. The coding region of a mouse ESP gene is separated into two exons, each encoding signal and mature sequences. Here, we report the origin and evolution of the ESP gene family. ESP genes were found only in the Muridea and Cricetidae families of rodents, suggesting a recent origin of ESP genes in the common ancestor of murids and cricetids. ESP genes show a great diversity in number, length, and sequence among different species as well as mouse strains. Some ESPs in rats and golden hamsters are expressed in the lacrimal gland and the salivary gland. We also found that a mature sequence of an ESP gene showed overall sequence similarity to the α-globin gene. The ancestral ESP gene seems to be generated by recombination of a retrotransposed α-globin gene with the signal-encoding exon of the CRISP2 gene located adjacent to the ESP gene cluster. This study provides an intriguing example of molecular tinkering in rapidly evolving species-specific proteinaceous pheromone genes.     

Evolutionarily plastic regions at human 3p21.3 coincide with tumor breakpoints identified by the "elimination test"

We have previously found with the microcell hybrid-based "elimination test" that human chromosome 3 transferred into murine or human tumor cells regularly lost certain 3p regions during tumor growth in SCID mice. The most common eliminated region, CER1, is approximately 2.4 Mb at 3p21.3. CER1 breakpoints were clustered in approximately 200-kb regions at both telomeric and centromeric borders. We have also shown, earlier, that tumor-related deletions often coincide with human/mouse synteny breakpoints on 3p12-p22. Here we describe the results of a comparative genomic analysis on the CER1 region in Caenorhabditis elegans, Drosophila melanogaster, Fugu rubripes, Gallus gallus, Mus musculus, Rattus norvegicus, and Canis familiaris. First, four independent synteny breaks were found within the CER1 telomeric breakpoint cluster region, comparing human, dog, and chicken genomes, and two independent synteny breaks within the CER1 centromeric breakpoint cluster region, comparing human, mouse, and chicken genomes, suggesting a nonrandom involvement of tumor breakpoint regions in chromosome evolution. Second, both CER1 breakpoint cluster regions show recent tandem duplications (seven Zn finger protein family genes at the telomeric and eight chemokine receptor genes at the centromeric side). Finally, all genes from these regions underwent horizontal evolution in mammals, with formation of new genes and expansion of gene families, which were displayed in the human genome as tandem gene duplications and pseudogene insertions. In contrast the CER1 middle region contained evolutionarily well-conserved solitary genes and a minimal amount of retroposed genes. The coincidence of evolutionary plasticity with CER1 breakpoints may suggest that regional structural instability is expressed in both evolutionary and cancer-associated chromosome rearrangements.     

Genomic organization and evolutionary analysis of <Emphasis Type="Italic">Ly49</Emphasis> genes encoding the rodent natural killer cell receptors: rapid evolution by repeated gene duplication

Ly49 genes regulate the cytotoxic activity of natural killer (NK) cells in rodents and provide important protection against virus-infected or tumor cells. About 15 Ly49 genes have been identified in mice, but only a few genes have been reported to date in rats. Here we studied all Ly49 genes in the entire rat genome sequence and identified 17 putative functional and 16 putative non-functional genes together with their genomic locations in a 1.8-Mb region of chromosome 4. Phylogenetic analysis of these genes indicated that the Ly49 gene family expanded rapidly in recent years, and this expansion was mediated by both tandem and genomic block duplication. The joint phylogenetic analysis of mouse and rat genes suggested that the most recent common ancestor of the two species had at least several Ly49 genes, but that the majority of current duplicate genes were generated after divergence of the two species. In both species Ly49 genes are apparently subject to birth-and-death evolution, but the birth and death rates of Ly49 genes are higher in rats than in mice. The rate of gene expansion in the Ly49 gene family in rats is one of the highest among all mammalian multigene families so far studied. The biochemical function of Ly49 genes is essentially the same as that of KIR genes in primates, but the molecular structures of the two groups of NK cell receptors are very different. A hypothesis was presented to explain the origin of the differential use of Ly49 and KIR genes in rodents and primates.     

The recent origin of the Sdic gene cluster in the melanogaster subgroup

Gene families are composed of closely related genes and are an important part of eukaryotic genomes. In the proximal region of the X chromosome of Drosophila melanogaster there is a cluster of four tandem Sdic genes, located between the gene Cdic and the gene AnnX. Sdic is a chimeric gene that encodes a novel protein with sperm-specific expression. It had been hypothesized that the Sdic gene cluster was formed after the split of D. melanogaster and D. simulans. To study the evolution of this cluster, the sequence of this region was studied in several Drosophilidae species. In all species analyzed, Sdic genes are absent and AnnX and Cdic are adjacent to each other. The results allowed the inference of the ancestral situation and the reconstruction of the evolution of the cluster, and confirm that the Sdic cluster was indeed formed in the lineage that gave rise to D. melanogaster, being one of the youngest gene clusters known. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Features of cytochrome P450 evolution]

It is shown that the process of mutation in the CYP2 family of the superfamily of P450 cytochromes is species-specific (man, rat, and mouse). It is also shown that, within one species (rat), different families (CYP2 and CYP11) have different mutation spectra, indicating a high specificity of the mutation process for the families of cytochrome genes. A similar specificity was demonstrated for five families (CYP1, CYP2, CYP6, CYP7, CYP11) as compared with globins and prions. The analysis of the evolutionary mutation pattern, and the pattern of pseudogenes and damaged alleles of the CYP21 family (found in patients with congenital adrenal hyperplasia) does not confirm the widely accepted hypothesis that mutations arising in pseudogenes are transduced to normal alleles of the CYP21 gene through gene conversion.     

The P450 superfamily: update on new sequences, gene mapping, and recommended nomenclature

We provide here a list of 154 P450 genes and seven putative pseudogenes that have been characterized as of October 20, 1990. These genes have been described in a total of 23 eukaryotes (including nine mammalian and one plant species) and six prokaryotes. Of 27 gene families so far described, 10 exist in all mammals. These 10 families comprise 18 subfamilies, of which 16 and 14 have been mapped in the human and mouse genomes, respectively; to date, each subfamily appears to represent a cluster of tightly linked genes. We propose here a modest revision of the initially proposed (Nebert et al., DNA 6, 1-11, 1987) and updated (Nebert et al., DNA 8, 1-13, 1989) nomenclature system based on evolution of the superfamily. For the gene we recommend that the italicized root symbol CYP for human (Cyp for mouse), representing cytochrome P450, be followed by an Arabic number denoting the family, a letter designating the subfamily (when two or more exist), and an Arabic numeral representing the individual gene within the subfamily. A hyphen should precede the final number in mouse genes. We suggest that the human nomenclature system be used for other species. This system is consistent with our earlier proposed nomenclature for P450 of all eukaryotes and prokaryotes, except that we are discouraging the future use of cumbersome Roman numerals.     

Why mice have lost genes for COL21A1, STK17A, GPR145 and AHRI: evidence for gene deletion at evolutionary breakpoints in the rodent lineage

The mouse genome has undergone extensive chromosome rearrangement relative to the human genome since these species last shared a common ancestor. One possible consequence of these rearrangements is the deletion of genes that are located within evolutionary breakpoint regions. In this article, we present evidence of four human genes (COL21A1, STK17A, GPR145 and ARHI) that are located in regions corresponding to evolutionary breakpoints in rodents and lack mouse and rat orthologues. We propose that "evolutionary breakpoint-associated gene deletion" is an unexpected consequence of evolutionary chromosome rearrangement, and we describe a novel mechanism through which genes can be lost during evolution.     

Patterns of gene duplication in the plant SKP1 gene family in angiosperms: evidence for multiple mechanisms of rapid gene birth

Gene duplication plays important roles in organismal evolution, because duplicate genes provide raw materials for the evolution of mechanisms controlling physiological and/or morphological novelties. Gene duplication can occur via several mechanisms, including segmental duplication, tandem duplication and retroposition. Although segmental and tandem duplications have been found to be important for the expansion of a number of multigene families, the contribution of retroposition is not clear. Here we show that plant SKP1 genes have evolved by multiple duplication events from a single ancestral copy in the most recent common ancestor (MRCA) of eudicots and monocots, resulting in 19 ASK (Arabidopsis SKP1-like) and 28 OSK (Oryza SKP1-like) genes. The estimated birth rates are more than ten times the average rate of gene duplication, and are even higher than that of other rapidly duplicating plant genes, such as type I MADS box genes, R genes, and genes encoding receptor-like kinases. Further analyses suggest that a relatively large proportion of the duplication events may be explained by tandem duplication, but few, if any, are likely to be due to segmental duplication. In addition, by mapping the gain/loss of a specific intron on gene phylogenies, and by searching for the features that characterize retrogenes/retrosequences, we show that retroposition is an important mechanism for expansion of the plant SKP1 gene family. Specifically, we propose that two and three ancient retroposition events occurred in lineages leading to Arabidopsis and rice, respectively, followed by repeated tandem duplications and chromosome rearrangements. Our study represents a thorough investigation showing that retroposition can play an important role in the evolution of a plant gene family whose members do not encode mobile elements.     

Differential evolution of voltage-gated sodium channels in tetrapods and teleost fishes

The voltage-gated sodium channel (SCN) alpha subunits are large proteins with central roles in the generation of action potentials. They consist of approximately 2,000 amino acids encoded by 24-27 exons. Previous evolutionary studies have been unable to reconcile the proposed gene duplication schemes with the species distribution and molecular phylogeny of the genes. We have carefully annotated the complete SCN gene sequences, correcting numerous database errors, for a broad range of vertebrate species and analyzed their phylogenetic relationships. We have also compared the chromosomal positions of the SCN genes relative to adjacent gene families. Our studies show that the ancestor of the vertebrates probably had a single sodium channel gene with two characteristic AT-AC introns, the second of which is unique to vertebrate SCN genes. This ancestral gene, located close to a HOX gene cluster, was quadrupled along with HOX in the two rounds of basal vertebrate tetraploidizations to generate the ancestors of the four channels SCN1A, SCN4A, SCN5A, and SCN8A. The third tetraploidization in the teleost fish ancestor doubled this set of genes and all eight are still present in at least three of four investigated teleost fish genomes. In tetrapods, the gene family expanded by local duplications before the radiation of amniotes, generating the cluster SCN5A, SCN10A, and SCN11A on one chromosome and the cluster SCN1A, SCN2A, SCN3A, and SCN9A on a different chromosome. In eutherian mammals, a tenth gene, SCN7A, arose in a local duplication in the SCN1A gene cluster. The SCN7A gene has undergone rapid evolution and has lost the ability to cause action potentials-instead, it functions as a sodium sensor. The three genes in the SCN5A cluster were translocated from the HOX-bearing chromosome in a mammalian ancestor along with several adjacent genes. This evolutionary scenario is supported by the adjacent TGF-β receptor superfamily (comprised of five distinct families) and the cysteine-serine-rich nuclear protein gene family as well as the HOX clusters. The independent expansions of the SCN repertoires in tetrapods and teleosts suggest that the functional diversification may differ between the two lineages.     

Rapid turnover and species-specificity of vomeronasal pheromone receptor genes in mice and rats

Pheromones are used by individuals of the same species to elicit behavioral or physiological changes, and they are perceived primarily by the vomeronasal organ (VNO) in terrestrial vertebrates. VNO pheromone receptors are encoded by the V1r and V2r gene superfamilies in mammals. A comparison of the V1r and V2r repertoires between closely related species can provide significant insights into the evolutionary genetic mechanisms responsible for species-specific pheromone communications. A total of 137 putatively functional V1r genes of 12 families were previously identified from the mouse genome. We report the identification of 95 putatively functional V1r genes from the draft rat genome sequence. These genes map primarily to four blocks in two chromosomes. The rat V1r genes can be phylogenetically grouped into 10 families, which are shared with mouse, and 2 new families, which are rat-specific. Even in many shared families, gene numbers differ between the two species, apparently due to frequent gene duplication and pseudogenization after the separation of the two species. Molecular dating suggests that most of the rat V1r families emerged before or during the radiation of mammalian orders, but many duplications within families occurred as recently as in the past 10 million years (MY). Our results show that the evolution of the V1r repertoire is characterized by exceptionally fast gene turnover via gains and losses of individual genes, suggesting rapid and substantial changes in pheromone communication between species.     

Rapid changes in the exon/intron structure of a mammalian thrombin inhibitor gene

The genomic organization of the heparin cofactor II (HCII) gene from rat and mouse was investigated and compared with their human counterpart. The genes share a common core structure consisting of five exons interrupted by four introns, but the mouse and rat gene reveal individual additional features. A unique differentially spliced exon is present in the 5'-untranslated region of the rat gene, which most probably has arisen de novo by point mutations in intronic sequences of the ancestor gene. In the mouse HCII gene, a novel intron/exon boundary has been created due to the presence of an additional DNA segment, which simultaneously provides a 3'-splice site and a polypyrimidine stretch leading to an alternatively used exon of increased size. Our data suggest that, in contrast to most other mammalian genes, the exon/intron pattern of the gene coding for HCII is in dynamic evolution.     

The cytochrome P450 genesis locus: the origin and evolution of animal cytochrome P450s

The neighbourhoods of cytochrome P450 (CYP) genes in deuterostome genomes, as well as those of the cnidarians Nematostella vectensis and Acropora digitifera and the placozoan Trichoplax adhaerens were examined to find clues concerning the evolution of CYP genes in animals. CYP genes created by the 2R whole genome duplications in chordates have been identified. Both microsynteny and macrosynteny were used to identify genes that coexisted near CYP genes in the animal ancestor. We show that all 11 CYP clans began in a common gene environment. The evidence implies the existence of a single locus, which we term the ‘cytochrome P450 genesis locus’, where one progenitor CYP gene duplicated to create a tandem set of genes that were precursors of the 11 animal CYP clans: CYP Clans 2, 3, 4, 7, 19, 20, 26, 46, 51, 74 and mitochondrial. These early CYP genes existed side by side before the origin of cnidarians, possibly with a few additional genes interspersed. The Hox gene cluster, WNT genes, an NK gene cluster and at least one ARF gene were close neighbours to this original CYP locus. According to this evolutionary scenario, the CYP74 clan originated from animals and not from land plants nor from a common ancestor of plants and animals. The CYP7 and CYP19 families that are chordate-specific belong to CYP clans that seem to have originated in the CYP genesis locus as well, even though this requires many gene losses to explain their current distribution. The approach to uncovering the CYP genesis locus overcomes confounding effects because of gene conversion, sequence divergence, gene birth and death, and opens the way to understanding the biodiversity of CYP genes, families and subfamilies, which in animals has been obscured by more than 600 Myr of evolution.     

Gene rearrangements in gekkonid mitochondrial genomes with shuffling,loss, and reassignment of tRNA genes

Background

Vertebrate mitochondrial genomes (mitogenomes) are 16–18 kbp double-stranded circular DNAs that encode a set of 37 genes. The arrangement of these genes and the major noncoding region is relatively conserved through evolution although gene rearrangements have been described for diverse lineages. The tandem duplication-random loss model has been invoked to explain the mechanisms of most mitochondrial gene rearrangements. Previously reported mitogenomic sequences for geckos rarely included gene rearrangements, which we explore in the present study.

Results

We determined seven new mitogenomic sequences from Gekkonidae using a high-throughput sequencing method. The Tropiocolotes tripolitanus mitogenome involves a tandem duplication of the gene block: tRNA^Arg, NADH dehydrogenase subunit 4L, and NADH dehydrogenase subunit 4. One of the duplicate copies for each protein-coding gene may be pseudogenized. A duplicate copy of the tRNA^Arg gene appears to have been converted to a tRNA^Gln gene by a C to T base substitution at the second anticodon position, although this gene may not be fully functional in protein synthesis. The Stenodactylus petrii mitogenome includes several tandem duplications of tRNA^Leu genes, as well as a translocation of the tRNA^Ala gene and a putative origin of light-strand replication within a tRNA gene cluster. Finally, the Uroplatus fimbriatus and U. ebenaui mitogenomes feature the apparent loss of the tRNA^Glu gene from its original position. Uroplatus fimbriatus appears to retain a translocated tRNA^Glu gene adjacent to the 5’ end of the major noncoding region.

Conclusions

The present study describes several new mitochondrial gene rearrangements from Gekkonidae. The loss and reassignment of tRNA genes is not very common in vertebrate mitogenomes and our findings raise new questions as to how missing tRNAs are supplied and if the reassigned tRNA gene is fully functional. These new examples of mitochondrial gene rearrangements in geckos should broaden our understanding of the evolution of mitochondrial gene arrangements.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-930) contains supplementary material, which is available to authorized users.     

Origins of flagellar gene operons and secondary flagellar systems

Forty-one flagellated species representing 11 bacterial phyla were used to investigate the origin of secondary flagellar systems and the structure and formation of flagellar gene operons over the course of bacterial evolution. Secondary (i.e., lateral) flagellar systems, which are harbored by five of the proteobacterial species considered, originated twice, once in the alphaproteobacterial lineage and again in the common ancestor of the Beta- and Gammaproteobacteria. The order and organization of flagellar genes have undergone extensive shuffling and rearrangement among lineages, and based on the phylogenetic distributions of flagellar gene complexes, the flagellar gene operons existed as small, usually two-gene units in the ancestor of Bacteria and have expanded through the recruitment of new genes and fusion of gene units. In contrast to the evolutionary trend towards larger flagellar gene complexes, operon structures have been highly disrupted through gene disassociation and rearrangements in the Epsilon- and Alphaproteobacteria. These results demonstrate that the genetic basis of this ancient and structurally conserved organelle has been subject to many lineage-specific modifications.     

A Revised Evolutionary History of the CYP1A Subfamily: Gene Duplication, Gene Conversion, and Positive Selection

Members of cytochrome P450 subfamily 1A (CYP1As) are involved in detoxification and bioactivation of common environmental pollutants. Understanding the functional evolution of these genes is essential to predicting and interpreting species differences in sensitivity to toxicity caused by such chemicals. The CYP1A gene subfamily comprises a single ancestral representative in most fish species and two paralogs in higher vertebrates, including birds and mammals. Phylogenetic analysis of complete coding sequences suggests that mammalian and bird paralog pairs (CYP1A1/2 and CYP1A4/5, respectively) are the result of independent gene duplication events. However, comparison of vertebrate genome sequences revealed that CYP1A genes lie within an extended region of conserved fine-scale synteny, suggesting that avian and mammalian CYP1A paralogs share a common genomic history. Algorithms designed to detect recombination between nucleotide sequences indicate that gene conversion has homogenized most of the length of the chicken CYP1A genes, as well as the 5′ end of mammalian CYP1As. Together, these data indicate that avian and mammalian CYP1A paralog pairs resulted from a single gene duplication event and that extensive gene conversion is responsible for the exceptionally high degree of sequence similarity between CYP1A4 and CYP1A5. Elevated nonsynonymous/synonymous substitution ratios within a putatively unconverted stretch of ∼250 bp suggests that positive selection may have reduced the effective rate of gene conversion in this region, which contains two substrate recognition sites. This work significantly alters our understanding of functional evolution in the CYP1A subfamily, suggesting that gene conversion and positive selection have been the dominant processes of sequence evolution. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Yves Van de Peer]     

Diversity and evolution of 5S rRNA gene family organization in Pythium

The 5S rRNA gene family organization among 87 species and varieties of Pythium was investigated to assess evolutionary stability of the two patterns detected and to determine which pattern is likely the ancestral state in the genus. Species with filamentous sporangia (Groups A-C according to the ITS phylogenetic tree for Pythium) had 5S genes linked to the rDNA repeat that were predominantly coded for on the DNA strand opposite to the one with the other rRNA genes (‘inverted’ orientation). A small group of species with contiguous sporangia (Group D) is related to Groups A-C but had unlinked 5S genes. The main group of species with spherical zoosporangia (Groups E-J) generally had unlinked 5S genes in tandem arrays. The six species in Group K, although they also have spherical sporangia, had linked genes on the same strand as the other rRNA genes ‘non-inverted’ and most of them had pairs of tandem 5S genes. The evolutionary stability of 5S sequence organization was compared with the stability of morphological characters as interpreted from a phylogeny based on ITS sequence analysis. Features of 5S sequence organization were found to be just as consistent within groups as were the morphological characters. To determine the ancestral type of 5S family organization, a survey of Phytophthora strains was conducted to supply an outgroup reference. The most parsimonious interpretation of the data in this survey yielded the tentative conclusion that the linked condition of the 5S sequences was ancestral.