Whole genome amplification of single cells: mathematical analysis of PEP and tagged PCR

Single-Cell Analysis in the Omics Era: Technologies and Applications in Cancer

Cancer molecular profiling obtained with conventional bulk sequencing describes average alterations obtained from the entire cellular population analyzed. In the era of precision medicine, this approach is unable to track tumor heterogeneity and cannot be exploited to unravel the biological processes behind clonal evolution. In the last few years, functional single-cell omics has improved our understanding of cancer heterogeneity. This approach requires isolation and identification of single cells starting from an entire population. A cell suspension obtained by tumor tissue dissociation or hematological material can be manipulated using different techniques to separate individual cells, employed for single-cell downstream analysis. Single-cell data can then be used to analyze cell-cell diversity, thus mapping evolving cancer biological processes. Despite its unquestionable advantages, single-cell analysis produces massive amounts of data with several potential biases, stemming from cell manipulation and pre-amplification steps. To overcome these limitations, several bioinformatic approaches have been developed and explored. In this work, we provide an overview of this entire process while discussing the most recent advances in the field of functional omics at single-cell resolution.

Confidence intervals for Markov chain transition probabilities based on next generation sequencing reads data

Background

Markov chains (MC) have been widely used to model molecular sequences. The estimations of MC transition matrix and confidence intervals of the transition probabilities from long sequence data have been intensively studied in the past decades. In next generation sequencing (NGS), a large amount of short reads are generated. These short reads can overlap and some regions of the genome may not be sequenced resulting in a new type of data. Based on NGS data, the transition probabilities of MC can be estimated by moment estimators. However, the classical asymptotic distribution theory for MC transition probability estimators based on long sequences is no longer valid.

Methods

In this study, we present the asymptotic distributions of several statistics related to MC based on NGS data. We show that, after scaling by the effective coverage d defined in a previous study by the authors, these statistics based on NGS data approximate to the same distributions as the corresponding statistics for long sequences.

Results

We apply the asymptotic properties of these statistics for finding the theoretical confidence regions for MC transition probabilities based on NGS short reads data. We validate our theoretical confidence intervals using both simulated data and real data sets, and compare the results with those by the parametric bootstrap method.

Conclusions

We find that the asymptotic distributions of these statistics and the theoretical confidence intervals of transition probabilities based on NGS data given in this study are highly accurate, providing a powerful tool for NGS data analysis.

Whole Genome Amplification and Branching Processes

Whole genome amplification is important for multipoint mapping by sperm or oocyte typing and genetic disease diagnosis. Polymerase chain reaction is not suitable for amplifying long DNA sequences. This paper studies a new technique, designated PEP-primer-extension-preamplification, for amplifying long DNA sequences using the theory of branching processes. A mathematical model for PEP is constructed and a closed formula for the expected target yield is obtained. A central limit theorem and a strong law of large numbers for the number of kth generation target sequences are proved.

Whole Genome Amplification and Branching Processes

Whole genome amplification is important for multipoint mapping by sperm or oocyte typing and genetic disease diagnosis. Polymerase chain reaction is not suitable for amplifying long DNA sequences. This paper studies a new technique, designated PEP-primer-extension-preamplification, for amplifying long DNA sequences using the theory of branching processes. A mathematical model for PEP is constructed and a closed formula for the expected target yield is obtained. A central limit theorem and a strong law of large numbers for the number of kth generation target sequences are proved.

Modeling genome coverage in single-cell sequencing

MOTIVATION

Single-cell DNA sequencing is necessary for examining genetic variation at the cellular level, which remains hidden in bulk sequencing experiments. But because they begin with such small amounts of starting material, the amount of information that is obtained from single-cell sequencing experiment is highly sensitive to the choice of protocol employed and variability in library preparation. In particular, the fraction of the genome represented in single-cell sequencing libraries exhibits extreme variability due to quantitative biases in amplification and loss of genetic material.

RESULTS

We propose a method to predict the genome coverage of a deep sequencing experiment using information from an initial shallow sequencing experiment mapped to a reference genome. The observed coverage statistics are used in a non-parametric empirical Bayes Poisson model to estimate the gain in coverage from deeper sequencing. This approach allows researchers to know statistical features of deep sequencing experiments without actually sequencing deeply, providing a basis for optimizing and comparing single-cell sequencing protocols or screening libraries.

