Genomic sequencing of single microbial cells from environmental samples

Recently developed techniques allow genomic DNA sequencing from single microbial cells [Lasken RS: Single-cell genomic sequencing using multiple displacement amplification. Curr Opin Microbiol 2007, 10:510-516]. Here, we focus on research strategies for putting these methods into practice in the laboratory setting. An immediate consequence of single-cell sequencing is that it provides an alternative to culturing organisms as a prerequisite for genomic sequencing. The microgram amounts of DNA required as template are amplified from a single bacterium by a method called multiple displacement amplification (MDA) avoiding the need to grow cells. The ability to sequence DNA from individual cells will likely have an immense impact on microbiology considering the vast numbers of novel organisms, which have been inaccessible unless culture-independent methods could be used. However, special approaches have been necessary to work with amplified DNA. MDA may not recover the entire genome from the single copy present in most bacteria. Also, some sequence rearrangements can occur during the DNA amplification reaction. Over the past two years many research groups have begun to use MDA, and some practical approaches to single-cell sequencing have been developed. We review the consensus that is emerging on optimum methods, reliability of amplified template, and the proper interpretation of 'composite' genomes which result from the necessity of combining data from several single-cell MDA reactions in order to complete the assembly. Preferred laboratory methods are considered on the basis of experience at several large sequencing centers where >70% of genomes are now often recovered from single cells. Methods are reviewed for preparation of bacterial fractions from environmental samples, single-cell isolation, DNA amplification by MDA, and DNA sequencing.     

Chitinases from uncultured marine microorganisms.

Our understanding of the degradation of organic matter will benefit from a greater appreciation for the genes encoding enzymes involved in the hydrolysis of biopolymers such as chitin, one of the most abundant polymers in nature. To isolate representative and abundant chitinase genes from uncultivated marine bacteria, we constructed libraries of genomic DNA isolated from coastal and estuarine waters. The libraries were screened for genes encoding proteins that hydrolyze a fluorogenic analogue of chitin, 4-methylumbelliferyl beta-D-N,N'-diacetylchitobioside (MUF-diNAG). The abundance of clones capable of MUF-diNAG hydrolysis was higher in the library constructed with DNA from the estuary than in that constructed with DNA from coastal waters, although the abundance of positive clones was also dependent on the method used to screen the library. Plaque assays revealed nine MUF-diNAG-positive clones of 75,000 screened for the estuarine sample and two clones of 750,000 for the coastal sample. A microtiter plate assay revealed approximately 1 positive clone for every 500 clones screened in the coastal library. The number of clones detected with the plaque assay was consistent with estimates of the portion of culturable bacteria that degrade chitin. Our results suggest that culture-dependent methods do not greatly underestimate the portion of marine bacterial communities capable of chitin degradation.     

Screening for novel lipolytic enzymes from uncultured soil microorganisms

The construction and screening of metagenomic libraries constitute a valuable resource for obtaining novel biocatalysts. In this work, we present the construction of a metagenomic library in Escherichia coli using fosmid and microbial DNA directly isolated from forest topsoil and screened for lipolytic enzymes. The library consisted of 33,700 clones with an average DNA insert size of 35 kb. Eight unique lipolytic active clones were obtained from the metagenomic library on the basis of tributyrin hydrolysis. Subsequently, secondary libraries in a high-copy-number plasmid were generated to select lipolytic subclones and to characterize the individual genes responsible for the lipolytic activity. DNA sequence analysis of six genes revealed that the enzymes encoded by the metagenomic genes for lipolytic activity were novel with 34–48% similarity to known enzymes. They had conserved sequences similar to those in the hormone-sensitive lipase family. Based on their deduced amino acid similarity, the six genes encoding lipolytic enzymes were further divided into three subgroups, the identities among which ranged from 33% to 45%. The six predicted gene products were successfully expressed in E. coli and secreted into the culture broth. Most of the secreted enzymes showed a catalytic activity for hydrolysis of p-nitrophenyl butyrate (C₄) but not p-nitrophenyl palmitate (C₁₆).     

