Late changes in spliceosomal introns define clades in vertebrate evolution

The promise and pitfalls of synteny in phylogenomics

Reconstructing the tree of life remains a central goal in biology. Early methods, which relied on small numbers of morphological or genetic characters, often yielded conflicting evolutionary histories, undermining confidence in the results. Investigations based on phylogenomics, which use hundreds to thousands of loci for phylogenetic inquiry, have provided a clearer picture of life's history, but certain branches remain problematic. To resolve difficult nodes on the tree of life, 2 recent studies tested the utility of synteny, the conserved collinearity of orthologous genetic loci in 2 or more organisms, for phylogenetics. Synteny exhibits compelling phylogenomic potential while also raising new challenges. This Essay identifies and discusses specific opportunities and challenges that bear on the value of synteny data and other rare genomic changes for phylogenomic studies. Synteny-based analyses of highly contiguous genome assemblies mark a new chapter in the phylogenomic era and the quest to reconstruct the tree of life.

Molecular Barcoding: A Tool to Guarantee Correct Seafood Labelling and Quality and Preserve the Conservation of Endangered Species

The recent increase in international fish trade leads to the need for improving the traceability of fishery products. In relation to this, consistent monitoring of the production chain focusing on technological developments, handling, processing and distribution via global networks is necessary. Molecular barcoding has therefore been suggested as the gold standard in seafood species traceability and labelling. This review describes the DNA barcoding methodology for preventing food fraud and adulteration in fish. In particular, attention has been focused on the application of molecular techniques to determine the identity and authenticity of fish products, to discriminate the presence of different species in processed seafood and to characterize raw materials undergoing food industry processes. In this regard, we herein present a large number of studies performed in different countries, showing the most reliable DNA barcodes for species identification based on both mitochondrial (COI, cytb, 16S rDNA and 12S rDNA) and nuclear genes. Results are discussed considering the advantages and disadvantages of the different techniques in relation to different scientific issues. Special regard has been dedicated to a dual approach referring to both the consumer's health and the conservation of threatened species, with a special focus on the feasibility of the different genetic and genomic approaches in relation to both scientific objectives and permissible costs to obtain reliable traceability.

Trans‐specific polymorphism and the convergent evolution of supertypes in major histocompatibility complex class II genes in darters ( Etheostoma )

Major Histocompatibility Complex (MHC) genes are one of the most polymorphic gene groups known in vertebrates. MHC genes also exhibit allelic variants that are shared among taxa, referred to as trans‐specific polymorphism (TSP). The role that selection plays in maintaining such high diversity within species, as well as TSP, is an ongoing discussion in biology. In this study, we used deep‐sequencing techniques to characterize MHC class IIb gene diversity in three sympatric species of darters. We found at least 5 copies of the MHC gene in darters, with 126 genetic variants encoding 122 unique amino acid sequences. We identified four supertypes based on the binding properties of proteins encoded by the sequences. Although each species had a unique pool of variants, many variants were shared between species pairs and across all three species. Phylogenetic analysis showed that the variants did not group together monophyletically based on species identity or on supertype. An expanded phylogenetic analysis showed that some darter alleles grouped together with alleles from other percid fishes. Our findings show that TSP occurs in darters, which suggests that balancing selection is acting at the genotype level. Supertypes, however, are most likely evolving convergently, as evidenced by the fact that alleles do not form monophyletic groups based on supertype. Our research demonstrates that selection may be acting differently on MHC genes at the genotype and supertype levels, selecting for the maintenance of high genotypic diversity while driving the convergent evolution of similar MHC phenotypes across different species.

OUP accepted manuscript

Homology of highly divergent genes often cannot be determined from sequence similarity alone. For example, we recently identified in the aphid Hormaphis cornu a family of rapidly evolving bicycle genes, which encode novel proteins implicated as plant gall effectors, and sequence similarity search methods yielded few putative bicycle homologs in other species. Coding sequence-independent features of genes, such as intron-exon boundaries, often evolve more slowly than coding sequences, however, and can provide complementary evidence for homology. We found that a linear logistic regression classifier using only structural features of bicycle genes identified many putative bicycle homologs in other species. Independent evidence from sequence features and intron locations supported homology assignments. To test the potential roles of bicycle genes in other aphids, we sequenced the genome of a second gall-forming aphid, Tetraneura nigriabdominalis and found that many bicycle genes are strongly expressed in the salivary glands of the gall forming foundress. In addition, bicycle genes are strongly overexpressed in the salivary glands of a non-gall forming aphid, Acyrthosiphon pisum, and in the non-gall forming generations of H. cornu. These observations suggest that Bicycle proteins may be used by multiple aphid species to manipulate plants in diverse ways. Incorporation of gene structural features into sequence search algorithms may aid identification of deeply divergent homologs, especially of rapidly evolving genes involved in host-parasite interactions.

Evolutionary History and Taxonomic Reappraisal of Coral Reef Rabbitfishes (Siganidae): Patterns of Lineage Diversification and Speciation

Simple Summary

Herbivorous fish are recognized as being ecologically important to the health and survival of coral reef ecosystems because they remove algal turfs growing on corals. Apart from being one of the major components of herbivorous fish communities, rabbitfish are also characterized by possessing rabbit-like mouths. A total of 29 species of rabbitfish are confined to a single genus, Siganus, fish that are highly sought after for the aquarium trade and for food by humans. Natural hybridization between some species that have parapatric distributions across the Indo-West Pacific region may have homogenized their genotypic and morphological features. Relatively little is known, however, about how environmental factors may affect phylogenetic relationships among these siganid species. Based on sequencing of eight siganid species collected from the South China Sea and meta-analysis of sequences from ten siganid species retrieved from the NCBI database, we applied an integrated morphological–molecular approach to elucidate phylogenetic relationships and demographic histories of these species. Our results highlight that diversification and speciation of siganid species were influenced by a series of paleo-climatic events, changes to natural geographical distributions, and associated environmental changes. The target species were differentiated by body shape, and two morphometric parameters, notably body depth and snout length. Our results provide considerable baseline knowledge for strategizing improvement of both breeding and conservation programs for rabbitfish.

Abstract

