Search

Abstract Background Mammalian testis is a highly complex and heterogeneous tissue. This complexity, which mostly derives from spermatogenic cells, is reflected at the transcriptional level, with the largest number of tissue-specific genes and long noncoding RNAs (lncRNAs) compared to other tissues, and one of the highest rates of alternative splicing. Although it is known that adequate alternative-splicing patterns and stage-specific isoforms are critical for successful spermatogenesis, so far only a very limited number of reports have addressed a detailed study of alternative splicing and isoforms along the different spermatogenic stages. Results In the present work, using highly purified stage-specific testicular cell populations, we detected 33,002 transcripts expressed throughout mouse spermatogenesis not annotated so far. These include both splice variants of already annotated genes, and of hitherto unannotated genes. Using conservative criteria, we uncovered 13,471 spermatogenic lncRNAs, which reflects the still incomplete annotation of lncRNAs. A distinctive feature of lncRNAs was their lower number of splice variants compared to protein-coding ones, adding to the conclusion that lncRNAs are, in general, less complex than mRNAs. Besides, we identified 2,794 unannotated transcripts with high coding potential (including some arising from yet unannotated genes), many of which encode unnoticed putative testis-specific proteins. Some of the most interesting coding splice variants were chosen, and validated through RT-PCR. Remarkably, the largest number of stage-specific unannotated transcripts are expressed during early meiotic prophase stages, whose study has been scarcely addressed in former transcriptomic analyses. Conclusions We detected a high number of yet unannotated genes and alternatively spliced transcripts along mouse spermatogenesis, hence showing that the transcriptomic diversity of the testis is considerably higher than previously reported. This is especially prominent for specific, underrepresented stages such as those of early meiotic prophase, and its unveiling may constitute a step towards the understanding of their key events.

BackgroundDuring its life cycle, the human pathogen Trypanosoma cruzi must quickly adapt to different environments, in which the variation in the gene expression of the regulatory U-rich RNA-binding protein 1 (TcUBP1) plays a crucial role. We have previously demonstrated that the overexpression of TcUBP1 in insect-dwelling epimastigotes orchestrates an RNA regulon to promote differentiation to infective forms.MethodsIn an attempt to generate TcUBP1 knockout parasites by using CRISPR-Cas9 technology, in the present study, we obtained a variant transcript that encodes a protein with 95% overall identity and a modified N-terminal sequence. The expression of this mutant protein, named TcUBP1mut, was notably reduced compared to that of the endogenous form found in normal cells. TcUBP1mut-knockdown epimastigotes exhibited normal growth and differentiation into infective metacyclic trypomastigotes and were capable of infecting mammalian cells.ResultsWe analyzed the RNA-Seq expression profiles of these parasites and identified 276 up- and 426 downregulated genes with respect to the wildtype control sample. RNA-Seq comparison across distinct developmental stages revealed that the transcriptomic profile of these TcUBP1mut-knockdown epimastigotes significantly differs not only from that of epimastigotes in the stationary phase but also from the gene expression landscape characteristic of infective forms. This is both contrary to and consistent with the results of our recent study involving TcUBP1-overexpressing cells.ConclusionTogether, our findings demonstrate that the genes exhibiting opposite changes under overexpression and knockdown conditions unveil key mRNA targets regulated by TcUBP1. These mostly encompass transcripts that encode for trypomastigote-specific surface glycoproteins and ribosomal proteins, supporting a role for TcUBP1 in determining the molecular characteristics of the infective stage.

Previous studies have shown that overexpression of the Trypanosoma cruzi U-rich RNA-binding protein 1 (TcUBP1) in insect-dwelling epimastigotes results in a gene expression pattern resembling that of the infective form of the pathogen. Here, we used CRISPR-Cas9-induced edition of TcUBP1 and full-length protein overexpression in epimastigote cells to monitor transcriptomic changes during the epimastigote-to-metacyclic trypomastigote stage transition of T. cruzi. This dataset includes the bioinformatics analysis of three different RNA-seq samples, each with three biological replicates, showing differential mRNA abundances. The current transcriptome report has the potential to shed light on the quantitative variances in the expression of significantup- or down-regulated mRNAs as a consequence of the levels of the UBP1 protein. Raw data files were deposited at the NCBI Sequence Read Archive - SRA at http://ncbi.nlm.nih.gov/Traces/sra/sra.cgi with accession numbers PRJNA907231 and PRJNA949967.

Abstract Background The common liver fluke Fasciola hepatica is the agent of a zoonosis with significant economic consequences in livestock production worldwide, and increasing relevance to human health in developing countries. Although flukicidal drugs are available, re-infection and emerging resistance are demanding new efficient and inexpensive control strategies. Understanding the molecular mechanisms underlying the host-parasite interaction provide relevant clues in this search, while enlightening the physiological adaptations to parasitism. Genomics and transcriptomics are still in their infancy in F. hepatica, with very scarce information available from the invasive newly excysted juveniles (NEJ). Here we provide an initial glimpse to the transcriptomics of the NEJ, the first stage to interact with the mammalian host. Results We catalogued more than 500 clusters generated from the analysis of F. hepatica juvenile expressed sequence tags (EST), several of them not detected in the adult stage. A set of putative F. hepatica specific transcripts, and a group of sequences conserved exclusively in flatworms were identified. These novel sequences along with a set of parasite transcripts absent in the host genomes are putative new targets for future anti-parasitic drugs or vaccine development. Comparisons of the F. hepatica sequences with other metazoans genomes or EST databases were consistent with the basal positioning of flatworms in the bilaterian phylogeny. Notably, GC content, codon usage and amino acid frequencies are remarkably different in Schistosomes to F. hepatica and other trematodes. Functional annotation of predicted proteins showed a general representation of diverse biological functions. Besides proteases and antioxidant enzymes expected to participate in the early interaction with the host, various proteins involved in gene expression, protein synthesis, cell signaling and mitochondrial enzymes were identified. Differential expression of secreted protease gene family members between juvenile and adult stages may respond to different needs during host colonization. Conclusion The knowledge of the genes expressed by the invasive stage of Fasciola hepatica is a starting point to unravel key aspects of this parasite's biology. The integration of the emerging transcriptomics, and proteomics data and the advent of functional genomics tools in this organism are positioning F. hepatica as an interesting model for trematode biology.

