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.

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.

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.

Identification of a gene encoding a flavoprotein involved in bile acid metabolism by the human gut bacterium Clostridium scindens ATCC 35704.

BACKGROUND: The multi-step bile acid 7α-dehydroxylating pathway by which a few species of Clostridium convert host primary bile acids to toxic secondary bile acids is of great importance to gut microbiome structure and host physiology and disease. While genes in the oxidative arm of the 7α-dehydroxylating pathway have been identified, genes in the reductive arm of the pathway are still obscure. METHODS: We identified a candidate flavoprotein-encoding gene predicted to metabolize steroids. This gene was cloned and overexpressed in E. coli and affinity purified. Reaction substrate and product were separated by thin layer chromatography and identified by liquid chromatograph mass spectrometry-ion trap-time of flight (LCMS-IT-TOF). Phylogenetic analysis of the amino acid sequence was performed. RESULTS: We report the identification of a gene encoding a flavoprotein (EDS08212.1) involved in secondary bile acid metabolism by Clostridium scindens ATCC 35704 and related species. Purified rEDS08212.1 catalyzed formation of a product from 3-dehydro-deoxycholic acid that UPLC-IT-TOF-MS analysis suggests loses 4amu. Our phylogeny identified this gene in other bile acid 7α-dehydroxylating bacteria. CONCLUSIONS: These data suggest formation of a product, 3-dehydro-4,6-deoxycholic acid, a recognized intermediate in the reductive arm of bile acid 7α-dehydroxylation pathway and the first report of a gene in the reductive arm of the bile acid 7α-dehydroxylating pathway.

Metabolism of Oxo-Bile Acids and Characterization of Recombinant 12α-Hydroxysteroid Dehydrogenases from Bile Acid 7α-Dehydroxylating Human Gut Bacteria.

Bile acids are important cholesterol-derived nutrient signaling hormones, synthesized in the liver, that act as detergents to solubilize dietary lipids. Bile acid 7α-dehydroxylating gut bacteria generate the toxic bile acids deoxycholic acid and lithocholic acid from host bile acids. The ability of these bacteria to remove the 7-hydroxyl group is partially dependent on 7α-hydroxysteroid dehydrogenase (HSDH) activity, which reduces 7-oxo-bile acids generated by other gut bacteria. 3α-HSDH has an important enzymatic activity in the bile acid 7α-dehydroxylation pathway. 12α-HSDH activity has been reported for the low-activity bile acid 7α-dehydroxylating bacterium Clostridium leptum; however, this activity has not been reported for high-activity bile acid 7α-dehydroxylating bacteria, such as Clostridium scindens, Clostridium hylemonae, and Clostridium hiranonis Here, we demonstrate that these strains express bile acid 12α-HSDH. The recombinant enzymes were characterized from each species and shown to preferentially reduce 12-oxolithocholic acid to deoxycholic acid, with low activity against 12-oxochenodeoxycholic acid and reduced activity when bile acids were conjugated to taurine or glycine. Phylogenetic analysis suggests that 12α-HSDH is widespread among Firmicutes, Actinobacteria in the Coriobacteriaceae family, and human gut ArchaeaIMPORTANCE 12α-HSDH activity has been established in the medically important bile acid 7α-dehydroxylating bacteria C. scindens, C. hiranonis, and C. hylemonae Experiments with recombinant 12α-HSDHs from these strains are consistent with culture-based experiments that show a robust preference for 12-oxolithocholic acid over 12-oxochenodeoxycholic acid. Phylogenetic analysis identified novel members of the gut microbiome encoding 12α-HSDH. Future reengineering of 12α-HSDH enzymes to preferentially oxidize cholic acid may provide a means to industrially produce the therapeutic bile acid ursodeoxycholic acid. In addition, a cholic acid-specific 12α-HSDH expressed in the gut may be useful for the reduction in deoxycholic acid concentration, a bile acid implicated in cancers of the gastrointestinal (GI) tract.

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.

The Evolutionary Loss of RNAi Key Determinants in Kinetoplastids as a Multiple Sporadic Phenomenon.

