Intron gain and loss in the evolution of the conserved eukaryotic recombination machinery

Is there any intron sliding in mammals?

BACKGROUND

Eukaryotic protein-coding genes consist of exons and introns. Exon-intron borders are conserved between species and thus their changes might be observed only on quite long evolutionary distances. One of the rarest types of change, in which intron relocates over a short distance, is called "intron sliding", but the reality of this event has been debated for a long time. The main idea of a search for intron sliding is to use the most accurate genome annotation and genome sequence, as well as high-quality transcriptome data. We applied them in a search for sliding introns in mammals in order to widen knowledge about the presence or absence of such phenomena in this group.

RESULTS

We didn't find any significant evidence of intron sliding in the primate group (human, chimpanzee, rhesus macaque, crab-eating macaque, green monkey, marmoset). Only one possible intron sliding event supported by a set of high quality transcriptomes was observed between EIF1AX human and sheep gene orthologs. Also, we checked a list of previously observed intron sliding events in mammals and showed that most likely they are artifacts of genome annotations and are not shown in subsequent annotation versions as well as are not supported by transcriptomic data.

CONCLUSIONS

We assume that intron sliding is indeed a very rare evolutionary event if it exists at all. Every case of intron sliding needs a lot of supportive data for detection and confirmation.

Introns: Good Day Junk Is Bad Day Treasure

Introns are ubiquitous in eukaryotic transcripts. They are often viewed as junk RNA but the huge energetic burden of transcribing, removing, and degrading them suggests a significant evolutionary advantage. Ostensibly, an intron functions within the host pre-mRNA to regulate its splicing, transport, and degradation. However, recent studies have revealed an entirely new class of trans-acting functions where the presence of intronic RNA in the cell impacts the expression of other genes in trans. Here, we review possible new mechanisms of intron functions, with a focus on the role of yeast introns in regulating the cell growth response to starvation.

GeMoMa: Homology-Based Gene Prediction Utilizing Intron Position Conservation and RNA-seq Data

GeMoMa is a homology-based gene prediction program that predicts gene models in target species based on gene models in evolutionary related reference species. GeMoMa utilizes amino acid sequence conservation, intron position conservation, and RNA-seq data to accurately predict protein-coding transcripts. Furthermore, GeMoMa supports the combination of predictions based on several reference species allowing to transfer high-quality annotation of different reference species to a target species. Here, we present a detailed description of GeMoMa modules and the GeMoMa pipeline and how they can be used on the command line to address particular biological problems.

Genomics and evolutionary aspect of calcium signaling event in calmodulin and calmodulin-like proteins in plants

BACKGROUND

Ca2+ ion is a versatile second messenger that operate in a wide ranges of cellular processes that impact nearly every aspect of life. Ca2+ regulates gene expression and biotic and abiotic stress responses in organisms ranging from unicellular algae to multi-cellular higher plants through the cascades of calcium signaling processes.

RESULTS

In this study, we deciphered the genomics and evolutionary aspects of calcium signaling event of calmodulin (CaM) and calmodulin like- (CML) proteins. We studied the CaM and CML gene family of 41 different species across the plant lineages. Genomic analysis showed that plant encodes more calmodulin like-protein than calmodulins. Further analyses showed, the majority of CMLs were intronless, while CaMs were intron rich. Multiple sequence alignment showed, the EF-hand domain of CaM contains four conserved D-x-D motifs, one in each EF-hand while CMLs contain only one D-x-D-x-D motif in the fourth EF-hand. Phylogenetic analysis revealed that, the CMLs were evolved earlier than CaM and later diversified. Gene expression analysis demonstrated that different CaM and CMLs genes were express differentially in different tissues in a spatio-temporal manner.

CONCLUSION

In this study we provided in detailed genome-wide identifications and characterization of CaM and CML protein family, phylogenetic relationships, and domain structure. Expression study of CaM and CML genes were conducted in Glycine max and Phaseolus vulgaris. Our study provides a strong foundation for future functional research in CaM and CML gene family in plant kingdom.

Using intron position conservation for homology-based gene prediction

