Intron length evolution in Drosophila

I present data on the evolution of intron lengths among 3 closely related Drosophila species, D. melanogaster, Drosophila simulans, and Drosophila yakuba. Using D. yakuba as an outgroup, I mapped insertion and deletion mutations in 148 introns (spanning approximately 30 kb) to the D. melanogaster and D. simulans lineages. Intron length evolution in the 2 sister species has been different: in D. melanogaster, X-linked introns have increased slightly in size, whereas autosomal ones have decreased slightly in size; in D. simulans, both X-linked and autosomal introns have decreased in size. To understand the possible evolutionary causes of these lineage- and chromosome-specific patterns of intron evolution, I studied insertion-deletion (indel) polymorphism and divergence in D. melanogaster. Small insertion mutations segregate at elevated frequencies and enjoy elevated probabilities of fixation, particularly on the X chromosome. In contrast, there is no detectable X chromosome effect on fixations in D. simulans. These findings suggest X chromosome-specific selection or biased gene conversion-gap repair favoring insertions in D. melanogaster but not in D. simulans. These chromosome- and lineage-specific patterns of indel substitution are not easily explained by existing general population genetic models of intron length evolution. Genomic data from D. melanogaster further suggest that the forces described here affect introns and intergenic regions similarly.     

The novel MER transposon-derived miRNAs in human genome

MicroRNAs (miRNAs) are small RNA molecules (~ 20–30 nucleotides) that generally act in gene silencing and translational repression through the RNA interference pathway. They generally originate from intergenic genomic regions, but some are found in genomic regions that have been characterized such as introns, exons, and transposable elements (TE). To identify the miRNAs that are derived from palindromic MERs, we analyzed MER paralogs in human genome. The structures of the palindromic MERs were similar to the hairpin structure of miRNA in humans. Three miRNAs derived from MER96 located on chromosome 3, and MER91C paralogs located on chromosome 8 and chromosome 17 were identified in HeLa, HCT116, and HEK293 cell lines. The interactions between these MER-derived miRNAs and AGO1, AGO2, and AGO3 proteins were validated by immunoprecipitation assays. The data suggest that miRNAs derived from transposable elements could widely affect various target genes in the human genome.     

Prevalence of intron gain over intron loss in the evolution of paralogous gene families

The mechanisms and evolutionary dynamics of intron insertion and loss in eukaryotic genes remain poorly understood. Reconstruction of parsimonious scenarios of gene structure evolution in paralogous gene families in animals and plants revealed numerous gains and losses of introns. In all analyzed lineages, the number of acquired new introns was substantially greater than the number of lost ancestral introns. This trend held even for lineages in which vertical evolution of genes involved more intron losses than gains, suggesting that gene duplication boosts intron insertion. However, dating gene duplications and the associated intron gains and losses based on the molecular clock assumption showed that very few, if any, introns were gained during the last ~100 million years of animal and plant evolution, in agreement with previous conclusions reached through analysis of orthologous gene sets. These results are generally compatible with the emerging notion of intensive insertion and loss of introns during transitional epochs in contrast to the relative quiet of the intervening evolutionary spans.     

Colonization of heterochromatic genes by transposable elements in Drosophila

As a further step toward understanding transposable element-host genome interactions, we investigated the molecular anatomy of introns from five heterochromatic and 22 euchromatic protein-coding genes of Drosophila melanogaster. A total of 79 kb of intronic sequences from heterochromatic genes and 355 kb of intronic sequences from euchromatic genes have been used in Blast searches against Drosophila transposable elements (TEs). The results show that TE-homologous sequences belonging to 19 different families represent about 50% of intronic DNA from heterochromatic genes. In contrast, only 0.1% of the euchromatic intron DNA exhibits homology to known TEs. Intraspecific and interspecific size polymorphisms of introns were found, which are likely to be associated with changes in TE-related sequences. Together, the enrichment in TEs and the apparent dynamic state of heterochromatic introns suggest that TEs contribute significantly to the evolution of genes located in heterochromatin.     

Analysis of transposon interruptions suggests selection for L1 elements on the X chromosome