We screened the genomes of a broad panel of kinetoplastid protists for genes encoding proteins associated with the RNA interference (RNAi) system using probes from the Argonaute (AGO1), Dicer1 (DCL1), and Dicer2 (DCL2) genes of Leishmania brasiliensis and Crithidia fasciculata. We identified homologs for all the three of these genes in the genomes of a subset of these organisms. However, several of these organisms lacked evidence for any of these genes, while others lacked only DCL2. The open reading frames encoding these putative proteins were structurally analyzed in silico. The alignments indicated that the genes are homologous with a high degree of confidence, and three-dimensional structural models strongly supported a functional relationship to previously characterized AGO1, DCL1, and DCL2 proteins. Phylogenetic analysis of these putative proteins showed that these genes, when present, evolved in parallel with other nuclear genes, arguing that the RNAi system genes share a common progenitor, likely across all Kinetoplastea. In addition, the genome segments bearing these genes are highly conserved and syntenic, even among those taxa in which they are absent. However, taxa in which these genes are apparently absent represent several widely divergent branches of kinetoplastids, arguing that these genes were independently lost at least six times in the evolutionary history of these organisms. The mechanisms responsible for the apparent coordinate loss of these RNAi system genes independently in several lineages of kinetoplastids, while being maintained in other related lineages, are currently unknown.

Merozoite surface protein-1 genetic diversity in Plasmodium malariae and Plasmodium brasilianum from Brazil.

BACKGROUND: The merozoite surface protein 1 (MSP1) gene encodes the major surface antigen of invasive forms of the Plasmodium erythrocytic stages and is considered a candidate vaccine antigen against malaria. Due to its polymorphisms, MSP1 is also useful for strain discrimination and consists of a good genetic marker. Sequence diversity in MSP1 has been analyzed in field isolates of three human parasites: P. falciparum, P. vivax, and P. ovale. However, the extent of variation in another human parasite, P. malariae, remains unknown. This parasite shows widespread, uneven distribution in tropical and subtropical regions throughout South America, Asia, and Africa. Interestingly, it is genetically indistinguishable from P. brasilianum, a parasite known to infect New World monkeys in Central and South America. METHODS: Specific fragments (1 to 5) covering 60 % of the MSP1 gene (mainly the putatively polymorphic regions), were amplified by PCR in isolates of P. malariae and P. brasilianum from different geographic origin and hosts. Sequencing of the PCR-amplified products or cloned PCR fragments was performed and the sequences were used to construct a phylogenetic tree by the maximum likelihood method. Data were computed to give insights into the evolutionary and phylogenetic relationships of these parasites. RESULTS: Except for fragment 4, sequences from all other fragments consisted of unpublished sequences. The most polymorphic gene region was fragment 2, and in samples where this region lacks polymorphism, all other regions are also identical. The low variability of the P. malariae msp1 sequences of these isolates and the identification of the same haplotype in those collected many years apart at different locations is compatible with a low transmission rate. We also found greater diversity among P. brasilianum isolates compared with P. malariae ones. Lastly, the sequences were segregated according to their geographic origins and hosts, showing a strong genetic and geographic structure. CONCLUSIONS: Our data show that there is a low level of sequence diversity and a possible absence of allelic dimorphism of MSP1 in these parasites as opposed to other Plasmodium species. P. brasilianum strains apparently show greater divergence in comparison to P. malariae, thus P. malariae could derive from P. brasilianum, as it has been proposed.

Clostridium scindens: a human gut microbe with a high potential to convert glucocorticoids into androgens.

Clostridium scindens American Type Culture Collection 35704 is capable of converting primary bile acids to toxic secondary bile acids, as well as converting glucocorticoids to androgens by side-chain cleavage. The molecular structure of the side-chain cleavage product of cortisol produced by C. scindens was determined to be 11ß-hydroxyandrost-4-ene-3,17-dione (11ß-OHA) by high-resolution mass spectrometry, (1)H and (13)C NMR spectroscopy, and X-ray crystallography. Using RNA-Seq technology, we identified a cortisol-inducible (≈ 1,000-fold) operon (desABCD) encoding at least one enzyme involved in anaerobic side-chain cleavage. The desC gene was cloned, overexpressed, purified, and found to encode a 20α-hydroxysteroid dehydrogenase (HSDH). This operon also encodes a putative "transketolase" (desAB) hypothesized to have steroid-17,20-desmolase/oxidase activity, and a possible corticosteroid transporter (desD). RNA-Seq data suggests that the two-carbon side chain of glucocorticords may feed into the pentose-phosphate pathway and are used as a carbon source. The 20α-HSDH is hypothesized to function as a metabolic "rheostat" controlling rates of side-chain cleavage. Phylogenetic analysis suggests this operon is rare in nature and the desC gene evolved from a gene encoding threonine dehydrogenase. The physiological effect of 11ß-OHAD on the host or other gut microbes is currently unknown.

