Blastocrithidia nonstop mitochondrial genome and its expression are remarkably insulated from nuclear codon reassignment.

The canonical stop codons of the nuclear genome of the trypanosomatid Blastocrithidia nonstop are recoded. Here, we investigated the effect of this recoding on the mitochondrial genome and gene expression. Trypanosomatids possess a single mitochondrion and protein-coding transcripts of this genome require RNA editing in order to generate open reading frames of many transcripts encoded as 'cryptogenes'. Small RNAs that can number in the hundreds direct editing and produce a mitochondrial transcriptome of unusual complexity. We find B. nonstop to have a typical trypanosomatid mitochondrial genetic code, which presumably requires the mitochondrion to disable utilization of the two nucleus-encoded suppressor tRNAs, which appear to be imported into the organelle. Alterations of the protein factors responsible for mRNA editing were also documented, but they have likely originated from sources other than B. nonstop nuclear genome recoding. The population of guide RNAs directing editing is minimal, yet virtually all genes for the plethora of known editing factors are still present. Most intriguingly, despite lacking complex I cryptogene guide RNAs, these cryptogene transcripts are stochastically edited to high levels.

Trypanosoma cruzi strain and starvation-driven mitochondrial RNA editing and transcriptome variability.

Trypanosoma cruzi is a unicellular protistan parasitic species that is comprised of strains and isolates exhibiting high levels of genetic and metabolic variability. In the insect vector, it is known to be highly responsive to starvation, a signal for progression to a life stage in which it can infect mammalian cells. Most mRNAs encoded in its mitochondrion require the targeted insertion and deletion of uridines to become translatable transcripts. This study defined differences in uridine-insertion/deletion RNA editing among three strains and established the mechanism whereby abundances of edited (and, thus, translatable) mitochondrial gene products increase during starvation. Our approach utilized our custom T-Aligner toolkit to describe transcriptome-wide editing events and reconstruct editing products from high-throughput sequencing data. We found that the relative abundance of mitochondrial transcripts and the proportion of mRNAs that are edited varies greatly between analyzed strains, a characteristic that could potentially impact metabolic capacity. Starvation typically led to an increase in overall editing activity rather than affecting a specific step in the process. We also determined that transcripts CR3, CR4, and ND3 produce multiple open reading frames that, if translated, would generate different proteins. Finally, we quantitated the inherent flexibility of editing in T. cruzi and found it to be higher relative to that in a related trypanosomatid lineage. Over time, new editing domains or patterns could prove advantageous to the organism and become more widespread within individual transcriptomes or among strains.

Shining the spotlight on the neglected: new high-quality genome assemblies as a gateway to understanding the evolution of Trypanosomatidae.

BACKGROUND: Protists of the family Trypanosomatidae (phylum Euglenozoa) have gained notoriety as parasites affecting humans, domestic animals, and agricultural plants. However, the true extent of the group's diversity spreads far beyond the medically and veterinary relevant species. We address several knowledge gaps in trypanosomatid research by undertaking sequencing, assembly, and analysis of genomes from previously overlooked representatives of this protistan group. RESULTS: We assembled genomes for twenty-one trypanosomatid species, with a primary focus on insect parasites and Trypanosoma spp. parasitizing non-human hosts. The assemblies exhibit sizes consistent with previously sequenced trypanosomatid genomes, ranging from approximately 18 Mb for Obscuromonas modryi to 35 Mb for Crithidia brevicula and Zelonia costaricensis. Despite being the smallest, the genome of O. modryi has the highest content of repetitive elements, contributing nearly half of its total size. Conversely, the highest proportion of unique DNA is found in the genomes of Wallacemonas spp., with repeats accounting for less than 8% of the assembly length. The majority of examined species exhibit varying degrees of aneuploidy, with trisomy being the most frequently observed condition after disomy. CONCLUSIONS: The genome of Obscuromonas modryi represents a very unusual, if not unique, example of evolution driven by two antidromous forces: i) increasing dependence on the host leading to genomic shrinkage and ii) expansion of repeats causing genome enlargement. The observed variation in somy within and between trypanosomatid genera suggests that these flagellates are largely predisposed to aneuploidy and, apparently, exploit it to gain a fitness advantage. High heterogeneity in the genome size, repeat content, and variation in chromosome copy numbers in the newly-sequenced species highlight the remarkable genome plasticity exhibited by trypanosomatid flagellates. These new genome assemblies are a robust foundation for future research on the genetic basis of life cycle changes and adaptation to different hosts in the family Trypanosomatidae.

