Molecular Systematics and Evolution of the Growth Hormone Introns in the Salmoninae

DNA sequence data was collected for the C and D introns in the duplicate growth hormone loci (GH1 and GH2) from Brachmystax lenok, two subspecies of Hucho hucho, Hucho (Parahucho) perryi, Salmo salar, Salmo trutta, Acantholingua ohridana (Salmothymus), six species of Salvelinus, eight species of Oncorhynchus including O. masou, and three outgroups including Thymallus thymallus, Coregonus artedi, and Coregonus clupeaformis. Phylogenetic analyses were performed using maximum parsimony and maximum likelihood (PAUP, version 4.08beta) with gaps as missing data and as a fifth base. B. lenok was basal in all of the trees and all of the other genera were monophyletic with the exception that A. ohridana always placed within Salmo, and H. hucho sp. often placed with B. lenok. The GH1 introns supported a sister relationship between Oncorhynchus and Salvelinus, while the combined GH2 introns were ambiguous at this node. This result contrasts with trees based on morphology and the ribosomal ITS1 sequences that support a sister relationship between Salmo and Oncorhynchus. The only estrogen response element (ERE) in the gene is found in the C intron and has mutated in GH2 in all of the species except B. lenok. The ERE element in GH1 has undergone another mutation in all of the species except for B. lenok, and members of the two genera Salvelinus and Oncorhynchus. Thus these latter two genera are the only ones with a difference in expression of GH1 and GH2 in the presence of estrogen. Differences in selective pressure on the introns in the duplicate genes in different taxa could account for the conflicting results obtained in the phylogenetic analysis.     

Interspecific Relationships among Charrs Based on Phylogenetic Analysis of Nuclear Growth Hormone Intron Sequences

DNA samples were obtained from six charr species: Salvelinus alpinus, S. malma, S. confluentus, S. leucomaenis, S. namaycush and S. fontinalis and two outgroups: Salmo salar and S. trutta. The C and D introns from the duplicate growth hormone loci (GH1 and GH2) were amplified using the polymerase chain reaction and sequenced from each of the species. Phylogenetic analyses with and without gaps were performed using the beta version of PAUP. The GH introns are evolving relatively slowly, so the data were combined with sequences from the ITS1 and ITS2 for a phylogenetic analysis of the charr species. The GH intron data and the combined data supported two groups of sister species among charrs: S. alpinus and S. malma, and S. confluentus and S. leucomaenis.     

Parallel Loss of Plastid Introns and Their Maturase in the Genus Cuscuta

Plastid genome content and arrangement are highly conserved across most land plants and their closest relatives, streptophyte algae, with nearly all plastid introns having invaded the genome in their common ancestor at least 450 million years ago. One such intron, within the transfer RNA trnK-UUU, contains a large open reading frame that encodes a presumed intron maturase, matK. This gene is missing from the plastid genomes of two species in the parasitic plant genus Cuscuta but is found in all other published land plant and streptophyte algal plastid genomes, including that of the nonphotosynthetic angiosperm Epifagus virginiana and two other species of Cuscuta. By examining matK and plastid intron distribution in Cuscuta, we add support to the hypothesis that its normal role is in splicing seven of the eight group IIA introns in the genome. We also analyze matK nucleotide sequences from Cuscuta species and relatives that retain matK to test whether changes in selective pressure in the maturase are associated with intron deletion. Stepwise loss of most group IIA introns from the plastid genome results in substantial change in selective pressure within the hypothetical RNA-binding domain of matK in both Cuscuta and Epifagus, either through evolution from a generalist to a specialist intron splicer or due to loss of a particular intron responsible for most of the constraint on the binding region. The possibility of intron-specific specialization in the X-domain is implicated by evidence of positive selection on the lineage leading to C. nitida in association with the loss of six of seven introns putatively spliced by matK. Moreover, transfer RNA gene deletion facilitated by parasitism combined with an unusually high rate of intron loss from remaining functional plastid genes created a unique circumstance on the lineage leading to Cuscuta subgenus Grammica that allowed elimination of matK in the most species-rich lineage of Cuscuta.     

EST analysis of Ostreococcus lucimarinus, the most compact eukaryotic genome, shows an excess of introns in highly expressed genes

Background

The genome of the pico-eukaryotic (bacterial-sized) prasinophyte green alga Ostreococcus lucimarinus has one of the highest gene densities known in eukaryotes, yet it contains many introns. Phylogenetic studies suggest this unusually compact genome (13.2 Mb) is an evolutionarily derived state among prasinophytes. The presence of introns in the highly reduced O. lucimarinus genome appears to be in opposition to simple explanations of genome evolution based on unidirectional tendencies, either neutral or selective. Therefore, patterns of intron retention in this species can potentially provide insights into the forces governing intron evolution.

