High-throughput automated DNA sequencing facility with fluorescent labels at the European Molecular Biology Laboratory

Sparsely cross-linked "nanogel" matrixes as fluid, mechanically stabilized polymer networks for high-throughput microchannel DNA sequencing

We have developed sparsely cross-linked "nanogels", subcolloidal polymer structures composed of covalently linked, linear polyacrylamide chains, as novel replaceable DNA sequencing matrixes for capillary and microchip electrophoresis. Nanogels were synthesized via inverse emulsion (water-in-oil) copolymerization of acrylamide and a low percentage (approximately 10(-4) mol %) of N,N-methylene bisacrylamide (Bis). Nanogels and nanogel networks were characterized by multiangle laser light scattering and rheometry, respectively, and tested for DNA sequencing in both capillaries and chips with four-color LIF detection. Typical nanogels have an average radius of approximately 230 nm, with approximately 75% of chains incorporating a Bis cross-linker. The properties and performance of nanogel matrixes are compared here to those of a linear polyacrylamide (LPA) network, matched for both polymer weight-average molar mass (M(w)) and the extent of interchain entanglements (c/c). At sequencing concentrations, the two matrixes have similar flow characteristics, important for capillary and microchip loading. However, because of the physical network stability provided by the internally cross-linked structure of the nanogels, substantially longer average read lengths are obtained under standard conditions with the nanogel matrix at a 98.5% accuracy of base-calling (for CE: 680 bases, an 18.7% improvement over LPA, with the best reads as long as 726 bases, compared to 568 bases for the LPA matrix). We further investigated the use of the nanogel matrixes in a high-throughput microfabricated DNA sequencing device consists of 96 separation channels densely fabricated on a 6-in. glass wafer. Again, preliminary DNA sequencing results show that the nanogel matrixes are capable of delivering significantly longer average read length, compared to an LPA matrix of comparable properties. Moreover, nanogel matrixes require 30% less polymer per unit volume than LPA. The addition of a small amount of low molar mass LPA or ultrahigh molar mass LPA to the optimized nanogel sequencing matrix further improves read length as well as the reproducibility of read length (RSD < 1.6%). This is the first report of a replaceable DNA sequencing matrix that provides better performance than LPA, in a side-by-side comparison of polymer matrixes appropriately matched for molar mass and the extent of interchain entanglements. These results could have significant implications for the improvement of microchip-based DNA sequencing technology.

DNA sequencing and genotyping in miniaturized electrophoresis systems

Advances in microchannel electrophoretic separation systems for DNA analyses have had important impacts on biological and biomedical sciences, as exemplified by the successes of the Human Genome Project (HGP). As we enter a new era in genomic science, further technological innovations promise to provide other far-reaching benefits, many of which will require continual increases in sequencing and genotyping efficiency and throughput, as well as major decreases in the cost per analysis. Since the high-resolution size- and/or conformation-based electrophoretic separation of DNA is the most critical step in many genetic analyses, continual advances in the development of materials and methods for microchannel electrophoretic separations will be needed to meet the massive demand for high-quality, low-cost genomic data. In particular, the development (and commercialization) of miniaturized genotyping platforms is needed to support and enable the future breakthroughs of biomedical science. In this review, we briefly discuss the major sequencing and genotyping techniques in which high-throughput and high-resolution electrophoretic separations of DNA play a significant role. We review recent advances in the development of technology for capillary electrophoresis (CE), including capillary array electrophoresis (CAE) systems. Most of these CE/CAE innovations are equally applicable to implementation on microfabricated electrophoresis chips. Major effort is devoted to discussing various key elements needed for the development of integrated and practical microfluidic sequencing and genotyping platforms, including chip substrate selection, microchannel design and fabrication, microchannel surface modification, sample preparation, analyte detection, DNA sieving matrices, and device integration. Finally, we identify some of the remaining challenges, and some of the possible routes to further advances in high-throughput DNA sequencing and genotyping technologies.

Random exploration of the Kluyveromyces lactis genome and comparison with that of Saccharomyces cerevisiae

