Genomes of Endotrypanum monterogeii from Panama and Zelonia costaricensis from Brazil: Expansion of Multigene Families in Leishmaniinae Parasites That Are Close Relatives of Leishmania spp.

The Leishmaniinae subfamily of the Trypanosomatidae contains both genus Zelonia (monoxenous) and Endotrypanum (dixenous). They are amongst the nearest known relatives of Leishmania, which comprises many human pathogens widespread in the developing world. These closely related lineages are models for the genomic biology of monoxenous and dixenous parasites. Herein, we used comparative genomics to identify the orthologous groups (OGs) shared among 26 Leishmaniinae species to investigate gene family expansion/contraction and applied two phylogenomic approaches to confirm relationships within the subfamily. The Endotrypanum monterogeii and Zelonia costaricensis genomes were assembled, with sizes of 29.9 Mb and 38.0 Mb and 9.711 and 12.201 predicted protein-coding genes, respectively. The genome of E. monterogeii displayed a higher number of multicopy cell surface protein families, including glycoprotein 63 and glycoprotein 46, compared to Leishmania spp. The genome of Z. costaricensis presents expansions of BT1 and amino acid transporters and proteins containing leucine-rich repeat domains, as well as a loss of ABC-type transporters. In total, 415 and 85 lineage-specific OGs were identified in Z. costaricensis and E. monterogeii. The evolutionary relationships within the subfamily were confirmed using the supermatrix (3384 protein-coding genes) and supertree methods. Overall, this study showed new expansions of multigene families in monoxenous and dixenous parasites of the subfamily Leishmaniinae.

Complete genome sequence of the archetype bile acid 7α-dehydroxylating bacterium, Clostridium scindens VPI12708, isolated from human feces, circa 1980.

Clostridium scindens strain VPI12708 serves as model organism to study bile acid 7α-dehydroxylating pathways. The closed circular genome of C. scindens VPI12708 was obtained by PacBio sequencing. The genome is composed of 3,983,052 bp, with 47.59% G + C, and 3,707 coding DNA sequences are predicted.

SIEVE: joint inference of single-nucleotide variants and cell phylogeny from single-cell DNA sequencing data.

We present SIEVE, a statistical method for the joint inference of somatic variants and cell phylogeny under the finite-sites assumption from single-cell DNA sequencing. SIEVE leverages raw read counts for all nucleotides and corrects the acquisition bias of branch lengths. In our simulations, SIEVE outperforms other methods in phylogenetic reconstruction and variant calling accuracy, especially in the inference of homozygous variants. Applying SIEVE to three datasets, one for triple-negative breast (TNBC), and two for colorectal cancer (CRC), we find that double mutant genotypes are rare in CRC but unexpectedly frequent in the TNBC samples.

Formation of secondary allo-bile acids by novel enzymes from gut Firmicutes.

The gut microbiome of vertebrates is capable of numerous biotransformations of bile acids, which are responsible for intestinal lipid digestion and function as key nutrient-signaling molecules. The human liver produces bile acids from cholesterol predominantly in the A/B-cis orientation in which the sterol rings are "kinked", as well as small quantities of A/B-trans oriented "flat" stereoisomers known as "primary allo-bile acids". While the complex multi-step bile acid 7α-dehydroxylation pathway has been well-studied for conversion of "kinked" primary bile acids such as cholic acid (CA) and chenodeoxycholic acid (CDCA) to deoxycholic acid (DCA) and lithocholic acid (LCA), respectively, the enzymatic basis for the formation of "flat" stereoisomers allo-deoxycholic acid (allo-DCA) and allo-lithocholic acid (allo-LCA) by Firmicutes has remained unsolved for three decades. Here, we present a novel mechanism by which Firmicutes generate the "flat" bile acids allo-DCA and allo-LCA. The BaiA1 was shown to catalyze the final reduction from 3-oxo-allo-DCA to allo-DCA and 3-oxo-allo-LCA to allo-LCA. Phylogenetic and metagenomic analyses of human stool samples indicate that BaiP and BaiJ are encoded only in Firmicutes and differ from membrane-associated bile acid 5α-reductases recently reported in Bacteroidetes that indirectly generate allo-LCA from 3-oxo-Δ4-LCA. We further map the distribution of baiP and baiJ among Firmicutes in human metagenomes, demonstrating an increased abundance of the two genes in colorectal cancer (CRC) patients relative to healthy individuals.

