The hidden genomic diversity of ciliated protists revealed by single-cell genome sequencing

Background Ciliated protists are a widely distributed, morphologically diverse, and genetically heterogeneous group of unicellular organisms, usually known for containing two types of nuclei: a transcribed polyploid macronucleus involved in gene expression and a silent diploid micronucleus responsible for transmission of genetic material during sexual reproduction and generation of the macronucleus. Although studies in a few species of culturable ciliated protists have revealed the highly dynamic nature of replicative and recombination events relating the micronucleus to the macronucleus, the broader understanding of the genomic diversity of ciliated protists, as well as their phylogenetic relationships and metabolic potential, has been hampered by the inability to culture numerous other species under laboratory conditions, as well as the presence of symbiotic bacteria and microalgae which provide a challenge for current sequencing technologies. Here, we optimized single-cell sequencing methods and associated data analyses, to effectively remove contamination by commensal bacteria, and generated high-quality genomes for a number of Euplotia species. Results We obtained eight high-quality Euplotia genomes by using single-cell genome sequencing techniques. The genomes have high genomic completeness, with sizes between 68 and 125 M and gene numbers between 14K and 25K. Through comparative genomic analysis, we found that there are a large number of gene expansion events in Euplotia genomes, and these expansions are closely related to the phenotypic evolution and specific environmental adaptations of individual species. We further found four distinct subgroups in the genus Euplotes, which exhibited considerable genetic distance and relative lack of conserved genomic syntenies. Comparative genomic analyses of Uronychia and its relatives revealed significant gene expansion associated with the ciliary movement machinery, which may be related to the unique and strong swimming ability. Conclusions We employed single-cell genomics to obtain eight ciliate genomes, characterized the underestimated genomic diversity of Euplotia, and determined the divergence time of representative species in this subclass for the first time. We also further investigated the extensive duplication events associated with speciation and environmental adaptation. This study provides a unique and valuable resource for understanding the evolutionary history and genetic diversity of ciliates.

Whole chromosome loss and genomic instability in mouse embryos after CRISPR-Cas9 genome editing

Karyotype alterations have emerged as on-target complications from CRISPR-Cas9 genome editing. However, the events that lead to these karyotypic changes in embryos after Cas9-treatment remain unknown. Here, using imaging and single-cell genome sequencing of 8-cell stage embryos, we track both spontaneous and Cas9-induced karyotype aberrations through the first three divisions of embryonic development. We observe the generation of abnormal structures of the nucleus that arise as a consequence of errors in mitosis, including micronuclei and chromosome bridges, and determine their contribution to common karyotype aberrations including whole chromosome loss that has been recently reported after editing in embryos. Together, these data demonstrate that Cas9-mediated germline genome editing can lead to unwanted on-target side effects, including major chromosome structural alterations that can be propagated over several divisions of embryonic development.

Three Draft Single-Cell Genome Sequences of Novel SAR324 Strains Isolated from the Abyssopelagic Southern Ocean

SAR324 is a ubiquitous and phylogenetically distinct clade of Deltaproteobacteria in marine environments. Here, we present three single-cell amplified genome sequences from the SAR324 lineage, obtained from the abyssopelagic zone of the Indian sector of the Southern Ocean.

Reconstructing and counting genomic fragments through tagmentation-based haploid phasing

Single-cell sequencing provides a new level of granularity in studying the heterogeneous nature of cancer cells. For some cancers, this heterogeneity is the result of copy number changes of genes within the cellular genomes. The ability to accurately determine such copy number changes is critical in tracing and understanding tumorigenesis. Current single-cell genome sequencing methodologies infer copy numbers based on statistical approaches followed by rounding decimal numbers to integer values. Such methodologies are sample dependent, have varying calling sensitivities which heavily depend on the sample’s ploidy and are sensitive to noise in sequencing data. In this paper we have demonstrated the concept of integer-counting by using a novel bioinformatic algorithm built on our library construction chemistry in order to detect the discrete nature of the genome.

A unified model of human hemoglobin switching through single-cell genome editing

