Noise in gene expression is a main determinant of phenotypic variability. Increasing experimental evidence suggests that genome-wide cellular constraints largely contribute to the heterogeneity observed in gene products. It is still unclear, however, which global factors affect gene expression noise and to what extent. Since eukaryotic gene expression is an energy demanding process, differences in the energy budget of each cell could determine gene expression differences. Here, we quantify the contribution of mitochondrial variability (a natural source of ATP variation) to global variability in gene expression. We find that changes in mitochondrial content can account for ∼50% of the variability observed in protein levels. This is the combined result of the effect of mitochondria dosage on transcription and translation apparatus content and activities. Moreover, we find that mitochondrial levels have a large impact on alternative splicing, thus modulating both the abundance and type of mRNAs. A simple mathematical model in which mitochondrial content simultaneously affects transcription rate and splicing site choice can explain the alternative splicing data. The results of this study show that mitochondrial content (and/or probably function) influences mRNA abundance, translation, and alternative splicing, which ultimately affects cellular phenotype.

Background: Although the genome sequence of the protozoan parasite Leishmania major was determined several years ago, the knowledge of its transcriptome was incomplete, both regarding the real number of genes and their primary structure. Results: Here, we describe the first comprehensive transcriptome analysis of a parasite from the genus Leishmania. Using high-throughput RNA sequencing (RNA-seq), a total of 10285 transcripts were identified, of which 1884 were considered novel, as they did not match previously annotated genes. In addition, our data indicate that current annotations should be modified for many of the genes. The detailed analysis of the transcript processing sites revealed extensive heterogeneity in the spliced leader (SL) and polyadenylation addition sites. As a result, around 50% of the genes presented multiple transcripts differing in the length of the UTRs, sometimes in the order of hundreds of nucleotides. This transcript heterogeneity could provide an additional source for regulation as the different sizes of UTRs could modify RNA stability and/or influence the efficiency of RNA translation. In addition, for the first time for the Leishmania major promastigote stage, we are providing relative expression transcript levels. Conclusions: This study provides a concise view of the global transcriptome of the L. major promastigote stage, providing the basis for future comparative analysis with other development stages or other Leishmania species.

Leishmaniasis is a serious medical issue in many countries around the World, but it remains largely neglected in terms of research investment for developing new control and treatment measures. No vaccines exist for human use, and the chemotherapeutic agents currently used are scanty. Furthermore, for some drugs, resistance and treatment failure are increasing to alarming levels. The aim of this work was to identify genomic and trancriptomic alterations associated with experimental resistance against the common drugs used against VL: trivalent antimony (SbIII, S line), amphotericin B (AmB, A line), miltefosine (MIL, M line) and paromomycin (PMM, P line). A total of 1006 differentially expressed transcripts were identified in the S line, 379 in the A line, 146 in the M line, and 129 in the P line. Also, changes in ploidy of chromosomes and amplification/deletion of particular regions were observed in the resistant lines regarding the parental one. A series of genes were identified as possible drivers of the resistance phenotype and were validated in both promastigotes and amastigotes from Leishmania donovani, Leishmania infantum and Leishmania major species. Remarkably, a deletion of the gene LinJ.36.2510 (coding for 24-sterol methyltransferase, SMT) was found to be associated with AmB-resistance in the A line. In the P line, a dramatic overexpression of the transcripts LinJ.27.T1940 and LinJ.27.T1950 that results from a massive amplification of the collinear genes was suggested as one of the mechanisms of PMM resistance. This conclusion was reinforced after transfection experiments in which significant PMM-resistance was generated in WT parasites over-expressing either gene LinJ.27.1940 (coding for a D-lactate dehydrogenase-like protein, D-LDH) or gene LinJ.27.1950 (coding for an aminotransferase of branched-chain amino acids, BCAT). This work allowed to identify new drivers, like SMT, the deletion of which being associated with resistance to AmB, and the tandem D-LDH-BCAT, the amplification of which being related to PMM resistance.

Chagas disease is a complex illness caused by the protozoan Trypanosoma cruzi displaying highly diverse clinical outcomes. In this sense, the genome sequence elucidation and comparison between strains may lead to disease understanding. Here, two new T. cruzi strains, have been sequenced, Y using Illumina and Bug2148 using PacBio, assembled, analyzed and compared with the T. cruzi annotated genomes available to date. The assembly stats from the new sequences show effective improvement of T. cruzi genome over the actual ones. Such as, the largest contig assembled (1.3 Mb in Bug2148) in de novo attempts and the highest mean assembly coverage (71X for Y). Our analysis reveals a new genomic expansion and greater complexity for those multi-copy gene families related to infection process and disease development, such as Trans-sialidases, Mucins and Mucin Associated Surface Proteins, among others. On one side, we demonstrate that multi-copy gene families are located near telomeric regions of the “chromosome-like” 1.3 Mb contig assembled of Bug2148, where they likely suffer high evolutive pressure. On the other hand, we identified several strain-specific single copy genes that might help to understand the differences in infectivity and physiology among strains. In summary, our results indicate that T. cruzi has a complex genomic architecture that may have promoted its evolution.