Concerted and birth-and-death evolution of multigene families

Molecular and cytological characterization of 5S rDNA inOryzaspecies: genomic organization and phylogenetic implications

We investigated the molecular characteristics and chromosomal organization of 5S rDNA in the genus Oryza , including diploid and tetraploid species. A phylogenetic tree of Oryza species was constructed based on the non-transcribed spacer sequences of 5S rDNA, and some novel relationships were discovered. Specifically, comparative sequence analysis of 5S rDNA in several wild rice species showed unique characteristics inconsistent with the model of concerted evolution: (1) multiple distinct 5S rDNA types were detected within a species, leading to intraspecific divergence of 5S rDNA; (2) multiple identical 5S rDNA types were shared among species, resulting in interspecies clustering of 5S rDNA types; and (3) intraspecific nucleotide diversity was detected within a 5S rDNA class. Our results obtained by fluorescence in situ hybridization revealed that each rice species studied contained only one 5S rDNA locus with two hybridization sites, which were located on either chromosome 7 or chromosome 11. These results suggest that different 5S rDNA classes within the rice genome were arranged together and that one pair of 5S rDNA loci from a diploid progenitor of the tetraploid species might have been lost during evolution. Taken together, our data show that 5S rDNA in rice species is more informative at the gene level than at the chromosome level.

Evolutionary origins of a novel host plant detoxification gene in butterflies

Chemical interactions between plants and their insect herbivores provide an excellent opportunity to study the evolution of species interactions on a molecular level. Here, we investigate the molecular evolutionary events that gave rise to a novel detoxifying enzyme (nitrile-specifier protein [NSP]) in the butterfly family Pieridae, previously identified as a coevolutionary key innovation. By generating and sequencing expressed sequence tags, genomic libraries, and screening databases we found NSP to be a member of an insect-specific gene family, which we characterized and named the NSP-like gene family. Members consist of variable tandem repeats, are gut expressed, and are found across Insecta evolving in a dynamic, ongoing birth-death process. In the Lepidoptera, multiple copies of single-domain major allergen genes are present and originate via tandem duplications. Multiple domain genes are found solely within the brassicaceous-feeding Pieridae butterflies, one of them being NSP and another called major allergen (MA). Analyses suggest that NSP and its paralog MA have a unique single-domain evolutionary origin, being formed by intragenic domain duplication followed by tandem whole-gene duplication. Duplicates subsequently experienced a period of relaxed constraint followed by an increase in constraint, perhaps after neofunctionalization. NSP and its ortholog MA are still experiencing high rates of change, reflecting a dynamic evolution consistent with the known role of NSP in plant-insect interactions. Our results provide direct evidence to the hypothesis that gene duplication is one of the driving forces for speciation and adaptation, showing that both within- and whole-gene tandem duplications are a powerful force underlying evolutionary adaptation.

Biological applications of the theory of birth-and-death processes

In this review, we discuss applications of the theory of birth-and-death processes to problems in biology, primarily, those of evolutionary genomics. The mathematical principles of the theory of these processes are briefly described. Birth-and-death processes, with some straightforward additions such as innovation, are a simple, natural and formal framework for modeling a vast variety of biological processes such as population dynamics, speciation, genome evolution, including growth of paralogous gene families and horizontal gene transfer and somatic evolution of cancers. We further describe how empirical data, e.g. distributions of paralogous gene family size, can be used to choose the model that best reflects the actual course of evolution among different versions of birth-death-and-innovation models. We conclude that birth-and-death processes, thanks to their mathematical transparency, flexibility and relevance to fundamental biological processes, are going to be an indispensable mathematical tool for the burgeoning field of systems biology.

Inferring natural selection operating on conservative and radical substitution at single amino acid sites

Natural selection operating on amino acid substitution at single amino acid sites can be detected by comparing the rates of synonymous (r(S)) and nonsynonymous (r(N)) nucleotide substitution at single codon sites. Amino acid substitutions can be classified as conservative or radical according to whether they retain the properties of the substituted amino acid. Here methods for comparing the rates of conservative (r(C)) and radical (r(R)) nonsynonymous substitution with r(S) at single codon sites were developed to detect natural selection operating on these substitutions at single amino acid sites. A method for comparing r(C) and r(R) at single codon sites was also developed to detect biases toward these substitutions at single amino acid sites. Charge was used as the property of the amino acids. In a computer simulation, false-positive rates of these methods were always < 5%, unless termination sites were included in the computation of the numbers of sites and estimates of transition/transversion rate ratio were highly biased. The frequency of detection of natural selection operating on conservative substitution was almost independent of the presence of natural selection operating on radical substitution, and vice versa. Natural selection operating specifically on conservative and radical substitution was detected more efficiently by comparing r(S) with r(C) and r(S) with r(R) than by comparing r(S) with r(N). These methods also appeared to be robust against the occurrence of recombination during evolution. In an analysis of class I human leukocyte antigen, negative selection operating on conservative substitution, but not positive selection operating on radical substitution, was observed at some of the codon sites with r(R) > r(C), suggesting that r(R) > r(C) may not necessarily be an indicator of positive selection operating on radical substitution.

T2prhd: a tool to study the patterns of repeat evolution

Background

The models developed to characterize the evolution of multigene families (such as the birth-and-death and the concerted models) have also been applied on the level of sequence repeats inside a gene/protein. Phylogenetic reconstruction is the method of choice to study the evolution of gene families and also sequence repeats in the light of these models. The characterization of the gene family evolution in view of the evolutionary models is done by the evaluation of the clustering of the sequences with the originating loci in mind. As the locus represents positional information, it is straightforward that in the case of the repeats the exact position in the sequence should be used, as the simple numbering according to repeat order can be misleading.

Results

We have developed a novel rapid visual approach to study repeat evolution, that takes into account the exact repeat position in a sequence. The "pairwise repeat homology diagram" visualizes sequence repeats detected by a profile HMM in a pair of sequences and highlights their homology relations inferred by a phylogenetic tree. The method is implemented in a Perl script (t2prhd) available for downloading at and is also accessible as an online tool at . The power of the method is demonstrated on the EGF-like and fibronectin-III-like (Fn-III) domain repeats of three selected mammalian Tenascin sequences.

Conclusion