Rabbitfish (Siganidae) are coral reef fish that are distributed across diverse habitats that include estuaries, mangroves, reefs, and even seaweed mats. Given their ecological diversity and natural widespread distributions across the Indo-Pacific region, we were interested to investigate the evolutionary history of this group and patterns of divergence that have contributed to their present-day distributions. In the present study, samples were collected from the South China Sea to study taxonomic and phylogenetic relationships, and divergence times. We investigated the taxonomic relationships among modern rabbitfish species, reconstructed their molecular phylogeny, and estimated divergence times among selected lineages based on a fragment of the mtDNA cytochrome oxidase I (COI) and sequences of the nuclear rhodopsin retrogene (RHO). Our results indicate that modern rabbitfish likely originated in the Indo-West Pacific during the late Eocene [37.4 million years ago (mya)], following which they diverged into three major clades during the Pliocene/Pleistocene. Subsequent diversification and origins of the majority of siganids may likely be associated with episodes of paleo-oceanographic events, including greenhouse and glaciation events (Eocene–Miocene) as well as major plate tectonic events (Pliocene–Pleistocene). Some modern siganid species may naturally hybridize with congeneric species where their geographical ranges overlap. A comprehensive taxonomic analysis revealed that the phylogeny of Siganidae (cladogenesis of Clades I, II, and III) is characterized by divergence in several external morphological characters and morphometric parameters. Our study demonstrates that morphological characteristics, geographical heterogeneity, and environmental change have contributed to siganids’ historical diversification.

Evolutionary history of teleost intron-containing and intron-less rhodopsin genes

Recent progress in whole genome sequencing has revealed that animals have various kinds of opsin genes for photoreception. Among them, most opsin genes have introns in their coding regions. However, it has been known for a long time that teleost retinas express intron-less rhodopsin genes, which are presumed to have been formed by retroduplication from an ancestral intron-containing rhodopsin gene. In addition, teleosts have an intron-containing rhodopsin gene (exo-rhodopsin) exclusively for pineal photoreception. In this study, to unravel the evolutionary origin of the two teleost rhodopsin genes, we analyzed the rhodopsin genes of non-teleost fishes in the Actinopterygii. The phylogenetic analysis of full-length sequences of bichir, sturgeon and gar rhodopsins revealed that retroduplication of the rhodopsin gene occurred after branching of the bichir lineage. In addition, analysis of the tissue distribution and the molecular properties of bichir, sturgeon and gar rhodopsins showed that the abundant and exclusive expression of intron-containing rhodopsin in the pineal gland and the short lifetime of its meta II intermediate, which leads to optimization for pineal photoreception, were achieved after branching of the gar lineage. Based on these results, we propose a stepwise evolutionary model of teleost intron-containing and intron-less rhodopsin genes.

Rhodopsin gene evolution in early teleost fishes

Rhodopsin mediates an essential step in image capture and is tightly associated with visual adaptations of aquatic organisms, especially species that live in dim light environments (e.g., the deep sea). The rh1 gene encoding rhodopsin was formerly considered a single-copy gene in genomes of vertebrates, but increasing exceptional cases have been found in teleost fish species. The main objective of this study was to determine to what extent the visual adaptation of teleosts might have been shaped by the duplication and loss of rh1 genes. For that purpose, homologous rh1/rh1-like sequences in genomes of ray-finned fishes from a wide taxonomic range were explored using a PCR-based method, data mining of public genetic/genomic databases, and subsequent phylogenomic analyses of the retrieved sequences. We show that a second copy of the fish-specific intron-less rh1 is present in the genomes of most anguillids (Elopomorpha), Hiodon alosoides (Osteoglossomorpha), and several clupeocephalan lineages. The phylogenetic analysis and comparisons of alternative scenarios for putative events of gene duplication and loss suggested that fish rh1 was likely duplicated twice during the early evolutionary history of teleosts, with one event coinciding with the hypothesized fish-specific genome duplication and the other in the common ancestor of the Clupeocephala. After these gene duplication events, duplicated genes were maintained in several teleost lineages, whereas some were secondarily lost in specific lineages. Alternative evolutionary schemes of rh1 and comparison with previous studies of gene evolution are also reviewed.

Rhodopsin gene copies in Japanese eel originated in a teleost-specific genome duplication

BACKGROUND

Gene duplication is considered important to increasing the genetic diversity in animals. In fish, visual pigment genes are often independently duplicated, and the evolutionary significance of such duplications has long been of interest. Eels have two rhodopsin genes (rho), one of which (freshwater type, fw-rho) functions in freshwater and the other (deep-sea type, ds-rho) in marine environments. Hence, switching of rho expression in retinal cells is tightly linked with eels' unique life cycle, in which they migrate from rivers or lakes to the sea. These rho genes are apparently paralogous, but the timing of their duplication is unclear due to the deep-branching phylogeny. The aim of the present study is to elucidate the evolutionary origin of the two rho copies in eels using comparative genomics methods.

RESULTS

In the present study, we sequenced the genome of Japanese eel Anguilla japonica and reconstructed two regions containing rho by de novo assembly. We found a single corresponding region in a non-teleostean primitive ray-finned fish (spotted gar) and two regions in a primitive teleost (Asian arowana). The order of ds-rho and the neighboring genes was highly conserved among the three species. With respect to fw-rho, which was lost in Asian arowana, the neighboring genes were also syntenic between Japanese eel and Asian arowana. In particular, the pattern of gene losses in ds-rho and fw-rho regions was the same as that in Asian arowana, and no discrepancy was found in any of the teleost genomes examined. Phylogenetic analysis supports mutual monophyly of these two teleostean synteny groups, which correspond to the ds-rho and fw-rho regions.

CONCLUSIONS

Syntenic and phylogenetic analyses suggest that the duplication of rhodopsin gene in Japanese eel predated the divergence of eel (Elopomorpha) and arowana (Osteoglossomorpha). Thus, based on the principle of parsimony, it is most likely that the rhodopsin paralogs were generated through a whole genome duplication in the ancestor of teleosts, and have remained till the present in eels with distinct functional roles. Our result indicates, for the first time, that teleost-specific genome duplication may have contributed to a gene innovation involved in eel-specific migratory life cycle.

Structure and vascularization of the ventricular myocardium in Holocephali: their evolutionary significance