Annotation of protein-coding genes is very important in bioinformatics and biology and has a decisive influence on many downstream analyses. Homology-based gene prediction programs allow for transferring knowledge about protein-coding genes from an annotated organism to an organism of interest.Here, we present a homology-based gene prediction program called GeMoMa. GeMoMa utilizes the conservation of intron positions within genes to predict related genes in other organisms. We assess the performance of GeMoMa and compare it with state-of-the-art competitors on plant and animal genomes using an extended best reciprocal hit approach. We find that GeMoMa often makes more precise predictions than its competitors yielding a substantially increased number of correct transcripts. Subsequently, we exemplarily validate GeMoMa predictions using Sanger sequencing. Finally, we use RNA-seq data to compare the predictions of homology-based gene prediction programs, and find again that GeMoMa performs well.Hence, we conclude that exploiting intron position conservation improves homology-based gene prediction, and we make GeMoMa freely available as command-line tool and Galaxy integration.

Using intron position conservation for homology-based gene prediction

Annotation of protein-coding genes is very important in bioinformatics and biology and has a decisive influence on many downstream analyses. Homology-based gene prediction programs allow for transferring knowledge about protein-coding genes from an annotated organism to an organism of interest.Here, we present a homology-based gene prediction program called GeMoMa. GeMoMa utilizes the conservation of intron positions within genes to predict related genes in other organisms. We assess the performance of GeMoMa and compare it with state-of-the-art competitors on plant and animal genomes using an extended best reciprocal hit approach. We find that GeMoMa often makes more precise predictions than its competitors yielding a substantially increased number of correct transcripts. Subsequently, we exemplarily validate GeMoMa predictions using Sanger sequencing. Finally, we use RNA-seq data to compare the predictions of homology-based gene prediction programs, and find again that GeMoMa performs well.Hence, we conclude that exploiting intron position conservation improves homology-based gene prediction, and we make GeMoMa freely available as command-line tool and Galaxy integration.

Pax2/5/8 and Pax6 alternative splicing events in basal chordates and vertebrates: a focus on paired box domain

Paired box transcription factors play important role in development and tissue morphogenesis. The number of Pax homologs varies among species studied so far, due to genome and gene duplications that have affected PAX family to a great extent. Based on sequence similarity and functional domains, four Pax classes have been identified in chordates, namely Pax1/9, Pax2/5/8, Pax3/7, and Pax4/6. Numerous splicing events have been reported mainly for Pax2/5/8 and Pax6 genes. Of significant interest are those events that lead to Pax proteins with presumed novel properties, such as altered DNA-binding or transcriptional activity. In the current study, a thorough analysis of Pax2/5/8 splicing events from cephalochordates and vertebrates was performed. We focused more on Pax2/5/8 and Pax6 splicing events in which the paired domain is involved. Three new splicing events were identified in Oryzias latipes, one of which seems to be conserved in Acanthomorphata. Using representatives from deuterostome and protostome phyla, a comparative analysis of the Pax6 exon-intron structure of the paired domain was performed, during an attempt to estimate the time of appearance of the Pax6(5a) mRNA isoform. As shown in our analysis, this splicing event is characteristic of Gnathostomata and is absent in the other chordate subphyla. Moreover, expression pattern of alternative spliced variants was compared between cephalochordates and fish species. In summary, our data indicate expansion of alternative mRNA variants in paired box region of Pax2/5/8 and Pax6 genes during the course of vertebrate evolution.

The splicing fate of plant SPO11 genes