Although pairwise repeat homology diagrams do not carry all the information provided by the phylogenetic tree, they allow a rapid and intuitive assessment of repeat evolution. We believe, that t2prhd is a helpful tool with which to study the pattern of repeat evolution. This method can be particularly useful in cases of large datasets (such as large gene families), as the command line interface makes it possible to automate the generation of pairwise repeat homology diagrams with the aid of scripts.

Evolutionary dynamics of the immunoglobulin heavy chain variable region genes in vertebrates

Immunoglobulin heavy chains are polypeptides encoded by four genes: variable (IGHV), joining (IGHJ), diversity (IGHD), and constant (IGHC) region genes. The number of IGHV genes varies from species to species. To understand the evolution of the IGHV multigene family, we identified and analyzed the IGHV sequences from 16 vertebrate species. The results show that the numbers of functional and nonfunctional IGHV genes among different species are positively correlated. The number of IGHV genes is relatively stable in teleosts, but the intragenomic sequence variation is generally higher in teleosts than in tetrapods. The IGHV genes in tetrapods can be classified into three phylogenetic clans (I, II, and III). The clan III and/or II genes are relatively abundant, whereas clan I genes exist in small numbers or are absent in most species. The genomic organization of clan I, II, and III IGHV genes varies considerably among species, but the entire IGHV locus seems to be conserved in the subtelomeric or near-centromeric region of chromosome. The presence or absence of specific IGHV clan members and the lineage-specific expansion and contraction of IGHV genes indicate that the IGHV locus continues to evolve in a species-specific manner. Our results suggest that the evolution of IGHV multigene family is more complex than previously thought and that several factors may act synergistically for the development of antibody repertoire.

Purifying selection and birth-and-death evolution in the class II hydrophobin gene families of the ascomycete Trichoderma/Hypocrea

BACKGROUND

Hydrophobins are proteins containing eight conserved cysteine residues that occur uniquely in mycelial fungi. Their main function is to confer hydrophobicity to fungal surfaces in contact with air or during attachment of hyphae to hydrophobic surfaces of hosts, symbiotic partners or themselves resulting in morphogenetic signals. Based on their hydropathy patterns and solubility characteristics, hydrophobins are divided into two classes (I and II), the latter being found only in ascomycetes.

RESULTS

We have investigated the mechanisms driving the evolution of the class II hydrophobins in nine species of the mycoparasitic ascomycetous genus Trichoderma/Hypocrea, using three draft sequenced genomes (H. jecorina = T. reesei, H. atroviridis = T. atroviride; H. virens = T. virens) an additional 14,000 ESTs from six other Trichoderma spp. (T. asperellum, H. lixii = T. harzianum, T. aggressivum var. europeae, T. longibrachiatum, T. cf. viride). The former three contained six, ten and nine members, respectively. Ten is the highest number found in any ascomycete so far. All the hydrophobins we examined had the conserved four beta-strands/one helix structure, which is stabilized by four disulfide bonds. In addition, a small number of these hydrophobins (HFBs)contained an extended N-terminus rich in either proline and aspartate, or glycine-asparagine. Phylogenetic analysis reveals a mosaic of terminal clades containing duplicated genes and shows only three reasonably supported clades. Calculation of the ratio of differences in synonymous vs. non-synonymous nucleotide substitutions provides evidence for strong purifying selection (KS/Ka >> 1). A genome database search for class II HFBs from other ascomycetes retrieved a much smaller number of hydrophobins (2-4) from each species, and most were from Sordariomycetes. A combined phylogeny of these sequences with those of Trichoderma showed that the Trichoderma HFBs mostly formed their own clades, whereas those of other Sordariomycetes occurred in shared clades.

CONCLUSION

Our study shows that the genus Trichoderma/Hypocrea has a proliferated arsenal of class II hydrophobins which arose by birth-and-death evolution followed by purifying selection.

Rapid rates of lineage-specific gene duplication and deletion in the alpha-globin gene family

Phylogeny reconstructions of the globin gene families have revealed that paralogous genes within species are often more similar to one another than they are to their orthologous counterparts in closely related species. This pattern has been previously attributed to mechanisms of concerted evolution such as interparalog gene conversion that homogenize sequence variation between tandemly duplicated genes and therefore create the appearance of recent common ancestry. Here we report a comparative genomic analysis of the alpha-globin gene family in mammals that reveal a surprisingly high rate of lineage-specific gene duplication and deletion via unequal crossing-over. Results of our analysis reveal that patterns of sequence similarity between paralogous alpha-like globin genes from the same species are only partly explained by concerted evolution between preexisting gene duplicates. In a number of cases, sequence similarity between paralogous sequences from the same species is attributable to recent ancestry between the products of de novo gene duplications. As a result of this surprisingly rapid rate of gene gain and loss, many mammals possess alpha-like globin genes that have no orthologous counterparts in closely related species. The resultant variation in gene copy number among species may represent an important source of regulatory variation that affects physiologically important aspects of blood oxygen transport and aerobic energy metabolism.

Molecular diversity at the plant-pathogen interface

Plants have evolved a robust innate immune system that exhibits striking similarities as well as significant differences with various metazoan innate immune systems. For example, plants are capable of perceiving pathogen-associated molecular patterns through pattern recognition receptors that bear structural similarities to animal Toll-like receptors. In addition, plants have evolved a second surveillance system based on cytoplasmic "NB-LRR" proteins (nucleotide-binding, leucine-rich repeat) that are structurally similar to animal nucleotide-binding and oligomerization domain (NOD)-like receptors. Plant NB-LRR proteins do not detect PAMPs; rather, they perceive effector proteins that pathogens secrete into plant cells to promote virulence. This review summarizes the current state of knowledge about the molecular functionality and evolution of these immune surveillance genes.

Evolutionary mechanisms underlying secondary metabolite diversity

The enormous chemical diversity and the broad range of biological activities of secondary metabolites raise many questions about their role in nature and the specific traits leading to their evolution. The answers to these questions will not only be of fundamental interest but may also provide lessons that could help to improve the screening protocols of pharmaceutical companies and strategies for rational secondary metabolite engineering. In this review, we try to dissect evolutionary principles leading to the emergence, distribution, diversification and selection of genes involved in secondary metabolite biosyntheses. We give an overview about recent insights into the evolution of the different types of polyketide synthases (PKS) in microorganisms and plants and highlight unique mechanisms underlying polyketide diversity. Although phylogenetic and experimental data have significantly increased our knowledge about the role and evolution of secondary metabolites in the last decades there is still much dissent about the impact of natural selection. In order to understand the evolution towards metabolic diversity we therefore need more thorough investigations of the ecological role of secondary metabolites in the future.