It was generally assumed that the ventricle of the primitive vertebrate heart was composed of trabeculated, or spongy, myocardium, supplied by oxygen-poor luminal blood. In addition, it was presumed that the mixed ventricular myocardium, consisting of a compacta and a spongiosa, and its supply through coronary arteries appeared several times throughout fish evolution. Recent work has suggested, however, that a fully vascularized, mixed myocardium may be the primitive condition in gnathostomes. The present study of the heart ventricles of four holocephalan species aimed to clarify this controversy. Our observations showed that the ventricular myocardium of Chimaera monstrosa and Harriotta raleighana consists of a very thin compacta overlying a widespread spongiosa. The ventricle of Hydrolagus affinis is composed exclusively of trabeculated myocardium. In these three species there is a well-developed coronary artery system. The main coronary artery trunks run along the outflow tract, giving off subepicardial ventricular arteries. The trabeculae of the spongiosa are irrigated by branches of the subepicardial arteries and by penetrating arterial vessels arising directly from the main coronary trunks at the level of the conoventricular junction. The ventricle of Rhinochimaera atlantica has only spongy myocardium supplied by luminal blood. Small coronary arterial vessels are present in the subepicardium, but they do not enter the myocardial trabeculae. The present findings show for the first time that in a wild living vertebrate species, specifically H. affinis, an extensive coronary artery system supplying the whole cardiac ventricle exists in the absence of a well-developed compact ventricular myocardium. This is consistent with the notion derived from experimental work that myocardial cell proliferation and coronary vascular growth rely on distinct developmental programs. Our observations, together with data in the literature on elasmobranchs, support the view that the mixed ventricular myocardium is primitive for chondrichthyans. The reduction or even lack of compacta in holocephali has to be regarded as a derived anatomical trait. Our findings also fit in with the view that the mixed myocardium was the primitive condition in gnathostomes, and that the absence of compact ventricular myocardium in different actinopterygian groups is the result of a repeated loss of such type of cardiac muscle during fish evolution.

Rapid genomic DNA changes in allotetraploid fish hybrids

Rapid genomic change has been demonstrated in several allopolyploid plant systems; however, few studies focused on animals. We addressed this issue using an allotetraploid lineage (4nAT) of freshwater fish originally derived from the interspecific hybridization of red crucian carp (Carassius auratus red var., ♀, 2n=100) × common carp (Cyprinus carpio L., ♂, 2n=100). We constructed a bacterial artificial chromosome (BAC) library from allotetraploid hybrids in the 20th generation (F20) and sequenced 14 BAC clones representing a total of 592.126 kb, identified 11 functional genes and estimated the guanine-cytosine content (37.10%) and the proportion of repetitive elements (17.46%). The analysis of intron evolution using nine orthologous genes across a number of selected fish species detected a gain of 39 introns and a loss of 30 introns in the 4nAT lineage. A comparative study based on seven functional genes among 4nAT, diploid F1 hybrids (2nF1) (first generation of hybrids) and their original parents revealed that both hybrid types (2nF1 and 4nAT) not only inherited genomic DNA from their parents, but also demonstrated rapid genomic DNA changes (homoeologous recombination, parental DNA fragments loss and formation of novel genes). However, 4nAT presented more genomic variations compared with their parents than 2nF1. Interestingly, novel gene fragments were found for the iqca1 gene in both hybrid types. This study provided a preliminary genomic characterization of allotetraploid F20 hybrids and revealed evolutionary and functional genomic significance of allopolyploid animals.

Analysis of CACTA transposases reveals intron loss as major factor influencing their exon/intron structure in monocotyledonous and eudicotyledonous hosts

Background

CACTA elements are DNA transposons and are found in numerous organisms. Despite their low activity, several thousand copies can be identified in many genomes. CACTA elements transpose using a ‘cut-and-paste’ mechanism, which is facilitated by a DDE transposase. DDE transposases from CACTA elements contain, despite their conserved function, different exon numbers among various CACTA families. While earlier studies analyzed the ancestral history of the DDE transposases, no studies have examined exon loss and gain with a view of mechanisms that could drive the changes.

Results

We analyzed 64 transposases from different CACTA families among monocotyledonous and eudicotyledonous host species. The annotation of the exon/intron boundaries showed a range from one to six exons. A robust multiple sequence alignment of the 64 transposases based on their protein sequences was created and used for phylogenetic analysis, which revealed eight different clades. We observed that the exon numbers in CACTA transposases are not specific for a host genome. We found that ancient CACTA lineages diverged before the divergence of monocotyledons and eudicotyledons. Most exon/intron boundaries were found in three distinct regions among all the transposases, grouping 63 conserved intron/exon boundaries.

Conclusions

We propose a model for the ancestral CACTA transposase gene, which consists of four exons, that predates the divergence of the monocotyledons and eudicotyledons. Based on this model, we propose pathways of intron loss or gain to explain the observed variation in exon numbers. While intron loss appears to have prevailed, a putative case of intron gain was nevertheless observed.

A novel model for development, organization, and function of gonadotropes in fish pituitary

The gonadotropins follicle-stimulating hormone (FSH) and luteinizing hormone (LH) are key regulators of the reproductive axis in vertebrates. Despite the high popularity of zebrafish as a model organism for studying reproductive functions, to date no transgenic zebrafish with labeled gonadotropes have been introduced. Using gonadotropin regulatory elements from tilapia, we generated two transgenic zebrafish lines with labeled gonadotropes. The tilapia and zebrafish regulatory sequences were highly divergent but several conserved elements allowed the tilapia promoters to correctly drive the transgenes in zebrafish pituitaries. FSH cells reacted to stimulation with gonadotropin releasing hormone by proliferating and showing increased transgene fluorescence, whereas estrogen exposure caused a decrease in cell number and transgene fluorescence. Transgene fluorescence reflected the expression pattern of the endogenous fshb gene. Ontogenetic expression of the transgenes followed typical patterns, with FSH cells appearing early in development, and LH cells appearing later and increasing dramatically in number with the onset of puberty. Our transgenic lines provide a powerful tool for investigating the development, anatomy, and function of the reproductive axis in lower vertebrates.

The vertebrate ancestral repertoire of visual opsins, transducin alpha subunits and oxytocin/vasopressin receptors was established by duplication of their shared genomic region in the two rounds of early vertebrate genome duplications

BACKGROUND

Vertebrate color vision is dependent on four major color opsin subtypes: RH2 (green opsin), SWS1 (ultraviolet opsin), SWS2 (blue opsin), and LWS (red opsin). Together with the dim-light receptor rhodopsin (RH1), these form the family of vertebrate visual opsins. Vertebrate genomes contain many multi-membered gene families that can largely be explained by the two rounds of whole genome duplication (WGD) in the vertebrate ancestor (2R) followed by a third round in the teleost ancestor (3R). Related chromosome regions resulting from WGD or block duplications are said to form a paralogon. We describe here a paralogon containing the genes for visual opsins, the G-protein alpha subunit families for transducin (GNAT) and adenylyl cyclase inhibition (GNAI), the oxytocin and vasopressin receptors (OT/VP-R), and the L-type voltage-gated calcium channels (CACNA1-L).