Complete minicircle genome of Leptomonas pyrrhocoris reveals sources of its non-canonical mitochondrial RNA editing events.

Uridine insertion/deletion (U-indel) editing of mitochondrial mRNA, unique to the protistan class Kinetoplastea, generates canonical as well as potentially non-productive editing events. While the molecular machinery and the role of the guide (g) RNAs that provide required information for U-indel editing are well understood, little is known about the forces underlying its apparently error-prone nature. Analysis of a gRNA:mRNA pair allows the dissection of editing events in a given position of a given mitochondrial transcript. A complete gRNA dataset, paired with a fully characterized mRNA population that includes non-canonically edited transcripts, would allow such an analysis to be performed globally across the mitochondrial transcriptome. To achieve this, we have assembled 67 minicircles of the insect parasite Leptomonas pyrrhocoris, with each minicircle typically encoding one gRNA located in one of two similar-sized units of different origin. From this relatively narrow set of annotated gRNAs, we have dissected all identified mitochondrial editing events in L. pyrrhocoris, the strains of which dramatically differ in the abundance of individual minicircle classes. Our results support a model in which a multitude of editing events are driven by a limited set of gRNAs, with individual gRNAs possessing an inherent ability to guide canonical and non-canonical editing.

A link between mitochondrial gene expression and life stage morphologies in Trypanosoma cruzi.

The protozoan Trypanosoma cruzi has a complicated dual-host life cycle, and starvation can trigger transition from the replicating insect stage to the mammalian-infectious nonreplicating insect stage (epimastigote to trypomastigote differentiation). Abundance of some mature RNAs derived from its mitochondrial genome increase during culture starvation of T. cruzi for unknown reasons. Here, we examine T. cruzi mitochondrial gene expression in the mammalian intracellular replicating life stage (amastigote), and uncover implications of starvation-induced changes in gene expression. Mitochondrial RNA levels in general were found to be lowest in actively replicating amastigotes. We discovered that mitochondrial respiration decreases during starvation in insect stage cells, despite the previously observed increases in mitochondrial mRNAs encoding electron transport chain (ETC) components. Surprisingly, T. cruzi epimastigotes in replete medium grow at normal rates when we genetically compromised their ability to perform insertion/deletion editing and thereby generate mature forms of some mitochondrial mRNAs. However, these cells, when starved, were impeded in the epimastigote to trypomastigote transition. Further, they experience a short-flagella phenotype that may also be linked to differentiation. We hypothesize a scenario where levels of mature RNA species or editing in the single T. cruzi mitochondrion are linked to differentiation by a yet-unknown signaling mechanism.

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

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.

Marked for translation: long A/U tails as an interface between completion of RNA editing and ribosome recruitment.

The mechanism by which ribosomes select translatable mRNAs from the complex mixture of incompletely edited mRNAs in trypanosome mitochondria has remained a mystery. In this issue of Molecular Cell, Aphasizheva and colleagues (Aphasizheva et al., 2011) reveal a role for long 3' A/U tails in signaling ribosome recruitment to a fully edited, translatable mRNA.

circTAIL-seq, a targeted method for deep analysis of RNA 3' tails, reveals transcript-specific differences by multiple metrics.

Post-transcriptionally added RNA 3' nucleotide extensions, or tails, impose numerous regulatory effects on RNAs, including effects on RNA turnover and translation. However, efficient methods for in-depth tail profiling of a transcript of interest are still lacking, hindering available knowledge particularly of tail populations that are highly heterogeneous. Here, we developed a targeted approach, termed circTAIL-seq, to quantify both major and subtle differences of heterogeneous tail populations. As proof-of-principle, we show that circTAIL-seq quantifies the differences in tail qualities between two selected Trypanosoma brucei mitochondrial transcripts. The results demonstrate the power of the developed method in identification, discrimination, and quantification of different tail states that the population of one transcript can possess. We further show that circTAIL-seq can detect the tail characteristics for variants of transcripts that are not easily detectable by conventional approaches, such as degradation intermediates. Our findings are not only well supported by previous knowledge, but they also expand this knowledge and provide experimental evidence for previous hypotheses. In the future, this approach can be used to determine changes in tail qualities in response to environmental or internal stimuli, or upon silencing of genes of interest in mRNA-processing pathways. In summary, circTAIL-seq is an effective tool for comparing nonencoded RNA tails, especially when the tails are extremely variable or transcript of interest is low abundance.