Key mechanisms of fetal hemoglobin (HbF) regulation and switching have been elucidated through studies of human genetic variation, including mutations in the HBG1/2 promoters, deletions in the β-globin locus, and variation impacting BCL11A. While this has led to substantial insights, there has not been a unified understanding of how these distinct genetically-nominated elements, as well as other key transcription factors such as ZBTB7A, collectively interact to regulate HbF. A key limitation has been the inability to model specific genetic changes in primary isogenic human hematopoietic cells to uncover how each of these act individually and in aggregate. Here, we describe a single-cell genome editing functional assay that enables specific mutations to be recapitulated individually and in combination, providing insights into how multiple mutation-harboring functional elements collectively contribute to HbF expression. In conjunction with quantitative modeling and chromatin capture analyses, we illustrate how these genetic findings enable a comprehensive understanding of how distinct regulatory mechanisms can synergistically modulate HbF expression.

High throughput single-cell genome sequencing gives insights into the generation and evolution of mosaic aneuploidy in Leishmania donovani

Leishmania, a unicellular eukaryotic parasite, is a unique model for aneuploidy and cellular heterogeneity, along with their potential role in adaptation to environmental stresses. Somy variation within clonal populations was previously explored in a small subset of chromosomes using fluorescence hybridization methods. This phenomenon, termed mosaic aneuploidy (MA) might have important evolutionary and functional implications, but remains under-explored due to technological limitations. Here, we applied and validated a high throughput single-cell genome sequencing method to study for the first time the extent and dynamics of whole karyotype heterogeneity in two Leishmania clonal populations representing different stages of MA evolution in vitro. We found that drastic changes in karyotypes quickly emerge in a population stemming from an almost euploid founder cell. This possibly involves polyploidization/hybridization at an early stage of population expansion, followed by assorted ploidy reduction. During further stages of expansion, MA increases by moderate and gradual karyotypic alterations. MA usually affected a defined subset of chromosomes, of which some display enrichment in snoRNA genes which could represent an adaptative benefit to the amplification of these chromosomes. Our data provide the first complete characterization of MA in Leishmania and pave the way for further functional studies.

Single-cell landscapes of primary glioblastomas and matched organoids and cell lines reveal variable retention of inter- and intra-tumor heterogeneity

SummaryGlioblastomas (GBMs) are aggressive primary malignant brain tumors characterized by extensive levels of inter- and intra-tumor genetic and phenotypic heterogeneity. Patient-derived organoids (PDOs) have recently emerged as useful models to study such heterogeneity. Here, we present bulk exome as well as single-cell genome and transcriptome profiles of primary IDH wild type GBMs from ten patients, including two recurrent tumors, as well as PDOs and brain tumor-initiating cell (BTIC) lines derived from these patients. We find that PDOs are genetically similar to and variably retain gene expression characteristics of their parent tumors. At the phenotypic level, PDOs appear to exhibit similar levels of transcriptional heterogeneity as their parent tumors, whereas BTIC lines tend to be enriched for cells in a more uniform transcriptional state. The datasets introduced here will provide a valuable resource to help guide experiments using GBM-derived organoids, especially in the context of studying cellular heterogeneity.

Recovery of high-quality assembled genomes via single-cell genome-guided binning of metagenome assembly

High-quality (HQ) reference genomes are essential for understanding the phylogeny and function of uncultured microbes in complex microbial ecosystems. However, existing metagenomic binners often fail to reconstruct a reasonable number of reliable HQ genomes owing to a lack of ideal binning guides. Here, we present a single-cell genome-guided binning of metagenomic assemblies (SIGMA) to reconstruct the HQ genomes of multiple strains from microbial communities at once. SIGMA generates self-reference sequences from the same sample by single-cell sequencing and uses them as guides to reconstruct metagenomic bins. The single-cell genome guide enabled precise binning and sequence integration and produced the largest number of HQ genomes from mock community and human microbiota samples in comparison with conventional binners. SIGMA can recover rRNA and tRNA genes and link plasmids to the host. This ability will contribute to understanding intraspecies diversity and distribution of mobile genetic elements in uncultured microbes.