Comparative analysis of capture methods for genomic profiling of circulating tumor cells in colorectal cancer.

The genomic profiling of circulating tumor cells (CTCs) in the bloodstream should provide clinically relevant information on therapeutic efficacy and help predict cancer survival. Here, we contrasted the genomic profiles of CTC pools recovered from metastatic colorectal cancer (mCRC) patients using different enrichment strategies (CellSearch, Parsortix, and FACS). Mutations inferred in the CTC pools differed depending on the enrichment strategy and, in all cases, represented a subset of the mutations detected in the matched primary tumor samples. However, the CTC pools from Parsortix, and in part, CellSearch, showed diversity estimates, mutational signatures, and drug-suitability scores remarkably close to those found in matching primary tumor samples. In addition, FACS CTC pools were enriched in apparent sequencing artifacts, leading to much higher genomic diversity estimates. Our results highlight the utility of CTCs to assess the genomic heterogeneity of individual tumors and help clinicians prioritize drugs in mCRC.

Clonality and timing of relapsing colorectal cancer metastasis revealed through whole-genome single-cell sequencing.

Recurrence of tumor cells following local and systemic therapy is a significant hurdle in cancer. Most patients with metastatic colorectal cancer (mCRC) will relapse, despite resection of the metastatic lesions. A better understanding of the evolutionary history of recurrent lesions is required to identify the spatial and temporal patterns of metastatic progression and expose the genetic and evolutionary determinants of therapeutic resistance. With this goal in mind, here we leveraged a unique single-cell whole-genome sequencing dataset from recurrent hepatic lesions of an mCRC patient. Our phylogenetic analysis confirms that the treatment induced a severe demographic bottleneck in the liver metastasis but also that a previously diverged lineage survived this surgery, possibly after migration to a different site in the liver. This lineage evolved very slowly for two years under adjuvant drug therapy and diversified again in a very short period. We identified several non-silent mutations specific to this lineage and inferred a substantial contribution of chemotherapy to the overall, genome-wide mutational burden. All in all, our study suggests that mCRC subclones can migrate locally and evade resection, keep evolving despite rounds of chemotherapy, and re-expand explosively.

Single-cell mtDNA heteroplasmy in colorectal cancer.

Human mitochondria can be genetically distinct within the same individual, a phenomenon known as heteroplasmy. In cancer, this phenomenon seems exacerbated, and most mitochondrial mutations seem to be heteroplasmic. How this genetic variation is arranged within and among normal and tumor cells is not well understood. To address this question, here we sequenced single-cell mitochondrial genomes from multiple normal and tumoral locations in four colorectal cancer patients. Our results suggest that single cells, both normal and tumoral, can carry various mitochondrial haplotypes. Remarkably, this intra-cell heteroplasmy can arise before tumor development and be maintained afterward in specific tumoral cell subpopulations. At least in the colorectal patients studied here, the somatic mutations in the single-cells do not seem to have a prominent role in tumorigenesis.

CellPhy: accurate and fast probabilistic inference of single-cell phylogenies from scDNA-seq data.

We introduce CellPhy, a maximum likelihood framework for inferring phylogenetic trees from somatic single-cell single-nucleotide variants. CellPhy leverages a finite-site Markov genotype model with 16 diploid states and considers amplification error and allelic dropout. We implement CellPhy into RAxML-NG, a widely used phylogenetic inference package that provides statistical confidence measurements and scales well on large datasets with hundreds or thousands of cells. Comprehensive simulations suggest that CellPhy is more robust to single-cell genomics errors and outperforms state-of-the-art methods under realistic scenarios, both in accuracy and speed. CellPhy is freely available at https://github.com/amkozlov/cellphy .

Endosymbiont Capture, a Repeated Process of Endosymbiont Transfer with Replacement in Trypanosomatids Angomonas spp.

