Splicing profile based protein categorization between human and mouse genomes by use of the DDBJ Web services

In one scenario of gene evolution, exon shuffling plays a fundamental role in increasing gene diversity. This paper is an appraisal of the biological relevance of categorising proteins by their splicing profiles (exon-intron structures). The central question is whether protein function is more correlated with splicing profiles than sequence similarity, or not. To approach this question, a splicing profile similarity (SPS) index, which measures relative exon length discrepancy, was devised. Arbitrary human proteins were compared, in terms of SPS and amino acid sequence similarity, to their 1) mouse orthologues and 2) human paralogues, which epitomise functional equivalence and non-equivalence, respectively, to methodically elucidate the global relationship between a) biological function, b) splicing profile similarity, and c) sequence similarity. Protein function is more correlated with splicing profile similarity than sequence similarity as demonstrated by the fact that human-mouse orthologues (HMOs) display significantly higher splicing profile similarity than do human-human paralogues (HHPs), despite the mutual sequence similarity between these two categories. This finding indicates that splicing profile-based protein categorisation is biologically meaningful.     

Domain shuffling and the evolution of vertebrates

The evolution of vertebrates has included a number of important events: the development of cartilage, the immune system, and complicated craniofacial structures. Here, we examine domain shuffling as one of the mechanisms that contributes novel genetic material required for vertebrate evolution. We mapped domain-shuffling events during the evolution of deuterostomes with a focus on how domain shuffling contributed to the evolution of vertebrate- and chordate-specific characteristics. We identified ∼1000 new domain pairs in the vertebrate lineage, including ∼100 that were shared by all seven of the vertebrate species examined. Some of these pairs occur in the protein components of vertebrate-specific structures, such as cartilage and the inner ear, suggesting that domain shuffling made a marked contribution to the evolution of vertebrate-specific characteristics. The evolutionary history of the domain pairs is traceable; for example, the Xlink domain of aggrecan, one of the major components of cartilage, was originally utilized as a functional domain of a surface molecule of blood cells in protochordate ancestors, and it was recruited by the protein of the matrix component of cartilage in the vertebrate ancestor. We also identified genes that were created as a result of domain shuffling in ancestral chordates. Some of these are involved in the functions of chordate structures, such as the endostyle, Reissner''s fiber of the neural tube, and the notochord. Our analyses shed new light on the role of domain shuffling, especially in the evolution of vertebrates and chordates.The question of how novel structures are created by changing genetic information is one of the most challenging issues in evolutionary biology. It is generally accepted that the morphological features of various multicellular animals are built on a common set of genes. Carroll et al. (2001) and Davidson (2006) proposed that the novel features emerged as a result of altered gene expression patterns. Novel genetic material has also contributed to the evolution of novel structures. Much attention has been focused on gene duplications as a mechanism for the evolution of novel genetic material. Particularly in the case of vertebrate evolution, whole genome duplications that occurred twice in ancestral vertebrates have been regarded as the main force driving the evolution of novel structures (Holland et al. 1994). As an additional mechanism, we focus here on domain shuffling (Chothia et al. 2003; Babushok et al. 2007). Several different molecular mechanisms for domain shuffling have been proposed. Since the domains are often correlated with exon boundaries, exon shuffling is believed to be one of the major forces driving domain shuffling (Liu and Grigoriev 2004). In addition, the introns at the boundaries of domains show a marked excess of symmetrical phase combinations, and, consequently, these domains may be inserted without changing the reading frame of ancestral genes (Kaessmann et al. 2002; Vibranovski et al. 2005). Some other mechanisms might have been involved in domain shuffling, such as the simple fusion of genes and recruitment of mobile elements (Babushok et al. 2007; Ekman et al. 2007).In contrast to thorough investigation of the mechanistic aspects of domain shuffling, the contribution of the new genes created by domain shuffling to the evolution of phenotype has not been sufficiently explored. Because domain shuffling occurs more frequently in metazoan lineages than in unicellular organisms, suggesting that domain shuffling played an important role in the evolution of multicellularity (Patthy 2003; Ekman et al. 2007), we reasoned that domain shuffling also contributed to the elaboration of metazoan body plans. In this study, we comprehensively investigated the domain shuffling events during deuterostome evolution. The complete amphioxus genome sequence provides sufficient deuterostome genomic sequences (Dehal et al. 2002; Sea Urchin Genome Sequencing Consortium 2006; Putnam et al. 2008) to map the domain shuffling events accurately. We listed gene models that were created by domain shuffling in the ancestors of vertebrates, and examined how domain shuffling contributed to the evolution of new genes for vertebrate-specific characteristics. The contribution of domain shuffling was also examined in the evolution of chordates.     

