Evgeny S. Gerasimov

Genome of Leptomonas pyrrhocoris: a high-quality reference for monoxenous trypanosomatids and new insights into evolution of Leishmania.

Many high-quality genomes are available for dixenous (two hosts) trypanosomatid species of the genera Trypanosoma, Leishmania, and Phytomonas, but only fragmentary information is available for monoxenous (single-host) trypanosomatids. In trypanosomatids, monoxeny is ancestral to dixeny, thus it is anticipated that the genome sequences of the key monoxenous parasites will be instrumental for both understanding the origin of parasitism and the evolution of dixeny. Here, we present a high-quality genome for Leptomonas pyrrhocoris, which is closely related to the dixenous genus Leishmania. The L. pyrrhocoris genome (30.4 Mbp in 60 scaffolds) encodes 10,148 genes. Using the L. pyrrhocoris genome, we pinpointed genes gained in Leishmania. Among those genes, 20 genes with unknown function had expression patterns in the Leishmania mexicana life cycle suggesting their involvement in virulence. By combining differential expression data for L. mexicana, L. major and Leptomonas seymouri, we have identified several additional proteins potentially involved in virulence, including SpoU methylase and U3 small nucleolar ribonucleoprotein IMP3. The population genetics of L. pyrrhocoris was also addressed by sequencing thirteen strains of different geographic origin, allowing the identification of 1,318 genes under positive selection. This set of genes was significantly enriched in components of the cytoskeleton and the flagellum.

Diversity of mitochondrial genome organization

In this review, we discuss types of mitochondrial genome structural organization (architecture), which includes the following characteristic features: size and the shape of DNA molecule, number of encoded genes, presence of cryptogenes, and editing of primary transcripts.

Gene Loss and Error-Prone RNA Editing in the Mitochondrion of Perkinsela, an Endosymbiotic Kinetoplastid

ABSTRACT Perkinsela is an enigmatic early-branching kinetoplastid protist that lives as an obligate endosymbiont inside Paramoeba (Amoebozoa). We have sequenced the highly reduced mitochondrial genome of Perkinsela, which possesses only six protein-coding genes (cox1, cox2, cox3, cob, atp6, and rps12), despite the fact that the organelle itself contains more DNA than is present in either the host or endosymbiont nuclear genomes. An in silico analysis of two Perkinsela strains showed that mitochondrial RNA editing and processing machineries typical of kinetoplastid flagellates are generally conserved, and all mitochondrial transcripts undergo U-insertion/deletion editing. Canonical kinetoplastid mitochondrial ribosomes are also present. We have developed software tools for accurate and exhaustive mapping of transcriptome sequencing (RNA-seq) reads with extensive U-insertions/deletions, which allows detailed investigation of RNA editing via deep sequencing. With these methods, we show that up to 50% of reads for a given edited region contain errors of the editing system or, less likely, correspond to alternatively edited transcripts. IMPORTANCE Uridine insertion/deletion-type RNA editing, which occurs in the mitochondrion of kinetoplastid protists, has been well-studied in the model parasite genera Trypanosoma, Leishmania, and Crithidia. Perkinsela provides a unique opportunity to broaden our knowledge of RNA editing machinery from an evolutionary perspective, as it represents the earliest kinetoplastid branch and is an obligatory endosymbiont with extensive reductive trends. Interestingly, up to 50% of mitochondrial transcripts in Perkinsela contain errors. Our study was complemented by use of newly developed software designed for accurate mapping of extensively edited RNA-seq reads obtained by deep sequencing. Uridine insertion/deletion-type RNA editing, which occurs in the mitochondrion of kinetoplastid protists, has been well-studied in the model parasite genera Trypanosoma, Leishmania, and Crithidia. Perkinsela provides a unique opportunity to broaden our knowledge of RNA editing machinery from an evolutionary perspective, as it represents the earliest kinetoplastid branch and is an obligatory endosymbiont with extensive reductive trends. Interestingly, up to 50% of mitochondrial transcripts in Perkinsela contain errors. Our study was complemented by use of newly developed software designed for accurate mapping of extensively edited RNA-seq reads obtained by deep sequencing.

