Search

Background & Aims: Polyploidy in hepatocytes has been proposed as a genetic mechanism to buffer against transcriptional dysregulation. Here, we aim to demonstrate the role of polyploidy in modulating gene regulatory networks in hepatocytes during ageing., Methods: We performed single-nucleus RNA sequencing in hepatocyte nuclei of different ploidy levels isolated from young and old wild-type mice. Changes in the gene expression and regulatory network were compared to three independent strains that were haploinsufficient for HNF4A, CEBPA or CTCF, representing non-deleterious perturbations. Phenotypic characteristics of the liver section were additionally evaluated histologically, whereas the genomic allele composition of hepatocytes was analysed by BaseScope., Results: We observed that ageing in wild-type mice results in nuclei polyploidy and a marked increase in steatosis. Haploinsufficiency of liver-specific master regulators (HFN4A or CEBPA) results in the enrichment of hepatocytes with tetraploid nuclei at a young age, affecting the genomic regulatory network, and dramatically suppressing ageing-related steatosis tissue wide. Notably, these phenotypes are not the result of subtle disruption to liver-specific transcriptional networks, since haploinsufficiency in the CTCF insulator protein resulted in the same phenotype. Further quantification of genotypes of tetraploid hepatocytes in young and old HFN4A-haploinsufficient mice revealed that during ageing, tetraploid hepatocytes lead to the selection of wild-type alleles, restoring non-deleterious genetic perturbations., Conclusions: Our results suggest a model whereby polyploidisation leads to fundamentally different cell states. Polyploid conversion enables pleiotropic buffering against age-related decline via non-random allelic segregation to restore a wild-type genome., Impact and Implications: The functional role of hepatocyte polyploidisation during ageing is poorly understood. Using single-nucleus RNA sequencing and BaseScope approaches, we have studied ploidy dynamics during ageing in murine livers with non-deleterious genetic perturbations. We have identified that hepatocytes present different cellular states and the ability to buffer ageing-associated dysfunctions. Tetraploid nuclei exhibit robust transcriptional networks and are better adapted to genomically overcome perturbations. Novel therapeutic interventions aimed at attenuating age-related changes in tissue function could be exploited by manipulation of ploidy dynamics during chronic liver conditions., (Copyright © 2024 The Authors. Published by Elsevier B.V. All rights reserved.)

External controls (eControls) leverage historical data to create non-randomized control arms. The lack of randomization can result in confounding between the experimental and eControl cohorts. To balance potentially confounding variables between the cohorts, one of the proposed methods is to match on prognostic scores. Still, the performance of prognostic scores to construct eControls in oncology has not been analyzed yet. Using an electronic health record-derived de-identified database, we constructed eControls using one of three methods: ROPRO, a state-of-the-art prognostic score, or either a propensity score composed of five (5Vars) or 27 covariates (ROPROvars). We compared the performance of these methods in estimating the overall survival (OS) hazard ratio (HR) of 11 recent advanced non-small cell lung cancer. The ROPRO eControls had a lower OS HR error (median absolute deviation (MAD), 0.072, confidence interval (CI): 0.036-0.185), than the 5Vars (MAD 0.081, CI: 0.025-0.283) and ROPROvars eControls (MAD 0.087, CI: 0.054-0.383). Notably, the OS HR errors for all methods were even lower in the phase III studies. Moreover, the ROPRO eControl cohorts included, on average, more patients than the 5Vars (6.54%) and ROPROvars cohorts (11.7%). The eControls matched with the prognostic score reproduced the controls more reliably than propensity scores composed of the underlying variables. Additionally, prognostic scores could allow eControls to be built on many prognostic variables without a significant increase in the variability of the propensity score, which would decrease the number of matched patients., (© 2023 Roche Diagnostics GmbH. Clinical Pharmacology & Therapeutics published by Wiley Periodicals LLC on behalf of American Society for Clinical Pharmacology and Therapeutics.)