Molecular cloning and functional characterization of a novel decarboxylase from uncultured microorganisms

The metagenomic library approach has been used successfully to isolate novel biocatalyst genes from uncultured microorganisms. We report the cloning of a novel decarboxylase gene by sequence-based screening of a plasmid metagenomic library constructed with DNA from alkaline polluted soils. The gene was named undec1 A and had an open reading frame of 1077 base pairs. It encoded a 359 amino acid polypeptide with a molecular mass of 38 kDa. The predicted protein had 58% similarity to a decarboxylase from Chlorobium phaeobacteroides BS1. The putative decarboxylase gene was subcloned into pETBlue-2 vector and overexpressed in Escherichia coli Tuner (DE3) pLac. The recombinant protein was purified to homogeneity. Functional characterization with liquid chromatography-mass spectrometry confirmed that the recombinant Undec1 A protein catalyzed the decarboxylation of L-cysteine to form cysteamine.     

Isolation of a gene encoding endoglucanase activity from uncultured microorganisms in buffalo rumen

A gene, umcel5N, was isolated from a metagenomic library constructed from the contents of buffalo rumen. Its putative product belongs to the glycosyl hydrolase family 5 and is most closely related to an endoglucanase (ABN54006.1) from Clostridium thermocellum with 44% identity and 60% similarity. Gene umcel5N was heterologously expressed in Escherichia coli. The purified recombinant Umcel5N hydrolyzed carboxymethyl cellulose with a rapid decrease in the viscosity of the solution but with little release of reducing sugars, suggesting an endo mode of action. The enzyme exhibited optimal activity toward p-nitrophenyl β-d-cellobioside at pH 5.5 and 55°C, and had a Km of 1.56 mM and a Vmax of 285.6 U/mg. Two glutamic acids (E144 and E285) of the wild-type Umcel5N were predicted as a proton donor and a nucleophile, respectively. Site-directed mutagenesis confirmed that they were required for the enzyme’s activity.     

Specific single-cell isolation and genomic amplification of uncultured microorganisms

We in this study describe a new method for genomic studies of individual uncultured prokaryotic organisms, which was used for the isolation and partial genome sequencing of a soil archaeon. The diversity of Archaea in a soil sample was mapped by generating a clone library using group-specific primers in combination with a terminal restriction fragment length polymorphism profile. Intact cells were extracted from the environmental sample, and fluorescent in situ hybridization probing with Cy3-labeled probes designed from the clone library was subsequently used to detect the organisms of interest. Single cells with a bright fluorescent signal were isolated using a micromanipulator and the genome of the single isolated cells served as a template for multiple displacement amplification (MDA) using the Phi29 DNA polymerase. The generated MDA product was afterwards used for 16S rRNA gene sequence analysis and shotgun-cloned for additional genomic analysis. Sequence analysis showed >99% 16S rRNA gene homology to soil crenarchaeotal clone SCA1170 and shotgun fragments had the closest match to a crenarchaeotal BAC clone previously retrieved from a soil sample. The system was validated using Methanothermobacter thermoautotrophicus as single-cell test organism, and the validation setup produced 100% sequence homology to the ten tested regions of the genome of this organism.     

Genomic sequencing of key genes in mouse pancreatic cancer cells