Trypanosomatids of the subfamily Strigomonadinae bear permanent intracellular bacterial symbionts acquired by the common ancestor of these flagellates. However, the cospeciation pattern inherent to such relationships was revealed to be broken upon the description of Angomonas ambiguus, which is sister to A. desouzai, but bears an endosymbiont genetically close to that of A. deanei. Based on phylogenetic inferences, it was proposed that the bacterium from A. deanei had been horizontally transferred to A. ambiguus. Here, we sequenced the bacterial genomes from two A. ambiguus isolates, including a new one from Papua New Guinea, and compared them with the published genome of the A. deanei endosymbiont, revealing differences below the interspecific level. Our phylogenetic analyses confirmed that the endosymbionts of A. ambiguus were obtained from A. deanei and, in addition, demonstrated that this occurred more than once. We propose that coinfection of the same blowfly host and the phylogenetic relatedness of the trypanosomatids facilitate such transitions, whereas the drastic difference in the occurrence of the two trypanosomatid species determines the observed direction of this process. This phenomenon is analogous to organelle (mitochondrion/plastid) capture described in multicellular organisms and, thereafter, we name it endosymbiont capture.

Completion of the gut microbial epi-bile acid pathway.

Bile acids are detergent molecules that solubilize dietary lipids and lipid-soluble vitamins. Humans synthesize bile acids with α-orientation hydroxyl groups which can be biotransformed by gut microbiota to toxic, hydrophobic bile acids, such as deoxycholic acid (DCA). Gut microbiota can also convert hydroxyl groups from the α-orientation through an oxo-intermediate to the ß-orientation, resulting in more hydrophilic, less toxic bile acids. This interconversion is catalyzed by regio- (C-3 vs. C-7) and stereospecific (α vs. ß) hydroxysteroid dehydrogenases (HSDHs). So far, genes encoding the urso- (7α-HSDH & 7ß-HSDH) and iso- (3α-HSDH & 3ß-HSDH) bile acid pathways have been described. Recently, multiple human gut clostridia were reported to encode 12α-HSDH, which interconverts DCA and 12-oxolithocholic acid (12-oxoLCA). 12ß-HSDH completes the epi-bile acid pathway by converting 12-oxoLCA to the 12ß-bile acid denoted epiDCA; however, a gene(s) encoding this enzyme has yet to be identified. We confirmed 12ß-HSDH activity in cultures of Clostridium paraputrificum ATCC 25780. From six candidate C. paraputrificum ATCC 25780 oxidoreductase genes, we discovered the first gene (DR024_RS09610) encoding bile acid 12ß-HSDH. Phylogenetic analysis revealed unforeseen diversity for 12ß-HSDH, leading to validation of two additional bile acid 12ß-HSDHs through a synthetic biology approach. By comparison to a previous phylogenetic analysis of 12α-HSDH, we identified the first potential C-12 epimerizing strains: Collinsella tanakaei YIT 12063 and Collinsella stercoris DSM 13279. A Hidden Markov Model search against human gut metagenomes located putative 12ß-HSDH genes in about 30% of subjects within the cohorts analyzed, indicating this gene is relevant in the human gut microbiome.

The Importance of Glycerophospholipid Production to the Mutualist Symbiosis of Trypanosomatids.

The symbiosis in trypanosomatids is a mutualistic relationship characterized by extensive metabolic exchanges between the bacterium and the protozoan. The symbiotic bacterium can complete host essential metabolic pathways, such as those for heme, amino acid, and vitamin production. Experimental assays indicate that the symbiont acquires phospholipids from the host trypanosomatid, especially phosphatidylcholine, which is often present in bacteria that have a close association with eukaryotic cells. In this work, an in-silico study was performed to find genes involved in the glycerophospholipid (GPL) production of Symbiont Harboring Trypanosomatids (SHTs) and their respective bacteria, also extending the search for trypanosomatids that naturally do not have symbionts. Results showed that most genes for GPL synthesis are only present in the SHT. The bacterium has an exclusive sequence related to phosphatidylglycerol production and contains genes for phosphatidic acid production, which may enhance SHT phosphatidic acid production. Phylogenetic data did not indicate gene transfers from the bacterium to the SHT nucleus, proposing that enzymes participating in GPL route have eukaryotic characteristics. Taken together, our data indicate that, differently from other metabolic pathways described so far, the symbiont contributes little to the production of GPLs and acquires most of these molecules from the SHT.

Phylogeny and temporal diversification of mosquitoes (Diptera: Culicidae) with an emphasis on the Neotropical fauna