Significant gene content variation characterizes the genomes of inbred mouse strains

The contribution to genetic diversity of genomic segmental copy number variations (CNVs) is less well understood than that of single-nucleotide polymorphisms (SNPs). While less frequent than SNPs, CNVs have greater potential to affect phenotype. In this study, we have performed the most comprehensive survey to date of CNVs in mice, analyzing the genomes of 42 Mouse Phenome Consortium priority strains. This microarray comparative genomic hybridization (CGH)-based analysis has identified 2094 putative CNVs, with an average of 10 Mb of DNA in 51 CNVs when individual mouse strains were compared to the reference strain C57BL/6J. This amount of variation results in gene content that can differ by hundreds of genes between strains. These genes include members of large families such as the major histocompatibility and pheromone receptor genes, but there are also many singleton genes including genes with expected phenotypic consequences from their deletion or amplification. Using a whole-genome association analysis, we demonstrate that complex multigenic phenotypes, such as food intake, can be associated with specific copy number changes.

Euglena light-harvesting complexes are encoded by multifarious polyprotein mRNAs that evolve in concert

Light-harvesting complexes (LHCs) are a superfamily of chlorophyll- and carotenoid-binding proteins that are responsible for the capture of light energy and its transfer to the photosynthetic reaction centers. Unlike those of most eukaryotes, the LHCs of Euglena gracilis are translated from large mRNAs, producing polyprotein precursors consisting of multiple concatenated LHC subunits that are separated by conserved decapeptide linkers. These precursors are posttranslationally targeted to the chloroplast and cleaved into individual proteins. We analyzed expressed sequence tags from Euglena to further characterize the structural features of the LHC polyprotein-coding genes and to examine the evolution of this multigene family. Of the 19 different LHC transcriptional units we detected, 17 encoded polyproteins composed of both tandem and nontandem repeats of LHC subunits; organizations that likely occurred through unequal crossing-over. Of the 2 nonpolyprotein-encoding LHC transcripts detected, 1 evolved from the truncation of a polyprotein-coding gene. Duplication of LHC polyprotein-coding genes was particularly important in the LHCI gene family where multiple paralogous sequences were detected. Intriguingly, several of the individual LHC-coding subunits both within and between transcriptional units appeared to be evolving in concert, suggesting that gene conversion has been a significant mechanism for LHC evolution in Euglena.

Species-specific evolution of immune receptor tyrosine based activation motif-containing CEACAM1-related immune receptors in the dog

BACKGROUND

Although the impact of pathogens on the evolution of the mammalian immune system is still under debate, proteins, which both regulate immune responses and serve as cellular receptors for pathogens should be at the forefront of pathogen-driven host evolution. The CEA (carcinoembryonic antigen) gene family codes for such proteins and indeed shows tremendous species-specific variation between human and rodents. Since little is known about the CEA gene family in other lineages of placental mammals, we expected to gain new insights into the evolution of the rapidly diverging CEA family by analyzing the CEA family of the dog.

RESULTS

Here we describe the complete CEA gene family in the dog. We found that the gene coding for the ITIM-bearing immunoregulatory molecule CEACAM1 gave rise to a recent expansion of the canine CEA gene family by gene duplication, similar to that previously found in humans and mice. However, while the murine and human CEACAMs (carcinoembryonic antigen-related cell adhesion molecules) are predominantly secreted and GPI-anchored, respectively, in the dog, most of the CEACAMs represent ITAM-bearing transmembrane proteins. One of these proteins, CEACAM28, exhibits nearly complete sequence identity with the ligand-binding N domain of CEACAM1, but antagonizing signaling motifs in the cytoplasmic tail. Comparison of nonsynonymous and synonymous substitutions indicates that the CEACAM28 N domain is under the strongest purifying selection of all canine CEACAM1-related CEACAMs. In addition, CEACAM28 shows a similar expression pattern in resting immune cells and tissues as CEACAM1. However, upon activation CEACAM28 mRNA and CEACAM1 mRNA are differentially regulated.

CONCLUSION

Thus, CEACAM1 and CEACAM28 are the first paired immune receptors identified within the CEA gene family, which are expressed on T cells and are most likely involved in the fine-tuning of T cell responses. The direction of gene conversion accompanied by purifying selection and expression in immune cells suggests the possibility that CEACAM28 evolved in response to selective pressure imposed by species-specific pathogens.

Birth and death of gene overlaps in vertebrates

Background

Between five and fourteen per cent of genes in the vertebrate genomes do overlap sharing some intronic and/or exonic sequence. It was observed that majority of these overlaps are not conserved among vertebrate lineages. Although several mechanisms have been proposed to explain gene overlap origination the evolutionary basis of these phenomenon are still not well understood. Here, we present results of the comparative analysis of several vertebrate genomes. The purpose of this study was to examine overlapping genes in the context of their evolution and mechanisms leading to their origin.

Results

Based on the presence and arrangement of human overlapping genes orthologs in rodent and fish genomes we developed 15 theoretical scenarios of overlapping genes evolution. Analysis of these theoretical scenarios and close examination of genomic sequences revealed new mechanisms leading to the overlaps evolution and confirmed that many of the vertebrate gene overlaps are not conserved. This study also demonstrates that repetitive elements contribute to the overlapping genes origination and, for the first time, that evolutionary events could lead to the loss of an ancient overlap.

Conclusion

Birth as well as most probably death of gene overlaps occurred over the entire time of vertebrate evolution and there wasn't any rapid origin or 'big bang' in the course of overlapping genes evolution. The major forces in the gene overlaps origination are transposition and exaptation. Our results also imply that origin of overlapping genes is not an issue of saving space and contracting genomes size.

X-tox: an atypical defensin derived family of immune-related proteins specific to Lepidoptera

