Intron-rich Genes Research Articles

Identification and characterization of NHX gene family for their role under salt stress in Vigna mungo.

In the current study, we have performed a comprehensive analysis of the Sodium Hydrogen Exchanger (NHX) gene family in Vigna mungo, and a total of 44 NHX genes were identified. A bimodal distribution based on domains, gene structure and phylogenetic analysis was evident. All intronpoor and intron-rich genes were clustered in clades I and II, respectively. Interestingly, all genes of subclade IIb were localized to vacuoles and possess only the NHX domain. The isoelectric point and trans-membrane domain analysis reflect the wide distribution of the NHX genes. Interestingly, Vm_NHX2 and Vm_NHX3 lacked trans-membrane domain but were found to interact with other NHX genes as well as vital salinity pathway genes, including calcium-mediated salt-responsive genes. The comparison of the mRNA sequences with that of V. marina, a halophytic species, reflects their independent evolution, majorly supporting the convergent evolution. The Ka/Ks ratio reflects the abundance of purifying selection supporting their conserved function during evolution. In our analysis, several abiotic stress and hormone-responsive elements and transcription factor binding sites were present in the promoter of the NHX genes. Further, the ion partitioning of a tolerant (K90) and a susceptible (K49) variety of V. mungo suggested that K90 managed the Na+/K+ ratio more affluently, which was also supported by profiling of superoxide radicals, hydrogen peroxide, phenol, peroxidase activity and superoxide dismutase activity. From the expression, we identified five candidate Vm_NHX genes, four of which, i.e. Vm_NHX16, Vm_NHX17, Vm_NHX29 and Vm_NHX33, were localized to the vacuolar and lysosomal membrane.

Comparative mitochondrial genome analyses reveal conserved gene arrangement but massive expansion/contraction in two closely related Exserohilum pathogens

Exserohilum turcicum and E. rostratum, two closely related fungal species, are both economically important pathogens but have quite different target hosts (specific to plants and cross-kingdom infection, respectively). In the present study, complete circular mitochondrial genomes of the two Exserohilum species were sequenced and de novo assembled, which mainly comprised the same set of 13 core protein-coding genes (PCGs), two rRNAs, and a certain number of tRNAs and unidentified open reading frames (ORFs). Comparative analyses indicated that these two fungi had significant mitogenomic collinearity and consistent mitochondrial gene arrangement, yet with vastly different mitogenome sizes, 264,948 bp and 64,620 bp, respectively. By contrast with the 17 introns containing 17 intronic ORFs (one-to-one) in the E. rostratum mitogenome, E. turcicum involved far more introns (70) and intronic ORFs (126), which was considered as the main contributing factors of their mitogenome expansion/contraction. Within the generally intron-rich gene cox1, a total of 18 and 10 intron position classes (Pcls) were identified separately in the two mitogenomes. Moreover, 16.16% and 10.85% ratios of intra-mitogenomic repetitive regions were detected in E. turcicum and E. rostratum, respectively. Based on the combined mitochondrial gene dataset, we established a well-supported topology of phylogeny tree of 98 ascomycetes, implying that mitogenomes may act as an effective molecular marker for fungal phylogenetic reconstruction. Our results served as the first report on mitogenomes in the genus Exserohilum, and would have significant implications in understanding the origin, evolution and pathogenic mechanisms of this fungal lineage.

Mutational Effects of Mobile Introns on the Mitochondrial Genomes of Metschnikowia Yeasts

It has been argued that DNA repair by homologous recombination in the context of endonuclease-mediated cleavage can cause mutations. To better understand this phenomenon, we examined homologous recombination following endonuclease cleavage in a native genomic context: the movement of self-splicing introns in the mitochondrial genomes of Metschnikowia yeasts. Self-splicing mitochondrial introns are mobile elements, which can copy and paste themselves at specific insertion sites in mitochondrial DNA using a homing endonuclease in conjunction with homologous recombination. Here, we explore the mutational effects of self-splicing introns by comparing sequence variation within the intron-rich cox1 and cob genes from 71 strains (belonging to 40 species) from the yeast genus Metschnikowia. We observed a higher density of single nucleotide polymorphisms around self-splicing-intron insertion sites. Given what is currently known about the movement of organelle introns, it is likely that their mutational effects result from the high binding affinity of endonucleases and their interference with repair machinery during homologous recombination (or, alternatively, via gene conversion occurring during the intron insertion process). These findings suggest that there are fitness costs to harbouring self-splicing, mobile introns and will help us better understand the risks associated with modern biotechnologies that use endonuclease-mediated homologous recombination, such as CRISPR-Cas9 gene editing.