The genome of the yeast Kluyveromyces lactis was explored by sequencing 588 short tags from two random genomic libraries (random sequenced tags, or RSTs), representing altogether 1.3% of the K. lactis genome. After systematic translation of the RSTs in all six possible frames and comparison with the complete set of proteins predicted from the Saccharomyces cerevisiae genomic sequence using an internally standardized threshold, 296 K.lactis genes were identified of which 292 are new. This corresponds to approximately 5% of the estimated genes of this organism and triples the total number of identified genes in this species. Of the novel K.lactis genes, 169 (58%) are homologous to S.cerevisiae genes of known or assigned functions, allowing tentative functional assignment, but 59 others (20%) correspond to S.cerevisiae genes of unknown function and previously without homolog among all completely sequenced genomes. Interestingly, a lower degree of sequence conservation is observed in this latter class. In nearly all instances in which the novel K.lactis genes have homologs in different species, sequence conservation is higher with their S.cerevisiae counterparts than with any of the other organisms examined. Conserved gene order relationships (synteny) between the two yeast species are also observed for half of the cases studied.

Optimization and automation of fluorescence-based DNA hybridization for high-throughput clone mapping

Large-scale hybridization-based genome mapping projects, such as the production of sequence-ready physical clone maps, call for robust and cheap DNA labeling techniques that are amenable to automation. We routinely use a high-throughput protocol based on fluorescence detection. DNA probes are labeled via polymerase chain reaction (PCR) amplification with primers that are digoxigenin-modified at their 5' ends. Alternatively, digoxigenin-labeled dUTP is incorporated in a random hexamer priming reaction. Hybridization takes place in small volumes by sandwiching the probe between filters and plastic sheets. A fluorescence signal is produced by the activity of alkaline phosphatase attached to an anti-digoxigenin antibody upon the addition of AttoPhos substrate. Signals are directly detected with a charge-coupled device (CCD) camera and scored by an image data analysis system. DNA filters can be reused at least 40 times without loss of data quality. Significant advantages compared to radioactive techniques are the reduced health risk, enabling highly parallel processing; the production of spot signals uniform in size and intensity, which is essential for efficient image analysis; and a cost reduction of about 70%.

DNA sequencing and analysis of 130 kb from yeast chromosome XV