We report here the isolation in Spodoptera frugiperda (Lepidoptera) of an immune-related protein (hereafter named Spod-11-tox), characterized by imperfectly conserved tandem repeats of 11 cysteine-stabilized alpha beta motifs (CS-alphabeta), the structural scaffold characteristic of invertebrate defensins and scorpion toxins. Spod-11-tox orthologs were only found in Lepidopteran species, suggesting that this new protein family (named X-tox) is specific to this insect order. Moreover, phylogenetic analysis suggests that X-tox proteins represent a new class of proteins restricted to Lepidoptera and likely derived from Lepidopteran defensins. In S. frugiperda, analysis of gene expression revealed that spod-11-tox is rapidly induced by infection. However, and conversely to what is known for most insect antimicrobial peptides (AMP), spod-11-tox is mainly expressed in blood cells. Moreover, recombinant Spod-11-tox produced in the Sf9 cell line does not show any antimicrobial activity. Altogether, these results suggest that although X-tox proteins are derived from defensins, they may play a different and still unknown role in Lepidoptera immune response.

Co-regulation of a large and rapidly evolving repertoire of odorant receptor genes

The olfactory system meets niche- and species-specific demands by an accelerated evolution of its odorant receptor repertoires. In this review, we describe evolutionary processes that have shaped olfactory and vomeronasal receptor gene families in vertebrate genomes. We emphasize three important periods in the evolution of the olfactory system evident by comparative genomics: the adaptation to land in amphibian ancestors, the decline of olfaction in primates, and the delineation of putative pheromone receptors concurrent with rodent speciation. The rapid evolution of odorant receptor genes, the sheer size of the repertoire, as well as their wide distribution in the genome, presents a developmental challenge: how are these ever-changing odorant receptor repertoires coordinated within the olfactory system? A central organizing principle in olfaction is the specialization of sensory neurons resulting from each sensory neuron expressing only ~one odorant receptor allele. In this review, we also discuss this mutually exclusive expression of odorant receptor genes. We have considered several models to account for co-regulation of odorant receptor repertoires, as well as discussed a new hypothesis that invokes important epigenetic properties of the system.

The fungus-specific HET domain mediates programmed cell death in Podospora anserina

Vegetative incompatibility is a programmed cell death reaction that occurs when fungal cells of unlike genotypes fuse. Genes defining vegetative incompatibility (het genes) are highly polymorphic, and most if not all incompatibility systems include a protein partner bearing the fungus-specific domain termed the HET domain. The nonallelic het-C/het-E incompatibility system is the best-characterized incompatibility system in Podospora anserina. Cell death is triggered by interaction of specific alleles of het-C, encoding a glycolipid transfer protein, and het-E, encoding a HET domain and a WD repeat domain involved in recognition. We show here that overexpression of the isolated HET domain from het-E results in cell death. This cell death is characterized by induction of autophagy, increased vacuolization, septation, and production of lipid droplets, which are hallmarks of cell death by incompatibility. In addition, the HET domain lethality is suppressed by the same mutations as vegetative incompatibility, but not by the inactivation of het-C. These results establish the HET domain as the mediator of cell death by incompatibility and lead to a modular conception of incompatibility systems whereby recognition is ensured by the variable regions of incompatibility proteins and cell death is triggered by the HET domain.

Unique nucleotide polymorphism of ankyrin gene cluster in Arabidopsis

The ankyrin (ANK) gene cluster is a part of a multigene family encoding ANK transmembrane proteins in Arabidopsis thaliana, and plays an important role in protein-protein interactions and in signal pathways. In contrast to other regions of a genome, the ANK gene cluster exhibits an extremely high level of DNA polymorphism in an approximately 5-kb region, without apparent decay. Phylogenetic analysis detects two clear, deeply differentiated haplotypes (dimorphism). The divergence between haplotypes of accession Col-0 and Ler-0 (Hap-C and Hap-L) is estimated to be 10.7%, approximately equal to the 10.5% average divergence between A. thaliana and A. lyrata. Sequence comparisons for the ANK gene cluster homologues in Col-0 indicate that the members evolve independently, and that the similarity among paralogues is lower than between alleles. Very little intralocus recombination or gene conversion is detected in ANK regions. All these characteristics of the ANK gene cluster are consistent with a tandem gene duplication and birth-and-death process. The possible mechanisms for and implications of this elevated nucleotide variation are also discussed, including the suggestion of balancing selection.

Gene Duplication: A Drive for Phenotypic Diversity and Cause of Human Disease

Gene duplication is one of the key factors driving genetic innovation, i.e., producing novel genetic variants. Although the contribution of whole-genome and segmental duplications to phenotypic diversity across species is widely appreciated, the phenotypic spectrum and potential pathogenicity of small-scale duplications in individual genomes are less well explored. This review discusses the nature of small-scale duplications and the phenotypes produced by such duplications. Phenotypic variation and disease phenotypes induced by duplications are more diverse and widespread than previously anticipated, and duplications are a major class of disease-related genomic variation. Pathogenic duplications particularly involve dosage-sensitive genes with both similar and dissimilar over- and underexpression phenotypes, and genes encoding proteins with a propensity to aggregate. Phenotypes related to human-specific copy number variation in genes regulating environmental responses and immunity are increasingly recognized. Small genomic duplications containing defense-related genes also contribute to complex common phenotypes.

A model of segmental duplication formation in Drosophila melanogaster

Segmental duplications (SDs) are low-copy repeats of DNA segments that have long been recognized to be involved in genome organization and evolution. But, to date, the mechanism of their formation remains obscure. We propose a model for SD formation that we name "duplication-dependent strand annealing" (DDSA). This model is a variant of the synthesis-dependent strand annealing (SDSA) model--a double-strand break (DSB) homologous repair model. DSB repair in Drosophila melanogaster genome usually occurs primarily through homologous repair, more preferentially through the SDSA model. The DDSA model predicts that after a DSB, the search for an ectopic homologous region--here a repeat--initiates the repair. As expected by the model, the analysis of SDs detected by a computational analysis of the D. melanogaster genome indicates a high enrichment in transposable elements at SD ends. It shows moreover a preferential location of SDs in heterochromatic regions. The model has the advantage of also predicting specific traces left during synthesis. The observed traces support the DDSA model as one model of formation of SDs in D. melanogaster genome. The analysis of these DDSA signatures suggests moreover a sequestration of the dissociated strand in the repair complex.

Evolution of mammalian chitinase(-like) members of family 18 glycosyl hydrolases