Understanding phylogenetic relationships within the family Culicidae informs mosquito evolution and may have public health implications as this family includes numerous species of medical and veterinary importance. We investigated the mitochondrial genomes of 102 mosquitoes, including six newly sequenced species, representing 21 genera with an emphasis on the Neotropical region. We estimated divergence times based on sequence data and three fossil calibration points, using Bayesian relaxed clock methods. Bayesian and maximum-likelihood phylogenetic analyses based on the DNA sequences of 13 PCGs of the 102 species provided robust support for the monophyly of the subfamily Anophelinae and the tribes Aedini, Culicini, Mansoniini and Sabethini. Despite the current genera of Anophelinae being consistently recovered as monophyletic, relationships among them proved to be quite variable depending on the method used (concatenated or partitioned) and the number of taxa sampled. Molecular divergence time estimates revealed that the two mosquito subfamilies, Anophelinae and Culicinae, diverged in the early Jurassic (approximately 197.5 Mya). However, most major lineages of these groups arose after the Cretaceous, coincident with the emergence of angiosperms and the expansion of mammals and birds. The diversification and worldwide distribution of Culicidae may also be determined in part by geographic isolation as a result of continental drift during the Cretaceous.

A collection of bacterial isolates from the pig intestine reveals functional and taxonomic diversity.

Our knowledge about the gut microbiota of pigs is still scarce, despite the importance of these animals for biomedical research and agriculture. Here, we present a collection of cultured bacteria from the pig gut, including 110 species across 40 families and nine phyla. We provide taxonomic descriptions for 22 novel species and 16 genera. Meta-analysis of 16S rRNA amplicon sequence data and metagenome-assembled genomes reveal prevalent and pig-specific species within Lactobacillus, Streptococcus, Clostridium, Desulfovibrio, Enterococcus, Fusobacterium, and several new genera described in this study. Potentially interesting functions discovered in these organisms include a fucosyltransferase encoded in the genome of the novel species Clostridium porci, and prevalent gene clusters for biosynthesis of sactipeptide-like peptides. Many strains deconjugate primary bile acids in in vitro assays, and a Clostridium scindens strain produces secondary bile acids via dehydroxylation. In addition, cells of the novel species Bullifex porci are coccoidal or spherical under the culture conditions tested, in contrast with the usual helical shape of other members of the family Spirochaetaceae. The strain collection, called 'Pig intestinal bacterial collection' (PiBAC), is publicly available at www.dsmz.de/pibac and opens new avenues for functional studies of the pig gut microbiota.

Bacterial steroid-17,20-desmolase is a taxonomically rare enzymatic pathway that converts prednisone to 1,4-androstanediene-3,11,17-trione, a metabolite that causes proliferation of prostate cancer cells.

The adrenal gland has traditionally been viewed as a source of "weak androgens"; however, emerging evidence indicates 11-oxy-androgens of adrenal origin are metabolized in peripheral tissues to potent androgens. Also emerging is the role of gut bacteria in the conversion of C21 glucocorticoids to 11-oxygenated C19 androgens. Clostridium scindens ATCC 35,704 is a gut microbe capable of converting cortisol into 11-oxy-androgens by cleaving the side-chain. The desA and desB genes encode steroid-17,20-desmolase. Our prior study indicated that the urinary tract bacterium, Propionimicrobium lymphophilum ACS-093-V-SCH5 encodes desAB and converts cortisol to 11ß-hydroxyandrostenedione. We wanted to determine how widespread this function occurs in the human microbiome. Phylogenetic and sequence similarity network analyses indicated that the steroid-17,20-desmolase pathway is taxonomically rare and located in gut and urogenital microbiomes. Two microbes from each of these niches, C. scindens and Propionimicrobium lymphophilum, respectively, were screened for activity against endogenous (cortisol, cortisone, and allotetrahydrocortisol) and exogenous (prednisone, prednisolone, dexamethasone, and 9-fluorocortisol) glucocorticoids. LC/MS analysis showed that both microbes were able to side-chain cleave all glucocorticoids, forming 11-oxy-androgens. Pure recombinant DesAB from C. scindens showed the highest activity against prednisone, a commonly prescribed glucocorticoid. In addition, 0.1 nM 1,4-androstadiene-3,11,17-trione, bacterial side-chain cleavage product of prednisone, showed significant proliferation relative to vehicle in androgen-dependent growth LNCaP prostate cancer cells after 24 h (2.3 fold; P <  0.01) and 72 h (1.6 fold; P < 0.01). Taken together, DesAB-expressing microbes may be an overlooked source of androgens in the body, potentially contributing to various disease states, such as prostate cancer.