AVAILABILITY AND IMPLEMENTATION

The method is available as part of the preseq software package. Source code is available at http://smithlabresearch.org/preseq.

CONTACT

andrewds@usc.edu

SUPPLEMENTARY INFORMATION

Supplementary material is available at Bioinformatics online.

The future is now: single-cell genomics of bacteria and archaea

Interest in the expanding catalog of uncultivated microorganisms, increasing recognition of heterogeneity among seemingly similar cells, and technological advances in whole-genome amplification and single-cell manipulation are driving considerable progress in single-cell genomics. Here, the spectrum of applications for single-cell genomics, key advances in the development of the field, and emerging methodology for single-cell genome sequencing are reviewed by example with attention to the diversity of approaches and their unique characteristics. Experimental strategies transcending specific methodologies are identified and organized as a road map for future studies in single-cell genomics of environmental microorganisms. Over the next decade, increasingly powerful tools for single-cell genome sequencing and analysis will play key roles in accessing the genomes of uncultivated organisms, determining the basis of microbial community functions, and fundamental aspects of microbial population biology.

A short and simple improved-primer extension preamplification (I-PEP) procedure for whole genome amplification (WGA) of bovine cells

Embryo transfer is a reproductive technique that has a major impact on the dissemination of economically important genes and the rate of genetic gain in breeding schemes. In recent years, there has been increasing interest in the use of sexed and genotyped embryos in commercial embryo transfer programs. Marker/gene assisted selection (MAS/GAS) projects can be performed in the pre-implantation stage through mass production of characterized embryos. Biopsy of a few cells in the morulla stage is essential for pre-implantation genetic diagnosis (PGD), in which sex determination, evaluation of disease genes, and genotyping for candidate genes are performed. Limited quantity of cells and low amount of DNA restrict the use of multiple molecular analyses in PGD programs. Recently, whole genome amplification (WGA) techniques promise to overcome this problem by providing sufficient input DNA for analysis. Among several techniques proposed for WGA, the primer extension pre-amplification (PEP) and the improved-primer extension pre-amplification (I-PEP) methods are the most commonly used. However, these methods are time-consuming and need more than 12 h amplification cycles. Since the time is a critical parameter in the successful characterized embryo transfer, the shortening of diagnosis time is highly desirable. In this study, we developed a short and simple I-PEP procedure (~3 h) and evaluated its performance for the amplification of bovine genomic DNA. We assessed short WGA procedure by polymerase chain reaction (PCR) amplification of 7 specific loci. The results indicated that the short procedure possesses enough sensitivity for the molecular genetic analysis of 1 input cell. Although the efficiency of the method was 100%, there was an inconsistency between genomic DNA (gDNA) and whole genome amplification product (wgaDNA) genotypes for kappa-casein locus; that is, however, most likely due to allele drop-out (ADO) or false homozigocity. The results of this study indicate that with the application of reliable methods, WGA-amplified bovine DNA will be a useful source for sexing and genotyping bovine embryos in several quantitative trait locus (QTL) markers.

dcDegenerate oligonucleotide primed-PCR for multilocus, genome-wide analysis from limited quantities of DNA

This study modified the degenerate oligonucleotide primed-polymerase chain reaction (DOP-PCR)-based whole genome amplification method for improvement of downstream genome-wide analysis of low copy number DNA samples (<or= 0.100 ng). Experiments involved altering the degeneracy of the DOP primer, nonspecific cycle number, and adding proofreading polymerases. Increasing the degeneracy of the primer and the number of cycles that use a low annealing temperature should improve the nonspecific amplification of the DOP-PCR reaction. The addition of proofreading enzymes should allow for longer amplification products, increasing the genome coverage of the reaction. Low-input DNA quantities were examined for the primer and the cycle number studies using standard DOP-PCR parameters. The optimized DOP-PCR technique was then implemented for the polymerase study. All DOP-PCR products were amplified by using a multiplex microsatellite amplification kit to evaluate products from multiple chromosomes, followed by separation and detection by capillary electrophoresis. The 10 N primer, 12 nonspecific cycles, and the addition of the DeepVent proofreading enzyme all significantly increased the number of short tandem repeat alleles successfully amplified. All modifications also lowered the rate of allele drop-in, or sporadic additional allele occurrence, when compared with DOP-PCR results published earlier. Further, an average of > 0.50 intralocus heterozygote peak ratios were observed for most DNA input quantities examined. These results show that modifications of the traditional DOP-PCR reaction (dcDOP-PCR) to include the use of a more degenerate primer (10 N), 12 nonspecific cycles, and a proofreading enzyme allows for a more complete, balanced chromosome amplification from limited and/or compromised clinical and biological samples.