Genomic study of the Ket: a Paleo-Eskimo-related ethnic group with significant ancient North Eurasian ancestry

The Kets, an ethnic group in the Yenisei River basin, Russia, are considered the last nomadic hunter-gatherers of Siberia, and Ket language has no transparent affiliation with any language family. We investigated connections between the Kets and Siberian and North American populations, with emphasis on the Mal’ta and Paleo-Eskimo ancient genomes, using original data from 46 unrelated samples of Kets and 42 samples of their neighboring ethnic groups (Uralic-speaking Nganasans, Enets, and Selkups). We genotyped over 130,000 autosomal SNPs, identified mitochondrial and Y-chromosomal haplogroups, and performed high-coverage genome sequencing of two Ket individuals. We established that Nganasans, Kets, Selkups, and Yukaghirs form a cluster of populations most closely related to Paleo-Eskimos in Siberia (not considering indigenous populations of Chukotka and Kamchatka). Kets are closely related to modern Selkups and to some Bronze and Iron Age populations of the Altai region, with all these groups sharing a high degree of Mal’ta ancestry. Implications of these findings for the linguistic hypothesis uniting Ket and Na-Dene languages into a language macrofamily are discussed.

Trypanosomatid mitochondrial RNA editing: Dramatically complex transcript repertoires revealed with a dedicated mapping tool

Abstract RNA editing by targeted insertion and deletion of uridine is crucial to generate translatable mRNAs from the cryptogenes of the mitochondrial genome of kinetoplastids. This type of editing consists of a stepwise cascade of reactions generally proceeding from 3′ to 5′ on a transcript, resulting in a population of partially edited as well as pre-edited and completely edited molecules for each mitochondrial cryptogene of these protozoans. Often, the number of uridines inserted and deleted exceed the number of nucleotides that are genome-encoded. Thus, analysis of kinetoplastid mitochondrial transcriptomes has proven frustratingly complex. Here we present our analysis of Leptomonas pyrrhocoris mitochondrial cDNA deep sequencing reads using T-Aligner, our new tool which allows comprehensive characterization of RNA editing, not relying on targeted transcript amplification and on prior knowledge of final edited products. T-Aligner implements a pipeline of read mapping, visualization of all editing states and their coverage, and assembly of canonical and alternative translatable mRNAs. We also assess T-Aligner functionality on a more challenging deep sequencing read input from Trypanosoma cruzi. The analysis reveals that transcripts of cryptogenes of both species undergo very complex editing that includes the formation of alternative open reading frames and whole categories of truncated editing products.

Selective amplification of maxicircle classes during the life cycle of Leishmania major.

The kinetoplast genome contains several thousands of minicircles of various sequence classes and several scores of maxicircles. We demonstrated that maxicircles are heterogeneous in clonal cultures of Leishmania major, and, therefore, probably heterogeneous (heteroplasmic) within the kinetoplast. Sequence heterogeneity was observed in a non-coding fragment upstream of the 12S rRNA gene. We identified about 20 stable variants of this fragment, which were composed of one to five non-identical repeats 200-300bp in length. Promastigote-to-amastigote and amastigote-to-promastigote differentiation was often accompanied by shifts in abundance of some maxicircle classes. Reversion to promastigote-specific maxicircle patterns was usually observed in the life cycle (promastigote-amastigote-promastigote), however there were many exceptions.

Evolution of the Genome 3D Organization: Comparison of Fused and Segregated Globin Gene Clusters

The genomes are folded in a complex three-dimensional (3D) structure. Some features of this organization are common for all eukaryotes, but little is known about its evolution. Here, we have studied the 3D organization and regulation of zebrafish globin gene domain and compared its organization and regulation with those of other vertebrate species. In birds and mammals, the α- and β-globin genes are segregated into separate clusters located on different chromosomes and organized into chromatin domains of different types, whereas in cold-blooded vertebrates, including Danio rerio, α- and β-globin genes are organized into common clusters. The major globin gene locus of Danio rerio is of particular interest as it is located in a genomic area that is syntenic in vertebrates and is controlled by a conserved enhancer. We have found that the major globin gene locus of Danio rerio is structurally and functionally segregated into two spatially distinct subloci harboring either adult or embryo-larval globin genes. These subloci demonstrate different organization at the level of chromatin domains and different modes of spatial organization, which appears to be due to selective interaction of the upstream enhancer with the sublocus harboring globin genes of the adult type. These data are discussed in terms of evolution of linear and 3D organization of gene clusters in vertebrates.