Rapid evolution and biogeographic spread in a colorectal cancer.

How and when tumoral clones start spreading to surrounding and distant tissues is currently unclear. Here we leveraged a model-based evolutionary framework to investigate the demographic and biogeographic history of a colorectal cancer. Our analyses strongly support an early monoclonal metastatic colonization, followed by a rapid population expansion at both primary and secondary sites. Moreover, we infer a hematogenous metastatic spread under positive selection, plus the return of some tumoral cells from the liver back to the colon lymph nodes. This study illustrates how sophisticated techniques typical of organismal evolution can provide a detailed, quantitative picture of the complex tumoral dynamics over time and space.

Investigation of the pathways related to intrinsic miltefosine tolerance in Leishmania (Viannia) braziliensis clinical isolates reveals differences in drug uptake.

In Brazil, cutaneous leishmaniasis is caused predominantly by L. (V.) braziliensis. The few therapeutic drugs available exhibit several limitations, mainly related to drug toxicity and reduced efficacy in some regions. Miltefosine (MF), the only oral drug available for leishmaniasis treatment, is not widely available and has not yet been approved for human use in Brazil. Our group previously reported the existence of differential susceptibility among L. (V.) braziliensis clinical isolates. In this work, we further characterized three of these isolates of L. (V.) braziliensis chosen because they exhibited the lowest and the highest MF half maximal inhibitory concentrations and were therefore considered less tolerant or more tolerant, respectively. Uptake of MF, and also of phosphocholine, were found to be significantly different in more tolerant parasites compared to the less sensitive isolate, which raised the hypothesis of differences in the MF transport complex Miltefosine Transporter (MT)-Ros3. Although some polymorphisms in those genes were found, they did not correlate with the drug susceptibility phenotype. Drug efflux and compartmentalization were similar in the isolates tested, and amphotericin B susceptibility was retained in MF tolerant parasites, suggesting that increased fitness was also not the basis of observed differences. Transcriptomic analysis revealed that Ros3 mRNA levels were upregulated in the sensitive strain compared to the tolerant ones. Increased mRNA abundance in more tolerant isolates was validated by quantitative PCR. Our results suggest that differential gene expression of the MT transporter complex is the basis of the differential susceptibility in these unselected, naturally occurring parasites.

Clostridium scindens ATCC 35704: Integration of Nutritional Requirements, the Complete Genome Sequence, and Global Transcriptional Responses to Bile Acids.

In the human gut, Clostridium scindens ATCC 35704 is a predominant bacterium and one of the major bile acid 7α-dehydroxylating anaerobes. While this organism is well-studied relative to bile acid metabolism, little is known about the basic nutrition and physiology of C. scindens ATCC 35704. To determine the amino acid and vitamin requirements of C. scindens, the leave-one-out (one amino acid group or vitamin) technique was used to eliminate the nonessential amino acids and vitamins. With this approach, the amino acid tryptophan and three vitamins (riboflavin, pantothenate, and pyridoxal) were found to be required for the growth of C. scindens In the newly developed defined medium, C. scindens fermented glucose mainly to ethanol, acetate, formate, and H2. The genome of C. scindens ATCC 35704 was completed through PacBio sequencing. Pathway analysis of the genome sequence coupled with transcriptome sequencing (RNA-Seq) under defined culture conditions revealed consistency with the growth requirements and end products of glucose metabolism. Induction with bile acids revealed complex and differential responses to cholic acid and deoxycholic acid, including the expression of potentially novel bile acid-inducible genes involved in cholic acid metabolism. Responses to toxic deoxycholic acid included expression of genes predicted to be involved in DNA repair, oxidative stress, cell wall maintenance/metabolism, chaperone synthesis, and downregulation of one-third of the genome. These analyses provide valuable insight into the overall biology of C. scindens which may be important in treatment of disease associated with increased colonic secondary bile acids.IMPORTANCEC. scindens is one of a few identified gut bacterial species capable of converting host cholic acid into disease-associated secondary bile acids such as deoxycholic acid. The current work represents an important advance in understanding the nutritional requirements and response to bile acids of the medically important human gut bacterium, C. scindens ATCC 35704. A defined medium has been developed which will further the understanding of bile acid metabolism in the context of growth substrates, cofactors, and other metabolites in the vertebrate gut. Analysis of the complete genome supports the nutritional requirements reported here. Genome-wide transcriptomic analysis of gene expression in the presence of cholic acid and deoxycholic acid provides a unique insight into the complex response of C. scindens ATCC 35704 to primary and secondary bile acids. Also revealed are genes with the potential to function in bile acid transport and metabolism.