It has been hypothesised that the massive accumulation of L1 transposable elements on the X chromosome is due to their function in X inactivation, and that the accumulation of Alu elements near genes is adaptive. We tested the possible selective advantage of these two transposable element (TE) families with a novel method, interruption analysis. In mammalian genomes, a large number of TEs interrupt other TEs due to the high overall abundance and age of repeats, and these interruptions can be used to test whether TEs are selectively neutral. Interruptions of TEs, which are beneficial for the host, are expected to be deleterious and underrepresented compared with neutral ones. We found that L1 elements in the regions of the X chromosome that contain the majority of the inactivated genes are significantly less frequently interrupted than on the autosomes, while L1s near genes that escape inactivation are interrupted with higher frequency, supporting the hypothesis that L1s on the X chromosome play a role in its inactivation. In addition, we show that TEs are less frequently interrupted in introns than in intergenic regions, probably due to selection against the expansion of introns, but the insertion pattern of Alus is comparable to other repeats.     

The Y chromosome as a target for acquired and amplified genetic material in evolution

The special properties of the Y chromosome stem form the fact that it is a non-recombining degenerate derivative of the X chromosome. The absence of homologous recombination between the X and the Y chromosome leads to gradual degeneration of various Y chromosome genes on an evolutionary timescale. The absence of recombination, however, also favors the accumulation of transposable elements on the Y chromosome during its evolution, as seen with both Drosophila and mammalian Y chromosomes. Alongside these processes, the acquisition and amplification of autosomal male benefit genes occur. This review will focus on recent studies that reveal the autosome-acquired genes on the Y chromosome of both Drosophila and humans. The evolution of the acquired and amplified genes on the Y chromosome is also discussed. Molecular and comparative analyses of Y-linked repeats in the Drosophila melanogaster genome demonstrate that there was a period of their degeneration followed by a period of their integration into RNAi silencing, which was beneficial for male fertility. Finally, the function of non-coding RNA produced by amplified Y chromosome genetic elements will be discussed.     

A glimpse of a putative pre-intron phase of eukaryotic evolution

Comparison of the exon-intron structures of ancient eukaryotic paralogs reveals the absence of conserved intron positions in these genes. This is in contrast to the conservation of intron positions in orthologous genes from even the most evolutionarily distant eukaryotes and in more recent paralogs. The lack of conserved intron positions in ancient paralogs probably reflects the origination of these genes during the earliest phase of eukaryotic evolution, which was characterized by concomitant invasion of genes by group II self-splicing elements (which were to become introns in the future) and extensive duplication of genes.     

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

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

Molecular analysis of an alcohol dehydrogenase (Adh) gene from chromosome 1 of wheat

We have cloned and determined the nucleotide sequence of a gene encoding alcohol dehydrogenase (Adh) from Triticum aestivum cv. Millewa. Southern analysis using cv. Chinese Spring nullisomic-tetrasomic and ditelosomic lines established that the cloned gene mapped to the long arm of chromosome 1A and does not correspond to any previously identified wheat Adh locus. Southern analysis also provided evidence for triplicate copies of this Adh gene on the homoeologous group 1 chromosomes, while Northern blots indicated that the homoeologous group 1 Adh genes, like several other plant Adh genes, are transcribed under anaerobic conditions. Sequence analysis indicates that the cloned gene has a structure similar to both monocot and dicot Adh genes with an open reading frame encoding a polypeptide of 379 amino acids. Sequences important for eucaryotic gene expression such as the TATA box, polyadenylation signal, and intron splice sites were found in the expected positions. The open reading frame is interrupted by 8 introns which are in identical positions with 8 of the 9 introns in maize and pea Adh genes, suggesting that during evolution there are processes occurring that result in the loss of introns. Sequence analysis also revealed that the cloned wheat Adh gene shared extensive homology with the barley Adh3 gene not only in the coding region but also in the noncoding regions. However, this homology is discontinuous as a result of a 1.8-kbp insertion (TLM), which is present in the cloned wheat Adh gene and absent in the barley Adh3 gene. Sequence analysis of this insertion reveals features characteristic of the short terminal inverted repeat class of eucaryotic transposable elements. We have no evidence for the transposition of the TLM element. However, Southern blots reveal multiple copies of sequences related to TLM in the wheat genome and in other closely related species, suggesting that transposition may once have played an important role in the evolution of the Gramineae family.     