Endosymbiosis in trypanosomatids: the genomic cooperation between bacterium and host in the synthesis of essential amino acids is heavily influenced by multiple horizontal gene transfers.

BACKGROUND: Trypanosomatids of the genera Angomonas and Strigomonas live in a mutualistic association characterized by extensive metabolic cooperation with obligate endosymbiotic Betaproteobacteria. However, the role played by the symbiont has been more guessed by indirect means than evidenced. Symbiont-harboring trypanosomatids, in contrast to their counterparts lacking symbionts, exhibit lower nutritional requirements and are autotrophic for essential amino acids. To evidence the symbiont's contributions to this autotrophy, entire genomes of symbionts and trypanosomatids with and without symbionts were sequenced here. RESULTS: Analyses of the essential amino acid pathways revealed that most biosynthetic routes are in the symbiont genome. By contrast, the host trypanosomatid genome contains fewer genes, about half of which originated from different bacterial groups, perhaps only one of which (ornithine cyclodeaminase, EC:4.3.1.12) derived from the symbiont. Nutritional, enzymatic, and genomic data were jointly analyzed to construct an integrated view of essential amino acid metabolism in symbiont-harboring trypanosomatids. This comprehensive analysis showed perfect concordance among all these data, and revealed that the symbiont contains genes for enzymes that complete essential biosynthetic routes for the host amino acid production, thus explaining the low requirement for these elements in symbiont-harboring trypanosomatids. Phylogenetic analyses show that the cooperation between symbionts and their hosts is complemented by multiple horizontal gene transfers, from bacterial lineages to trypanosomatids, that occurred several times in the course of their evolution. Transfers occur preferentially in parts of the pathways that are missing from other eukaryotes. CONCLUSION: We have herein uncovered the genetic and evolutionary bases of essential amino acid biosynthesis in several trypanosomatids with and without endosymbionts, explaining and complementing decades of experimental results. We uncovered the remarkable plasticity in essential amino acid biosynthesis pathway evolution in these protozoans, demonstrating heavy influence of horizontal gene transfer events, from Bacteria to trypanosomatid nuclei, in the evolution of these pathways.

Cloning and characterization of bifunctional enzyme farnesyl diphosphate/geranylgeranyl diphosphate synthase from Plasmodium falciparum.

BACKGROUND: Isoprenoids are the most diverse and abundant group of natural products. In Plasmodium falciparum, isoprenoid synthesis proceeds through the methyl erythritol diphosphate pathway and the products are further metabolized by farnesyl diphosphate synthase (FPPS), turning this enzyme into a key branch point of the isoprenoid synthesis. Changes in FPPS activity could alter the flux of isoprenoid compounds downstream of FPPS and, hence, play a central role in the regulation of a number of essential functions in Plasmodium parasites. METHODS: The isolation and cloning of gene PF3D7_18400 was done by amplification from cDNA from mixed stage parasites of P. falciparum. After sequencing, the fragment was subcloned in pGEX2T for recombinant protein expression. To verify if the PF3D7_1128400 gene encodes a functional rPfFPPS protein, its catalytic activity was assessed using the substrate [4-14C] isopentenyl diphosphate and three different allylic substrates: dimethylallyl diphosphate, geranyl diphosphate or farnesyl diphosphate. The reaction products were identified by thin layer chromatography and reverse phase high-performance liquid chromatography. To confirm the product spectrum formed of rPfFPPS, isoprenic compounds were also identified by mass spectrometry. Apparent kinetic constants KM and Vmax for each substrate were determined by Michaelis-Menten; also, inhibition assays were performed using risedronate. RESULTS: The expressed protein of P. falciparum FPPS (rPfFPPS) catalyzes the synthesis of farnesyl diphosphate, as well as geranylgeranyl diphosphate, being therefore a bifunctional FPPS/geranylgeranyl diphosphate synthase (GGPPS) enzyme. The apparent KM values for the substrates dimethylallyl diphosphate, geranyl diphosphate and farnesyl diphosphate were, respectively, 68 ± 5 µM, 7.8 ± 1.3 µM and 2.06 ± 0.4 µM. The protein is expressed constitutively in all intra-erythrocytic stages of P. falciparum, demonstrated by using transgenic parasites with a haemagglutinin-tagged version of FPPS. Also, the present data demonstrate that the recombinant protein is inhibited by risedronate. CONCLUSIONS: The rPfFPPS is a bifunctional FPPS/GGPPS enzyme and the structure of products FOH and GGOH were confirmed mass spectrometry. Plasmodial FPPS represents a potential target for the rational design of chemotherapeutic agents to treat malaria.

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.

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.

