rKOMICS:An R package for processing mitochondrial minicircle assemblies in population-scale genome projects

Abstract Background The advent of population-scale genome projects has revolutionized our biological understanding of parasitic protozoa. However, while hundreds to thousands of nuclear genomes of parasitic protozoa have been generated and analyzed, information about the diversity, structure and evolution of their mitochondrial genomes remains fragmentary, mainly because of their extraordinary complexity. Indeed, unicellular flagellates of the order Kinetoplastida contain structurally the most complex mitochondrial genome of all eukaryotes, organized as a giant network of homogeneous maxicircles and heterogeneous minicircles. We recently developed KOMICS, an analysis toolkit that automates the assembly and circularization of the mitochondrial genomes of Kinetoplastid parasites. While this tool overcomes the limitation of extracting mitochondrial assemblies from Next-Generation Sequencing datasets, interpreting and visualizing the genetic (dis)similarity within and between samples remains a time-consuming process. Results Here, we present a new analysis toolkit—rKOMICS—to streamline the analyses of minicircle sequence diversity in population-scale genome projects. rKOMICS is a user-friendly R package that has simple installation requirements and that is applicable to all 27 trypanosomatid genera. Once minicircle sequence alignments are generated, rKOMICS allows to examine, summarize and visualize minicircle sequence diversity within and between samples through the analyses of minicircle sequence clusters. We showcase the functionalities of the (r)KOMICS tool suite using a whole-genome sequencing dataset from a recently published study on the history of diversification of the Leishmania braziliensis species complex in Peru. Analyses of population diversity and structure highlighted differences in minicircle sequence richness and composition between Leishmania subspecies, and between subpopulations within subspecies. Conclusion The rKOMICS package establishes a critical framework to manipulate, explore and extract biologically relevant information from mitochondrial minicircle assemblies in tens to hundreds of samples simultaneously and efficiently. This should facilitate research that aims to develop new molecular markers for identifying species-specific minicircles, or to study the ancestry of parasites for complementary insights into their evolutionary history

Unveiling drug-tolerant and persister-like cells in Leishmania braziliensis lines derived from patients with cutaneous leishmaniasis

IntroductionResistance against anti-Leishmania drugs (DR) has been studied for years, giving important insights into long-term adaptations of these parasites to drugs, through genetic modifications. However, microorganisms can also survive lethal drug exposure by entering into temporary quiescence, a phenomenon called drug tolerance (DT), which is rather unexplored in Leishmania.MethodsWe studied a panel of nine Leishmania braziliensis strains highly susceptible to potassium antimonyl tartrate (PAT), exposed promastigotes to lethal PAT pressure, and compared several cellular and molecular parameters distinguishing DT from DR.Results and discussionWe demonstrated in vitro that a variable proportion of cells remained viable, showing all the criteria of DT and not of DR: i) signatures of quiescence, under drug pressure: reduced proliferation and significant decrease of rDNA transcription; ii) reversibility of the phenotype: return to low IC50 after removal of drug pressure; and iii) absence of significant genetic differences between exposed and unexposed lineages of each strain and absence of reported markers of DR. We found different levels of quiescence and DT among the different L. braziliensis strains. We provide here a new in-vitro model of drug-induced quiescence and DT in Leishmania. Research should be extended in vivo, but the current model could be further exploited to support R&D, for instance, to guide the screening of compounds to overcome the quiescence resilience of the parasite, thereby improving the therapy of leishmaniasis

Analysis of cell hyperplasia and parietal cell dysfunction induced by Ostertagia ostertagi infection

Infections in cattle with the gastric nematode Ostertagia ostertagi are associated with decreased acid secretion and profound physio-morphological changes of the gastric mucosa. The purpose of the current study was to investigate the mechanisms triggering these pathophysiological changes. O. ostertagi infection resulted in a marked cellular hyperplasia, which can be explained by increased transcriptional levels of signaling molecules related to the homeostasis of gastric epithelial cells such as HES1, WNT5A, FGF10, HB-EGF, AREG, ADAM10 and ADAM17. Intriguingly, histological analysis indicated that the rapid rise in the gastric pH, observed following the emergence of adult worms, cannot be explained by a loss of parietal cells, as a decrease in the number of parietal cells was only observed following a long term infection of several weeks, but is likely to be caused by an inhibition of parietal cell activity. To investigate whether this inhibition is caused by a direct effect of the parasites, parietal cells were co-cultured with parasite Excretory/Secretory products (ESP) and subsequently analyzed for acid production. The results indicate that adult ESP inhibited acid secretion, whereas ESP from the L4 larval stages did not alter parietal cell function. In addition, our data show that the inhibition of parietal cell activity could be mediated by a marked upregulation of inflammatory factors, which are partly induced by adult ESP in abomasal epithelial cells. In conclusion, this study shows that the emergence of adult O. ostertagi worms is associated with marked cellular changes that can be partly triggered by the worm's Excretory/secretory antigens