RESULTS

Sequence-based phylogenies and analyses of conserved synteny show that the above-mentioned gene families, and many neighboring gene families, expanded in the early vertebrate WGDs. This allows us to deduce the following evolutionary scenario: The vertebrate ancestor had a chromosome containing the genes for two visual opsins, one GNAT, one GNAI, two OT/VP-Rs and one CACNA1-L gene. This chromosome was quadrupled in 2R. Subsequent gene losses resulted in a set of five visual opsin genes, three GNAT and GNAI genes, six OT/VP-R genes and four CACNA1-L genes. These regions were duplicated again in 3R resulting in additional teleost genes for some of the families. Major chromosomal rearrangements have taken place in the teleost genomes. By comparison with the corresponding chromosomal regions in the spotted gar, which diverged prior to 3R, we could time these rearrangements to post-3R.

CONCLUSIONS

We present an extensive analysis of the paralogon housing the visual opsin, GNAT and GNAI, OT/VP-R, and CACNA1-L gene families. The combined data imply that the early vertebrate WGD events contributed to the evolution of vision and the other neuronal and neuroendocrine functions exerted by the proteins encoded by these gene families. In pouched lamprey all five visual opsin genes have previously been identified, suggesting that lampreys diverged from the jawed vertebrates after 2R.

Characterization of polymorphic microsatellite markers and genetic diversity in wild bronze featherback, Notopterus notopterus (Pallas, 1769)

Six polymorphic microsatellite DNA loci were identified in the primitive fish, bronze featherback, Notopterus notopterus for the first time and demonstrated significant population genetic structure. Out of the six primers, one primer (NN90) was specific to N. notopterus (microsatellite sequence within the RAG1 gene) and five primers were product of successful cross-species amplification. Sixty-four primers available from 3 fish species of order Osteoglossiformes and families Notopteridae and Osteoglossidae were tested to amplify homologous microsatellite loci in N. notopterus. Fifteen primer pairs exhibited successful cross-priming PCR product. However, polymorphism was detected only at five loci. To assess the significance of these six loci (including NN90) in population genetic study, 215 samples of N. notopterus from five rivers, viz Satluj, Gomti, Yamuna, Brahmaputra and Mahanadi were analyzed. The five sample sets displayed different diversity levels and observed heterozygosity ranged from 0.6036 to 0.7373. Significant genotype heterogeneity (P < 0.0001) and high FST (0.2205) over all loci indicated that the samples are not drawn from the same genepool. The identified microsatellite loci are promising for use in fine-scale population structure analysis of N. notopterus.

Can long-range PCR be used to amplify genetically divergent mitochondrial genomes for comparative phylogenetics? A case study within spiders (Arthropoda: Araneae)

The development of second generation sequencing technology has resulted in the rapid production of large volumes of sequence data for relatively little cost, thereby substantially increasing the quantity of data available for phylogenetic studies. Despite these technological advances, assembling longer sequences, such as that of entire mitochondrial genomes, has not been straightforward. Existing studies have been limited to using only incomplete or nominally intra-specific datasets resulting in a bottleneck between mitogenome amplification and downstream high-throughput sequencing. Here we assess the effectiveness of a wide range of targeted long-range PCR strategies, encapsulating single and dual fragment primer design approaches to provide full mitogenomic coverage within the Araneae (Spiders). Despite extensive rounds of optimisation, full mitochondrial genome PCR amplifications were stochastic in most taxa, although 454 Roche sequencing confirmed the successful amplification of 10 mitochondrial genomes out of the 33 trialled species. The low success rates of amplification using long-Range PCR highlights the difficulties in consistently obtaining genomic amplifications using currently available DNA polymerases optimised for large genomic amplifications and suggests that there may be opportunities for the use of alternative amplification methods.

Body plan convergence in the evolution of skates and rays (Chondrichthyes: Batoidea)

Skates, rays and allies (Batoidea) comprise more than half of the species diversity and much of the morphological disparity among chondrichthyan fishes, the sister group to all other jawed vertebrates. While batoids are morphologically well characterized and have an excellent fossil record, there is currently no consensus on the interrelationships of family-level taxa. Here we construct a resolved, robust and time-calibrated batoid phylogeny using mitochondrial genomes, nuclear genes, and fossils, sampling densely across taxa. Data partitioning schemes, biases in the sequence data, and the relative informativeness of each fossil are explored. The molecular phylogeny is largely congruent with morphology crownward in the tree, but the branching orders of major batoid groups are mostly novel. Body plan convergence appears to be widespread in batoids. A depressed, rounded pectoral disk supported to the snout tip by fin radials, common to skates and stingrays, is indicated to have been derived independently by each group, while the long, spiny rostrum of sawfishes similarly appears to be convergent with that of sawsharks, which are not batoids. The major extant batoid lineages are inferred to have arisen relatively rapidly from the Late Triassic into the Jurassic, with long stems followed by subsequent radiations in each group around the Cretaceous/Tertiary boundary. The fossil record indicates that batoids were affected with disproportionate severity by the end-Cretaceous extinction event.

Opsin gene duplication and divergence in ray-finned fish

Opsin gene sequences were first reported in the 1980s. The goal of that research was to test the hypothesis that human opsins were members of a single gene family and that variation in human color vision was mediated by mutations in these genes. While the new data supported both hypotheses, the greatest contribution of this work was, arguably, that it provided the data necessary for PCR-based surveys in a diversity of other species. Such studies, and recent whole genome sequencing projects, have uncovered exceptionally large opsin gene repertoires in ray-finned fishes (taxon, Actinopterygii). Guppies and zebrafish, for example, have 10 visual opsin genes each. Here we review the duplication and divergence events that have generated these gene collections. Phylogenetic analyses revealed that large opsin gene repertories in fish have been generated by gene duplication and divergence events that span the age of the ray-finned fishes. Data from whole genome sequencing projects and from large-insert clones show that tandem duplication is the primary mode of opsin gene family expansion in fishes. In some instances gene conversion between tandem duplicates has obscured evolutionary relationships among genes and generated unique key-site haplotypes. We mapped amino acid substitutions at so-called key-sites onto phylogenies and this exposed many examples of convergence. We found that dN/dS values were higher on the branches of our trees that followed gene duplication than on branches that followed speciation events, suggesting that duplication relaxes constraints on opsin sequence evolution. Though the focus of the review is opsin sequence evolution, we also note that there are few clear connections between opsin gene repertoires and variation in spectral environment, morphological traits, or life history traits.