Whole-Genome Amplification by Improved Primer Extension Preamplification PCR (I-PEP-PCR)

INTRODUCTIONPCR-based whole-genome amplification (WGA) has the goal of generating microgram quantities of genome-representative DNA from picogram or nanogram amounts of starting material. This amplification should introduce little, or ideally no, representational bias. In contrast to other techniques for WGA, PCR-based methods are generally less affected by DNA quality and are more applicable to DNA extracted from various sources (fixed and fresh tissues). Primer extension preamplification PCR (PEP-PCR), in contrast to degenerate oligonucleotide primed PCR (DOP-PCR), uses totally degenerate 15-mer PCR primers. An additional difference is that in PEP-PCR, the number of potential priming sites is orders of magnitude larger. The effectiveness of PEP-PCR has been increased by several alterations. The improved PEP (I-PEP) PCR approach, described in this protocol, uses a DNA polymerase cocktail that includes Taq DNA polymerase (to carry out the primer extension as in a traditional PCR) and a proofreading DNA polymerase (to provide 3'-to-5'-exonuclease activity, excising misincorporated nucleotides that slow the progression of Taq DNA polymerase). The result is far more efficient WGA, with increased fidelity due to the removal of the misincorporated nucleotides. Similar to DOP-PCR, PEP-PCR generates a smear of DNA fragments that are visible on an agarose gel.

Prediction of recurrence of nonfunctioning pituitary tumours by loss of heterozygosity analysis

BACKGROUND

Postsurgical regrowth or recurrence of nonfunctioning pituitary adenomas (NFAs) is not uncommon and often requires further surgery or radiotherapy (DXT). Routine postoperative DXT increases the incidence of hypopituitarism, which is associated with increased morbidity and mortality. Identification of genetic abnormalities in the tumour tissue, which can predict recurrence, may allow targeting DXT to the most appropriate patients.

DESIGN AND METHODS

We have performed loss of heterozygosity (LOH) analysis on 96 NFAs of which 43 (45%) were recurrent and 53 (55%) were nonrecurrent tumours. Analysis of all tumours was performed on the surgical specimen obtained at the time of first surgery. All tumours underwent allelotyping across nine highly informative microsatellite markers selected on the basis of high LOH frequency in an earlier study involving genome-wide allelotyping. LOH frequency across all microsatellite markers as well as across individual markers was compared between the two cohorts of tumours.

RESULTS

LOH frequency in tumours that subsequently recurred was significantly higher across all microsatellite markers as compared to tumours that did not recur (P < 0.05). Allelic loss across one or more microsatellite marker was significantly higher in recurrent tumours (30/43) as compared to their nonrecurrent counterparts (17/53) (P < 0.01). On Poisson regression analysis, the higher LOH frequency in recurrent tumours was independent of the invasiveness of tumours determined radiologically. In addition, LOH at the microsatellite markers D1S215 and D1S459 was significantly higher in tumours that recurred as compared to tumours that did not (32%vs. 3% and 27%vs. 2%, respectively; P < 0.01 for both). No significant difference in LOH frequency between the two tumour groups was evident at the other markers. No association could be demonstrated between the frequency and pattern of LOH and the time to manifest recurrence.

CONCLUSIONS

We have shown that it may be possible to predict recurrence of NFAs by LOH analysis of the initial tumour specimen at predefined microsatellite markers, especially on chromosome 1q. This merits further prospective study.

Whole genome amplification of DNA from laser capture-microdissected tissue for high-throughput single nucleotide polymorphism and short tandem repeat genotyping