Comparative immune responses against Psoroptes ovis in two cattle breeds with different susceptibility to mange

The sheep scab mite, Psoroptes ovis, is a major problem in the beef cattle industry, especially in Belgian Blue (BB) cattle. This breed is naturally more predisposed to psoroptic mange but reasons for this high susceptibility remain unknown. Different immune responses could be a potential cause; thus in this study, the cutaneous immune response and in vitro cellular immune response after antigen re-stimulation were examined in naturally infested BB. Cytokine production in the skin and in circulating re-stimulated peripheral blood mononuclear cells (PBMC) demonstrated a mixed pro-inflammatory Th2/Th17 profile, with transcription of IL-4, IL-13, IL-6 and IL-17. Strong IL-17 up-regulation in the skin of BB was associated with an influx of eosinophils and other immune cells, potentially leading towards more severe symptoms. Virtually no changes in cutaneous IFN-gamma transcription were detected, while there was substantial IFN-gamma up-regulation in re-stimulated PBMC from infested and uninfested animals, potentially indicating a role of this pro-inflammatory cytokine in the innate immune response. In Holstein-Friesian (HF) cattle, generally more resistant to P.ovis infection, a largely similar immunologic response was observed. Differences between HF and BB were the lack of cutaneous IL-17 response in infested HF and low transcription levels of IFN-gamma and high IL-10 transcription in re-stimulated PBMC from both infested and uninfested animals. Further research is needed to identify potential cell sources and biological functions for these cytokines and to fully unravel the basis of this different breed susceptibility to P. ovis

Genome diversity of Leishmania aethiopica

Leishmania aethiopica is a zoonotic Old World parasite transmitted by Phlebotomine sand flies and causing cutaneous leishmaniasis in Ethiopia and Kenya. Despite a range of clinical manifestations and a high prevalence of treatment failure, L. aethiopica is one of the most neglected species of the Leishmania genus in terms of scientific attention. Here, we explored the genome diversity of L. aethiopica by analyzing the genomes of twenty isolates from Ethiopia. Phylogenomic analyses identified two strains as interspecific hybrids involving L. aethiopica as one parent and L. donovani and L. tropica respectively as the other parent. High levels of genome-wide heterozygosity suggest that these two hybrids are equivalent to F1 progeny that propagated mitotically since the initial hybridization event. Analyses of allelic read depths further revealed that the L. aethiopica - L. tropica hybrid was diploid and the L. aethiopica - L. donovani hybrid was triploid, as has been described for other interspecific Leishmania hybrids. When focusing on L. aethiopica, we show that this species is genetically highly diverse and consists of both asexually evolving strains and groups of recombining parasites. A remarkable observation is that some L. aethiopica strains showed an extensive loss of heterozygosity across large regions of the nuclear genome, which likely arose from gene conversion/mitotic recombination. Hence, our prospection of L. aethiopica genomics revealed new insights into the genomic consequences of both meiotic and mitotic recombination in Leishmania

Deep kinetoplast genome analyses result in a novel molecular assay for detecting trypanosoma brucei gambiense-specific minicircles

