Datasets of Iso-Seq transcripts for decoding transcriptome complexity in four Leishmania species.

The Iso-Seq technology, based on PacBio sequencing, enables the generation of high-quality, full-length transcripts, providing insights into transcriptome complexity. In this study, total RNA from promastigotes of four Leishmania species (Leishmania braziliensis, Leishmania donovani, Leishmania infantum and Leishmania major) was sequenced using Single Molecule, Real-Time (SMRT) Sequencing (PacBio) methodology. The Iso-seq transcripts were categorized as either complete or truncated according to the presence or absence of the Spliced-Leader (SL) sequence at their 5'-end, respectively. Moreover, only transcripts having a poly-A+ at their 3'-end were considered. Supplied datasets represent valuable information that may help to uncover novel transcripts and alternative splicing events in a parasite that regulates its gene expression at the post-transcriptional level. A better knowledge of gene expression regulation in Leishmania will open avenues for the development of new drugs to treat leishmaniasis, a devastating disease that has worldwide distribution. Additionally, the bioinformatics pipeline followed here may guide the analysis of Iso-Seq data derived from related trypanosomatids like Trypanosoma cruzi (Chagas disease agent) and Trypanosoma brucei (sleeping disease). © 2023 The Authors. Published by Elsevier Inc. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).

Refinement of Leishmania donovani Genome Annotations in the Light of Ribosome-Protected mRNAs Fragments (Ribo-Seq Data).

Advances in next-generation sequencing methodologies have facilitated the assembly of an ever-increasing number of genomes. Gene annotations are typically conducted via specialized software, but the most accurate results require additional manual curation that incorporates insights derived from functional and bioinformatic analyses (e.g., transcriptomics, proteomics, and phylogenetics). In this study, we improved the annotation of the Leishmania donovani (strain HU3) genome using publicly available data from the deep sequencing of ribosome-protected mRNA fragments (Ribo-Seq). As a result of this analysis, we uncovered 70 previously non-annotated protein-coding genes and improved the annotation of around 600 genes. Additionally, we present evidence for small upstream open reading frames (uORFs) in a significant number of transcripts, indicating their potential role in the translational regulation of gene expression. The bioinformatics pipelines developed for these analyses can be used to improve the genome annotations of other organisms for which Ribo-Seq data are available. The improvements provided by these studies will bring us closer to the ultimate goal of a complete and accurately annotated L. donovani genome and will enhance future transcriptomics, proteomics, and genetics studies.

Cognitive decline in diabetic mice predisposed to Alzheimer's disease is greater than in wild type.

In this work, we tested the hypothesis that the development of dementia in individuals with type 2 diabetes (T2DM) requires a genetic background of predisposition to neurodegenerative disease. As a proof of concept, we induced T2DM in middle-aged hAPP NL/F mice, a preclinical model of Alzheimer's disease. We show that T2DM produces more severe behavioral, electrophysiological, and structural alterations in these mice compared with wild-type mice. Mechanistically, the deficits are not paralleled by higher levels of toxic forms of Aß or by neuroinflammation but by a reduction in γ-secretase activity, lower levels of synaptic proteins, and by increased phosphorylation of tau. RNA-seq analysis of the cerebral cortex of hAPP NL/F and wild-type mice suggests that the former could be more susceptible to T2DM because of defects in trans-membrane transport. The results of this work, on the one hand, confirm the importance of the genetic background in the severity of the cognitive disorders in individuals with T2DM and, on the other hand, suggest, among the involved mechanisms, the inhibition of γ-secretase activity.

Leishmania infantum (JPCM5) Transcriptome, Gene Models and Resources for an Active Curation of Gene Annotations.

Leishmania infantum is one of the causative agents of visceral leishmaniases, the most severe form of leishmaniasis. An improved assembly for the L. infantum genome was published five years ago, yet delineation of its transcriptome remained to be accomplished. In this work, the transcriptome annotation was attained by a combination of both short and long RNA-seq reads. The good agreement between the results derived from both methodologies confirmed that transcript assembly based on Illumina RNA-seq and further delimitation according to the positions of spliced leader (SAS) and poly-A (PAS) addition sites is an adequate strategy to annotate the transcriptomes of Leishmania, a procedure previously used for transcriptome annotation in other Leishmania species and related trypanosomatids. These analyses also confirmed that the Leishmania transcripts boundaries are relatively slippery, showing extensive heterogeneity at the 5'- and 3'-ends. However, the use of RNA-seq reads derived from the PacBio technology (referred to as Iso-Seq) allowed the authors to uncover some complex transcription patterns occurring at particular loci that would be unnoticed by the use of short RNA-seq reads alone. Thus, Iso-Seq analysis provided evidence that transcript processing at particular loci would be more dynamic than expected. Another noticeable finding was the observation of a case of allelic heterozygosity based on the existence of chimeric Iso-Seq reads that might be generated by an event of intrachromosomal recombination. In addition, we are providing the L. infantum gene models, including both UTRs and CDS regions, that would be helpful for undertaking whole-genome expression studies. Moreover, we have built the foundations of a communal database for the active curation of both gene/transcript models and functional annotations for genes and proteins.