The MHC Big Bang

Summary: The human Major Histocompatibility Complex (MHC) shares similarities with three other chromosome regions in human. This could be the vestige of ancestral large scale duplications. We discuss here the possibility i) that these duplications occurred during two rounds of tetraploidization supposed to have taken place during chordate evolution before the jawed vertebrate radiation, and ii) that one of the quadruplicate regions, relaxed of functional constraints, gave rise to the vertebrate MHC by a quick round of gene cis-duplication and cis-exon shuffling. These different rounds of cis-duplications and exon shufflings allowed the emergence of new genes participating in novel biological functions i.e. adaptive immune responses, Cis-duplications and cis-exon shufflings are ongoing processes in the evolution of some of these genes in this region as they have occurred and were fixed at different times and in different lineages during vertebrate evolution. In contrast, other genes within the MHC have remained stable since the emergence of jawed vertebrates.     

Origin and evolution of new exons in rodents

Gene number difference among organisms demonstrates that new gene origination is a fundamental biological process in evolution. Exon shuffling has been universally observed in the formation of new genes. Yet to be learned are the ways new exons originate and evolve, and how often new exons appear. To address these questions, we identified 2695 newly evolved exons in the mouse and rat by comparing the expressed sequences of 12,419 orthologous genes between human and mouse, using 743,856 pig ESTs as the outgroup. The new exon origination rate is about 2.71 x 10(-3) per gene per million years. These new exons have markedly accelerated rates both of nonsynonymous substitutions and of insertions/deletions (indels). A much higher proportion of new exons have K(a)/K(s) ratios >1 (where K(a) is the nonsynonymous substitution rate and K(s) is the synonymous substitution rate) than do the old exons shared by human and mouse, implying a role of positive selection in the rapid evolution. The majority of these new exons have sequences unique in the genome, suggesting that most new exons might originate through "exonization" of intronic sequences. Most of the new exons appear to be alternative exons that are expressed at low levels.     

Gene duplications as a recurrent theme in the evolution of the human pseudoautosomal region 1: isolation of the gene ASMTL

We have isolated a novel gene, ASMTL (acetylserotonin methytransferase- like ), in the pseudoautosomal region (PAR1) on the human sex chromosomes. ASMTL represents a unique fusion product of two different full-length genes of different evolutionary origin and function. One part is homologous to the bacterial maf/orfE genes. The other part shows significant homology to the entire open reading frame of the previously described pseudoautosomal gene ASMT, encoding the enzyme catalysing the last step in the synthesis of melatonin. We have also detected the identity of one exon (1A) of ASMT to exon 3 in yet another pseudoautosomal gene, XE7. The data presented suggest that exon duplication and exon shuffling as well as gene fusion may represent common characteristics in the pseudoautosomal region.      

The mouse neurological mutant flailer expresses a novel hybrid gene derived by exon shuffling between Gnb5 and Myo5a

Exon shuffling is thought to be an important mechanism for evolution of new genes. Here we show that the mouse neurological mutation flailer (flr) expresses a novel gene that combines the promoter and first two exons of guanine nucleotide binding protein beta 5 (Gnb5) with the C-terminal exons of the closely linked Myosin 5A (MyoVA) gene (Myo5a). The flailer protein, which is expressed predominantly in brain, contains the N-terminal 83 amino acids of Gnb5 fused in-frame with the C-terminal 711 amino acids of MyoVA, including the globular tail domain that binds organelles for intracellular transport. Biochemical and genetic studies indicate that the flailer protein competes with wild-type MyoVA in vivo, preventing the localization of smooth endoplasmic reticulum vesicles in the dendritic spines of cerebellar Purkinje cells. The flailer protein thus has a dominant-negative mechanism of action with a recessive mode of inheritance due to the dependence of competitive binding on the ratio between mutant and wild-type proteins. The chromosomal arrangement of Myo5a upstream of Gnb5 is consistent with non-homologous recombination as the mutational mechanism. To our knowledge, flailer is the first example of a mammalian mutation caused by germ line exon shuffling between unrelated genes.     

Complement factor H: sequence analysis of 221 kb of human genomic DNA containing the entire fH, fHR-1 and fHR-3 genes