Ribosome biogenesis requires a highly diverged XRN family 5'->3' exoribonuclease for rRNA processing in Trypanosoma brucei.

Although biogenesis of ribosomes is a crucial process in all organisms and is thus well conserved, Trypanosoma brucei ribosome biogenesis, of which maturation of rRNAs is an early step, has multiple points of divergence. Our aim was to determine whether in the processing of the pre-rRNA precursor molecule, 5'→3' exoribonuclease activity in addition to endonucleolytic cleavage is necessary in T. brucei as in other organisms. Our approach initiated with the bioinformatic identification of a putative 5'→3' exoribonuclease, XRNE, which is highly diverged from the XRN2/Rat1 enzyme responsible for rRNA processing in other organisms. Tagging this protein in vivo allowed us to classify XRNE as nucleolar by indirect immunofluorescence and identify by copurification interacting proteins, many of which were ribosomal proteins, ribosome biogenesis proteins, and/or RNA processing proteins. To determine whether XRNE plays a role in ribosome biogenesis in procyclic form cells, we inducibly depleted the protein by RNA interference. This resulted in the generation of aberrant preprocessed 18S rRNA and 5' extended 5.8S rRNA, implicating XRNE in rRNA processing. Polysome profiles of XRNE-depleted cells demonstrated abnormal features including an increase in ribosome small subunit abundance, a decrease in large subunit abundance, and defects in polysome assembly. Furthermore, the 5' extended 5.8S rRNA in XRNE-depleted cells was observed in the large subunit, monosomes, and polysomes in this gradient. Therefore, the function of XRNE in rRNA processing, presumably due to exonucleolytic activity very early in ribosome biogenesis, has consequences that persist throughout all biogenesis stages.

Functional characterization of two paralogs that are novel RNA binding proteins influencing mitochondrial transcripts of Trypanosoma brucei.

A majority of Trypanosoma brucei proteins have unknown functions, a consequence of its independent evolutionary history within the order Kinetoplastida that allowed for the emergence of several unique biological properties. Among these is RNA editing, needed for expression of mitochondrial-encoded genes. The recently discovered mitochondrial RNA binding complex 1 (MRB1) is composed of proteins with several functions in processing organellar RNA. We characterize two MRB1 subunits, referred to herein as MRB8170 and MRB4160, which are paralogs arisen from a large chromosome duplication occurring only in T. brucei. As with many other MRB1 proteins, both have no recognizable domains, motifs, or orthologs outside the order. We show that they are both novel RNA binding proteins, possibly representing a new class of these proteins. They associate with a similar subset of MRB1 subunits but not directly with each other. We generated cell lines that either individually or simultaneously target the mRNAs encoding both proteins using RNAi. Their dual silencing results in a differential effect on moderately and pan-edited RNAs, suggesting a possible functional separation of the two proteins. Cell growth persists upon RNAi silencing of each protein individually in contrast to the dual knockdown. Yet, their apparent redundancy in terms of cell viability is at odds with the finding that only one of these knockdowns results in the general degradation of pan-edited RNAs. While MRB8170 and MRB4160 share a considerable degree of conservation, our results suggest that their recent sequence divergence has led to them influencing mitochondrial mRNAs to differing degrees.

Role of PIP39 in oxidative stress response appears conserved in kinetoplastids.

Kinetoplastids, a group of flagellated protists that are often insect intestinal parasites, encounter various sources of oxidative stress. Such stressors include reactive oxygen species, both internally produced within the protist, and induced externally by host immune responses. This investigation focuses on the role of a highly conserved aspartate-based protein phosphatase, PTP-Interacting protein (PIP39) in managing oxidative stress. In addition to its well accepted role in a Trypanosoma brucei life stage transition, there is evidence of PIP39 participation in the T. brucei oxidative stress response. To examine whether this latter PIP39 role may exist more broadly, we aimed to elucidate PIP39's contribution to redox homeostasis in the monoxenous parasite Leptomonas seymouri. Utilizing CRISPR-Cas9-mediated elimination of PIP39 in conjunction with oxidative stress assays, we demonstrate that PIP39 is required for cellular tolerance to oxidative stress in L. seymouri, positing it as a putative regulatory node for adaptive stress responses. We propose that future analysis of L. seymouri PIP39 enzymatic activity, regulation, and potential localization to a specialized organelle termed a glycosome will contribute to a deeper understanding of the molecular mechanisms by which protozoan parasites adapt to oxidative environments. Our study also demonstrates success at using gene editing tools developed for Leishmania for the related L. seymouri.