The function of a cell is defined by its intrinsic characteristics and its niche: the tissue microenvironment in which it dwells. Here we combine single-cell and spatial transcriptomics data to discover cellular niches within eight regions of the human heart. We map cells to microanatomical locations and integrate knowledge-based and unsupervised structural annotations. We also profile the cells of the human cardiac conduction system 1 . The results revealed their distinctive repertoire of ion channels, G-protein-coupled receptors (GPCRs) and regulatory networks, and implicated FOXP2 in the pacemaker phenotype. We show that the sinoatrial node is compartmentalized, with a core of pacemaker cells, fibroblasts and glial cells supporting glutamatergic signalling. Using a custom CellPhoneDB.org module, we identify trans-synaptic pacemaker cell interactions with glia. We introduce a druggable target prediction tool, drug2cell, which leverages single-cell profiles and drug-target interactions to provide mechanistic insights into the chronotropic effects of drugs, including GLP-1 analogues. In the epicardium, we show enrichment of both IgG + and IgA + plasma cells forming immune niches that may contribute to infection defence. Overall, we provide new clarity to cardiac electro-anatomy and immunology, and our suite of computational approaches can be applied to other tissues and organs., (© 2023. The Author(s).)

Single-cell technologies have transformed our understanding of human tissues. Yet, studies typically capture only a limited number of donors and disagree on cell type definitions. Integrating many single-cell datasets can address these limitations of individual studies and capture the variability present in the population. Here we present the integrated Human Lung Cell Atlas (HLCA), combining 49 datasets of the human respiratory system into a single atlas spanning over 2.4 million cells from 486 individuals. The HLCA presents a consensus cell type re-annotation with matching marker genes, including annotations of rare and previously undescribed cell types. Leveraging the number and diversity of individuals in the HLCA, we identify gene modules that are associated with demographic covariates such as age, sex and body mass index, as well as gene modules changing expression along the proximal-to-distal axis of the bronchial tree. Mapping new data to the HLCA enables rapid data annotation and interpretation. Using the HLCA as a reference for the study of disease, we identify shared cell states across multiple lung diseases, including SPP1 + profibrotic monocyte-derived macrophages in COVID-19, pulmonary fibrosis and lung carcinoma. Overall, the HLCA serves as an example for the development and use of large-scale, cross-dataset organ atlases within the Human Cell Atlas., (© 2023. The Author(s).)

Human infection with the intestinal nematode Strongyloides stercoralis is persistent unless effectively treated, and potentially fatal in immunosuppressed individuals. Epidemiological data are lacking, partially due to inadequate diagnosis. A rapid antigen detection test is a priority for population surveillance, validating cure after treatment, and for screening prior to immunosuppression. We used a targeted analysis of open access 'omics' data sets and used online predictors to identify S. stercoralis proteins that are predicted to be present in infected stool, Strongyloides-specific, and antigenic. Transcriptomic data from gut and non-gut dwelling life cycle stages of S. stercoralis revealed 328 proteins that are differentially expressed. Strongyloides ratti proteomic data for excreted and secreted (E/S) proteins were matched to S. stercoralis, giving 1,057 orthologues. Five parasitism-associated protein families (SCP/TAPS, prolyl oligopeptidase, transthyretin-like, aspartic peptidase, acetylcholinesterase) were compared phylogenetically between S. stercoralis and outgroups, and proteins with least homology to the outgroups were selected. Proteins that overlapped between the transcriptomic and proteomic datasets were analysed by multiple sequence alignment, epitope prediction and 3D structure modelling to reveal S. stercoralis candidate peptide/protein coproantigens. We describe 22 candidates from seven genes, across all five protein families for further investigation as potential S. stercoralis diagnostic coproantigens, identified using open access data and freely-available protein analysis tools. This powerful approach can be applied to many parasitic infections with 'omic' data to accelerate development of specific diagnostic assays for laboratory or point-of-care field application., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Marlais et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

The increasing availability of large-scale single-cell atlases has enabled the detailed description of cell states. In parallel, advances in deep learning allow rapid analysis of newly generated query datasets by mapping them into reference atlases. However, existing data transformations learned to map query data are not easily explainable using biologically known concepts such as genes or pathways. Here we propose expiMap, a biologically informed deep-learning architecture that enables single-cell reference mapping. ExpiMap learns to map cells into biologically understandable components representing known 'gene programs'. The activity of each cell for a gene program is learned while simultaneously refining them and learning de novo programs. We show that expiMap compares favourably to existing methods while bringing an additional layer of interpretability to integrative single-cell analysis. Furthermore, we demonstrate its applicability to analyse single-cell perturbation responses in different tissues and species and resolve responses of patients who have coronavirus disease 2019 to different treatments across cell types., (© 2023. The Author(s).)