Assembly of a Large Collection of Maxicircle Sequences and Their Usefulness for Leishmania Taxonomy and Strain Typing.

Parasites of medical importance, such as Leishmania and Trypanosoma, are characterized by the presence of thousands of circular DNA molecules forming a structure known as kinetoplast, within the mitochondria. The maxicircles, which are equivalent to the mitochondrial genome in other eukaryotes, have been proposed as a promising phylogenetic marker. Using whole-DNA sequencing data, it is also possible to assemble maxicircle sequences as shown here and in previous works. In this study, based on data available in public databases and using a bioinformatics workflow previously reported by our group, we assembled the complete coding region of the maxicircles for 26 prototypical strains of trypanosomatid species. Phylogenetic analysis based on this dataset resulted in a robust tree showing an accurate taxonomy of kinetoplastids, which was also able to discern between closely related Leishmania species that are usually difficult to discriminate by classical methodologies. In addition, we provide a dataset of the maxicircle sequences of 60 Leishmania infantum field isolates from America, Western Europe, North Africa, and Eastern Europe. In agreement with previous studies, our data indicate that L. infantum parasites from Brazil are highly homogeneous and closely related to European strains, which were transferred there during the discovery of America. However, this study showed the existence of different L. infantum populations/clades within the Mediterranean region. A maxicircle signature for each clade has been established. Interestingly, two L. infantum clades were found coexisting in the same region of Spain, one similar to the American strains, represented by the Spanish JPCM5 reference strain, and the other, named "non-JPC like", may be related to an important leishmaniasis outbreak that occurred in Madrid a few years ago. In conclusion, the maxicircle sequence emerges as a robust molecular marker for phylogenetic analysis and species typing within the kinetoplastids, which also has the potential to discriminate intraspecific variability.

The Astonishing Large Family of HSP40/DnaJ Proteins Existing in Leishmania.

Abrupt environmental changes are faced by Leishmania parasites during transmission from a poikilothermic insect vector to a warm-blooded host. Adaptation to harsh environmental conditions, such as nutrient deprivation, hypoxia, oxidative stress and heat shock needs to be accomplished by rapid reconfiguration of gene expression and remodeling of protein interaction networks. Chaperones play a central role in the maintenance of cellular homeostasis, and they are responsible for crucial tasks such as correct folding of nascent proteins, protein translocation across different subcellular compartments, avoiding protein aggregates and elimination of damaged proteins. Nearly one percent of the gene content in the Leishmania genome corresponds to members of the HSP40 family, a group of proteins that assist HSP70s in a variety of cellular functions. Despite their expected relevance in the parasite biology and infectivity, little is known about their functions or partnership with the different Leishmania HSP70s. Here, we summarize the structural features of the 72 HSP40 proteins encoded in the Leishmania infantum genome and their classification into four categories. A review of proteomic data, together with orthology analyses, allow us to postulate cellular locations and possible functional roles for some of them. A detailed study of the members of this family would provide valuable information and opportunities for drug discovery and improvement of current treatments against leishmaniasis.

Gene Annotation and Transcriptome Delineation on a De Novo Genome Assembly for the Reference Leishmania major Friedlin Strain.

Leishmania major is the main causative agent of cutaneous leishmaniasis in humans. The Friedlin strain of this species (LmjF) was chosen when a multi-laboratory consortium undertook the objective of deciphering the first genome sequence for a parasite of the genus Leishmania. The objective was successfully attained in 2005, and this represented a milestone for Leishmania molecular biology studies around the world. Although the LmjF genome sequence was done following a shotgun strategy and using classical Sanger sequencing, the results were excellent, and this genome assembly served as the reference for subsequent genome assemblies in other Leishmania species. Here, we present a new assembly for the genome of this strain (named LMJFC for clarity), generated by the combination of two high throughput sequencing platforms, Illumina short-read sequencing and PacBio Single Molecular Real-Time (SMRT) sequencing, which provides long-read sequences. Apart from resolving uncertain nucleotide positions, several genomic regions were reorganized and a more precise composition of tandemly repeated gene loci was attained. Additionally, the genome annotation was improved by adding 542 genes and more accurate coding-sequences defined for around two hundred genes, based on the transcriptome delimitation also carried out in this work. As a result, we are providing gene models (including untranslated regions and introns) for 11,238 genes. Genomic information ultimately determines the biology of every organism; therefore, our understanding of molecular mechanisms will depend on the availability of precise genome sequences and accurate gene annotations. In this regard, this work is providing an improved genome sequence and updated transcriptome annotations for the reference L. major Friedlin strain.