Toward the global understanding of plant meiosis, it seems to be essential to decipher why all as yet sequenced plants need or at least encode for two different meiotic SPO11 genes. This is in contrast to mammals and fungi, where only one SPO11 is present. Both SPO11 in Arabidopsis thaliana are essential for the initiation of double strand breaks (DSBs) during the meiotic prophase. In nearly all eukaryotic organisms DSB induction during prophase I by SPO11 leads to meiotic DSB repair, thereby ensuring the formation of a necessary number of crossovers (CO) as physical connections between the homologous chromosomes. We aim to investigate the specific functions and evolution of both SPO11 genes in land plants. Therefore, we identified and cloned the respective orthologous genes from Brassica rapa, Carica papaya, Oryza sativa, and Physcomitrella patens. In parallel we determined the full length cDNA sequences of SPO11-1 and -2 from all of these plants by RT-PCR. During these experiments we observed that the analyzed plants exhibit a pattern of alternative splicing products of both SPO11 mRNAs. Such an aberrant splicing has previously been described for Arabidopsis and therefore seems to be conserved throughout evolution. Most of the splicing forms of SPO11-1 and -2 seem to be non-functional as they either showed intron retention (IR) or shortened exons. However, the positional distribution and number of alternative splicing events vary strongly between the different plants. The cDNAs showed in most cases premature termination codons (PTCs) due to frameshift. Nevertheless, in some cases we found alternatively spliced but functional cDNAs. These findings let us suggest that alternative splicing of SPO11 depends on the respective gene sequence and on the plant species. Therefore, this conserved mechanism might play a role concerning regulation of SPO11.

Progenitor-derivative relationships of Hordeum polyploids (Poaceae, Triticeae) inferred from sequences of TOPO6, a nuclear low-copy gene region

Polyploidization is a major mechanism of speciation in plants. Within the barley genus Hordeum, approximately half of the taxa are polyploids. While for diploid species a good hypothesis of phylogenetic relationships exists, there is little information available for the polyploids (4×, 6×) of Hordeum. Relationships among all 33 diploid and polyploid Hordeum species were analyzed with the low-copy nuclear marker region TOPO6 for 341 Hordeum individuals and eight outgroup species. PCR products were either directly sequenced or cloned and on average 12 clones per individual were included in phylogenetic analyses. In most diploid Hordeum species TOPO6 is probably a single-copy locus. Most sequences found in polyploid individuals phylogenetically cluster together with sequences derived from diploid species and thus allow the identification of parental taxa of polyploids. Four groups of sequences occurring only in polyploid taxa are interpreted as footprints of extinct diploid taxa, which contributed to allopolyploid evolution. Our analysis identifies three key species involved in the evolution of the American polyploids of the genus. (i) All but one of the American tetraploids have a TOPO6 copy originating from the Central Asian diploid H. roshevitzii, the second copy clustering with different American diploid species. (ii) All hexaploid species from the New World have a copy of an extinct close relative of H. californicum and (iii) possess the TOPO6 sequence pattern of tetraploid H. jubatum, each with an additional copy derived from different American diploids. Tetraploid H. bulbosum is an autopolyploid, while the assumed autopolyploid H. brevisubulatum (4×, 6×) was identified as allopolyploid throughout most of its distribution area. The use of a proof-reading DNA polymerase in PCR reduced the proportion of chimerical sequences in polyploids in comparison to Taq polymerase.

Two extinct diploid progenitors were involved in allopolyploid formation in the Hordeum murinum (Poaceae: Triticeae) taxon complex

Wall barley (Hordeum murinum) occurs with three subspecies, naturally distributed from southern Central Asia through the Mediterranean region to northwestern Europe, but now is an invasive weed in many parts of the world. Subspecies glaucum is diploid, while subspp. murinum and leporinum are tetraploids, the latter also occurring with a hexaploid cytotype. Earlier analyses were inconclusive regarding auto- or allopolyploid origins of subspp. murinum and leporinum. We analyzed the phylogeny of the taxon group using amplified fragment length polymorphisms (AFLP), sequences of cloned PCR products of the nuclear ribosomal DNA internal transcribed spacer region (ITS), a part of the nuclear single-copy gene topoisomerase 6 (Topo6) spanning two introns, and sequences of the chloroplast trnL-F region together with length variation at six chloroplast microsatellite loci, including multiple individuals of each subspecies and cytotype, covering the entire natural distribution area of the species. Phylogenetic analyses with all used markers differentiate diploid and polyploids. Sequences of both nuclear regions indicated that diploid subsp. glaucum was involved in tetraploid formation together with a now extinct species belonging to the same Hordeum genome group (Xu). Furthermore, AFLP and ITS analyses suggest that a third, though closely related extinct taxon contributed to hexaploid formation. No method was able to discern tetraploid subspp. murinum and leporinum, which we attribute to the young age of subsp. murinum. None of the used molecular markers revealed a strong geographic pattern of genetic variation that would allow comprehensive phylogeographic analysis, most probably due to the very effective seed dispersal of the taxa.