Family 18 of glycosyl hydrolases encompasses chitinases and so-called chi-lectins lacking enzymatic activity due to amino acid substitutions in their active site. Both types of proteins widely occur in mammals although these organisms lack endogenous chitin. Their physiological function(s) as well as evolutionary relationships are still largely enigmatic. An overview of all family members is presented and their relationships are described. Molecular phylogenetic analyses suggest that both active chitinases (chitotriosidase and AMCase) result from an early gene duplication event. Further duplication events, followed by mutations leading to loss of chitinase activity, allowed evolution of the chi-lectins. The homologous genes encoding chitinase(-like) proteins are clustered in two distinct loci that display a high degree of synteny among mammals. Despite the shared chromosomal location and high homology, individual genes have evolved independently. Orthologs are more closely related than paralogues, and calculated substitution rate ratios indicate that protein-coding sequences underwent purifying selection. Substantial gene specialization has occurred in time, allowing for tissue-specific expression of pH optimized chitinases and chi-lectins. Finally, several family 18 chitinase-like proteins are present only in certain lineages of mammals, exemplifying recent evolutionary events in the chitinase protein family.

Extensive gains and losses of olfactory receptor genes in mammalian evolution

Odor perception in mammals is mediated by a large multigene family of olfactory receptor (OR) genes. The number of OR genes varies extensively among different species of mammals, and most species have a substantial number of pseudogenes. To gain some insight into the evolutionary dynamics of mammalian OR genes, we identified the entire set of OR genes in platypuses, opossums, cows, dogs, rats, and macaques and studied the evolutionary change of the genes together with those of humans and mice. We found that platypuses and primates have <400 functional OR genes while the other species have 800–1,200 functional OR genes. We then estimated the numbers of gains and losses of OR genes for each branch of the phylogenetic tree of mammals. This analysis showed that (i) gene expansion occurred in the placental lineage each time after it diverged from monotremes and from marsupials and (ii) hundreds of gains and losses of OR genes have occurred in an order-specific manner, making the gene repertoires highly variable among different orders. It appears that the number of OR genes is determined primarily by the functional requirement for each species, but once the number reaches the required level, it fluctuates by random duplication and deletion of genes. This fluctuation seems to have been aided by the stochastic nature of OR gene expression.

Characterization and evolution of the cell cycle-associated mob domain-containing proteins in eukaryotes

The MOB family includes a group of cell cycle-associated proteins highly conserved throughout eukaryotes, whose founding members are implicated in mitotic exit and co-ordination of cell cycle progression with cell polarity and morphogenesis. Here we report the characterization and evolution of the MOB domain-containing proteins as inferred from the 43 eukaryotic genomes so far sequenced. We show that genes for Mob-like proteins are present in at least 41 of these genomes, confirming the universal distribution of this protein family and suggesting its prominent biological function. The phylogenetic analysis reveals five distinct MOB domain classes, showing a progressive expansion of this family from unicellular to multicellular organisms, reaching the highest number in mammals. Plant Mob genes appear to have evolved from a single ancestor, most likely after the loss of one or more genes during the early stage of Viridiplantae evolutionary history. Three of the Mob classes are widespread among most of the analyzed organisms. The possible biological and molecular function of Mob proteins and their role in conserved signaling pathways related to cell proliferation, cell death and cell polarity are also presented and critically discussed.

Large-scale, lineage-specific expansion of a bric-a-brac/tramtrack/broad complex ubiquitin-ligase gene family in rice

Selective ubiquitination of proteins is directed by diverse families of ubiquitin-protein ligases (or E3s) in plants. One important type uses Cullin-3 as a scaffold to assemble multisubunit E3 complexes containing one of a multitude of bric-a-brac/tramtrack/broad complex (BTB) proteins that function as substrate recognition factors. We previously described the 80-member BTB gene superfamily in Arabidopsis thaliana. Here, we describe the complete BTB superfamily in rice (Oryza sativa spp japonica cv Nipponbare) that contains 149 BTB domain-encoding genes and 43 putative pseudogenes. Amino acid sequence comparisons of the rice and Arabidopsis superfamilies revealed a near equal repertoire of putative substrate recognition module types. However, phylogenetic comparisons detected numerous gene duplication and/or loss events since the rice and Arabidopsis BTB lineages split, suggesting possible functional specialization within individual BTB families. In particular, a major expansion and diversification of a subset of BTB proteins containing Meprin and TRAF homology (MATH) substrate recognition sites was evident in rice and other monocots that likely occurred following the monocot/dicot split. The MATH domain of a subset appears to have evolved significantly faster than those in a smaller core subset that predates flowering plants, suggesting that the substrate recognition module in many monocot MATH-BTB E3s are diversifying to ubiquitinate a set of substrates that are themselves rapidly changing. Intriguing possibilities include pathogen proteins attempting to avoid inactivation by the monocot host.

Structural evolution of nrDNA ITS in Pinaceae and its phylogenetic implications

Nuclear ribosomal DNA (nrDNA) has been considered as an important tool for inferring phylogenetic relationships at many taxonomic levels. In comparison with its fast concerted evolution in angiosperms, nrDNA is symbolized by slow concerted evolution and substantial ITS region length variation in gymnosperms, particularly in Pinaceae. Here we studied structure characteristics, including subrepeat composition, size, GC content and secondary structure, of nrDNA ITS regions of all Pinaceae genera. The results showed that the ITS regions of all taxa studied contained subrepeat units, ranging from 2 to 9 in number, and these units could be divided into two types, longer subrepeat (LSR) without the motif (5'-GGCCACCCTAGTC) and shorter subrepeat (SSR) with the motif. Phylogenetic analyses indicate that the homology of some SSRs still can be recognized, providing important informations for the evolutionary history of nrDNA ITS and phylogeny of Pinaceae. In particular, the adjacent tandem SSRs are not more closely related to one another than they are to remote SSRs in some genera, which may imply that multiple structure variations such as recombination have occurred in the ITS1 region of these groups. This study also found that GC content in the ITS1 region is relevant to its sequence length and subrepeat number, and could provide some phylogenetic information, especially supporting the close relationships among Picea, Pinus, and Cathaya. Moreover, several characteristics of the secondary structure of Pinaceae ITS1 were found as follows: (1) the structure is dominated by several extended hairpins; (2) the configuration complexity is positively correlated with subrepeat number; (3) paired subrepeats often partially overlap at the conserved motif (5'-GGCCACCCTAGTC), and form a long stem, while other subrepeats fold onto itself, leaving part of the conserved motif exposed in hairpin loops.

Large-scale, lineage-specific expansion of a bric-a-brac/tramtrack/broad complex ubiquitin-ligase gene family in rice