Pancreatic cancer is a multiple genetic disorder with many mutations identified during the progression. Two mouse pancreatic cancer cell lines were established which showed different phenotype in vivo: a non-metastatic cell line, Panc02, and a highly metastatic cell line, Panc02-H7, a derivative of Panc02. In order to investigate whether the genetic mutations of key genes in pancreatic cancer such as KRAS, TP53 (p53), CDKN2A (p16), SMAD4, ZIP4, and PDX-1 contribute to the phenotypic difference of these two mouse pancreatic cancer cells, we sequenced the exonic regions of these key genes in both cell lines and in the normal syngeneic mouse pancreas and compared them with the reference mouse genome sequence. The exons of KRAS, SMAD4, CDKN2A (p16), TP53 (p53), ZIP4, and PDX-1 genes were amplified and the genotype of these genes was determined by Sanger sequencing. The sequences were analyzed with Sequencher software. A mutation in SMAD4 was identified in both cell lines. This homozygote G to T mutation in the first position of codon 174 (GAA) generated a stop codon resulting in the translation of a truncated protein. Further functional analysis indicates that different TGF-β/SMAD signaling pathways were involved in those two mouse cell lines, which may explain the phonotypic difference between the two cells. A single nucleotide polymorphism (SNP) in KRAS gene (TAT to TAC at codon 32) was also identified in the normal pancreas DNA of the syngenic mouse and in both derived tumoral Panc02 and Panc02-H7 cells. No mutation or SNP was found in CDKN2A (p16), TP53 (p53), ZIP4, and PDX-1 genes in these two cell lines. The absence of mutations in genes such as KRAS, TP53, and CDKN2A, which are considered as key genes in the development of human pancreatic cancer suggests that SMAD4 might play a central and decisive role in mouse pancreatic cancer. These results also suggest that other mechanisms are involved in the substantial phenotypic difference between these two mouse pancreatic cancer cell lines. Further studies are warranted to elucidate the molecular pathways that lead to the aggressive metastatic potential of Panc02-H7.     

Genomic content of uncultured Bacteroidetes from contrasting oceanic provinces in the North Atlantic Ocean

Bacteroidetes are widespread in marine systems where they play a crucial role in organic matter degradation. Whole genome analysis of several strains has revealed a broad glycolytic and proteolytic potential. In this study, we used a targeted metagenomic approach to investigate the degradation capabilities of distinct Bacteroidetes clades from two contrasting regions of the North Atlantic Ocean, the Polar Biome (BPLR) and the North Atlantic Subtropical (NAST). We present here the analysis of 76 Bacteroidetes fosmids, of which 28 encode the 16S rRNA gene as phylogenetic marker, and their comparison to complete Bacteroidetes genomes. Almost all of the 16S rRNA harbouring fosmids belonged to clades that we previously identified in BPLR and NAST. The majority of sequenced fosmids could be assigned to Bacteroidetes affiliated with the class Flavobacteria. We also present novel genomic information on the classes Cytophagia and Sphingobacteria, suggesting a capability of the latter for attachment to algal surfaces. In our fosmid set we identified a larger potential for polysaccharide degradation and cell surface attachment in the phytoplankton-rich BPLR. Particularly, two flavobacterial fosmids, one affiliated with the genus Polaribacter, showed a whole armoury of enzymes that likely function in degradation of sulfated polysaccharides known to be major constituents of phytoplankton cell walls. Genes involved in protein and peptidoglycan degradation, although present in both fosmid sets, seemed to have a slight preponderance in NAST. This study provides support for the hypothesis of a distinct specialization among marine Bacteroidetes for the degradation of certain types of polymers.     

Cloning and identification of cellulase genes from uncultured microorganisms in pulp sediments from paper mill effluent]

The metagenomic DNA of pulp sediments from paper mill effluent was extracted and purified. The 16S rDNA was amplified using the purified metagenomic DNA as template and a 16S rDNA library was prepared. Sequence analysis of 16S rDNA clones showed that diverse of uncultured bacteria inhabit in this environment, which can be classified into 4 clusters as Spirochaetes, Proteobacteria, Bacteroidetes and Firmicutes. A metagenomic library containing 10000 clones was constructed into cosmid vector, and the capacity of inserted DNA of which was 3.53 x 10(8) bp. Functional screening of the library resulted in isolation of two independent clones expressing endoglucanase activity, three independent clones expressing exoglucanase activity and two independent clones expressing beta-glucosidase activity. One clone expressing strongest enzyme activity from each activity category was chosen to be further analyzed. Three novel cellulase genes designated as umcel5L, umcel5M and umbgl3D were identified by subcloning, sequencing and expression. The umcel5L encodes an endoglucanase belonging to glycosyl hydrolase family 5, which is most related to an endoglucanase from Bradyrhizobium japonicum at 43% identity and 59% similarity. The umcel5M encodes a cellodextrinase belonging to glycosyl hydrolase family 5, which is most similar to a cellodextrinase from Fibrobacter succinogenes at 48% identity and 69% similarity. The umbgl3D encodes a putative beta-glucosidase belonging to glycosyl hydrolase family 3, which shares highest homology with a beta-glucosidase from Thermotoga maritima at 46% identity and 61% similarity. It is the first time to reveal the bacterial diversity of pulp sediments from paper mill effluent and clone novel cellulase genes from the bacteria by culture-independent method.     