Complement factor H (fH) is a member of a family of proteins involved in the regulation of complement activation (RCA). These proteins share a common structural motif, the Short Consensus Repeat (SCR), which is structurally conserved among related genes and between phylogenetically divergent species. fH is composed of 20 such SCRs and a variety of biological functions have been localised to specific SCR domains. The majority of individual SCRs identified are encoded by single exons, and processes such as gene conversion, duplication and exon shuffling have been implicated in the evolution and genomic radiation of SCR-encoding genes. We have analysed two GenBank sequence entries relating to two overlapping PAC clones sequenced at the Sanger Centre which contain the entire human fH gene and two adjacent fH-related (fHR) genes, fHR-1 and fHR-3. Here, we report the detailed analysis of the assembled 221 kb of contiguous, ungapped genomic sequence from human chromosome 1q32, in part employing the RUMMAGE-DP automated annotation tool. Genomic duplications involving fH and fHR exons were identified and Alu/L1 repeat dating established that the duplications occurred after the separation of rodent and primate lineages. The analysis indicates that retrotransposition as well as single and multiple exon duplication events are likely to have been involved in SCR radiation and RCA gene evolution, facilitated by conservation of splice-phasing and the single-exon, single-SCR nature of the encoded domains.     

How are exons encoding transmembrane sequences distributed in the exon–intron structure of genes?

The exon–intron structure of eukaryotic genes raises a question about the distribution of transmembrane regions in membrane proteins. Were exons that encode transmembrane regions formed simply by inserting introns into preexisting genes or by some kind of exon shuffling? To answer this question, the exon‐per‐gene distribution was analyzed for all genes in 40 eukaryotic genomes with a particular focus on exons encoding transmembrane segments. In 21 higher multicellular eukaryotes, the percentage of multi‐exon genes (those containing at least one intron) within all genes in a genome was high (>70%) and with a mean of 87%. When genes were grouped by the number of exons per gene in higher eukaryotes, good exponential distributions were obtained not only for all genes but also for the exons encoding transmembrane segments, leading to a constant ratio of membrane proteins independent of the exon‐per‐gene number. The positional distribution of transmembrane regions in single‐pass membrane proteins showed that they are generally located in the amino or carboxyl terminal regions. This nonrandom distribution of transmembrane regions explains the constant ratio of membrane proteins to the exon‐per‐gene numbers because there are always two terminal (i.e., the amino and carboxyl) regions – independent of the length of sequences.     

Identification of peptide ligands to the chemokine receptor CCR5 and their maturation by gene shuffling

The determination of protein-protein interactions and their role in diverse pathophysiological processes is a promising approach to the identification of molecules of therapeutic potential. This paper describes the identification of peptidic CCR5 receptor ligands as potential drug leads against HIV-1 infection using in vitro evolution based on phage display. A phage-displayed peptide library was used to select for anti-CCR5 peptide. Further in vitro evolution of the peptide by exon shuffling was performed to identify peptides with optimized characteristics for CCR5 receptor. This peptide inhibited HIV coreceptor activity in a cell fusion assay with an IC50 of 5 microM. It did not exhibit either agonistic or antagonistic activity on CCR5 in the concentration range used. To our knowledge, this is a first report that describes the identification of peptide ligands specific to the CCR5 receptor from a phage-displayed library and the maturation of the selected peptide sequence by gene shuffling.     

Phylogenetically Older Introns Strongly Correlate With Module Boundaries in Ancient Proteins

The hypothesis that some (but not all) introns were used to construct ancient genes by exon shuffling of modules at the earliest stages of evolution is supported by the finding of an excess of phase-zero intron positions in the boundary regions of such modules in 276 ancient proteins (defined as common to eukaryotes and prokaryotes). Here we show further that as phase-zero intron positions are shared by distant taxa, and thus are truly phylogenetically ancient, their excess in the boundaries becomes greater, rising to an 80% excess if shared by four out of the five taxa: vertebrates, invertebrates, fungi, plants, and protists.     

Genomic analysis of the nuclear receptor family: new insights into structure, regulation, and evolution from the rat genome

Completion of the Rattus norvegicus genome sequence enabled a global inventory and analysis of the nuclear receptors (NRs) in three mammalian species. Forty-nine NR members were found in mouse, 48 in human. Forty-seven were found in the rat, with gaps at the locations expected for the other two. Pairwise comparisons of their distribution in rat, mouse, and human identified 11 syntenic NR gene blocks, including three small clusters of two or three closely related genes, each spanning 40 kb to 1700 kb. The exon structure of the ligand-binding domain suggests that exon shuffling has played a role in the evolution of this family. An invariant splice junction in all members of the NR family except LXRbeta suggests a functional role for the intron. The ligand-binding domains of PXR and CAR are among the most divergent in the family. Their higher nucleotide substitution rates may be related to the central role played by these two NRs in the metabolism of the foreign compounds and may have resulted from limited positive selection.     