Selective ubiquitination of proteins is directed by diverse families of ubiquitin-protein ligases (or E3s) in plants. One important type uses Cullin-3 as a scaffold to assemble multisubunit E3 complexes containing one of a multitude of bric-a-brac/tramtrack/broad complex (BTB) proteins that function as substrate recognition factors. We previously described the 80-member BTB gene superfamily in Arabidopsis thaliana. Here, we describe the complete BTB superfamily in rice (Oryza sativa spp japonica cv Nipponbare) that contains 149 BTB domain-encoding genes and 43 putative pseudogenes. Amino acid sequence comparisons of the rice and Arabidopsis superfamilies revealed a near equal repertoire of putative substrate recognition module types. However, phylogenetic comparisons detected numerous gene duplication and/or loss events since the rice and Arabidopsis BTB lineages split, suggesting possible functional specialization within individual BTB families. In particular, a major expansion and diversification of a subset of BTB proteins containing Meprin and TRAF homology (MATH) substrate recognition sites was evident in rice and other monocots that likely occurred following the monocot/dicot split. The MATH domain of a subset appears to have evolved significantly faster than those in a smaller core subset that predates flowering plants, suggesting that the substrate recognition module in many monocot MATH-BTB E3s are diversifying to ubiquitinate a set of substrates that are themselves rapidly changing. Intriguing possibilities include pathogen proteins attempting to avoid inactivation by the monocot host.

Recent duplication and positive selection of the GAGE gene family

We report that the GAGE gene family of human Cancer/testis antigen (CTA) genes is likely to be in an early stage of its evolution. Members of this gene family are tandemly arranged on the X chromosome only in human, chimpanzee and macaque genomes and share a very high similarity. Phylogenetic trees show that the GAGE gene family began to duplicate after the split of human and chimpanzee. The estimated ages of the duplication events range from 4 million years ago to the present. The Ka/Ks values between the duplicates are significantly greater than 1, indicating that the mutation rate is higher in coding regions than non-coding regions of the genes, which suggests that the GAGE gene family is under positive selection. These findings indicate that the GAGE gene family may be a newly formed gene family undergoing rapid functional evolution.

Estimating meiotic gene conversion rates from population genetic data

Gene conversion plays an important part in shaping genetic diversity in populations, yet estimating the rate at which it occurs is difficult because of the short lengths of DNA involved. We have developed a new statistical approach to estimating gene conversion rates from genetic variation, by extending an existing model for haplotype data in the presence of crossover events. We show, by simulation, that when the rate of gene conversion events is at least comparable to the rate of crossover events, the method provides a powerful approach to the detection of gene conversion and estimation of its rate. Application of the method to data from the telomeric X chromosome of Drosophila melanogaster, in which crossover activity is suppressed, indicates that gene conversion occurs approximately 400 times more often than crossover events. We also extend the method to estimating variable crossover and gene conversion rates and estimate the rate of gene conversion to be approximately 1.5 times higher than the crossover rate in a region of human chromosome 1 with known recombination hotspots.

Screen for ISG15-crossreactive deubiquitinases

Background

The family of ubiquitin-like molecules (UbLs) comprises several members, each of which has sequence, structural, or functional similarity to ubiquitin. ISG15 is a homolog of ubiquitin in vertebrates and is strongly upregulated following induction by type I interferon. ISG15 can be covalently attached to proteins, analogous to ubiquitination and with actual support of ubiquitin conjugating factors. Specific proteases are able to reverse modification with ubiquitin or UbLs by hydrolyzing the covalent bond between their C-termini and substrate proteins. The tail regions of ubiquitin and ISG15 are identical and we therefore hypothesized that promiscuous deubiquitinating proteases (DUBs) might exist, capable of recognizing both ubiquitin and ISG15.

Results

We have cloned and expressed 22 human DUBs, representing the major clades of the USP protease family. Utilizing suicide inhibitors based on ubiquitin and ISG15, we have identified USP2, USP5 (IsoT1), USP13 (IsoT3), and USP14 as ISG15-reactive proteases, in addition to the bona fide ISG15-specific protease USP18 (UBP43). USP14 is a proteasome-associated DUB, and its ISG15 isopeptidase activity increases when complexed with the proteasome.

Conclusions

By evolutionary standards, ISG15 is a newcomer among the UbLs and it apparently not only utilizes the conjugating but also the deconjugating machinery of its more established relative ubiquitin. Functional overlap between these two posttranslational modifiers might therefore be more extensive than previously appreciated and explain the rather innocuous phenotype of ISG15 null mice.

The new mutation theory of phenotypic evolution

Recent studies of developmental biology have shown that the genes controlling phenotypic characters expressed in the early stage of development are highly conserved and that recent evolutionary changes have occurred primarily in the characters expressed in later stages of development. Even the genes controlling the latter characters are generally conserved, but there is a large component of neutral or nearly neutral genetic variation within and between closely related species. Phenotypic evolution occurs primarily by mutation of genes that interact with one another in the developmental process. The enormous amount of phenotypic diversity among different phyla or classes of organisms is a product of accumulation of novel mutations and their conservation that have facilitated adaptation to different environments. Novel mutations may be incorporated into the genome by natural selection (elimination of preexisting genotypes) or by random processes such as genetic and genomic drift. However, once the mutations are incorporated into the genome, they may generate developmental constraints that will affect the future direction of phenotypic evolution. It appears that the driving force of phenotypic evolution is mutation, and natural selection is of secondary importance.

Gene Duplication and the Evolution of Vertebrate Skeletal Mineralization

The mineralized skeleton is a critical innovation that evolved early in vertebrate history. The tissues found in dermal skeletons of ancient vertebrates are similar to the dental tissues of modern vertebrates; both consist of a highly mineralized surface hard tissue, enamel or enameloid, more resilient body dentin, and basal bone. Many proteins regulating mineralization of these tissues are evolutionarily related and form the secretory calcium-binding phosphoprotein (SCPP) family. We hypothesize here the duplication histories of SCPP genes and their common ancestors, SPARC and SPARCL1. At around the same time that Paleozoic jawless vertebrates first evolved mineralized skeleton, SPARCL1 arose from SPARC by whole genome duplication. Then both before and after the split of ray-finned fish and lobe-finned fish, tandem gene duplication created two types of SCPP genes, each residing on the opposite side of SPARCL1. One type was subsequently used in surface tissue and the other in body tissue. In tetrapods, these two types of SCPP genes were separated by intrachromosomal rearrangement. While new SCPP genes arose by duplication, some old genes were eliminated from the genome. As a consequence, phenogenetic drift occurred: while mineralized skeleton is maintained by natural selection, the underlying genetic basis has changed.