Variability of the Pr77 sequence of L1Tc retrotransposon among six T. cruzi strains belonging to different discrete typing units (DTUs).

All trypanosomatid genomes are colonized by non-LTR retrotransposons which exhibit a highly conserved 77-nt sequence at their 5' ends, known as the Pr77-hallmark (Pr77). The wide distribution of Pr77 is expected to be related to the gene regulation processes in these organisms as it has promoter and HDV-like ribozyme activities at the DNA and RNA levels, respectively. The identification of Pr77 hallmark-bearing retrotransposons and the study of the associations of mobile elements with relevant genes have been analyzed in the genomes of six strains of Trypanosoma cruzi belonging to different discrete typing units (DTUs) and with different geographical origins and host/vectors. The genomes have been sequenced, assembled and annotated. BUSCO analyses indicated a good quality for the assemblies that were used in comparative analyses. The results show differences among the six genomes in the copy number of genes related to virulence processes, the abundance of retrotransposons bearing the Pr77 sequence and the presence of the Pr77 hallmarks not associated with retroelements. The analyses also show frequent associations of Pr77-bearing retrotransposons and single Pr77 hallmarks with genes coding for trans-sialidases, RHS, MASP or hypothetical proteins, showing variable proportion depending on the type of retroelement, gene class and parasite strain. These differences in the genomic distribution of active retroelements and other Pr77-containing elements have shaped the genome architecture of these six strains and might be contributing to the phenotypic variability existing among them.

The Experimental Proteome of Leishmania infantum Promastigote and Its Usefulness for Improving Gene Annotations.

Leishmania infantum causes visceral leishmaniasis (kala-azar), the most severe form of leishmaniasis, which is lethal if untreated. A few years ago, the re-sequencing and de novo assembling of the L. infantum (JPCM5 strain) genome was accomplished, and now we aimed to describe and characterize the experimental proteome of this species. In this work, we performed a proteomic analysis from axenic cultured promastigotes and carried out a detailed comparison with other Leishmania experimental proteomes published to date. We identified 2352 proteins based on a search of mass spectrometry data against a database built from the six-frame translated genome sequence of L. infantum. We detected many proteins belonging to organelles such as glycosomes, mitochondria, or flagellum, as well as many metabolic enzymes and many putative RNA binding proteins and molecular chaperones. Moreover, we listed some proteins presenting post-translational modifications, such as phosphorylations, acetylations, and methylations. On the other hand, the identification of peptides mapping to genomic regions previously annotated as non-coding allowed for the correction of annotations, leading to the N-terminal extension of protein sequences and the uncovering of eight novel protein-coding genes. The alliance of proteomics, genomics, and transcriptomics has resulted in a powerful combination for improving the annotation of the L. infantum reference genome.

Liquid biopsy by NGS: Differential presence of exons (DPE) is related to metastatic potential in a colon-cancer model in the rat.

Differential presence of exons (DPE) is a method of interpretation of exome sequencing, which has been proposed to design a predictive algorithm with clinical value in patients with colorectal cancer (CRC). The goal of the present study was to examine the reproducibility in a rat model of metastatic colon cancer. DHD/K12-TRb cells were injected in syngenic immunocompetent BD-IX rats. Cells were from two stocks with low and normal metastatic potential, and injected into two separate groups of rats. Five to ten weeks after injection, blood samples were taken prior euthanasia and whole exome sequencing performed. Through DPE analysis, we identified a set of exons whose differential presence in plasma allowed us to compare both groups of tumor-bearing animals. A verification test was performed to confirm that the algorithm was able to classify extracted samples into their corresponding groups of origin. The highest mean probability was 0.8954. In conclusion, the DPE analysis in tumor-bearing animals was able to discriminate between different disease status, which fully supports previous results in CRC patients.

Trypanosoma cruzi Ikiakarora (TcIII) Draft Genome Sequence.