The emergence and evolution of intron-poor and intronless genes in intron-rich plant gene families.

Eukaryotic genes can be classified into intronless (no introns), intron-poor (three or fewer introns per gene) or intron-rich. Early eukaryotic genes were mostly intron-rich, and their alternative splicing into multiple transcripts, giving rise to different proteins, might have played pivotal roles in adaptation and evolution. Interestingly, extant plant genomes contain many gene families with one or sometimes few sub-families with genes that are intron-poor or intronless, and it remains unknown when and how these intron-poor or intronless genes have originated and evolved, and what their possible functions are. In this study, we identified 33 such gene families that contained intronless and intron-poor sub-families. Intronless genes seemed to have first emerged in early land plant evolution, while intron-poor sub-families seemed first to have appeared in green algae. In contrast to intron-rich genes, intronless genes in intron-poor sub-families occurred later, and were subject to stronger functional constraints. Based on RNA-seq analyses in Arabidopsis and rice, intronless or intron-poor genes in AP2, EF-hand_7, bZIP, FAD_binding_4, STE_STE11, CAMK_CAMKL-CHK1 and C2 gene families were more likely to play a role in response to drought and salt stress, compared with intron-rich genes in the same gene families, whereas intronless genes in the B_lectin and S_locus_glycop gene family were more likely to participate in epigenetic processes and plant development. Understanding the origin and evolutionary trajectory, as well as the potential functions, of intronless and intron-poor sub-families provides further insight into plant genome evolution and the functional divergence of genes.

Myosin II sequences for Lethocerusindicus

We present the genomic and expressed myosin II sequences from the giant waterbug, Lethocerus indicus. The intron rich gene appears relatively ancient and contains six regions of mutually exclusive exons that are alternatively spliced. Alternatively spliced regions may be involved in the asymmetric myosin dimer structure known as the interacting heads motif, as well as stabilizing the interacting heads motif within the thick filament. A lack of negative charge in the myosin S2 domain may explain why Lethocerus thick filaments display a perpendicular interacting heads motif, rather than one folded back to contact S2, as is seen in other thick filament types such as those from tarantula.

Recent gene multiplication and evolution of a juvenile hormone esterase-related gene in a lepidopteran pest

Juvenile hormone esterase (JHE) is a carboxylesterase that plays critical roles in regulating larval to adult transition by hydrolyzing the key developmental and reproductive hormone, juvenile hormone (JH). In the current work, we have cloned and sequenced a superfamily of juvenile hormone esterase related genes in Sesamia nonagrioides (JHERs). These seemed to have been recently multiplied from a common ancestral gene and consequently were inherited in the resulting populations as intron-less and intron-rich genes (loci). We sequenced three JHER genes (one intron-rich and two intron-less) and four cDNAs encoding for juvenile hormone esterase related sequences. Three cDNAs presented nucleotide deletions similar to alternative splicing events when compared with the introns of the intron-rich gene. The exons of the intron-rich gene were >98% identical with one of the intron-less gene and the homologous sequences of all the four cDNAs. Moreover, the second intron-less gene seemed to be almost identical with one of the four cDNAs. The fourth cDNA contained an extensive (in-frame) deletion inside its ORF. This mRNA seems to be encoded by another gene which's deletion was generated by homologous recombination. Interestingly, our data revealed differential expression patterns for the four cDNAs. This study provides an initial assessment of the diversity of JHER genes in a population of Sesamia and presents this species as an attractive model to study the diversification of JHE-like esterase genes and their functional consequences.

Multifunctional polyketide synthase genes identified by genomic survey of the symbiotic dinoflagellate, Symbiodinium minutum.