Comparative genome analysis of Pseudogymnoascus spp. reveals primarily clonal evolution with small genome fragments exchanged between lineages

BackgroundPseudogymnoascus spp. is a wide group of fungi lineages in the family Pseudorotiaceae including an aggressive pathogen of bats P. destructans. Although several lineages of P. spp. were shown to produce ascospores in culture, the vast majority of P. spp. demonstrates no evidence of sexual reproduction. P. spp. can tolerate a wide range of different temperatures and salinities and can survive even in permafrost layer. Adaptability of P. spp. to different environments is accompanied by extremely variable morphology and physiology.ResultsWe sequenced genotypes of 14 strains of P. spp., 5 of which were extracted from permafrost, 1 from a cryopeg, a layer of unfrozen ground in permafrost, and 8 from temperate surface environments. All sequenced genotypes are haploid. Nucleotide diversity among these genomes is very high, with a typical evolutionary distance at synonymous sites dS ≈ 0.5, suggesting that the last common ancestor of these strains lived >50Mya. The strains extracted from permafrost do not form a separate clade. Instead, each permafrost strain has close relatives from temperate environments.We observed a strictly clonal population structure with no conflicting topologies for ~99% of genome sequences. However, there is a number of short (~100–10,000 nt) genomic segments with the total length of 67.6 Kb which possess phylogenetic patterns strikingly different from the rest of the genome. The most remarkable case is a MAT-locus, which has 2 distinct alleles interspersed along the whole-genome phylogenetic tree.ConclusionsPredominantly clonal structure of genome sequences is consistent with the observations that sexual reproduction is rare in P. spp. Small number of regions with noncanonical phylogenies seem to arise due to some recombination events between derived lineages of P. spp., with MAT-locus being transferred on multiple occasions. All sequenced strains have heterothallic configuration of MAT-locus.

Minicircle kinetoplast genome of insect trypanosomatid Leptomonas pyrrhocoris

We present here the structure of a minicircle population based on transcriptome sequencing of Leptomonas pyrrhocoris. We show that minicircle DNA molecules are dimeric. As in dixenous species, the entire molecule of minicircle DNA is transcribed. This is the first minicircle transcriptome of monoxenous trypanosomatid species determined using NGS technology.

Gene expression to mitochondrial metabolism: Variability among cultured Trypanosoma cruzi strains

The insect-transmitted protozoan parasite Trypanosoma cruzi experiences changes in nutrient availability and rate of flux through different metabolic pathways across its life cycle. The species encompasses much genetic diversity of both the nuclear and mitochondrial genomes among isolated strains. The genetic or expression variation of both genomes are likely to impact metabolic responses to environmental stimuli, and even steady state metabolic function, among strains. To begin formal characterization these differences, we compared aspects of metabolism between genetically similar strains CL Brener and Tulahuen with less similar Esmeraldo and Sylvio X10 strains in a culture environment. Epimastigotes of all strains took up glucose at similar rates. However, the degree of medium acidification that could be observed when glucose was absent from the medium varied by strain, indicating potential differences in excreted metabolic byproducts. Our main focus was differences related to electron transport chain function. We observed differences in ATP-coupled respiration and maximal respiratory capacity, mitochondrial membrane potential, and mitochondrial morphology between strains, despite the fact that abundances of two nuclear-encoded proteins of the electron transport chain are similar between strains. RNA sequencing reveals strain-specific differences in abundances of mRNAs encoding proteins of the respiratory chain but also other metabolic processes. From these differences in metabolism and mitochondrial phenotypes we have generated tentative models for the differential metabolic fluxes or differences in gene expression that may underlie these results.