Genomic comparison of Trypanosoma conorhini and Trypanosoma rangeli to Trypanosoma cruzi strains of high and low virulence.

BACKGROUND: Trypanosoma conorhini and Trypanosoma rangeli, like Trypanosoma cruzi, are kinetoplastid protist parasites of mammals displaying divergent hosts, geographic ranges and lifestyles. Largely nonpathogenic T. rangeli and T. conorhini represent clades that are phylogenetically closely related to the T. cruzi and T. cruzi-like taxa and provide insights into the evolution of pathogenicity in those parasites. T. rangeli, like T. cruzi is endemic in many Latin American countries, whereas T. conorhini is tropicopolitan. T. rangeli and T. conorhini are exclusively extracellular, while T. cruzi has an intracellular stage in the mammalian host. RESULTS: Here we provide the first comprehensive sequence analysis of T. rangeli AM80 and T. conorhini 025E, and provide a comparison of their genomes to those of T. cruzi G and T. cruzi CL, respectively members of T. cruzi lineages TcI and TcVI. We report de novo assembled genome sequences of the low-virulent T. cruzi G, T. rangeli AM80, and T. conorhini 025E ranging from ~ 21-25 Mbp, with ~ 10,000 to 13,000 genes, and for the highly virulent and hybrid T. cruzi CL we present a ~ 65 Mbp in-house assembled haplotyped genome with ~ 12,500 genes per haplotype. Single copy orthologs of the two T. cruzi strains exhibited ~ 97% amino acid identity, and ~ 78% identity to proteins of T. rangeli or T. conorhini. Proteins of the latter two organisms exhibited ~ 84% identity. T. cruzi CL exhibited the highest heterozygosity. T. rangeli and T. conorhini displayed greater metabolic capabilities for utilization of complex carbohydrates, and contained fewer retrotransposons and multigene family copies, i.e. trans-sialidases, mucins, DGF-1, and MASP, compared to T. cruzi. CONCLUSIONS: Our analyses of the T. rangeli and T. conorhini genomes closely reflected their phylogenetic proximity to the T. cruzi clade, and were largely consistent with their divergent life cycles. Our results provide a greater context for understanding the life cycles, host range expansion, immunity evasion, and pathogenesis of these trypanosomatids.

The reduced genome of Candidatus Kinetoplastibacterium sorsogonicusi, the endosymbiont of Kentomonas sorsogonicus (Trypanosomatidae): loss of the haem-synthesis pathway.

Trypanosomatids of the genera Angomonas and Strigomonas (subfamily Strigomonadinae) have long been known to contain intracellular beta-proteobacteria, which provide them with many important nutrients such as haem, essential amino acids and vitamins. Recently, Kentomonas sorsogonicus, a divergent member of Strigomonadinae, has been described. Herein, we characterize the genome of its endosymbiont, Candidatus Kinetoplastibacterium sorsogonicusi. This genome is completely syntenic with those of other known Ca. Kinetoplastibacterium spp., but more reduced in size (~742 kb, compared with 810-833 kb, respectively). Gene losses are not concentrated in any hot-spots but are instead distributed throughout the genome. The most conspicuous loss is that of the haem-synthesis pathway. For long, removing haemin from the culture medium has been a standard procedure in cultivating trypanosomatids isolated from insects; continued growth was considered as an evidence of endosymbiont presence. However, we demonstrate that, despite bearing the endosymbiont, K. sorsogonicus cannot grow in culture without haem. Thus, the traditional test cannot be taken as a reliable criterion for the absence or presence of endosymbionts in trypanosomatid flagellates. It remains unclear why the ability to synthesize such an essential compound was lost in Ca. K. sorsogonicusi, whereas all other known bacterial endosymbionts of trypanosomatids retain them.