BackgroundDinoflagellates are unicellular marine and freshwater eukaryotes. They possess large nuclear genomes (1.5–245 gigabases) and produce structurally unique and biologically active polyketide secondary metabolites. Although polyketide biosynthesis is well studied in terrestrial and freshwater organisms, only recently have dinoflagellate polyketides been investigated. Transcriptomic analyses have characterized dinoflagellate polyketide synthase genes having single domains. The Genus Symbiodinium, with a comparatively small genome, is a group of major coral symbionts, and the S. minutum nuclear genome has been decoded.ResultsThe present survey investigated the assembled S. minutum genome and identified 25 candidate polyketide synthase (PKS) genes that encode proteins with mono- and multifunctional domains. Predicted proteins retain functionally important amino acids in the catalytic ketosynthase (KS) domain. Molecular phylogenetic analyses of KS domains form a clade in which S. minutum domains cluster within the protist Type I PKS clade with those of other dinoflagellates and other eukaryotes. Single-domain PKS genes are likely expanded in dinoflagellate lineage. Two PKS genes of bacterial origin are found in the S. minutum genome. Interestingly, the largest enzyme is likely expressed as a hybrid non-ribosomal peptide synthetase-polyketide synthase (NRPS-PKS) assembly of 10,601 amino acids, containing NRPS and PKS modules and a thioesterase (TE) domain. We also found intron-rich genes with the minimal set of catalytic domains needed to produce polyketides. Ketosynthase (KS), acyltransferase (AT), and acyl carrier protein (ACP) along with other optional domains are present. Mapping of transcripts to the genome with the dinoflagellate-specific spliced leader sequence, supports expression of multifunctional PKS genes. Metabolite profiling of cultured S. minutum confirmed production of zooxanthellamide D, a polyhydroxy amide polyketide and other unknown polyketide secondary metabolites.ConclusionThis genomic survey demonstrates that S. minutum contains genes with the minimal set of catalytic domains needed to produce polyketides and provides evidence of the modular nature of Type I PKS, unlike monofunctional Type I PKS from other dinoflagellates. In addition, our study suggests that diversification of dinoflagellate PKS genes comprises dinoflagellate-specific PKS genes with single domains, multifunctional PKS genes with KS domains orthologous to those of other protists, and PKS genes of bacterial origin.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-2195-8) contains supplementary material, which is available to authorized users.

Alternative splicing acting as a bridge in evolution.

Alternative splicing (AS) regulates diverse cellular and developmental functions through alternative protein structures of different isoforms. Alternative exons dominate AS in vertebrates; however, very little is known about the extent and function of AS in lower eukaryotes. To understand the role of introns in gene evolution, we examined AS from a green algal and five fungal genomes using a novel EST-based gene-modeling algorithm (COMBEST). AS from each genome was classified with COMBEST that maps EST sequences to genomes to build gene models. Various aspects of AS were analyzed through statistical methods. The interplay of intron 3n length, phase, coding property, and intron retention (RI) were examined with Chi-square testing. With 3 to 834 times EST coverage, we identified up to 73% of AS in intron-containing genes and found preponderance of RI among 11 types of AS. The number of exons, expression level, and maximum intron length correlated with number of AS per gene (NAG), and intron-rich genes suppressed AS. Genes with AS were more ancient, and AS was conserved among fungal genomes. Among stopless introns, non-retained introns (NRI) avoided, but major RI preferred 3n length. In contrast, stop-containing introns showed uniform distribution among 3n, 3n+1, and 3n+2 lengths. We found a clue to the intron phase enigma: it was the coding function of introns involved in AS that dictates the intron phase bias. Majority of AS is non-functional, and the extent of AS is suppressed for intron-rich genes. RI through 3n length, stop codon, and phase bias bridges the transition from functionless to functional alternative isoforms.

Intron gain, a dominant evolutionary process supporting high levels of gene expression in rice

Presence of introns in eukaryotic genes, their evolution and biological function has been a subject of considerable debate ever since their discovery in 1977. To understand the effect of number of introns on the structural and functional characteristics of rice genes, we carried out whole genome analysis of the relationship of the number of introns per gene with predicted cDNA sequence (CDS) length, average exon length and gene expression patterns. There was a direct correlation between the number of introns and the average CDS length among the expressed rice genes, as determined by expressed sequence tags (EST) support. The percentage of expressed genes in groups of rice genes representing different intron numbers showed a significant positive correlation with the number of introns providing evidence for higher level of expression for intron-rich genes. This was further supported by higher abundance of ESTs for the intron-rich genes in the rice EST database. Higher number of introns may be providing post-transcriptional stability to the mRNA leading to higher expression levels. Here we first report the detailed genome wide analysis of distribution pattern of introns in rice that provides important insight in to understanding the evolution, structure and expression of genes in plant species. Particularly, the complex gene structure and functional advantage of the intron containing genes supports the gain of intron theory for the evolution of eukaryotic genes.