Two distinct MUS81-EME1 complexes from Arabidopsis process Holliday junctions

The MUS81 endonuclease complex has been shown to play an important role in the repair of stalled or blocked replication forks and in the processing of meiotic recombination intermediates from yeast to humans. This endonuclease is composed of two subunits, MUS81 and EME1. Surprisingly, unlike other organisms, Arabidopsis (Arabidopsis thaliana) has two EME1 homologs encoded in its genome. AtEME1A and AtEME1B show 63% identity on the protein level. We were able to demonstrate that, after expression in Escherichia coli, each EME1 protein can assemble with the unique AtMUS81 to form a functional endonuclease. Both complexes, AtMUS81-AtEME1A and AtMUS81-AtEME1B, are not only able to cleave 3'-flap structures and nicked Holliday junctions (HJs) but also, with reduced efficiency, intact HJs. While the complexes have the same cleavage patterns with both nicked DNA substrates, slight differences in the processing of intact HJs can be detected. Our results are in line with an involvement of both MUS81-EME1 endonuclease complexes in DNA recombination and repair processes in Arabidopsis.

Genome-wide analysis of recombination machinery for spliceosomal introns gain

What caused spliceosomal introns gain remains an unsolved problem. To this, defining what spliceosomal introns arise from is critical. Here, the introns density of the genomes is calculated for four species, indicating:(1) sex chromosomes in mammals have lower intron densities, (2) despite that, the proportion of UTRs (untranslated regions) with introns in sex chromosomes is higher than other ones, and (3) AT content of introns is more similar to that of intergenic regions when these regions comprise the majority of a chromosome, and more similar to that of exons, when exons are the majority of the chromosome. On the other hand, introns have been clearly demonstrated to invade genetic sequences in recent times while sex chromosomes evolved from a pair of autosomes within the last 300 millions years. One main difference between sex chromosomes and autosomes in mammalian is that sex chromosomes recombination stopped. Thus, recombination might be the main determinant for eukaryotes gaining spliceosomal introns. To further prove that and avoid giving weak signal, the whole genomes from eight eukaryotic species are analyzed and present strong signal for above the trend (3) in three species (t-test, P = 0.55 for C. elegans, P = 0.72 for D. melanogaster and P = 0.83 for A. thaliana). These results suggest that the genome-wide coincidence as above (3) can only be caused by the large-scale random unequal crossover in eukaryote meiosis, which might have fueled spliceosomal introns but hardly occurred in prokaryotes.

A homolog of ScRAD5 is involved in DNA repair and homologous recombination in Arabidopsis

Rad5 is the key component in the Rad5-dependent error-free branch of postreplication repair in yeast (Saccharomyces cerevisiae). Rad5 is a member of the Snf2 ATPase/helicase family, possessing as a characteristic feature, a RING-finger domain embedded in the Snf2-helicase domain and a HIRAN domain. Yeast mutants are sensitive to DNA-damaging agents and reveal differences in homologous recombination. By sequence comparisons we were able to identify two homologs (AtRAD5a and AtRAD5b) in the Arabidopsis thaliana genome, sharing about 30% identity and 45% similarity to yeast Rad5. AtRad5a and AtRad5b have the same kind of domain organization with a higher degree of similarity to each other than to ScRad5. Surprisingly, both genes differ in function: whereas two independent mutants of Atrad5a are hypersensitive to the cross-linking agents mitomycin C and cis-platin and to a lesser extent to the methylating agent, methyl methane sulfonate, the Atrad5b mutants did not exhibit any sensitivity to all DNA-damaging agents tested. An Atrad5a/Atrad5b double mutant resembles the sensitivity phenotype of the Atrad5a single mutants. Moreover, in contrast to Atrad5b, the two Atrad5a mutants are deficient in homologous recombination after treatment with the double-strand break-inducing agent bleomycin. Our results suggest that the RAD5-dependent error-free branch of postreplication repair is conserved between yeast and plants, and that AtRad5a might be functionally homologous to ScRad5.