On the origin of new genes in Drosophila

Several mechanisms have been proposed to account for the origination of new genes. Despite extensive case studies, the general principles governing this fundamental process are still unclear at the whole-genome level. Here, we unveil genome-wide patterns for the mutational mechanisms leading to new genes and their subsequent lineage-specific evolution at different time nodes in the Drosophila melanogaster species subgroup. We find that (1) tandem gene duplication has generated approximately 80% of the nascent duplicates that are limited to single species (D. melanogaster or Drosophila yakuba); (2) the most abundant new genes shared by multiple species (44.1%) are dispersed duplicates, and are more likely to be retained and be functional; (3) de novo gene origination from noncoding sequences plays an unexpectedly important role during the origin of new genes, and is responsible for 11.9% of the new genes; (4) retroposition is also an important mechanism, and had generated approximately 10% of the new genes; (5) approximately 30% of the new genes in the D. melanogaster species complex recruited various genomic sequences and formed chimeric gene structures, suggesting structure innovation as an important way to help fixation of new genes; and (6) the rate of the origin of new functional genes is estimated to be five to 11 genes per million years in the D. melanogaster subgroup. Finally, we survey gene frequencies among 19 globally derived strains for D. melanogaster-specific new genes and reveal that 44.4% of them show copy number polymorphisms within a population. In conclusion, we provide a panoramic picture for the origin of new genes in Drosophila species.     

Signatures of domain shuffling in the human genome

To elucidate the role of exon shuffling in shaping the complexity of the human genome/proteome, we have systematically analyzed intron phase distributions in the coding sequence of human protein domains. We found that introns at the boundaries of domains show high excess of symmetrical phase combinations (i.e., 0-0, 1-1, and 2-2), whereas nonboundary introns show no excess symmetry. This suggests that exon shuffling has primarily involved rearrangement of structural and functional domains as a whole. Furthermore, we found that domains flanked by phase 1 introns have dramatically expanded in the human genome due to domain shuffling and that 1-1 symmetrical domains and domain families are nonrandomly distributed with respect to their age. The predominance and extracellular location of 1-1 symmetrical domains among domains specific to metazoans suggests that they are associated with the rise of multicellularity. On the other hand, 0-0 symmetrical domains tend to be over-represented among ancient protein domains that are shared between the eukaryotic and prokaryotic kingdoms, which is compatible with the suggestion of primordial domain shuffling in the progenote. To see whether the human data reflect general genomic patterns of metazoans, similar analyses were done for the nematode Caenorhabditis elegans. Although the C. elegans data generally concur with the human patterns, we identified fewer intron-bounded domains in this organism, consistent with the lower complexity of C. elegans genes. [The following individuals kindly provided reagents, samples, or unpublished information as indicated in the paper: Z. Gu and R. Stevens.]     

Multi-alignment of orthologous genome regions in five species provides new insights into the evolutionary make-up of mammalian genomes

Evidence has shown that bacterial genomes have undergone random shuffling of genomic elements consisting of one to two genes. In order to delineate such genome-shuffling events in mammals, we constructed a high-resolution map of Sus scrofa chromosome 3 (SSC3) with a total of 116 genes/markers. Alignment of this pig map to orthologous regions in human, dog, mouse and rat led to the identification of 31 provisional conserved ancestral blocks (CABs) in these five species. Among them, only 3 CABs (<10%) had one gene, indicating that one-gene shuffling is not frequent in mammals. The sizes of CABs vary significantly within a species, but each may be relatively consistent in different species with a scale to species-genome evolution. The type and frequency of rearrangement events that takes place, either intra- or interchromosomal, depends on the evolutionary regions and species under comparison. Characterization of 36 tentative breakpoint regions flanking these 31 CABs indicated that they occupied ∼43 Mb in length and featured genome deserts, gene duplications, and birth/death of species-specific genes in humans. Identification of CABs provides an alternative for further determination of the evolutionary make-up of mammalian genomes.     

Variety of antimicrobial peptides in the Bombina maxima toad and evidence of their rapid diversification