A segmental genomic duplication generates a functional intron

An intron is an extended genomic feature whose function requires multiple constrained positions—donor and acceptor splice sites, a branch point, a polypyrimidine tract and suitable splicing enhancers—that may be distributed over hundreds or thousands of nucleotides. New introns are therefore unlikely to emerge by incremental accumulation of functional sub-elements. Here we demonstrate that a functional intron can be created de novo in a single step by a segmental genomic duplication. This experiment recapitulates in vivo the birth of an intron that arose in the ancestral jawed vertebrate lineage nearly half-a-billion years ago.

The appearance of a new intron that splits an exon without disrupting the corresponding peptide sequence is a rare event in vertebrate genomes. Hellsten et al. demonstrate that, under certain circumstances, a functional intron can be produced in a single step by segmental genomic duplication.

Spliceosomal intron insertions in genome compacted ray-finned fishes as evident from phylogeny of MC receptors, also supported by a few other GPCRs

BACKGROUND

Insertions of spliceosomal introns are very rare events during evolution of vertebrates and the mechanisms governing creation of novel intron(s) remain obscure. Largely, gene structures of melanocortin (MC) receptors are characterized by intron-less architecture. However, recently a few exceptions have been reported in some fishes. This warrants a systematic survey of MC receptors for understanding intron insertion events during vertebrate evolution.

METHODOLOGY/PRINCIPAL FINDINGS

We have compiled an extended list of MC receptors from different vertebrate genomes with variations in fishes. Notably, the closely linked MC2Rs and MC5Rs from a group of ray-finned fishes have three and one intron insertion(s), respectively, with conserved positions and intron phase. In both genes, one novel insertion was in the highly conserved DRY motif at the end of helix TM3. Further, the proto-splice site MAG↑R is maintained at intron insertion sites in these two genes. However, the orthologs of these receptors from zebrafish and tetrapods are intron-less, suggesting these introns are simultaneously created in selected fishes. Surprisingly, these novel introns are traceable only in four fish genomes. We found that these fish genomes are severely compacted after the separation from zebrafish. Furthermore, we also report novel intron insertions in P2Y receptors and in CHRM3. Finally, we report ultrasmall introns in MC2R genes from selected fishes.

CONCLUSIONS/SIGNIFICANCE

The current repository of MC receptors illustrates that fishes have no MC3R ortholog. MC2R, MC5R, P2Y receptors and CHRM3 have novel intron insertions only in ray-finned fishes that underwent genome compaction. These receptors share one intron at an identical position suggestive of being inserted contemporaneously. In addition to repetitive elements, genome compaction is now believed to be a new hallmark that promotes intron insertions, as it requires rapid DNA breakage and subsequent repair processes to gain back normal functionality.

Adaptation of pineal expressed teleost exo-rod opsin to non-image forming photoreception through enhanced Meta II decay

Photoreception by vertebrates enables both image-forming vision and non-image-forming responses such as circadian photoentrainment. Over the recent years, distinct non-rod non-cone photopigments have been found to support circadian photoreception in diverse species. By allowing specialization to this sensory task a selective advantage is implied, but the nature of that specialization remains elusive. We have used the presence of distinct rod opsin genes specialized to either image-forming (retinal rod opsin) or non-image-forming (pineal exo-rod opsin) photoreception in ray-finned fish (Actinopterygii) to gain a unique insight into this problem. A comparison of biochemical features for these paralogous opsins in two model teleosts, Fugu pufferfish (Takifugu rubripes) and zebrafish (Danio rerio), reveals striking differences. While spectral sensitivity is largely unaltered by specialization to the pineal environment, in other aspects exo-rod opsins exhibit a behavior that is quite distinct from the cardinal features of the rod opsin family. While they display a similar thermal stability, they show a greater than tenfold reduction in the lifetime of the signaling active Meta II photoproduct. We show that these features reflect structural changes in retinal association domains of helices 3 and 5 but, interestingly, not at either of the two residues known to define these characteristics in cone opsins. Our findings suggest that the requirements of non-image-forming photoreception have lead exo-rod opsin to adopt a characteristic that seemingly favors efficient bleach recovery but not at the expense of absolute sensitivity.

Cytogenetic analysis in Polypterus ornatipinnis (Actinopterygii, Cladistia, Polypteridae) and 5S rDNA

Polypteridae is a family of archaic freshwater African fish that constitute an interesting subject for the study of the karyological evolution in vertebrates, on account of their primitive morphological characters and peculiar relationships with lower Osteichthyans. In this paper, a cytogenetic analysis on twenty specimens of both sexes of Polypterus ornatipinnis the ornate "bichir", coming from the Congo River basin, was performed by using both classical and molecular techniques. The karyotypic formula (2n=36; FN=72) was composed of 26 M+10 SM. The Alu I banding, performed to characterize heterochromatin in this species, was mainly centromeric. Both the chromosome location of the ribosomal 5S and 18S rRNA genes were examined by using Ag-NOR, classical C-banding, CMA(3) staining and FISH. CMA(3) marked all centromerical regions and showed the presence of two GC rich regions on the p arm of the chromosome pair n°1 and on the q arm of the pair n°14. Staining with Ag-NOR marked the only telomeric region of the chromosome n°1 p arm. After PCR, the 5S rDNA in this species was cloned, sequenced and analyzed. In the 665bp 5S rDNA sequence of P.ornatipinnis, a conserved 120bp gene region for the 5S rDNA was identified, followed by a non-transcribed variable spacer (NTS) which included simple repeats, microsatellites and a fragment of a non-LTR retrotransposon R-TEX. FISH with 5S rDNA marked the subtelomeric region of the q arm of the chromosome pair n°14, previously marked by CMA(3). FISH with 18S rDNA marked the telomeric region of the p arm of the pair n°1, previously marked both by Ag-NOR and CMA(3). The (GATA)(7) repeats marked the telomeric regions of all chromosome pairs, with the exclusion of the n°1, n°3 and n°14; hybridization with telomeric probes (TTAGGG)(n) showed signals at the end of all chromosomes. Karyotype evolution in Polypterus genus was finally discussed, including the new data obtained.

Distinguishing species of European sturgeons Acipenser spp. using microsatellite allele sequences