Methodology/Principal Findings

Here we studied intron features and levels of expression in O. lucimarinus using expressed sequence tags (ESTs) to annotate the current genome assembly. ESTs were assembled into unigene clusters that were mapped back to the O. lucimarinus Build 2.0 assembly using BLAST and the level of gene expression was inferred from the number of ESTs in each cluster. We find a positive correlation between expression levels and both intron number (R = +0.0893, p = <0.0005) and intron density (number of introns/kb of CDS; R = +0.0753, p = <0.005).

Conclusions/Significance

In a species with a genome that has been recently subjected to a great reduction of non-coding DNA, these results imply the existence of selective/functional roles for introns that are principally detectable in highly expressed genes. In these cases, introns are likely maintained by balancing the selective forces favoring their maintenance with other mutational and/or selective forces acting on genome size.     

GROUP I INTRON VERSUS ITS SEQUENCES IN PHYLOGENY OF CETRARIOID LICHENS

Abstract: Phylogenetic trees based on group I intron sequences and on internal transcribed spacer (ITS) sequences of mycobiont ribosomal genes were calculated and compared. Eight cetrarioid and four non-cetrarioid species of the Parmeliaceae were compared. The phylogeny based on group I intron sequences is partly congruent with the ITS sequence phylogeny. Group I intron sequences are presumably less informative for infragenic studies. The introns have a length of 214–233 nucleotides, and differ at up to 33% of the bases between species. All introns analysed are located between the positions 1516 and 1517 of the fungal 18S ribosomal RNA gene. Cetrarioid lichens form a non-homogeneous group within theParmeliaceae according to both group I intron and ITS sequences.     

Frequent Gain and Loss of Introns in Fungal Cytochrome b Genes

In this study, all available cytochrome b (Cyt b) genes from the GOBASE database were compiled and the evolutionary dynamics of the Cyt b gene introns was assessed. Cyt b gene introns were frequently present in the fungal kingdom and some lower plants, but generally absent or rare in Chromista, Protozoa, and Animalia. Fungal Cyt b introns were found at 35 positions in Cyt b genes and the number of introns varied at individual positions from a single representative to 32 different introns at position 131, showing a wide and patchy distribution. Many homologous introns were present at the same position in distantly related species but absent in closely related species, suggesting that introns of the Cyt b genes were frequently lost. On the other hand, highly similar intron sequences were observed in some distantly related species rather than in closely related species, suggesting that these introns were gained independently, likely through lateral transfers. The intron loss-and-gain events could be mediated by transpositions that might have occurred between nuclear and mitochondria. Southern hybridization analysis confirmed that some introns contained repetitive sequences and might be transposable elements. An intron gain in Botryotinia fuckeliana prevented the development of QoI fungicide resistance, suggesting that intron loss-and-gain events were not necessarily beneficial to their host organisms.     

Analysis of Fungal Genomes Reveals Commonalities of Intron Gain or Loss and Functions in Intron-Poor Species

Previous evolutionary reconstructions have concluded that early eukaryotic ancestors including both the last common ancestor of eukaryotes and of all fungi had intron-rich genomes. By contrast, some extant eukaryotes have few introns, underscoring the complex histories of intron–exon structures, and raising the question as to why these few introns are retained. Here, we have used recently available fungal genomes to address a variety of questions related to intron evolution. Evolutionary reconstruction of intron presence and absence using 263 diverse fungal species supports the idea that massive intron reduction through intron loss has occurred in multiple clades. The intron densities estimated in various fungal ancestors differ from zero to 7.6 introns per 1 kb of protein-coding sequence. Massive intron loss has occurred not only in microsporidian parasites and saccharomycetous yeasts, but also in diverse smuts and allies. To investigate the roles of the remaining introns in highly-reduced species, we have searched for their special characteristics in eight intron-poor fungi. Notably, the introns of ribosome-associated genes RPL7 and NOG2 have conserved positions; both intron-containing genes encoding snoRNAs. Furthermore, both the proteins and snoRNAs are involved in ribosome biogenesis, suggesting that the expression of the protein-coding genes and noncoding snoRNAs may be functionally coordinated. Indeed, these introns are also conserved in three-quarters of fungi species. Our study shows that fungal introns have a complex evolutionary history and underappreciated roles in gene expression.     

Novel Introner-Like Elements in fungi Are Involved in Parallel Gains of Spliceosomal Introns