Reverse transcriptase and intron number evolution.

Introns are universal in eukaryotic genomes and play important roles in transcriptional regulation, mRNA export to the cytoplasm, nonsense-mediated decay as both a regulatory and a splicing quality control mechanism, R-loop avoidance, alternative splicing, chromatin structure, and evolution by exon-shuffling. Sixteen complete fungal genomes were used 13 of which were sequenced and annotated by JGI. Ustilago maydis, Cryptococcus neoformans, and Coprinus cinereus (also named Coprinopsis cinerea) were from the Broad Institute. Gene models from JGI-annotated genomes were taken from the GeneCatalog track that contained the best representative gene models. Varying fractions of the GeneCatalog were manually curated by external users. For clarity, we used the JGI unique database identifier. The last common ancestor of eukaryotes (LECA) has an estimated 6.4 coding exons per gene (EPG) and evolved into the diverse eukaryotic life forms, which is recapitulated by the development of a stem cell. We found a parallel between the simulated reverse transcriptase (RT)-mediated intron loss and the comparative analysis of 16 fungal genomes that spanned a wide range of intron density. Although footprints of RT (RTF) were dynamic, relative intron location (RIL) to the 5'-end of mRNA faithfully traced RT-mediated intron loss and revealed 7.7 EPG for LECA. The mode of exon length distribution was conserved in simulated intron loss, which was exemplified by the shared mode of 75 nt between fungal and Chlamydomonas genomes. The dominant ancient exon length was corroborated by the average exon length of the most intron-rich genes in fungal genomes and consistent with ancient protein modules being ~25 aa. Combined with the conservation of a protein length of 400 aa, the earliest ancestor of eukaryotes could have 16 EPG. During earlier evolution, Ascomycota's ancestor had significantly more 3'-biased RT-mediated intron loss that was followed by dramatic RTF loss. There was a down trend of EPG from more conserved to less conserved genes. Moreover, species-specific genes have higher exon-densities, shorter exons, and longer introns when compared to genes conserved at the phylum level. However, intron length in species-specific genes became shorter than that of genes conserved in all species after genomes experiencing drastic intron loss. The estimated EPG from the most frequent exon length is more than double that from the RIL method. This implies significant intron loss during the very early period of eukaryotic evolution. De novo gene-birth contributes to shorter exons, longer introns, and higher exon-density in species-specific genes relative to conserved genes.

Genome of Acanthamoeba castellanii highlights extensive lateral gene transfer and early evolution of tyrosine kinase signaling

BackgroundThe Amoebozoa constitute one of the primary divisions of eukaryotes, encompassing taxa of both biomedical and evolutionary importance, yet its genomic diversity remains largely unsampled. Here we present an analysis of a whole genome assembly of Acanthamoeba castellanii (Ac) the first representative from a solitary free-living amoebozoan.ResultsAc encodes 15,455 compact intron-rich genes, a significant number of which are predicted to have arisen through inter-kingdom lateral gene transfer (LGT). A majority of the LGT candidates have undergone a substantial degree of intronization and Ac appears to have incorporated them into established transcriptional programs. Ac manifests a complex signaling and cell communication repertoire, including a complete tyrosine kinase signaling toolkit and a comparable diversity of predicted extracellular receptors to that found in the facultatively multicellular dictyostelids. An important environmental host of a diverse range of bacteria and viruses, Ac utilizes a diverse repertoire of predicted pattern recognition receptors, many with predicted orthologous functions in the innate immune systems of higher organisms.ConclusionsOur analysis highlights the important role of LGT in the biology of Ac and in the diversification of microbial eukaryotes. The early evolution of a key signaling facility implicated in the evolution of metazoan multicellularity strongly argues for its emergence early in the Unikont lineage. Overall, the availability of an Ac genome should aid in deciphering the biology of the Amoebozoa and facilitate functional genomic studies in this important model organism and environmental host.

Identification and Functional Analysis of the erh1+ Gene Encoding Enhancer of Rudimentary Homolog from the Fission Yeast Schizosaccharomyces pombe