Trypanosoma cruzi is a digenetic unicellular parasite that alternates between a blood-sucking insect and a mammalian, host causing Chagas disease or American trypanosomiasis. In the insect gut, the parasite differentiates from the non-replicative trypomastigote forms that arrive upon blood ingestion to the non-infective replicative epimastigote forms. Epimastigotes develop into infective non-replicative metacyclic trypomastigotes in the rectum and are delivered via the feces. In addition to these parasite stages, transitional forms have been reported. The insect-feeding behavior, characterized by few meals of large blood amounts followed by long periods of starvation, impacts the parasite population density and differentiation, increasing the transitional forms while diminishing both epimastigotes and metacyclic trypomastigotes. To understand the molecular changes caused by nutritional restrictions in the insect host, mid-exponentially growing axenic epimastigotes were cultured for more than 30 days without nutrient supplementation (prolonged starvation). We found that the parasite population in the stationary phase maintains a long period characterized by a total RNA content three times smaller than that of exponentially growing epimastigotes and a distinctive transcriptomic profile. Among the transcriptomic changes induced by nutrient restriction, we found differentially expressed genes related to managing protein quality or content, the reported switch from glucose to amino acid consumption, redox challenge, and surface proteins. The contractile vacuole and reservosomes appeared as cellular components enriched when ontology term overrepresentation analysis was carried out, highlighting the roles of these organelles in starving conditions possibly related to their functions in regulating cell volume and osmoregulation as well as metabolic homeostasis. Consistent with the quiescent status derived from nutrient restriction, genes related to DNA metabolism are regulated during the stationary phase. In addition, we observed differentially expressed genes related to the unique parasite mitochondria. Finally, our study identifies gene expression changes that characterize transitional parasite forms enriched by nutrient restriction. The analysis of the here-disclosed regulated genes and metabolic pathways aims to contribute to the understanding of the molecular changes that this unicellular parasite undergoes in the insect vector.

Introduction:Trypanosoma cruzi, Trypanosoma brucei, and Leishmania spp., commonly referred to as TriTryps, are a group of protozoan parasites that cause important human diseases affecting millions of people belonging to the most vulnerable populations worldwide. Current treatments have limited efficiencies and can cause serious side effects, so there is an urgent need to develop new control strategies. Presently, the identification and prioritization of appropriate targets can be aided by integrative genomic and computational approaches.Methods: In this work, we conducted a genome-wide multidimensional data integration strategy to prioritize drug targets. We included genomic, transcriptomic, metabolic, and protein structural data sources, to delineate candidate proteins with relevant features for target selection in drug development.Results and Discussion: Our final ranked list includes proteins shared by TriTryps and covers a range of biological functions including essential proteins for parasite survival or growth, oxidative stress-related enzymes, virulence factors, and proteins that are exclusive to these parasites. Our strategy found previously described candidates, which validates our approach as well as new proteins that can be attractive targets to consider during the initial steps of drug discovery.

In Trypanosoma cruzi DNA is packaged into chromatin by octamers of histone proteins that form nucleosomes. Transcription of protein coding genes in trypanosomes is constitutive producing polycistronic units and gene expression is primarily regulated post-transcriptionally. However, chromatin organization influences DNA dependent processes. Hence, determining nucleosome position is of uppermost importance to understand the peculiarities found in trypanosomes. To map nucleosomes genome-wide in several organisms, digestion of chromatin with micrococcal nuclease followed by deep sequencing has been applied. Nonetheless, the special requirements for cell manipulation and the uniqueness of the chromatin organization in trypanosomes entails a customized analytical approach. In this work, we adjusted this broadly used method to the hybrid reference strain, CL Brener. Particularly, we implemented an exhaustive and thorough computational workflow to overcome the difficulties imposed by this complex genome. We tested the performance of two aligners, Bowtie2 and HISAT2, and discuss their advantages and caveats. Specifically, we highlight the relevance of using the whole genome as a reference instead of the commonly used Esmeraldo-like haplotype to avoid spurious alignments. Additionally, we show that using the whole genome refines the average nucleosome representation, but also the quality of mapping for every region represented. Moreover, we show that the average nucleosome organization around trans-splicing acceptor site described before, is not just an average since the same chromatin pattern is detected for most of the represented regions. In addition, we extended the study to a non-hybrid strain applying the experimental and analytical approach to Sylvio-X10 strain. Furthermore, we provide a source code for the construction of 2D plots and heatmaps which are easy to adapt to any T. cruzi strain.

We developed a genomic surveillance program for real-time monitoring of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOCs) in Uruguay. We report on a PCR method for SARS-CoV-2 VOCs, the surveillance workflow, and multiple independent introductions and community transmission of the SARS-CoV-2 P.1 VOC in Uruguay.

Functional genomics studies through transcriptomics, translatomics and proteomics have become increasingly important tools to understand the molecular basis of biological systems in the last decade. In most cases, when these approaches are applied to the nervous system, they are centered in cell bodies or somatodendritic compartments, as these are easier to isolate and, at least in vitro, contain most of the mRNA and proteins present in all neuronal compartments. However, key functional processes and many neuronal disorders are initiated by changes occurring far away from cell bodies, particularly in axons (axopathologies) and synapses (synaptopathies). Both neuronal compartments contain specific RNAs and proteins, which are known to vary depending on their anatomical distribution, developmental stage and function, and thus form the complex network of molecular pathways required for neuron connectivity. Modifications in these components due to metabolic, environmental, and/or genetic issues could trigger or exacerbate a neuronal disease. For this reason, detailed profiling and functional understanding of the precise changes in these compartments may thus yield new insights into the still intractable molecular basis of most neuronal disorders. In the case of synaptic dysfunctions or synaptopathies, they contribute to dozens of diseases in the human brain including neurodevelopmental (i.e., autism, Down syndrome, and epilepsy) as well as neurodegenerative disorders (i.e., Alzheimer’s and Parkinson’s diseases). Histological, biochemical, cellular, and general molecular biology techniques have been key in understanding these pathologies. Now, the growing number of omics approaches can add significant extra information at a high and wide resolution level and, used effectively, can lead to novel and insightful interpretations of the biological processes at play. This review describes current approaches that use transcriptomics, translatomics and proteomic related methods to analyze the axon and presynaptic elements, focusing on the relationship that axon and synapses have with neurodegenerative diseases.