Trypanosoma cruzi shows a genetic diversity that has been associated with the variability of clinical manifestations, geographical distribution, and preferential parasite-vector interactions. In an effort to better understand this genetic variability, here, the draft genome of T. cruzi strain Ikiakarora (discrete typing unit TcIII), which has been associated with the sylvatic cycle, is reported.

Draft Genome Sequence of the Trypanosoma cruzi B. M. López Strain (TcIa), Isolated from a Colombian Patient.

Trypanosoma cruzi parasite strains are classified into six lineages (discrete typing units TcI to TcVI). The broad genetic diversity of T. cruzi strains has an influence on the development of the host response and pathogenesis, as well as drug susceptibility. Here, the draft genome of the T. cruzi B. M. López strain (TcIa) is reported.

Proteins interacting with Leishmania major PUF1: A proteomic dataset.

Leishmania parasites must deal with stressful environmental conditions (thermal, nutritional and oxidative) along their digenetic life cycles. This requires drastic changes in gene expression, which in this parasite occurs mainly through post-transcriptional mechanisms involving RNA binding proteins (RBPs). PUF proteins, a class of RBPs existing in most eukaryotic organisms, might play too an essential role modulating the fate of mRNAs and regulating gene expression in Leishmania parasites. A proteomic approach to identify putative protein partners (interactome) of the Leishmania major PUF1 protein was performed. The PUF1 interactome was characterized by co-immunoprecipitation using L. major cellular extracts and an anti-LiPUF1 antibody, and a subsequent analysis of the co-immunoprecipitated proteins by mass-spectrometry, identifying 90 LmPUF1 candidate partners. Remarkably, many of the identified proteins are other RBPs and/or putative P bodies and mRNA-exporting machinery components. Raw mass spectrometry data are available via ProteomeXchange with identifier PXD022581.

Leishmania Mitochondrial Genomes: Maxicircle Structure and Heterogeneity of Minicircles.

The mitochondrial DNA (mtDNA), which is present in almost all eukaryotic organisms, is a useful marker for phylogenetic studies due to its relative high conservation and its inheritance manner. In Leishmania and other trypanosomatids, the mtDNA (also referred to as kinetoplast DNA or kDNA) is composed of thousands of minicircles and a few maxicircles, catenated together into a complex network. Maxicircles are functionally similar to other eukaryotic mtDNAs, whereas minicircles are involved in RNA editing of some maxicircle-encoded transcripts. Next-generation sequencing (NGS) is increasingly used for assembling nuclear genomes and, currently, a large number of genomic sequences are available. However, most of the time, the mitochondrial genome is ignored in the genome assembly processes. The aim of this study was to develop a pipeline to assemble Leishmania minicircles and maxicircle DNA molecules, exploiting the raw data generated in the NGS projects. As a result, the maxicircle molecules and the plethora of minicircle classes for Leishmania major, Leishmania infantum and Leishmania braziliensis have been characterized. We have observed that whereas the heterogeneity of minicircle sequences existing in a single cell hampers their use for Leishmania typing and classification, maxicircles emerge as an extremely robust genetic marker for taxonomic studies within the clade of kinetoplastids.

Analysis by RNA-seq of transcriptomic changes elicited by heat shock in Leishmania major.

Besides their medical relevance, Leishmania is an adequate model for studying post-transcriptional mechanisms of gene expression. In this microorganism, mRNA degradation/stabilization mechanisms together with translational control and post-translational modifications of proteins are the major drivers of gene expression. Leishmania parasites develop as promastigotes in sandflies and as amastigotes in mammalians, and during host transmission, the parasite experiences a sudden temperature increase. Here, changes in the transcriptome of Leishmania major promastigotes after a moderate heat shock were analysed by RNA-seq. Several of the up-regulated transcripts code for heat shock proteins, other for proteins previously reported to be amastigote-specific and many for hypothetical proteins. Many of the transcripts experiencing a decrease in their steady-state levels code for transporters, proteins involved in RNA metabolism or translational factors. In addition, putative long noncoding RNAs were identified among the differentially expressed transcripts. Finally, temperature-dependent changes in the selection of the spliced leader addition sites were inferred from the RNA-seq data, and particular cases were further validated by RT-PCR and Northern blotting. This study provides new insights into the post-transcriptional mechanisms by which Leishmania modulate gene expression.

Complete assembly of the Leishmania donovani (HU3 strain) genome and transcriptome annotation.