Five microsatellite markers were analysed and their alleles were sequenced for the three sturgeon species that lived in western Europe: the European sturgeon Acipenser sturio, the Atlantic sturgeon Acipenser oxyrinchus and the Adriatic sturgeon Acipenser naccarii. A total of 94 different allele sequences were obtained. Fixed mutations in the flanking regions or in the core repeat of microsatellites provided a clear distinction between the different species. Comparison of allele sequences also provided some insights into microsatellites and the evolution of Acipenser species. These nuclear markers can be used to solve species determination problems, and combined with mitochondrial markers, will be useful to identify introgression and hybridization among the three species. Moreover, because they are short and with a limited allele size range, they are particularly suited for analysis of museum specimens or archaeological remains.

The globin gene family of the cephalochordate amphioxus: implications for chordate globin evolution

BACKGROUND

The lancelet amphioxus (Cephalochordata) is a close relative of vertebrates and thus may enhance our understanding of vertebrate gene and genome evolution. In this context, the globins are one of the best studied models for gene family evolution. Previous biochemical studies have demonstrated the presence of an intracellular globin in notochord tissue and myotome of amphioxus, but the corresponding gene has not yet been identified. Genomic resources of Branchiostoma floridae now facilitate the identification, experimental confirmation and molecular evolutionary analysis of its globin gene repertoire.

RESULTS

We show that B. floridae harbors at least fifteen paralogous globin genes, all of which reveal evidence of gene expression. The protein sequences of twelve globins display the conserved characteristics of a functional globin fold. In phylogenetic analyses, the amphioxus globin BflGb4 forms a common clade with vertebrate neuroglobins, indicating the presence of this nerve globin in cephalochordates. Orthology is corroborated by conserved syntenic linkage of BflGb4 and flanking genes. The kinetics of ligand binding of recombinantly expressed BflGb4 reveals that this globin is hexacoordinated with a high oxygen association rate, thus strongly resembling vertebrate neuroglobin. In addition, possible amphioxus orthologs of the vertebrate globin X lineage and of the myoglobin/cytoglobin/hemoglobin lineage can be identified, including one gene as a candidate for being expressed in notochord tissue. Genomic analyses identify conserved synteny between amphioxus globin-containing regions and the vertebrate β-globin locus, possibly arguing against a late transpositional origin of the β-globin cluster in vertebrates. Some amphioxus globin gene structures exhibit minisatellite-like tandem duplications of intron-exon boundaries ("mirages"), which may serve to explain the creation of novel intron positions within the globin genes.

CONCLUSIONS

The identification of putative orthologs of vertebrate globin variants in the B. floridae genome underlines the importance of cephalochordates for elucidating vertebrate genome evolution. The present study facilitates detailed functional studies of the amphioxus globins in order to trace conserved properties and specific adaptations of respiratory proteins at the base of chordate evolution.

Size polymorphism in alleles of the myoglobin gene from biomphalaria mollusks

Introns are common among all eukaryotes, while only a limited number of introns are found in prokaryotes. Globin and globin-like proteins are widely distributed in nature, being found even in prokaryotes and a wide range of patterns of intron-exon have been reported in several eukaryotic globin genes. Globin genes in invertebrates show considerable variation in the positions of introns; globins can be found without introns, with only one intron or with three introns in different positions. In this work we analyzed the introns in the myoglobin gene from Biomphalaria glabrata, B. straminea and B. tenagophila. In the Biomphalaria genus, the myoglobin gene has three introns; these were amplified by PCR and analyzed by PCR-RFLP. Results showed that the size (number or nucleotides) and the nucleotide sequence of the coding gene of the myoglobin are variable in the three species. We observed the presence of size polymorphisms in intron 2 and 3; this characterizes a homozygous/heterozygous profile and it indicates the existence of two alleles which are different in size in each species of Biomphalaria. This polymorphism could be explored for specific identification of Biomphalaria individuals.

Intron-loss evolution of hatching enzyme genes in Teleostei

Background

Hatching enzyme, belonging to the astacin metallo-protease family, digests egg envelope at embryo hatching. Orthologous genes of the enzyme are found in all vertebrate genomes. Recently, we found that exon-intron structures of the genes were conserved among tetrapods, while the genes of teleosts frequently lost their introns. Occurrence of such intron losses in teleostean hatching enzyme genes is an uncommon evolutionary event, as most eukaryotic genes are generally known to be interrupted by introns and the intron insertion sites are conserved from species to species. Here, we report on extensive studies of the exon-intron structures of teleostean hatching enzyme genes for insight into how and why introns were lost during evolution.

Results

We investigated the evolutionary pathway of intron-losses in hatching enzyme genes of 27 species of Teleostei. Hatching enzyme genes of basal teleosts are of only one type, which conserves the 9-exon-8-intron structure of an assumed ancestor. On the other hand, otocephalans and euteleosts possess two types of hatching enzyme genes, suggesting a gene duplication event in the common ancestor of otocephalans and euteleosts. The duplicated genes were classified into two clades, clades I and II, based on phylogenetic analysis. In otocephalans and euteleosts, clade I genes developed a phylogeny-specific structure, such as an 8-exon-7-intron, 5-exon-4-intron, 4-exon-3-intron or intron-less structure. In contrast to the clade I genes, the structures of clade II genes were relatively stable in their configuration, and were similar to that of the ancestral genes. Expression analyses revealed that hatching enzyme genes were high-expression genes, when compared to that of housekeeping genes. When expression levels were compared between clade I and II genes, clade I genes tends to be expressed more highly than clade II genes.

Conclusions

Hatching enzyme genes evolved to lose their introns, and the intron-loss events occurred at the specific points of teleostean phylogeny. We propose that the high-expression hatching enzyme genes frequently lost their introns during the evolution of teleosts, while the low-expression genes maintained the exon-intron structure of the ancestral gene.

Polypteridae (Actinopterygii: Cladistia) and DANA-SINEs insertions