The ERH gene encodes a highly conserved small nuclear protein with a unique amino acid sequence and three-dimensional structure but unknown function. The gene is present in animals, plants, and protists but to date has only been found in few fungi. Here we report that ERH homologs are also present in all four species from the genus Schizosaccharomyces, S. pombe, S. octosporus, S. cryophilus, and S. japonicus, which, however, are an exception in this respect among Ascomycota and Basidiomycota. The ERH protein sequence is moderately conserved within the genus (58% identity between S. pombe and S. japonicus), but the intron-rich genes have almost identical intron-exon organizations in all four species. In S. pombe, erh1+ is expressed at a roughly constant level during vegetative growth and adaptation to unfavorable conditions such as nutrient limitation and hyperosmotic stress caused by sorbitol. Erh1p localizes preferentially to the nucleus with the exception of the nucleolus, but is also present in the cytoplasm. Cells lacking erh1+ have an aberrant cell morphology and a comma-like shape when cultured to the stationary phase, and exhibit a delayed recovery from this phase followed by slower growth. Loss of erh1+ in an auxotrophic background results in enhanced arrest in the G1 phase following nutritional stress, and also leads to hypersensitivity to agents inducing hyperosmotic stress (sorbitol), inhibiting DNA replication (hydroxyurea), and destabilizing the plasma membrane (SDS); this hypersensitivity can be abolished by expression of S. pombe erh1+ and, to a lesser extent, S. japonicus erh1+ or human ERH. Erh1p fails to interact with the human Ciz1 and PDIP46/SKAR proteins, known molecular partners of human ERH. Our data suggest that in Schizosaccharomyces sp. erh1+ is non-essential for normal growth and Erh1p could play a role in response to adverse environmental conditions and in cell cycle regulation.

Origin and evolution of spliceosomal introns

Evolution of exon-intron structure of eukaryotic genes has been a matter of long-standing, intensive debate. The introns-early concept, later rebranded ‘introns first’ held that protein-coding genes were interrupted by numerous introns even at the earliest stages of life's evolution and that introns played a major role in the origin of proteins by facilitating recombination of sequences coding for small protein/peptide modules. The introns-late concept held that introns emerged only in eukaryotes and new introns have been accumulating continuously throughout eukaryotic evolution. Analysis of orthologous genes from completely sequenced eukaryotic genomes revealed numerous shared intron positions in orthologous genes from animals and plants and even between animals, plants and protists, suggesting that many ancestral introns have persisted since the last eukaryotic common ancestor (LECA). Reconstructions of intron gain and loss using the growing collection of genomes of diverse eukaryotes and increasingly advanced probabilistic models convincingly show that the LECA and the ancestors of each eukaryotic supergroup had intron-rich genes, with intron densities comparable to those in the most intron-rich modern genomes such as those of vertebrates. The subsequent evolution in most lineages of eukaryotes involved primarily loss of introns, with only a few episodes of substantial intron gain that might have accompanied major evolutionary innovations such as the origin of metazoa. The original invasion of self-splicing Group II introns, presumably originating from the mitochondrial endosymbiont, into the genome of the emerging eukaryote might have been a key factor of eukaryogenesis that in particular triggered the origin of endomembranes and the nucleus. Conversely, splicing errors gave rise to alternative splicing, a major contribution to the biological complexity of multicellular eukaryotes. There is no indication that any prokaryote has ever possessed a spliceosome or introns in protein-coding genes, other than relatively rare mobile self-splicing introns. Thus, the introns-first scenario is not supported by any evidence but exon-intron structure of protein-coding genes appears to have evolved concomitantly with the eukaryotic cell, and introns were a major factor of evolution throughout the history of eukaryotes. This article was reviewed by I. King Jordan, Manuel Irimia (nominated by Anthony Poole), Tobias Mourier (nominated by Anthony Poole), and Fyodor Kondrashov. For the complete reports, see the Reviewers’ Reports section.

A Detailed History of Intron-rich Eukaryotic Ancestors Inferred from a Global Survey of 100 Complete Genomes