Spliceosomal introns are key components of the eukaryotic gene structure. Although they contributed to the emergence of eukaryotes, their origin remains elusive. In fungi, they might originate from the multiplication of invasive introns named Introner-Like Elements (ILEs). However, so far ILEs have been observed in six fungal species only, including Fulvia fulva and Dothistroma septosporum (Dothideomycetes), arguing against ILE insertion as a general mechanism for intron gain. Here, we identified novel ILEs in eight additional fungal species that are phylogenetically related to F. fulva and D. septosporum using PCR amplification with primers derived from previously identified ILEs. The ILE content appeared unique to each species, suggesting independent multiplication events. Interestingly, we identified four genes each containing two gained ILEs. By analysing intron positions in orthologues of these four genes in Ascomycota, we found that three ILEs had inserted within a 15 bp window that contains regular spliceosomal introns in other fungal species. These three positions are not the result of intron sliding because ILEs are newly gained introns. Furthermore, the alternative hypothesis of an inferred ancestral gain followed by independent losses contradicts the observed degeneration of ILEs. These observations clearly indicate three parallel intron gains in four genes that were randomly identified. Our findings suggest that parallel intron gain is a phenomenon that has been highly underestimated in ILE-containing fungi, and likely in the whole fungal kingdom.     

Sense-antisense gene overlap is probably a cause for retaining the few introns in <Emphasis Type="Italic">Giardia</Emphasis> genome and the implications

Background

It is widely accepted that the last eukaryotic common ancestor and early eukaryotes were intron-rich and intron loss dominated subsequent evolution, thus the presence of only very few introns in some modern eukaryotes must be the consequence of massive loss. But it is striking that few eukaryotes were found to have completely lost introns. Despite extensive research, the causes of massive intron losses remain elusive. Actually the reverse question -- how the few introns can be retained under the evolutionary selection pressure of intron loss -- is equally significant but was rarely studied, except that it was conjectured that the essential functions of some introns prevent their loss. The situation that extremely few (eight) spliceosome-mediated cis-spliced introns present in the relatively simple genome of Giardia lamblia provides an excellent opportunity to explore this question.

Results

Our investigation found three types of distribution patterns of the few introns in the intron-containing genes: ancient intron in ancient gene, later-evolved intron in ancient gene, and later-evolved intron in later-evolved gene, which can reflect to some extent the dynamic evolution of introns in Giardia. Without finding any special features or functional importance of these introns responsible for their retention, we noticed and experimentally verified that some intron-containing genes form sense-antisense gene pairs with transcribable genes on their complementary strands, and that the introns just reside in the overlapping regions.

Conclusions

In Giardia’s evolution, despite constant evolutionary selection pressure of intron loss, intron gain can still occur in both ancient and later-evolved genes, but only a few introns are retained; at least the evolutionary retention of some of the introns might not be due to the functional constraint of the introns themselves but the causes outside of introns, such as the constraints imposed by other genomic functional elements overlapping with the introns. These findings can not only provide some clues to find new genomic functional elements -- in the areas overlapping with introns, but suggest that “functional constraint” of introns may not be necessarily directly associated with intron loss and gain, and that the real functions are probably still outside of our current knowledge.

Reviewers

This article was reviewed by Mikhail Gelfand, Michael Gray, and Igor Rogozin.

     

Separate Introns Gained within Short and Long Soluble Peridinin-Chlorophyll a-Protein Genes during Radiation of Symbiodinium (Dinophyceae) Clade A and B Lineages

Here we document introns in two Symbiodinium clades that were most likely gained following divergence of this genus from other peridinin-containing dinoflagellate lineages. Soluble peridinin-chlorophyll a-proteins (sPCP) occur in short and long forms in different species. Duplication and fusion of short sPCP genes produced long sPCP genes. All short and long sPCP genes characterized to date, including those from free living species and Symbiodinium sp. 203 (clade C/type C2) are intronless. However, we observed that long sPCP genes from two Caribbean Symbiodinium clade B isolates each contained two introns. To test the hypothesis that introns were gained during radiation of clade B, we compared sPCP genomic and cDNA sequences from 13 additional distinct Caribbean and Pacific Symbiodinium clade A, B, and F isolates. Long sPCP genes from all clade B/B1 and B/B19 descendants contain orthologs of both introns. Short sPCP genes from S. pilosum (A/A2) and S. muscatinei (B/B4) plus long sPCP genes from S. microadriaticum (A/A1) and S. kawagutii (F/F1) are intronless. Short sPCP genes of S. microadriaticum have a third unique intron. Symbiodinium clade B long sPCP sequences are useful for assessing divergence among B1 and B19 descendants. Phylogenetic analyses of coding sequences from four dinoflagellate orders indicate that introns were gained independently during radiation of Symbiodinium clades A and B. Long sPCP introns were present in the most recent common ancestor of Symbiodinium clade B core types B1 and B19, which apparently diverged sometime during the Miocene. The clade A short sPCP intron was either gained by S. microadriaticum or possibly by the ancestor of Symbiodinium types A/A1, A3, A4 and A5. The timing of short sPCP intron gain in Symbiodinium clade A is less certain. But, all sPCP introns were gained after fusion of ancestral short sPCP genes, which we confirm as occurring once in dinoflagellate evolution.     