Genome-wide screening of genetic alterations between normal and cancer cells, as well as among subgroups of tumors, is important for establishing molecular mechanism and classification of cancer. Gene silencing through loss of heterozygosity is widely observed in cancer cells and detectable by analyzing allelic loss of single nucleotide polymorphism and/or short tandem repeat markers. To use minute quantities of DNA that are available through laser capture microdissection (LCM) of cancer cells, a whole genome amplification method that maintains locus and allele balance is essential. We have successfully used a ø29 polymerase-based isothermal whole genome amplification method to amplify LCM DNA using a proteinase K lysis procedure coupled with a pooling strategy. Through single nucleotide polymorphism and short tandem repeat genotype analysis we demonstrate that using pooled DNA from two or three separate amplification reactions significantly reduces any allele bias introduced during amplification. This strategy is especially effective when using small quantities of source DNA. Although a convenient alkaline lysis DNA extraction procedure provided satisfactory results from using 1500 to 3000 LCM cells, proteinase K digestion was superior for lower cell numbers. Accurate genotyping is achieved with as few as 100 cells when both proteinase K extraction and pooling are applied.

Laser microdissection of small tissue samples--application to chronic pancreatitis tissues

Laser microdissection is considered to be the gold standard of tissue sampling, especially if a defined small tissue area consisting of single or few cells within a heterogeneous tissue compartment is of interest. This sophisticated technique offers the opportunity of rapid and contamination-free tissue sampling for RNA- or DNA-based molecular genetic studies. We have applied laser microdissection to a molecular genetic study of pancreatic intraductal lesions (PanINs) in tissues of chronic pancreatitis, where an exact microdissection of small ducts within a dense fibrous tissue is of paramount importance for following analysis. From nine patients suffering from chronic pancreatitis, formalin-fixed, paraffin-embedded tissue specimens were laser microdissected, and a total of 202 normal ducts and PanINs of grade PanIN-1A to grade PanIN-2 were harvested. After whole genome amplification by improved primer extension and preamplification PCR (I-PEP-PCR), microsatellite-PCR based loss of heterozygosity analysis (LOH) of the tumor suppressor gene loci TP53, p16INK4, and DPC4 was performed. One of 85 informative duct lesions (1.2%) had LOH of TP53, 1 of 76 duct lesions (1.3%) had LOH of DPC4, and 2/29 duct lesions (6.9%) showed LOH of p16INK4. Microsatellite instability (MSI) was seen in 2 of 178 duct lesions (1.1%). Immunohistochemical staining of p53 protein and DPC4 protein revealed no aberrant expression. These preliminary data indicate that LOH of tumor suppressor genes, important in pancreatic cancer genesis or MSI, can be found in chronic pancreatitis tissues, but their incidence is low.

Genome-wide amplification and allelotyping of sporadic pituitary adenomas identify novel regions of genetic loss

Through the use of a candidate gene approach, several previous studies have identified loss of heterozygosity (LOH) at putative tumor-suppressor gene (TSG) loci in sporadic pituitary tumors. This study reports a genome-wide allelotyping by use of 122 microsatellite markers in a large cohort of tumors, consisting of somatotrophinomas and non-functioning adenomas. Samples were first subject to prior whole genome amplification by primer extension pre-amplification (PEP) to circumvent limitations imposed by insufficient DNA for whole-genome analysis with this number of microsatellite markers. The overall mean frequency of loss in invasive tumors was significantly higher than that in their non-invasive counterparts (7 vs. 3% somatotrophinomas; 6 vs. 3% non-functioning adenomas, respectively). Analysis of the mean frequency of LOH, across all markers to individual chromosomal arms, identified 13 chromosomal arms in somatotrophinomas and 10 in non-functioning tumors, with LOH greater than the 99% upper confidence interval calculated for the rate of overall random allelic loss. In the majority of cases, these losses were more frequent in invasive tumors than in their non-invasive counterparts, suggesting these to be markers of tumor progression. Other regions showed similar frequencies of LOH in both invasive and non-invasive tumors, implying these to be early changes in pituitary tumorigenesis. This genome-wide study also revealed chromosomal regions where losses were frequently associated with an individual marker, for example, chromosome arm 1q (LOH > 30%). In some cases, these losses were subtype-specific and were found at a higher frequency in invasive tumors than in their non-invasive counterparts. Identification of these regions of loss provides the first preliminary evidence for the location of novel putative TSGs involved in pituitary tumorigenesis that are, in some cases, subtype-specific. This investigation provides an unbiased estimate of global aberrations in sporadic pituitary tumors as assessed by LOH analysis. The identification of multiple "hotspots" throughout the genome may be a reflection of an unstable chromatin structure that is susceptible to a deletion or epigenetic-mediated gene-silencing events.

Progress on methods of gene detection in preimplantation embryos

