DC-SIGN(+) Macrophages Control the Induction of Transplantation Tolerance.

Tissue effector cells of the monocyte lineage can differentiate into different cell types with specific cell function depending on their environment. The phenotype, developmental requirements, and functional mechanisms of immune protective macrophages that mediate the induction of transplantation tolerance remain elusive. Here, we demonstrate that costimulatory blockade favored accumulation of DC-SIGN-expressing macrophages that inhibited CD8(+) T cell immunity and promoted CD4(+)Foxp3(+) Treg cell expansion in numbers. Mechanistically, that simultaneous DC-SIGN engagement by fucosylated ligands and TLR4 signaling was required for production of immunoregulatory IL-10 associated with prolonged allograft survival. Deletion of DC-SIGN-expressing macrophages in vivo, interfering with their CSF1-dependent development, or preventing the DC-SIGN signaling pathway abrogated tolerance. Together, the results provide new insights into the tolerogenic effects of costimulatory blockade and identify DC-SIGN(+) suppressive macrophages as crucial mediators of immunological tolerance with the concomitant therapeutic implications in the clinic.

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.

Global variability in gene expression and alternative splicing is modulated by mitochondrial content.

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.

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.

ChimPipe: accurate detection of fusion genes and transcription-induced chimeras from RNA-seq data.

BACKGROUND: Chimeric transcripts are commonly defined as transcripts linking two or more different genes in the genome, and can be explained by various biological mechanisms such as genomic rearrangement, read-through or trans-splicing, but also by technical or biological artefacts. Several studies have shown their importance in cancer, cell pluripotency and motility. Many programs have recently been developed to identify chimeras from Illumina RNA-seq data (mostly fusion genes in cancer). However outputs of different programs on the same dataset can be widely inconsistent, and tend to include many false positives. Other issues relate to simulated datasets restricted to fusion genes, real datasets with limited numbers of validated cases, result inconsistencies between simulated and real datasets, and gene rather than junction level assessment. RESULTS: Here we present ChimPipe, a modular and easy-to-use method to reliably identify fusion genes and transcription-induced chimeras from paired-end Illumina RNA-seq data. We have also produced realistic simulated datasets for three different read lengths, and enhanced two gold-standard cancer datasets by associating exact junction points to validated gene fusions. Benchmarking ChimPipe together with four other state-of-the-art tools on this data showed ChimPipe to be the top program at identifying exact junction coordinates for both kinds of datasets, and the one showing the best trade-off between sensitivity and precision. Applied to 106 ENCODE human RNA-seq datasets, ChimPipe identified 137 high confidence chimeras connecting the protein coding sequence of their parent genes. In subsequent experiments, three out of four predicted chimeras, two of which recurrently expressed in a large majority of the samples, could be validated. Cloning and sequencing of the three cases revealed several new chimeric transcript structures, 3 of which with the potential to encode a chimeric protein for which we hypothesized a new role. Applying ChimPipe to human and mouse ENCODE RNA-seq data led to the identification of 131 recurrent chimeras common to both species, and therefore potentially conserved. CONCLUSIONS: ChimPipe combines discordant paired-end reads and split-reads to detect any kind of chimeras, including those originating from polymerase read-through, and shows an excellent trade-off between sensitivity and precision. The chimeras found by ChimPipe can be validated in-vitro with high accuracy.

The need for transparency and good practices in the qPCR literature.

Two surveys of over 1,700 publications whose authors use quantitative real-time PCR (qPCR) reveal a lack of transparent and comprehensive reporting of essential technical information. Reporting standards are significantly improved in publications that cite the Minimum Information for Publication of Quantitative Real-Time PCR Experiments (MIQE) guidelines, although such publications are still vastly outnumbered by those that do not.

Expanded CTG repeats trigger miRNA alterations in Drosophila that are conserved in myotonic dystrophy type 1 patients.

Myotonic dystrophy type 1 (DM1) is caused by the expansion of CTG repeats in the 3' untranslated region of the DMPK gene. Several missplicing events and transcriptional alterations have been described in DM1 patients. A large number of these defects have been reproduced in animal models expressing CTG repeats alone. Recent studies have also reported miRNA dysregulation in DM1 patients. In this work, a Drosophila model was used to investigate miRNA transcriptome alterations in the muscle, specifically triggered by CTG expansions. Twenty miRNAs were differentially expressed in CTG-expressing flies. Of these, 19 were down-regulated, whereas 1 was up-regulated. This trend was confirmed for those miRNAs conserved between Drosophila and humans (miR-1, miR-7 and miR-10) in muscle biopsies from DM1 patients. Consistently, at least seven target transcripts of these miRNAs were up-regulated in DM1 skeletal muscles. The mechanisms involved in dysregulation of miR-7 included a reduction of its primary precursor both in CTG-expressing flies and in DM1 patients. Additionally, a regulatory role for Muscleblind (Mbl) was also suggested for miR-1 and miR-7, as these miRNAs were down-regulated in flies where Mbl had been silenced. Finally, the physiological relevance of miRNA dysregulation was demonstrated for miR-10, since over-expression of this miRNA in Drosophila extended the lifespan of CTG-expressing flies. Taken together, our results contribute to our understanding of the origin and the role of miRNA alterations in DM1.

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/).