Sequence and structure evolved separately in a ribosomal ubiquitin variant

Encoded by a multigene family, ubiquitin is expressed in the form of three precursor proteins, two of which are fusions to the ribosomal subunits S27a and L40. Ubiquitin assists in ribosome biogenesis and also functions as a post-translational modifier after its release from S27a or L40. However, several species do not conserve the ribosomal ubiquitin domains. We report here the solution structure of a distant variant of ubiquitin, found at the N-terminus of S27a in Giardia lamblia, referred to as GlUb(S27a). Despite the considerable evolutionary distance that separates ubiquitin from GlUb(S27a), the structure of GlUb(S27a) is largely identical to that of ubiquitin. The variant domain remains attached to S27a and is part of the assembled holoribosome. Thus, conservation of tertiary structure suggests a role of this variant as a chaperone, while conservation of the primary structure--necessary for ubiquitin's function as a post-translational modifier--is no longer required. Based on these observations, we propose a model to explain the origin of the widespread ubiquitin superfold in eukaryotes.

Intragenomic variation and evolution of the internal transcribed spacer of the rRNA operon in bacteria

Variation in the internal transcribed spacer (ITS) of the rRNA (rrn) operon is increasingly used to infer population-level diversity in bacterial communities. However, intragenomic ITS variation may skew diversity estimates that do not correct for multiple rrn operons within a genome. This study characterizes variation in ITS length, tRNA composition, and intragenomic nucleotide divergence across 155 Bacteria genomes. On average, these genomes encode 4.8 rrn operons (range: 2-15) and contain 2.4 unique ITS length variants (range: 1-12) and 2.8 unique sequence variants (range: 1-12). ITS variation stems primarily from differences in tRNA gene composition, with ITS regions containing tRNA-Ala + tRNA-Ile (48% of sequences), tRNA-Ala or tRNA-Ile (10%), tRNA-Glu (11%), other tRNAs (3%), or no tRNA genes (27%). Intragenomic divergence among paralogous ITS sequences grouped by tRNA composition ranges from 0% to 12.11% (mean: 0.94%). Low divergence values indicate extensive homogenization among ITS copies. In 78% of alignments, divergence is <1%, with 54% showing zero variation and 81% containing at least two identical sequences. ITS homogenization occurs over relatively long sequence tracts, frequently spanning the entire ITS, and is largely independent of the distance (basepairs) between operons. This study underscores the potential contribution of interoperon ITS variation to bacterial microdiversity studies, as well as unequivocally demonstrates the pervasiveness of concerted evolution in the rrn gene family.

Non-neutral evolution of the major histocompatibility complex class II gene DRB1 in the sac-winged bat Saccopteryx bilineata

The immune genes of the major histocompatibility complex (MHC) are classical examples for high levels of genetic diversity and non-neutral evolution. This is particularly true for the regions containing the antigen-binding sites as, for instance, in the exon 2 of the MHC class II gene DRB. We surveyed, for the first time in the order Chiroptera, the genetic diversity within this exon in the sac-winged bat Saccopteryx bilineata. We detected 11 alleles among 85 bats, of which 79 were sampled in one population. Pairwise comparisons revealed that interallelic sequence differences ranged between 3 and 22%, although nucleotide substitutions were not evenly distributed along the exon sequence. This was most probably the result of intragenic recombination. High levels of sequence divergence and significantly more nonsynonymous than synonymous substitutions (d_N/d_S>1) suggest long-term balancing selection. Thus, the data are consistent with the hypothesis that recombination gives rise to new alleles at the DRB locus of the sac-winged bat, and these are maintained in the population through balancing selection. In this respect, the sac-winged bat closely resembles other mammalian species.

Genomic organization, nucleotide sequence analysis of the core histone genes cluster in Chlamys farreri and molecular evolution assessment of the H2A and H2B

This work represents the nucleotide sequence of the core histone gene cluster from scallop Chlamys farreri. The tandemly repeated unit of 5671 bp containing a copy of the four core histone genes H4, H2B, H2A and H3 was amplified and identified by the techniques of homology cloning and genomic DNA walking. All the histone genes in the cluster had the structures in their 3' flanking region which related to the evolution of histone gene expression patterns throughout the cell cycle, including two different termination signals, the hairpin structure and at least one AATAAA polyadenylation signal. In their 5' region, the transcription initiation sites with a conserved sequence of 5'-PyATTCPu-3' known as the CAP site were present in all genes except to H2B, generally 37-45 bp upstream of the start code. Canonical TATA and CAAT boxes were identified only in certain histone genes. In the case of the promoters of H2B and H2A genes, there was a 5'-GATCC-3' element, which had been found to be essential to start transcription at the appropriate site. After this element, in the promoter of H2B, there was another sequence, 5'-GGATCGAAACGTTC-3', which was similar to the consensus sequence of 5'-GGAATAAACGTATTC-3' corresponding to the H2B-specific promoter element. The presence of enhancer sequences (5'-TGATATATG-3') was identified from the H4 and H3 genes, matching perfectly with the consensus sequence defined for histone genes. There were several slightly more complex repetitive DNA in the intergene regions. The presence of the series of conserved sequences and reiterated sequences was consistent with the view that mollusc histone gene cluster arose by duplicating of an ancestral precursor histone gene, the birth-and-death evolution model with strong purifying selection enabled the histone cluster less variation and more conserved function. Meanwhile, the H2A and the H2B were demonstrated to be potential good marks for phylogenetic analysis. All the results will be contributed to the characterization of repeating histone gene families in molluscs.

Finely orchestrated movements: evolution of the ribosomal RNA genes

Evolution of the tandemly repeated ribosomal RNA (rRNA) genes is intriguing because in each species all units within the array are highly uniform in sequence but that sequence differs between species. In this review we summarize the origins of the current models to explain this process of concerted evolution, emphasizing early studies of recombination in yeast and more recent studies in Drosophila and mammalian systems. These studies suggest that unequal crossover is the major driving force in the evolution of the rRNA genes with sister chromatid exchange occurring more often than exchange between homologs. Gene conversion is also believed to play a role; however, direct evidence for its involvement has not been obtained. Remarkably, concerted evolution is so well orchestrated that even transposable elements that insert into a large fraction of the rRNA genes appear to have little effect on the process. Finally, we summarize data that suggest that recombination in the rDNA locus of higher eukaryotes is sufficiently frequent to monitor changes within a few generations.

