Revealing gene expression heterogeneity in a clonal population of Tetrahymena thermophila through single-cell RNA sequencing.

We performed single-cell RNA sequencing (scRNA-seq) on a population of 5,000 Tetrahymena thermophila, using the 10x Genomics 3' gene expression analysis, to investigate gene expression variability within this clonal population. Initially, we estimated the 3'-untranslated regions (3' UTRs), which were absent in existing annotation files but are crucial for the 10x Genomics 3' gene expression analysis, using the peaks2utr method. This allowed us to create a modified annotation file, which was then utilized in our scRNA-seq analysis. Our analysis revealed significant gene expression variability within the population, even after removing the effect of cell phase-related features. This variability predominantly appeared in six distinct clusters. Through gene ontology and KEGG pathway enrichment analyses, we identified that these were primarily associated with ribosomal proteins, proteins specific to mitochondria, proteins involved in peroxisome-specific carbon metabolism, cytoskeletal proteins, motor proteins, and immobilized antigens.

Comparative genomics hints at dispensability of multiple essential genes in two Escherichia coli L-form strains.

Despite the critical role of bacterial cell walls in maintaining cell shapes, certain environmental stressors can induce the transition of many bacterial species into a wall-deficient state called L-form. Long-term induced Escherichia coli L-forms lose their rod shape and usually hold significant mutations that affect cell division and growth. Besides this, the genetic background of L-form bacteria is still poorly understood. In the present study, the genomes of two stable L-form strains of E. coli (NC-7 and LWF+) were sequenced and their gene mutation status was determined and compared with their parental strains. Comparative genomic analysis between two L-forms reveals both unique adaptions and common mutated genes, many of which belong to essential gene categories not involved in cell wall biosynthesis, indicating that L-form genetic adaptation impacts crucial metabolic pathways. Missense variants from L-forms and Lenski's long-term evolution experiment (LTEE) were analyzed in parallel using an optimized DeepSequence pipeline to investigate predicted mutation effects (α) on protein functions. We report that the two L-form strains analyzed display a frequency of 6-10% (0% for LTEE) in mutated essential genes where the missense variants have substantial impact on protein functions (α<0.5). This indicates the emergence of different survival strategies in L-forms through changes in essential genes during adaptions to cell wall deficiency. Collectively, our results shed light on the detailed genetic background of two E. coli L-forms and pave the way for further investigations of the gene functions in L-form bacterial models.

Improvement of solubility of phospholipase D from Streptomyces antibioticus in recombinant Escherichia coli and its application for the enzymatic synthesis of a non-natural plasmalogen.

Plasmalogens are a subclass of glycerophospholipids that have a vinyl-ether bond at the sn-1 position and are thought to have several physiological functions. The creation of non-natural plasmalogens with functional groups is desired for the establishment of the prevention of diseases caused by the depletion of plasmalogens. Phospholipase D (PLD) has both hydrolysis and transphosphatidylation activities. In particular, PLD from Streptomyces antibioticus has been investigated extensively due to its high transphosphatidylation activity. However, it has been difficult to stably express recombinant PLD in Escherichia coli and to express it as a soluble protein. In this study, we used the E. coli strain, SoluBL21™, and achieved stable PLD expression from the T7 promoter and increased soluble fraction in the cell. We also improved the purification method of PLD using His-tag at the C terminus. We obtained PLD with ∼730 mU mg-1 protein of specific activity, and the yield was ∼420 mU l-1 culture, corresponding to 76 mU per gram of wet cells. Finally, we synthesized a non-natural plasmalogen with 1,4-cyclohexanediol bound to the phosphate group at the sn-3 position by transphosphatidylation of the purified PLD. This method will contribute to the expansion of the chemical structure library of non-natural plasmalogens.

Construction of a mini-RNA replicon in Escherichia coli.

How the ribonucleic acid (RNA) world transited to the deoxyribonucleic acid (DNA) world has remained controversial in evolutionary biology. At a certain time point in the transition from the RNA world to the DNA world, 'RNA replicons', in which RNAs produce proteins to replicate their coding RNA, and 'DNA replicons', in which DNAs produce RNA to synthesize proteins that replicate their coding DNA, can be assumed to coexist. The coexistent state of RNA replicons and DNA replicons is desired for experimental approaches to determine how the DNA world overtook the RNA world. We constructed a mini-RNA replicon in Escherichia coli. This mini-RNA replicon encoded the ß subunit, one of the subunits of the Qß replicase derived from the positive-sense single-stranded Qß RNA phage and is replicated by the replicase in E. coli. To maintain the mini-RNA replicon persistently in E. coli cells, we employed a system of α complementation of LacZ that was dependent on the Qß replicase, allowing the cells carrying the RNA replicon to grow in the lactose minimal medium selectively. The coexistent state of the mini-RNA replicon and DNA replicon (E. coli genome) was successively synthesized. The coexistent state can be used as a starting system to experimentally demonstrate the transition from the RNA-protein world to the DNA world, which will contribute to progress in the research field of the origin of life.

Lymphatic Absorption of Microbial Plasmalogens in Rats.

Plasmalogens, functional glycerophospholipids with biological roles in the human body, are associated with various diseases. Although a variety of saturated and/or unsaturated fatty acids in plasmalogens are presumed to have different functions in the human body, there are limited reports validating such functions of plasmalogens. In this study, we focused on the bacterial plasmalogen derived from Selenomonas ruminantium subsp. lactilytica (NBRC No. 103574) with different main species of hydrocarbon chains at the sn-1 position and shorter fatty acids at the sn-2 position than animal plasmalogens. Optimum culture conditions of S. ruminantium for high-yield production of plasmalogens, such as pH and the concentration of caproic acid, were investigated under anaerobic conditions using a 2-L scale jar fermenter. The obtained plasmalogen mainly consisted of the ethanolamine plasmalogen (PlsEtn). The molar ratios of PlsEtn species obtained from S. ruminantium, at sn-1/sn-2 positions, were p16:1/14:0 (68.4%), p16:1/16:1 (29.2%), p16:1/16:0 (0.7%), p16:1/15:0 (0.3%), and p17:1/14:0 (0.3%). Subsequently, duodenal infusion of the emulsion carrying the lipid extracted from S. ruminantium was carried out in lymph duct-cannulated rats. In the lymphatic plasmalogen of rats, the level of PlsEtns with molar ratios p16:1/14:0 and p16:1/16:1, the main species of plasmalogens from S. ruminantium, increased gradually until 3-4 h after lipid injection and then gradually decreased. In addition, the level of PlsEtns with p16:1/20:4 and p16:1/22:6 rapidly increased, peaking at 1-1.5 h and 1.5-2 h after lipid injection, respectively. The increase in the number of PlsEtns with p16:1/20:4 and p16:1/22:6 suggested that 20:4 and 22:6, the main fatty acids at the sn-2 position in the rat lymphatic plasmalogen, were preferentially re-esterified at the sn-2 position, regardless of the types of hydrocarbon chains at the sn-1 position. Thus, we showed that bacterial PlsEtns with "unnatural" structures against rats could be absorbed into the lymph. Our findings provide insights into the association between the chemical structure of plasmalogens and their biological functions in humans.

Host selection-producing variations in the genome of hop stunt viroid.

A random mutant pool of hop stunt viroid (HSVd) was created by shuffling cDNA fragments prepared from three natural HSVd variants obtained from grapevine, citrus, and plum. It was used to infect five host plant species: hop, cucumber, grapevine, peach, and citrus. After infection, progenies having variations characteristic for grapevine and citrus HSVd variants have been preferentially enriched in the homologous plant species, suggesting that strong but different selection pressures affected the genomic RNA when HSVd-infected either grapevine or citrus. In the progeny propagated in cucumber, hop, and peach, variations characteristic to grapevine, citrus, and plum HSVd variants were detected simultaneously as a blend. Accordingly, we showed that at least some of the host-specific variations found in HSVd variants isolated from host plant species, e.g., grapevine and citrus, seemed to have arisen from positive host selection pressures. The HSVd-grapevine variant was found to be ideally adaptable not only to grapevine but to various host plants as well.

Microbial Diversity in the Phyllosphere and Rhizosphere of an Apple Orchard Managed under Prolonged "Natural Farming" Practices.

Microbial diversity in an apple orchard cultivated with natural farming practices for over 30 years was compared with conventionally farmed orchards to analyze differences in disease suppression. In this long-term naturally farmed orchard, major apple diseases were more severe than in conventional orchards but milder than in a short-term natural farming orchard. Among major fungal species in the phyllosphere, we found that Aureobasidium pullulans and Cryptococcus victoriae were significantly less abundant in long-term natural farming, while Cladosporium tenuissimum predominated. However, diversity of fungal species in the phyllosphere was not necessarily the main determinant in the disease suppression observed in natural farming; instead, the maintenance of a balanced, constant selection of fungal species under a suitable predominant species such as C. tenuissimum seemed to be the important factors. Analysis of bacteria in the phyllosphere revealed Pseudomonas graminis, a potential inducer of plant defenses, predominated in long-term natural farming in August. Rhizosphere metagenome analysis showed that Cordyceps and Arthrobotrys, fungal genera are known to include insect- or nematode-infecting species, were found only in long-term natural farming. Among soil bacteria, the genus Nitrospira was most abundant, and its level in long-term natural farming was more than double that in the conventionally farmed orchard.

The Single-Stranded RNA Bacteriophage Qß Adapts Rapidly to High Temperatures: An Evolution Experiment.

Single-stranded (ss)RNA viruses are thought to evolve rapidly due to an inherently high mutation rate. However, it remains unclear how ssRNA viruses adapt to novel environments and/or how many and what types of substitutions are needed to facilitate this evolution. In this study, we followed the adaptation of the ssRNA bacteriophage Qß using thermally adapted Escherichia coli as a host, which can efficiently grow at temperatures between 37.2 and 45.3 °C. This made it possible to evaluate Qß adaptation to the highest known temperature that supports growth, 45.3 °C. We found that Qß was capable of replication at this temperature; within 114 days (~1260 generations), we detected more than 34 novel point mutations in the genome of the thermally adapted Qß population, representing 0.8% of the total Qß genome. In addition, we returned the 45.3 °C-adapted Qß populations to 37.2 °C and passaged them for 8 days (~124 generations). We found that the reverse-adapted Qß population showed little to no decrease in fitness. These results indicate that Qß can evolve in response to increasing temperatures in a short period of time with the accumulation of point mutations.

Complete genomic sequence of Pseudomonas lactis bacteriophage HU1 isolated from raw cow's milk.

The lytic cold-active phage HU1, a member of the family Podoviridae, infects Pseudomonas lactis and was first isolated from raw cow's milk. In this study, we used deep sequencing to determine and analyze the DNA genome sequence of HU1. We identified a 42,551-base-pair genome comprising double-stranded DNA, with 69 predicted open reading frames and a GC content of 56.4%. A whole-genome comparison did not identify HU1 as a member of any previously reported cluster of Pseudomonas phages. By contrast, HU1 was most similar to AF, which infects P. putida, with nucleotide sequence alignment coverage of 24%. These results suggest that HU1 is a novel Pseudomonas phage.

A Lytic Bacteriophage for Controlling Pseudomonas lactis in Raw Cow's Milk.

The control of bacterial growth during milk processing is crucial for the quality maintenance of commercial milk and milk products. During a period of cold storage prior to heat treatments, some psychrotrophic bacteria grow and produce extracellular heat-resistant lipases and proteases that cause product defects. The use of lytic bacteriophages (phages) that infect and kill bacteria could be a useful tool for suppressing bacterial growth during this cold storage phase. In this study, we isolated a Pseudomonas lactis strain and a phage from raw cow's milk. Quantitative characterization of the phage was used to elucidate whether this phage was active under low temperatures and neutral pH and whether it was inactivated during pasteurization. Phage titer determination was possible under conditions ranging from pH 4 to 9 and from 3°C to 25°C; the phage was inactivated under pasteurization conditions at 63°C for 30 min. Furthermore, we showed that this phage reduced viable bacterial cell counts in both skim and whole milk. The results of this study represent the potential uses of phages for controlling psychrotrophic bacterial growth in raw cow's milk during cold storage.IMPORTANCE Suppression of bacterial growth in raw milk under cold storage is crucial for the quality control of commercially supplied milk. The use of lytic phages as low-cost microbicides is an attractive prospect. Due to strict host specificities, phages must be isolated from the raw milk where the host bacteria are growing. We first isolated the P. lactis bacterial strain and then the phage infecting that strain. Partial phage genomic analysis showed that this is a newly isolated phage, different from any previously reported. This study reports a lytic phage for P. lactis, and we have presented evidence here that this phage reduced viable bacterial cell counts not only in rich medium but also in skim and whole milk. As a result, we have concluded that the phage reported in this study would be useful in milk processing.

Influence of adaptive mutations, from thermal adaptation experiments, on the infection cycle of RNA bacteriophage Qß.

A population's growth rate is determined by multiple 'life history traits'. To quantitatively determine which life history traits should be improved to allow a living organism to adapt to an inhibitory environment is an important issue. Previously, we conducted thermal adaptation experiments on the RNA bacteriophage Qß using three independent replicates and reported that all three end-point populations could grow at a temperature (43.6°C) that inhibited the growth of the ancestral strain. Even though the fitness values of the endpoint populations were almost the same, their genome sequence was not, indicating that the three thermally adapted populations may have different life history traits. In this study, we introduced each mutation observed in these three end-point populations into the cDNA of the Qß genome and prepared three different mutants. Quantitative analysis showed that they tended to increase their fitness by increasing the adsorption rate to their host, shortening their latent period (i.e., the duration between phage infection and progeny release), and increasing the burst size (i.e., the number of progeny phages per infected cell), but all three mutants decreased their thermal stability. However, the degree to which these traits changed differed. The mutant with the least mutations showed a smaller decrease in thermal stability, the largest adsorption rate to the host, and the shortest latent period. These results indicated that several different adaptive routes exist by which Qß can adapt to higher temperatures, even though Qß is a simple RNA bacteriophage with a small genome size, encoding only four genes.

Characterization of a single mutation in TraQ in a strain of Escherichia coli partially resistant to Qß infection.

Bacteria and virulent bacteriophages are in a prey-predator relationship. Experimental models under simplified conditions with the presence of bacteria and bacteriophages have been used to elucidate the mechanisms that have enabled both prey and predator to coexist over long periods. In experimental coevolution conducted with Escherichia coli and the virulent RNA bacteriophage Qß in serial transfer, both coexisted for at least for 54 days, during which time they continued to change genetically and phenotypically. By day 16, an E. coli strain partially resistant to Qß appeared and caused an approximately 10(4)-fold decrease in Qß amplification. Whole-genome analysis of this strain suggested that a single mutation in TraQ was responsible for the partially resistant phenotype. TraQ interacts with propilin, encoded by the traA gene and a precursor of pilin, which is a component of the F pilus. The present study was performed to elucidate the mechanism underlying the coexistence of E. coli and Qß by investigating how a mutation in TraQ altered the physiological state of E. coli, and thus the amplification of Qß. Overexpression of wild-type TraQ in the partially resistant E. coli strain resulted in recovery of both TraA protein content, including propilin and pilin, and Qß amplification to levels comparable to those observed in the susceptible strain. Intriguingly, overexpression of the mutant TraQ in the partially resistant strains also increased the levels of TraA protein and Qß amplification, but these increases were smaller than those observed in the wild-type strain or the partially resistant strain expressing wild-type TraQ. The results of this study represent an example of how E. coli can become partially resistant to RNA bacteriophage infection via changes in a protein involved in maturation of a receptor rather than in the receptor itself and of how E. coli can stably coexist with virulent RNA bacteriophages.

Contribution of silent mutations to thermal adaptation of RNA bacteriophage Qß.

UNLABELLED: Changes in protein function and other biological properties, such as RNA structure, are crucial for adaptation of organisms to novel or inhibitory environments. To investigate how mutations that do not alter amino acid sequence may be positively selected, we performed a thermal adaptation experiment using the single-stranded RNA bacteriophage Qß in which the culture temperature was increased from 37.2°C to 41.2°C and finally to an inhibitory temperature of 43.6°C in a stepwise manner in three independent lines. Whole-genome analysis revealed 31 mutations, including 14 mutations that did not result in amino acid sequence alterations, in this thermal adaptation. Eight of the 31 mutations were observed in all three lines. Reconstruction and fitness analyses of Qß strains containing only mutations observed in all three lines indicated that five mutations that did not result in amino acid sequence changes but increased the amplification ratio appeared in the course of adaptation to growth at 41.2°C. Moreover, these mutations provided a suitable genetic background for subsequent mutations, altering the fitness contribution from deleterious to beneficial. These results clearly showed that mutations that do not alter the amino acid sequence play important roles in adaptation of this single-stranded RNA virus to elevated temperature. IMPORTANCE: Recent studies using whole-genome analysis technology suggested the importance of mutations that do not alter the amino acid sequence for adaptation of organisms to novel environmental conditions. It is necessary to investigate how these mutations may be positively selected and to determine to what degree such mutations that do not alter amino acid sequences contribute to adaptive evolution. Here, we report the roles of these silent mutations in thermal adaptation of RNA bacteriophage Qß based on experimental evolution during which Qß showed adaptation to growth at an inhibitory temperature. Intriguingly, four synonymous mutations and one mutation in the untranslated region that spread widely in the Qß population during the adaptation process at moderately high temperature provided a suitable genetic background to alter the fitness contribution of subsequent mutations from deleterious to beneficial at a higher temperature.

Adaptation of a cyanobacterium to a biochemically rich environment in experimental evolution as an initial step toward a chloroplast-like state.

Chloroplasts originated from cyanobacteria through endosymbiosis. The original cyanobacterial endosymbiont evolved to adapt to the biochemically rich intracellular environment of the host cell while maintaining its photosynthetic function; however, no such process has been experimentally demonstrated. Here, we show the adaptation of a model cyanobacterium, Synechocystis sp. PCC 6803, to a biochemically rich environment by experimental evolution. Synechocystis sp. PCC 6803 does not grow in a biochemically rich, chemically defined medium because several amino acids are toxic to the cells at approximately 1 mM. We cultured the cyanobacteria in media with the toxic amino acids at 0.1 mM, then serially transferred the culture, gradually increasing the concentration of the toxic amino acids. The cells evolved to show approximately the same specific growth rate in media with 0 and 1 mM of the toxic amino acid in approximately 84 generations and evolved to grow faster in the media with 1 mM than in the media with 0 mM in approximately 181 generations. We did not detect a statistically significant decrease in the autotrophic growth of the evolved strain in an inorganic medium, indicating the maintenance of the photosynthetic function. Whole-genome resequencing revealed changes in the genes related to the cell membrane and the carboxysome. Moreover, we quantitatively analyzed the evolutionary changes by using simple mathematical models, which evaluated the evolution as an increase in the half-maximal inhibitory concentration (IC50) and estimated quantitative characteristics of the evolutionary process. Our results clearly demonstrate not only the potential of a model cyanobacterium to adapt to a biochemically rich environment without a significant decrease in photosynthetic function but also the properties of its evolutionary process, which sheds light of the evolution of chloroplasts at the initial stage.

Quantitative comparison of the RNA bacteriophage Qß infection cycle in rich and minimal media.

As bacteriophages are dependent on the host for multiplication, their infection cycle is expected to be influenced by the host's physiological state. To elucidate how and which steps of the bacteriophage infection cycle are influenced by changes in the physiological state of the host, we quantitatively compared the infection cycle of lytic RNA bacteriophage Qß in Escherichia coli cultured in rich and minimal media. The adsorption rate constants in both media were almost the same. A difference of 15 min in the latent period and an approximately twofold increase in the rate of phage release were observed, although approximately 10(5) molecules of coat proteins, equivalent to approximately 600-1000 phage particles, accumulated in an infected cell prior to burst. Addition of Mg(2+) to minimal medium markedly affected the Qß infection cycle, and these results suggest that Mg(2+) is required for the stages of the infectious cycle after adsorption.

Ongoing phenotypic and genomic changes in experimental coevolution of RNA bacteriophage Qß and Escherichia coli.

According to the Red Queen hypothesis or arms race dynamics, coevolution drives continuous adaptation and counter-adaptation. Experimental models under simplified environments consisting of bacteria and bacteriophages have been used to analyze the ongoing process of coevolution, but the analysis of both parasites and their hosts in ongoing adaptation and counter-adaptation remained to be performed at the levels of population dynamics and molecular evolution to understand how the phenotypes and genotypes of coevolving parasite-host pairs change through the arms race. Copropagation experiments with Escherichia coli and the lytic RNA bacteriophage Qß in a spatially unstructured environment revealed coexistence for 54 days (equivalent to 163-165 replication generations of Qß) and fitness analysis indicated that they were in an arms race. E. coli first adapted by developing partial resistance to infection and later increasing specific growth rate. The phage counter-adapted by improving release efficiency with a change in host specificity and decrease in virulence. Whole-genome analysis indicated that the phage accumulated 7.5 mutations, mainly in the A2 gene, 3.4-fold faster than in Qß propagated alone. E. coli showed fixation of two mutations (in traQ and csdA) faster than in sole E. coli experimental evolution. These observations suggest that the virus and its host can coexist in an evolutionary arms race, despite a difference in genome mutability (i.e., mutations per genome per replication) of approximately one to three orders of magnitude.

Adaptation by stochastic switching of a monostable genetic circuit in Escherichia coli.

Stochastic switching is considered as a cost-saving strategy for adaptation to environmental challenges. We show here that stochastic switching of a monostable circuit can mediate the adaptation of the engineered OSU12-hisC Escherichia coli strain to histidine starvation. In this strain, the hisC gene was deleted from the His operon and placed under the control of a monostable foreign promoter. In response to histidine depletion, the OSU12-hisC population shifted to a higher HisC expression level, which is beneficial under starving conditions but is not favoured by the monostable circuit. The population shift was accompanied by growth recovery and was reversible upon histidine addition. A weak directionality in stochastic switching of hisC was observed in growing microcolonies under histidine-free conditions. Directionality and fate decision were in part dependent on the initial cellular status. Finally, microarray analysis indicated that OSU12-hisC reorganized its transcriptome to reach the appropriate physiological state upon starvation. These findings suggest that bacteria do not necessarily need to evolve signalling mechanisms to control gene expression appropriately, even for essential genes.

Cooperative adaptation to establishment of a synthetic bacterial mutualism.

To understand how two organisms that have not previously been in contact can establish mutualism, it is first necessary to examine temporal changes in their phenotypes during the establishment of mutualism. Instead of tracing back the history of known, well-established, natural mutualisms, we experimentally simulated the development of mutualism using two genetically-engineered auxotrophic strains of Escherichia coli, which mimic two organisms that have never met before but later establish mutualism. In the development of this synthetic mutualism, one strain, approximately 10 hours after meeting the partner strain, started oversupplying a metabolite essential for the partner's growth, eventually leading to the successive growth of both strains. This cooperative phenotype adaptively appeared only after encountering the partner strain but before the growth of the strain itself. By transcriptome analysis, we found that the cooperative phenotype of the strain was not accompanied by the local activation of the biosynthesis and transport of the oversupplied metabolite but rather by the global activation of anabolic metabolism. This study demonstrates that an organism has the potential to adapt its phenotype after the first encounter with another organism to establish mutualism before its extinction. As diverse organisms inevitably encounter each other in nature, this potential would play an important role in the establishment of a nascent mutualism in nature.

Probabilistic transition from unstable predator-prey interaction to stable coexistence of Dictyostelium discoideum and Escherichia coli.

Predator-prey interactions have been found at all levels within ecosystems. Despite their ecological ubiquity and importance, the process of transition to a stable coexistent state has been poorly verified experimentally. To investigate the stabilization process of predator-prey interactions, we previously constructed a reproducible experimental predator-prey system between Dictyostelium discoideum and Escherichia coli, and showed that the phenotypically changed E. coli contributed to stabilization of the system. In the present study, we focused on the transition to stable coexistence of both species after the phenotypic change in E. coli. Analysis of E. coli cells isolated from co-culture plates as single colony enabled us to readily identify the appearance of phenotypically changed E. coli that differed in colony morphology and growth rate. It was also demonstrated that two types of viscous colony, i.e., the dense-type and sparse-type, differing in spatial distribution of both species emerged probabilistically and all of the viscous colonies maintained stably were of the sparse-type. These results suggest that the phenotypically changed E. coli may produce two types of viscous colonies probabilistically. The difference in spatial distribution would affect localized interactions between both species and then cause probabilistic stabilization of predator-prey interactions.

Single-cell isolation and cloning of Tetrahymena thermophila cells with a fluorescence-activated cell sorter.

We developed a method for cloning cells of the ciliate Tetrahymena thermophila in chemically defined medium (CDM) using a fluorescence-activated cell sorter (FACS). Although T. thermophila is a model unicellular eukaryote, two major technical difficulties remain in its cloning. First, T. thermophila fails to proliferate from low density in CDM, particularly if the inoculum contains single cells. Second, general cloning methods are time consuming and have low throughput. Here, we modified the CDM by addition of bovine serum albumin that helped growth from an inoculum with a density of 10 cell/ml (1 cell/100 µl). In addition, we applied a FACS for isolation of single cells. We showed that it is possible to separate cell populations based on the presence or absence of phagocytosed fluorescent beads and to isolate single cells in a modified CDM by FACS. Our techniques allow the direct isolation of single cells and facilitate the establishment of clonal strains.