Kinetoplast Genome of Leishmania spp. Is under Strong Purifying Selection.

Instability is an intriguing characteristic of many protist genomes, and trypanosomatids are not an exception in this respect. Some regions of trypanosomatid genomes evolve fast. For instance, the trypanosomatid mitochondrial (kinetoplast) genome consists of fairly conserved maxicircle and minicircle molecules that can, nevertheless, possess high nucleotide substitution rates between closely related strains. Recent experiments have demonstrated that rapid laboratory evolution can result in the non-functionality of multiple genes of kinetoplast genomes due to the accumulation of mutations or loss of critical genomic components. An example of a loss of critical components is the reported loss of entire minicircle classes in Leishmania tarentolae during laboratory cultivation, which results in an inability to generate some correctly encoded genes. In the current work, we estimated the evolutionary rates of mitochondrial and nuclear genome regions of multiple natural Leishmania spp. We analyzed synonymous and non-synonymous substitutions and, rather unexpectedly, found that the coding regions of kinetoplast maxicircles are among the most variable regions of both genomes. In addition, we demonstrate that synonymous substitutions greatly predominate among maxicircle coding regions and that most maxicircle genes show signs of purifying selection. These results imply that maxicircles in natural Leishmania populations remain functional despite their high mutation rate.

Circular mitochondrial-encoded mRNAs are a distinct subpopulation of mitochondrial mRNA in Trypanosoma brucei.

Since the first identification of circular RNA (circRNA) in viral-like systems, reports of circRNAs and their functions in various organisms, cell types, and organelles have greatly expanded. Here, we report the first evidence of circular mRNA in the mitochondrion of the eukaryotic parasite, Trypanosoma brucei . While using a circular RT-PCR technique developed to sequence mRNA tails of mitochondrial transcripts, we found that some mRNAs are circularized without an in vitro circularization step normally required to produce PCR products. Starting from total in vitro circularized RNA and in vivo circRNA, we high-throughput sequenced three transcripts from the 3' end of the coding region, through the 3' tail, to the 5' start of the coding region. We found that fewer reads in the circRNA libraries contained tails than in the total RNA libraries. When tails were present on circRNAs, they were shorter and less adenine-rich than the total population of RNA tails of the same transcript. Additionally, using hidden Markov modelling we determined that enzymatic activity during tail addition is different for circRNAs than for total RNA. Lastly, circRNA UTRs tended to be shorter and more variable than those of the same transcript sequenced from total RNA. We propose a revised model of Trypanosome mitochondrial tail addition, in which a fraction of mRNAs is circularized prior to the addition of adenine-rich tails and may act as a new regulatory molecule or in a degradation pathway.

Circular mitochondrial-encoded mRNAs are a distinct subpopulation of mitochondrial mRNA in Trypanosoma brucei.

Since the first identification of circular RNA (circRNA) in viral-like systems, reports of circRNAs and their functions in various organisms, cell types, and organelles have greatly expanded. Here, we report the first evidence, to our knowledge, of circular mRNA in the mitochondrion of the eukaryotic parasite, Trypanosoma brucei. While using a circular RT-PCR technique developed to sequence mRNA tails of mitochondrial transcripts, we found that some mRNAs are circularized without an in vitro circularization step normally required to produce PCR products. Starting from total in vitro circularized RNA and in vivo circRNA, we high-throughput sequenced three transcripts from the 3' end of the coding region, through the 3' tail, to the 5' start of the coding region. We found that fewer reads in the circRNA libraries contained tails than in the total RNA libraries. When tails were present on circRNAs, they were shorter and less adenine-rich than the total population of RNA tails of the same transcript. Additionally, using hidden Markov modelling we determined that enzymatic activity during tail addition is different for circRNAs than for total RNA. Lastly, circRNA UTRs tended to be shorter and more variable than those of the same transcript sequenced from total RNA. We propose a revised model of Trypanosome mitochondrial tail addition, in which a fraction of mRNAs is circularized prior to the addition of adenine-rich tails and may act as a new regulatory molecule or in a degradation pathway.

A novel member of the RNase D exoribonuclease family functions in mitochondrial guide RNA metabolism in Trypanosoma brucei.