The role of introns in the conservation of the metabolic genes of Arabidopsis thaliana

In Arabidopsis thaliana, primary metabolic genes (PMGs) are more evolutionarily conserved and intron-rich than secondary metabolic genes. We observed that PMGs are more primitive and pan-taxonomically persistent as compared to secondary (SMGs) and non-metabolic genes (NMGs). This difference in primitiveness and persistence is primarily correlated with intron number and is independent of gene expression level. We propose a twofold explanation behind higher intron enrichment in PMGs. Firstly, introns might increase protein versatility amongst PMGs through alternative splicing, providing selective advantage of PMGs and making them more persistent across diverse plant taxa. Also, multifunctional PMGs may acquire functional domains by increasing the intronic burden. Additionally, single nucleotide polymorphisms (SNPs) accumulate at a higher rate in introns as compared to exons. Moreover, a strong negative correlation between cumulative exonic SNPs density and intron number indicates that introns may protect the exonic regions against the deleterious effect of these mutations, making them more conserved.     

Test of intron predictions reveals novel splice sites, alternatively spliced mRNAs and new introns in meiotically regulated genes of yeast

Correct identification of all introns is necessary to discern the protein-coding potential of a eukaryotic genome. The existence of most of the spliceosomal introns predicted in the genome of Saccharomyces cerevisiae remains unsupported by molecular evidence. We tested the intron predictions for 87 introns predicted to be present in non-ribosomal protein genes, more than a third of all known or suspected introns in the yeast genome. Evidence supporting 61 of these predictions was obtained, 20 predicted intron sequences were not spliced and six predictions identified an intron-containing region but failed to specify the correct splice sites, yielding a successful prediction rate of <80%. Alternative splicing has not been previously described for this organism, and we identified two genes (YKL186C/MTR2 and YML034W) which encode alternatively spliced mRNAs; YKL186C/MTR2 produces at least five different spliced mRNAs. One gene (YGR225W/SPO70) has an intron whose removal is activated during meiosis under control of the MER1 gene. We found eight new introns, suggesting that numerous introns still remain to be discovered. The results show that correct prediction of introns remains a significant barrier to understanding the structure, function and coding capacity of eukaryotic genomes, even in a supposedly simple system like yeast.     

Extensive structural conservation exists among several homologs of two Euglena chloroplast group II introns

82 of the 155 chloroplast introns in Euglena gracilis have been categorized as group II introns. Because they are shorter and more divergent than group II introns from other organisms, the assignment of these Euglena introns to the group II class has been questioned. In the current study, two homologs of E. gracilispetB intron 1 and four homologs of psbC intron 2 have been isolated from related species and characterized. Based on a comparative sequence analysis of intron homologs, the intron core and four of the six helical domains present in the canonical group II intron structural model are conserved in E. gracilispetB intron 1 and psbC intron 2 and all of their homologs. Distal portions of domain I, which are involved in most of the tertiary interactions, are less well conserved than the central core.     

Cloning,sequencing and expression of the two genes encoding the mitochondrial single-stranded DNA-binding protein in Xenopus laevis

In Xenopus laevis the single-stranded DNA binding protein imported into the mitochondria consists of two highly related polypeptides. The establishment of the genomic nucleotide sequences reveals that they are encoded by two different genes, XLSSB1 and XLSSB2. The deduced amino acid sequence is identical to the direct amino acid sequence determined by Edman degradation of the mitochondrial polypeptides [Ghrir, R., Lecaer, J.P., Dufresne, C. and Gueride, M. (1991) Primary structure of the two variants of Xenopus laevis mtSSB, a mitochondrial DNA binding protein. Arch. Biochem. Biophys. 291, 395–400]. Both genes are organized in seven exons and six introns, the sequence of the peptide leader is interrupted by an intervening sequence (intron 2). The exon/intron junctions are in exactly conserved positions, splitting the same codon. A high level of identity is observed between corresponding introns of the two genes over part or most of their lengths. Structural features of intronic sequences reveal multiple rearrangements and exchanges during the evolution of X. laevis species. A CCAAT box and the potential regulatory elements NRF-2 and Sp 1 are observed in the 5′-flanking region of both genes. During oogenesis, XLSSB gene expression is correlated with the replicative activity of the mitochondrial DNA.