SINE sequences are interspersed throughout virtually all eukaryotic genomes and greatly outnumber the other repetitive elements. These sequences are of increasing interest for phylogenetic studies because of their diagnostic power for establishing common ancestry among taxa, once properly characterized. We identified and characterized a peculiar family of composite tRNA-derived short interspersed SINEs, DANA-SINEs, associated with mutational activities in Danio rerio, in a group of species belonging to one of the most basal bony fish families, the Polypteridae, in order to investigate their own inner specific phylogenetic relationships. DANA sequences were identified, sequenced and then localized, by means of fluorescent in situ hybridization (FISH), in six Polypteridae species (Polypterus delhezi, P. ornatipinnis, P. palmas, P. buettikoferi P. senegalus and Erpetoichthys calabaricus) After cloning, the sequences obtained were aligned for phylogenetic analysis, comparing them with three Dipnoan lungfish species (Protopterus annectens, P. aethiopicus, Lepidosiren paradoxa), and Lethenteron reissneri (Petromyzontidae)was used as outgroup. The obtained overlapping MP, ML and NJ tree clustered together the species belonging to the two taxonomically different Osteichthyans groups: the Polypteridae, by one side, and the Protopteridae by the other, with the monotypic genus Erpetoichthys more distantly related to the Polypterus genus comprising three distinct groups: P. palmas and P. buettikoferi, P. delhezi and P. ornatipinnis and P. senegalus. In situ hybridization with DANA probes marked along the whole chromosome arms in the metaphases of all the Polypteridae species examined.

Some novel intron positions in conserved Drosophila genes are caused by intron sliding or tandem duplication

BACKGROUND

Positions of spliceosomal introns are often conserved between remotely related genes. Introns that reside in non-conserved positions are either novel or remnants of frequent losses of introns in some evolutionary lineages. A recent gain of such introns is difficult to prove. However, introns verified as novel are needed to evaluate contemporary processes of intron gain.

RESULTS

We identified 25 unambiguous cases of novel intron positions in 31 Drosophila genes that exhibit near intron pairs (NIPs). Here, a NIP consists of an ancient and a novel intron position that are separated by less than 32 nt. Within a single gene, such closely-spaced introns are very unlikely to have coexisted. In most cases, therefore, the ancient intron position must have disappeared in favour of the novel one. A survey for NIPs among 12 Drosophila genomes identifies intron sliding (migration) as one of the more frequent causes of novel intron positions. Other novel introns seem to have been gained by regional tandem duplications of coding sequences containing a proto-splice site.

CONCLUSIONS

Recent intron gains sometimes appear to have arisen by duplication of exonic sequences and subsequent intronization of one of the copies. Intron migration and exon duplication together may account for a significant amount of novel intron positions in conserved coding sequences.

Characterization of the neuroligin gene family expression and evolution in zebrafish

Neuroligins constitute a family of transmembrane proteins localized at the postsynaptic side of both excitatory and inhibitory synapses of the central nervous system. They are involved in synaptic function and maturation and recent studies have linked mutations in specific human Neuroligins to mental retardation and autism. We isolated the human Neuroligin homologs in Danio rerio. Next, we studied their gene structures and we reconstructed the evolution of the Neuroligin genes across vertebrate phyla. Using reverse-transcriptase polymerase chain reaction, we analyzed the expression and alternative splicing pattern of each gene during zebrafish embryonic development and in different adult organs. By in situ hybridization, we analyzed the temporal and spatial expression pattern during embryonic development and larval stages and we found that zebrafish Neuroligins are expressed throughout the nervous system. Globally, our results indicate that, during evolution, specific subfunctionalization events occurred within paralogous members of this gene family in zebrafish.

Nonsense-mediated decay enables intron gain in Drosophila

Intron number varies considerably among genomes, but despite their fundamental importance, the mutational mechanisms and evolutionary processes underlying the expansion of intron number remain unknown. Here we show that Drosophila, in contrast to most eukaryotic lineages, is still undergoing a dramatic rate of intron gain. These novel introns carry significantly weaker splice sites that may impede their identification by the spliceosome. Novel introns are more likely to encode a premature termination codon (PTC), indicating that nonsense-mediated decay (NMD) functions as a backup for weak splicing of new introns. Our data suggest that new introns originate when genomic insertions with weak splice sites are hidden from selection by NMD. This mechanism reduces the sequence requirement imposed on novel introns and implies that the capacity of the spliceosome to recognize weak splice sites was a prerequisite for intron gain during eukaryotic evolution.

The surprising observation 30 years ago that genes are interrupted by non-coding introns changed our view of gene architecture. Intron number varies dramatically among species; ranging from nine introns/gene in humans to less than one in some simple eukyarotes. Here we ask where new introns come from and how they are maintained in a population. We find that novel introns do not arise from pre-existing introns, although the mechanisms that generate novel introns remain unclear. We also show that novel introns carry only weak signals for their identification and removal, and therefore depend on nonsense-mediated decay (NMD). NMD maintains RNA quality control by degrading transcripts that have not been spliced properly. We propose that NMD shelters novel introns from natural selection. This increases the likelihood that a novel intron will rise in frequency and be maintained within a population, thus increasing the rate of intron gain.

Base composition, selection, and phylogenetic significance of indels in the recombination activating gene-1 in vertebrates

BACKGROUND

The Recombination Activating Proteins, RAG1 and RAG2, play a crucial role in the immune response in vertebrates. Among the nuclear markers currently used for phylogenetic purposes, Rag1 has especially enjoyed enormous popularity, since it successfully contributed to elucidating the relationships among and within a large variety of vertebrate lineages. We here report on a comparative investigation of the genetic variation, base composition, presence of indels, and selection in Rag1 in different vertebrate lineages (Actinopterygii, Amphibia, Aves, Chondrichthyes, Crocodylia, Lepidosauria, Mammalia, and Testudines) through the analysis of 582 sequences obtained from Genbank. We also analyze possible differences between distinct parts of the gene with different type of protein functions.

RESULTS

In the vertebrate lineages studied, Rag1 is over 3 kb long. We observed a high level of heterogeneity in base composition at the 3(rd )codon position in some of the studied vertebrate lineages and in some specific taxa. This result is also paralleled by taxonomic differences in the GC content at the same codon position. Moreover, positive selection occurs at some sites in Aves, Lepidosauria and Testudines. Indels, which are often used as phylogenetic characters, are more informative across vertebrates in the 5' than in the 3'-end of the gene. When the entire gene is considered, the use of indels as phylogenetic character only recovers one major vertebrate clade, the Actinopterygii. However, in numerous cases insertions or deletions are specific to a monophyletic group.

CONCLUSIONS