Antimicrobial peptides secreted by the skin of many amphibians play an important role in innate immunity. From two skin cDNA libraries of two individuals of the Chinese red belly toad (Bombina maxima), we identified 56 different antimicrobial peptide cDNA sequences, each of which encodes a precursor peptide that can give rise to two kinds of antimicrobial peptides, maximin and maximin H. Among these cDNA, we found that the mean number of nucleotide substitution per non-synonymous site in both the maximin and maximin H domains significantly exceed the mean number of nucleotide substitution per synonymous site, whereas the same pattern was not observed in other structural regions, such as the signal and propiece peptide regions, suggesting that these antimicrobial peptide genes have been experiencing rapid diversification driven by Darwinian selection. We cloned and sequenced seven genes amplified from skin or liver genomic DNA. These genes have three exons and share the same gene structure, in which both maximin and maximin H are encoded by the third exon. This suggests that alternative splicing and somatic recombination are less likely to play a role in creating the diversity of maximins and maximin Hs. The gene trees based on different domain regions revealed that domain shuffling or gene conversion among these genes might have happened frequently.     

Pathogenomics of mobile genetic elements of toxigenic bacteria

The growing knowledge of genetic diversity and whole genome organization in bacteria shows that pathogenicity islands (PAIs) represent a subtype of a more general genetic element, termed genomic island (GEI), which is widespread among pathogenic and non-pathogenic microbes. These findings mirror the importance of horizontal gene transfer, genome reduction and recombination events as fundamental mechanisms involved in evolution of bacterial variants. GEIs are part of the flexible gene pool and carry selfish genes, but also determinants which may be beneficial under certain conditions thus increasing bacterial fitness and consequently their survival or transmission. In this review, we focus on the role of mobile genetic elements that may also contain toxin-encoding genes for genome variability and evolution of bacteria.     

Differential accumulation of mutations localized in particular domains of the mucin genes expressed in the vertebrate host stage of Trypanosoma cruzi

The surface of Trypanosoma cruzi is covered by mucin-type glycoproteins involved in parasite protection, attachment and immunoevasion. The gene family coding for the mucins expressed by the parasite in the vertebrate host, named TcMUC, is composed of several hundred members and presents high variability. The genes encoding mucins expressed in the insect-dwelling parasite stages are part of a much more homogeneous family, named TcSMUG. Here, we addressed the organization and evolution of physically linked T. cruzi mucin genes by sequencing large chromosomal fragments containing these genes. Specific accumulation of mutations was restricted to particular domains of TcMUC genes, showing that these regions have, or have had, an accelerated evolution rate. Sequence analysis of several TcMUC genes allowed for the identification of members sharing features of TcMUC I and II, thus evidencing that one group of genes was generated from the other. The highly conserved intergenic regions of both TcMUC and TcSMUG families contained TG-rich microsatellites that were not present in unrelated genes in the cosmids, suggesting a role for homologous recombination in shuffling and/or amplification of T. cruzi mucin genes. The comparison of putative homologous TcMUC II genes from different strains of T. cruzi showed that their central variable domains are conserved. This conservation was always higher at the DNA level suggesting positive selection in these particular regions of TcMUC II genes.     

Novel immune-type receptor genes

Summary: Novel immune-type receptor ( NITR ) genes, which initially were identified in the Southern pufferfish ( Spheroides nephelus ), encode products which consist of an extracellular variable (V) and V-like C2 (V/C2) domain, a transmembrane region, and a cytoplasmic tail, which typically possesses an immunoreceptor tyrosine-based inhibition motif (ITIM). Multiple NITR genes have been identified in close, contiguous chromosomal linkage. The V regions of NITRs resemble prototypic forms defined for immunoglobulin (Ig) and T-cell antigen receptor (TCR), are present in multiple families and exhibit regionalized variation in sequence, which also occurs in Ig and TCR. Comparisons of exons encoding transmembrane and cytoplasmic regions of multiple NITRs suggest that exon shuffling has factored in the diversification of the NITR gene complex. Zebrafish ( Danio rerio ) NITRs exhibit many of these characteristics. NITRs that have been identified in additional species of bony fish demonstrate additional variation in the number of extracellular domains as well as in the presence of intramembranous charged residues, cytoplasmic tails and ITIMs. The presence in NITRs of V regions that are related closely to those found in Ig and TCR, as well as regulatory motifs and other structural features that are characteristic of immune inhibitory receptors encoded at the leukocyte receptor cluster, suggests that the NITRs are representative of an integral stage in the evolution of innate and adaptive immune function.
This research was supported by grants AI23338 to GWL and GM20231 to JAY from the National Institutes of Health as well as a grant from The Pediatric Cancer Foundation, Inc. to GWL.