Cryptosporidium as a testbed for single cell genome characterization of unicellular eukaryotes.

BACKGROUND: Infectious disease involving multiple genetically distinct populations of pathogens is frequently concurrent, but difficult to detect or describe with current routine methodology. Cryptosporidium sp. is a widespread gastrointestinal protozoan of global significance in both animals and humans. It cannot be easily maintained in culture and infections of multiple strains have been reported. To explore the potential use of single cell genomics methodology for revealing genome-level variation in clinical samples from Cryptosporidium-infected hosts, we sorted individual oocysts for subsequent genome amplification and full-genome sequencing. RESULTS: Cells were identified with fluorescent antibodies with an 80 % success rate for the entire single cell genomics workflow, demonstrating that the methodology can be applied directly to purified fecal samples. Ten amplified genomes from sorted single cells were selected for genome sequencing and compared both to the original population and a reference genome in order to evaluate the accuracy and performance of the method. Single cell genome coverage was on average 81 % even with a moderate sequencing effort and by combining the 10 single cell genomes, the full genome was accounted for. By a comparison to the original sample, biological variation could be distinguished and separated from noise introduced in the amplification. CONCLUSIONS: As a proof of principle, we have demonstrated the power of applying single cell genomics to dissect infectious disease caused by closely related parasite species or subtypes. The workflow can easily be expanded and adapted to target other protozoans, and potential applications include mapping genome-encoded traits, virulence, pathogenicity, host specificity and resistance at the level of cells as truly meaningful biological units.

A pan-genomic approach reveals novel Sulfurimonas clade in the ferruginous meromictic Lake Pavin.

The permanently anoxic waters in meromictic lakes create suitable niches for the growth of bacteria using sulphur metabolisms like sulphur oxidation. In Lake Pavin, the anoxic water mass hosts an active cryptic sulphur cycle that interacts narrowly with iron cycling, however the metabolisms of the microorganisms involved are poorly known. Here we combined metagenomics, single-cell genomics, and pan-genomics to further expand our understanding of the bacteria and the corresponding metabolisms involved in sulphur oxidation in this ferruginous sulphide- and sulphate-poor meromictic lake. We highlighted two new species within the genus Sulfurimonas that belong to a novel clade of chemotrophic sulphur oxidisers exclusive to freshwaters. We moreover conclude that this genus holds a key-role not only in limiting sulphide accumulation in the upper part of the anoxic layer but also constraining carbon, phosphate and iron cycling.

scDual-Seq of Toxoplasma gondii-infected mouse BMDCs reveals heterogeneity and differential infection dynamics.

Dendritic cells and macrophages are integral parts of the innate immune system and gatekeepers against infection. The protozoan pathogen, Toxoplasma gondii, is known to hijack host immune cells and modulate their immune response, making it a compelling model to study host-pathogen interactions. Here we utilize single cell Dual RNA-seq to parse out heterogeneous transcription of mouse bone marrow-derived dendritic cells (BMDCs) infected with two distinct genotypes of T. gondii parasites, over multiple time points post infection. We show that the BMDCs elicit differential responses towards T. gondii infection and that the two parasite lineages distinctly manipulate subpopulations of infected BMDCs. Co-expression networks define host and parasite genes, with implications for modulation of host immunity. Integrative analysis validates previously established immune pathways and additionally, suggests novel candidate genes involved in host-pathogen interactions. Altogether, this study provides a comprehensive resource for characterizing host-pathogen interplay at high-resolution.

Genomes from uncultivated prokaryotes: a comparison of metagenome-assembled and single-amplified genomes.

BACKGROUND: Prokaryotes dominate the biosphere and regulate biogeochemical processes essential to all life. Yet, our knowledge about their biology is for the most part limited to the minority that has been successfully cultured. Molecular techniques now allow for obtaining genome sequences of uncultivated prokaryotic taxa, facilitating in-depth analyses that may ultimately improve our understanding of these key organisms. RESULTS: We compared results from two culture-independent strategies for recovering bacterial genomes: single-amplified genomes and metagenome-assembled genomes. Single-amplified genomes were obtained from samples collected at an offshore station in the Baltic Sea Proper and compared to previously obtained metagenome-assembled genomes from a time series at the same station. Among 16 single-amplified genomes analyzed, seven were found to match metagenome-assembled genomes, affiliated with a diverse set of taxa. Notably, genome pairs between the two approaches were nearly identical (average 99.51% sequence identity; range 98.77-99.84%) across overlapping regions (30-80% of each genome). Within matching pairs, the single-amplified genomes were consistently smaller and less complete, whereas the genetic functional profiles were maintained. For the metagenome-assembled genomes, only on average 3.6% of the bases were estimated to be missing from the genomes due to wrongly binned contigs. CONCLUSIONS: The strong agreement between the single-amplified and metagenome-assembled genomes emphasizes that both methods generate accurate genome information from uncultivated bacteria. Importantly, this implies that the research questions and the available resources are allowed to determine the selection of genomics approach for microbiome studies.