Trypanosomatids are protozoan parasites that cause devastating vector-borne human diseases. Gene expression regulation of these organisms depends on post-transcriptional control in responding to diverse environments while going through multiple developmental stages of their complex life cycles. In this scenario, non-coding RNAs (ncRNAs) are excellent candidates for a very efficient, quick, and economic strategy to regulate gene expression. The advent of high throughput RNA sequencing technologies show the presence and deregulation of small RNA fragments derived from canonical ncRNAs. This review seeks to depict the ncRNA landscape in trypanosomatids, focusing on the small RNA fragments derived from functional RNA molecules observed in RNA sequencing studies. Small RNA fragments derived from canonical ncRNAs (tsRNAs, snsRNAs, sdRNAs, and sdrRNAs) were identified in trypanosomatids. Some of these RNAs display changes in their levels associated with different environments and developmental stages, demanding further studies to determine their functional characterization and potential roles. Nevertheless, a comprehensive and detailed ncRNA annotation for most trypanosomatid genomes is still needed, allowing better and more extensive comparative and functional studies.

Uruguay is one of the few countries in the Americas that successfully contained the coronavirus disease 19 (COVID-19) epidemic during the first half of 2020. Nevertheless, the intensive human mobility across the dry border with Brazil is a major challenge for public health authorities. We aimed to investigate the origin of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) strains detected in Uruguayan localities bordering Brazil as well as to measure the viral flux across this ∼1,100 km uninterrupted dry frontier. Using complete SARS-CoV-2 genomes from the Uruguayan–Brazilian bordering region and phylogeographic analyses, we inferred the virus dissemination frequency between Brazil and Uruguay and characterized local outbreak dynamics during the first months (May–July) of the pandemic. Phylogenetic analyses revealed multiple introductions of SARS-CoV-2 Brazilian lineages B.1.1.28 and B.1.1.33 into Uruguayan localities at the bordering region. The most probable sources of viral strains introduced to Uruguay were the Southeast Brazilian region and the state of Rio Grande do Sul. Some of the viral strains introduced in Uruguayan border localities between early May and mid-July were able to locally spread and originated the first outbreaks detected outside the metropolitan region. The viral lineages responsible for Uruguayan urban outbreaks were defined by a set of between four and 11 mutations (synonymous and non-synonymous) with respect to the ancestral B.1.1.28 and B.1.1.33 viruses that arose in Brazil, supporting the notion of a rapid genetic differentiation between SARS-CoV-2 subpopulations spreading in South America. Although Uruguayan borders have remained essentially closed to non-Uruguayan citizens, the inevitable flow of people across the dry border with Brazil allowed the repeated entry of the virus into Uruguay and the subsequent emergence of local outbreaks in Uruguayan border localities. Implementation of coordinated bi-national surveillance systems is crucial to achieve an efficient control of the SARS-CoV-2 spread across this kind of highly permeable borderland regions around the world.

Treatment for Chagas disease, a parasitosis caused by Trypanosoma cruzi, has always been based on two drugs, nifurtimox and benznidazole, despite the toxic side effects described after prolonged prescription. In this work, we study a new prospective antitrypanosomal drug based on vanadium, here named VIVO(5Brsal)(aminophen). We found a good IC50 value, (3.76 ± 0.08) μM, on CL Brener epimastigotes. The analysis of cell death mechanism allowed us to rule out the implication of a mechanism based on early apoptosis or necrosis. Recovery assays revealed a trypanostatic effect, accompanied by cell shape and motility alterations. An uptake mostly associated with the insoluble fraction of the parasites was deduced through vanadium determinations. Concordantly, no drastic changes of the parasite transcriptome were detected after 6 h of treatment. Instead, proteomic analysis uncovered the modulation of proteins involved in different processes such as energy and redox metabolism, transport systems, detoxifying pathways, ribosomal protein synthesis, and proteasome protein degradation. Overall, the results here presented lead us to propose that VIVO(5Brsal)(aminophen) exerts a trypanostatic effect on T. cruzi affecting parasite insoluble proteins.

Abstract Background Fasciola hepatica is the main agent of fasciolosis, a zoonotic disease affecting livestock worldwide, and an emerging food-borne disease in humans. Even when effective treatments are available, drugs are costly and can result in tolerance, liver damage and normally they do not prevent reinfection. Drug-resistant strains in livestock have been reported in various countries and, more worryingly, drug resistance in human cases has emerged in South America. The present study aims to characterize the transcriptome of two South American resistant isolates, the Cajamarca isolate from Peru, resistant to both triclabendazole and albendazole (TCBZR/ABZR) and the Rubino isolate from Uruguay, resistant to ABZ (TCBZS/ABZR), and compare them to a sensitive strain (Cenapa, Mexico, TCBZS/ABZS) to reveal putative molecular mechanisms leading to drug resistance. Results We observed a major reduction in transcription in the Cajamarca TCBZR/ABZR isolate in comparison to the other isolates. While most of the differentially expressed genes are still unannotated, several trends could be detected. Specific reduction in the expression levels of cytoskeleton proteins was consistent with a role of tubulins as putative targets of triclabendazole (TCBZ). A marked reduction of adenylate cyclase might be underlying pleiotropic effects on diverse metabolic pathways of the parasite. Upregulation of GST mu isoforms suggests this detoxifying mechanism as one of the strategies associated with resistance. Conclusions Our results stress the value of transcriptomic approaches as a means of providing novel insights to advance the understanding of drug mode of action and drug resistance. The results provide evidence for pleiotropic variations in drug-resistant isolates consistent with early observations of TCBZ and ABZ effects and recent proteomic findings.