Genomic sequence analysis and characterization of Sneathia amnii sp. nov.

BACKGROUND: Bacteria of the genus Sneathia are emerging as potential pathogens of the female reproductive tract. Species of Sneathia, which were formerly grouped with Leptotrichia, can be part of the normal microbiota of the genitourinary tracts of men and women, but they are also associated with a variety of clinical conditions including bacterial vaginosis, preeclampsia, preterm labor, spontaneous abortion, post-partum bacteremia and other invasive infections. Sneathia species also exhibit a significant correlation with sexually transmitted diseases and cervical cancer. Because Sneathia species are fastidious and rarely cultured successfully in vitro; and the genomes of members of the genus had until now not been characterized, very little is known about the physiology or the virulence of these organisms. RESULTS: Here, we describe a novel species, Sneathia amnii sp. nov, which closely resembles bacteria previously designated "Leptotrichia amnionii". As part of the Vaginal Human Microbiome Project at VCU, a vaginal isolate of S. amnii sp. nov. was identified, successfully cultured and bacteriologically cloned. The biochemical characteristics and virulence properties of the organism were examined in vitro, and the genome of the organism was sequenced, annotated and analyzed. The analysis revealed a reduced circular genome of ~1.34 Mbp, containing ~1,282 protein-coding genes. Metabolic reconstruction of the bacterium reflected its biochemical phenotype, and several genes potentially associated with pathogenicity were identified. CONCLUSIONS: Bacteria with complex growth requirements frequently remain poorly characterized and, as a consequence, their roles in health and disease are unclear. Elucidation of the physiology and identification of genes putatively involved in the metabolism and virulence of S. amnii may lead to a better understanding of the role of this potential pathogen in bacterial vaginosis, preterm birth, and other issues associated with vaginal and reproductive health.

On the structural plasticity of the human genome: chromosomal inversions revisited.

With the aid of novel and powerful molecular biology techniques, recent years have witnessed a dramatic increase in the number of studies reporting the involvement of complex structural variants in several genomic disorders. In fact, with the discovery of Copy Number Variants (CNVs) and other forms of unbalanced structural variation, much attention has been directed to the detection and characterization of such rearrangements, as well as the identification of the mechanisms involved in their formation. However, it has long been appreciated that chromosomes can undergo other forms of structural changes - balanced rearrangements - that do not involve quantitative variation of genetic material. Indeed, a particular subtype of balanced rearrangement - inversions - was recently found to be far more common than had been predicted from traditional cytogenetics. Chromosomal inversions alter the orientation of a specific genomic sequence and, unless involving breaks in coding or regulatory regions (and, disregarding complex trans effects, in their close vicinity), appear to be phenotypically silent. Such a surprising finding, which is difficult to reconcile with the classical interpretation of inversions as a mechanism causing subfertility (and ultimately reproductive isolation), motivated a new series of theoretical and empirical studies dedicated to understand their role in human genome evolution and to explore their possible association to complex genetic disorders. With this review, we attempt to describe the latest methodological improvements to inversions detection at a genome wide level, while exploring some of the possible implications of inversion rearrangements on the evolution of the human genome.

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 .

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.

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.

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.