Adaptive evolution of multiple-variable exons and structural diversity of drug-metabolizing enzymes

BACKGROUND

The human genome contains a large number of gene clusters with multiple-variable-first exons, including the drug-metabolizing UDP glucuronosyltransferase (UGT1) and I-branching beta-1,6-N-acetylglucosaminyltransferase (GCNT2, also known as IGNT) clusters, organized in a tandem array, similar to that of the protocadherin (PCDH), immunoglobulin (IG), and T-cell receptor (TCR) clusters. To gain insight into the evolutionary processes that may have shaped their diversity, we performed comprehensive comparative analyses for vertebrate multiple-variable-first-exon clusters.

RESULTS

We found that there are species-specific variable-exon duplications and mutations in the vertebrate Ugt1, Gcnt2, and Ugt2a clusters and that their variable and constant genomic organizations are conserved and vertebrate-specific. In addition, analyzing the complete repertoires of closely-related Ugt2 clusters in humans, mice, and rats revealed extensive lineage-specific duplications. In contrast to the Pcdh gene clusters, gene conversion does not play a predominant role in the evolution of the vertebrate Ugt1, Gcnt2 and Ugt2 gene clusters. Thus, their tremendous diversity is achieved through "birth-and-death" evolution. Comparative analyses and homologous modeling demonstrated that vertebrate UGT proteins have similar three-dimensional structures each with N-terminal and C-terminal Rossmann-fold domains binding acceptor and donor substrates, respectively. Molecular docking experiments identified key residues in donor and acceptor recognition and provided insight into the catalytic mechanism of UGT glucuronidation, suggesting the human UGT1A1 residue histidine 39 (H39) as a general base and the residue aspartic acid 151 (D151) as an important electron-transfer helper. In addition, we identified four hypervariable regions in the N-terminal Rossmann domain that form an acceptor-binding pocket. Finally, analyzing patterns of nonsynonymous and synonymous nucleotide substitutions identified codon sites that are subject to positive Darwinian selection at the molecular level. These diversified residues likely play an important role in recognition of myriad xenobiotics and endobiotics.

CONCLUSION

Our results suggest that enormous diversity of vertebrate multiple variable first exons is achieved through birth-and-death evolution and that adaptive evolution of specific codon sites enhances vertebrate UGT diversity for defense against environmental agents. Our results also have interesting implications regarding the staggering molecular diversity required for chemical detoxification and drug clearance.

Comparative genomic analysis of a mammalian beta-defensin gene cluster

Comparative genomic analyses have yielded valuable insights into conserved and divergent aspects of gene function, regulation, and evolution. Herein, we describe the characterization of a mouse beta-defensin gene cluster locus on chromosome 2F6. In addition, we present the evolutionary analysis of this cluster and its human, rhesus, and rat orthologs. Expression analysis in mouse revealed the occurrence of defensin cluster transcripts in multiple tissues, with the highest abundance in the urogenital tract. Molecular evolutionary analysis suggests that this cluster originated by a series of duplication events, and by positive selection occurring even after the rodent-primate split. In addition, the constraints analysis showed higher positive selection in rodents than in primates, especially distal to the six-cysteine array. Positive selection in the evolution of these defensins may relate not only to the evolving enhancement of ancestral host defense but also to functional innovations in reproduction. The multiplicity of defensins and their preferential overexpression in the urogenital tract indicate that defensins function in the protection and maintenance of fertility.

Evolution of hyperactive, repetitive antifreeze proteins in beetles

Some organisms that experience subzero temperatures, such as insects, fish, bacteria, and plants, synthesize antifreeze proteins (AFPs) that adsorb to surfaces of nascent ice crystals and inhibit their growth. Although some AFPs are globular and nonrepetitive, the majority are repetitive in both sequence and structure. In addition, they are frequently encoded by tandemly arrayed, multigene families. AFP isoforms from the mealworm beetle, Tenebrio molitor, are extremely potent and inhibit ice growth at temperatures below -5 degrees C. They contain a 12-amino acid repeat with the sequence TCTxSxxCxxAx, each of which makes up one coil of the beta-helix structure. TxT motifs are arrayed to form the ice-binding surface in all three known insect AFPs: the homologous AFPs from the two beetles, T. molitor and Dendroides canadensis, and the nonhomologous AFP from the spruce budworm, Choristoneura fumiferana. In this study, we have obtained the cDNA and genomic sequences of additional T. molitor isoforms. They show variation in the number of repeats (from 6 to 10) which can largely be explained by recombination at various TCT motifs. In addition, phylogenetic comparison of the AFPs from the two beetles suggests that gene loss and amplification may have occurred after the divergence of these species. In contrast to a previous study suggesting that T. molitor genes have undergone positive Darwinian selection (selection for heterogeneity), we propose that the higher than expected ratio of nonsynonymous-to-synonymous substitutions might result from selection for higher AT content in the third codon position.

Evolutionary dynamics of olfactory receptor genes in Drosophila species

Olfactory receptor (OR) genes are of vital importance for animals to find food, identify mates, and avoid dangers. In mammals, the number of OR genes is large and varies extensively among different orders, whereas, in insects, the extent of interspecific variation appears to be small, although only a few species have been studied. To understand the evolutionary changes of OR genes, we identified all OR genes from 12 Drosophila species, of which the evolutionary time is roughly equivalent to that of eutherian mammals. The results showed that all species examined have similar numbers ( approximately 60) of functional OR genes. Phylogenetic analysis indicated that the ancestral species also had similar numbers of genes, but there were frequent gains and losses of genes that occurred in each evolutionary lineage. It appears that tandem duplication and random inactivation of duplicate genes are the major factors of gene number change. However, chromosomal rearrangements have contributed to the establishment of genome-wide distribution of OR genes. These results suggest that the repertoire of OR genes in Drosophila has been quite stable compared with the mammalian genes. The difference in evolutionary pattern between Drosophila and mammals can be explained partly by the differences of gene expression mechanisms and partly by the environmental and behavioral differences.