Genomic DNA amplification from a single bacterium

Genomic DNA was amplified about 5 billion-fold from single, flow-sorted bacterial cells by the multiple displacement amplification (MDA) reaction, using phi 29 DNA polymerase. A 662-bp segment of the 16S rRNA gene could be accurately sequenced from the amplified DNA. MDA methods enable new strategies for studying non-culturable microorganisms.     

Genomic sequencing by ligation-mediated PCR

Genomic sequencing permits studies of in vivo DNA methylation and protein-DNA interactions, but its use has been limited due to the complexity of the mammalian genome. Ligation-mediated PCR (LMPCR) is a sensitive genomic sequencing procedure that generates high quality, reproducible sequence ladders starting with only 1 μg of uncloned mammalian DNA per reaction. This genomic sequencing procedure can be adapted for various methylation, in vivo footprinting and DNA adduct mapping procedures. We provide a detailed protocol for genomic sequencing by LMPCR and discuss the principles and applications of the method.     

Genomic sequencing: the complexity conundrum

The concept of 'organismal complexity' has had a chequered career in genetics, with no rigorous operational definition available for the term. The recent finding that Drosophila melanogaster has more than four thousand fewer genes than the nematode forces a re-examination of whether gene number, in itself, can be taken as any real guide to complexity.     

Genomic and phenotypic analyses of microorganisms isolated from the sediments of Lake Baikal

Genetic and biochemical methods and morphological examination were used to study microorganisms isolated from samples of deep drilling of the Lake Baikal bottom sediments. Based on blot hybridization patterns, the strains investigated were divided into several groups according to the degree of homology of their genomic DNA. Morphological, biochemical, and ultrastructural characteristics of bacterial strains are described, and their compliance with the genomic analysis data is demonstrated.     

PCR and sequencing from a single pollen grain

In order to eliminate the laborious step of DNA extraction preceding all studies within the field of plant molecular biology we attempted to do PCR amplifications directly on pollen grains. Successful PCR amplification was obtained in reactions including a single pollen grain from Hordeum vulgare or Secale strictum. Both the plastid gene encoding ribulose-1,5-biphosphate carboxylase/oxygenase (rbcL) and the nuclear-encoded internal transcribed spacer regions (ITS) and the 5.8S rDNA region were amplified and sequenced to verify PCR amplification.     

Automated identification of single nucleotide polymorphisms from sequencing data

The single nucleotide polymorphism (SNP) is the difference of the DNA sequence between individuals and provides abundant information about genetic variation. Large scale discovery of high frequency SNPs is being undertaken using various methods. However, the publicly available SNP data sometimes need to be verified. If only a particular gene locus is concerned, locus-specific polymerase chain reaction amplification may be useful. Problem of this method is that the secondary peak has to be measured. We have analyzed trace data from conventional sequencing equipment and found an applicable rule to discern SNPs from noise. The rule is applied to multiply aligned sequences with a trace and the peak height of the traces are compared between samples. We have developed software that integrates this function to automatically identify SNPs. The software works accurately for high quality sequences and also can detect SNPs in low quality sequences. Further, it can determine allele frequency, display this information as a bar graph and assign corresponding nucleotide combinations. It is also designed for a person to verify and edit sequences easily on the screen. It is very useful for identifying de novo SNPs in a DNA fragment of interest.