The liver is a complex and heterogenous tissue responsible for carrying out many critical physiological functions, such as the maintenance of energy homeostasis and the metabolism of xenobiotics, among others. These tasks are performed through tight coordination between hepatic parenchymal and non-parenchymal cells. Additionally, various metabolic activities are confined to specific areas of the hepatic lobule-a phenomenon called liver zonation. Recent advances in single-cell sequencing technologies have empowered researchers to investigate tissue heterogeneity at a single-cell resolution. In many complex tissues, including the liver, harsh enzymatic and/or mechanical dissociation protocols can negatively affect the viability or the quality of the single-cell suspensions needed to comprehensively characterize this organ in health and disease. This paper describes a robust and reproducible protocol for isolating nuclei from frozen, archived liver tissues. This method yields high-quality nuclei that are compatible with downstream, single-cell omics approaches, including single-nucleus RNA-seq, assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq), as well as multimodal omics (joint RNA-seq and ATAC-seq). This method has been successfully used for the isolation of nuclei from healthy and diseased human, mouse, and non-human primate frozen liver samples. This approach allows the unbiased isolation of all the major cell types in the liver and, therefore, offers a robust methodology for studying the liver at the single-cell resolution.

Protozoa and fungi are known to have extraordinarily diverse mechanisms of genetic exchange. However, the presence and epidemiological relevance of genetic exchange in Trypanosoma cruzi , the agent of Chagas disease, has been controversial and debated for many years. Field studies have identified both predominantly clonal and sexually recombining natural populations. Two of six natural T. cruzi lineages (TcV and TcVI) show hybrid mosaicism, using analysis of single-gene locus markers. The formation of hybrid strains in vitro has been achieved and this provides a framework to study the mechanisms and adaptive significance of genetic exchange. Using whole genome sequencing of a set of experimental hybrids strains, we have confirmed that hybrid formation initially results in tetraploid parasites. The hybrid progeny showed novel mutations that were not attributable to either (diploid) parent showing an increase in amino acid changes. In long-term culture, up to 800 generations, there was a variable but gradual erosion of progeny genomes towards triploidy, yet retention of elevated copy number was observed at several core housekeeping loci. Our findings indicate hybrid formation by fusion of diploid T. cruzi , followed by sporadic genome erosion, but with substantial potential for adaptive evolution, as has been described as a genetic feature of other organisms, such as some fungi., Competing Interests: GM, MY, EG, LM, BA The authors declare no competing interests, ML, CT, MM No competing interests declared, (© 2022, Matos et al.)

Trypanosoma cruzi , a zoonotic kinetoplastid protozoan parasite, is the causative agent of American trypanosomiasis (Chagas disease). Having a very plastic, repetitive and complex genome, the parasite displays a highly diverse repertoire of surface molecules, with pivotal roles in cell invasion, immune evasion and pathogenesis. Before 2016, the complexity of the genomic regions containing these genes impaired the assembly of a genome at chromosomal level, making it impossible to study the structure and function of the several thousand repetitive genes encoding the surface molecules of the parasite. We here describe the genome assembly of the Sylvio X10/1 genome sequence, which since 2016 has been used as a reference genome sequence for T. cruzi clade I (TcI), produced using high coverage PacBio single-molecule sequencing. It was used to analyze deep Illumina sequence data from 34 T. cruzi TcI isolates and clones from different geographic locations, sample sources and clinical outcomes. Resolution of the surface molecule gene distribution showed the unusual duality in the organization of the parasite genome, a synteny of the core genomic region with related protozoa flanked by unique and highly plastic multigene family clusters encoding surface antigens. The presence of abundant interspersed retrotransposons in these multigene family clusters suggests that these elements are involved in a recombination mechanism for the generation of antigenic variation and evasion of the host immune response on these TcI strains. The comparative genomic analysis of the cohort of TcI strains revealed multiple cases of such recombination events involving surface molecule genes and has provided new insights into T. cruzi population structure., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Talavera-López, Messenger, Lewis, Yeo, Reis-Cunha, Matos, Bartholomeu, Calzada, Saldaña, Ramírez, Guhl, Ocaña-Mayorga, Costales, Gorchakov, Jones, Nolan, Teixeira, Carrasco, Bottazzi, Hotez, Murray, Grijalva, Burleigh, Grisard, Miles and Andersson.)