Characterization of the gene Mre11 and evidence of silencing after polyploidization in Triticum

The MRE11 protein is a component of the highly conserved MRN complex, along with RAD50 and NBS1. This complex is crucial in the repair of breaks in double stranded DNA, and is involved in many other cell processes. The present paper reports the molecular characterization of Mre11 gene in all three genomes of wheat, making use of the diploid species Triticum monococcum (genome A) and Aegilops Tauschii (genome D), the tetraploid T. turgidum (genomes A and B), and the hexaploid T. aestivum (genomes A, B and D). The genomic sequences characterized ranged from 4,662 to 4,766 bp in length; the cDNA corresponding to the processed mRNA was 2,440-2,510 bp long. In all cases, Mre11 coded for a highly conserved protein of 699 amino acids with a structure involving 22 exons. Mre11 expression was determined by real-time PCR in all the species analysed. The tetraploid species showed an expression similar to that of the diploid Ae. tauschii and lower than that of T. monococcum. Stronger expression was detected in the hexaploid T. aestivum. The SSCP technique was modified by introducing fluorescent labelling to the procedure in order to analyse the expression of the different Mre11 genes (i.e., those belonging to the different genomes) in the polyploid species. In both polyploids, the Mre11 gene belonging to the B genome was the least expressed. This probably reflects a first step in the process of silencing duplicate genes after polyploidization.

Overexpression of putative topoisomerase 6 genes from rice confers stress tolerance in transgenic Arabidopsis plants

DNA topoisomerase 6 (TOP6) belongs to a novel family of type II DNA topoisomerases present, other than in archaebacteria, only in plants. Here we report the isolation of full-length cDNAs encoding putative TOP6 subunits A and B from rice (Oryza sativa ssp. indica), preserving all the structural domains conserved among archaebacterial TOP6 homologs and eukaryotic meiotic recombination factor SPO11. OsTOP6A1 was predominantly expressed in prepollinated flowers. The transcript abundance of OsTOP6A2, OsTOP6A3 and OsTOP6B was also higher in prepollinated flowers and callus. The expression of OsTOP6A2, OsTOP6A3 and OsTOP6B was differentially regulated by the plant hormones, auxin, cytokinin, and abscisic acid. Yeast two-hybrid analysis revealed that the full-length OsTOP6B protein interacts with both OsTOP6A2 and OsTOP6A3, but not with OsTOP6A1. The nuclear localization of OsTOP6A3 and OsTOP6B was established by the transient expression of their beta-glucuronidase fusion proteins in onion epidermal cells. Overexpression of OsTOP6A3 and OsTOP6B in transgenic Arabidopsis plants conferred reduced sensitivity to the stress hormone, abscisic acid, and tolerance to high salinity and dehydration. Moreover, the stress tolerance coincided with enhanced induction of many stress-responsive genes in transgenic Arabidopsis plants. In addition, microarray analysis revealed that a large number of genes are expressed differentially in transgenic plants. Taken together, our results demonstrate that TOP6 genes play a crucial role in stress adaptation of plants by altering gene expression.

The role of AtMUS81 in DNA repair and its genetic interaction with the helicase AtRecQ4A

The endonuclease MUS81 has been shown in a variety of organisms to be involved in DNA repair in mitotic and meiotic cells. Homologues of the MUS81 gene exist in the genomes of all eukaryotes, pointing to a conserved role of the protein. However, the biological role of MUS81 varies between different eukaryotes. For example, while loss of the gene results in strongly impaired fertility in Saccharomyces cerevisiae and nearly complete sterility in Schizosaccharomyces pombe, it is not essential for meiosis in mammals. We identified a functional homologue (AtMUS81/At4g30870) in the genome of Arabidopsis thaliana and isolated a full-length cDNA of this gene. Analysing two independent T-DNA insertion lines of AtMUS81, we found that they are sensitive to the mutagens MMS and MMC. Both mutants have a deficiency in homologous recombination in somatic cells but only after induction by genotoxic stress. In contrast to yeast, no meiotic defect of AtMUS81 mutants was detectable and the mutants are viable. Crosses with a hyperrecombinogenic mutant of the AtRecQ4A helicase resulted in synthetic lethality in the double mutant. Thus, the nuclease AtMUS81 and the helicase AtRecQ4A seem to be involved in two alternative pathways of resolution of replicative DNA structures in somatic cells.

