Four layer multi-omics reveals molecular responses to aneuploidy in Leishmania.

Aneuploidy causes system-wide disruptions in the stochiometric balances of transcripts, proteins, and metabolites, often resulting in detrimental effects for the organism. The protozoan parasite Leishmania has an unusually high tolerance for aneuploidy, but the molecular and functional consequences for the pathogen remain poorly understood. Here, we addressed this question in vitro and present the first integrated analysis of the genome, transcriptome, proteome, and metabolome of highly aneuploid Leishmania donovani strains. Our analyses unambiguously establish that aneuploidy in Leishmania proportionally impacts the average transcript- and protein abundance levels of affected chromosomes, ultimately correlating with the degree of metabolic differences between closely related aneuploid strains. This proportionality was present in both proliferative and non-proliferative in vitro promastigotes. However, as in other Eukaryotes, we observed attenuation of dosage effects for protein complex subunits and in addition, non-cytoplasmic proteins. Differentially expressed transcripts and proteins between aneuploid Leishmania strains also originated from non-aneuploid chromosomes. At protein level, these were enriched for proteins involved in protein metabolism, such as chaperones and chaperonins, peptidases, and heat-shock proteins. In conclusion, our results further support the view that aneuploidy in Leishmania can be adaptive. Additionally, we believe that the high karyotype diversity in vitro and absence of classical transcriptional regulation make Leishmania an attractive model to study processes of protein homeostasis in the context of aneuploidy and beyond.

Acquired non-thermal plasma resistance mediates a shift towards aerobic glycolysis and ferroptotic cell death in melanoma.

AIMS: To gain insights into the underlying mechanisms of NTP therapy sensitivity and resistance, using the first-ever NTP-resistant cell line derived from sensitive melanoma cells (A375). METHODS: Melanoma cells were exposed to NTP and re-cultured for 12 consecutive weeks before evaluation against the parental control cells. Whole transcriptome sequencing analysis was performed to identify differentially expressed genes and enriched molecular pathways. Glucose uptake, extracellular lactate, media acidification, and mitochondrial respiration was analyzed to determine metabolic changes. Cell death inhibitors were used to assess the NTP-induced cell death mechanisms, and apoptosis and ferroptosis was further validated via Annexin V, Caspase 3/7, and lipid peroxidation analysis. RESULTS: Cells continuously exposed to NTP became 10 times more resistant to NTP compared to the parental cell line of the same passage, based on their half-maximal inhibitory concentration (IC50). Sequencing and metabolic analysis indicated that NTP-resistant cells had a preference towards aerobic glycolysis, while cell death analysis revealed that NTP-resistant cells exhibited less apoptosis but were more vulnerable to lipid peroxidation and ferroptosis. CONCLUSIONS: A preference towards aerobic glycolysis and ferroptotic cell death are key physiological changes in NTP-resistance cells, which opens new avenues for further, in-depth research into other cancer types.

A new Plasmodium vivax reference genome for South American isolates.

BACKGROUND: Plasmodium vivax is the second most important cause of human malaria worldwide, and accounts for the majority of malaria cases in South America. A high-quality reference genome exists for Papua Indonesia (PvP01) and Thailand (PvW1), but is lacking for South America. A reference genome specifically for South America would be beneficial though, as P. vivax is a genetically diverse parasite with geographical clustering. RESULTS: This study presents a new high-quality assembly of a South American P. vivax isolate, referred to as PvPAM (P. vivax Peruvian AMazon). The genome was obtained from a low input patient sample from the Peruvian Amazon and sequenced using PacBio technology, resulting in a highly complete assembly with 6497 functional genes. Telomeric ends were present in 17 out of 28 chromosomal ends, and additional (sub)telomeric regions are present in 12 unassigned contigs. A comparison of multigene families between PvPAM and the PvP01 genome revealed remarkable variation in vir genes, and the presence of merozoite surface proteins (MSP) 3.6 and 3.7. Three dhfr and dhps drug resistance associated mutations are present in PvPAM, similar to those found in other Peruvian isolates. Mapping of publicly available South American whole genome sequencing (WGS) data to PvPAM resulted in significantly fewer variants and truncated reads compared to the use of PvP01 or PvW1 as reference genomes. To minimize the number of core genome variants in non-South American samples, PvW1 is most suited for Southeast Asian isolates, both PvPAM and PvW1 are suited for South Asian isolates, and PvPAM is recommended for African isolates. Interestingly, non-South American samples still contained the least subtelomeric variants when mapped to PvPAM, indicating high quality of the PvPAM subtelomeric regions. CONCLUSIONS: Our findings show that the PvPAM reference genome more accurately represents South American P. vivax isolates in comparison to PvP01 and PvW1. In addition, PvPAM has a high level of completeness, and contains a similar number of annotated genes as PvP01 or PvW1. The PvPAM genome therefore will be a valuable resource to improve future genomic analyses on P. vivax isolates from the South American continent.