Rag1 is a phylogenetic marker of undoubted quality. Our study points to the need of carrying out a preliminary investigation on the base composition and the possible existence of sites under selection of this gene within the groups studied to avoid misleading resolution. The gene shows highly heterogeneous base composition, which affects some taxa in particular and contains sites under positive selection in some vertebrate lineages in the 5'-end. The first part of the gene (5'-end) is more variable than the second (3'-end), and less affected by a heterogeneous base composition. However, in some vertebrate lineages the 5'-end of the gene is not yet widely used for phylogenetic studies.

Tandem repeats modify the structure of human genes hosted in segmental duplications

BACKGROUND

Recently duplicated genes are often subject to genomic rearrangements that can lead to the development of novel gene structures. Here we specifically investigated the effect of variations in internal tandem repeats (ITRs) on the gene structure of human paralogs located in segmental duplications.

RESULTS

We found that around 7% of the primate-specific genes located within duplicated regions of the genome contain variable tandem repeats. These genes are members of large groups of recently duplicated paralogs that are often polymorphic in the human population. Half of the identified ITRs occur within coding exons and may be either kept or spliced out from the mature transcript. When ITRs reside within exons, they encode variable amino acid repeats. When located at exon-intron boundaries, ITRs can generate alternative splicing patterns through the formation of novel introns.

CONCLUSIONS

Our study shows that variation in the number of ITRs impacts on recently duplicated genes by modifying their coding sequence, splicing pattern, and tissue expression. The resulting effect is the production of a variety of primate-specific proteins, which mostly differ in number and sequence of amino acid repeats.

Whole-genome phylogeny of mammals: evolutionary information in genic and nongenic regions

Ten complete mammalian genome sequences were compared by using the "feature frequency profile" (FFP) method of alignment-free comparison. This comparison technique reveals that the whole nongenic portion of mammalian genomes contains evolutionary information that is similar to their genic counterparts--the intron and exon regions. We partitioned the complete genomes of mammals (such as human, chimp, horse, and mouse) into their constituent nongenic, intronic, and exonic components. Phylogenic species trees were constructed for each individual component class of genome sequence data as well as the whole genomes by using standard tree-building algorithms with FFP distances. The phylogenies of the whole genomes and each of the component classes (exonic, intronic, and nongenic regions) have similar topologies, within the optimal feature length range, and all agree well with the evolutionary phylogeny based on a recent large dataset, multispecies, and multigene-based alignment. In the strictest sense, the FFP-based trees are genome phylogenies, not species phylogenies. However, the species phylogeny is highly related to the whole-genome phylogeny. Furthermore, our results reveal that the footprints of evolutionary history are spread throughout the entire length of the whole genome of an organism and are not limited to genes, introns, or short, highly conserved, nongenic sequences that can be adversely affected by factors (such as a choice of sequences, homoplasy, and different mutation rates) resulting in inconsistent species phylogenies.

Multiple gains of spliceosomal introns in a superfamily of vertebrate protease inhibitor genes

Background

Intron gains reportedly are very rare during evolution of vertebrates, and the mechanisms underlying their creation are largely unknown. Previous investigations have shown that, during metazoan radiation, the exon-intron patterns of serpin superfamily genes were subject to massive changes, in contrast to many other genes.

Results

Here we investigated intron dynamics in the serpin superfamily in lineages pre- and postdating the split of vertebrates. Multiple intron gains were detected in a group of ray-finned fishes, once the canonical groups of vertebrate serpins had been established. In two genes, co-occurrence of non-standard introns was observed, implying that intron gains in vertebrates may even happen concomitantly or in a rapidly consecutive manner. DNA breakage/repair processes associated with genome compaction are introduced as a novel factor potentially favoring intron gain, since all non-canonical introns were found in a lineage of ray-finned fishes that experienced genomic downsizing.

Conclusion

Multiple intron acquisitions were identified in serpin genes of a lineage of ray-finned fishes, but not in any other vertebrates, suggesting that insertion rates for introns may be episodically increased. The co-occurrence of non-standard introns within the same gene discloses the possibility that introns may be gained simultaneously. The sequences flanking the intron insertion points correspond to the proto-splice site consensus sequence MAG↑N, previously proposed to serve as intron insertion site. The association of intron gains in the serpin superfamily with a group of fishes that underwent genome compaction may indicate that DNA breakage/repair processes might foster intron birth.

Beyond linear sequence comparisons: the use of genome-level characters for phylogenetic reconstruction

The first whole genomes to be compared for phylogenetic inference were those of mitochondria, which provided the first sets of genome-level characters for phylogenetic reconstruction. Most powerful among these characters has been the comparisons of the relative arrangements of genes, which has convincingly resolved numerous branch points, including those that had remained recalcitrant even to very large molecular sequence comparisons. Now the world faces a tsunami of complete nuclear genome sequences. In addition to the tremendous amount of DNA sequence that is becoming available for comparison, there is also a potential for many more genome-level characters to be developed, including the relative positions of introns, the domain structures of proteins, gene family membership, the presence of particular biochemical pathways, aspects of DNA replication or transcription, and many others. These characters can be especially convincing owing to their low likelihood of reverting to a primitive condition or occurring independently in separate lineages, thereby reducing the occurrence of homoplasy. The comparisons of organelle genomes pioneered the way for using such features for phylogenetic reconstructions, and it is almost certainly true, as ever more genomic sequence becomes available, that further use of genome-level characters will play a big role in outlining the relationships among major animal groups.

Alternative splicing: a missing piece in the puzzle of intron gain

Spliceosomal introns, a hallmark of eukaryotic gene organization, were an unexpected discovery. After three decades, crucial issues such as when and how introns first appeared in evolution remain unsettled. An issue yet to be answered is how intron positions arise de novo. Phylogenetic investigations concur that intron positions continue to emerge, at least in some lineages. Yet genomic scans for the sources of introns occupying new positions have been fruitless. Two alternative solutions to this paradox are: (i) formation of new intron positions halted before the recent past and (ii) it continues to occur, but through processes different from those generally assumed. One process generally dismissed is intron sliding--the relocation of a preexisting intron over short distances--because of supposed associated deleterious effects. The puzzle of intron gain arises owing to a pervasive operational definition of introns, which sees them as precisely demarcated segments of the genome separated from the neighboring nonintronic DNA by unmovable limits. Intron homology is defined as position homology. Recent studies of pre-mRNA processing indicate that this assumption needs to be revised. We incorporate recent advances on the evolutionarily frequent process of alternative splicing, by which exons of primary transcripts are spliced in different patterns, into a new model of intron sliding that accounts for the diversity of intron positions. We posit that intron positional diversity is driven by two overlapping processes: (i) background process of continuous relocation of preexisting introns by sliding and (ii) spurts of extensive gain/loss of new intron sequences.