Protein translation is a key step in gene expression. The development of Ribosome Profiling has allowed the global analysis of this process at sub-codon resolution. In the last years the method has been applied to several models ranging from bacteria to mammalian cells yielding a surprising amount of insight on the mechanism and the regulation of translation. In this review we describe the key aspects of the experimental protocol and comment on the main conclusions raised in different models. Keywords: Ribo-seq, Translation, Translatome, Transcriptome, Ribosome profiling

Uruguay controlled the viral dissemination during the first nine months of the SARS-CoV-2 pandemic. Unfortunately, towards the end of 2020, the number of daily new cases exponentially increased. Herein, we analyzed the country-wide genetic diversity of SARS-CoV-2 between November 2020 and April 2021. We identified that the most prevalent viral variant during the first epidemic wave in Uruguay (December 2020–February 2021) was a B.1.1.28 sublineage carrying Spike mutations Q675H + Q677H, now designated as P.6, followed by lineages P.2 and P.7. P.6 probably arose around November 2020, in Montevideo, Uruguay’s capital department, and rapidly spread to other departments, with evidence of further local transmission clusters; it also spread sporadically to the USA and Spain. The more efficient dissemination of lineage P.6 with respect to P.2 and P.7 and the presence of mutations (Q675H and Q677H) in the proximity of the key cleavage site at the S1/S2 boundary suggest that P.6 may be more transmissible than other lineages co-circulating in Uruguay. Although P.6 was replaced by the variant of concern (VOC) P.1 as the predominant lineage in Uruguay since April 2021, the monitoring of the concurrent emergence of Q675H + Q677H in VOCs should be of worldwide interest.

Abstract Background Trypanosoma cruzi and Trypanosoma brucei are protozoan parasites causing Chagas disease and African sleeping sickness, displaying unique features of cellular and molecular biology. Remarkably, no canonical signals for RNA polymerase II promoters, which drive protein coding genes transcription, have been identified so far. The secondary structure of DNA has long been recognized as a signal in biological processes and more recently, its involvement in transcription initiation in Leishmania was proposed. In order to study whether this feature is conserved in trypanosomatids, we undertook a genome wide search for intrinsic DNA curvature in T. cruzi and T. brucei. Results Using a region integrated intrinsic curvature (RIIC) scoring that we previously developed, a non-random distribution of sequence-dependent curvature was observed. High RIIC scores were found to be significantly correlated with transcription start sites in T. cruzi, which have been mapped in divergent switch regions, whereas in T. brucei, the high RIIC scores correlated with sites that have been involved not only in RNA polymerase II initiation but also in termination. In addition, we observed regions with high RIIC score presenting in-phase tracts of Adenines, in the subtelomeric regions of the T. brucei chromosomes that harbor the variable surface glycoproteins genes. Conclusions In both T. cruzi and T. brucei genomes, a link between DNA conformational signals and gene expression was found. High sequence dependent curvature is associated with transcriptional regulation regions. High intrinsic curvature also occurs at the T. brucei chromosome subtelomeric regions where the recombination processes involved in the evasion of the immune host system take place. These findings underscore the relevance of indirect DNA readout in these ancient eukaryotes.

Abstract Background Small non-coding RNAs, including miRNAs, and gene silencing mediated by RNA interference have been described in free-living and parasitic lineages of flatworms, but only few key factors of the small RNA pathways have been exhaustively investigated in a limited number of species. The availability of flatworm draft genomes and predicted proteomes allowed us to perform an extended survey of the genes involved in small non-coding RNA pathways in this phylum. Results Overall, findings show that the small non-coding RNA pathways are conserved in all the analyzed flatworm linages; however notable peculiarities were identified. While Piwi genes are amplified in free-living worms they are completely absent in all parasitic species. Remarkably all flatworms share a specific Argonaute family (FL-Ago) that has been independently amplified in different lineages. Other key factors such as Dicer are also duplicated, with Dicer-2 showing structural differences between trematodes, cestodes and free-living flatworms. Similarly, a very divergent GW182 Argonaute interacting protein was identified in all flatworm linages. Contrasting to this, genes involved in the amplification of the RNAi interfering signal were detected only in the ancestral free living species Macrostomum lignano. We here described all the putative small RNA pathways present in both free living and parasitic flatworm lineages. Conclusion These findings highlight innovations specifically evolved in platyhelminths presumably associated with novel mechanisms of gene expression regulation mediated by small RNA pathways that differ to what has been classically described in model organisms. Understanding these phylum-specific innovations and the differences between free living and parasitic species might provide clues to adaptations to parasitism, and would be relevant for gene-silencing technology development for parasitic flatworms that infect hundreds of million people worldwide.

Trypanosoma cruzi, the protozoan parasite that causes Chagas' disease, exhibits peculiar biological features. Among them, the presence of a unique mitochondrion is remarkable. Even though the mitochondrial DNA constitutes up to 25% of total cellular DNA, the structure and functionality of the mitochondrion are dependent on the expression of the nuclear genome. As in other eukaryotes, specific peptide signals have been proposed to drive the mitochondrial localization of a subset of trypanosomatid proteins. However, there are mitochondrial proteins encoded in the nuclear genome that lack of a peptide signal. In other eukaryotes, alternative protein targeting to subcellular organelles via mRNA localization has also been recognized and specific mRNA localization towards the mitochondria has been described. With the aim of seeking for mitochondrial localization signals in T. cruzi, we developed a strategy to build a comprehensive database of nuclear genes encoding predicted mitochondrial proteins (MiNT) in the TriTryps (T. cruzi, T. brucei and L. major). We found that approximately 15% of their nuclear genome encodes mitochondrial products. In T. cruzi the MiNT database reaches 1438 genes and a conserved peptide signal, M(L/F) R (R/S) SS, named TryM-TaPe is found in 60% of these genes, suggesting that the canonical mRNA guidance mechanism is present. In addition, the search for compositional signals in the transcripts of T. cruzi MiNT genes produce a list, being worth to note a conserved non-translated element represented by the consensus sequence DARRVSG. Taking into account its reported interaction with the T. brucei TRRM3 protein which is enriched in the mitochondrial membrane fraction, we here suggest a putative zip code role for this element. Globally, here we provide an inventory of the mitochondrial proteins in T. cruzi and give evidence for the existence of both peptide and mRNA signals specific to nuclear encoded mitochondrial proteins.