Angiotensin-converting enzyme 2 (ACE2) and accessory proteases (TMPRSS2 and CTSL) are needed for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) cellular entry, and their expression may shed light on viral tropism and impact across the body. We assessed the cell-type-specific expression of ACE2, TMPRSS2 and CTSL across 107 single-cell RNA-sequencing studies from different tissues. ACE2, TMPRSS2 and CTSL are coexpressed in specific subsets of respiratory epithelial cells in the nasal passages, airways and alveoli, and in cells from other organs associated with coronavirus disease 2019 (COVID-19) transmission or pathology. We performed a meta-analysis of 31 lung single-cell RNA-sequencing studies with 1,320,896 cells from 377 nasal, airway and lung parenchyma samples from 228 individuals. This revealed cell-type-specific associations of age, sex and smoking with expression levels of ACE2, TMPRSS2 and CTSL. Expression of entry factors increased with age and in males, including in airway secretory cells and alveolar type 2 cells. Expression programs shared by ACE2 + TMPRSS2 + cells in nasal, lung and gut tissues included genes that may mediate viral entry, key immune functions and epithelial-macrophage cross-talk, such as genes involved in the interleukin-6, interleukin-1, tumor necrosis factor and complement pathways. Cell-type-specific expression patterns may contribute to the pathogenesis of COVID-19, and our work highlights putative molecular pathways for therapeutic intervention.

Cardiovascular disease is the leading cause of death worldwide. Advanced insights into disease mechanisms and therapeutic strategies require a deeper understanding of the molecular processes involved in the healthy heart. Knowledge of the full repertoire of cardiac cells and their gene expression profiles is a fundamental first step in this endeavour. Here, using state-of-the-art analyses of large-scale single-cell and single-nucleus transcriptomes, we characterize six anatomical adult heart regions. Our results highlight the cellular heterogeneity of cardiomyocytes, pericytes and fibroblasts, and reveal distinct atrial and ventricular subsets of cells with diverse developmental origins and specialized properties. We define the complexity of the cardiac vasculature and its changes along the arterio-venous axis. In the immune compartment, we identify cardiac-resident macrophages with inflammatory and protective transcriptional signatures. Furthermore, analyses of cell-to-cell interactions highlight different networks of macrophages, fibroblasts and cardiomyocytes between atria and ventricles that are distinct from those of skeletal muscle. Our human cardiac cell atlas improves our understanding of the human heart and provides a valuable reference for future studies.

We investigated SARS-CoV-2 potential tropism by surveying expression of viral entry-associated genes in single-cell RNA-sequencing data from multiple tissues from healthy human donors. We co-detected these transcripts in specific respiratory, corneal and intestinal epithelial cells, potentially explaining the high efficiency of SARS-CoV-2 transmission. These genes are co-expressed in nasal epithelial cells with genes involved in innate immunity, highlighting the cells' potential role in initial viral infection, spread and clearance. The study offers a useful resource for further lines of inquiry with valuable clinical samples from COVID-19 patients and we provide our data in a comprehensive, open and user-friendly fashion at www.covid19cellatlas.org.

Expression Atlas is EMBL-EBI's resource for gene and protein expression. It sources and compiles data on the abundance and localisation of RNA and proteins in various biological systems and contexts and provides open access to this data for the research community. With the increased availability of single cell RNA-Seq datasets in the public archives, we have now extended Expression Atlas with a new added-value service to display gene expression in single cells. Single Cell Expression Atlas was launched in 2018 and currently includes 123 single cell RNA-Seq studies from 12 species. The website can be searched by genes within or across species to reveal experiments, tissues and cell types where this gene is expressed or under which conditions it is a marker gene. Within each study, cells can be visualized using a pre-calculated t-SNE plot and can be coloured by different features or by cell clusters based on gene expression. Within each experiment, there are links to downloadable files, such as RNA quantification matrices, clustering results, reports on protocols and associated metadata, such as assigned cell types., (© The Author(s) 2019. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Although the spleen is broadly accepted as the major lymphoid organ involved in generating immune responses to the erythrocytic stages of the malaria parasite, Plasmodium, human splenic tissue is not readily available in most cases. As a result, most studies of malaria in humans rely on peripheral blood to assess cellular immune responses to malaria. The suitability of peripheral blood as a proxy for splenic immune responses is however unknown. Here, we have simultaneously analysed the transcriptomes of whole blood and spleen over 12 days of erythrocytic stage Plasmodium chabaudi infection in C57BL/6 mice. Using both unsupervised and directed approaches, we compared gene expression between blood and spleen over the course of infection. Taking advantage of publicly available datasets, we used machine learning approaches to infer cell proportions and cell-specific gene expression signatures from our whole tissue transcriptome data. Our findings demonstrate that spleen and blood are quite dissimilar, sharing only a small amount of transcriptional information between them, with transcriptional differences in both cellular composition and transcriptional activity. These results suggest that while blood transcriptome data may be useful in defining surrogate markers of protection and pathology, they should not be used to predict specific immune responses occurring in lymphoid organs.

