Daniel C. Jeffares

Comparative genomic analysis of three Leishmania species that cause diverse human disease.

Leishmania parasites cause a broad spectrum of clinical disease. Here we report the sequencing of the genomes of two species of Leishmania: Leishmania infantum and Leishmania braziliensis. The comparison of these sequences with the published genome of Leishmania major reveals marked conservation of synteny and identifies only ∼200 genes with a differential distribution between the three species. L. braziliensis, contrary to Leishmania species examined so far, possesses components of a putative RNA-mediated interference pathway, telomere-associated transposable elements and spliced leader–associated SLACS retrotransposons. We show that pseudogene formation and gene loss are the principal forces shaping the different genomes. Genes that are differentially distributed between the species encode proteins implicated in host-pathogen interactions and parasite survival in the macrophage.

Genome variation and evolution of the malaria parasite Plasmodium falciparum

Infections with the malaria parasite Plasmodium falciparum result in more than 1 million deaths each year worldwide. Deciphering the evolutionary history and genetic variation of P. falciparum is critical for understanding the evolution of drug resistance, identifying potential vaccine candidates and appreciating the effect of parasite variation on prevalence and severity of malaria in humans. Most studies of natural variation in P. falciparum have been either in depth over small genomic regions (up to the size of a small chromosome) or genome wide but only at low resolution. In an effort to complement these studies with genome-wide data, we undertook shotgun sequencing of a Ghanaian clinical isolate (with fivefold coverage), the IT laboratory isolate (with onefold coverage) and the chimpanzee parasite P. reichenowi (with twofold coverage). We compared these sequences with the fully sequenced P. falciparum 3D7 isolate genome. We describe the most salient features of P. falciparum polymorphism and adaptive evolution with relation to gene function, transcript and protein expression and cellular localization. This analysis uncovers the primary evolutionary changes that have occurred since the P. falciparum–P. reichenowi speciation and changes that are occurring within P. falciparum.NOTE: In the original version of this paper, the authors failed to acknowledge that sequencing of the P. falciparum IT laboratory isolate was funded by a European Union 6th Framework Program grant to the BioMalPar Consortium (contract number LSHP-LT-2004-503578). This error has been corrected in the PDF version of the article.

Rapidly regulated genes are intron poor

We show that genes with rapidly changing expression levels in response to stress contain significantly lower intron densities in yeasts, thale cress and mice. Therefore, we propose that introns can delay regulatory responses and are selected against in genes whose transcripts require rapid adjustment for survival of environmental challenges. These findings could provide an explanation for the apparent extensive intron loss during the evolution of some eukaryotic lineages.

EARLY EVOLUTION: PROKARYOTES, THE NEW KIDS ON THE BLOCK

Prokaryotes are generally assumed to be the oldest existing form of life on earth. This assumption, however, makes it difficult to understand certain aspects of the transition from earlier stages in the origin of life to more complex ones, and it does not account for many apparently ancient features in the eukaryotes. From a model of the RNA world, based on relic RNA species in modern organisms, one can infer that there was an absolute requirement for a high-accuracy RNA replicase even before proteins evolved. In addition, we argue here that the ribosome (together with the RNAs involved in its assembly) is so large that it must have had a prior function before protein synthesis. A model that connects and equates these two requirements (high-accuracy RNA replicase and prior function of the ribosome) can explain many steps in the origin of life while accounting for the observation that eukaryotes have retained more vestiges of the RNA world. The later derivation of prokaryote RNA metabolism and genome structure can be accounted for by the two complementary mechanisms of r-selection and thermoreduction.

The genomic and phenotypic diversity of Schizosaccharomyces pombe

Natural variation within species reveals aspects of genome evolution and function. The fission yeast Schizosaccharomyces pombe is an important model for eukaryotic biology, but researchers typically use one standard laboratory strain. To extend the usefulness of this model, we surveyed the genomic and phenotypic variation in 161 natural isolates. We sequenced the genomes of all strains, finding moderate genetic diversity (π = 3 × 10−3 substitutions/site) and weak global population structure. We estimate that dispersal of S. pombe began during human antiquity (∼340 BCE), and ancestors of these strains reached the Americas at ∼1623 CE. We quantified 74 traits, finding substantial heritable phenotypic diversity. We conducted 223 genome-wide association studies, with 89 traits showing at least one association. The most significant variant for each trait explained 22% of the phenotypic variance on average, with indels having larger effects than SNPs. This analysis represents a rich resource to examine genotype-phenotype relationships in a tractable model.

An Overview of the Introns-First Theory

We review the introns-first hypothesis a decade after it was first proposed. It is that exons emerged from non-coding regions interspersed between RNA genes in an early RNA world, and is a subcomponent of a more general ‘RNA-continuity’ hypothesis. The latter is that some RNA-based systems, especially in RNA processing, are ‘relics’ that can be traced back either to the RNA world that preceded both DNA and encoded protein synthesis or to the later ribonucleoprotein (RNP) world (before DNA took over the main coding role). RNA-continuity is based on independent evidence—in particular, the relative inefficiency of RNA catalysis compared with protein catalysis—and leads to a wide range of predictions, ranging from the origin of the ribosome, the spliceosome, small nucleolar RNAs, RNases P and MRP, and mRNA, and it is consistent with the wide involvement of RNA-processing and regulation of RNA in modern eukaryotes. While there may still be cause to withhold judgement on intron origins, there is strong evidence against introns being uncommon in the last eukaryotic common ancestor (LECA), and expanding only within extant eukaryotic groups—the ‘very-late’ intron invasion model. Similarly, it is clear that there are selective forces on numbers and positions of introns; their existence may not always be neutral. There is still a range of viable alternatives, including introns first, early, and ‘latish’ (i.e. well established in LECA), and regardless of which is ultimately correct, it pays to separate out various questions and to focus on testing the predictions of sub-theories.