Protein-coding genes in eukaryotes are interrupted by introns, but intron densities widely differ between eukaryotic lineages. Vertebrates, some invertebrates and green plants have intron-rich genes, with 6–7 introns per kilobase of coding sequence, whereas most of the other eukaryotes have intron-poor genes. We reconstructed the history of intron gain and loss using a probabilistic Markov model (Markov Chain Monte Carlo, MCMC) on 245 orthologous genes from 99 genomes representing the three of the five supergroups of eukaryotes for which multiple genome sequences are available. Intron-rich ancestors are confidently reconstructed for each major group, with 53 to 74% of the human intron density inferred with 95% confidence for the Last Eukaryotic Common Ancestor (LECA). The results of the MCMC reconstruction are compared with the reconstructions obtained using Maximum Likelihood (ML) and Dollo parsimony methods. An excellent agreement between the MCMC and ML inferences is demonstrated whereas Dollo parsimony introduces a noticeable bias in the estimations, typically yielding lower ancestral intron densities than MCMC and ML. Evolution of eukaryotic genes was dominated by intron loss, with substantial gain only at the bases of several major branches including plants and animals. The highest intron density, 120 to 130% of the human value, is inferred for the last common ancestor of animals. The reconstruction shows that the entire line of descent from LECA to mammals was intron-rich, a state conducive to the evolution of alternative splicing.

Characterization of Maltase Clusters in the Genus Drosophila

To reveal evolutionary history of maltase gene family in the genus Drosophila, we undertook a bioinformatics study of maltase genes from available genomes of 12 Drosophila species. Molecular evolution of a closely related glycoside hydrolase, the α-amylase, in Drosophila has been extensively studied for a long time. The α-amylases were even used as a model of evolution of multigene families. On the other hand, maltase, i.e., the α-glucosidase, got only scarce attention. In this study, we, therefore, investigated spatial organization of the maltase genes in Drosophila genomes, compared the amino acid sequences of the encoded enzymes and analyzed the intron/exon composition of orthologous genes. We found that the Drosophila maltases are more numerous than previously thought (ten instead of three genes) and are localized in two clusters on two chromosomes (2L and 2R). To elucidate the approximate time line of evolution of the clusters, we estimated the order and dated duplication of all the 10 genes. Both clusters are the result of ancient series of subsequent duplication events, which took place from 352 to 61million years ago, i.e., well before speciation to extant Drosophila species. Also observed was a remarkable intron/exon composition diversity of particular maltase genes of these clusters, probably a result of independent intron loss after duplication of intron-rich gene ancestor, which emerged well before speciation in a common ancestor of all extant Drosophila species.

A general framework for optimization of probes for gene expression microarray and its application to the fungus Podospora anserina

BackgroundThe development of new microarray technologies makes custom long oligonucleotide arrays affordable for many experimental applications, notably gene expression analyses. Reliable results depend on probe design quality and selection. Probe design strategy should cope with the limited accuracy of de novo gene prediction programs, and annotation up-dating. We present a novel in silico procedure which addresses these issues and includes experimental screening, as an empirical approach is the best strategy to identify optimal probes in the in silico outcome.FindingsWe used four criteria for in silico probe selection: cross-hybridization, hairpin stability, probe location relative to coding sequence end and intron position. This latter criterion is critical when exon-intron gene structure predictions for intron-rich genes are inaccurate. For each coding sequence (CDS), we selected a sub-set of four probes. These probes were included in a test microarray, which was used to evaluate the hybridization behavior of each probe. The best probe for each CDS was selected according to three experimental criteria: signal-to-noise ratio, signal reproducibility, and representative signal intensities. This procedure was applied for the development of a gene expression Agilent platform for the filamentous fungus Podospora anserina and the selection of a single 60-mer probe for each of the 10,556 P. anserina CDS.ConclusionsA reliable gene expression microarray version based on the Agilent 44K platform was developed with four spot replicates of each probe to increase statistical significance of analysis.

Molecular characterization of beta-tubulin from Phakopsora pachyrhizi, the causal agent of Asian soybean rust

β-tubulins are structural components of microtubules and the targets of benzimidazole fungicides used to control many diseases of agricultural importance. Intron polymorphisms in the intron-rich genes of these proteins have been used in phylogeographic investigations of phytopathogenic fungi. In this work, we sequenced 2764 nucleotides of the β-tubulin gene (Pp tubB) in samples of Phakopsora pachyrhizi collected from seven soybean fields in Brazil. Pp tubB contained an open reading frame of 1341 nucleotides, including nine exons and eight introns. Exon length varied from 14 to 880 nucleotides, whereas intron length varied from 76 to 102 nucleotides. The presence of only four polymorphic sites limited the usefulness of Pp tubB for phylogeographic studies in P. pachyrhizi. The gene structures of Pp tubB and orthologous β-tubulin genes of Melampsora lini and Uromyces viciae-fabae were highly conserved. The amino acid substitutions in β-tubulin proteins associated with the onset of benzimidazole resistance in model organisms, especially at His 6 , Glu 198 and Phe 200 , were absent from the predicted sequence of the P. pachyrhizi β-tubulin protein.