Leishmania donovani is a unicellular parasite that causes visceral leishmaniasis, a fatal disease in humans. In this study, a complete assembly of the genome of L. donovani is provided. Apart from being the first published genome of this strain (HU3), this constitutes the best assembly for an L. donovani genome attained to date. The use of a combination of sequencing platforms enabled to assemble, without any sequence gap, the 36 chromosomes for this species. Additionally, based on this assembly and using RNA-seq reads derived from poly-A + RNA, the transcriptome for this species, not yet available, was delineated. Alternative SL addition sites and heterogeneity in the poly-A addition sites were commonly observed for most of the genes. After a complete annotation of the transcriptome, 2,410 novel transcripts were defined. Additionally, the relative expression for all transcripts present in the promastigote stage was determined. Events of cis-splicing have been documented to occur during the maturation of the transcripts derived from genes LDHU3_07.0430 and LDHU3_29.3990. The complete genome assembly and the availability of the gene models (including annotation of untranslated regions) are important pieces to understand how differential gene expression occurs in this pathogen, and to decipher phenotypic peculiarities like tissue tropism, clinical disease, and drug susceptibility.

Complete and de novo assembly of the Leishmania braziliensis (M2904) genome

Leishmania braziliensis is the etiological agent of American mucosal leishmaniasis, one of the most severe clinical forms of leishmaniasis. Here, we report the assembly of the L. braziliensis (M2904) genome into 35 continuous chromosomes. Also, the annotation of 8395 genes is provided. The public availability of this information will contribute to a better knowledge of this pathogen and help in the search for vaccines and novel drug targets aimed to control the disease caused by this Leishmania species.

Complete and de novo assembly of the Leishmania braziliensis (M2904) genome.

Leishmania braziliensis is the etiological agent of American mucosal leishmaniasis, one of the most severe clinical forms of leishmaniasis. Here, we report the assembly of the L. braziliensis (M2904) genome into 35 continuous chromosomes. Also, the annotation of 8395 genes is provided. The public availability of this information will contribute to a better knowledge of this pathogen and help in the search for vaccines and novel drug targets aimed to control the disease caused by this Leishmania species.

Genomic and transcriptomic alterations in Leishmania donovani lines experimentally resistant to antileishmanial drugs.

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.

Deep-Sequencing Analysis of the Dynamics of HIV-1 Quasiespecies in Naive Patients during a Short Exposure to Maraviroc.

In this study, we have characterized quasispecies dynamics and the evolution of viral tropism in naive HIV-1-infected patients treated with a short course of maraviroc monotherapy (ClinicalTrials.gov registration no. NCT01060618) independently of the tropism of the infecting virus. We randomly selected 20 patients infected with viruses displaying different basal tropisms-10 carrying R5 and 10 carrying dual/mixed X4 (DM/X4) viruses-at recruitment as determined by phenotypic assay (Trofile). Evolution of viral quasiespecies at the end of treatment was determined by ultradeep sequencing of the V3 region using a 454 Life Sciences Platform and geno2pheno (g2p) algorithm for viral tropism prediction. The false-positive rate (FPR) that defines the probability of classifying an R5 virus falsely as X4 was set at 10%. X4-specific HIV-1 viral load (VL) was calculated from sequences with an FPR of <3.75%. Virological response as defined as >1-log10 copies/ml reduction in VL was detected in 70% of patients independently of the basal tropism of the infecting virus. Viral tropism remained stable, and nonsignificant differences in FPR values before and after treatment were found for the majority of patients in both tropism groups. Only three patients (one with R5 and two with DM/X4 viruses) showed an increased (>1 log) X4 VL, and one patient harboring a DM/X4-tropic virus displayed a significant reduction in FPR values at the end of treatment. Fast changes in the composition of viral populations were observed in all patients after 10 days of maraviroc (MVC) monotherapy treatment, and a complete replacement of viral quasiespecies was found in 3/10 patients carrying R5-using viruses and 4/10 patients carrying DM/X4-using viruses.IMPORTANCE Initiation of treatment with maraviroc requires previous determination of viral tropism by genotypic or phenotypic methods because of the risk of treatment failure and selection of DM/X4-tropic variants. In this study, we confirm previous work showing that the virologic response to maraviroc is independent of basal tropism. By deep-sequencing analysis, we determined that fast changes in viral populations were due to the emergence of minority variants in some patients whereas in others generation of new strains was detected. The risk of DM/X4 selection was very low as FPR values remained stable, and only one patient showed a detrimental switch to DM/X4 variants. Our data show that some DM/X4 viruses are sensitive to maraviroc treatment probably because only a low proportion of DM/X4 viruses use preferentially the X4 receptor and contain authentically maraviroc-resistant viruses that are not accurately detected by current assays.