Evolutionarily Stable Association of Intronic snoRNAs and microRNAs with Their Host Genes

Small nucleolar RNAs (snoRNAs) and microRNAs (miRNAs) are integral to a range of processes, including ribosome biogenesis and gene regulation. Some are intron encoded, and this organization may facilitate coordinated coexpression of host gene and RNA. However, snoRNAs and miRNAs are known to be mobile, so intron-RNA associations may not be evolutionarily stable. We have used genome alignments across 11 mammals plus chicken to examine positional orthology of snoRNAs and miRNAs and report that 21% of annotated snoRNAs and 11% of miRNAs are positionally conserved across mammals. Among RNAs traceable to the bird–mammal common ancestor, 98% of snoRNAs and 76% of miRNAs are intronic. Comparison of the most evolutionarily stable mammalian intronic snoRNAs with those positionally conserved among primates reveals that the former are more overrepresented among host genes involved in translation or ribosome biogenesis and are more broadly and highly expressed. This stability is likely attributable to a requirement for overlap between host gene and intronic snoRNA expression profiles, consistent with an ancestral role in ribosome biogenesis. In contrast, whereas miRNA positional conservation is comparable to that observed for snoRNAs, intronic miRNAs show no obvious association with host genes of a particular functional category, and no statistically significant differences in host gene expression are found between those traceable to mammalian or primate ancestors. Our results indicate evolutionarily stable associations of numerous intronic snoRNAs and miRNAs and their host genes, with probable continued diversification of snoRNA function from an ancestral role in ribosome biogenesis.

LaSSO, a strategy for genome-wide mapping of intronic lariats and branch points using RNA-seq

Both canonical and alternative splicing of RNAs are governed by intronic sequence elements and produce transient lariat structures fastened by branch points within introns. To map precisely the location of branch points on a genomic scale, we developed LaSSO (Lariat Sequence Site Origin), a data-driven algorithm which utilizes RNA-seq data. Using fission yeast cells lacking the debranching enzyme Dbr1, LaSSO not only accurately identified canonical splicing events, but also pinpointed novel, but rare, exon-skipping events, which may reflect aberrantly spliced transcripts. Compromised intron turnover perturbed gene regulation at multiple levels, including splicing and protein translation. Notably, Dbr1 function was also critical for the expression of mitochondrial genes and for the processing of self-spliced mitochondrial introns. LaSSO showed better sensitivity and accuracy than algorithms used for computational branch-point prediction or for empirical branch-point determination. Even when applied to a human data set acquired in the presence of debranching activity, LaSSO identified both canonical and exon-skipping branch points. LaSSO thus provides an effective approach for defining high-resolution maps of branch-site sequences and intronic elements on a genomic scale. LaSSO should be useful to validate introns and uncover branch-point sequences in any eukaryote, and it could be integrated into RNA-seq pipelines.

Selected Schizosaccharomyces pombe Strains Have Characteristics That Are Beneficial for Winemaking.

At present, wine is generally produced using Saccharomyces yeast followed by Oenococus bacteria to complete malolactic fermentation. This method has some unsolved problems, such as the management of highly acidic musts and the production of potentially toxic products including biogenic amines and ethyl carbamate. Here we explore the potential of the fission yeast Schizosaccharomyces pombe to solve these problems. We characterise an extensive worldwide collection of S. pombe strains according to classic biochemical parameters of oenological interest. We identify three genetically different S. pombe strains that appear suitable for winemaking. These strains compare favourably to standard Saccharomyces cerevisiae winemaking strains, in that they perform effective malic acid deacidification and significantly reduce levels of biogenic amines and ethyl carbamate precursors without the need for any secondary bacterial malolactic fermentation. These findings indicate that the use of certain S. pombe strains could be advantageous for winemaking in regions where malic acid is problematic, and these strains also show superior performance with respect to food safety.

Diversification of genes encoding mei2-like RNA binding proteins in plants

A predominantly plant-based family of genes encoding RNA binding proteins is defined by the presence of a highly conserved RNA binding motif first described in the mei2 gene of the fission yeast Schizosaccharomyces pombe. In silico analyses reveal nine mei2-like genes in Arabidopsis thaliana and six in Oryza sativa. These predicted genes group into four distinct clades, based on overall sequence similarity and subfamily-specific sequence elements. In situ analysis show that Arabidopsis genes from one of these clades, TEL1 and TEL2, are specifically expressed in central zone of the shoot apical meristem and the quiescent center of the root apical meristem, suggesting that they may somehow function to maintain indeterminacy in these tissues. By contrast, members of two sister clades, AML1 through AML5, are expressed more broadly, a trend that was confirmed by Q-PCR analysis. mei2-like transcripts with similar sequences showed similar expression patterns, suggesting functional redundancy within the four clades. Phenotypic analyses of lines that contain T-DNA insertions to individual mei2-like genes reveal no obvious phenotypes, further suggesting redundant activities for these gene products.