Most signals involved in post-transcriptional regulatory networks are located in the untranslated regions (UTRs) of the mRNAs. Therefore, to deepen our understanding of gene expression regulation, delimitation of these regions with high accuracy is needed. The trypanosomatid lineage includes a variety of parasitic protozoans causing a significant worldwide burden on human health. Given their peculiar mechanisms of gene expression, these organisms depend on post-transcriptional regulation as the main level of gene expression control. In this context, the definition of the UTR regions becomes of key importance. We have developed UTR-mini-exon (UTRme), a graphical user interface (GUI) stand-alone application to identify and annotate 5′ and 3′ UTR regions in a highly accurate way. UTRme implements a multiple scoring system tailored to address the issue of false positive UTR assignment that frequently arise because of the characteristics of the intergenic regions. Even though it was developed for trypanosomatids, the tool can be used to predict 3′ sites in any eukaryote and 5′ UTRs in any organism where trans-splicing occurs (such as the model organism C. elegans). UTRme offers a way for non-bioinformaticians to precisely determine UTRs from transcriptomic data. The tool is freely available via the conda and github repositories.

Trypanosoma cruzi is the protozoan parasite causing American trypanosomiasis or Chagas disease, a neglected parasitosis with important human health impact in Latin America. The efficacy of current therapy is limited, and its toxicity is high. Since parasite proliferation is a fundamental target for rational drug design, we sought to progress into its understanding by applying a genome-wide approach. Treating a TcI linage strain with hydroxyurea, we isolated epimastigotes in late G1, S and G2/M cell cycle stages at 70% purity. The sequencing of each phase identified 305 stage-specific transcripts (1.5-fold change, p≤0.01), coding for conserved cell cycle regulated proteins and numerous proteins whose cell cycle dependence has not been recognized before. Comparisons with the parasite T. brucei and the human host reveal important differences. The meta-analysis of T. cruzi transcriptomic and ribonomic data indicates that cell cycle regulated mRNAs are subject to sub-cellular compartmentalization. Compositional and structural biases of these genes- including CAI, GC content, UTR length, and polycistron position- may contribute to their regulation. To discover nucleotide motifs responsible for the co-regulation of cell cycle regulated genes, we looked for overrepresented motifs at their UTRs and found a variant of the cell cycle sequence motif at the 3' UTR of most of the S and G2 stage genes. We additionally identified hairpin structures at the 5' UTRs of a high proportion of the transcripts, suggesting that periodic gene expression might also rely on translation initiation in T. cruzi. In summary, we report a comprehensive list of T. cruzi cell cycle regulated genes, including many previously unstudied proteins, we show evidence favoring a multi-step control of their expression, and we identify mRNA motifs that may mediate their regulation. Our results provide novel information of the T. cruzi proliferative proteins and the integrated levels of their gene expression control.

Trypanosoma cruzi

It is generally accepted that the presence of ORFs in the 5′ untranslated region of eukaryotic transcripts modulates the production of proteins by controlling the translation initiation rate of the main CDS. In trypanosomatid parasites, which almost exclusively depend on post-transcriptional mechanisms to regulate gene expression, translation has been identified as a key step. However, the mechanisms of control of translation are not fully understood. In the present work, we have annotated the 5′UTRs of the Trypanosoma cruzi genome both in epimastigotes and metacyclic trypomastigotes and, using a stringent classification approach, we identified putative regulatory uORFs in about 9% of the analyzed 5′UTRs. The translation efficiency (TE) and translational levels of transcripts containing putative repressive uORFs were found to be significantly reduced. These findings are supported by the fact that proteomic methods only identify a low number of proteins coded by transcripts containing repressive uORF. We additionally show that AUG is the main translation initiator codon of repressive uORFs in T. cruzi. Interestingly, the decrease in TE is more pronounced when the uORFs overlaps the main CDS. In conclusion, we show that the presence of the uORF and features such as initiation codon and/or location of the uORFs may be acting to fine tune translation levels in these parasites.

In the present work, we performed a comparative genome-wide analysis of 22 species representative of the main clades and lifestyles of the phylum Platyhelminthes. We selected a set of 700 orthologous genes conserved in all species, measuring changes in GC content, codon, and amino acid usage in orthologous positions. Values of 3rd codon position GC spanned over a wide range, allowing to discriminate two distinctive clusters within freshwater turbellarians, Cestodes and Trematodes respectively. Furthermore, a hierarchical clustering of codon usage data differs remarkably from the phylogenetic tree. Additionally, we detected a synonymous codon usage bias that was more dramatic in extreme GC-poor or GC-rich genomes, i.e., GC-poor Schistosomes preferred to use AT-rich terminated synonymous codons, while GC-rich M. lignano showed the opposite behavior. Interestingly, these biases impacted the amino acidic usage, with preferred amino acids encoded by codons following the GC content trend. These are associated with non-synonymous substitutions at orthologous positions. The detailed analysis of the synonymous and non-synonymous changes provides evidence for a two-hit mechanism where both mutation and selection forces drive the diverse coding strategies of flatworms.

In trypanosomatids, the RNA polymerase I (RNAPI)-dependent promoters controlling the ribosomal RNA (rRNA) genes have been well identified. Although the RNAPI transcription machinery recognizes the DNA conformation instead of the DNA sequence of promoters, no conformational study has been reported for these promoters. Here we present the in silico analysis of the intrinsic DNA curvature of the rRNA gene core promoters in Trypanosoma brucei, Trypanosoma cruzi, and Leishmania major. We found that, in spite of the absence of sequence conservation, these promoters hold conformational properties similar to other eukaryotic rRNA promoters. Our results also indicated that the intrinsic DNA curvature pattern is conserved within the Leishmania genus and also among strains of T. cruzi and T. brucei. Furthermore, we analyzed the impact of point mutations on the intrinsic curvature and their impact on the promoter activity. Furthermore, we found that the core promoters of protein-coding genes transcribed by RNAPI in T. brucei show the same conserved conformational characteristics. Overall, our results indicate that DNA intrinsic curvature of the rRNA gene core promoters is conserved in these ancient eukaryotes and such conserved curvature might be a requirement of RNAPI machinery for transcription of not only rRNA genes but also protein-coding genes.