We have determined the nucleotide sequence of 129,524 bases of yeast (Saccharomyces cerevisiae) chromosome XV. Sequence analysis revealed the presence of 59 non-overlapping open reading frames (ORFs) of length > 300 bp, three tRNA genes, four delta elements and one Ty-element. Among the 21 previously known yeast genes (36% of all ORFs in this fragment) were nucleoporin (NUP1), ras protein (RAS1), RNA polymerase III (RPC1) and elongation factor 2 (EF2). Further, 31 ORFs (53% of the total) were found to be homologous to known protein or DNA sequences, or sequence patterns. For seven ORFs (11% of the total) no homology was found. Among the most interesting protein identification in this DNA fragment are an inositol polyphosphatase, the second gene of this type found in yeast (homologous to the human OCRL gene involved in Lowe's syndrome), a new ADP ribosylation factor of the arf6 subfamily, the first protein containing three C2 domains, and an ORF similar to a Bacillus subtilis cell-cycle related protein. For each ORF detailed sequence analysis was carried out, with a full consideration of its biological function and pointing out key regions of interest for further functional analysis.

Evaluation of the ALF DNA sequencer for high-speed sizing of short tandem repeat alleles

DNA typing of short tandem repeat (STR) loci with automated real-time analysis of the fluorescent polymerase chain reaction (PCR) products has given forensic DNA analysis a new dimension. In the present work, the ALF DNA sequencer was evaluated for automated size determination of tetra-nucleotide STRs at high speed. Short gel plates, with a well-to-laser distance of 10 cm, allowed for the analysis of four STR loci (HUMvWF, HUMTH01, D21S11 and HPRT) in one gel lane in less than 75 min. Allele size determination was done with two external allelic ladders for each locus. Lane-to-lane variations were overcome by the inclusion in each lane of two fluorescent PCR products of constant size (123 and 375 bp) that migrated below and above the multiplex of the four STR loci. The accuracy of sizing and allele detection within and between different gels was high (99.89%) for all four STR systems investigated and the gels could be reloaded without a decrease in accuracy of the allele size estimation. This way, the throughput of the system was increased, which is of interest for linkage studies, gene mapping, and population diversity studies.

Multiple affinity domains for the detection, purification and immobilization of recombinant proteins

Affinity systems based on specific molecular recognition are valuable tools for detection, purification and immobilization of recombinant proteins. Here, novel multipartite affinity fusion vectors were assembled and investigated to allow flexible binding and elution conditions. The rationale for the assembly of different combinations of affinity domains was to take advantage of the wide variety of molecular interactions of these domains for purification, solubilization, detection and immobilization. In total, seven different affinity tags representing five different types of tag-ligand interactions were studied: (i) monoclonal antibodies-peptides (T7-tag and FLAG peptide); (ii) streptavidin-peptide (Strep-tag); (iii) hexahistidyl-metal ions (His6-tag; (iv) bacterial receptors-serum proteins (staphyloccal protein A-Fc and streptococcal protein G-serum albumin); (v) streptavidin-biotin (in vivo biotinylated peptide). Selected tags were evaluated for the production and purification of Escherichia coli DNA polymerase I (Klenow fragment). On the basis of the results, a vector (pAff2c) was assembled using a novel combination of affinity domains: (i) an in vivo biotinylated peptide; (ii) a His6 sequence, and (iii) a highly soluble serum albumin binding region. Using these three affinities, a wide variety of conditions can be employed for both the binding and the elution steps.

Sequencing and analysis of 51 kb on the right arm of chromosome XV from Saccharomyces cerevisiae reveals 30 open reading frames

We have sequenced a region of 51 kb of the right arm from chromosome XV of Saccharomyces cerevisiae. The sequence contains 30 open reading frames (ORFs) of more than 100 amino acid residues. Thirteen new genes have been identified. Thirteen ORFs correspond to known yeast genes. One delta element and one tRNA gene were identified. Upstream of the RPO31 gene, encoding the largest subunit of RNA polymerase III, lies a Abf1p binding site. The nucleotide sequence data reported in this paper are available in the EMBL, GenBank and DDBJ nucleotide sequence databases under the Accession Number X90518.

Sequencing 39,350 bp of Saccharomyces cerevisiae chromosome XII utilizing ordered shotgun libraries

A 39,350 bp cosmid containing DNA of Saccharomyces cerevisiae chromosome XII was sequenced by making use of ordered sub-clones of 1 kb insert-length selected from a physical clone map. In a first analysis, 96 clones were sequenced from both ends (10 gels) with two standard sequencing primers covering 91% of the total sequence (49% double-stranded). After selection of another eight clones six gaps of a total of 1.8 kb and several single-stranded stretches remained. These gaps were closed by 86 primer walks leading to an overall redundancy of 4.4 per base and a total of 292 sequencing reactions. The number of walking primers can be reduced significantly by more uniform clone lengths and longer sequencing reads, thus, the total amount of sequencing reactions can approach the minimum value achieved with primer walking strategies, with only very few walking primers needed for gap closure.

A modified screening method for pcDNA-1 expression libraries which is applicable to both surface and intracellular antigens. Cloning of a colon carcinoma antigen

A modification of the COS cell transient expression system for isolating cDNAs encoding antibody defined molecules was used to clone the carcinoma antigen defined by monoclonal antibody CORC-6. Bacteria were transformed with cDNA prepared from colon carcinoma cells and immediately divided into aliquots containing approximately 500 independent clones. The aliquots were grown for 8 h before frozen. Screening was performed by isolating plasmid from each aliquot, transfecting COS-7 cells, and staining the fixed cells with monoclonal antibodies. The aliquot giving rise to antigen expressing cells was then subdivided until a single antigen encoding colony was obtained. This approach allows a library to be rescreened and is equally applicable to cytoplasmic and surface localized antigens. Using this approach, the antigen defined by CORC-6 was identified as carcinoembryonic antigen.

Fluorescence-based classification of microsatellites using a single-wavelength semiautomatic sequencer: genotype assignment and identity tests by analysis of comigrating peak profiles

Microsatellite analysis, based on fluorescein labeling and reading through a semiautomatic single wavelength sequencer, is described. Pairs of labeled polymerase chain reaction (PCR) samples, mixed in equimolar proportion, were electrophoresed and the specific peaks read in a single gel lane. Identity was asserted when peaks overlapped in a unique fluorescent signal which, compared with individual sample profiles, had a twofold intensity. Classification was achieved by blending individual PCR products to 'locus specific allelic ladders' (composite samples containing a repertory of fragments allelic to a given locus) and by noticing the specific peak enhancement. The resulting protocol of analysis assigned no size and classified allelic forms by tandem repeat number. Applied to a large repertory of PCR products and compared with manual electrophoresis, this protocol proved to be reliable and reduced times and costs of genotype analysis. Analysis of comigrating peak profiles is highly objective and provides convincing evidence for diagnostics and identity tests.

Efficient low redundancy large-scale DNA sequencing at EMBL

An efficient low redundancy DNA sequencing strategy should allow high accuracy determination of the consensus sequence on both strands of a DNA fragment from a minimal number of sequencing reactions with minimal overlap. At EMBL we developed a directed strategy for cosmid-scale sequencing based on primer walking, whereas most other sequencing projects of this scale rely on the random 'shotgun' strategy. In our strategy, highly accurate raw data are obtained from automated double-stranded Sanger dideoxy sequencing with inexpensive walking primers (8 to 10 $ per primer), T7 DNA polymerase and internal labelling by fluorescein-15- dATP on A.L.F. DNA sequencers (Pharmacia Biotech). The use of 60-cm long glass plates enables reading length of up to 1000 bases. Comparing various random and directed sequencing strategies in the course of the European Community yeast genome sequencing project on cosmids from chromosomes IX, XI and XV, primer walking was found to be the strategy resulting in the lowest possible redundancy of 2.6 to 2.8. Future development of the sequencing strategy is based on the new EMBL 2-dye sequencing device for simultaneous sequencing on both strands, and implementation of an initial limited random sequencing phase to reduce the number of walking primers required by a factor of 3, while still maintaining a low redundancy of approx. 3.

High-resolution physical mapping of a 250-kb region of human chromosome 11q24 by genomic sequence sampling (GSS)

A physical map of the region of human chromosome 11q24 containing the FLI1 gene, disrupted by the t(11;22) translocation in Ewing sarcoma and primitive neuroectodermal tumors, was analyzed by genomic sequence sampling. Using a 4- to 5-fold coverage chromosome 11-specific library, 22 region-specific cosmid clones were identified by phenol emulsion reassociation hybridization, with a 245-kb yeast artificial chromosome clone containing the FLI1 gene, and by directed "walking" techniques. Cosmid contigs were constructed by individual clone fingerprinting using restriction enzyme digestion and assembly with the Genome Reconstruction and AsseMbly (GRAM) computer algorithm. The relative orientation and spacing of cosmid contigs with respect to the chromosome was determined by the structural analysis of cosmid clones and by direct visual in situ hybridization mapping. Each cosmid clone in the contig was subjected to "one-pass" end sequencing, and the resulting ordered sequence fragments represent approximately 5% of the complete DNA sequence, making the entire region accessible by PCR amplification. The sequence samples were analyzed for putative exons, repetitive DNAs, and simple sequence repeats using a variety of computer algorithms. Based upon the computer predictions, Southern and Northern blot experiments led to the independent identification and localization of the FLI1 gene as well as a previously unknown gene located in this region of chromosome 11q24. This approach to high-resolution physical analysis of human chromosomes allows the assembly of detailed sequence-based maps and provides a tool for further structural and functional analysis of the genome.

Nucleotide sequence and analysis of the centromeric region of yeast chromosome IX

We have determined the nucleotide sequence of a cosmid (pIX338) containing the centromere region of yeast (Saccharomyces cerevisiae) chromosome IX. The complete nucleotide sequence of 33.8 kb was obtained by using an efficient directed sequencing strategy in combination with automated DNA sequencing on the A.L.F. DNA sequencer. Sequence analysis revealed the presence of 17 open reading frames (ORFs), four of them previously known yeast genes (sly12, pan1, sts1 and prl1), a tRNA gene and the centromere motif. Exhaustive database searches detected sequence homologues of known function for as many as 14 of the 17 ORFs. These include a mammalian tyrosine kinase substrate; the Escherichia coli cell cycle protein MinD; the human inositol polyphosphate-5-phosphatase (gene OCRL) involved in Lowe's syndrome, a developmental disorder; and helicases, for which the new yeast member defines a distinct DEAD/H-box subfamily. A surprisingly large fraction of the ORFs (at least six out of 17) in the centromeric region are apparently involved in RNA or DNA binding.

DMSO resolves certain compressions and signal dropouts in fluorescent dye labeled primer-based DNA sequencing reactions

Automated base calling algorithms are more sensitive to the quality of the DNA sequencing data than are the labor intensive visual methods of base calling. To improve this quality, data from DNA sequencing reactions have been compared in order to determine the effects of the inclusion of dimethyl sulfoxide (DMSO). Inclusion of 10% DMSO into the reaction cocktail resolves at least one type of sequence compression. This compression may be due to the lack of ability in T7 DNA polymerase to read through certain sequences correctly. The poor quality of these data is seen as radioactive bands or fluorescent signal peaks that have an abnormal alignment, either in the wrong order or as single bands/peaks. The inclusion of DMSO also resolves sequences where the peak signal is absent or severely diminished, leading to a "gap" in the chromatogram profile. DMSO is better than deaza-dITP for resolving certain compressions. Addition of DMSO is a cheaper and more efficient method for high-throughput DNA sequencing than repeating reactions with base analogs.

Human protein NEFA, a novel DNA binding/EF-hand/leucine zipper protein. Molecular cloning and sequence analysis of the cDNA, isolation and characterization of the protein

The cDNA libraries constructed from the human acute lymphoblastic leukemia cell line KM3 in the expression vector lambda gt11, were screened with the anti-CALLA (common acute lymphoblastic leukemia antigen) mAb (monoclonal antibody) J5. The selected J5-positive clone I containing a partial cDNA insert was isolated and sequenced. For completing the cDNA sequence the cDNA libraries were further screened by hybridization with the DIG (digoxigenin)-labelled DNA probe derived from clone I, the 5'-end region was analysed by 5'-RACE (rapid amplification of cDNA ends) using a sequence specific primer. In total a 1639 bp cDNA sequence was determined. The cDNA sequence contains a 1260 bp open reading frame and the untranslated 3'- and 5'-end sides. The 420 residue amino acid sequence, deduced from the cDNA sequence, unexpectedly differs fundamentally from CALLA (CD10) although clones I and II were J5-positive in immuno screening. The mature protein corresponding to the cDNA was isolated and characterized from the KM3 cells using polyclonal antisera raised against the in vitro expressed polypeptide from clone I. The protein is expressed on plasma membrane, in cytosol and is secreted into culture medium, its relative molecular mass was determined to be 55 kDa on SDS-PAGE. The deduced amino acid sequence from cDNA was confirmed by peptide sequences. The new protein contains a basic amino acid rich putative DNA binding domain (b) with a potential nuclear targeting signal, two helix-loop-helix (HLH) motif regions, concurrently EF-hand motifs, an acidic amino acid rich region (a) between the EF-hands, and a leucine zipper (Z) motif. This DNA binding protein therefore is characterized by a linked motif "b/HLH/a/HLH/Z". The protein was designated NEFA: DNA binding/EF-hand/acidic amino acid rich region.

Genomic sequence sampling: a strategy for high resolution sequence-based physical mapping of complex genomes

We present a simple and efficient method for constructing high resolution physical maps of large regions of genomic DNA based upon sampled sequencing. The physical map is constructed by ordering high density cosmid contigs and determining a sequence fragment from each end of every clone. The resulting map, which contains 30-50% of the complete DNA sequence, allows the identification of many genes and makes possible PCR amplification of virtually any part of the genome. We apply this strategy to the automated analysis of the genome of the primitive eukaryote Giardia lamblia and evaluate its applicability to the physical mapping and DNA sequencing of the human genome.

Purification and characterization of two putative HLA class II associated proteins: PHAPI and PHAPII

In addition to their well defined role in presentation of processed antigen on the cell surface, class II molecules are able to transduce signals into the cell after binding of ligands. The cytoplasmic regions of class II molecules might function as docking sites for as yet unidentified proteins that are components of this signalling pathway. Here we report on two putative HLA class II associated proteins (PHAPI and PHAPII) which have been purified from the cytosolic fraction of the human lymphoblastoid B-cell line H2LCL using an affinity matrix composed of the synthetic biotinylated cytoplasmic region of the DR2 alpha chain immobilized on avidin agarose. The sequence obtained for PHAPI revealed a novel primary structure with a leucine/isoleucine rich N-terminal region. Protein data and the cDNA sequence obtained for PHAPII agree with the cDNA sequence of SET that has been described recently. Both PHAPI and PHAPII have an extended highly acidic C-terminal region. Based on their primary structure we speculate that PHAPI and PHAPII are involved in the generation of intracellular signalling events that lead to regulation of transcriptional activity after binding of a ligand to HLA class II molecules.

How far have we progressed toward automated electrophoresis in sieving media of the twenty-first century?

The automation of electrophoresis in polymeric sieving media requires (i) an objective definition of the conditions (the polymer, its concentration, solvent, buffer, pH, ionic strength, temperature) under which a particular separation proceeds most effectively; (ii) apparatus capable of zone detection, acquisition by computer, evaluation (migration distance, zone width and area) and a print-out of the number of components, their size and net charge, and the polymer conditions under which each component separates most effectively from its two neighboring zones. Both of these prerequisites of automation have been met to a first approximation at this time and, after further maturation, assembly and streamlining should be able to fill the need of the coming century for a more efficient, nonarbitrary and cost-effective mode of macromolecular and cellular particle separation. (i) The realization of the qualitative equivalence of polymer solutions and gels has greatly increased our options in the choice of sieving media. That choice can be made objectively by correlating separation efficiency (S), particle size (R) and intrinsic viscosity (eta o) of the polymer. (S) is a function of the slope, KR(R), of the Ferguson plot [log(mobility) vs. polymer concentration], or with nonlinear plots (DNA, agarose) KR(R,T). KR is at present most easily derived from transverse pore gradient gels or by conducting capillary zone electrophoresis (CZE) at multiple polymer concentrations. Pore gradient CZE appears promising. CZE also defines the free mobility unequivocally. Computer programs exist to generate KR from migration distances (times), and optimal S and polymer concentration for a particular R from KR.(ABSTRACT TRUNCATED AT 250 WORDS)

Sequencing and analysis of 51.6 kilobases on the left arm of chromosome XI from Saccharomyces cerevisiae reveals 23 open reading frames including the FAS1 gene

We have sequenced two segments containing a total of 51.6 kb of the left arm from chromosome XI of Saccharomyces cerevisiae. The first segment of 38.5 kb contains 18 open reading frames (ORFs) of more than 100 amino acid residues. Five ORFs encode known yeast genes, including the fatty acid synthase gene (FAS1). Three new yeast genes were discovered with homologies to non-yeast genes and ten new genes without homologies to any known sequences. The second segment of 13 kb contains five ORFs with two known yeast genes and three unknown genes. The sequences from cosmid pUKG041 were obtained entirely with the walking primer strategy resulting in a very low overall sequence redundancy of 2.8 and an average reading length of 443 bases.

Identification of predominant T-cell receptor rearrangements by temperature-gradient gel electrophoresis and automated DNA sequencing

To assess the diversity of T-cell receptor (TCR) gene rearrangements in uncloned lymphocytes we used a three-stage strategy that allows the detection of a restricted TCR repertoire and the identification of the predominant, rearranged sequence(s). We have analyzed in parallel T cells obtained from a renal cell carcinoma infiltrate that specifically lyse the autologous tumor after in vitro culture and T cells from autologous peripheral blood. First, DNA amplification by the polymerase chain reaction (PCR) was performed with a number of oligonucleotide primers specific for several TCR V alpha gene families. All V alpha primers displayed specific amplification products in the peripheral blood, while a restricted TCR repertoire was present in the tumor-infiltrating lymphocytes. Subsequently, positively amplified PCR products were run in a temperature-gradient gel electrophoresis. A limited number of bands corresponding to predominant homo- and heteroduplexes was only found in the tumor-infiltrating lymphocytes. The presence of a low number of rearranged TCR sequences in these samples was confirmed by automated single-stranded DNA sequencing using a single fluorescent dye. These results support the broad application of this strategy for targeted sequencing of those PCR products carrying predominant DNA templates without previous DNA cloning.

The Caenorhabditis elegans genome sequencing project: first steps in automation

Novel biochemical processes using magnetic particles have been automated to provide a robotic system to perform processes for large-scale shotgun sequencing.