The World Health Organization targeted Trypanosoma brucei gambiense (Tbg) human African trypanosomiasis for elimination of transmission by 2030. Sensitive molecular markers that specifically detect Tbg type 1 (Tbg1) parasites will be important tools to assist in reaching this goal. We aim at improving molecular diagnosis of Tbg1 infections by targeting the abundant mitochondrial minicircles within the kinetoplast of these parasites. Using Next-Generation Sequencing of total cellular DNA extracts, we assembled and annotated the kinetoplast genome and investigated minicircle sequence diversity in 38 animal- and human-infective trypanosome strains. Computational analyses recognized a total of 241 Minicircle Sequence Classes as Tbg1-specific, of which three were shared by the 18 studied Tbg1 strains. We developed a minicircle-based assay that is applicable on animals and as specific as the TgsGP-based assay, the current golden standard for molecular detection of Tbg1. The median copy number of the targeted minicircle was equal to eight, suggesting our minicircle-based assay may be used for the sensitive detection of Tbg1 parasites. Annotation of the targeted minicircle sequence indicated that it encodes genes essential for the survival of the parasite and will thus likely be preserved in natural Tbg1 populations, the latter ensuring the reliability of our novel diagnostic assay

Genomic population structure associated with repeated escape of Salmonella enterica ATCC14028s from the laboratory into nature

Salmonella enterica serovar Typhimurium strain ATCC14028s is commercially available from multiple national type culture collections, and has been widely used since 1960 for quality control of growth media and experiments on fitness (“laboratory evolution”). ATCC14028s has been implicated in multiple cross-contaminations in the laboratory, and has also caused multiple laboratory infections and one known attempt at bioterrorism. According to hierarchical clustering of 3002 core gene sequences, ATCC14028s belongs to HierCC cluster HC20_373 in which most internal branch lengths are only one to three SNPs long. Many natural Typhimurium isolates from humans, domesticated animals and the environment also belong to HC20_373, and their core genomes are almost indistinguishable from those of laboratory strains. These natural isolates have infected humans in Ireland and Taiwan for decades, and are common in the British Isles as well as the Americas. The isolation history of some of the natural isolates confirms the conclusion that they do not represent recent contamination by the laboratory strain, and 10% carry plasmids or bacteriophages which have been acquired in nature by HGT from unrelated bacteria. We propose that ATCC14028s has repeatedly escaped from the laboratory environment into nature via laboratory accidents or infections, but the escaped micro-lineages have only a limited life span. As a result, there is a genetic gap separating HC20_373 from its closest natural relatives due to a divergence between them in the late 19th century followed by repeated extinction events of escaped HC20_373

Genome analysis of triploid hybrid Leishmania parasite from the Neotropics

We discovered a hybrid Leishmania parasite in Costa Rica that is genetically similar to hybrids from Panama. Genome analyses demonstrated the hybrid is triploid and identified L. braziliensis and L. guyanensis-related strains as parents. Our findings highlight the existence of poorly sampled Leishmania (Viannia) variants infectious to humans

Reconciling biodiversity and carbon stock conservation in an Afrotropical forest landscape

Protecting aboveground carbon stocks in tropical forests is essential for mitigating global climate change and is assumed to simultaneously conserve biodiversity. Although the relationship between tree diversity and carbon stocks is generally positive, the relationship remains unclear for consumers or decomposers. We assessed this relationship for multiple trophic levels across the tree of life (10 organismal groups, 3 kingdoms) in lowland rainforests of the Congo Basin. Comparisons across regrowth and old-growth forests evinced the expected positive relationship for trees, but not for other organismal groups. Moreover, differences in species composition between forests increased with difference in carbon stock. These variable associations across the tree of life contradict the implicit assumption that maximum co-benefits to biodiversity are associated with conservation of forests with the highest carbon storage. Initiatives targeting climate change mitigation and biodiversity conservation should include both old-growth and regenerating forests to optimally benefit biodiversity and carbon storage

The phylodynamics of SARS-CoV-2 during 2020 in Finland

Finland has had a low incidence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV2) infections as compared to most European countries. Here we report the origins and turnover of SARS-CoV-2 lineages circulating in Finland in 2020. SARS-CoV-2 introduced to Finland in January 2020 and spread rapidly across southern Finland during spring. We observed rapid turnover among Finnish lineages during this period. Clade 20C became the most prevalent among sequenced cases and was replaced by other strains in fall 2020. Bayesian phylogeographic reconstructions suggested 42 independent introductions into Finland during spring 2020, mainly from Italy, Austria, and Spain, which might have been the source for a third of cases. The investigations of the original introductions of SARS-CoV-2 to Finland during the early stages of the pandemic and of the subsequent lineage dynamics could be utilized to assess the role of transboundary movements and effects of early intervention and public health measures.Peer reviewe