The advent of the polymerase chain reaction (PCR) and the development of fluorescence in situ hybridization (FISH) have had a tremendous impact on preimplantation genetic diagnosis (PGD). While PCR is a powerful tool in detecting genetic diseases or molecular markers affecting quantitative trait loci, the main use of FISH is screening for chromosomal aberrations. This presentation reviews the recent progress in preimplantation genetic diagnosis with an emphasis on bovine embryos. In particular the importance of biopsy size and strategies to avoid PCR contamination are discussed. Alternative DNA amplification and detection methods as well as methods to meet the challenge of multiple locus detection for marker assisted selection are presented.

Spontaneous and thiazolidinedione-induced B6C3F1 mouse hemangiosarcomas exhibit low ras oncogene mutation frequencies

Hemangiosarcomasare uncommon malignant endothelial cell tumors in humans and experimental animal species. The mechanisms giving rise to these tumors are poorly understood even though the histotypes are comparable between humans and rodents. Activating mutations in cellular ras protooncogenes have been detected in sporadic and chemically induced human and rodent hemangiosarcomas. Ras activation significantly modulates tumor angiogenesis, suggesting that mutations in ras genes might be causally related to vascular tumorigenesis. To more clearly define the role of ras in experimental vascular tumorigenesis, mutations in the Ki- and Ha-ras genes were characterized in 63 hemangiosarcomas that arose unexpectedly in control and treated B6C3F1 mice during a two-year carcinogenicity study of the thiazolidinedione troglitazone. DNA was extracted from paraffin sections of mouse hemangiosarcomas, control liver, or positive control hepatocellular carcinomas with defined mutations in the Ki- or Ha-ras genes. Exons 1 and 2 of the Ki- and Ha-ras genes were independently amplified using primer extension preamplification/locus-specific heminested PCR, and PCR amplicons were directly sequenced to identify mutations in codons 12, 13, or 61. Activating mutations were detected in 3 of 63 hemangiosarcomas: a single G-->A transition in the second position of Ki-ras codon 13 in a tumor from a treated animal and two G-->T transversions in the second position of Ha-ras codon 13, one in a single tumor from a control animal and one in a tumor from a treated animal. These mutations are consistent with endogenous mutagenesis arising from oxidative DNA damage. The low frequency of mutation (<5%) indicates that ras mutations did not contribute significantly to hemangiosarcoma development and suggests that mutational ras activation may not be a necessary step in vascular tumorigenesis in mice.

Freezer anthropology: new uses for old blood

Archived blood fractions (plasma, settled red cells, white cells) have proved to be a rich and valuable source of DNA for human genetic studies. Large numbers of such samples were collected between 1960 and the present for protein and blood group studies, many of which are languishing in freezers or have already been discarded. More are discarded each year because the usefulness of these samples is not widely understood. Data from DNA derived from 10-35-year-old blood samples have been used to address the peopling of the New World and of the Pacific. Mitochondrial DNA haplotypes from studies using this source DNA support a single wave of migration into the New World (or a single source population for the New World), and that Mongolia was the likely source of the founding population. Data from Melanesia have shown that Polynesians are recent immigrants into the Pacific and did not arise from Melanesia.

Genotypic analysis of flow-sorted and microdissected head and neck squamous lesions by whole-genome amplification

To investigate the utility of primer extension preamplification (PEP) in the genetic analysis of head and neck squamous tumorigenesis, microsatellite analysis was performed on matched deoxyribonucleic acid (DNA) samples extracted from 32 flow-sorted and microdissected specimens before and after PEP. Eighteen fresh and nine archival specimens were taken from invasive carcinomas, and five specimens were obtained from microdissected archival premalignant squamous epithelial lesions. Identical microsatellite patterns were observed in 276 (87%) of the 319 paired PEP and non-PEP genotypes with sufficient DNA. Overall, 13 (4%) of the PEP and 28 (8.8%) of the non-PEP fresh tissue samples failed specific microsatellite amplification. All 14 PEP-archival specimens were successfully amplified. Sorted cells showed a higher incidence (42.8%) of loss of heterozygosity (LOH) in both PEP and non-PEP samples compared with their unsorted counterparts. The results of this study indicate that (a) PEP is a simple and reliable technique for enhancing the DNA yield from small specimens; (b) flow sorting, in certain cases, improves the interpretation of genetic results; and (c) PEP may be used to compensate for PCR failure of unamplified DNA specimens in these lesions.