Origination of the Split Structure of Spliceosomal Genes from Random Genetic Sequences

The mechanism by which protein-coding portions of eukaryotic genes came to be separated by long non-coding stretches of DNA, and the purpose for this perplexing arrangement, have remained unresolved fundamental biological problems for three decades. We report here a plausible solution to this problem based on analysis of open reading frame (ORF) length constraints in the genomes of nine diverse species. If primordial nucleic acid sequences were random in sequence, functional proteins that are innately long would not be encoded due to the frequent occurrence of stop codons. The best possible way that a long protein-coding sequence could have been derived was by evolving a split-structure from the random DNA (or RNA) sequence. Results of the systematic analyses of nine complete genome sequences presented here suggests that perhaps the major underlying structural features of split-genes have evolved due to the indigenous occurrence of split protein-coding genes in primordial random nucleotide sequence. The results also suggest that intron-rich genes containing short exons may have been the original form of genes intrinsically occurring in random DNA, and that intron-poor genes containing long exons were perhaps derived from the original intron-rich genes.

From Stop to Start: Tandem Gene Arrangement, Copy Number and Trans-Splicing Sites in the Dinoflagellate Amphidinium carterae

Dinoflagellate genomes present unique challenges including large size, modified DNA bases, lack of nucleosomes, and condensed chromosomes. EST sequencing has shown that many genes are found as many slightly different variants implying that many copies are present in the genome. As a preliminary survey of the genome our goal was to obtain genomic sequences for 47 genes from the dinoflagellate Amphidinium carterae. A PCR approach was used to avoid problems with large insert libraries. One primer set was oriented inward to amplify the genomic complement of the cDNA and a second primer set would amplify outward between tandem repeats of the same gene. Each gene was also tested for a spliced leader using cDNA as template. Almost all (14/15) of the highly expressed genes (i.e. those with high representation in the cDNA pool) were shown to be in tandem arrays with short intergenic spacers, and most were trans-spliced. Only two moderately expressed genes were found in tandem arrays. A polyadenylation signal was found in genomic copies containing the sequence AAAAG/C at the exact polyadenylation site and was conserved between species. Four genes were found to have a high intron density (>5 introns) while most either lacked introns, or had only one to three. Actin was selected for deeper sequencing of both genomic and cDNA copies. Two clusters of actin copies were found, separated from each other by many non-coding features such as intron size and sequence. One intron-rich gene was selected for genomic walking using inverse PCR, and was not shown to be in a tandem repeat. The first glimpse of dinoflagellate genome indicates two general categories of genes in dinoflagellates, a highly expressed tandem repeat class and an intron rich less expressed class. This combination of features appears to be unique among eukaryotes.

Extremely Intron-Rich Genes in the Alveolate Ancestors Inferred with a Flexible Maximum-Likelihood Approach

Chromalveolates are a large, diverse supergroup of unicellular eukaryotes that includes Apicomplexa, dinoflagellates, ciliates (three lineages that form the alveolate branch), heterokonts, haptophytes, and cryptomonads (three lineages comprising the chromist branch). All sequenced genomes of chromalveolates have relatively low intron density in protein-coding genes, and few intron positions are shared between chromalveolate lineages. In contrast, genes of different chromalveolates share many intron positions with orthologous genes from other eukaryotic supergroups, in particular, the intron-rich orthologs from animals and plants. Reconstruction of the history of intron gain and loss during the evolution of chromalveolates using a general and flexible maximum-likelihood approach indicates that genes of the ancestors of chromalveolates and, particularly, alveolates had unexpectedly high intron densities. It is estimated that the chromalveolate ancestor had, approximately, two-third of the human intron density, whereas the intron density in the genes of the alveolate ancestor is estimated to be slightly greater than the human intron density. Accordingly, it is inferred that the evolution of chromalveolates was dominated by intron loss. The conclusion that ancestral chromalveolate forms had high intron densities is unexpected because all extant unicellular eukaryotes have relatively few introns and are thought to be unable to maintain numerous introns due to intense purifying selection in their, typically, large populations. It is suggested that, at early stages of evolution, chromalveolates went through major population bottlenecks that were accompanied by intron invasion.