Genetic structure and evolution of the Vps25 family, a yeast ESCRT-II component

Background

Vps25p is the product of yeast gene VPS25 and is found in an endosomal sorting complex required for transport (ESCRT)-II, along with Vps22p and Vps36p. This complex is essential for sorting of ubiquitinated biosynthetic and endosomal cargoes into endosomes.

Results

We found that VPS25 is a highly conserved and widely expressed eukaryotic gene, with single orthologs in chromalveolate, excavate, amoebozoan, plant, fungal and metazoan species. Two paralogs were found in Trichomonas vaginalis. An ortholog was strikingly absent from the Encephalitozoon cuniculi genome. Intron positions were analyzed in VPS25 from 36 species. We found evidence for five ancestral VPS25 introns, intron loss, and single instances of intron gain (a Paramecium species) and intron slippage (Theileria species). Processed pseudogenes were identified in four mammalian genomes, with a notable absence in the mouse genome. Two retropseudogenes were found in the chimpanzee genome, one more recently inserted, and one evolving from a common primate ancestor. The amino acid sequences of 119 Vps25 orthologs are aligned, compared with the known secondary structure of yeast Vps25p, and used to carry out phylogenetic analysis. Residues in two amino-terminal PPXY motifs (motif I and II), involved in dimerization of Vps25p and interaction with Vps22p and Vps36p, were closely, but not absolutely conserved. Specifically, motif I was absent in Vps25 homologs of chromalveolates, euglenozoa, and diplomonads. A highly conserved carboxy-terminal lysine was identified, which suggests Vps25 is ubiquitinated. Arginine-83 of yeast Vps25p involved in Vps22p interaction was highly, but not absolutely, conserved. Human tissue expression analysis showed universal expression.

Conclusion

We have identified 119 orthologs of yeast Vps25p. Expression of mammalian VPS25 in a wide range of tissues, and the presence in a broad range of eukaryotic species, indicates a basic role in eukaryotic cell function. Intron splice site positions were highly conserved across all major eukaryotic species, suggesting an ancestral origin. Amino acid sequence analysis showed the consensus for the amino-terminal proline-rich motifs is P- [WP]-X-[YF] for motif I (when present) and P-P-[FYL]-[FY] for motif II, and that Vps25 may be ubiquitinated.

The RecQ gene family in plants

RecQ helicases are conserved throughout all kingdoms of life regarding their overall structure and function. They are 3'-5' DNA helicases resolving different recombinogenic DNA structures. The RecQ helicases are key factors in a number of DNA repair and recombination pathways involved in the maintenance of genome integrity. In eukaryotes the number of RecQ genes and the structure of RecQ proteins vary strongly between organisms. Therefore, they have been named RecQ-like genes. Knockouts of several RecQ-like genes cause severe diseases in animals or harmful cellular phenotypes in yeast. Until now the largest number of RecQ-like genes per organism has been found in plants. Arabidopsis and rice possess seven different RecQ-like genes each. In the almost completely sequenced genome of the moss Physcomitrella patens at least five RecQ-like genes are present. One of the major present and future research aims is to define putative plant-specific functions and to assign their roles in DNA repair and recombination pathways in relation to RecQ genes from other eukaryotes. Regarding their intron positions, the structures of six RecQ-like genes of dicots and monocots are virtually identical indicating a conservation over a time scale of 150 million years. In contrast to other eukaryotes one gene (RecQsim) exists exclusively in plants. It possesses an interrupted helicase domain but nevertheless seems to have maintained the RecQ function. Owing to a recent gene duplication besides the AtRecQl4A gene an additional RecQ-like gene (AtRecQl4B) exists in the Brassicaceae only. Genetic studies indicate that a AtRecQl4A knockout results in sensitivity to mutagens as well as an hyper-recombination phenotype. Since AtRecQl4B was still present, both genes must have non-redundant roles. Analysis of plant RecQ-like genes will not only increase the knowledge on DNA repair and recombination, but also on the evolution and radiation of protein families.