A Proteogenomic Approach to Unravel New Proteins Encoded in the Leishmania donovani (HU3) Genome.

The high-throughput proteomics data generated by increasingly more sensible mass spectrometers greatly contribute to our better understanding of molecular and cellular mechanisms operating in live beings. Nevertheless, proteomics analyses are based on accurate genomic and protein annotations, and some information may be lost if these resources are incomplete. Here, we show that most proteomics data may be recovered by interconnecting genomics and proteomics approaches (i.e., following a proteogenomic strategy), resulting, in turn, in an improvement of gene/protein models. In this study, we generated proteomics data from Leishmania donovani (HU3 strain) promastigotes that allowed us to detect 1908 proteins in this developmental stage on the basis of the currently annotated proteins available in public databases. However, when the proteomics data were searched against all possible open reading frames existing in the L. donovani genome, twenty new protein-coding genes could be annotated. Additionally, 43 previously annotated proteins were extended at their N-terminal ends to accommodate peptides detected in the proteomics data. Also, different post-translational modifications (phosphorylation, acetylation, methylation, among others) were found to occur in a large number of Leishmania proteins. Finally, a detailed comparative analysis of the L. donovani and Leishmania major experimental proteomes served to illustrate how inaccurate conclusions can be raised if proteomes are compared solely on the basis of the listed proteins identified in each proteome. Finally, we have created data entries (based on freely available repositories) to provide and maintain updated gene/protein models. Raw data are available via ProteomeXchange with the identifier PXD051920.

Intron retention and transcript chimerism conserved across mammals: Ly6g5b and Csnk2b-Ly6g5b as examples.

BACKGROUND: Alternative splicing (AS) is a major mechanism for modulating gene expression of an organism, allowing the synthesis of several structurally and functionally distinct mRNAs and protein isoforms from a unique gene. Related to AS is the Transcription Induced Chimerism (TIC) or Tandem Chimerism, by which chimeric RNAs between adjacent genes can be found, increasing combinatorial complexity of the proteome. The Ly6g5b gene presents particular behaviours in its expression, involving an intron retention event and being capable to form RNA chimera transcripts with the upstream gene Csnk2b. We wanted to characterise these events more deeply in four tissues in six different mammals and analyse their protein products. RESULTS: While canonical Csnk2b isoform was widely expressed, Ly6g5b canonical isoform was less ubiquitous, although the Ly6g5b first intron retained transcript was present in all the tissues and species analysed. Csnk2b-Ly6g5b chimeras were present in all the samples analysed, but with restricted expression patterns. Some of these chimeric transcripts maintained correct structural domains from Csnk2b and Ly6g5b. Moreover, we found Csnk2b, Ly6g5b, and Csnk2b-Ly6g5b transcripts that present exon skipping, alternative 5' and 3' splice site and intron retention events. These would generate truncated or aberrant proteins whose role remains unknown. Some chimeric transcripts would encode CSNK2B proteins with an altered C-terminus, which could affect its biological function broadening its substrate specificity. Over-expression of human CSNK2B, LY6G5B, and CSNK2B-LY6G5B proteins, show different patterns of post-translational modifications and cell distribution. CONCLUSIONS: Ly6g5b intron retention and Csnk2b-Ly6g5b transcript chimerism are broadly distributed in tissues of different mammals.

The transcriptome of Leishmania major in the axenic promastigote stage: transcript annotation and relative expression levels by RNA-seq.

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.

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.

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.

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.

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.

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.

A FRET-based assay for characterization of alternative splicing events using peptide nucleic acid fluorescence in situ hybridization.

We describe a quantitative method for detecting RNA alternative splicing variants that combines in situ hybridization of fluorescently labeled peptide nucleic acid (PNA) probes with confocal microscopy Förster resonance energy transfer (FRET). The use of PNA probes complementary to sequences flanking a given splice junction allows to specifically quantify, within the cell, the RNA isoform generating such splice junction by FRET measure. As a proof of concept we analyzed two alternative splicing events originating from lymphocyte antigen 6 (LY6) complex, locus G5B (LY6G5B) pre-mRNA. These are characterized by the removal of the first intron (Fully Spliced Isoform, FSI) or by retention of such intron (Intron-Retained Isoform, IRI). The use of PNA probe pairs labeled with donor (Cy3) and acceptor (Cy5) fluorophores, suitable to FRET, flanking FSI and IRI specific splice junctions specifically detected both mRNA isoforms in HeLa cells. We have observed that the method works efficiently with probes 5-11 nt apart. The data supports that this FRET-based PNA fluorescence in situ hybridization (FP-FISH) method offers a conceptually new approach for characterizing at the subcellular level not only splice variant isoform structure, location and dynamics but also potentially a wide variety of close range RNA-RNA interactions.

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.

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.