Recombination drives vertebrate genome contraction

Selective and/or neutral processes may govern variation in DNA content and, ultimately, genome size. The observation in several organisms of a negative correlation between recombination rate and intron size could be compatible with a neutral model in which recombination is mutagenic for length changes. We used whole-genome data on small insertions and deletions within transposable elements from chicken and zebra finch to demonstrate clear links between recombination rate and a number of attributes of reduced DNA content. Recombination rate was negatively correlated with the length of introns, transposable elements, and intergenic spacer and with the rate of short insertions. Importantly, it was positively correlated with gene density, the rate of short deletions, the deletion bias, and the net change in sequence length. All these observations point at a pattern of more condensed genome structure in regions of high recombination. Based on the observed rates of small insertions and deletions and assuming that these rates are representative for the whole genome, we estimate that the genome of the most recent common ancestor of birds and lizards has lost nearly 20% of its DNA content up until the present. Expansion of transposable elements can counteract the effect of deletions in an equilibrium mutation model; however, since the activity of transposable elements has been low in the avian lineage, the deletion bias is likely to have had a significant effect on genome size evolution in dinosaurs and birds, contributing to the maintenance of a small genome. We also demonstrate that most of the observed correlations between recombination rate and genome contraction parameters are seen in the human genome, including for segregating indel polymorphisms. Our data are compatible with a neutral model in which recombination drives vertebrate genome size evolution and gives no direct support for a role of natural selection in this process.     

Investigation of loss and gain of introns in the compact genomes of pufferfishes (Fugu and Tetraodon)

We have investigated intron evolution in the compact genomes of 2 closely related species of pufferfishes, Fugu rubripes and Tetraodon nigroviridis, that diverged about 32 million years ago (MYA). Analysis of 148,028 aligned intron positions in 13,547 gene pairs using human as an outgroup identified 57 and 24 intron losses in Tetraodon and fugu lineages, respectively, and no gain in either lineage. For comparison, we analyzed 144,545 intron positions in 12,866 orthologous pairs of genes in human and mouse that diverged about 61 MYA using dog as an outgroup and identified 51 intron losses in mouse and 3 losses in human and no gain. The rate of intron loss in Tetraodon is higher than that in fugu, mouse, and human but lower than the previous estimates for other eukaryotes. The introns lost in pufferfishes and mammals are significantly shorter than the mean size of introns in the genome. One intron deleted in fugu and another in Tetraodon have left behind 6 and 3 nucleotides, respectively, suggesting that they were lost due to genomic deletions. Such losses of introns are likely to be the result of a higher rate of DNA deletions experienced by the genomes of pufferfishes compared with mammals. The shorter generation time of Tetraodon compared with fugu, and the rich diversity and higher activity of transposable elements in pufferfishes compared with mammals, may be responsible for the higher rate of intron loss in Tetraodon. Our findings indicate that overall very little intron turnover has occurred in pufferfishes and mammals during recent evolution and that intron gain is an extremely rare event in vertebrate evolution.     

Survey of group I and group II introns in 29 sequenced genomes of the Bacillus cereus group: insights into their spread and evolution

Group I and group II introns are different catalytic self-splicing and mobile RNA elements that contribute to genome dynamics. In this study, we have analyzed their distribution and evolution in 29 sequenced genomes from the Bacillus cereus group of bacteria. Introns were of different structural classes and evolutionary origins, and a large number of nearly identical elements are shared between multiple strains of different sources, suggesting recent lateral transfers and/or that introns are under a strong selection pressure. Altogether, 73 group I introns were identified, inserted in essential genes from the chromosome or newly described prophages, including the first elements found within phages in bacterial plasmids. Notably, bacteriophages are an important source for spreading group I introns between strains. Furthermore, 77 group II introns were found within a diverse set of chromosomal and plasmidic genes. Unusual findings include elements located within conserved DNA metabolism and repair genes and one intron inserted within a novel retroelement. Group II introns are mainly disseminated via plasmids and can subsequently invade the host genome, in particular by coupling mobility with host cell replication. This study reveals a very high diversity and variability of mobile introns in B. cereus group strains.     