Food borne trematodes (FBTs) are an assemblage of platyhelminth parasites transmitted through the food chain, four of which are recognized as neglected tropical diseases (NTDs). Fascioliasis stands out among the other NTDs due to its broad and significant impact on both human and animal health, as Fasciola sp., are also considered major pathogens of domesticated ruminants. Here we present a reference genome sequence of the common liver fluke, Fasciola hepatica isolated from sheep, complementing previously reported isolate from cattle. A total of 14,642 genes were predicted from the 1.14 GB genome of the liver fluke. Comparative genomics indicated that F. hepatica Oregon and related food-borne trematodes are metabolically less constrained than schistosomes and cestodes, taking advantage of the richer millieux offered by the hepatobiliary organs. Protease families differentially expanded between diverse trematodes may facilitate migration and survival within the heterogeneous environments and niches within the mammalian host. Surprisingly, the sequencing of Oregon and Uruguay F. hepatica isolates led to the first discovery of an endobacteria in this species. Two contigs from the F. hepatica Oregon assembly were joined to complete the 859,205 bp genome of a novel Neorickettsia endobacterium (nFh) closely related to the etiological agents of human Sennetsu and Potomac horse fevers. Immunohistochemical studies targeting a Neorickettsia surface protein found nFh in specific organs and tissues of the adult trematode including the female reproductive tract, eggs, the Mehlis' gland, seminal vesicle, and oral suckers, suggesting putative routes for fluke-to-fluke and fluke-to-host transmission. The genomes of F. hepatica and nFh will serve as a resource for further exploration of the biology of F. hepatica, and specifically its newly discovered trans-kingdom interaction with nFh and the impact of both species on disease in ruminants and humans.

Leishmania major is a flagellated protozoan parasite of medical importance. Like other members of the Trypanosomatidae family, it possesses unique mechanisms of gene expression such as constitutive polycistronic transcription of directional gene clusters, gene amplification, mRNA trans-splicing, and extensive editing of mitochondrial transcripts. The molecular signals underlying most of these processes remain under investigation. In order to investigate the role of DNA secondary structure signals in gene expression, we carried out a genome-wide in silico analysis of the intrinsic DNA curvature. The L. major genome revealed a lower frequency of high intrinsic curvature regions as well as inter- and intra- chromosomal distribution heterogeneity, when compared to prokaryotic and eukaryotic organisms. Using a novel method aimed at detecting region-integrated intrinsic curvature (RIIC), high DNA curvature was found to be associated with regions implicated in transcription initiation. Those include divergent strand-switch regions between directional gene clusters and regions linked to markers of active transcription initiation such as acetylated H3 histone, TRF4 and SNAP50. These findings suggest a role for DNA curvature in transcription initiation in Leishmania supporting the relevance of DNA secondary structures signals.

Genome sequences are available for the human blood flukes, Schistosoma japonicum, S. mansoni and S. haematobium. Functional genomic approaches could aid in identifying the role and importance of these newly described schistosome genes. Transgenesis is established for functional genomics in model species, which can lead to gain- or loss-of-functions, facilitate vector-based RNA interference, and represents an effective forward genetics tool for insertional mutagenesis screens. Progress toward routine transgenesis in schistosomes might be expedited if germ cells could be reliably localized in cultured schistosomes. Vasa, a member of the ATP-dependent DEAD-box RNA helicase family, is a prototypic marker of primordial germ cells and the germ line in the Metazoa. Using bioinformatics, 33 putative DEAD-box RNA helicases exhibiting conserved motifs that characterize helicases of this family were identified in the S. mansoni genome. Moreover, three of the helicases exhibited vasa-like sequences; phylogenetic analysis confirmed the three vasa-like genes-termed Smvlg1, Smvlg2, and Smvlg3-were members of the Vasa/PL10 DEAD-box subfamily. Transcripts encoding Smvlg1, Smvlg2, and Smvlg3 were cloned from cDNAs from mixed sex adult worms, and quantitative real time PCR revealed their presence in developmental stages of S. mansoni with elevated expression in sporocysts, adult females, eggs, and miracidia, with strikingly high expression in the undeveloped egg. Whole mount in situ hybridization (WISH) analysis revealed that Smvlg1, Smvlg2 and Smvlg3 were transcribed in the posterior ovary where the oocytes mature. Germ cell specific expression of schistosome vasa-like genes should provide an informative landmark for germ line transgenesis of schistosomes, etiologic agents of major neglected tropical diseases.