Genetic analysis using genomic representations

Analysis of the genetic changes in human tumors is often problematical because of the presence of normal stroma and the limited availability of pure tumor DNA. However, large amounts of highly reproducible "representations" of tumor and normal genomes can be made by PCR from nanogram amounts of restriction endonuclease cleaved DNA that has been ligated to oligonucleotide adaptors. We show here that representations are useful for many types of genetic analyses, including measuring relative gene copy number, loss of heterozygosity, and comparative genomic hybridization. Representations may be prepared even from sorted nuclei from fixed and archived tumor biopsies.

Genomic analysis of single cells from human basal cell cancer using laser-assisted capture microscopy

In this study, we show that direct mutational analysis of genomic DNA can be performed on single somatic cells extracted from a frozen, immunohistochemically stained tissue section using laser-assisted capture microscopy. Eighty-nine single tumor cells were separately dissected from one case of human basal cell cancer (BCC) and p53 mutations were analyzed by direct semi-automated sequencing of PCR fragments. Amplification was obtained for at least one of the two analyzed exons from approximately 50% of the single tumor cells. Identical p53 mutations were found in widely spread areas of the tumor, suggesting a clonal proliferation originating from one cell. Interestingly, comparison between results of immunohistochemistry and genetic analysis of the single cells revealed the same p53 mutations irrespective of the p53 immunoreactivity. We propose that this approach has a great potential to allow investigation of genotypic differences in single cells and more specifically to resolve important and fundamental questions determining cancer heterogeneity.

Whole genome amplification using a degenerate oligonucleotide primer allows hundreds of genotypes to be performed on less than one nanogram of genomic DNA

Genetic analysis of limiting quantities of genomic DNA play an important role in DNA forensics, paleoarcheology, genetic disease diagnosis, genetic linkage analysis, and genetic diversity studies. We have tested the ability of degenerate oligonucleotide primed polymerase chain reaction (DOP-PCR) to amplify picogram quantities of human genomic DNA for the purpose of increasing the amount of template for genotyping with microsatellite repeat markers. DNA was uniformly amplified at a large number of typable loci throughout the human genome with starting template DNAs from as little as 15 pg to as much as 400 ng. A much greater-fold enrichment was seen for the smaller genomic DOP-PCRs. All markers tested were amplified from starting genomic DNAs in the range of 0.6-40 ng with amplifications of 200- to 600-fold. The DOP-PCR-amplified genomic DNA was an excellent and reliable template for genotyping with microsatellites, which give distinct bands with no increase in stutter artifact on di-, tri-, and tetranucleotide repeats. There appears to be equal amplification of genomic DNA from 55 of 55 tested discrete microsatellites implying near complete coverage of the human genome. Thus, DOP-PCR appears to allow unbiased, hundreds-fold whole genome amplification of human genomic DNA for genotypic analysis.

Use of tagged random hexamer amplification (TRHA) to clone and sequence minute quantities of DNA--application to a 180 kb plasmid isolated from Sphingomonas F199

We have developed a novel method to clone and sequence minute quantities of DNA. The method was applied to sequence a 180 kb plasmid pNL1. The first step was the production of a size distributed population of DNA molecules that were derived from the 180 kb plasmid pNL1. The first step was accomplished by a random synthesis reaction using Klenow fragment and random hexamers tagged with a T7 primer at the primer 5'-end (T7-dN6, 5'-GTAATACGACTCACTATAGGGCNNNNNN-3'. In the second step, Klenow-synthesized molecules were amplified by PCR using T7 primer (5'-GTAATACGACTCACTATAGGGC-3'). With a hundred nanograms starting plasmid DNA from pNL1, we were able to generate Klenow-synthesized molecules with sizes ranging from 28 bp to >23 kb which were detectable on an agarose gel. The Klenow-synthesized molecules were then used as templates for standard PCR with T7 primer. PCR products of sizes ranging from 0.3 to 1.3 kb were obtained for cloning and sequencing. From the same Klenow-synthesized molecules, we were also able to generate PCR products with sizes up to 23 kb by long range PCR. A total 232.5 kb sequences were obtained from 593 plasmid clones and over twenty putative genes were identified. Sequences from these 593 clones were assembled into 62 contigs and 99 individual sequence fragments with a total unique sequence of 86.3 kb.