Universal tag arrays in forensic SNP analysis.

Microarray-based single nucleotide polymorphism (SNP) genotyping enables simultaneous and rapid detection of a large number of markers and is thus an attractive method for forensic individual acid identification. This assay relies on a one-color detection system and minisequencing in solution before hybridization to universal tag arrays. The minisequencing reaction is based on incorporation of a fluorescent dideoxynucleotide to a primer containing a tag-sequence flanking the position to be interrogated. This one-color system detects C and T polymorphisms in separate reactions on multiple polymerase chain reaction targets with the fluorophore TAMRA coupled to the respective dideoxynucleotide. After incorporation, tagged primer sequences are hybridized through their complementary sequence on the array, and positive signals are detected by a confocal laser-scanner.

A DNA microarray system for forensic SNP analysis.

Forensic DNA analysis is routinely performed using polymorphic short tandem repeat (STR) markers. However, for degraded or minute DNA samples, analysis of autosomal single nucleotide polymorphisms (SNPs) in short fragments might be more successful. Furthermore, sequencing of mitochondrial DNA (mtDNA) is often performed on highly degraded or scarce samples due to the high copy number of mtDNA in each cell. Due to the increasing number of complete mtDNA genome sequences available, the limited discrimination power of an mtDNA analysis, may be increased by analysis of coding region polymorphisms in addition to the non-coding variation. Since sequence analysis of the coding region would require more material than generally present in forensic samples, an alternative SNP analysis approach is required. We have developed a one-colour microarray-based SNP detection system for limited forensic materials. The method is based on minisequencing in solution prior to hybridisation to universal tag-arrays. In a first outline of a forensic chip, a combination of 12 nuclear and 21 mitochondrial SNP markers are analysed simultaneously. The mitochondrial markers on the chip are polymorphisms within the hypervariable region as well as in the coding region. Even though the number of markers in the current system is limited, it can easily be extended to yield a greater power of discrimination. When fully developed, microarray analysis provides a promising system for efficient sensitive SNP analysis of forensic samples in the future.

Forensic casework analysis using the HVI/HVII mtDNA linear array assay.

The mitochondrial hypervariable regions I and II have proven to be a useful target for analysis of forensic materials, in which the amount of DNA is limited or highly degraded. Conventional mitochondrial DNA (mtDNA) sequencing can be time-consuming and expensive, limitations that can be minimized using a faster and less expensive typing assay. We have evaluated the exclusion capacity of the linear array mtDNA HVI/HVII region-sequence typing assay (Roche Applied Science) in 16 forensic cases comprising 90 samples. Using the HVI/HVII mtDNA linear array, 56% of the samples were excluded and thus less than half of the samples require further sequencing due to a match or inconclusive results. Of all the samples that were excluded by sequence analysis, 79% could be excluded using the HVI/HVII linear array alone. Using the HVI/HVII mtDNA linear array assay, we demonstrate the potential to decrease sequencing efforts substantially and thereby reduce the cost and the turn-around time in casework analysis.

Forensic analysis of autosomal STR markers using Pyrosequencing.

Short tandem repeats (STRs) are highly variable, and therefore routinely used in forensic investigations for a DNA-based individual identification. The routine assay is commonly performed by size separation using capillary electrophoresis, but alternative technologies can also be used. In this study, a Pyrosequencing assay was developed for analysis of STR markers useful in forensic DNA analysis. The assay was evaluated for 10 different STR loci (CSF1PO, TH01, TPOX, D3S1358, D5S818, D7S820, D8S1179, D13S317, D16S539 and Penta E) and a total of 114 Swedish individuals were genotyped. This genotyping strategy reveal the actual sequence and variant alleles were seen at several loci, providing additional information compared to fragment size analysis. At the D13S317 locus a T/A SNP located in the last repeat unit was observed in 92% of the genotypes. Moreover, an upstream flanking SNP at locus D7S820, a SNP within the repeats at D3S1358 and D8S1179 and a deletion in the flanking region at locus D5S818 were observed. The Pyrosequencing method was first developed for SNP typing and sequencing of shorter DNA fragments but the method also provides an alternative method for STR analysis of less complex repeats. This assay is suitable for investigation of new markers, a rapid compilation of population data and for confirmation of variant and new alleles.