Variants in Bedaquiline-Candidate-Resistance Genes: Prevalence in Bedaquiline-Naive Patients, Effect on MIC, and Association with Mycobacterium tuberculosis Lineage.

Studies have shown that variants in bedaquiline-resistance genes can occur in isolates from bedaquiline-naive patients. We assessed the prevalence of variants in all bedaquiline-candidate-resistance genes in bedaquiline-naive patients, investigated the association between these variants and lineage, and the effect on phenotype. We used whole-genome sequencing to identify variants in bedaquiline-resistance genes in isolates from 509 bedaquiline treatment naive South African tuberculosis patients. A phylogenetic tree was constructed to investigate the association with the isolate lineage background. Bedaquiline MIC was determined using the UKMYC6 microtiter assay. Variants were identified in 502 of 509 isolates (98.6%), with the highest (85%) prevalence of variants in the Rv0676c (mmpL5) gene. We identified 36 unique variants, including 19 variants not reported previously. Only four isolates had a bedaquiline MIC equal to or above the epidemiological cut-off value of 0.25 µg/mL. Phylogenetic analysis showed that 14 of the 15 variants observed more than once occurred monophyletically in one Mycobacterium tuberculosis (sub)lineage. The bedaquiline MIC differed between isolates belonging to lineage 2 and 4 (Fisher's exact test, P = 0.0004). The prevalence of variants in bedaquiline-resistance genes in isolates from bedaquiline-naive patients is high, but very few (<2%) isolates were phenotypically resistant. We found an association between variants in bedaquiline resistance genes and Mycobacterium tuberculosis (sub)lineage, resulting in a lineage-dependent difference in bedaquiline phenotype. Future studies should investigate the impact of the presence of variants on bedaquiline-resistance acquisition and treatment outcome.

Nanopore Sequencing for Mycobacterium tuberculosis: a Critical Review of the Literature, New Developments, and Future Opportunities.

The next-generation, short-read sequencing technologies that generate comprehensive, whole-genome data with single nucleotide resolution have already advanced tuberculosis diagnosis, treatment, surveillance, and source investigation. Their high costs, tedious and lengthy processes, and large equipment remain major hurdles for research use in high tuberculosis burden countries and implementation into routine care. The portable next-generation sequencing devices developed by Oxford Nanopore Technologies (ONT) are attractive alternatives due to their long-read sequence capability, compact low-cost hardware, and continued improvements in accuracy and throughput. A systematic review of the published literature demonstrated limited uptake of ONT sequencing in tuberculosis research and clinical care. Of the 12 eligible articles presenting ONT sequencing data on at least one Mycobacterium tuberculosis sample, four addressed software development for long-read ONT sequencing data with potential applications for M. tuberculosis. Only eight studies presented results of ONT sequencing of M. tuberculosis, of which five performed whole-genome and three did targeted sequencing. Based on these findings, we summarize the standard processes, reflect on the current limitations of ONT sequencing technology, and the research needed to overcome the main hurdles. The low capital cost, portable nature and continued improvement in the performance of ONT sequencing make it an attractive option for sequencing for research and clinical care, but limited data are available on its application in the tuberculosis field. Important research investment is needed to unleash the full potential of ONT sequencing for tuberculosis research and care.

Sorted B cell transcriptomes point towards actively regulated B cell responses during ongoing chronic hepatitis B infections.

The natural course of chronic hepatitis B virus (HBV) infections follows distinct clinical disease phases, characterized by fluctuating levels of serum HBV DNA and ALT. The immune cells and their features that govern these clinical disease transitions remain unknown. In the current study, we performed RNA sequencing on purified B cells from blood (n = 42) and liver (n = 10) of healthy controls and chronic HBV patients. We found distinct gene expression profiles between healthy controls and chronic HBV patients, as evidenced by 190 differentially expressed genes (DEG), but also between the clinical phenotypes of a chronic HBV infection (17-110 DEG between each phase). Numerous immune pathways, including the B cell receptor pathway were upregulated in liver B cells when compared to peripheral B cells. Further investigation of the detected DEG suggested an activation of B cells during HBeAg seroconversion and an active regulation of B cell signalling in the liver.

Ferroptosis Induction in Multiple Myeloma Cells Triggers DNA Methylation and Histone Modification Changes Associated with Cellular Senescence.

Disease relapse and therapy resistance remain key challenges in treating multiple myeloma. Underlying (epi-)mutational events can promote myelomagenesis and contribute to multi-drug and apoptosis resistance. Therefore, compounds inducing ferroptosis, a form of iron and lipid peroxidation-regulated cell death, are appealing alternative treatment strategies for multiple myeloma and other malignancies. Both ferroptosis and the epigenetic machinery are heavily influenced by oxidative stress and iron metabolism changes. Yet, only a limited number of epigenetic enzymes and modifications have been identified as ferroptosis regulators. In this study, we found that MM1 multiple myeloma cells are sensitive to ferroptosis induction and epigenetic reprogramming by RSL3, irrespective of their glucocorticoid-sensitivity status. LC-MS/MS analysis revealed the formation of non-heme iron-histone complexes and altered expression of histone modifications associated with DNA repair and cellular senescence. In line with this observation, EPIC BeadChip measurements of significant DNA methylation changes in ferroptotic myeloma cells demonstrated an enrichment of CpG probes located in genes associated with cell cycle progression and senescence, such as Nuclear Receptor Subfamily 4 Group A member 2 (NR4A2). Overall, our data show that ferroptotic cell death is associated with an epigenomic stress response that might advance the therapeutic applicability of ferroptotic compounds.

Reflections on a journey: a retrospective of the ISCB Student Council symposium series.

This article describes the motivation, origin and evolution of the student symposia series organised by the ISCB Student Council. The meeting series started thirteen years ago in Madrid and has spread to four continents. The article concludes with the highlights of the most recent edition of annual Student Council Symposium held in conjunction with the 25th Conference on Intelligent Systems for Molecular Biology and the 16th European Conference on Computational Biology, in Prague, in July 2017.

On the feasibility of mining CD8+ T cell receptor patterns underlying immunogenic peptide recognition.

Current T cell epitope prediction tools are a valuable resource in designing targeted immunogenicity experiments. They typically focus on, and are able to, accurately predict peptide binding and presentation by major histocompatibility complex (MHC) molecules on the surface of antigen-presenting cells. However, recognition of the peptide-MHC complex by a T cell receptor (TCR) is often not included in these tools. We developed a classification approach based on random forest classifiers to predict recognition of a peptide by a T cell receptor and discover patterns that contribute to recognition. We considered two approaches to solve this problem: (1) distinguishing between two sets of TCRs that each bind to a known peptide and (2) retrieving TCRs that bind to a given peptide from a large pool of TCRs. Evaluation of the models on two HIV-1, B*08-restricted epitopes reveals good performance and hints towards structural CDR3 features that can determine peptide immunogenicity. These results are of particular importance as they show that prediction of T cell epitope and T cell epitope recognition based on sequence data is a feasible approach. In addition, the validity of our models not only serves as a proof of concept for the prediction of immunogenic T cell epitopes but also paves the way for more general and high-performing models.

Experimental resistance to drug combinations in Leishmania donovani: metabolic and phenotypic adaptations.

Together with vector control, chemotherapy is an essential tool for the control of visceral leishmaniasis (VL), but its efficacy is jeopardized by growing resistance and treatment failure against first-line drugs. To delay the emergence of resistance, the use of drug combinations of existing antileishmanial agents has been tested systematically in clinical trials for the treatment of visceral leishmaniasis (VL). In vitro, Leishmania donovani promastigotes are able to develop experimental resistance to several combinations of different antileishmanial drugs after 10 weeks of drug pressure. Using an untargeted liquid chromatography-mass spectrometry (LC-MS) metabolomics approach, we identified metabolic changes in lines that were experimentally resistant to drug combinations and their respective single-resistant lines. This highlighted both collective metabolic changes (found in all combination therapy-resistant [CTR] lines) and specific ones (found in certain CTR lines). We demonstrated that single-resistant and CTR parasite cell lines show distinct metabolic adaptations, which all converge on the same defensive mechanisms that were experimentally validated: protection against drug-induced and external oxidative stress and changes in membrane fluidity. The membrane fluidity changes were accompanied by changes in drug uptake only in the lines that were resistant against drug combinations with antimonials, and surprisingly, drug accumulation was higher in these lines. Together, these results highlight the importance and the central role of protection against oxidative stress in the different resistant lines. Ultimately, these phenotypic changes might interfere with the mode of action of all drugs that are currently used for the treatment of VL and should be taken into account in drug development.

Metabolic adaptations of Leishmania donovani in relation to differentiation, drug resistance, and drug pressure.

Antimonial (sodium stibogluconate, SSG) resistance and differentiation have been shown to be closely linked in Leishmania donovani, with SSG-resistant strains showing an increased capacity to generate infectious (metacyclic) forms. This is the first untargeted LC-MS metabolomics study which integrated both phenomena in one experimental design and provided insights into metabolic differences between three clinical L. donovani strains with a similar genetic background but different SSG-susceptibilities. We performed this analysis at different stages during promastigote growth and in the absence or presence of drug pressure. When comparing SSG-resistant and SSG-sensitive strains, a number of metabolic changes appeared to be constitutively present in all growth stages, pointing towards a clear link with SSG-resistance, whereas most metabolic changes were only detected in the stationary stage. These changes reflect the close intertwinement between SSG-resistance and an increased metacyclogenesis in resistant parasites. The metabolic changes suggest that SSG-resistant parasites have (i) an increased capacity for protection against oxidative stress; (ii) a higher fluidity of the plasma membrane; and (iii) a metabolic survival kit to better endure infection. These changes were even more pronounced in a resistant strain kept under Sb(III) drug pressure.

Selective whole-genome sequencing of Plasmodium parasites directly from blood samples by nanopore adaptive sampling.

IMPORTANCE: Malaria is caused by parasites of the genus Plasmodium, and reached a global disease burden of 247 million cases in 2021. To study drug resistance mutations and parasite population dynamics, whole-genome sequencing of patient blood samples is commonly performed. However, the predominance of human DNA in these samples imposes the need for time-consuming laboratory procedures to enrich Plasmodium DNA. We used the Oxford Nanopore Technologies' adaptive sampling feature to circumvent this problem and enrich Plasmodium reads directly during the sequencing run. We demonstrate that adaptive nanopore sequencing efficiently enriches Plasmodium reads, which simplifies and shortens the timeline from blood collection to parasite sequencing. In addition, we show that the obtained data can be used for monitoring genetic markers, or to generate nearly complete genomes. Finally, owing to its inherent mobility, this technology can be easily applied on-site in endemic areas where patients would benefit the most from genomic surveillance.

Malate production, sugar metabolism, and redox homeostasis in the leaf growth zone of Rye (Secale cereale) increase stress tolerance to aluminum stress: A biochemical and genome-wide transcriptional study.

Global soil acidification is increasing, enlarging aluminum (Al) availability in soils, leading to reductions in plant growth. This study investigates the effect of Al stress on the leaf growth zones of Rye (Secale cereale, cv Beira). Kinematic analysis showed that the effect of Al on leaf growth rates was mainly due to a reduced cell production rate in the meristem. Transcriptomic analysis identified 2272 significantly (log2fold > |0.5| FDR < 0.05) differentially expressed genes (DEGs) for Al stress. There was a downregulation in several DEGs associated with photosynthetic processes and an upregulation in genes for heat/light response, and H2O2 production in all leaf zones. DEGs associated with heavy metals and malate transport were increased, particularly, in the meristem. To determine the putative function of these processes in Al tolerance, we performed biochemical analyses comparing the tolerant Beira with an Al sensitive variant RioDeva. Beira showed improved sugar metabolism and redox homeostasis, specifically in the meristem compared to RioDeva. Similarly, a significant increase in malate and citrate production, which are known to aid in Al detoxification in plants, was found in Beira. This suggests that Al tolerance in Rye is linked to its ability for Al exclusion from the leaf meristem.

Persistent T cell unresponsiveness associated with chronic visceral leishmaniasis in HIV-coinfected patients.

A large proportion of HIV-coinfected visceral leishmaniasis (VL-HIV) patients exhibit chronic disease with frequent VL recurrence. However, knowledge on immunological determinants underlying the disease course is scarce. We longitudinally profiled the circulatory cellular immunity of an Ethiopian HIV cohort that included VL developers. We show that chronic VL-HIV patients exhibit high and persistent levels of TIGIT and PD-1 on CD8+/CD8- T cells, in addition to a lower frequency of IFN-γ+ TIGIT- CD8+/CD8- T cells, suggestive of impaired T cell functionality. At single T cell transcriptome and clonal resolution, the patients show CD4+ T cell anergy, characterised by a lack of T cell activation and lymphoproliferative response. These findings suggest that PD-1 and TIGIT play a pivotal role in VL-HIV chronicity, and may be further explored for patient risk stratification. Our findings provide a strong rationale for adjunctive immunotherapy for the treatment of chronic VL-HIV patients to break the recurrent disease cycle.

A Lean Reverse Vaccinology Pipeline with Publicly Available Bioinformatic Tools.

Reverse vaccinology (RV) marked an outstanding improvement in vaccinology employing bioinformatics tools to extract effective features from protein sequences to drive the selection of potential vaccine candidates (Rappuoli, Curr Opin Microbiol 3(5):445-450, 2000). Pioneered by Rino Rappuoli and first used against serogroup B meningococcus, since then, it has been used on several other bacterial vaccines, varying during time the adopted bioinformatics tools. Based on our experience in the field of RV and following an extensive literature review, we consolidate a lean RV pipeline of publicly available bioinformatic tools whose usage is described in this contribution. The protein features, whose extraction is reported in this contribution, can be also the input in a matrix format for machine learning-based approaches.

IL-22-Activated MUC13 Impacts on Colonic Barrier Function through JAK1/STAT3, SNAI1/ZEB1 and ROCK2/MAPK Signaling.

Overexpression of the transmembrane mucin MUC13, as seen in inflammatory bowel diseases (IBD), could potentially impact barrier function. This study aimed to explore how inflammation-induced MUC13 disrupts epithelial barrier integrity by affecting junctional protein expression in IBD, thereby also considering the involvement of MUC1. RNA sequencing and permeability assays were performed using LS513 cells transfected with MUC1 and MUC13 siRNA and subsequently stimulated with IL-22. In vivo intestinal permeability and MUC13-related signaling pathways affecting barrier function were investigated in acute and chronic DSS-induced colitis wildtype and Muc13-/- mice. Finally, the expression of MUC13, its regulators and other barrier mediators were studied in IBD and control patients. Mucin knockdown in intestinal epithelial cells affected gene expression of several barrier mediators in the presence/absence of inflammation. IL-22-induced MUC13 expression impacted barrier function by modulating the JAK1/STAT3, SNAI1/ZEB1 and ROCK2/MAPK signaling pathways, with a cooperating role for MUC1. In response to DSS, MUC13 was protective during the acute phase whereas it caused more harm upon chronic colitis. The pathways accounting for the MUC13-mediated barrier dysfunction were also altered upon inflammation in IBD patients. These novel findings indicate an active role for aberrant MUC13 signaling inducing intestinal barrier dysfunction upon inflammation with MUC1 as collaborating partner.

ISCB Student Council Symposium 2021, a virtual global venue: challenges and lessons learned.

Since 2004, the ISCB Student Council has been organizing different symposia worldwide, gathering together the community of young computational biologists. Due to the coronavirus disease 2019 (COVID-19) pandemic situation, the world scientific community was forced to cancel in-person meetings for almost two years, imposing the adoption of virtual formats instead. After the successful editions of our continental symposia in 2020 in the USA, Latin America, and Europe, we organized our flagship global event, the Student Council Symposium (SCS) 2021, trying to apply all previous lessons learned and to exploit the advantages that virtuality has to offer.

Fractionated irradiation of MCF7 breast cancer cells rewires a gene regulatory circuit towards a treatment-resistant stemness phenotype.

Radiotherapy is the standard of care for breast cancer. However, surviving radioresistant cells can repopulate following treatment and provoke relapse. Better understanding of the molecular mechanisms of radiation resistance may help to improve treatment of radioresistant tumours. To emulate radiation therapy at the cellular level, we exposed MCF7 breast cancer cells to daily radiation doses of 2 Gy up to an accumulated dose of 20 Gy. Fractionally irradiated cells (FIR20) displayed increased clonogenic survival and population doubling time as compared with age-matched sham-irradiated cells and untreated parental MCF7 cells. RNA-sequencing revealed a core signature of 229 mRNAs and 7 circular RNAs of which the expression was significantly altered in FIR20 cells. Dysregulation of several top genes was mirrored at the protein level. The FIR20 cell transcriptome overlapped significantly with canonical radiation response signatures and demonstrated a remarkable commonality with radiation and endocrine therapy resistance expression profiles, suggesting crosstalk between both acquired resistance pathways, as indicated by reduced sensitivity to tamoxifen cytotoxicity of FIR20 cells. Using predictive analyses and functional enrichment, we identified a gene-regulatory network that promotes stemness and inflammatory signalling in FIR20 cells. We propose that these phenotypic traits render breast cancer cells more radioresistant but may at the same time serve as potential targets for combination therapies.

Band 3-mediated Plasmodium vivax invasion is associated with transcriptional variation in PvTRAg genes.

The Plasmodium vivax reticulocyte invasion process is still poorly understood, with only a few receptor-ligand interactions identified to date. Individuals with the Southeast Asian ovalocytosis (SAO) phenotype have a deletion in the band 3 protein on the surface of erythrocytes, and are reported to have a lower incidence of clinical P. vivax malaria. Based on this observation, band 3 has been put forward as a receptor for P. vivax invasion, although direct proof is still lacking. In this study, we combined functional ex vivo invasion assays and transcriptome sequencing to uncover a band 3-mediated invasion pathway in P. vivax and potential band 3 ligands. Invasion by P. vivax field isolates was 67%-71% lower in SAO reticulocytes compared with non-SAO reticulocytes. Reticulocyte invasion was decreased by 40% and 27%-31% when blocking with an anti-band 3 polyclonal antibody and a PvTRAg38 peptide, respectively. To identify new band 3 receptor candidates, we mRNA-sequenced schizont-stage isolates used in the invasion assays, and observed high transcriptional variability in multigene and invasion-related families. Transcriptomes of isolates with low or high dependency on band 3 for invasion were compared by differential expression analysis, which produced a list of band 3 ligand candidates with high representation of PvTRAg genes. Our ex vivo invasion assays have demonstrated that band 3 is a P. vivax invasion receptor and confirm previous in vitro studies showing binding between PvTRAg38 and band 3, although the lower and variable inhibition levels observed suggest the involvement of other ligands. By coupling transcriptomes and invasion phenotypes from the same isolates, we identified a list of band 3 ligand candidates, of which the overrepresented PvTRAg genes are the most promising for future research.