Background: During its life cycle, the human pathogen Trypanosoma cruzi must quickly adapt to different environments, in which the variation in the gene expression of the regulatory U-rich RNA-binding protein 1 (TcUBP1) plays a crucial role. We have previously demonstrated that the overexpression of TcUBP1 in insect-dwelling epimastigotes orchestrates an RNA regulon to promote differentiation to infective forms., Methods: In an attempt to generate TcUBP1 knockout parasites by using CRISPR-Cas9 technology, in the present study, we obtained a variant transcript that encodes a protein with 95% overall identity and a modified N-terminal sequence. The expression of this mutant protein, named TcUBP1mut, was notably reduced compared to that of the endogenous form found in normal cells. TcUBP1mut-knockdown epimastigotes exhibited normal growth and differentiation into infective metacyclic trypomastigotes and were capable of infecting mammalian cells., Results: We analyzed the RNA-Seq expression profiles of these parasites and identified 276 up- and 426 downregulated genes with respect to the wildtype control sample. RNA-Seq comparison across distinct developmental stages revealed that the transcriptomic profile of these TcUBP1mut-knockdown epimastigotes significantly differs not only from that of epimastigotes in the stationary phase but also from the gene expression landscape characteristic of infective forms. This is both contrary to and consistent with the results of our recent study involving TcUBP1-overexpressing cells., Conclusion: Together, our findings demonstrate that the genes exhibiting opposite changes under overexpression and knockdown conditions unveil key mRNA targets regulated by TcUBP1. These mostly encompass transcripts that encode for trypomastigote-specific surface glycoproteins and ribosomal proteins, supporting a role for TcUBP1 in determining the molecular characteristics of the infective stage., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Sabalette et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Infectious bovine keratoconjunctivitis (IBK) is an ocular disease that affects bovines and has significant economic and health effects worldwide. Gram negative bacteria Moraxella bovis and Moraxella bovoculi are its main etiological agents. Antimicrobial therapy against IBK is often difficult in beef and dairy herds and, although vaccines are commercially available, their efficacy is variable and dependent on local strains. The aim of this study was to analyze for the first time the genomes of Uruguayan clinical isolates of M. bovis and M. bovoculi. The genomes were de novo assembled and annotated; the genetic basis of fimbrial synthesis was analyzed and virulence factors were identified. A 94% coverage in the reference genomes of both species, and more than 80% similarity to the reference genomes were observed. The mechanism of fimbrial phase variation in M. bovis was detected, and the tfpQ orientation of these genes confirmed, in an inversion region of approximately 2.18kb. No phase variation was determined in the fimbrial gene of M. bovoculi. When virulence factors were compared between strains, it was observed that fimbrial genes have 36.2% sequence similarity. In contrast, the TonB-dependent lactoferrin/transferrin receptor exhibited the highest percentage of amino acid similarity (97.7%) between strains, followed by cytotoxins MbxA/MbvA and the ferric uptake regulator. The role of these virulence factors in the pathogenesis of IBK and their potential as vaccine components should be explored., (Copyright © 2024 Asociación Argentina de Microbiología. Publicado por Elsevier España, S.L.U. All rights reserved.)

Previous studies have shown that overexpression of the Trypanosoma cruzi U-rich RNA-binding protein 1 (TcUBP1) in insect-dwelling epimastigotes results in a gene expression pattern resembling that of the infective form of the pathogen. Here, we used CRISPR-Cas9-induced edition of TcUBP1 and full-length protein overexpression in epimastigote cells to monitor transcriptomic changes during the epimastigote-to-metacyclic trypomastigote stage transition of T. cruzi . This dataset includes the bioinformatics analysis of three different RNA-seq samples, each with three biological replicates, showing differential mRNA abundances. The current transcriptome report has the potential to shed light on the quantitative variances in the expression of significant up- or down-regulated mRNAs as a consequence of the levels of the UBP1 protein. Raw data files were deposited at the NCBI Sequence Read Archive - SRA at http://ncbi.nlm.nih.gov/Traces/sra/sra.cgi with accession numbers PRJNA907231 and PRJNA949967., (© 2024 The Authors.)

Trypanosoma brucei occupies distinct niches throughout its life cycle, within both the mammalian and tsetse fly hosts. The immunological and biochemical complexity and variability of each of these environments require a reshaping of the protein landscape of the parasite both to evade surveillance and face changing metabolic demands. In kinetoplastid protozoa, including T. brucei , posttranscriptional control mechanisms are the primary means of gene regulation, and these are often mediated by RNA-binding proteins. DRBD18 is a T. brucei RNA-binding protein that reportedly interacts with ribosomal proteins and translation factors. Here, we tested a role for DRBD18 in translational control. We validate the DRBD18 interaction with translating ribosomes and the translation initiation factor, eIF3a. We further show that DRBD18 depletion by RNA interference leads to altered polysomal profiles with a specific depletion of heavy polysomes. Ribosome profiling analysis reveals that 101 transcripts change in translational efficiency (TE) upon DRBD18 depletion: 41 exhibit decreased TE and 60 exhibit increased TE. A further 66 transcripts are buffered, that is, changes in transcript abundance are compensated by changes in TE such that the total translational output is expected not to change. In DRBD18-depleted cells, a set of transcripts that codes for procyclic form-specific proteins is translationally repressed while, conversely, transcripts that code for bloodstream form- and metacyclic form-specific proteins are translationally enhanced. RNA immunoprecipitation/qRT-PCR indicates that DRBD18 associates with members of both repressed and enhanced cohorts. These data suggest that DRBD18 contributes to the maintenance of the procyclic state through both positive and negative translational control of specific mRNAs., (© 2023 Ciganda et al.; Published by Cold Spring Harbor Laboratory Press for the RNA Society.)