RNA turnover and RNA editing are essential for regulation of mitochondrial gene expression in Trypanosoma brucei. RNA turnover is controlled in part by RNA 3' adenylation and uridylation status, with trans-acting factors also impacting RNA homeostasis. However, little is known about the mitochondrial degradation machinery or its regulation in T. brucei. We have identified a mitochondrial exoribonuclease, TbRND, whose expression is highly up-regulated in the insect proliferative stage of the parasite. TbRND shares sequence similarity with RNase D family enzymes but differs from all reported members of this family in possessing a CCHC zinc finger domain. In vitro, TbRND exhibits 3' to 5' exoribonuclease activity, with specificity toward uridine homopolymers, including the 3' oligo(U) tails of guide RNAs (gRNAs) that provide the sequence information for RNA editing. Several lines of evidence generated from RNAi-mediated knockdown and overexpression cell lines indicate that TbRND functions in gRNA metabolism in vivo. First, TbRND depletion results in gRNA tails extended by 2-3 nucleotides on average. Second, overexpression of wild type but not catalytically inactive TbRND results in a substantial decrease in the total gRNA population and a consequent inhibition of RNA editing. The observed effects on the gRNA population are specific as rRNAs, which are also 3'-uridylated, are unaffected by TbRND depletion or overexpression. Finally, we show that gRNA binding proteins co-purify with TbRND. In summary, TbRND is a novel 3' to 5' exoribonuclease that appears to have evolved a function highly specific to the mitochondrion of trypanosomes.

Using a Community-Based Participatory Approach to Address Gender Equity in Academic Medicine: The Center for Women in Medicine and Science at the University of Minnesota.

Many medical schools are instituting gender equity initiatives to address long-standing inequities (e.g., salary, leadership positions, resource distribution) between women and men in academic medicine. However, few theory-driven models exist with built-in metrics to assess the impact of gender equity initiatives. The authors describe the theory- and metric-driven process used to create the Center for Women in Medicine and Science (CWIMS) at the University of Minnesota (UMN) Medical School. An innovative theory-driven approach using community-based participatory research (CBPR) was used to create and organize CWIMS. CBPR acknowledges community members (e.g., faculty members, staff), academic organizational representatives (e.g., department heads, center directors), and administrative leaders (e.g., deans) as equal contributors in carrying out all aspects of gender equity work. CBPR values collaborative approaches that empower faculty, promote co-learning and co-creation of initiatives among all university partners, and build upon already existing community strengths and resources. Four CWIMS action groups were created using CBPR principles. The action groups are retention and recruitment; mentoring; salary, resource, and leadership equity; and strategic communications and collaborations. Faculty members across all medical school departments joined these 4 action groups to co-create and carry out all CWIMS gender equity initiatives. The process of developing the CWIMS center and action groups, the CBPR theoretical model guiding the approach, the initiatives developed by the action groups and metrics created, and the outcomes achieved to date are described. In addition, 4 lessons learned from the development of the CWIMS-use of theoretically driven and evidence-based models is key to building a sustainable organization; bottom-up and top-down engagement of partners is crucial for sustainability; passion and innovation are critical for long-term momentum; and not all faculty members and leaders will be enthusiastic about gender equity issues-are shared for the benefit of other medical schools wanting to develop similar centers.

Reintegrating Biology Through the Nexus of Energy, Information, and Matter.

Information, energy, and matter are fundamental properties of all levels of biological organization, and life emerges from the continuous flux of matter, energy, and information. This perspective piece defines and explains each of the three pillars of this nexus. We propose that a quantitative characterization of the complex interconversions between matter, energy, and information that comprise this nexus will help us derive biological insights that connect phenomena across different levels of biological organization. We articulate examples from multiple biological scales that highlight how this nexus approach leads to a more complete understanding of the biological system. Metrics of energy, information, and matter can provide a common currency that helps link phenomena across levels of biological organization. The propagation of energy and information through levels of biological organization can result in emergent properties and system-wide changes that impact other hierarchical levels. Deeper consideration of measured imbalances in energy, information, and matter can help researchers identify key factors that influence system function at one scale, highlighting avenues to link phenomena across levels of biological organization and develop predictive models of biological systems.

Mitochondrial RNA editing in Trypanoplasma borreli: New tools, new revelations.