The immune system encompasses a large degree of phenotypic diversity and plasticity in its cell types, and more is to be uncovered. We argue that large, multiomic datasets of single-cell resolution, in conjunction with improved computational methods, will be essential to resolving immune cell identity. Existing datasets, combined with 'big data' methodologies, can serve as a platform to support future studies in immunology. Technical and analytical advances in multiomics and spatial integration can provide a reference for gene regulation and cellular interactions in spatially structured tissue contexts. We posit that these developments may allow guided functional studies of immune cell populations and lay the groundwork for informed cell engineering and precision medicine., (Copyright © 2019 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Human lungs enable efficient gas exchange and form an interface with the environment, which depends on mucosal immunity for protection against infectious agents. Tightly controlled interactions between structural and immune cells are required to maintain lung homeostasis. Here, we use single-cell transcriptomics to chart the cellular landscape of upper and lower airways and lung parenchyma in healthy lungs, and lower airways in asthmatic lungs. We report location-dependent airway epithelial cell states and a novel subset of tissue-resident memory T cells. In the lower airways of patients with asthma, mucous cell hyperplasia is shown to stem from a novel mucous ciliated cell state, as well as goblet cell hyperplasia. We report the presence of pathogenic effector type 2 helper T cells (T H 2) in asthmatic lungs and find evidence for type 2 cytokines in maintaining the altered epithelial cell states. Unbiased analysis of cell-cell interactions identifies a shift from airway structural cell communication in healthy lungs to a T H 2-dominated interactome in asthmatic lungs.

Following publication of the original article [1], an error was reported in Table 1. The data repository links in the 4th column were incorrect. In this Correction, the corrected version of Table 1 is shown. The original publication of this article has been corrected.

Objectives: Plasmodium berghei ANKA infection in mice is a model for human cerebral malaria, the most severe complication of Plasmodium falciparum infection. Responses of brain microglia have been little investigated, and may contribute to the pathogenesis of cerebral malaria. We showed previously that microglia are activated in P. berghei infections, and that Type 1-Interferon signaling is important for activation. This dataset compares transcriptomic profiles of brain microglia of infected mice in the presence and absence of Type 1 interferon signaling, with the aim of identifying genes in microglia in this pathway during experimental cerebral malaria., Data Description: We documented global gene expression from microglial RNA from uninfected and P berghei-infected wild-type C57BL/6 and IFNA Receptor Knock-out mice using Illumina Beadarrays. Principal component analysis showed 4 groups of samples corresponding to naïve wild-type, naïve IFNA Receptor knock-out, infected wild-type, and IFNA Receptor knock-out mice. Differentially-expressed genes of microglia from the two groups of infected mice are documented. Gene set enrichment analysis showing the top 500 genes assigned to Reactome pathways from infected IFNA Receptor knock-out versus naïve, and infected WT versus naïve has been generated. These data will be useful for those interested in microglia cells, and in experimental cerebral malaria.

Salamanders exhibit an extraordinary ability among vertebrates to regenerate complex body parts. However, scarce genomic resources have limited our understanding of regeneration in adult salamanders. Here, we present the ~20 Gb genome and transcriptome of the Iberian ribbed newt Pleurodeles waltl, a tractable species suitable for laboratory research. We find that embryonic stem cell-specific miRNAs mir-93b and mir-427/430/302, as well as Harbinger DNA transposons carrying the Myb-like proto-oncogene have expanded dramatically in the Pleurodeles waltl genome and are co-expressed during limb regeneration. Moreover, we find that a family of salamander methyltransferases is expressed specifically in adult appendages. Using CRISPR/Cas9 technology to perturb transcription factors, we demonstrate that, unlike the axolotl, Pax3 is present and necessary for development and that contrary to mammals, muscle regeneration is normal without functional Pax7 gene. Our data provide a foundation for comparative genomic studies that generate models for the uneven distribution of regenerative capacities among vertebrates.

Cerebral malaria is a pathology involving inflammation in the brain. There are many immune cell types activated during this process, but there is little information on the response of microglia, in this severe complication. We examined microglia by genome wide transcriptomic analysis in a model of experimental cerebral malaria (ECM), in which C57BL/6 mice are infected with Plasmodium berghei ANKA. Thousands of transcripts were differentially expressed in microglia at two different time points during infection. Proliferation of microglia was a dominant feature before the onset of ECM, and supporting this, we observed an increase in numbers of these cells in the brain. When cerebral malaria symptoms were manifest, genes involved in immune responses and chemokine production were upregulated, which were possibly driven by Type I Interferon. Consistent with this hypothesis, in vitro culture of a microglial cell line with Interferon-β, but not infected red blood cells, resulted in production of several of the chemokines shown to be upregulated in the gene expression analysis. It appears that these responses are associated with ECM, as microglia from mice infected with a mutant P. berghei parasite (ΔDPAP3), which does not cause ECM, did not show the same level of activation or proliferation.

Kinetoplastids are an ancestral group of protists that contains free-living species and parasites with distinct mechanisms in response to stress. Here, we compared genes involved in antioxidant defense (AD), proposing an evolution model among trypanosomatids. All genes were identified in Bodo saltans , suggesting that AD mechanisms have evolved prior to adaptation for parasitic lifestyles. While most of the monoxenous and dixenous parasites revealed minor differences from B. saltans , the endosymbiont-bearing species have an increased number of genes. The absence of these genes was mainly observed in the extracellular parasites of the genera Phytomonas and Trypanosoma . In trypanosomes, a distinction was observed between stercorarian and salivarian parasites, except for Trypanosoma rangeli . Our analyses indicate that the variability of AD among trypanosomatids at the genomic level is not solely due to the geographical isolation, being mainly related to specific adaptations of their distinct biological cycles within insect vectors and to a parasitism of a wide range of hosts., Competing Interests: Authors disclose no potential conflicts of interest.

Trypanosoma cruzi affects millions of people worldwide. Clinical variability of Chagas disease can be due to the genetic variability of this parasite, requiring further genome studies. Here we report the genome sequence of the T. cruzi Dm28c clone (TcI), a strain related to the sylvatic cycle of the parasite.

Conifers have dominated forests for more than 200 million years and are of huge ecological and economic importance. Here we present the draft assembly of the 20-gigabase genome of Norway spruce (Picea abies), the first available for any gymnosperm. The number of well-supported genes (28,354) is similar to the >100 times smaller genome of Arabidopsis thaliana, and there is no evidence of a recent whole-genome duplication in the gymnosperm lineage. Instead, the large genome size seems to result from the slow and steady accumulation of a diverse set of long-terminal repeat transposable elements, possibly owing to the lack of an efficient elimination mechanism. Comparative sequencing of Pinus sylvestris, Abies sibirica, Juniperus communis, Taxus baccata and Gnetum gnemon reveals that the transposable element diversity is shared among extant conifers. Expression of 24-nucleotide small RNAs, previously implicated in transposable element silencing, is tissue-specific and much lower than in other plants. We further identify numerous long (>10,000 base pairs) introns, gene-like fragments, uncharacterized long non-coding RNAs and short RNAs. This opens up new genomic avenues for conifer forestry and breeding.

Background: Trypanosoma cruzi marinkellei is a bat-associated parasite of the subgenus Schizotrypanum and it is regarded as a T. cruzi subspecies. Here we report a draft genome sequence of T. c. marinkellei and comparison with T. c. cruzi. Our aims were to identify unique sequences and genomic features, which may relate to their distinct niches., Results: The T. c. marinkellei genome was found to be ~11% smaller than that of the human-derived parasite T. c. cruzi Sylvio X10. The genome size difference was attributed to copy number variation of coding and non-coding sequences. The sequence divergence in coding regions was ~7.5% between T. c. marinkellei and T. c. cruzi Sylvio X10. A unique acetyltransferase gene was identified in T. c. marinkellei, representing an example of a horizontal gene transfer from eukaryote to eukaryote. Six of eight examined gene families were expanded in T. c. cruzi Sylvio X10. The DGF gene family was expanded in T. c. marinkellei. T. c. cruzi Sylvio X10 contained ~1.5 fold more sequences related to VIPER and L1Tc elements. Experimental infections of mammalian cell lines indicated that T. c. marinkellei has the capacity to invade non-bat cells and undergo intracellular replication., Conclusions: Several unique sequences were identified in the comparison, including a potential subspecies-specific gene acquisition in T. c. marinkellei. The identified differences reflect the distinct evolutionary trajectories of these parasites and represent targets for functional investigation.