In Trypanosoma cruzi DNA is packaged into chromatin by octamers of histone proteins that form nucleosomes. Transcription of protein coding genes in trypanosomes is constitutive producing polycistronic units and gene expression is primarily regulated post-transcriptionally. However, chromatin organization influences DNA dependent processes. Hence, determining nucleosome position is of uppermost importance to understand the peculiarities found in trypanosomes. To map nucleosomes genome-wide in several organisms, digestion of chromatin with micrococcal nuclease followed by deep sequencing has been applied. Nonetheless, the special requirements for cell manipulation and the uniqueness of the chromatin organization in trypanosomes entails a customized analytical approach. In this work, we adjusted this broadly used method to the hybrid reference strain, CL Brener. Particularly, we implemented an exhaustive and thorough computational workflow to overcome the difficulties imposed by this complex genome. We tested the performance of two aligners, Bowtie2 and HISAT2, and discuss their advantages and caveats. Specifically, we highlight the relevance of using the whole genome as a reference instead of the commonly used Esmeraldo-like haplotype to avoid spurious alignments. Additionally, we show that using the whole genome refines the average nucleosome representation, but also the quality of mapping for every region represented. Moreover, we show that the average nucleosome organization around trans-splicing acceptor site described before, is not just an average since the same chromatin pattern is detected for most of the represented regions. In addition, we extended the study to a non-hybrid strain applying the experimental and analytical approach to Sylvio-X10 strain. Furthermore, we provide a source code for the construction of 2D plots and heatmaps which are easy to adapt to any T. cruzi strain., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Beati et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Trypanosomes regulate gene expression mainly by using posttranscriptional mechanisms. Key factors responsible for carrying out this regulation are RNA-binding proteins, affecting subcellular localization, translation, and/or transcript stability. Trypanosoma cruzi U-rich RNA-binding protein 1 (TcUBP1) is a small protein that modulates the expression of several surface glycoproteins of the trypomastigote infective stage of the parasite. Its mRNA targets are known, but the impact of its overexpression at the transcriptome level in the insect-dwelling epimastigote cells has not yet been investigated. Thus, in the present study, by using a tetracycline-inducible system, we generated a population of TcUBP1-overexpressing parasites and analyzed its effect by RNA-Seq methodology. This allowed us to identify 793 up- and 371 downregulated genes with respect to the wildtype control sample. Among the upregulated genes, it was possible to identify members coding for the TcS superfamily, MASP, MUCI/II, and protein kinases, whereas among the downregulated transcripts, we found mainly genes coding for ribosomal, mitochondrial, and synthetic pathway proteins. RNA-Seq comparison with two previously published datasets revealed that the expression profile of this TcUBP1-overexpressing replicative epimastigote form resembles the transition to the infective metacyclic trypomastigote stage. We identified novel cis-regulatory elements in the 3'-untranslated region of the affected transcripts and confirmed that UBP1m, a signature TcUBP1 binding element previously characterized in our laboratory, is enriched in the list of stabilized genes. We can conclude that the overall effect of TcUBP1 overexpression on the epimastigote transcriptome is mainly the stabilization of mRNAs coding for proteins that are important for parasite infection., Competing Interests: Conflict of interest The authors declare that they have no conflicts of interest with the contents of this article., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Trypanosoma cruzi is a digenetic unicellular parasite that alternates between a blood-sucking insect and a mammalian, host causing Chagas disease or American trypanosomiasis. In the insect gut, the parasite differentiates from the non-replicative trypomastigote forms that arrive upon blood ingestion to the non-infective replicative epimastigote forms. Epimastigotes develop into infective non-replicative metacyclic trypomastigotes in the rectum and are delivered via the feces. In addition to these parasite stages, transitional forms have been reported. The insect-feeding behavior, characterized by few meals of large blood amounts followed by long periods of starvation, impacts the parasite population density and differentiation, increasing the transitional forms while diminishing both epimastigotes and metacyclic trypomastigotes. To understand the molecular changes caused by nutritional restrictions in the insect host, mid-exponentially growing axenic epimastigotes were cultured for more than 30 days without nutrient supplementation (prolonged starvation). We found that the parasite population in the stationary phase maintains a long period characterized by a total RNA content three times smaller than that of exponentially growing epimastigotes and a distinctive transcriptomic profile. Among the transcriptomic changes induced by nutrient restriction, we found differentially expressed genes related to managing protein quality or content, the reported switch from glucose to amino acid consumption, redox challenge, and surface proteins. The contractile vacuole and reservosomes appeared as cellular components enriched when ontology term overrepresentation analysis was carried out, highlighting the roles of these organelles in starving conditions possibly related to their functions in regulating cell volume and osmoregulation as well as metabolic homeostasis. Consistent with the quiescent status derived from nutrient restriction, genes related to DNA metabolism are regulated during the stationary phase. In addition, we observed differentially expressed genes related to the unique parasite mitochondria. Finally, our study identifies gene expression changes that characterize transitional parasite forms enriched by nutrient restriction. The analysis of the here-disclosed regulated genes and metabolic pathways aims to contribute to the understanding of the molecular changes that this unicellular parasite undergoes in the insect vector., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Smircich, Pérez-Díaz, Hernández, Duhagon and Garat.)

Introduction: Trypanosoma cruzi , Trypanosoma brucei , and Leishmania spp. , commonly referred to as TriTryps, are a group of protozoan parasites that cause important human diseases affecting millions of people belonging to the most vulnerable populations worldwide. Current treatments have limited efficiencies and can cause serious side effects, so there is an urgent need to develop new control strategies. Presently, the identification and prioritization of appropriate targets can be aided by integrative genomic and computational approaches. Methods: In this work, we conducted a genome-wide multidimensional data integration strategy to prioritize drug targets. We included genomic, transcriptomic, metabolic, and protein structural data sources, to delineate candidate proteins with relevant features for target selection in drug development. Results and Discussion: Our final ranked list includes proteins shared by TriTryps and covers a range of biological functions including essential proteins for parasite survival or growth, oxidative stress-related enzymes, virulence factors, and proteins that are exclusive to these parasites. Our strategy found previously described candidates, which validates our approach as well as new proteins that can be attractive targets to consider during the initial steps of drug discovery., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Rivara-Espasandín, Palumbo, Sosa, Radío, Turjanski, Sotelo-Silveira, Fernandez Do Porto and Smircich.)

Motivation: The use of high precision for representing quality scores in nanopore sequencing data makes these scores hard to compress and, thus, responsible for most of the information stored in losslessly compressed FASTQ files. This motivates the investigation of the effect of quality score information loss on downstream analysis from nanopore sequencing FASTQ files., Results: We polished de novo assemblies for a mock microbial community and a human genome, and we called variants on a human genome. We repeated these experiments using various pipelines, under various coverage level scenarios and various quality score quantizers. In all cases, we found that the quantization of quality scores causes little difference (or even sometimes improves) on the results obtained with the original (non-quantized) data. This suggests that the precision that is currently used for nanopore quality scores may be unnecessarily high, and motivates the use of lossy compression algorithms for this kind of data. Moreover, we show that even a non-specialized compressor, such as gzip, yields large storage space savings after the quantization of quality scores., Availability and Supplementary Information: Quantizers are freely available for download at: https://github.com/mrivarauy/QS-Quantizer., (© The Author(s) 2022. Published by Oxford University Press.)