Evolution and Survival on Eutherian Sex Chromosomes

Since the two eutherian sex chromosomes diverged from an ancestral autosomal pair, the X has remained relatively gene-rich, while the Y has lost most of its genes through the accumulation of deleterious mutations in nonrecombining regions. Presently, it is unclear what is distinctive about genes that remain on the Y chromosome, when the sex chromosomes acquired their unique evolutionary rates, and whether X-Y gene divergence paralleled that of paralogs located on autosomes. To tackle these questions, here we juxtaposed the evolution of X and Y homologous genes (gametologs) in eutherian mammals with their autosomal orthologs in marsupial and monotreme mammals. We discovered that genes on the X and Y acquired distinct evolutionary rates immediately following the suppression of recombination between the two sex chromosomes. The Y-linked genes evolved at higher rates, while the X-linked genes maintained the lower evolutionary rates of the ancestral autosomal genes. These distinct rates have been maintained throughout the evolution of X and Y. Specifically, in humans, most X gametologs and, curiously, also most Y gametologs evolved under stronger purifying selection than similarly aged autosomal paralogs. Finally, after evaluating the current experimental data from the literature, we concluded that unique mRNA/protein expression patterns and functions acquired by Y (versus X) gametologs likely contributed to their retention. Our results also suggest that either the boundary between sex chromosome strata 3 and 4 should be shifted or that stratum 3 should be divided into two strata.     

Intron size, abundance, and distribution within untranslated regions of genes

Most research concerning the evolution of introns has largely considered introns within coding sequences (CDSs), without regard for introns located within untranslated regions (UTRs) of genes. Here, we directly determined intron size, abundance, and distribution in UTRs of genes using full-length cDNA libraries and complete genome sequences for four species, Arabidopsis thaliana, Drosophila melanogaster, human, and mouse. Overall intron occupancy (introns/exon kbp) is lower in 5' UTRs than CDSs, but intron density (intron occupancy in regions containing introns) tends to be higher in 5' UTRs than in CDSs. Introns in 5' UTRs are roughly twice as large as introns in CDSs, and there is a sharp drop in intron size at the 5' UTR-CDS boundary. We propose a mechanistic explanation for the existence of selection for larger intron size in 5' UTRs, and outline several implications of this hypothesis. We found introns to be randomly distributed within 5' UTRs, so long as a minimum required exon size was assumed. Introns in 3' UTRs were much less abundant than in 5' UTRs. Though this was expected for human and mouse that have intron-dependent nonsense-mediated decay (NMD) pathways that discourage the presence of introns within the 3' UTR, it was also true for A. thaliana and D. melanogaster, which may lack intron-dependent NMD. Our findings have several implications for theories of intron evolution and genome evolution in general.     

Generation of a database containing discordant intron positions in eukaryotic genes (MIDB)

MOTIVATION: Intron sliding is the relocation of intron-exon boundaries over short distances and is often also referred to as intron slippage or intron migration or intron drift. We have generated a database containing discordant intron positions in homologous genes (MIDB--Mismatched Intron DataBase). Discordant intron positions are those that are either closely located in homologous genes (within a window of 10 nucleotides) or an intron position that is present in one gene but not in any of its homologs. The MIDB database aims at systematically collecting information about mismatched introns in the genes from GenBank and organizing it into a form useful for understanding the genomics and dynamics of introns thereby helping understand the evolution of genes. RESULTS: Intron displacement or sliding is critically important for explaining the present distribution of introns among orthologous and paralogous genes. MIDB allows examining of intron movements and allows mapping of intron positions from homologous proteins onto a single sequence. The database is of potential use for molecular biologists in general and for researchers who are interested in gene evolution and eukaryotic gene structure. Partial analysis of this database allowed us to identify a few putative cases of intron sliding. AVAILABILITY: http://intron.bic.nus.edu.sg/midb/midb.html