Rates of intron loss and gain: implications for early eukaryotic evolution

We study the intron-exon structures of 684 groups of orthologs from seven diverse eukaryotic genomes and provide maximum likelihood estimates for rates and numbers of intron losses and gains in these same genes for a variety of lineages. Rates of intron loss vary from approximately 2 x 10(-9) to 2 x 10(-10) per year. Rates of gain vary from 6 x 10(-13) to 4 x 10(-12) per possible intron insertion site per year. There is an inverse correspondence between rates of intron loss and gain, leading to a 20-fold variation among lineages in the ratio of the rates of the two processes. The observed rates of intron gain are insufficient to explain the large number of introns estimated to have been present in the plant-animal ancestor, suggesting that introns present in early eukaryotes may have been created by a fundamentally different process than more recently gained introns.

A phylogenomic approach to reconstructing the diversification of serine proteases in fungi

Using a phylogenomic approach with 10 fungi of very different virulence and habitat, we determined that there was substantial diversification of subtilase-type proteases early in ascomycete history (with subsequent loss in many lineages) but with no comparable diversification of trypsins. Patterns of intron loss and the degree of divergence between paralogues demonstrated that the proliferation of proteinase K subtilases and subtilisin type subtilases seen in pathogenic ascomycetes (Metarhizium anisopliae, Magnaporthe grisea, Fusarium graminearum) occurred after the basidiomycete/ascomycete split but predated radiation of ascomycete lineages. This suggests that the early ascomycetes had a lifestyle that selected for multiple proteases, whereas the current disparity in gene numbers between ascomycete lineages results from retention of genes in at least some pathogens that have been lost in other lineages (yeasts, Aspergillus nidulans, Neurospora crassa). A similar prevailing trend towards lineage specific gene loss of trypsins in saprophytes and some pathogens suggests that their phylogenetic breadth will have been much wider in early fungi than currently.

Patterns of intron gain and loss in fungi

Little is known about the patterns of intron gain and loss or the relative contributions of these two processes to gene evolution. To investigate the dynamics of intron evolution, we analyzed orthologous genes from four filamentous fungal genomes and determined the pattern of intron conservation. We developed a probabilistic model to estimate the most likely rates of intron gain and loss giving rise to these observed conservation patterns. Our data reveal the surprising importance of intron gain. Between about 150 and 250 gains and between 150 and 350 losses were inferred in each lineage. We discuss one gene in particular (encoding 1-phosphoribosyl-5-pyrophosphate synthetase) that displays an unusually high rate of intron gain in multiple lineages. It has been recognized that introns are biased towards the 5' ends of genes in intron-poor genomes but are evenly distributed in intron-rich genomes. Current models attribute this bias to 3' intron loss through a poly-adenosine-primed reverse transcription mechanism. Contrary to standard models, we find no increased frequency of intron loss toward the 3' ends of genes. Thus, recent intron dynamics do not support a model whereby 5' intron positional bias is generated solely by 3'-biased intron loss.

Epizootiologic survey for Babesia microti among small wild mammals in northeastern Eurasia and a geographic diversity in the beta-tubulin gene sequences