The kinetoplastids are unicellular flagellates that derive their name from the 'kinetoplast', a region within their single mitochondrion harboring its organellar genome of high DNA content, called kinetoplast (k) DNA. Some protein products of this mitochondrial genome are encoded as cryptogenes; their transcripts require editing to generate an open reading frame. This happens through RNA editing, whereby small regulatory guide (g)RNAs direct the proper insertion and deletion of one or more uridines at each editing site within specific transcript regions. An accurate perspective of the kDNA expansion and evolution of their unique uridine insertion/deletion editing across kinetoplastids has been difficult to achieve. Here, we resolved the kDNA structure and editing patterns in the early-branching kinetoplastid Trypanoplasma borreli and compare them with those of the well-studied trypanosomatids. We find that its kDNA consists of circular molecules of about 42 kb that harbor the rRNA and protein-coding genes, and 17 different contigs of approximately 70 kb carrying an average of 23 putative gRNA loci per contig. These contigs may be linear molecules, as they contain repetitive termini. Our analysis uncovered a putative gRNA population with unique length and sequence parameters that is massive relative to the editing needs of this parasite. We validated or determined the sequence identity of four edited mRNAs, including one coding for ATP synthase 6 that was previously thought to be missing. We utilized computational methods to show that the T. borreli transcriptome includes a substantial number of transcripts with inconsistent editing patterns, apparently products of non-canonical editing. This species utilizes the most extensive uridine deletion compared to other studied kinetoplastids to enforce amino acid conservation of cryptogene products, although insertions still remain more frequent. Finally, in three tested mitochondrial transcriptomes of kinetoplastids, uridine deletions are more common in the raw mitochondrial reads than aligned to the fully edited, translationally competent mRNAs. We conclude that the organization of kDNA across known kinetoplastids represents variations on partitioned coding and repetitive regions of circular molecules encoding mRNAs and rRNAs, while gRNA loci are positioned on a highly unstable population of molecules that differ in relative abundance across strains. Likewise, while all kinetoplastids possess conserved machinery performing RNA editing of the uridine insertion/deletion type, its output parameters are species-specific.

Probabilistic models of biological enzymatic polymerization.

In this study, hierarchies of probabilistic models are evaluated for their ability to characterize the untemplated addition of adenine and uracil to the 3' ends of mitochondrial mRNAs of the human pathogen Trypanosoma brucei, and for their generative abilities to reproduce populations of these untemplated adenine/uridine "tails". We determined the most ideal Hidden Markov Models (HMMs) for this biological system. While our HMMs were not able to generatively reproduce the length distribution of the tails, they fared better in reproducing nucleotide composition aspects of the tail populations. The HMMs robustly identified distinct states of nucleotide addition that correlate to experimentally verified tail nucleotide composition differences. However they also identified a surprising subclass of tails among the ND1 gene transcript populations that is unexpected given the current idea of sequential enzymatic action of untemplated tail addition in this system. Therefore, these models can not only be utilized to reflect biological states that we already know about, they can also identify hypotheses to be experimentally tested. Finally, our HMMs supplied a way to correct a portion of the sequencing errors present in our data. Importantly, these models constitute rare simple pedagogical examples of applied bioinformatic HMMs, due to their binary emissions.

The Remarkable Metabolism of Vickermania ingenoplastis: Genomic Predictions.

A recently redescribed two-flagellar trypanosomatid Vickermania ingenoplastis is insensitive to the classical inhibitors of respiration and thrives under anaerobic conditions. Using genomic and transcriptomic data, we analyzed its genes of the core metabolism and documented that subunits of the mitochondrial respiratory complexes III and IV are ablated, while those of complexes I, II, and V are all present, along with an alternative oxidase. This explains the previously reported conversion of glucose to acetate and succinate by aerobic fermentation. Glycolytic pyruvate is metabolized to acetate and ethanol by pyruvate dismutation, whereby a unique type of alcohol dehydrogenase (shared only with Phytomonas spp.) processes an excess of reducing equivalents formed under anaerobic conditions, leading to the formation of ethanol. Succinate (formed to maintain the glycosomal redox balance) is converted to propionate by a cyclic process involving three enzymes of the mitochondrial methyl-malonyl-CoA pathway, via a cyclic process, which results in the formation of additional ATP. The unusual structure of the V. ingenoplastis genome and its similarity with that of Phytomonas spp. imply their relatedness or convergent evolution. Nevertheless, a critical difference between these two trypanosomatids is that the former has significantly increased its genome size by gene duplications, while the latter streamlined its genome.