We previously reported that small wild rodents in Japan harbor two types of novel Babesia microti-like parasites (designated as Hobetsu and Kobe types), but not the type commonly found in the northeastern United States (U.S. type) where human babesiosis is endemic. To determine whether these new types of parasites are distributed in places surrounding Japan, an epizootiologic survey was undertaken in three geographically distant areas in northeastern Eurasia; South Korea, Vladivostok in Russia, and Xinjiang in China. Blood samples were collected from a total of 387 animals comprising 24 species. DNAs extracted from the samples were tested by nested PCR targeting babesial nuclear small-subunit rRNA gene (rDNA), which revealed that small rodents harboring B. microti exist in all three survey areas. Sequence analysis showed that all PCR-positive samples had rDNA sequences virtually identical to that of U.S.-type B. microti. However, when beta-tubulin gene sequences were compared, evident geographic variations were seen. By use of primers specific for each of the beta-tubulin genes of Kobe-, Hobetsu-, and U.S.-type parasites, a type-specific PCR was developed. Parasite with Hobetsu- or Kobe-type sequence was not detected from any of the three survey areas. These findings suggest that U.S.-type B. microti is widely distributed among small wild mammals in temperate zones of not only North America, but also Eurasia, whereas that Hobetsu- and Kobe-type parasites may be uniquely distributed in Japan.

Recent Intron Gain in Elongation Factor-1α of Colletid Bees (Hymenoptera: Colletidae)

We discovered the presence of a unique spliceosomal intron in the F1 copy of elongation factor-1alpha (EF-1alpha) restricted to the bee family Colletidae (Hymenoptera: Apoidae). The intron ranges in size from 101 to 1044 bp and shows no positional sliding. Our data also demonstrate the complete absence of this intron from exemplars representing all other bee families, as well as from close hymenopteran relatives. A review of the literature finds that this intron is likewise absent from all other arthropods for which data are available. This provides unambiguous evidence for a relatively recent intron insertion event in the colletid common ancestor and, at least in this specific instance, lends support to the introns-late hypothesis. The comparative distribution of this novel intron also supports the monophyly of Colletidae and the exclusion of the Stenotritidae from this family, providing an example of the potential of some introns to act as robust markers of shared descent.

Gene structure and evolution of Tieg3, a new member of the Tieg family of proteins

TGF beta-inducible immediate early gene, Tieg, belongs to the superfamily of Sp1-like transcription factors containing three C(2)H(2)-zinc finger DNA binding motifs close to the C-terminus. So far, Tieg1 and Tieg2 have been identified in human and mouse. We identified Tieg3, a new member of the Tieg protein family by screening a mouse cDNA library. Tieg3 has almost all the known features of the Tieg protein family: it shares a highly conserved C(2)H(2) zinc finger DNA binding domain and is 96% identical to Tieg2 and 86% to Tieg1, respectively. In addition, the three repression domains at the N-terminus, R1, R2 and R3 are conserved in all the Tiegs. Similar to Tieg1 and Tieg2, Tieg3 mRNA is up-regulated in response to TGF beta 1 treatment and can perform the Sp1 sites mediated repression of transcription. A 4 kilobase (kb) long transcript of mouse Tieg3 can be detected using Northern-blot analysis. The gene of mouse Tieg3 contains four exons. Due to the amino acid sequence similarity, mouse Tieg2 is regarded as an orthologue of human Tieg2. However, the mouse Tieg3 gene is localized in a conserved segment on mouse chromosome 12 corresponding to human Tieg2 on chromosome 2 with the same gene order. An interesting explanation for this apparent contradiction might be a homologous recombination leading to loci exchange between the mouse Tieg3 and Tieg2.

Intron requirement for AFP gene expression in Trichoderma viride

The 430 bp ORF of the Aspergillus giganteus antifungal protein (AFP) gene, containing two small introns, was fused between the promoter and the terminator of the Aspergillus nidulans trpC gene. The AFP gene in this vector produced detectable levels of spliced mRNA in Trichoderma viride. In contrast, in the same vector configuration, its 285 bp intronless derivative showed no accumulation of mRNA when transformed into T. viride. Such expression results were confirmed at the protein level. This fact demonstrated that the introns were required for AFP gene expression in T. viride. This is thought to be a novel phenomenon found in filamentous fungi. Although the mechanism of splicing in filamentous fungi might be similar to that in other eukaryotes, little is known of how it affects expression. This study suggests that the small introns in filamentous fungal genes may not only act as intervening elements, but may also play crucial roles in gene expression by affecting mRNA accumulation. Furthermore, it may provide new evidence for intron-dependent evolution.