Nucleic acid compositions conferring altered metabolic characteristics

FIELD OF THE INVENTION

This invention relates to deoxyribonucleic acid (DNA) and amino acid sequences that confer altered metabolic characteristics in plants.

BACKGROUND OF THE INVENTION

Plants are photosynthetic organisms able to fix inorganic carbon (CO.sub.2) in organic matter via energy from light and minerals contained in water. All carbons fixed primarily via the pentose/triose phosphate cycle are converted in numerous anabolic pathways necessary to sustain life (primary metabolism). To survive plants must adapt to their environment and synthesize an extremely wide range of organic compounds required to interact with the elements of their microenvironment (secondary metabolism). To capture the biochemical diversity of this particular kingdom both primary and secondary metabolism have to be taken into account. The primary metabolism is represented by the biosynthesis of building blocks of macromolecules such as amino acids, fatty acids, carbohydrates, and sterols.

Each of these groups of compounds is of economic importance. Fatty acids can also be used as a raw material for industrial applications in a variety of products, including soaps, lubricants, paints, detergents, adhesives, and plasticizers. Furthermore, fatty acids are the major components of edible oils. For example, fatty acid compounds are involved in building blocks for protection (cell membrane, epicuticular polymers), storage of energy in the plant seeds and as secondary messengers in the plant cell. As another example, carbohydrates are intermediates in the biosynthesis of energy reserves (starch, cellulose) and building blocks of the cell wall giving the plant shape and structure. The carbohydrates are the carbon skeletons of many biosynthetic reactions. As such, the ability to alter carbohydrate metabolism could lead to many improvements in plants, including increased transport and accumulation of starch by accumulation of hexose phosphate that could improve starch yield in the seed and the plant; alterations in the cell wall for better resistance to pest and drought; better digestibility for forage plants; and better processivity for pulp production in paper industry (e.g., less lignin and hemicellulose).

The advent of modern biology, particularly molecular biology and genetics, has opened up new avenues for altering the production of compounds of economic importance by plants. Scientists have focused on utilizing recombinant DNA (rDNA) methods, that allow new varieties of plants to be produced much faster than by conventional breeding. rDNA techniques allow the introduction of genes from distantly related species or even from different biological kingdoms into crop plants, conferring traits that provide significant agronomic advantages. Furthermore, detailed knowledge of the traits being introduced, such as cellular function and localization, can lead to less variability in offspring, and fine-tuning of secondary effects (e.g., permitting variation from what is customarily observed). After a trait has been introduced into a plant by transgenic methods, conventional breeding can be used to hybridize the transgenic line with useful varieties and elite germplasms, resulting in crops containing numerous advantageous properties.

Most efforts to engineer plants with specific traits thus far have been based on the rational design paradigm of transforming a plant with a gene of known function with the intent of introducing a known trait. As agricultural biotechnology hurtles into the genomics and post-genomics era, the massive amounts of genetic and functional data being generated are being used to direct the search for genes that can be utilized with recombinant methods. However, if the use of this information is limited to the rational design paradigm, the identification of genes with truly profound effects on the production of desired compounds by plants could be extremely time-consuming and slow.

Accordingly, what is needed in the art are methods for rapidly screening and identifying gene sequences and polypeptide sequences of previously unknown function whose expression causes altered metabolic characteristics in biological systems, including, but not limited to, plants.

SUMMARY OF THE INVENTION

This invention relates to deoxyribonucleic acid (DNA) and amino acid sequences that confer altered metabolic characteristics in plants. In some embodiments, the present invention provides polynucleotides and polypeptides that confer altered metabolic characteristics when expressed in plants. The present invention is not limited to the alteration of amounts or levels of any particular metabolite. Indeed, the alteration of the levels or amounts of a variety of metabolites is contemplated, including, but not limited to acids, fatty acids, amino acids, hydroxy fatty acids, branched fatty acids, carbohydrates, hydrocarbons, glycerides, phenols, strerols, oxygenated terpenes, and other isoprenoids, alcohols, alkenes and alkynes. The present invention is not limited to any particular polypeptide or polynucleotide sequences that confer altered metabolic characteristics. Indeed, a variety of such sequences are contemplated. Accordingly, in some embodiments the present invention provides an isolated nucleic acid selected from the group consisting of SEQ ID NOs: 1-7554 and nucleic acid sequences that hybridize to any thereof under conditions of low stringency, wherein expression of the isolated nucleic acid in a plant results in a altered metabolic characteristic.

In some embodiments, the present invention provides an isolated nucleic acid selected from SEQ ID NOs: 162, 212, 3781, 3970, 3990, 492, 3796, 3975, and 4028, wherein expression of the nucleic acid in a plant results in altered acid metabolism. In other embodiments, the present invention provides an isolated nucleic acid selected from SEQ ID NOs: 4049, 210, 4045, 229, 3825, 4015, 3835, 4039, 1048 and 1106, wherein expression the nucleic acid in a plant results in altered alcohol metabolism. In still other embodiments, the present invention provides an isolated nucleic acid selected from SEQ ID NOs: 7548, 283, 3957, 3734, 3739, 3797, 7516, 3762, 4020 and 1062, wherein expression of the nucleic acid in a plant results in altered fatty acid metabolism. In further embodiments, the present invention provides an isolated nucleic acid of selected from SEQ ID NOs: 1148, 4147, 273, 281, 299, 3920, 450, 7463 and 4074, wherein expression of the nucleic acid in a plant results in altered branched fatty acid metabolism.

In still further embodiments, the present invention provides an isolated nucleic acid selected from SEQ ID NOs: 258, 456, 3859, 3817, 4018, 3848, 3862, 4008 and 1000, wherein expression of the nucleic acid in a plant results in altered alkaloid or other base metabolism. In some embodiments, the present invention provides an isolated nucleic acid selected from SEQ ID NOs: 372, 3714, 3717, 3963, 3775, 3757, 7462, 3743, 3744 and 7480, wherein expression of the nucleic acid in a plant results in altered amino acid metabolism.

In some other embodiments, the present invention provides an isolated nucleic acid selected from SEQ ID NOs: 7404, 180, 181, 225, 231, 366, 3983, 3833, 1121 and 1062, wherein expression of the nucleic acid in a plant results in altered ester metabolism. In some further embodiments, the present invention provides an isolated nucleic acid selected from SEQ ID NOs: 3773, 583, 3821, 7403, 988, 1002, 1007 and 1129, wherein expression of the nucleic acid in a plant results in altered glyceride metabolism. In still other embodiments, the present invention provides an isolated nucleic acid selected from SEQ ID NOs: 150, 7410, 175, 7553, 619, 1078, 1122 and 1124, wherein expression of the nucleic acid in a plant results in altered phenolic compound metabolism.

In further embodiments, the present invention provides an isolated nucleic acid selected from SEQ ID NOs: 3891, 7545, 7551, 4121, 157, 159, 7411, 3792, 3799 and 3997, wherein expression of the nucleic acid in a plant results in altered carbohydrate metabolism. In other embodiments, the present invention provides an isolated nucleic acid selected from SEQ ID NOs: 7405, 7406, 173, 183, 220, 227, 3778, 3803, 3847 and 1005, wherein expression of the nucleic acid in a plant results in altered sterol, oxygenated terpene, or isoprenoid metabolism. In still other embodiments, the present invention provides an isolated nucleic acid selected from SEQ ID NOs: 7408, 351, 378, 3864, 4103, 996, 1006 and 1098, wherein expression of the nucleic acid in a plant results in altered alkene or alkyne metabolism.

In further embodiments, the present invention provides an isolated nucleic acid selected from SEQ ID NOs: 177, 7442, 4038, 3836, 3855, 1012, 1015, 1119 and 1024, wherein expression of the nucleic acid in a plant results in altered hydrocarbon metabolism. In still further embodiments, the present invention provides an isolated nucleic acid selected from SEQ ID NOs: 360, 4001, 3703, 7399, 645, 3849 and 7552, wherein expression of the nucleic acid in a plant results in altered ketone or quinone metabolism.

In further preferred embodiments, the present invention provides vectors comprising the foregoing polynucleotide sequences. In still further embodiments, the foregoing sequences are operably linked to an exogenous promoter, most preferably a plant promoter. However, the present invention is not limited to the use of any particular promoter. Indeed, the use of a variety of promoters is contemplated, including, but not limited to, 35S, 19S, heat shock, and Rubisco promoters, subgenomic promoters such as the CaMV promoter and TMV coat protein promoter, and dual promoters systems such as DHSPES (see U.S. Pat. No. 6,303,848, incorporated herein by reference). In some embodiments, the nucleic acid sequences of the present invention are arranged in sense orientation, while in other embodiments, the nucleic acid sequences are arranged in the vector in antisense orientation. In still further embodiments, the present invention provides a plant comprising one of the foregoing nucleic acid sequences or vectors, as well as seeds, leaves, roots, stems and fruit from the plant. In some particularly preferred embodiments, the present invention provides at least one of the foregoing sequences for use in conferring altered metabolism in a plant.

In still other embodiments, the present invention provides processes for making a transgenic plant comprising providing a vector as described above and a plant, and transfecting the plant with the vector. In other preferred embodiments, the present invention provides processes for providing an altered metabolic characteristic in a plant or population of plants comprising providing a vector as described above and a plant, and transfecting the plant with the vector under conditions such that an altered metabolic characteristic is conferred by expression of the isolated nucleic acid from the vector. In still further embodiments, the present invention provides an isolated nucleic acid selected from the group consisting of SEQ ID NOs: 1-7554 and nucleic acid sequences that hybridize to any thereof under conditions of low stringency for use in producing a plant with altered metabolism. In other embodiments, the present invention provides an isolated nucleic acid, composition or vector substantially as described herein in any of the examples or claims.

DESCRIPTION OF THE DRAWINGS

FIG. 1 presents the contig sequences corresponding to SEQ ID NOs:1-1165, 3703-4153, and 7389-7554.

FIG. 2 presents homologous sequences 1166-3702 and 4154-7388.

FIG. 3 is a table of BLAST search results from public databases.

FIG. 4 is a table of BLAST search results from the Derwent.TM. amino acid database.

FIG. 5 is a table of BLAST search results from the Derwent.TM. nucleotide database.

FIG. 6 provides a summary of the metabolic alterations caused by expression of the indicated sequences. NQ=present in reference, but not detected or below the limit of quantification in the sample.

FIGS. 7a-d summarizes the gas chromatography flame ionization detection (GC/FD) parameters used to analyze metabolite samples.

FIG. 8 provides a list of SEQ IDs, grouped by the respective fractionation chemistry, observed to confer altered metabolic characteristic in plants based exclusively upon a pattern recognition automated data analysis technique (ADA).

FIG. 9 provides tables exemplifying the functional correlation between sequences that share homology.

DEFINITIONS

Before the present proteins (including their fragments and peptides), nucleotide sequences, and methods are described, it should be noted that this invention is not limited to the particular methodology, protocols, cell lines, vectors, and reagents described herein as these may vary. It should also be understood that the terminology used herein is for the purpose of describing particular aspects of the invention, and is not intended to limit its scope.

It must be noted that as used herein and in the appended claims, the singular forms "a", "an", and "the" include plural reference unless the context clearly dictates otherwise. Thus, for example, reference to "a host cell" includes a plurality of such host cells, reference to the "antibody" is a reference to one or more antibodies and equivalents thereof known to those skilled in the art, and so forth.

Unless defined otherwise, all technical and scientific terms used herein have the same meanings as commonly understood by one of ordinary skill in the art to which this invention belongs. Although any methods and materials similar or equivalent to those described herein can be used in the practice or testing of the present invention, the preferred methods, devices, and materials are now described. All publications mentioned herein are incorporated herein by reference for the purpose of describing and disclosing the cell lines, vectors, and methodologies that are reported in the publications that might be used in connection with the invention. Nothing herein is to be construed as an admission that the invention is not entitled to antedate such disclosure by virtue of prior invention.

"Acylate", as used herein, refers to the introduction of an acyl group into a molecule, (for example, acylation).

"Adjacent", as used herein, refers to a position in a nucleotide sequence immediately 5' or 3' to a defined sequence.

"Agonist", as used herein, refers to a molecule that, when bound to a polypeptide (for example, a polypeptide encoded by a nucleic acid of the present invention), increases the biological or immunological activity of the polypeptide. Agonists may include proteins, nucleic acids, carbohydrates, or any other molecules that bind to the protein.

"Allele" or "allelic sequence", as used herein, refers an alternative form of the gene that may result from at least one mutation in the nucleic acid sequence.

"Altered", as used herein, refers to modification in the metabolic profile compared to a reference or control where the amount of biochemical and/or chemical compound is increased or decreased.

"Alterations" in a polynucleotide (for example, a polypeptide encoded by a nucleic acid of the present invention), as used herein, comprise any deletions, insertions, and point mutations in the polynucleotide sequence. Included within this definition are alterations to the genomic DNA sequence that encodes the polypeptide.

"Amino acid sequence", as used herein, refers to an oligopeptide, peptide, polypeptide, or protein sequence, and fragments or portions thereof, and to naturally occurring or synthetic molecules. "Amino acid sequence" and like terms, such as "polypeptide" or "protein" as recited herein are not meant to limit the amino acid sequence to the complete, native amino acid sequence associated with the recited protein molecule.

"Amplification", as used herein, refers to the production of additional copies of a nucleic acid sequence and is generally carried out using polymerase chain reaction (PCR) technologies well known in the art (Dieffenbach, C. W. and G. S. Dveksler (1995) PCR Primer, a Laboratory Manual, Cold Spring Harbor Press, Plainview, N.Y.).

"Antibody" refers to intact molecules as well as fragments thereof that are capable of specific binding to a epitopic determinant. Antibodies that bind a polypeptide (for example, a polypeptide encoded by a nucleic acid of the present invention) can be prepared using intact polypeptides or fragments as the immunizing antigen. These antigens may be conjugated to a carrier protein, if desired.

"Antigenic determinant", "determinant group", or "epitope of an antigenic macromolecule", as used herein, refer to any region of the macromolecule with the ability or potential to elicit, and combine with, specific antibody. Determinants exposed on the surface of the macromolecule are likely to be immunodominant, that is, more immunogenic than other (immunorecessive) determinants that are less exposed, while some (for example, those within the molecule) are non-immunogenic (immunosilent). As used herein, "antigenic determinant" refers to that portion of a molecule that makes contact with a particular antibody (for example, an epitope). When a protein or fragment of a protein is used to immunize a host animal, numerous regions of the protein may induce the production of antibodies that bind specifically to a given region or three-dimensional structure on the protein; these regions or structures are referred to as antigenic determinants. An antigenic determinant may compete with the intact antigen (the immunogen used to elicit the immune response) for binding to an antibody.

"Antisense", as used herein, refers to a deoxyribonucleotide sequence whose sequence of deoxyribonucleotide residues is in reverse 5' to 3' orientation in relation to the sequence of deoxyribonucleotide residues in a sense strand of a DNA duplex. A "sense strand" of a DNA duplex refers to a strand in a DNA duplex that is transcribed by a cell in its natural state into a "sense mRNA". Thus an "antisense" sequence is a sequence having the same sequence as the non-coding strand in a DNA duplex. The term "antisense RNA" refers to a RNA transcript that is complementary to all or part of a target primary transcript or mRNA and that blocks the expression of a target gene by interfering with the processing, transport and/or translation of its primary transcript or mRNA. The complementarity of an antisense RNA may be with any part of the specific gene transcript, for example, at the 5' non-coding sequence, 3' non-coding sequence, introns, or the coding sequence. In addition, as used herein, antisense RNA may contain regions of ribozyme sequences that increase the efficacy of antisense RNA to block gene expression. "Ribozyme" refers to a catalytic RNA and includes sequence-specific endoribonucleases.

"Anti-sense inhibition", as used herein, refers to a type of gene regulation based on cytoplasmic, nuclear, or organelle inhibition of gene expression due to the presence in a cell of an RNA molecule complementary to at least a portion of the mRNA being translated. It is specifically contemplated that DNA molecules may be from either an RNA virus or mRNA from the host cell genome or from a DNA virus.

"Antagonist" or "inhibitor", as used herein, refer to a molecule that, when bound to a polypeptide (for example, a polypeptide encoded by a nucleic acid of the present invention), decreases the biological or immunological activity of the polypeptide. Antagonists and inhibitors may include proteins, nucleic acids, carbohydrates, or any other molecules that bind to the polypeptide.

"Biologically active", as used herein, refers to a molecule having the structural, regulatory, or biochemical functions of a naturally occurring molecule.

"Biological material", as used herein, refers to: a portion or portions of one or more cells, organs, or organisms; a whole cell, organelle, organ, or organism; or a group of cells, organelles, organs, or organisms. For example, if the organism(s) supplying the biological material is a garden variety carrot, a single leaf of one carrot plant could be used, or one or more whole carrot plant(s) could be used, or partial or whole taproots from a number of different individuals could be used, or mitochondria extracted from the crown of one carrot plant could be used.

"Cell culture", as used herein, refers to a proliferating mass of cells that may be in either an undifferentiated or differentiated state.

"Chimeric plasmid", as used herein, refers to any recombinant plasmid formed (by cloning techniques) from nucleic acids derived from organisms that do not normally exchange genetic information (for example, Escherichia coli and Saccharomyces cerevisiae).

"Chimeric sequence" or "chimeric gene", as used herein, refer to a nucleotide sequence derived from at least two heterologous parts. The sequence may comprise DNA or RNA.

"Chromatogram", as used herein, refers to an electronic and/or graphic record of data representing the absolutely or relatively quantitative detection of a plurality of separated chemical species obtained or derived from a group of metabolites, whether or not such separation has been performed by chromatography or some other method (e.g., electrophoresis).

"Control chromatogram", as used herein, refers to an individual chromatogram, or an average chromatograrn based on multiple individual chromatograms or a mathematical model based on multiple individual chromatograms, of chemical species obtained from a group of metabolites extracted from "control" biological material.

"Subject chromatogram", as used herein, refers to an individual chromatogram, or an average or model chromatogram based on multiple individual chromatograms, of chemical species obtained from a group of metabolites extracted from "subject" biological material. In either case, a model chromatogram may contain data including, e.g.: peak migration distance (or elution time) ranges and averages; peak height and peak area ranges and averages; and other parameters.

"Chromatographic data", as used herein, refers to chromatograms (e.g., including, but not limited to, total ion chromatograms or chromatograms generated from flame ionization detection) corresponding to individual biological or reference samples. Data such as retention time, retention index, peak areas, and peak areas normalized to internal standards can be extracted from total ion chromatograms to generate "peak tables."

"Coding sequence", as used herein, refers to a deoxyribonucleotide sequence that, when transcribed and translated, results in the formation of a cellular polypeptide or a ribonucleotide sequence that, when translated, results in the formation of a cellular polypeptide.

"Compatible", as used herein, refers to the capability of operating with other components of a system. A vector or plant viral nucleic acid that is compatible with a host is one that is capable of replicating in that host. A coat protein that is compatible with a viral nucleotide sequence is one capable of encapsidating that viral sequence.

"Coding region", as used herein, refers to that portion of a gene that codes for a protein. The term "non-coding region" refers to that portion of a gene that is not a coding region.

"Complementary" or "complementarity", as used herein, refer to the Watson-Crick base-pairing of two nucleic acid sequences. For example, for the sequence 5'-AGT-3' binds to the complementary sequence 3'-TCA-5'. Complementarity between two nucleic acid sequences may be "partial", in which only some of the bases bind to their complement, or it may be complete as when every base in the sequence binds to it's complementary base. The degree of complementarity between nucleic acid strands has significant effects on the efficiency and strength of hybridization between nucleic acid strands.

"Contig", as used herein, refers to a nucleic acid sequence that is derived from the contiguous assembly of two or more nucleic acid sequences.

"Control biological material" and "subject biological material", as used herein, both refer to biological material taken from (cultivated/domesticated or uncultivated/non-domesticated wild-type or genetically modified) individual(s) of any taxonomic category or categories, i.e. kingdom, phylum, subphylum, class, subclass, order, suborder, family, subfamily, genus, subgenus, species, subspecies, variety, breed, or strain. The "control" and "subject" biological material may be, and typically are, taken from individual(s) of the same taxonomic category, preferably from the same species, subspecies, variety, breed, or strain. However, when comparison between different types of organisms is desired, the "control" and "subject" biological material may be taken from individual(s) of different taxonomic categories. The "control" and "subject" biological materials differ from each other in at least one way. This difference may be that the "control" and "subject" biological materials were obtained from individual(s) of different taxonomic categories. Alternatively, or additionally, they may be different parts of the same organ(s), they may be different organelles or different groups of organelles, different cells or different groups of cells, different organs or different groups of organs, or different whole organisms or different groups of whole organisms. The difference may be that the organisms providing the biological materials are identical, but for, e.g., their growth stages.

"Correlates with expression of a polynucleotide", as used herein, indicates that the detection of the presence of ribonucleic acid that is similar to a nucleic acid (for example, SEQ ID NOs:1-7554) and is indicative of the presence of mRNA encoding a polypeptide (for example, a polypeptide encoded by a nucleic acid of the present invention) in a sample and thereby correlates with expression of the transcript from the polynucleotide encoding the protein.

"Customized reporting", as used herein, refers to the modification of a preliminary analyst report to generate an interim report (e.g., including, but not limited to, a modified analyst report and a cross-referenced modified analyst report) and a final report. In some embodiments, modifications include, but are not limited to, substitution of underivatized compound names for derivatized compound names and generation of a hit score. In other embodiments, customized reporting includes data mining of databases to generate biochemical profiling and genetic expression information and/or reports.

"Data analysis and reporting software", as used herein, refers to software configured for the analysis of spectroscopic and chromatographic data corresponding to biological subject and reference samples. Data analysis and reporting software is configured to perform data reduction, two-dimensional peak matching, quantitative peak differentiation, peak identification, querying, data mining, and customized reporting functions.

"Data reduction", as used herein, refers to the process of organizing, compiling, and normalizing data (for example, chromatographic and spectroscopic data). In some embodiments, data reduction includes the normalization of raw chromatogram peak areas and the generation of peak tables. In some embodiments, data reduction also includes the process of filtering peaks based on their normalized area. This step removes peaks that are considered to be background.

"Data sorting", as used herein, refers to the generation of a preliminary analyst report. In some embodiment, the preliminary analyst report can include equivalence value, retention time, retention index, normalized peak area, peak identification status, compound name or other unique identifier, compound identification number (e.g., a CAS number), mass spectral library name, ID number, MS-XCR value, relative % change, notes, and other information about the biological sample.

"Data mining", as used herein, refers to the process of querying and mining databases to analyze and to obtain information (e.g., to use in the generation of customized reports of information pertaining to biochemical profiling and gene function and expression).

"Deletion", as used herein, refers to a change made in either an amino acid or nucleotide sequence resulting in the absence of one or more amino acids or nucleotides, respectively.

"Encapsidation", as used herein, refers to the process during virion assembly in which nucleic acid becomes incorporated in the viral capsid or in a head/capsid precursor (for example, in certain bacteriophages).

"Exon", as used herein, refers to a polynucleotide sequence in a nucleic acid that encodes information for protein synthesis and that is copied and spliced together with other such sequences to form messenger RNA.

"Expression", as used herein, is meant to incorporate transcription, reverse transcription, and translation.

"Expressed sequence tag (EST)" as used herein, refers to relatively short single-pass DNA sequences obtained from one or more ends of cDNA clones and RNA derived therefrom. They may be present in either the 5' or the 3' orientation. ESTs have been shown to be useful for identifying particular genes.

"Fractionated biological sample", as used herein, refers to a biological sample that has been fractionated into two or more fractions based on one or more properties of the sample. For example, in some embodiments, leaf extracts are fractionated based on extraction with organic solvents.

"Industrial crop", as used herein, refers to crops grown primarily for consumption by humans or animals or use in industrial processes (for example, as a source of fatty acids for manufacturing or sugars for producing alcohol). It will be understood that either the plant or a product produced from the plant (for example, sweeteners, oil, flour, or meal) can be consumed. Examples of food crops include, but are not limited to, corn, soybean, rice, wheat, oilseed rape, cotton, oats, barley, and potato plants.

"Foreign gene", as used herein, refers to any sequence that is not native to the organism.

"Fusion protein", as used herein, refers to a protein containing amino acid sequences from each of two distinct proteins; it is formed by the expression of a recombinant gene in which two coding sequences have been joined together such that their reading frames are in phase. Hybrid genes of this type may be constructed in vitro in order to label the product of a particular gene with a protein that can be more readily assayed (for example, a gene fused with lacZ in E. coli to obtain a fusion protein with .beta.-galactosidase activity). Alternatively, a protein may be linked to a signal peptide to allow its secretion by the cell. The products of certain viral oncogenes are fusion proteins.

"Gene", as used herein, refers to a discrete nucleic acid sequence responsible for a discrete cellular product. The term "gene", as used herein, refers not only to the nucleotide sequence encoding a specific protein, but also to any adjacent 5' and 3' non-coding nucleotide sequence involved in the regulation of expression of the protein encoded by the gene of interest. These non-coding sequences include terminator sequences, promoter sequences, upstream activator sequences, regulatory protein binding sequences, and the like. These non-coding sequence gene regions may be readily identified by comparison with previously identified eukaryotic non-coding sequence gene regions. Furthermore, the person of average skill in the art of molecular biology is able to identify the nucleotide sequences forming the non-coding regions of a gene using well-known techniques such as a site-directed mutagenesis, sequential deletion, promoter probe vectors, and the like.

"Genetically modified" and "genetically unmodified" when used in relation to subject biological material and control biological material, respectively, refer to the fact that the subject biological material has been treated to produce a genetic modification thereof, whereas the control biological material has not received that particular genetic modification. In this context, the term "genetically unmodified" does not imply that the "control" biological material must be, e.g., a naturally-occurring, wild-type plant; rather, both the control and subject biological materials may be (but need not be) the result of, e.g., hybridization, selection, or genetic engineering.

"Growth cycle", as used herein, is meant to include the replication of a nucleus, an organelle, a cell, or an organism.

"Heterologous", as used herein, refers to the association of a molecular or genetic element associated with a distinctly different type of molecular or genetic element.

"HIT" and/or "hit", as used herein, refers to the result of a test or series of tests that meets a defined criteria for each test. "Hit Detection", as used herein, refers to the process of determing a hit using a mathemetical or statistical model.

"Host", as used herein, refers to a cell, tissue or organism capable of replicating a vector or plant viral nucleic acid and that is capable of being infected by a virus containing the viral vector or plant viral nucleic acid. This term is intended to include prokaryotic and eukaryotic cells, organs, tissues or organisms, where appropriate.

The term "homolog" as in a "homolog" of a given nucleic acid sequence, as used herein, refers to a nucleic acid sequence (for example, a nucleic acid sequence from another organism), that shares a given degree of "homology" with the nucleic acid sequence.

"Homology", as used herein, refers to a degree of complementarity. There may be partial homology or complete homology (identity). A partially complementary sequence is one that at least partially inhibits a completely complementary sequence from hybridizing to a target nucleic acid and is referred to using the functional term "substantially homologous". The inhibition of hybridization of the completely complementary sequence to the target sequence may be examined using a hybridization assay (Southern or Northern blot, solution hybridization and the like) under conditions of low stringency. A substantially homologous sequence or probe will compete for and inhibit the binding (the hybridization) of a completely homologous sequence to a target under conditions of low stringency. This is not to say that conditions of low stringency are such that non-specific binding is permitted; low stringency conditions require that the binding of two sequences to one another be a specific (selective) interaction. The absence of non-specific binding may be tested by the use of a second target that lacks even a partial degree of complementarity (for example, less than about 30% identity); in the absence of non-specific binding the probe will not hybridize to the second non-complementary target.

Numerous equivalent conditions may be employed to comprise low stringency conditions; factors such as the length and nature (DNA, RNA, base composition) of the probe and nature of the target (DNA, RNA, base composition, present in solution or immobilized, etc.) and the concentration of the salts and other components (for example, the presence or absence of formamide, dextran sulfate, polyethylene glycol) are considered and the hybridization solution may be varied to generate conditions of low stringency hybridization different from, but equivalent to, the above listed conditions. In addition, the art knows conditions that promote hybridization under conditions of high stringency (for example, increasing the temperature of the hybridization and/or wash steps, the use of formamide in the hybridization solution, etc.).

When used in reference to a double-stranded nucleic acid sequence such as a cDNA or genomic clone, the term "substantially homologous" refers to any probe that can hybridize to either or both strands of the double-stranded nucleic acid sequence under conditions of low stringency as described above.

A gene may produce multiple RNA species that are generated by differential splicing of the primary RNA transcript. cDNAs that are splice variants of the same gene will contain regions of sequence identity or complete homology (representing the presence of the same exon or portion of the same exon on both cDNAs) and regions of complete non-identity (for example, representing the presence of exon "A" on cDNA 1 wherein cDNA 2 contains exon "B" instead). Because the two cDNAs contain regions of sequence identity they will both hybridize to a probe derived from the entire gene or portions of the gene containing sequences found on both cDNAs; the two splice variants are therefore substantially homologous to such a probe and to each other.

When used in reference to a single-stranded nucleic acid sequence, the term "substantially homologous" refers to any probe that can hybridize (it is the complement of) the single-stranded nucleic acid sequence under conditions of low stringency as described above.

The term "hybridization" is used in reference to the pairing of complementary nucleic acids. Hybridization and the strength of hybridization (for example, the strength of the association between the nucleic acids) is impacted by such factors as the degree of complementary between the nucleic acids, stringency of the conditions involved, the T.sub.m of the formed hybrid, and the G:C ratio within the nucleic acids.

"Hybridization complex", as used herein, refers to a complex formed between nucleic acid strands by virtue of hydrogen bonding, stacking or other non-covalent interactions between bases. A hybridization complex may be formed in solution or between nucleic acid sequences present in solution and nucleic acid sequences immobilized on a solid support (for example, membranes, filters, chips, pins or glass slides to which cells have been fixed for in situ hybridization).

"Immunologically active", as used herein, refers to the capability of a natural, recombinant, or synthetic polypeptide, or any oligopeptide thereof, to bind with specific antibodies and induce a specific immune response in appropriate animals or cells.

"Induction" and the terms "induce", "induction" and "inducible", as used herein, refer generally to a gene and a promoter operably linked thereto which is in some manner dependent upon an external stimulus, such as a molecule, in order to actively transcribed and/or translate the gene.

"Infection", as used herein, refers to the ability of a virus to transfer its nucleic acid to a host or introduce viral nucleic acid into a host, wherein the viral nucleic acid is replicated, viral proteins are synthesized, and new viral particles assembled. In this context, the terms "transmissible" and "infective" are used interchangeably herein.

"Insertion" or "addition", as used herein, refers to the replacement or addition of one or more nucleotides or amino acids, to a nucleotide or amino acid sequence, respectively.

"In cis", as used herein, indicates that two sequences are positioned on the same strand of RNA or DNA.

"In trans", as used herein, indicates that two sequences are positioned on different strands of RNA or DNA.

"Intron", as used herein, refers to a polynucleotide sequence in a nucleic acid that does not encode information for protein synthesis and is removed before translation of messenger RNA.

"Isolated", as used herein, refers to a polypeptide or polynucleotide molecule separated not only from other peptides, DNAs, or RNAs, respectively, that are present in the natural source of the macromolecule. "Isolated" and "purified" do not encompass either natural materials in their native state or natural materials that have been separated into components (for example, in an acrylamide gel) but not obtained either as pure substances or as solutions.

"Kinase", as used herein, refers to an enzyme (for example, hexokinase and pyruvate kinase) that catalyzes the transfer of a phosphate group from one substrate (commonly ATP) to another.

"Marker" or "genetic marker", as used herein, refer to a genetic locus that is associated with a particular, usually readily detectable, genotype or phenotypic characteristic (for example, an antibiotic resistance gene).

"Metabolic characteristics", as used herein, refers to a biochemical/chemical trait/metabolite that is genetically expressed in a biological system. "Altered metabolic characteristic", as used herein, refers to the production of a given metabolite that has been altered (for example, increased or decreased) in a biological system, especially plants. Examples of metabolites that can be altered in a plant include, but are not limited to, acids, fatty acids, amino acids, hydroxy fatty acids, branched fatty acids, carbohydrates, hydrocarbons, glycerides, phenols, strerols, oxygenated terpenes, and other isoprenoids, alcohols, ketones, quinones, alkenes and alkynes.

"Metabolome", as used herein, indicates the complement of relatively low molecular weight molecules that is present in a plant, plant part, or plant sample, or in a suspension or extract thereof. Examples of such molecules include, but are not limited to: acids and related compounds; mono-, di-,and tri-carboxylic acids (saturated, unsaturated, aliphatic and cyclic, aryl, alkaryl); aldo-acids, keto-acids; lactone forms; gibberellins; abscisic acid; alcohols, polyols, derivatives, and related compounds; ethyl alcohol, benzyl alcohol, methanol; propylene glycol, glycerol, phytol; inositol, furfuryl alcohol, menthol; aldehydes, ketones, quinones, derivatives, and related compounds; acetaldehyde, butyraldehyde, benzaldehyde, acrolein, furfural, glyoxal; acetone, butanone; anthraquinone; carbohydrates; mono-, di-, tri-saccharides; alkaloids, amines, and other bases; pyridines (including nicotinic acid, nicotinamide); pyrimidines (including cytidine, thymine); purines (including guanine, adenine, xanthines/hypoxanthines, kinetin); pyrroles; quinolines (including isoquinolines); morphinans, tropanes, cinchonans; nucleotides, oligonucleotides, derivatives, and related compounds; guanosine, cytosine, adenosine, thymidine, inosine; amino acids, oligopeptides, derivatives, and related compounds; esters; phenols and related compounds; heterocyclic compounds and derivatives; pyrroles, tetrapyrroles (corrinoids and porphines/porphyrins, w/w/o metal-ion); flavonoids; indoles; lipids (including fatty acids and triglycerides), derivatives, and related compounds; carotenoids, phytoene; and sterols, isoprenoids including terpenes.

"Modulate", as used herein, refers to a change or an alteration in the biological activity of a polypeptide (for example, a polypeptide encoded by a nucleic acid of the present invention). Modulation may be an increase or a decrease in protein activity, a change in binding characteristics, or any other change in the biological, functional or immunological properties of the polypeptide.

"Movement protein", as used herein, refers to a noncapsid protein required for cell to cell movement of replicons or viruses in plants.

"Multigene family", as used herein, refers to a set of genes descended by duplication and variation from some ancestral gene. Such genes may be clustered together on the same chromosome or dispersed on different chromosomes. Examples of multigene families include those that encode the histones, hemoglobins, immunoglobulins, histocompatibility antigens, actins, tubulins, keratins, collagens, heat shock proteins, salivary glue proteins, chorion proteins, cuticle proteins, yolk proteins, and phaseolins.

"Non-native", as used herein, refers to any RNA sequence that promotes production of subgenomic mRNA including, but not limited to, 1) plant viral promoters such as ORSV and brome mosaic virus, 2) viral promoters from other organisms such as human Sindbis viral promoter, and 3) synthetic promoters.

"Nucleic acid sequence", as used herein, refers to a polymer of nucleotides in which the 3' position of one nucleotide sugar is linked to the 5' position of the next by a phosphodiester bridge. In a linear nucleic acid strand, one end typically has a free 5' phosphate group, the other a free 3' hydroxyl group. Nucleic acid sequences may be used herein to refer to oligonucleotides, or polynucleotides, and fragments or portions thereof, and to DNA or RNA of genomic or synthetic origin that may be single- or double-stranded, and represent the sense or antisense strand.

"Polypeptide", as used herein, refers to an amino acid sequence obtained from any species and from any source whether natural, synthetic, semi-synthetic, or recombinant.

"Principal component analysis", as used herein, refers to algorithms designed to represent large and complex data sets by linear combinations of the original variables. These linear combinations of variables are extracted to maximize the explained variance and are mutually orthogonal. Principal component analysis is designed to represent large complex data sets by linear combinations of the original variables that maximize the explained variance and are mutually orthogonal.

"Oil-producing species", as used herein, refers to plant species that produce and store triacylglycerol in specific organs, primarily in seeds. Such species include soybean (Glycine max), rapeseed and canola (including Brassica napus and B. campestris), sunflower (Helianthus annus), cotton (Gossypium hirsutum), corn (Zea mays), cocoa (Theobroma cacao), safflower (Carthamus tinctorius), oil palm (Elaeis guineensis), coconut palm (Cocos nucifera), flax (Linum usitatissimum), castor (Ricinus communis) and peanut (Arachis hypogaea). The group also includes non-agronomic species that are useful in developing appropriate expression vectors such as tobacco, rapid cycling Brassica species, and Arabidopsis thaliana, and wild species that may be a source of unique fatty acids.

"Operably linked", as used herein, refers to a juxtaposition of components, particularly nucleotide sequences, such that the normal function of the components can be performed. Thus, a coding sequence that is operably linked to regulatory sequences refers to a configuration of nucleotide sequences wherein the coding sequences can be expressed under the regulatory control, that is, transcriptional and/or translational control, of the regulatory sequences.

"Origin of assembly", as used herein, refers to a sequence where self-assembly of the viral RNA and the viral capsid protein initiates to form virions.

"Ortholog", as used herein, refers to genes that have evolved from an ancestral locus.

"Outlier peak", as used herein, indicates a peak of a chromatogram of a test sample, or the relative or absolute detected response data, or amount or concentration data thereof. An outlier peak: 1) may have a significantly different peak height or area as compared to a like chromatogram of a control sample; or 2) be an additional or missing peak as compared to a like chromatogram of a control sample.

"Overexpression", as used herein, refers to the production of a gene product in transgenic organisms that exceeds levels of production in normal or non-transformed organisms.

"Cosuppression", as used herein, refers to the expression of a foreign gene that has substantial homology to an endogenous gene resulting in the suppression of expression of both the foreign and the endogenous gene. As used herein, the term "altered levels" refers to the production of gene product(s) in transgenic organisms in amounts or portions that differ from that of normal or non-transformed organisms.

"Peak identification", as used herein, refers to the characterization and identification of a chemical compound represented by a given chromatographic peak. In some embodiments, the chemical compound corresponding to a given peak is identified by searching mass spectral libraries. In other embodiments, the chemical compounds are identified by searching additional libraries or databases (for example, biotechnology databases).

"Quantitative peak differentiation", as used herein, refers to the process of confirming matched peaks by calculating their relative quantitative differentiation, which is expressed as a percent change of the subject peak area relative to the area of the reference peak. A predetermined threshold for change is used to confirm that the peaks are of significant biological alteration.

"Plant", as used herein, refers to any plant and progeny thereof. The term also includes parts of plants, including seed, cuttings, tubers, fruit, flowers, etc.

"Plant cell", as used herein, refers to the structural and physiological unit of plants, consisting of a protoplast and the cell wall.

"Plant organ", as used herein, refers to a distinct and visibly differentiated part of a plant, such as root, stem, leaf or embryo.

"Plant tissue", as used herein, refers to any tissue of a plant in planta or in culture. This term is intended to include a whole plant, plant cell, plant organ, protoplast, cell culture, or any group of plant cells organized into a structural and functional unit.

"Portion", as used herein, with regard to a protein ("a portion of a given protein") refers to fragments of that protein. The fragments may range in size from four amino acid residues to the entire amino acid sequence minus one amino acid (10 nucleotides, 20, 30, 40, 50, 100, 200, etc.). A "portion" is preferably at least 25 nucleotides, more preferably at least 50 nucleotides, and even more preferably at least 100 nucleotides.

"Positive-sense inhibition", as used herein, refers to a type of gene regulation based on cytoplasmic inhibition of gene expression due to the presence in a cell of an RNA molecule substantially homologous to at least a portion of the mRNA being translated.

"Production cell", as used herein, refers to a cell, tissue or organism capable of replicating a vector or a viral vector, but which is not necessarily a host to the virus. This term is intended to include prokaryotic and eukaryotic cells, organs, tissues or organisms, such as bacteria, yeast, fungus, and plant tissue.

"Promoter", as used herein, refers to the 5'-flanking, non-coding sequence adjacent a coding sequence that is involved in the initiation of transcription of the coding sequence.

"Protoplast", as used herein, refers to an isolated plant cell without cell walls, having the potency for regeneration into cell culture or a whole plant.

"Purified", as used herein, when referring to a peptide or nucleotide sequence, indicates that the molecule is present in the substantial absence of other biological macromolecular, for example, polypeptides, polynucleic acids, and the like of the same type. The term "purified" as used herein preferably means at least 95% by weight, more preferably at least 99.8% by weight, of biological macromolecules of the same type present (but water, buffers, and other small molecules, especially molecules having a molecular weight of less than 1000 can be present).

"Pure", as used herein, preferably has the same numerical limits as "purified" immediately above. "Substantially purified", as used herein, refers to nucleic or amino acid sequences that are removed from their natural environment, isolated or separated, and are at least 60% free, preferably 75% free, and most preferably 90% free from other components with which they are naturally associated.

"Recombinant plant viral nucleic acid", as used herein, refers to a plant viral nucleic acid that has been modified to contain non-native nucleic acid sequences. These non-native nucleic acid sequences may be from any organism or purely synthetic, however, they may also include nucleic acid sequences naturally occurring in the organism into which the recombinant plant viral nucleic acid is to be introduced.

"Recombinant plant virus", as used herein, refers to a plant virus containing a recombinant plant viral nucleic acid.

"Reference sample", as used herein, refers to a sample taken from an individual receiving treatment that is not believed to alter the chemistry thereof.

"Regulated", as used herein, refers to an alteration that occurs in an expressed metabolite in a biological system. "Up-regulated", as used herein, refers to an increase in a give metabolite level relative to a control or reference. "Down-regulated", as used herein, refers to a decrease in a given metabolite level relative to a control or reference.

"Regulatory region" or "regulatory sequence", as used herein, in reference to a specific gene refers to the non-coding nucleotide sequences within that gene that are necessary or sufficient to provide for the regulated expression of the coding region of a gene. Thus the term regulatory region includes promoter sequences, regulatory protein binding sites, upstream activator sequences, and the like. Specific nucleotides within a regulatory region may serve multiple functions. For example, a specific nucleotide may be part of a promoter and participate in the binding of a transcriptional activator protein.

"Replication origin", as used herein, refers to the minimal terminal sequences in linear viruses that are necessary for viral replication.

"Replicon", as used herein, refers to an arrangement of RNA sequences generated by transcription of a transgene that is integrated into the host DNA that is capable of replication in the presence of a helper virus. A replicon may require sequences in addition to the replication origins for efficient replication and stability.

"Sample", as used herein, is used in its broadest sense. A biological sample suspected of containing nucleic acid encoding a polypeptide (for example, a polypeptide encoded by a nucleic acid of the present invention) or fragments thereof may comprise a tissue, a cell, an extract from cells, chromosomes isolated from a cell (for example, a spread of metaphase chromosomes), genomic DNA (in solution or bound to a solid support such as for Southern analysis), RNA (in solution or bound to a solid support such as for northern analysis), cDNA (in solution or bound to a solid support), and the like.

"Site-directed mutagenesis", as used herein, refers to the in-vitro induction of mutagenesis at a specific site in a given target nucleic acid molecule.

"Subgenomic promoter", as used herein, refers to a promoter of a subgenomic mRNA of a viral nucleic acid.

"Subject sample", as used herein, refers to a sample taken from an individual that has been treated in order to alter the chemistry thereof.

"T.sub.m" is used in reference to the "melting temperature". The melting temperature is the temperature at which a population of double-stranded nucleic acid molecules becomes half-dissociated into single strands. The equation for calculating the T.sub.m of nucleic acids is well known in the art. As indicated by standard references, a simple estimate of the T.sub.m value may be calculated by the equation: T.sub.m=81.5+0.41(% G+C), when a nucleic acid is in aqueous solution at 1 M NaCl [See for example, Anderson and Young, Quantitative Filter Hybridization, in Nucleic Acid Hybridization (1985)]. Other references include more sophisticated computations that take structural as well as sequence characteristics into account for the calculation of T.sub.m.

"Stringency" is used in reference to the conditions of temperature, ionic strength, and the presence of other compounds such as organic solvents, under which nucleic acid hybridizations are conducted. Those skilled in the art will recognize that "stringency" conditions may be changed by varying the parameters just described either individually or in concert. With "high stringency" conditions, nucleic acid base pairing will occur only between nucleic acid fragments that have a high frequency of complementary base sequences (for example, hybridization under "high stringency" conditions may occur between homologs with about 85-100% identity, preferably about 70-100% identity). With medium stringency conditions, nucleic acid base pairing will occur between nucleic acids with an intermediate frequency of complementary base sequences (for example, hybridization under "medium stringency" conditions may occur between homologs with about 50-70% identity). Thus, conditions of "weak" or "low" stringency are often required with nucleic acids that are derived from organisms that are genetically diverse, as the frequency of complementary sequences is usually less.

"High stringency conditions" when used in reference to nucleic acid hybridization comprise conditions equivalent to binding or hybridization at 42.degree. C. in a solution consisting of 5.times.SSPE (43.8 g/l NaCl, 6.9 g/l NaH.sub.2PO.sub.4 H.sub.2O and 1.85 g/l EDTA, pH adjusted to 7.4 with NaOH), 0.5% SDS, 5.times. Denhardt's reagent and 100 .mu.g/mL denatured salmon sperm DNA followed by washing in a solution comprising 0.1.times.SSPE, 1.0% SDS at 42.degree. C. when a probe of about 500 nucleotides in length is employed.

"Medium stringency conditions" when used in reference to nucleic acid hybridization comprise conditions equivalent to binding or hybridization at 42.degree. C. in a solution consisting of 5.times.SSPE (43.8 g/l NaCl, 6.9 g/l NaH.sub.2PO.sub.4 H.sub.2O and 1.85 g/l EDTA, pH adjusted to 7.4 with NaOH), 0.5% SDS, 5.times. Denhardt's reagent and 100 .mu.g/mL denatured salmon sperm DNA followed by washing in a solution comprising 1.1.times.SSPE, 1.0% SDS at 42.degree. C. when a probe of about 500 nucleotides in length is employed. "Low stringency conditions" when used in a reference to nucleic acid hybridization comprise conditions equivalent to binding or hybridization at 42.degree. C. in a solution consisting of 5.times.SSPE (43.8 g/l NaCl, 6.9 g/l NaH.sub.2PO.sub.4 H.sub.2O and 1.85 g/l EDTA, pH adjusted to 7.4 with NaOH), 0.1% SDS, 5.times. Denhardt's reagent [50.times. Denhardt's contains per 500 mL: 5 g Ficoll (Type 400, Pharamcia), 5 g BSA (Fraction V; Sigma)] and 100 .mu.g/mL denatured salmon sperm DNA followed by washing in a solution comprising 5.times.SSPE, 0.1% SDS at 42.degree. C. when a probe of about 500 nucleotides in length is employed.

"Substitution", as used herein, refers to a change made in an amino acid of nucleotide sequence that results in the replacement of one or more amino acids or nucleotides by different amino acids or nucleotides, respectively.

"Symptom", as used herein refers to a visual condition resulting from the action of the GENEWARE.TM. (trademark of Large Scale Biology Corporation) vector or the clone insert. The GENEWARE.TM. vector is described in U.S. application Ser. No. 09/008,186 (incorporated herein by reference).

"Systemic infection", as used herein, denotes infection throughout a substantial part of an organism including mechanisms of spread other than mere direct cell inoculation but rather including transport from one infected cell to additional cells either nearby or distant.

"Transcription", as used herein, refers to the production of an RNA molecule by RNA polymerase as a complementary copy of a DNA sequence.

"Transformation", as used herein, describes a process by which exogenous DNA enters and changes a recipient cell. It may occur under natural or artificial conditions using various methods well known in the art. Transformation may rely on any known method for the insertion of foreign nucleic acid sequences into a prokaryotic or eukaryotic host cell. The method is selected based on the host cell being transformed and may include, but is not limited to, viral infection, electroporation, lipofection, and particle bombardment. Such "transformed" cells include stably transformed cells in which the inserted DNA is capable of replication either as an autonomously replicating plasmid or as part of the host chromosome. They also include cells that transiently express the inserted DNA or RNA for limited periods of time.

"Transfection", as used herein, refers to the introduction of foreign nucleic acid into eukaryotic cells. Transfection may be accomplished by a variety of means known to the art including calcium phosphate-DNA co-precipitation, DEAE-dextran-mediated transfection, polybrene-mediated transfection, electroporation, microinjection, liposome fusion, lipofection, protoplast fusion, retroviral infection, and biolistics. Transfection may, for example, result in cells in which the inserted nucleic acid is capable of replication either as an autonomously replicating molecule or as part of the host chromosome, or cells that transiently express the inserted nucleic acid for limited periods of time.

"Transgenic plant", as used herein, refers to a plant that contains a foreign nucleotide sequence inserted into either its nuclear genome or organellar genome.

"Transgene", as used herein, refers to the DNA sequence coding for the replicon that is inserted into the host DNA.

"Two-dimensional peak matching", as used herein, refers to the pairing or matching of peaks in reference and subject biological samples. Peaks are first paired based on their retention index. A match is then confirmed by spectral matching.

"Unmatched peak", as used herein, refers to a peak reported in the chromatographic and/or spectroscopic data corresponding to reference biological sample but missing from chromatographic and/or spectroscopic data corresponding to subject biological sample, based upon the criteria for quantitation and reporting, or a peak reported in chromatographic and/or spectroscopic data corresponding to subject biological sample but missing from chromatographic and/or spectroscopic data corresponding to reference biological sample, based upon criteria for quantitation and reporting.

"Variants" of a polypeptide (for example, a polypeptide encoded by a nucleic acid of the present invention), as used herein, refers to a sequence resulting when a polypeptide is modified by one or more amino acids. The variant may have "conservative" changes, wherein a substituted amino acid has similar structural or chemical properties, for example, replacement of leucine with isoleucine. More rarely, a variant may have "nonconservative" changes, for example, replacement of a glycine with a tryptophan. Variants may also include sequences with amino acid deletions or insertions, or both. Guidance in determining which amino acid residues may be substituted, inserted, or deleted without abolishing biological or immunological activity may be found using computer programs well known in the art.

"Vector", as used herein, refers to a self-replicating DNA or RNA molecule that transfers a nucleic acid segment between cells.

"Virion", as used herein, refers to a particle composed of viral RNA and viral capsid protein.

"Virus", as used herein, refers to an infectious agent composed of a nucleic acid encapsidated in a protein. A virus may be a mono-, di-, tri- or multi-partite virus.

DESCRIPTION OF THE INVENTION

I. Identification of Nucleotide and Amino Acid Sequences

The invention is based on the discovery of deoxyribonucleic acid (DNA) and amino acid sequences that confer altered metabolic characteristics when expressed in plants. In particular, the present invention encompasses the nucleic acid sequences encoded by SEQ ID NOs:1-1165, 3703-4153 and 7389-7554 and variants and portions thereof. These sequences are contiguous sequences prepared from a database of 5' single pass sequences and are thus referred to as contig sequences.

Nucleic acids of the present invention were identified in clones generated from a variety of cDNA libraries. The cDNA libraries were constructed in the GENEWARE.TM. vector. The GENEWARE.TM. vector is described in U.S. application Ser. No. 09/008,186 (incorporated herein by reference). Each of the complete set of clones from the GENEWARE.TM. library was used to prepare an infectious viral unit. An infectious unit corresponding to each clone was used to inoculate Nicotiana benthamiana (a dicotyledonous plant). The plants were grown under identical conditions and a phenotypic analysis of each plant was carried out. The altered metabolic characteristic was observed in the plants that had been infected by an infectious unit created from the nucleic acids of the present invention.

Following the identification of the altered metabolic characteristic in plant samples, further analyses of the sequences were carried out. In particular, the nucleotide sequences of the present invention were analyzed using bioinformatics methods as described below.

II. Bioinformatics Methods

A. Phred, Phrap and Consed

Phred, Phrap and Consed are a set of programs that read DNA sequencer traces, make base calls, assemble the shotgun DNA sequence data and analyze the sequence regions that are likely to contribute to errors. Phred is the initial program used to read the sequencer trace data, call the bases and assign quality values to the bases. Phred uses a Fourier-based method to examine the base traces generated by the sequencer. The output files from Phred are written in FASTA, phd or scf format. Phrap is used to assemble contiguous sequences from only the highest quality portion of the sequence data output by Phred. Phrap is amenable to high-throughput data collection. Finally, Consed is used as a finishing tool to assign error probabilities to the sequence data. Detailed descriptions of the Phred, Phrap and Consed software and its use can be found in the following references: Ewing et al., Genome Res., 8:175 [1998]; Ewing and Green, Genome Res. 8:186 [1998]; Gordon et al., Genome Res. 8: 195 [1998].

B. BLAST

The BLAST set of programs may be used to compare the large numbers of sequences and obtain homologies to known protein families. These homologies provide information regarding the function of newly sequenced genes. Detailed descriptions of the BLAST software and its uses can be found in the following references Altschul et al., J. Mol. Biol., 215:403 [1990]; Altschul, J. Mol. Biol. 219:555 [1991].

Generally, BLAST performs sequence similarity searching and is divided into 5 basic subroutines: (1) BLASTP compares an amino acid sequence to a protein sequence database; (2) BLASTS compares a nucleotide sequence to a nucleic acid sequence database; (3) BLASTX compares translated protein sequences done in 6 frames to a protein sequence database; (4) TBLASTN compares a protein sequence to a nucleotide sequence database that is translated into all 6 reading frames; (5) TBLASTX compares the 6 frame translated protein sequence to the 6-frame translation of a nucleotide sequence database. Subroutines (3)-(5) may be used to identify weak similarities in nucleic acid sequence.

The BLAST program is based on the High Segment Pair (HSP), two sequence fragments of arbitrary but equal length whose alignment is locally maximized and whose alignment meets or exceeds a cutoff threshold. BLAST determines multiple HSP sets statistically using sum statistics. The score of the HSP is then related to its expected chance of frequency of occurrence, E. The value, E, is dependent on several factors such as the scoring system, residue composition of sequences, length of query sequence and total length of database. In the output file will be listed these E values, typically in a histogram format, which are useful in determining levels of statistical significance at the user s predefined expectation threshold. Finally, the Smallest Sum Probability, P(N) is the probability of observing the shown matched sequences by chance alone and is typically in the range of 0-1.

BLAST measures sequence similarity using a matrix of similarity scores for all possible pairs of residues and these specify scores for aligning pairs of amino acids. The matrix of choice for a specific use depends on several factors: the length of the query sequence and whether or not a close or distant relationship between sequences is suspected. Several matrices are available including PAM40, PAM120, PAM250, BLOSUM 62 and BLOSUM 50. Altschul et al. (1990) found PAM120 to be the most broadly sensitive matrix (for example point accepted mutation matrix per 100 residues). However, in some cases the PAM120 matrix may not find short but strong or long but weak similarities between sequences. In these cases, pairs of PAM matrices may be used, such as PAM40 and PAM 250, and the results compared. Typically, PAM 40 is used for database searching with a query of 9-21 residues long, while PAM 250 is used for lengths of 47-123.

The BLOSUM (Blocks Substitution Matrix) series of matrices are constructed based on percent identity between two sequence segments of interest. Thus, the BLOSUM62 matrix is based on a matrix of sequence segments in which the members are less than 62% identical. BLOSUM62 shows very good performance for BLAST searching. However, other BLOSUM matrices, like the PAM matrices, may be useful in other applications. For example, BLOSUM45 is particularly strong in profile searching.

C. FASTA

The FASTA suite of programs permits the evaluation of DNA and protein similarity based on local sequence alignment. The FASTA search algorithm utilizes Smith/Waterman- and Needleman/Wunsch-based optimization methods. These algorithms consider all of the alignment possibilities between the query sequence and the library in the highest-scoring sequence regions. The search algorithm proceeds in four basic steps: 1. The identities or pairs of identities between the two DNA or protein sequences are determined. The ktup parameter, as set by the user, is operative and determines how many consecutive sequence identities are required to indicate a match. 2. The regions identified in step I are re-scored using a PAM or BLOSUM matrix. This allows conservative replacements and runs of identities shorter than that specified by ktup to contribute to the similarity score. 3. The region with the single best scoring initial region is used to characterize pairwise similarity and these scores are used to rank the library sequences. 4. The highest scoring library sequences are aligned using the Smith-Waterman algorithm. This final comparison takes into account the possible alignments of the query and library sequence in the highest scoring region.

Further detailed description of the FASTA software and its use can be found in the following reference: Pearson and Lipman, Proc. Natl. Acad. Sci., 85: 2444 [1988].

D. Pfam

Despite the large number of different protein sequences determined through genomics-based approaches, relatively few structural and functional domains are known. Pfam is a computational method that utilizes a collection of multiple alignments and profile hidden Markov models of protein domain families to classify existing and newly found protein sequences into structural families. Detailed descriptions of the Pfam software and its uses can be found in the following references: Sonhammer et al., Proteins: Structure, Function and Genetics, 28:405 [1997]; Sonhammer et al., Nucleic Acids Res., 26:320 [1998]; Bateman et al., Nucleic Acids Res., 27: 260 [1999].

Pfam 3.1, the latest version, includes 54% of proteins in SWISS_PROT [For a recent reference see: Barker W. C., Garavelli J. S., Hou Z., Huang H., Ledley R. S., McGarvey P. B., Mewes H.-W., Orcutt B. C., Pfeiffer F., Tsugita A., Vinayaka C. R., Xiao C., Yeh L. S., Wu C.; Nucleic Acids Res. 29:29-32(2001)] and SP-TrEMBL-5 (A supplement to SWISS_PROT) as a match to the database and includes expectation values for matches. Pfam consists of parts A and B. Pfam-A contains a hidden Markov model and includes curated families. Pfam-B uses the Domainer program to cluster sequence segments not included in Pfam-A. Domainer uses pairwise homology data from BLASTP to construct aligned families.

Alternative protein family databases that may be used include PRINTS and BLOCKS. Both are based on a set of ungapped blocks of aligned residues. However, these programs typically contain short conserved regions whereas Pfam represents a library of complete domains that facilitates automated annotation. Comparisons of Pfam profiles may also be performed using genomic and EST (An abbreviation for expressed sequence tag which are defined as single pass sequencing of cDNAs usually 5') data with the programs, Genewise and ESTwise, respectively. Both of these programs allow for introns and frame shifting errors.

E. BLOCKS

The determination of sequence relationships between unknown sequences and those that have been categorized can be problematic because background noise increases with the number of sequences, especially at a low level of similarity detection. One recent approach to this problem has been tested that efficiently detects and confirms weak or distant relationships among protein sequences based on a database of blocks. The BLOCKS database provides multiple alignments of sequences and contains blocks or protein motifs found in known families of proteins.

Other programs such as PRINTS [The PRINTS database of protein fingerprints prepared under the supervision of Terri Attwood at the University of Manchester.Reference: Attwood T. K., Croning M. D. R., Flower D. R., Lewis A. P., Mabey J. E., Scordis P., Selley J. N. and Wright W.; Nucleic Acids Res. 28:225-227(2000)] and Prodom also provide alignments, however, the BLOCKS database differs in the manner in which the database was constructed. Construction of the BLOCKS database [S. Henikoff & J. G. Henikoff, "Protein family classification based on searching a database of blocks", Genomics 19:97-107 (1994). S. Henikoff, J. G. Henikoff, W. J. Alford & S. Pietrokovski, "Automated construction and graphical presentation of protein blocks from unaligned sequences", Gene-COMBIS, Gene 163 (1995) GC 17-26. S. Pietrokovski, "Searching Databases of Conserved Seqeuence Regions by Aligning Protein Multiple-Alignments", NAR 24:3836-3845 (1996)] proceeds as follows: one starts with a group of sequences that presumably have one or motifs in common, such as those from the PROSITE database [Hofmann K., Bucher P., Falquet L. and Bairoch A.; Nucleic Acids Res. 27:215-219(1999)]. The PROTOMAT program [S Henikoff & J G Henikoff, "Automated assembly of protein blocks for database searching", NAR (1991) 19:6565-6572) using the MOTIF algorithm (H O Smith, et al, "Finding sequence motifs in groups of functionally related proteins", PNAS (1990) 87:826-830] then uses a motif finding program to scan sequences for similarity looking for spaced triplets of amino acids. The located blocks are then entered into the MOTOMAT program [The first step (PROTOMAT) finds candidate alignments and the second step (MOTOMAT) extends the alignments, then sorts them in such a way that a best set is chosen] for block assembly. Weights are computed for all sequences. Following construction of a BLOCKS database one can use BLIMPS [this is a tool to search BLOCKS. I believe the reference is Henikoff S, Henikoff J G, Alford W J and Pietroskouski S (1995) Gene 163 GC17-26] to performs searches of the BLOCKS database. Detailed descriptions of the construction and use of a BLOCKS database can be found in the following references: Henikoff, S. and Henikoff, J. G., Genomics, 19:97 [1994]; Henikoff, J. G. and Henikoff, S., Meth. Enz., 266:88 [1996].

F. PRINTS

The PRINTS database of protein family fingerprints can be used in addition to BLOCKS and PROSITE. These databases are considered to be secondary databases because they diagnose the relationship between sequences that yield function information. Presently, however, it is not recommended that these databases be used alone. Rather, it is strongly suggested that these pattern databases be used in conjunction with each other so that a direct comparison of results can be made to analyze their robustness.

Generally, these programs utilize pattern recognition to discover motifs within protein sequences. However, PRINTS goes one step further, it takes into account not simply single motifs but several motifs simultaneously that might characterize a family signature. Other programs, such as PROSITE, rely on pattern recognition but are limited by the fact that query sequences must match them exactly. Thus, sequences that vary slightly will be missed. In contrast, the PRINTS database fingerprinting approach is capable of identifying distant relatives due to its reliance on the fact that sequences do not have to match the query exactly. Instead they are scored according to how well they fit each motif in the signature. Another advantage of PRINTS is that it allows the user to search both PRINTS and PROSITE simultaneously. A detailed description of the use of PRINTS can be found in the following reference: Attwood et al., Nucleic Acids Res. 25: 212 [1997].

III. Nucleic Acid Sequences, Including Related, Variant, Modified and Extended Sequences

This invention encompasses nucleic acids, polypeptides encoded by the nucleic acid sequences, and variants that retain at least one biological or other functional activity of the polynucleotide or polypeptide of interest. A preferred polynucleotide variant is one having at least 80%, and more preferably 90%, sequence identity to the sequence of interest. A most preferred polynucleotide variant is one having at least 95% sequence identity to the polynucleotide of interest.

In particularly preferred embodiments, the invention encompasses the polynucleotides comprising a polynucleotide encoded by SEQ ID NOs:1-7554. In particularly preferred embodiments, the nucleic acids are operably linked to an exogenous promoter (and in most preferred embodiments to a plant promoter) or present in a vector.

It will be appreciated by those skilled in the art that as a result of the degeneracy of the genetic code, a multitude of nucleotide sequences encoding a given polypeptide (for example, a polypeptide encoded by a nucleic acid of the present invention), some bearing minimal homology to the nucleotide sequences of any known and naturally occurring gene, may be produced. Thus, the invention contemplates each and every possible variation of nucleotide sequence that could be made by selecting combinations based on possible codon choices. These combinations are made in accordance with the standard triplet genetic code as applied to the nucleotide sequence of the naturally occurring polypeptide, and all such variations are to be considered as being specifically disclosed.

Although nucleotide sequences that encode a given polypeptide (for example, a polypeptide encoded by a nucleic acid of the present invention) and its variants are preferably capable of hybridizing to the nucleotide sequence of the naturally occurring polypeptide under appropriately selected conditions of stringency, it may be advantageous to produce nucleotide sequences encoding the polypeptide or its derivatives possessing a substantially different codon usage. Codons may be selected to increase the rate at which expression of the peptide occurs in a particular prokaryotic or eukaryotic host in, accordance with the frequency with which particular codons are utilized by the host. Other reasons for substantially altering the nucleotide sequence encoding a polypeptide and its derivatives without altering the encoded amino acid sequences include the production of RNA transcripts having more desirable properties, such as a greater half-life, than transcripts produced from the naturally occurring sequence.

The invention also encompasses production of DNA sequences, or portions thereof, that encode a polynucleotide and its variants, entirely by synthetic chemistry. After production, the synthetic sequence may be inserted into any of the many available expression vectors and cell systems using reagents that are well known in the art. Moreover, synthetic chemistry may be used to introduce mutations into a sequence encoding a polynucleotide of the present invention or any portion thereof.

Also encompassed by the invention are polynucleotide sequences that are capable of hybridizing to SEQ ID NOs:1-7554 under various conditions of stringency (for example, conditions ranging from low to high stringency). Hybridization conditions are based on the melting temperature T.sub.m of the nucleic acid binding complex or probe, as taught in Wahl and Berger, Methods Enzymol., 152:399 [1987] and Kimmel, Methods Enzymol., 152:507 [1987], and may be used at a defined stringency.

Modified nucleic acid sequences encoding a polynucleotide of the present invention include deletions, insertions, or substitutions of different nucleotides resulting in a polynucleotide that encodes the same polypeptide or a functionally equivalent polynucleotide or polypeptide. The encoded protein may also contain deletions, insertions, or substitutions of amino acid residues that produce a silent change and result in a functionally equivalent polypeptide. Deliberate amino acid substitutions may be made on the basis of similarity in polarity, charge, solubility, hydrophobicity, hydrophilicity, and/or the amphipathic nature of the residues as long as the biological activity of the polypeptide is retained. For example, negatively charged amino acids may include aspartic acid and glutamic acid; positively charged amino acids may include lysine and arginine; and amino acids with uncharged polar head groups having similar hydrophilicity values may include leucine, isoleucine, and valine; glycine and alanine; asparagine and glutamine; serine and threonine; phenylalanine and tyrosine.

Also included within the scope of the present invention are alleles of the genes encoding polypeptides. As used herein, an "allele" or "allelic sequence" is an alternative form of the gene that may result from at least one mutation in the nucleic acid sequence. Alleles may result in modified mRNAs or polypeptides whose structure or function may or may not be modified. Any given gene may have none, one, or many allelic forms. Common mutational changes that give rise to alleles are generally ascribed to natural deletions, additions, or substitutions of nucleotides. Each of these types of changes may occur alone, or in combination with the others, one or more times in a given sequence.

Methods for DNA sequencing that are well known and generally available in the art may be used to practice any embodiments of the invention. The methods may employ such enzymes as the Klenow fragment of DNA polymerase I, SEQUENASE (US Biochemical Corporation; Cleveland, Ohio), TAQ polymerase (U.S. Biochemical Corporation, Cleveland, Ohio), thermostable T7 polymerase (Amersham Pharmacia Biotech; Chicago, Ill.), or combinations of recombinant polymerases and proofreading exonucleases such as the ELONGASE amplification system (Life Technologies, Inc.; Rockville, Md.). Preferably, the process is automated with machines such as the MICROLAB 2200 (Hamilton Company; Reno, Nev.), PTC200 DNA Engine thermal cycler (MJ Research; Watertown, Mass.) and the ABI 377 DNA sequencer (Perkin Elmer).

The nucleic acid sequences encoding a polynucleotide of the present invention may be extended utilizing a partial nucleotide sequence and employing various methods known in the art to detect upstream sequences such as promoters and regulatory elements. For example, one method that may be employed, "restriction-site" PCR, uses universal primers to retrieve unknown sequence adjacent to a known locus [Sarkar, PCR Methods Applic. 2:318 (1993)]. In particular, genomic DNA is first amplified in the presence of primer to linker sequence and a primer specific to the known region. The amplified sequences are then subjected to a second round of PCR with the same linker primer and another specific primer internal to the first one. Products of each round of PCR are transcribed with an appropriate RNA polymerase and sequenced using reverse transcriptase.

Inverse PCR may also be used to amplify or extend sequences using divergent primers based on a known region [Triglia et al., Nucleic Acids Res. 16:8186 (1988)]. The primers may be designed using OLIGO 4.06 primer analysis software (National Biosciences Inc.; Plymouth, Minn.), or another appropriate program, to be 22-30 nucleotides in length, to have a GC content of 50% or more, and to anneal to the target sequence at temperatures about 68-72.degree. C. The method uses several restriction enzymes to generate a suitable fragment in the known region of a gene. The fragment is then circularized by intramolecular ligation and used as a PCR template.

Another method that may be used is capture PCR which involves PCR amplification of DNA fragments adjacent to a known sequence in human and yeast artificial chromosome DNA [Lagerstrom et al., PCR Methods Applic. 1:111 (1991)]. In this method, multiple restriction enzyme digestions and ligations may also be used to place an engineered double-stranded sequence into an unknown portion of the DNA molecule before performing PCR.

Another method that may be used to retrieve unknown sequences is that of Parker et al., Nucleic Acids Res., 19:3055 [1991]. Additionally, one may use PCR, nested primers, and PROMOTERFINDER DNA Walking Kits libraries (Clontech; Palo Alto, Calif.) to walk in genomic DNA. This process avoids the need to screen libraries and is useful in finding intron/exon junctions.

When screening for full-length cDNAs, it is preferable to use libraries that have been size-selected to include larger cDNAs. Also, random-primed libraries are preferable, in that they will contain more sequences that contain the 5' regions of genes. Use of a randomly primed library may be especially preferable for situations in which an oligo d(T) library does not yield a full-length cDNA. Genomic libraries may be useful for extension of sequence into the 5' and 3' non-transcribed regulatory regions.

Capillary electrophoresis systems that are commercially available (for example, from PE Biosystems, Inc.; Foster City, Calif.) may be used to analyze the size or confirm the nucleotide sequence of sequencing or PCR products. In particular, capillary sequencing may employ flowable polymers for electrophoretic separation, four different fluorescent dyes (one for each nucleotide) that are laser activated, and detection of the emitted wavelengths by a charge coupled device camera. Output/light intensity may be converted to electrical signal using appropriate software (for example, GENOTYPER and SEQUENCE NAVIGATOR from PE Biosystems; Foster City, Calif.) and the entire process from loading of samples to computer analysis and electronic data display may be computer controlled. Capillary electrophoresis is especially preferable for the sequencing of small pieces of DNA that might be present in limited amounts in a particular sample.

It is contemplated that the nucleic acids disclosed herein can be utilized as starting nucleic acids for directed evolution. In some embodiments, artificial-evolution is performed by random mutagenesis (for example, by utilizing error-prone PCR to introduce random mutations into a given coding sequence). This method requires that the frequency of mutation be finely tuned. As a general rule, beneficial mutations are rare, while deleterious mutations are common. This is because the combination of a deleterious mutation and a beneficial mutation often results in an inactive enzyme. The ideal number of base substitutions for a targeted gene is usually between 1.5 and 5 [Moore and Arnold, Nat. Biotech., 14, 458-67 (1996); Leung et al., Technique, 1:11-15 (1989); Eckert and Kunkel, PCR Methods Appl., 1:17-24 (1991); Caldwell and Joyce, PCR Methods Appl., 2:28-33 (1992); and Zhao and Arnold, Nuc. Acids. Res., 25:1307-08 (1997)]. After mutagenesis, the resulting clones are selected for desirable activity. Successive rounds of mutagenesis and selection are often necessary to develop enzymes with desirable properties. It should be noted that only the useful mutations are carried over to the next round of mutagenesis.

In other embodiments of the present invention, the polynucleotides of the present invention are used in gene shuffling or sexual PCR procedures (for example, Smith, Nature, 370:324-25 [1994]; U.S. Pat. Nos. 5,837,458; 5,830,721; 5,811,238; and 5,733,731, each of which is herein incorporated by reference). Gene shuffling involves random fragmentation of several mutant DNAs followed by their reassembly by PCR into full length molecules. Examples of various gene shuffling procedures include, but are not limited to, assembly following DNase treatment, the staggered extension process (STEP), and random priming in vitro recombination. In the DNase mediated method, DNA segments isolated from a pool of positive mutants are cleaved into random fragments with DNaseI and subjected to multiple rounds of PCR with no added primer. The lengths of random fragments approach that of the uncleaved segment as the PCR cycles proceed, resulting in mutations in present in different clones becoming mixed and accumulating in some of the resulting sequences. Multiple cycles of selection and shuffling have led to the functional enhancement of several enzymes [Stemmer, Nature, 370:398-91 (1994); Stemmer, Proc. Natl. Acad. Sci. USA, 91, 10747-51 (1994); Crameri et al., Nat. Biotech., 14:315-19 (1996); Zhang et al., Proc. Natl. Acad. Sci. USA, 94:4504-09 (1997); and Crameri et al., Nat. Biotech., 15:436-38 (1997)].

IV. Vectors, Engineering, and Expression of Sequences

In another embodiment of the invention, the polynucleotide sequences of the present invention and fragments and portions thereof, may be used in recombinant DNA molecules to direct expression of an mRNA or polypeptide in appropriate host cells. Due to the inherent degeneracy of the genetic code, other DNA sequences that encode substantially the same or a functionally equivalent mRNA or amino acid sequence may be produced and these sequences may be used to clone and express polypeptides (for example, a polypeptide encoded by a nucleic acid of the present invention).

As will be understood by those of skill in the art, it may be advantageous to produce nucleotide sequences possessing non-naturally occurring codons. For example, codons preferred by a particular prokaryotic or eukaryotic host can be selected to increase the rate of protein expression or to produce a recombinant RNA transcript having desirable properties, such as a half-life that is longer than that of a transcript generated from the naturally occurring sequence.

The nucleotide sequences of the present invention can be engineered using methods generally known in the art in order to alter the polypeptide sequences for a variety of reasons, including but not limited to, alterations that modify the cloning, processing, and/or expression of the gene product. DNA shuffling by random fragmentation and PCR reassembly of gene fragments and synthetic oligonucleotides may be used to engineer the nucleotide sequences. For example, site-directed mutagenesis may be used to insert new restriction sites, alter glycosylation patterns, change codon preference, produce splice variants, or introduce mutations, and so forth.

In another embodiment of the invention, natural, modified, or recombinant nucleic acid sequences encoding a polypeptide may be ligated to a heterologous sequence to encode a fusion protein. For example, to screen peptide libraries for inhibitors of the polypeptides activity (for example, enzymatic activity), it may be useful to encode a chimeric protein that can be recognized by a commercially available antibody. A fusion protein may also be engineered to contain a cleavage site located between the polypeptide encoding sequence and the heterologous protein sequence, so that the polypeptide of interest may be cleaved and purified away from the heterologous moiety.

In another embodiment, sequences encoding a polypeptide (for example, a polypeptide encoded by a nucleic acid of the present invention) may be synthesized, in whole or in part, using chemical methods well known in the art [See for example, Caruthers et al., Nucl. Acids Res. Symp. Ser. 215 (1980); Horn et al., Nucl. Acids Res. Symp. Ser. 225 (1980)]. Alternatively, the protein itself may be produced using chemical methods to synthesize the amino acid sequence of the polypeptide of interest (for example, a polypeptide encoded by a nucleic acid of the present invention), or a portion thereof. For example, peptide synthesis can be performed using various solid-phase techniques [Roberge et al., Science 269:202 (1995)] and automated synthesis may be achieved, for example, using the ABI 431A peptide synthesizer (PE Corporation, Norwalk, Conn.).

The newly synthesized peptide may be substantially purified by preparative high performance liquid chromatography [See for example, Creighton, T. (1983) Proteins, Structures and Molecular Principles, WH Freeman and Co., New York, N.Y.]. The composition of the synthetic peptides may be confirmed by amino acid analysis or sequencing (for example, the Edman degradation procedure; or Creighton, supra). Additionally, the amino acid sequence of the polypeptide of interest or any part thereof, may be changed during direct synthesis and/or combined using chemical methods with sequences from other proteins, or any part thereof, to produce a variant polypeptide.

In order to express a biologically active polypeptide (for example, a polypeptide encoded by a nucleic acid of the present invention) or RNA, the nucleotide sequences encoding the polypeptide or functional equivalents, may be inserted into appropriate expression vector, that is, a vector that contains the necessary elements for the transcription and translation of the inserted coding sequence.

Methods that are well known to those skilled in the art may be used to construct expression vectors containing sequences encoding polypeptides (for example, a polypeptide encoded by a nucleic acid of the present invention) and appropriate transcriptional and translational control elements. These methods include in vitro recombinant DNA techniques, synthetic techniques, and in vivo genetic recombination. Such techniques are described in Sambrook. et al. (1989) Molecular Cloning, A Laboratory Manual, Cold Spring Harbor Press, Plainview, N.Y., and Ausubel, F. M. et al. (1989) Current Protocols in Molecular Biology, John Wiley & Sons, New York, N.Y.

A variety of expression vector/host systems may be utilized to contain and express sequences encoding a polypeptide of interest. These include, but are not limited to, microorganisms such as bacteria transformed with recombinant bacteriophage, plasmid, or cosmid DNA expression vectors; yeast transformed with yeast expression vectors; insect cell systems infected with virus expression vectors (for example, baculovirus); plant cell systems transformed with virus expression vectors (for example, cauliflower mosaic virus, CaMV; tobacco mosaic virus, TMV; brome mosaic virus) or with bacterial expression vectors (for example, Ti or pBR322 plasmids); or animal cell systems.

The "control elements" or "regulatory sequences" are those non-translated regions of the vector (for example, enhancers, promoters, 5' and 3' untranslated regions) that interact with host cellular proteins to carry out transcription and translation. Such elements may vary in their strength and specificity. Depending on the vector system and host utilized, any number of suitable transcription and translation elements, including constitutive and inducible promoters, may be used. For example, when cloning in bacterial systems, inducible promoters such as the hybrid lacZ promoter of the BLUESCRIPT phagemid (Stratagene; LaJolla, Calif.) or PSPORT1 plasmid (Life Technologies, Inc.; Rockville, Md.) and the like may be used. The baculovirus polyhedrin promoter may be used in insect cells. Promoters or enhancers derived from the genomes of plant cells (for example, heat shock, RUBISCO; and storage protein genes) or from plant viruses (for example, viral promoters or leader sequences) may be cloned into the vector. In mammalian cell systems, promoters from mammalian genes or from mammalian viruses are preferable. If it is necessary to generate a cell line that contains multiple copies of the sequence encoding a polypeptide, vectors based on SV40 or EBV may be used with an appropriate selectable marker.

In bacterial systems, a number of expression vectors may be selected depending upon the use intended for the polypeptide of interest. For example, when large quantities of the polypeptide are needed for the induction of antibodies, vectors that direct high level expression of fusion proteins that are readily purified may be used. Such vectors include, but are not limited to, the multifunctional E. coli cloning and expression vectors such as BLUESCRIPT phagemid (Stratagene; La Jolla, Calif.), in which the sequence encoding the polypeptide of interest may be ligated into the vector in frame with sequences for the amino-terminal Met and the subsequent 7 residues of beta-galactosidase so that a hybrid protein is produced; pIN vectors (Van Heeke and Schuster, J. Biol. Chem. 264:5503 [1989]; and the like. pGEMX vectors (Promega Corporation; Madison, Wis.) may also be used to express foreign polypeptides as fusion proteins with glutathione S-transferase (GST). In general, such fusion proteins are soluble and can easily be purified from lysed cells by adsorption to glutathione-agarose beads followed by elution in the presence of free glutathione. Proteins made in such systems may be designed to include heparin, thrombin, or factor XA protease cleavage sites so that the cloned polypeptide of interest can be released from the GST moiety at will.

In the yeast Saccharomyces cerevisiae, a number of vectors containing constitutive or inducible promoters such as alpha factor, alcohol oxidase, and PGH may be used. For reviews, See for example, Ausubel et al. (supra) and Grant et al., Methods Enzymol. 153:516 [1987].

In cases where plant expression vectors are used, the expression of sequences encoding polypeptides may be driven by any of a number of promoters. In a preferred embodiment, plant vectors are created using a recombinant plant virus containing a recombinant plant viral nucleic acid, as described in PCT publication WO 96/40867. Subsequently, the recombinant plant viral nucleic acid that contains one or more non-native nucleic acid sequences may be transcribed or expressed in the infected tissues of the plant host and the product of the coding sequences may be recovered from the plant, as described in WO 99/36516.

An important feature of this embodiment is the use of recombinant plant viral nucleic acids that contain one or more non-native subgenomic promoters capable of transcribing or expressing adjacent nucleic acid sequences in the plant host and that result in replication and local and/or systemic spread in a compatible plant host. The recombinant plant viral nucleic acids have substantial sequence homology to plant viral nucleotide sequences and may be derived from an RNA, DNA, cDNA or a chemically synthesized RNA or DNA. A partial listing of suitable viruses is described below.

The first step in producing recombinant plant viral nucleic acids according to this particular embodiment is to modify the nucleotide sequences of the plant viral nucleotide sequence by known conventional techniques such that one or more non-native subgenomic promoters are inserted into the plant viral nucleic acid without destroying the biological, function of the plant viral nucleic acid. The native coat protein coding sequence may be deleted in some embodiments, placed under the control of a non-native subgenomic promoter in other embodiments, or retained in a further embodiment. If it is deleted or otherwise inactivated, a non-native coat protein gene is inserted under control of one of the non-native subgenomic promoters, or optionally under control of the native coat protein gene subgenomic promoter. The non-native coat protein is capable of encapsidating the recombinant plant viral nucleic acid to produce a recombinant plant virus. Thus, the recombinant plant viral nucleic acid contains a coat protein coding sequence, that may be native or a nonnative coat protein coding sequence, under control of one of the native or non-native subgenomic promoters. The coat protein is involved in the systemic infection of the plant host.

Some of the viruses that meet this requirement include viruses from the tobamovirus group such as Tobacco Mosaic virus (TMV), Ribgrass Mosaic Virus (RGM), Cowpea Mosaic virus (CMV), Alfalfa Mosaic virus (AMV), Cucumber Green Mottle Mosaic virus Watermelon Strain (CGMMV-W) and Oat Mosaic virus (OMV) and viruses from the brome mosaic virus group such as Brome Mosaic virus (BMV), Broad Bean Mottle virus and Cowpea Chlorotic mottle virus. Additional suitable viruses include Rice Necrosis virus (RNV), and geminiviruses such as Tomato Golden Mosaic virus (TGMV), Cassava Latent virus (CLV) and Maize Streak virus (MSV). However, the invention should not be construed as limited to using these particular viruses, but rather the method of the present invention is contemplated to include all plant viruses at a minimum.

Other embodiments of plant vectors used for the expression of sequences encoding polypeptides include, for example, viral promoters such as the 35S and 19S promoters of CaMV used alone or in combination with the omega leader sequence from TMV [Takamatsu, EMBO J. 6:307 (1987)]. Alternatively, plant promoters such as the small subunit of RUBISCO or heat shock promoters may be used [Coruzzi et al., EMBO J. 3:1671 (1984); Broglie et al., Science 224:838 (1984); and Winter et al., Results Probl. Cell Differ. 17:85 (1991)]. These constructs can be introduced into plant cells by direct DNA transformation or pathogen-mediated transfection. Such techniques are described in a number of generally available reviews (See for example, Hobbs, S. or Murry, L. E. in McGraw Hill Yearbook of Science and Technology (1992) McGraw Hill, New York, N.Y.; pp. 191-196.

The present invention further provides transgenic plants comprising the polynucleotides of the present invention. In some preferred embodiments, Agrobacterium mediated transfection is utilized to create transgenic plants. Since most dicotyledonous plant are natural hosts for Agrobacterium, almost every dicotyledonous plant may be transformed by Agrobacterium in vitro. Although monocotyledonous plants, and in particular, cereals and grasses, are not natural hosts to Agrobacterium, work to transform them using Agrobacterium has also been carried out [Hooykas-Van Slogteren et al. (1984) Nature 311:763-764]. Plant genera that may be transformed by Agrobacterium include Arabidopsis, Chrysanthemum, Dianthus, Gerbera, Euphorbia, Pelaronium, Ipomoea, Passiflora, Cyclamen, Malus, Prunus, Rosa, Rubus, Populus, Santalum, Allium, Lilium, Narcissus, Ananas, Arachis, Phaseolus and Pisum.

For transformation with Agrobacterium, disarmed Agrobacterium cells are transformed with recombinant Ti plasmids of Agrobacterium tumefaciens or Ri plasmids of Agrobacterium rhizogenes (such as those described in U.S. Pat. No. 4,940,838, the entire contents of which are herein incorporated by reference). The nucleic acid sequence of interest is then stably integrated into the plant genome by infection with the transformed Agrobacterium strain. For example, heterologous nucleic acid sequences have been introduced into plant tissues using the natural DNA transfer system of Agrobacterium tumefaciens and Agrobacterium rhizogenes bacteria [for review, see Klee et al. (1987) Ann. Rev. Plant Phys. 38:467-486].

There are three common methods to transform plant cells with Agrobacterium. The first method is co-cultivation of Agrobacterium with cultured isolated protoplasts. This method requires an established culture system that allows culturing protoplasts and plant regeneration from cultured protoplasts. The second method is transformation of cells or tissues with Agrobacterium. This method requires (a) that the plant cells or tissues can be transformed by Agrobacterium and (b) that the transformed cells or tissues can be induced to regenerate into whole plants. The third method is transformation of seeds, apices or meristems with Agrobacterium. This method requires micropropagation.

The efficiency of transformation by Agrobacterium may be enhanced by using a number of methods known in the art. For example, the inclusion of a natural wound response molecule such as acetosyringone (AS) to the Agrobacterium culture has been shown to enhance transformation efficiency with Agrobacterium tumefaciens [Shahla et al., (1987) Plant Molec. Biol. 8:291-298]. Alternatively, transformation efficiency may be enhanced by wounding the target tissue to be transformed. Wounding of plant tissue may be achieved, for example, by punching, maceration, bombardment with microprojectiles, etc. [See e.g., Bidney et al., (1992) Plant Molec. Biol. 18:301-313].

In still further embodiments, the plant cells are transfected with vectors via particle bombardment (i.e., with a gene gun). Particle mediated gene transfer methods are known in the art, are commercially available, and include, but are not limited to, the gas driven gene delivery instrument descried in McCabe, U.S. Pat. No. 5,584,807, the entire contents of which are herein incorporated by reference. This method involves coating the nucleic acid sequence of interest onto heavy metal particles, and accelerating the coated particles under the pressure of compressed gas for delivery to the target tissue.

Other particle bombardment methods are also available for the introduction of heterologous nucleic acid sequences into plant cells. Generally, these methods involve depositing the nucleic acid sequence of interest upon the surface of small, dense particles of a material such as gold, platinum, or tungsten. The coated particles are themselves then coated onto either a rigid surface, such as a metal plate, or onto a carrier sheet made of a fragile material such as mylar. The coated sheet is then accelerated toward the target biological tissue. The use of the flat sheet generates a uniform spread of accelerated particles that maximizes the number of cells receiving particles under uniform conditions, resulting in the introduction of the nucleic acid sample into the target tissue.

An insect system may also be used to express polypeptides (for example, a polypeptide encoded by a nucleic acid of the present invention). For example, in one such system, Autographa californica nuclear polyhedrosis virus (AcNPV) is used as a vector to express foreign genes in Spodoptera frugiperda cells or in Trichoplusia larvae. The sequences encoding a polypeptide of interest may be cloned into a non-essential region of the virus, such as the polyhedrin gene, and placed under control of the polyhedrin promoter. Successful insertion of the nucleic acid sequence encoding the polypeptide of interest will render the polyhedrin gene inactive and produce recombinant virus lacking coat protein. The recombinant viruses may then be used to infect, for example, S. frugiperda cells or Trichoplusia larvae in which the polypeptide may be expressed [Engelhard et al., Proc. Nat. Acad. Sci. 91:3224 (1994)].

In mammalian host cells, a number of viral-based expression systems may be utilized. In cases where an adenovirus is used as an expression vector, sequences encoding polypeptides may be ligated into an adenovirus transcription/translation complex consisting of the late promoter and tripartite leader sequence. Insertion in a non-essential E1 or E3 region of the viral genome may be used to obtain a viable virus that is capable of expressing the polypeptide in infected host cells [Logan and Shenk, Proc. Natl. Acad. Sci., 81:3655 (1984)]. In addition, transcription enhancers, such as the Rous sarcoma virus (RSV) enhancer, may be used to increase expression in mammalian host cells.

Specific initiation signals may also be used to achieve more efficient translation of sequences encoding the polypeptide of interest. Such signals include the ATG initiation codon and adjacent sequences. In cases where sequences encoding the polypeptide of interest, its initiation codon, and upstream sequences are inserted into the appropriate expression vector, no additional transcriptional or translational control signals may be needed. However, in cases where only coding sequence, or a portion thereof, is inserted, exogenous translational control signals including the ATG initiation codon should be provided. Furthermore, the initiation codon should be in the correct reading frame to ensure translation of the entire insert. Exogenous translational elements and initiation codons may be of various origins, both natural and synthetic. The efficiency of expression may be enhanced by the inclusion of enhancers that are appropriate for the particular cell system that is used, such as those described in the literature [Scharf et al., Results Probl. Cell Differ., 20:125 (1994)].

In addition, a host cell strain may be chosen for its ability to modulate the expression of the inserted sequences or to process the expressed protein in the desired fashion. Such modifications of the polypeptide include, but are not limited to, acetylation, carboxylation, glycosylation, phosphorylation, lipidation, and acylation. Post-translational processing that cleaves a "prepro" form of the protein may also be used to facilitate correct insertion, folding and/or function. Different host cells such as CHO, HeLa, MDCK, HEK293, and WI38, that have specific cellular machinery and characteristic mechanisms for such post-translational activities, may be chosen to ensure the correct modification and processing of the foreign protein.

For long-term, high-yield production of recombinant proteins, stable expression is preferred. For example, cell lines that stably express the polypeptide of interest (for example, a polypeptide encoded by a nucleic acid of the present invention) may be transformed using expression vectors that may contain viral origins of replication and/or endogenous expression elements and a selectable marker gene on the same or on a separate vector. Following the introduction of the vector, cells may be allowed to grow for 1-2 days in an enriched media before they are switched to selective media. The purpose of the selectable marker is to confer resistance to selection, and its presence allows growth and recovery of cells that successfully express the introduced sequences. Resistant clones of stably transformed cells may be proliferated using tissue culture techniques appropriate to the cell type.

Any number of selection systems may be used to recover transformed cell lines. These include, but are not limited to, the herpes simplex virus thymidine kinase [Wigler et al., Cell 11:223 (1977)] and adenine phosphoribosyltransferase [Lowy et al., Cell 22:817 (1980)] genes that can be employed in tk.sup.- or aprt.sup.- cells, respectively. Also, antimetabolite, antibiotic, or herbicide resistance can be used as the basis for selection; for example, dhfr, which confers resistance to methotrexate [Wigler et al., Proc. Natl. Acad. Sci., 77:3567 (1980)]; npt, which confers resistance to the aminoglycosides neomycin and G-418 [Colbere-Garapin et al., J. Mol. Biol., 150:1 (1981)]; and als or pat, which confer resistance to chlorsulfuron and phosphinotricin acetyltransferase, respectively (Murry, supra). Additional selectable genes have been described, for example, trpB, which allows cells to utilize indole in place of tryptophan, or hisD, which allows cells to utilize histinol in place of histidine [Hartman and Mulligan, Proc. Natl. Acad. Sci., 85:8047 (1988)]. Recently, the use of visible markers has gained popularity with such markers as anthocyanins, .alpha.-glucuronidase and its substrate GUS, and luciferase and its substrate luciferin, being widely used not only to identify transformants, but also to quantify the amount of transient or stable protein expression attributable to a specific vector system [Rhodes et al., Methods Mol. Biol., 55:121 (1995)].

Although the presence/absence of marker gene expression suggests that the gene of interest is also present, its presence and expression may need to be confirmed. For example, if the sequence encoding a polypeptide is inserted within a marker gene sequence, recombinant cells containing sequences encoding the polypeptide can be identified by the absence of marker gene function. Alternatively, a marker gene can be placed in tandem with a sequence encoding the polypeptide under the control of a single promoter. Expression of the marker gene in response to induction or selection usually indicates expression of the tandem gene as well.

Alternatively, host cells that contain the nucleic acid sequence encoding the polypeptide of interest (for example, a polypeptide encoded by a nucleic acid of the present invention) and express the polypeptide may be identified by a variety of procedures known to those of skill in the art. These procedures include, but are not limited to, DNA-DNA or DNA-RNA hybridizations and protein bioassay or immunoassay techniques that include membrane, solution, or chip based technologies for the detection and/or quantification of nucleic acid or protein.

The presence of polynucleotide sequences encoding a polypeptide of interest (for example, a polypeptide encoded by a nucleic acid of the present invention) can be detected by DNA-DNA or DNA-RNA hybridization or amplification using probes or portions or fragments of polynucleotides encoding the polypeptide. Nucleic acid amplification based assays involve the use of oligonucleotides or oligomers based on the sequences encoding the polypeptide to detect transformants containing DNA or RNA encoding the polypeptide. As used herein "oligonucleotides" or "oligomers" refer to a nucleic acid sequence of at least about 10 nucleotides and as many as about 60 nucleotides, preferably about 15 to 30 nucleotides, and more preferably about 20-25 nucleotides, that can be used as a probe or amplimer.

A variety of protocols for detecting and measuring the expression of a polypeptide (for example, a polypeptide encoded by a nucleic acid of the present invention), using either polyclonal or monoclonal antibodies specific for the protein are known in the art. Examples include enzyme-linked immunosorbent assay (ELISA), radioimmunoassay (RIA), and fluorescence activated cell sorting (FACS). A two-site, monoclonal-based immunoassay utilizing monoclonal antibodies reactive to two non-interfering epitopes on the polypeptide is preferred, but a competitive binding assay may be employed. These and other assays are described, among other places, in Hampton et al., 1990; Serological Methods, a Laboratory Manual, APS Press, St Paul, Minn. and Maddox et al., J. Exp. Med., 158:1211 (1983).

A wide variety of labels and conjugation techniques are known by those skilled in the art and may be used in various nucleic acid and amino acid assays. Means for producing labeled hybridization or PCR probes for detecting sequences related to polynucleotides encoding a polypeptide of interest include oligonucleotide labeling, nick translation, end-labeling or PCR amplification using a labeled nucleotide. Alternatively, the sequences encoding the polypeptide, or any portions thereof may be cloned into a vector for the production of an mRNA probe. Such vectors are known in the art, are commercially available, and may be used to synthesize RNA probes in vitro by addition of an appropriate RNA polymerase such as T7, T3, or SP6 and labeled nucleotides. These procedures may be conducted using a variety of commercially available kits from Pharmacia & Upjohn (Kalamazoo, Mich.), Promega Corporation (Madison, Wis.) and U.S. Biochemical Corp. (Cleveland, Ohio). Suitable reporter molecules or labels, that may be used, include radionuclides, enzymes, fluorescent, chemiluminescent, or chromogenic agents as well as substrates, cofactors, inhibitors, magnetic particles, and the like.

Host cells transformed with nucleotide sequences encoding a polypeptide of interest may be cultured under conditions suitable for the expression and recovery of the protein from cell culture. The protein produced by a recombinant cell may be secreted or contained intracellularly depending on the sequence and/or the vector used. As will be understood by those of skill in the art, expression vectors containing polynucleotides that encode the polypeptide of interest (for example, a polypeptide encoded by a nucleic acid of the present invention) may be designed to contain signal sequences that direct secretion of the polypeptide through a prokaryotic or eukaryotic cell membrane. Other recombinant constructions may be used to join sequences encoding the polypeptide to nucleotide sequence encoding a polypeptide domain that will facilitate purification of soluble proteins. Such purification facilitating domains include, but are not limited to, metal chelating peptides such as histidine-tryptophan modules that allow purification on immobilized metals, protein A domains that allow purification on immobilized immunoglobulin, and the domain utilized in the FLAGS extension/affinity purification system (Immunex Corp., Seattle, Wash.). The inclusion of cleavable linker sequences such as those specific for Factor XA or enterokinase (available from Invitrogen; San Diego, Calif.) between the purification domain and the polypeptide of interest may be used to facilitate purification. One such expression vector provides for expression of a fusion protein containing the polypeptide of interest and a nucleic acid encoding 6 histidine residues preceding a thioredoxin or an enterokinase cleavage site. The histidine residues facilitate purification on IMIAC (immobilized metal ion affinity chromatography) as described in Porath et al., Prot. Exp. Purif., 3:263 [1992] while the enterokinase cleavage site provides a means for purifying the polypeptide from the fusion protein. A discussion of vectors that contain fusion proteins is provided in Kroll et al., DNA Cell Biol., 12:441 (1993).

In addition to recombinant production, fragments of the polypeptide of interest may be produced by direct peptide synthesis using solid-phase techniques [Merrifield, J. Am. Chem. Soc., 85:2149 (1963)]. Protein synthesis may be performed using manual techniques or by automation. Automated synthesis may be achieved, for example, using the Applied Biosystems 431A peptide synthesizer (Perkin Elmer). Various fragments of the polypeptide may be chemically synthesized separately and combined using chemical methods to produce the full-length molecule.

V. Alteration of Gene Expression

It is contemplated that the polynucleotides of the present invention (for example, SEQ ID NOs:1-7554) may be utilized to either increase or decrease the level of corresponding mRNA and/or protein in transfected cells as compared to the levels in wild-type cells. Accordingly, in some embodiments, expression in plants by the methods described above leads to the overexpression of the polypeptide of interest in transgenic plants, plant tissues, or plant cells. The present invention is not limited to any particular mechanism. Indeed, an understanding of a mechanism is not required to practice the present invention. However, it is contemplated that overexpression of the polynucleotides of the present invention will alter the expression of the gene comprising the nucleic acid sequence of the present invention. In some embodiments, more than one of SEQ ID NOs:1-7554 are expressed in a given plant. The sequences may be contained in the same vector or in different vectors. The sequences can influence the same metabolic trait (e.g., fatty acid metabolism or one of the other traits discussed in more detail below) or multiple metabolic traits (e.g., fatty acid and carbohydrate metabolism).

In other embodiments of the present invention, the polynucleotides are utilized to decrease the level of the protein or mRNA of interest in transgenic plants, plant tissues, or plant cells as compared to wild-type plants, plant tissues, or plant cells. One method of reducing protein expression utilizes expression of antisense transcripts (for example, U.S. Pat. Nos. 6,031,154; 5,453,566; 5,451,514; 5,859,342; and 4,801,340, each of which is incorporated herein by reference). Antisense RNA has been used to inhibit plant target genes in a tissue-specific manner [for example, Van der Krol et al., Biotechniques 6:958-976 (1988)]. Antisense inhibition has been shown using the entire cDNA sequence as well as a partial cDNA sequence [for example, Sheehy et al., Proc. Natl. Acad. Sci. USA 85:8805-8809 (1988); Cannon et al., Plant Mol. Biol. 15:39-47 (1990)]. There is also evidence that 3' non-coding sequence fragment and 5' coding sequence fragments, containing as few as 41 base-pairs of a 1.87 kb cDNA, can play important roles in antisense inhibition [Ch'ng et al., Proc. Natl. Acad. Sci. USA 86:10006-10010 (1989)].

Accordingly, in some embodiments, the nucleic acids of the present invention (for example, SEQ ID NOs:1-7554, and fragments and variants thereof) are oriented in a vector and expressed so as to produce antisense transcripts. To accomplish this, a nucleic acid segment from the desired gene is cloned and operably linked to a promoter such that the antisense strand of RNA will be transcribed. The expression cassette is then transformed into plants and the antisense strand of RNA is produced. The nucleic acid segment to be introduced generally will be substantially identical to at least a portion of the endogenous gene or genes to be repressed. The sequence, however, need not be perfectly identical to inhibit expression. The vectors of the present invention can be designed such that the inhibitory effect applies to other proteins within a family of genes exhibiting homology or substantial homology to the target gene.

Furthermore, for antisense suppression, the introduced sequence also need not be full length relative to either the primary transcription product or fully processed mRNA. Generally, higher homology can be used to compensate for the use of a shorter sequence. Furthermore, the introduced sequence need not have the same intron or exon pattern, and homology of non-coding segments may be equally effective. Normally, a sequence of between about 30 or 40 nucleotides and up to about the full length of the coding region should be used, although a sequence of at least about 100 nucleotides is preferred, a sequence of at least about 200 nucleotides is more preferred, and a sequence of at least about 500 nucleotides is especially preferred.

Catalytic RNA molecules or ribozymes can also be used to inhibit expression of the target gene or genes. It is possible to design ribozymes that specifically pair with virtually any target RNA and cleave the phosphodiester backbone at a specific location, thereby functionally inactivating the target RNA. In carrying out this cleavage, the ribozyme is not itself changed, and is thus capable of recycling and cleaving other molecules, making it a true enzyme. The inclusion of ribozyme sequences within antisense RNAs confers RNA-cleaving activity upon them, thereby increasing the activity of the constructs.

A number of classes of ribozymes have been identified. One class of ribozymes is derived from a number of small circular RNAs that are capable of self-cleavage and replication in plants. The RNAs replicate either alone (viroid RNAs) or with a helper virus (satellite RNAs). Examples include RNAs from avocado sunblotch viroid and the satellite RNAs from tobacco ringspot virus, lucerne transient streak virus, velvet tobacco mottle virus, Solanum nodiflorum mottle virus and subterranean clover mottle virus. The design and use of target RNA-specific ribozymes is described in Haseloff, et al., Nature 334:585-591 (1988).

Another method of reducing protein expression utilizes the phenomenon of cosuppression or gene silencing (for example, U.S. Pat. Nos. 6,063,947; 5,686,649; and 5,283,184; each of which is incorporated herein by reference). The phenomenon of cosuppression has also been used to inhibit plant target genes in a tissue-specific manner. Cosuppression of an endogenous gene using a full-length cDNA sequence as well as a partial cDNA sequence (730 bp of a 1770 bp cDNA) are known [for example, Napoli et al., Plant Cell 2:279-289 [1990]; van der Krol et al., Plant Cell 2:291-299 (1990); Smith et al., Mol. Gen. Genetics 224:477-481 (1990]). Accordingly, in some embodiments the nucleic acids (for example, SEQ ID NOs: 1-7554, and fragments and variants thereof) from one species of plant are expressed in another species of plant to effect cosuppression of a homologous gene. Generally, where inhibition of expression is desired, some transcription of the introduced sequence occurs. The effect may occur where the introduced sequence contains no coding sequence per se, but only intron or untranslated sequences homologous to sequences present in the primary transcript of the endogenous sequence. The introduced sequence generally will be substantially identical to the endogenous sequence intended to be repressed. This minimal identity will typically be greater than about 65%, but a higher identity might exert a more effective repression of expression of the endogenous sequences. Substantially greater identity of more than about 80% is preferred, though about 95% to absolute identity would be most preferred. As with antisense regulation, the effect should apply to any other proteins within a similar family of genes exhibiting homology or substantial homology.

For cosuppression, the introduced sequence in the expression cassette, needing less than absolute identity, also need not be full length, relative to either the primary transcription product or fully processed mRNA. This may be preferred to avoid concurrent production of some plants that are overexpressers. A higher identity in a shorter than full length sequence compensates for a longer, less identical sequence. Furthermore, the introduced sequence need not have the same intron or exon pattern, and identity of non-coding segments will be equally effective. Normally, a sequence of the size ranges noted above for antisense regulation is used.

VI. Expression of Sequences Producing Altered Metabolic Characteristics

The present invention provides nucleic sequences involved in providing altered metabolic characteristics in plants. Plants transformed with viral vectors comprising the nucleic acid sequences of the present invention were screened for an altered metabolic characteristic. The results are presented in FIG. 6. Accordingly, in some embodiments, the present invention provides nucleic acid sequences that produce an altered metabolic characteristic when expressed in a plant (SEQ ID NOs:1-1165, 3703-4153 and 7389-7554, FIG. 1). The present invention is not limited to the particular nucleic acid sequences listed. Indeed, the present invention encompasses nucleic acid sequences (including sequences of the same, shorter, and longer lengths) that hybridize to the listed nucleic sequences under conditions ranging from low to high stringency and that also cause the altered metabolic characteristic. These sequences are conveniently identified by insertion into GENEWARE.TM. vectors and expression in plants as detailed in the examples.

The present invention is not limited to any particular mechanism. Indeed, an understanding of a mechanism is not required to practice the present invention. However, it is contemplated that the expression of genes comprising the nucleic acid sequences of the present invention effect biochemical pathways that lead to the alteration of metabolic characteristics of the present invention. For example, the expression of genetic function that effects the acyl acetylglycerol pathway that leads to the production of the altered esters such as 2-steroyl-1-acetylglycerol and 2-eicosanoyl-1-acetylglycerol. The present invention is not limited to alterations of any particular metabolic pathway. Indeed, the alteration of a variety of metabolic pathways is contemplated, including, but not limited to the pathways involved in the production of the following compounds: acids, fatty acids, branched fatty acids, alcohols, alkaloids, aleknes, alkynes, amino acids, carbohydrates, esters, glycerol, phenols, sterols, terpenes, isoprenoids, ketones, and quinones.

In some embodiments, the sequences are operably linked to a plant promoter or provided in a vector as described in more detail above. This present invention also contemplates plants transformed or transfected with these sequences, as well as seeds, fruit, leaves, stems and roots from such transfected plants. Furthermore, the sequences can be expressed in either sense or antisense orientation. In particularly preferred embodiments, the sequences are at least 30 nucleotides in length up to the length of the full-length of the corresponding gene. It is contemplated that sequences of less than full length (for example, greater than about 30 nucleotides) are useful for down regulation of gene expression via antisense or cosuppression. Suitable sequences are selected by chemically synthesizing the sequences, cloning into GENEWARE.TM. expression vectors, expressing in plants, and selecting plants with an altered metabolic characteristic.

In some embodiments, the present invention provides methods and compositions for increasing, decreasing, or otherwise altering the production of acids in plants. Examples of acids that can be altered according to the present invention include, but are not limited to, fumaric acid, malic acid, carbamic acid, glyceric acid, citric acid, ketoglutaric acid, quinic acid, shikimic acid, and sugar acids, such as, gluconic acid and galacturonic acid. The alterations in metabolic profiles are preferably accomplished by expressing, in a plant, one or more of the nucleic acid sequences in FIG. 1 or 2 (or sequences that hybridize thereto) shown to alter the production of acids in plants. In preferred embodiments, expression in plants of the sequences that hybridize to the preceding sequences also results in an increase, decrease, or alteration of acid production in a plant.

In some embodiments, the present invention provides methods and compositions for increasing, decreasing, or otherwise altering the production of fatty acids in plants. Examples of fatty acids that can be altered according to the present invention include, but are not limited to, tetradecanoic acid, hexadecanoic acid, heptadecanoic acid, octadecanoic acid, eicosanoic acid, docosanoic acid, 9-hexadecenoic acid, 6-octadecenoic acid, 9-octadecanoic acid, 7,10-hexadecadienoic acid, 9,12-octadecadienoic acid, 9,12,15-octadecatrienoic acid, and 7,10,13-docosatrienic acid. The alterations in metabolic profiles are preferably accomplished by expressing, in a plant, one or more of the nucleic acid sequences in FIG. 1 or 2 (or sequences that hybridize thereto) shown to alter the production of fatty acids in plants. In preferred embodiments, expression in plants of the sequences that hybridize to the preceding sequences also results in an increase, decrease, or alteration of fatty acid production in a plant.

In some embodiments, the present invention provides methods and compositions for increasing, decreasing, or otherwise altering the production of branched fatty acids in plants. Examples of branched fatty acids that can be altered according to the present invention include, but are not limited to, 14-methylhexadecanoic acid, 16-methylheptadecanoic acid, 17-methylheptadecanoic acid, 20-methylheneicosanoic acid, and 3,7,11,15-tetramethylhexadecanoic acid. The alterations in metabolic profiles are preferably accomplished by expressing, in a plant, one or more of the nucleic acid sequences in FIG. 1 or 2 (or sequences that hybridize thereto) shown to alter the production of branched fatty acids in plants. In preferred embodiments, expression in plants of the sequences that hybridize to the preceding sequences also results in an increase, decrease, or alteration of branched fatty acid production in a plant.

In some embodiments, the present invention provides methods and compositions for increasing, decreasing, or otherwise altering the production of alcohols in plants. Examples of alcohols that can be altered according to the present invention include, but are not limited to, octadecanol, phytol, valereneol, and sugar alcohols, such as, inositol and mannitol. The alterations in metabolic profiles are preferably accomplished by expressing, in a plant, one or more of the nucleic acid sequences in FIG. 1 or 2 (or sequences that hybridize thereto) shown to alter the production of alcohols in plants. In preferred embodiments, expression in plants of the sequences that hybridize to the preceding sequences also results in an increase, decrease, or alteration of alcohol production in a plant.

In some embodiments, the present invention provides methods and compositions for increasing, decreasing, or otherwise altering the production of alkaloids and other bases in plants. Examples of alkaloids and other bases that can be altered according to the present invention include, but are not limited to, nicotine, nornicotine, and 1,4-butanediamine. The alterations in metabolic profiles are preferably accomplished by expressing, in a plant, one or more of the nucleic acid sequences in FIG. 1 or 2 (or sequences that hybridize thereto) shown to alter the production of alkaloids (or other bases) in plants. In preferred embodiments, expression in plants of the sequences that hybridize to the preceding sequences also results in an increase, decrease, or alteration of alkaloid and other base production in a plant.

In some embodiments, the present invention provides methods and compositions for increasing, decreasing, or otherwise altering the production of alkenes and alkynes (unsaturated hydrocarbons) in plants. Examples of alkenes and alkynes that can be altered according to the present invention include, but are not limited to, limonene, dimethylcyclooctadiene, 4-methyldecene, eicosene, tetramethylhexadecene, dehydroisolongifolene, and squalene. The alterations in metabolic profiles are preferably accomplished by expressing, in a plant, one or more of the nucleic acid sequences in FIG. 1 or 2 (or sequences that hybridize thereto) shown to alter the production of alkenes and alkynes in plants. In preferred embodiments, expression in plants of the sequences that hybridize to the preceding sequences also results in an increase, decrease, or alteration of alkene and alkyne production in a plant.

In some embodiments, the present invention provides methods and compositions for increasing, decreasing, or otherwise altering the production of amino acids and related compounds in plants. Examples of amino acids and related compounds that can be altered according to the present invention include, but are not limited to, proline, glycine, alanine, serine, aspartic acid, glutamic acid, lysine, tyrosine, phenylalanine, valine, threonine, arginine, glutamine, tryptophan, isoleucine, and 5-oxo-proline. The alterations in metabolic profiles are preferably accomplished by expressing, in a plant, one or more of the nucleic acid sequences in FIG. 1 or 2 (or sequences that hybridize thereto) shown to alter the production of amino acids in plants. In preferred embodiments, expression in plants of the sequences that hybridize to the preceding sequences also results in an increase, decrease, or alteration of amino acid and related compounds production in a plant.

In some embodiments, the present invention provides methods and compositions for increasing, decreasing, or otherwise altering the production of carbohydrates in plants. Examples of carbohydrates that can be altered according to the present invention include, but are not limited to, arabinose, xylose, glucose, fructose, galactose, fructose, mannose, rhamnose, and sucrose. The alterations in metabolic profiles are preferably accomplished by expressing, in a plant, one or more of the nucleic acid sequences in FIG. 1 or 2 (or sequences that hybridize thereto) shown to alter the production of carbohydrates in plants. In preferred embodiments, expression of the sequences that hybridize to the preceding sequences also results in an increase, decrease, or alteration of carbohydrate production in a plant.

In some embodiments, the present invention provides methods and compositions for increasing, decreasing, or otherwise altering the production of esters in plants. Examples of esters that can be altered according to the present invention include, but are not limited to, acylates, such as 2-steroyl-1-acetylglycerol and 2-eicosanoyl-1-acetylglycerol. The alterations in metabolic profiles are preferably accomplished by expressing, in a plant, one or more of the nucleic acid sequences in FIG. 1 or 2 (or sequences that hybridize thereto) shown to alter the production of esters in plants. In preferred embodiments, expression in plants of the sequences that hybridize to the preceding sequences also results in an increase, decrease, or alteration of ester production in a plant.

In some embodiments, the present invention provides methods and compositions for increasing, decreasing, or otherwise altering the production of glycerides in plants. Examples of glycerides that can be altered according to the present invention include, but are not limited to, glycerol palmitate, and glycerol linoleate, and glyceryl linolenate. The alterations in metabolic profiles are preferably accomplished by expressing, in a plant, one or more of the nucleic acid sequences in FIG. 1 or 2 (or sequences that hybridize thereto) shown to alter the production of glycerides in plants. In preferred embodiments, expression in plants of the sequences that hybridize to the preceding sequences also results in an increase, decrease, or alteration of glyceride production in a plant.

In some embodiments, the present invention provides methods and compositions for increasing, decreasing, or otherwise altering the production of hydrocarbons (saturated) in plants. Examples of hydrocarbons that can be altered according to the present invention include, but are not limited to, eicosane, hentriacontane, 2-methyloctacosane, 3-methylnonacosane, 2-methyltriacontane, 3-methylhentriacontane, and 2-methyldotriacontane. The alterations in metabolic profiles are preferably accomplished by expressing, in a plant, one or more of the nucleic acid sequences in FIG. 1 or 2 (or sequences that hybridize thereto) shown to alter the production of hydrocarbons in plants. In preferred embodiments, expression in plants of the sequences that hybridize to the preceding sequences also results in an increase, decrease, or alteration of hydrocarbon production in a plant.

In some embodiments, the present invention provides methods and compositions for increasing, decreasing, or otherwise altering the production of phenols and related compounds in plants. Examples of phenols and related compounds that can be altered according to the present invention include, but are not limited to, caffeic acid and chlorogenic acid. The alterations in metabolic profiles are preferably accomplished by expressing, in a plant, one or more of the nucleic acid sequences in FIG. 1 or 2 (or sequences that hybridize thereto) shown to alter the production of phenols (and related compounds) in plants. In preferred embodiments, expression in plants of the sequences that hybridize to the preceding sequences also results in an increase, decrease, or alteration of phenol and related compounds in a plant.

In some embodiments, the present invention provides methods and compositions for increasing, decreasing, or otherwise altering the production of sterols, oxygenated terpenes, and other isoprenoids in plants. Examples of sterols, oxygenated terpenes, and other isoprenoids that can be altered according to the present invention include, but are not limited to, solanesol, cycloartenol, alpha-tocopherol, alpha-tocopherol quinone, beta-tocopherol, gamma-tocopherol, stigmastenol, cycloartenol, stigmastatriene, campesterol, cholesterol, sitosterol, stigmasterol, methylene-lophenol, methylene-cycloartenol, dimethylergostadienol, fucosterol, ergostenone, fucosterol, stigmastadienol, and lanosterol. The alterations in metabolic profiles are preferably accomplished by expressing, in a plant, one or more of the nucleic acid sequences in FIG. 1 or 2 (or sequences that hybridize thereto) shown to alter the production of sterols, oxygenated terpenes, or other isoprenoids, in plants. In preferred embodiments, expression in plants of the sequences that hybridize to the preceding sequences also results in an increase, decrease, or alteration of sterol, oxygenated terpene, and other isoprenoid production in a plant.

In some embodiments, the present invention provides methods and compositions for increasing, decreasing, or otherwise altering the production of ketones and quinones in plants. Examples of ketones and quinones that can be altered according to the present invention include, but are not limited to, 3-phytolmenadione and alpha-tocopherol quinone. The alterations in metabolic profiles are preferably accomplished by expressing, in a plant, one or more of the nucleic acid sequences in FIG. 1 or 2 (or sequences that hybridize thereto) shown to alter the production of ketones and quinones. In preferred embodiments, expression in plants of the sequences that hybridize to the preceding sequences also results in an increase, decrease, or alteration of ketone and quinone production in a plant.

VII. Identification of Homologs to Sequences

The present invention also provides homologs and variants of the sequences described above, but which may not hybridize to the sequences described above under conditions ranging from low to high stringency. In some preferred embodiments, the homologous and variant sequences are operably linked to an exogenous promoter. FIG. 3 provides BLASTX search results from publicly available databases. The relevant sequences are identified by Accession numbers in these databases. FIG. 4 contains the top BLASTX hits (identified by Accession number) versus all the amino acid sequences in the Derwent.TM. biweekly database. FIG. 5 contains the top BLASTN hits (identified by Accession number) versus all the nucleotide sequences in the Derwent.TM. biweekly database.

In some embodiments, the present invention comprises homologous nucleic acid sequences (SEQ ID NOs:1166-3702 and 4154-7388) identified by screening an internal database with SEQ ID NOs.1-1165, 3703-4153 and 7389-7554 at a confidence level of Pz<1.00E-20. These sequences are provided in FIG. 2. The headers list the sequence identifier for the sequence that produced the actual phenotypic hit first and the sequence identifier for the homologous contig second.

As will be understood by those skilled in the art, the present invention is not limited to the particular sequences of the homologs described above. Indeed, the present invention encompasses portions, fragments, and variants of the homologs as described above. Such variants, portions, and fragments can be produced and identified as described in Section III above. In particularly preferred embodiments, the present invention provides sequences that hybridize to SEQ ID NOs:1166-3702 and 4154-7388 under conditions ranging from low to high stringency. In other preferred embodiments, the present invention provides nucleic acid sequences that inhibit the binding of SEQ ID NOs:1166-3702 and 4154-7388 to their complements under conditions ranging from low to high stringency. Furthermore, as described above in Section IV, the homologs can be incorporated into vectors for expression in a variety of hosts, including transgenic plants.

Homolog contigs, FIG. 2 (as described in Example 16, Section D: Identification of Homologous Sequences) are formed from individual sequence runs belonging to clones whose sequences share a predefined level of nucleotide identity to each other and are presumably independent isolates of a single gene sequence. This list of clones composing any one homolog contig are the actual entities that are screened. If clones sharing homology to a particular hit sequence, FIG. 1, perform a very similar or identical function within the organism, then these clones should also result in metabolic alterations when tested in the metabolic screen. The data contained in FIG. 9 are provided as examples of the metabolic phenotype correlation between homolog clones. FIG. 9 shows the correlation between the homolog sequence pairs and the metabolic alterations observed in this invention. Biochemicals, common to both clones, are listed with the corresponding chemical and biochemical classes. The alterations, up-regulated or down-regulated, observed for these biochemicals are also reported.

EXAMPLES

Example 1

Construction of Tissue-Specific N. benthamiana cDNA Libraries

A. mRNA Isolation: Leaf, root, flower, meristem, and pathogen-challenged leaf cDNA libraries were constructed. Total RNA samples from 10.5 .mu.g of the above tissues were isolated by TRIZOL reagent (Life Technologies, Inc.; Rockville, Md.). The typical yield of total RNA was 1 mg PolyA.sup.+RNA and was purified from total RNA by DYNABEADS oligo (T).sub.25. Purified mRNA was quantified by UV absorbance at OD.sub.260 The typical yield of mRNA was 2% of total RNA. The purity was also determined by the ratio of OD.sub.260/OD.sub.280. The integrity of the samples had OD values of 1.8-2.0.

B. cDNA Synthesis: cDNA was synthesized from mRNA using the SUPERSCRIPT plasmid system (Life Technologies, Inc.; Rockville, Md.) with cloning sites of NotI at the 3' end and SalI at the 5' end. After fractionation through a gel column to eliminate adapter fragments and short sequences, cDNA was cloned into both GENEWARE.TM. vector p1057 NP and phagemid vector PSPORT in the multiple cloning region between NotI and XhoI sites. Over 20,000 recombinants were obtained for all of the tissue-specific libraries.

C. Library Analysis: The quality of the libraries was evaluated by checking the insert size and percentage from representative 24 clones. Overall, the average insert size was above 1 kb, and the recombinant percentage was >95%.

Example 2

Construction of Normalized N. benthamiana cDNA Library in GENEWARE.TM. Vectors

A. cDNA synthesis. A pooled RNA source from the tissues described above was used to construct a normalized cDNA library. Total RNA samples were pooled in equal amounts first, then polyA+RNA was isolated by DYNABEADS oligo (dT).sub.25. The first strand cDNA was synthesized by the Smart III system (Clontech; Palo Alto, Calif.). During the synthesis, adapter sequences with Sfi1a and Sfi1b sites were introduced by the polyA priming at the 3' end and 5' end by the template switch mechanism (Clontech; Palo Alto, Calif.). Eight .mu.g first strand cDNA was synthesized from 24 .mu.g mRNA. The yield and size were determined by UV absorbance and agarose gel electrophoresis.

B. Construction of Genomic DNA driver. Genomic DNA driver was constructed by immobilizing biotinylated DNA fragments onto streptavidin-coated magnetic beads. Fifty .mu.g genomic DNA was digested by EcoR1 and BamH1 followed by fill-in reaction using biotin-21-dUTP. The biotinylated fragments were denatured by boiling and immobilized onto DYNABEADS by the conjugation of streptavidin and biotin.

C. Normalization Procedure. Six .mu.g of the first strand cDNA was hybridized to 1 .mu.g of genomic DNA driver in 100 .mu.l of hybridization buffer (6.times.SSC, 0.1% SDS, 1.times. Denhardt's buffer) for 48 hours at 65.degree. C. with constant rotation. After hybridization, the cDNA bound on genomic DNA beads was washed 3 times by 20 .mu.l 1.times.SSC/0.1% SDS at 65.degree. C. for 15 min and one time by 0.1.times.SSC at room temperature. The cDNA bound to the beads was then eluted in 10 .mu.l of fresh-made 0.1N NaOH from the beads and purified by using a QIAGEN DNA purification column (QIAGEN GmbH; Hilden, Germany), which yielded 110 ng of normalized cDNA fragments. The normalized first strand cDNA was converted to double strand cDNA in 4 cycles of PCR with Smart primers annealed to the 3' and 5' end adapter sequences.

D. Evaluation of normalization efficiency. Ninety-six non-redundant cDNA clones selected from a randomly sequenced pool of 500 clones of a previously constructed whole seedling library were used to construct a nylon array. One hundred ng of the normalized cDNA fragments versus the non-normalized fragments were radioactively labeled by .sup.32P and hybridized to DNA array nylon filters. The hybridization images and intensity data were acquired by a PHOSPHORIMAGER (Amersham Pharmacia Biotech; Chicago, Ill.). Since the 96 clones on the nylon arrays represent different abundance classes of genes, the variance of hybridization intensity among these genes on the filter were measured by standard deviation before and after normalization. The results indicated that by using this type of normalization approach, a 1000-fold reduction in variance among this set of genes could be achieved.

E. Cloning of normalized cDNA into GENEWARE.TM. vector. The normalized cDNA fragments were digested by Sfi1 endonuclease, which recognizes 8-bp sites with variable sequences in the middle 4 nucleotides. After size fractionation, the cDNA was ligated into GENEWARE.TM. vector p1057 NP in antisense orientation and transformed into DH5.alpha. cells. Over 50,000 recombinants were obtained for this normalized library. The percentage of insert and size were evaluated by Sfi digestion of randomly picked 96 clones followed by electrophoresis on 1% of agarose gel. The average insert size was 1.5 kb, and the percentage of insert was 98% with vector only insertions of >2%.

F Sequence analysis of normalized cDNA library. Two plates of 96 randomly picked clones have been sequenced from the 5' end of cDNA inserts. One hundred ninety-two quality sequences were obtained after trimming of vector sequences and other standard quality checking and filtering procedure, and subjected to BLASTX search in DNA and protein databases. Over 40% of these sequences had no hit in the databases. Clustering analysis was conducted based on accession numbers of BLASTX matches among the 112 sequences that had hits in the databases. Only three genes (tumor-related protein, citrin, and rubit) appeared twice. All other members in this group appeared only once. This was a strong indication that this library is well-normalized. Sequence analysis also revealed that 68% of these 192 sequences had putative open reading frames using the ORF finder program (as described above), indicating possible full-length cDNA.

Example 3

Rice cDNA Library Construction in GENEWARE.TM. Vectors

Oryzae sativa var. Indica IR-7 was grown in greenhouses under standard conditions (12/12 photoperiod, 29.degree. C. daytime temp., 24.degree. C. night temp.). The following types of tissue were harvested, immediately frozen on dry ice and stored at -80.degree. C.: young leaves (20 days post sowing), mature leaves and panicles (122 days post sowing). Mature and immature root tissue (either 122 or 20 days post sowing) was harvested, rinsed in ddH.sub.2O to remove soil, frozen on dry ice and stored at -80.degree. C.

The following standard method (Life Technologies) was used for generation of cDNA and cloning. High quality total RNA was purified from target tissues using Trizol (LTI) reagent. mRNA was purified by binding to oligo (dT) and subsequent elution. Quality of mRNA samples is essential to cDNA library construction and was monitored spectrophotmetrically and via gel electrophoresis. 2-5 .mu.g of mRNA was primed with an oligo (dT)-NotI primer and cDNA was synthesized (no isotope was used in cDNA synthesis). SalI adaptors were ligated to the cDNA, which was then subjected to digestion with NotI. Restriction fragments were fractionated based on size and the first 10 fractions were measured for DNA quantity and quality. Fractions 6 to 9 were used for ligations. 100 ng of GENEWARE.TM. vector was ligated to 20 ng synthesized cDNA. Following ligations, the mixtures were kept at -20.degree. C. For transformation, 1 .mu.l to 10 .mu.l ligation reaction mixture was added to 100 .mu.l of competent E. coli cells (strain DH5.alpha.) and transformed using the heat shock method. After transformation, 900 .mu.l SOC medium was added to the culture and it was incubated at 37.degree. C. for 60 minutes. Transformation reactions were plated out on 22.times.22 cm LB/Amp agar plates and incubated overnight at 37.degree. C.

Example 4

Poppy cDNA Library Construction in GENEWARE.TM. Vectors

A. Plant Growth. A wild population of Papaver rhoeas resistant to auxin 2,4-Dichlorophenoxyacetic acid (2,4-D) was identified from a location in Spain and seed was collected. The seed was germinated and yielded a morphologically heterogeneous population. Leaf shape varied from deeply to shallowly indented. Latex color in some individuals was pure white when freshly cut, slowly changing to light orange then brown. Latex in other individuals was bright yellow or orange and rapidly changed to dark brown upon exposure to air. A single plant (PR4) with the white latex phenotype was used to generate the library.

B. RNA extraction. Approximately 1.5 g of leaves and stems were collected and frozen on liquid nitrogen. The tissue was ground to a fine powder and transferred to a 50 mL conical polypropylene screw cap centrifuge tube. Ten mL of TRIZOL reagent (Life Technologies, Inc.; Rockville, Md.) was added and vortexed at high speed for several minutes of short intervals until an aqueous mixture was attained. Two mL of chloroform was added and the suspension was again vortexed at high speed for several minutes. The tube was centrifuged 15 minutes at 3100 rpm in a tabletop centrifuge (GP Centrifuge, Beckman Coulter, Inc; Fullerton, Calif.) for resolution of the phases. The aqueous supernatant was then carefully transferred to diethylpyrocarbonate (DEPC)-treated 1.5 mL microtubes and total RNA was precipitated with 0.6 volumes of isopropanol. To facilitate precipitation, the solution was allowed to stand 10 minutes at room temperature after thorough mixing. Following centrifugation for 10 minutes at 8000 rpm in a microcentrifuge (model 5415C, Eppendorf AG, Hamburg), the pellet of total RNA was washed with 70% ethanol, briefly dried and resuspended in 200 .mu.L DEPC-treated deionized water. A 10 .mu.L aliquot was examined by non-denaturing agarose gel electrophoresis.

C. cDNA synthesis. To generate cDNA, approximately 50 .mu.g of total RNA was primed with 250 pmole of first strand oligo (TAIL: 5'-GAG-GAT-GTT-AAT-TAA-GCG-GCC-GCT-GCA-G(T).sub.23-3')(SEQ ID NO:7555) in a volume of 250 .mu.L using 1000 units of Superscript reverse transcriptase (Life Technologies, Inc.; Rockville, Md.) for 90 minutes at 42.degree. C. Phenol extraction was performed by adding an equal volume of phenol:chloroform:isoamyl alcohol (25:24:1 v/v), vortexing thoroughly, and centrifuging 5 minutes at 14,000 rpm in an Eppendorf microfuge. The aqueous supernatant phase was transferred to a fresh microfuge tube and the first strand cDNA:mRNA hybrids were precipitated with ethanol by adding 0.1 volume of 3 M sodium acetate and 2 volumes of absolute ethanol. After 5 minutes at room temperature, the tube was centrifuged 15 minutes at 14,000 rpm. The pellet was washed with 80% ethanol, dried briefly and resuspended in 100 .mu.L TE buffer (10 mM TrisCl, 1 mM EDTA, pH 8.0). After adding 10 .mu.L Klenow buffer (RE buffer 2, Life Technologies, Inc.; Rockville, Md.) and dNTPs (Life Technologies, Inc.; Rockville, Md.) to a final concentration of I MM, second strand cDNA was generated by adding 10 units of Klenow enzyme (Life Technologies, Inc.; Rockville, Md.), 2 units of RNase H (Life Technologies, Inc.; Rockville, Md.) and incubating at 37.degree. C. for 2 hrs. The buffer was adjusted with .beta.-nicotinamide adenine dinucleotide .beta.-NAD) by addition of E. coli ligase buffer (Life Technologies, Inc.; Rockville, Md.) and adenosine triphosphate (ATP, Sigma Chemical Company, St. Louis, Mo.) added to a final concentration of 0.6 mM. Double stranded phosphorylated cDNA was generated by addition of 10 units of E. coli DNA ligase (Life Technologies, Inc.; Rockville, Md.), 10 units of T4 polynucleotide kinase (Life Technologies, Inc.; Rockville, Md.) and incubating for 20 minutes at ambient temperature.

The double stranded cDNA was isolated through phenol extraction and ethanol precipitation, as described above. The pellet was washed with 80% ethanol, dried briefly and resuspended in a minimal volume of TE. The resuspended pellet was ligated overnight at 16.degree. C. with 50 pmole of kinased AP3-AP4 adapter (AP-3: 5'-GAT-CTT-AAT-TAA-GTC-GAC-GAA-TTC-3'/AP-4: 5'-GAA-TTC-GTCGAC-TTA-ATT-AA-3 ')(SEQ ID NOs:7556-7557) and 2 units of T4 DNA ligase (Life Technologies, Inc.; Rockville, Md.). Ligation products were amplified by 20 cycles of PCR using AP-3 primer and examined by agarose gel electrophoresis.

Expanded adapter-ligated cDNA was digested overnight at 37.degree. C. with PacI and NotI restriction endonucleases. The GENEWARE.TM. vector pBSG1056 (Large Scale Biology Corporation, Vacaville, Calif.) was similarly treated. Digested cDNA and vector were electrophoresed a short distance through low-melting temperature agarose. After visualizing with ethidium bromide and excising the appropriate fraction(s), the fragments were then isolated by melting the agarose and quickly diluting 5: I with TE buffer to keep from solidifying. The diluted fractions were mixed in the appropriate ratio (approximately 10:1 vector:insert ratio) and ligated overnight at 16.degree. C. using T4 DNA ligase. Characterization of the ligation revealed an average insert size of 1.27 kb. The ligation was transferred to LSBC, Inc. where large scale arraying was carried out. Random sequencing of nearly 100 clones indicated that about 40% of the inserts had full length open frames.

Example 5

ABRC Library Construction in GENEWARE.TM. Expression Vectors

Expressed sequence tag (EST) clones were obtained from the Arabidopsis Biological Resource Center (ABRC; The Ohio State University, Columbus, Ohio 43210). These clones originated from Michigan State University (from the labs of Dr. Thomas Newman of the DOE Plant Research Laboratory and Dr. Chris Somerville, Carnegie Institution of Washington) and from the Centre National de la Recherche Scientifique Project (CNRS project; donated by the Groupement De Recherche 1003, Centre National de la Recherche Scientifique, Dr. Bernard Lescure et al.). The clones were derived from cDNA libraries isolated from various tissues of Arabidopsis thaliana var Columbia. A clone set of 11,982 clones was received as glycerol stocks arrayed in 96 well plates, each with an ABRC identifier and associated EST sequence.

An ORF finding algorithm was performed on the EST clone set to find potential full-length genes. Approximately 3,200 full-length genes were found and used to make GENEWARE.TM. constructs in the sense orientation. Five thousand of the remaining clones (not full-length) were used to make GENEWARE.TM. constructs in the antisense orientation.

Full-length clones used to make constructs in the sense orientation were grown and DNA was isolated using Qiagen (Qiagen Inc.; Valencia, Calif. 91355) mini-preps. Each clone was digested with NotI and Sse 8387 eight base pair enzymes. The resultant fragments were individually isolated and then combined. The combined fragments were ligated into pGTN P/N vector (with polylinker extending from PstI to NotI--5' to 3'). For each set of 96 original clones approximately 192 colonies were picked from the pooled GENEWARE.TM. ligations, grown until confluent in deep-well 96-well plates, DNA prepped and sequenced. The ESTs matching the ABRC data was bioinformatically checked by BLAST and a list of missing clones was generated. Pools of clones found to be missing were prepared and subjected to the same process. The entire process resulted in greater than 3,000 full-length sense clones.

The negative sense clones were processed in the same manner, but ligated into pGTN N/P vector (with polylinker extending from NotI to PstI--5' to 3'). For each set of 96 original clones approximately 192 colonies were picked from the pooled GENEWARE.TM. ligations and DNA prepped. The DNA from the GENEWARE.TM. ligations was subjected to RFLP analysis using TaqI 4 base cutter. Novel patterns were identified for each set. The RFLP method was applied and only applicable for comparison within a single ABRC plate. This procedure resulted in greater than 6,000 negative sense clones.

The identified clones were re-arrayed, transcribed, encapsidated and used to inoculate plants.

Example 6

Regulatory Factors cDNA Library Construction in GENEWARE.TM. Vectors

Transcription factors represent a class of genes that regulate and control many aspects of plant physiology, including growth, development, metabolism and response to the environment. In order to analyze a collection of regulatory factor genes, the PCR-based methods described below were used to construct a library of such genes from Arabidopsis thaliana and Saccharomyces cerevisiae. In addition, clones containing genes corresponding to regulatory factors from N. benthamiana, Oryzae sativa and Papaver rhoeas were selected, based on cDNA sequence, from the libraries generated in GENEWARE.TM. vectors as described above.

A. Regulatory Factor Gene Targeting. Publicly accessible databases of genome sequence include data on a wide range of organisms, from microbes to human. Many of these databases include annotation along with gene sequences that predict function of the genes based on either experimental data or homology to characterized genes. The MIPS (Munich Information Center for Protein Sequences) database contains sequence information and annotation for both Arabidopsis thaliana and Saccharomyces cerevisiae genomes. Based on this annotation, open reading frame sequences of predicted yeast and Arabidopsis transcription factors were downloaded from MIPS and used for PCR primer design.

B. PCR Primer Design

18-20 base pairs of nucleotide sequences at both ends of each downloaded ORF were extracted and used to design the gene-specific portion of individual primers. In addition, flanking sequence and restriction sites were added to the ends of primers as shown in the following example:

TABLE-US-00001 SEQ ID NO:7558 5' primer GCCTTAATTAACTGCAGC atgtcgggtcgtgaagatgaag PacI ------- PstI 5' gene-specific sequence SEQ ID NO:7559 3' primer TTGATATCTAGAGCGGCCGCTTA tcatgtttcatcatcgaaatcatca EcoRV NotI ------ 3' gene-specific sequence XbaI

C. Arabidopsis and Yeast Template Preparation. Total RNA was isolated from flowers and apical meristems of the Arabidopsis ecotype Columbia using the Qiagen RNA-easy kit (Cat. no. 75162). mRNA was subsequently isolated from total RNA using the MACS mRNA isolation kit from Miltenyl Biotec (cat. no. 751-02). First strand cDNA was synthesized from 10 .mu.g of mRNA in the presence of Superscript II reverse transcriptase (Gibco BRL cDNA synthesis kit; cat. no. 18248-013) and NotI primer (5'-GACTAGTTCTAGATCGCGAGCGGCCGCCC(T).sub.30VN-3')(SEQ ID NO:7560). The second strand was synthesized based on the manufacturers instructions. This cDNA was diluted 1:5 prior to DNA amplification.

Since most yeast genes do not contain introns, genomic DNA was used directly as a template for PCR. Genomic DNA from S. cerevisiae S288C was obtained directly from Research Genetics (ResGen, an Invitrogen company, Huntsville Ala., catalog #40802).

D. PCR Amplification. 1 .mu.l of template DNA was subjected to PCR using the Hi Fi Platinum (hot start) DNA polymerase (Gibco-BRL cat. no. 11304-011) and gene-specific primers for each ORF. Each 50 .mu.l reaction contained: 5 .mu.l 10.times. buffer, 1 .mu.l of 10 mM dNTP, 2 .mu.l of 50 mM MgSO.sub.4, 1 .mu.l of template cDNA, 10 pmoles of each primer and 0.2 unit of Platinum Hi Fi DNA polymerase. PCR reactions were carried out in a MJ Research (Model PTC 200) thermal cycler programmed with the following conditions: 3 min at 95.degree. C. 30 cycles [95.degree. C. 30 sec., 50.degree. C. 30 sec., 72.degree. C. 3 min.] 72.degree. C. 10 min. Following PCR, reactions were stored at -20.degree. C. until ready for ligation.

D. Subcloning ORFs into GENEWARE.TM. Vectors. To minimize cost and the labor involved in cloning of individual ORF, PCR products containing different ORFs were cloned into the GENEWARE.TM. vectors as pooled DNAs. 30-75 PCR products were pooled, digested with PacI and NotI and purified from an agarose gel. Purified DNA was subsequently ligated into the GENEWARE.TM. vector (5PN-Cap digested with PacI and NotI). Single colonies were selected, grown and their DNA analyzed for the presence of insert. Inserts were gel purified and sequenced, and the sequence compared to the MIP protein database to confirm that they covered the complete ORF. Unique sequences representing various related genes were selected to cover different genes within a multi-gene family. The efficiency of pooled cloning ranged from 30-50% (i.e., 30-50 clones were identified from analysis of 100 pooled PCR products). Following sequence identification of the clones, PCR products that were not represented in the first round of cloning were subsequently pooled together and subjected to a second round.

Example 7

Other Libraries: Regulatory Gene Selection

For each of the cDNA libraries generated from N. benthamiana, Oryzae sativa and Papaver Rhoeas, a unigene set of clones was established. Following basic library construction, all DNA sequences were subjected to BLASTN analysis against each other. Sequences that showed perfect homology across a minimum of 50 base pairs were clustered together. At this level each cluster putatively represents a unique gene. The size of cluster varies depends on the size and complexity of sequence population (sequenced library). A cluster may have only one sequence member, or consist of hundreds of member sequences. The clone with 5'-most sequence in a cluster was then selected to represent the gene. A collection of all the 5'-most sequences or clones was established as the unigene set for that particular library. In the example illustrated below, 4 EST sequences were clustered, representing a putative gene. The EST Seq I contained the most sequence information toward the 5'-end, indicating that this clone had the longest insert relative to other cluster members. This process allows removal of redundant clones and selection of the longest and most-likely full-length clones for subsequent screens.

##STR00001##

Based on the analysis of the sequence, and annotations of each unigene from each library, all clones that were homologous to known regulatory genes/transcription factors were targets for selection. Depending on the level of homology, some of the clones represented well characterized regulatory genes; however, many of the selected clones had only a modest level of homology to known genes or genes of very distantly related organisms. It is believed that this selection process can increase the probability of gene discovery, and by eliminating non-relevant clones, increase screen efficiency.

Example 8

Trichoderma cDNA Library Construction in GENEWARE.TM. Vectors

A. Growth and Induction of Trichoderma harzianum rifai 1295-22. Cultures of Trichoderma harzianum rifai 1295-22 were obtained from ATCC (cat.# 20847) and propagated on PDA. Liquid cultures were inoculated and induced using a protocol derived from Vasseur et al. (Microbiology 141:767-774, 1995) and Cortes et al. (Mol. Gen. Genet. 260:218-225, 1998): agar-grown cells were used to inoculate a 100 mL culture in PDB and grown 48 hours at 29.degree. C. with agitation. Mycelia were harvested by centrifugation, transferred to Minimal Media (MM) +0.2% glucose, and incubated overnight at 29.degree. C. with agitation. Mycelia were harvested again by centrifugation, washed with MM, resuspended in MM and incubated 2 hours at 29.degree. C. with agitation. Mycelia were harvested again by centrifugation, divided into 2 aliquots, and used to inoculate 1)125 mL MM+0.2% glucose or 2) 125 mL MM +lmg/mL elicitor. Elicitor is a preparation of cell walls from Rhizoctonia solani grown in liquid culture and isolated according to Goldman et al. (Mol. Gen. Genet. 234:481-488, 1992). Induced and uninduced cultures were incubated at 29.degree. C. with agitation, harvested after 24 and 48 hours by filtration and immediately frozen in liquid nitrogen. Aliquots were assayed for induction using 2-D gel SDS-PAGE to compare induced and uninduced cultures. Both induced and uninduced (24 hours) tissue was used for subsequent RNA isolation and library construction.

B. RNA Isolation and Library Construction. mRNA isolation was accomplished by magnetically labeling polyA.sup.+RNA with oligo (dT) microbeads and selecting the magnetically labeled RNA over a column. The purified polyA.sup.+RNA was then used for cDNA synthesis using a modified version of the full-length enrichment reactions (cap-capture method) described by Seki et al. (Plant J. 15:707-720, 1998). Specifically, isolated mRNA was primed with NotI-oligo d(T) primer to synthesize the first strand cDNA. After the synthesis reaction, a biotin group was chemically introduced to the diol residue of the cap structure of the mRNA molecule. RNase I treatment was then used to digest the mRNA/cDNA hybrids, followed by binding of streptavidin magnetic beads. After this step, the full-length cDNAs were then removed from the beads by RNaseH and tailed with oligo dG by terminal transferase or used directly in the 2.sup.nd strand synthesis. For the oligo dG tailed samples, the second strand cDNAs were then synthesized with PacI-oligo dC primers and DNA polymerase. Additional modifications to the published procedure include: addition of trehalose and BSA as enzyme stabilizers in the reverse transcriptase reaction, a temperature of 50 to 60.degree. C. for the first strand cDNA synthesis reaction, high stringency binding and washing conditions for capturing biotinylated cap-RNA/cDNA hybrids and substitution of the cDNA poly (dG) tailing step with a Sal-I linker ligation. The cDNA was size-fractionated over a column and the largest 2-3 fractions were collected and used to ligate with GENEWARE.TM. vector pBSG1057. The ligation reaction was transformed into E. coli DH5.alpha. and plated, the transformation efficiency was calculated and the DNA from the transformants was subjected to the quality control steps described below: 1. cDNA synthesis/cloning: The cloning efficiency must be greater than 8.times.10.sup.5 cfu/.mu.g. 2. Restriction enzyme digestion and sequencing: 500 to 1,000 transformants were picked and DNA isolated. cDNA inserts were digested out by appropriate restriction enzymes and checked by gel electrophoresis. The average insert size was calculated from 100 random clones. If the average size was >0.9 kbp, the DNA preps were then passed on to the sequencing group to obtain 5'-end sequences. Those sequences were used to further evaluate the of the library. Libraries that did not meet QC standards, such as high vector background (>5%), low full-length percentage (<60%), or short average insert size (<0.7kbp), were discarded, and the entire procedure repeated.

C. Library Subtraction. The induced Trichoderma library in GENEWARE.TM. was constructed as above and a large number of clones were arrayed on a nylon membrane at high density (HD array). Based on the genomic size and expression levels of S. cerevisiae, 18,000 colonies were imprinted to provide 3-fold coverage of the expressed genes. Freshly grown colonies were plated out and picked into 384 well plates and then imprinted on Nylon membranes in 3.times.3 format at duplicated locations. First strand cDNAs to use as probes were synthesized from mRNAs isolated from both induced and uninduced tissue and used to hybridize the HD arrays. The intensity of each clone after hybridization was quantitated by phosphoimage scanning. The locations of all 18,000 spots were tracked by Array Vision software, which also determined the local background and calculated the signal/noise ratio for every clone on the membrane. The data generated were then converted to Excel format and analyzed to obtain the fold of induction or down-regulation. Based on the measured noise level, a 5-fold increase or decrease, relative to controls, was used as a cutoff value. Clones displaying .gtoreq.5-fold induction or reduction on duplicated samples were chosen. These clones were robotically re-arrayed using a Qbot device (see below, Colony Array) DNA was prepped as described below and sequenced. Based on the clustering results, 5'-most unigenes were selected and re-arrayed using the procedures described for the Poppy library above: the total number of clones that were selected was 1,019 for the up-regulated library (Th03), and 851 for the down-regulated library (Th04). These clones were prepared as described below (DNA Preparation, Transcription, Inoculation) and tested in a functional genomic screens for modified visual phenotypes.

Example 9

Colony Array

A. Colony Array--Picking. Ligations were transformed into E. coli DH5.alpha. cells and plated onto 22.times.22 cm Genetix "Q Trays" prepared with 200 mL agar, Amp.sup.100. A Qbot device (Genetix, Inc., Christchurch, Dorset UK) fitted with a 96 pin picking head was used to pick and transfer desired colonies into 384-well plates according to the manufacturers specifications and picking program SB384.SC1, with the following parameters:

Source

Container: Genetix bioassay tray

Color: White

Agar Volume: 200 mL

Destination

Container: Hotel (9 High)

Plate: Genetix 384 well plate

Time In Wells (sec): 2

Max Plates to use: # of 384 well plates

1.sup.st Plate: 1

Dips to Inoculate: 10

Well Offset: 1

Head

Head: 96 Pin Picking Head

First Picking Pin: 1

Pin Order: A1-H1, H2-A2 . . . (snaking)

Sterilizing

Qbot Bath #1

Bath Cycles: 4

Seconds in Dryer: 10

Wait After Drying: 10

(approximate picking time: 8 hrs /20,000 colonies)

Following picking, 384 well plates containing bacterial inoculum were grown in a HiGro chamber fitted with O.sub.2 at 30.degree. C., speed 6.5 for 12-14 hours. Following growth, plates were replicated using the Qbot with the following parameters, 2 replication runs per plate:

Source

Container: Hotel (9 High)

Plate: Genetix Plate 384 Well

Plates to replicate: 24

Start plate No.: 1

No. of copies: 1

Destination

Container: Universal Dest Plate Holder

Plate: Genetix Plate 384 Well

No. of Dips: 5

Head

Head: 384 Pin Gravity Gridding Head

Sterilizing

Qbot Bath #1

Bath cycles: 4

Seconds in Dryer: 10

Wait After Drying: 10

Airpore tape was placed over the replicated 384 well plates and the replicated plates were grown in the HiGro as above for 18-20 hours, sealed with foil tapes and stored at -80.degree. C.

B. Colony Array--Gridding. Membrane filters were soaked in LB/Ampicillin for 10 minutes. Filters were aligned onto fresh 22.times.22 cm agar plates and allowed to dry on the plates 30 min. in a Laminar flowhood. Plates and filters were placed in the Qbot and UV sterilized for 20 minutes. Following sterilization, plates/filters were gridded from 384 well plates using the Qbot according to the manufacturers specifications with the following parameters:

Gridding Routine

Name: 3.times.3

Source

Container: Hotel (9 High)

Plate: Genetix Plate 384 Well

Max Plates: 8

Inking time (ms): 1000

Destination

Filter holder: Qtray

Gridding Pattern: 3.times.3, non-duplicate, 8

Field Order: front 6 fields

No. Filters: up to 15

Max stamps per ink: 1

Max stamps per spot: 1

Stamp time (ms): 1000

No. Fields in Filter: 2

No. Identical Fields: 2

Stamps between sterilize: 1

Head: 384 pin gravity gridding head

Pin Height Adjustment: No change

Qbot Bath #1

Bath cycles: 4

Dry time: 10 (Seconds)

Wait After Drying: 10 (Seconds)

C. Plate Rearray. 384 well plates were rearrayed into deep 96 well block format using the Qbot according to the manufacturers instructions and the following rearray parameters X2 per plate:

Source

Container: Hotel (9 High)

Plate: Genetix Plate 384 Well

1.sup.st Plate: 1

Destination

Container: Universal Dest Plate Holder

Plate: Beckman 96 Deep Well Plate

1.sup.st plate: 1

Dips to Inoculate: 5

Well offset: 1

Max plates to use: 12 (or less)

Time in wells (sec): 2

Qbot Bath #1

Head: 96 pin picking head

First Picking Pin: 1

Pin Order: A1-H1, A2-H2, A3-H3 . . .

Bath cycles: 4

Sec. In dryer: 10

Wait after drying: 10

Following rearray, the 96-well blocks were covered with airpore tape and placed in incubator shakers at 37.degree. C., 500 rpm for a total of 24 hours. Plates were removed and used for DNA preparation.

Example 10

DNA Preparation

Plasmid DNA was prepared in a 96-well block format using a Qiagen Biorobot 9600 instrument (Qiagen; Valencia; Calif.) according to the manufacturers specifications. In this 96-well block format, 900 .mu.l of cell lysate was transferred to the Qiaprep filter and vacuumed 5 min. at 600 mbar. Following this vacuum, the filter was discarded and the Qiaprep Prep-Block was vacuumed for 2 min at 600 mbar. After adding buffer, samples were centrifuged for 5 min at 600 rpm (Eppendorf benchtop centrifuge fitted with 96-wp rotor) and subsequently washed .times.2 with PE. Elution was carried out for 1 minute, followed by a 5 min. centrifugation at 6000 rpm. Final volume of DNA product was approximately 75 .mu.l.

Example 11

Generation of Raw Sequence Data and Filtering Protocols

High-throughput sequencing was carried out using the PCT200 and TETRAD PCR machines (MJ Research; Watertown, Mass.) in 96-well plate format in combination with two ABI 377 automated DNA sequencers (PE Corporation; Norwalk, Conn.). The throughput at present is six 96-well plates per day. The quality of sequence data is improved by filtering the raw sequence output from sequencer. One criteria is to make sure that the unreadable bases are less than 10% of the total number of bases for any sequence and that there are no more than ten consecutive Ns in the middle part of the sequence (40-450). The sequences that pass these tests are defined as being of high quality. The second step for improving the quality of a sequence is to remove the vectors from the sequence. There are two advantages of this process. First, when locating the vector sequence, its position can be used to align to the input sequence. The quality of the sequence can be evaluated by the alignment between the vector sequence and the target sequence. Second, the removal of the vector sequence greatly improves the signal-to-noise ratio and makes the analysis of the resulting database search much easier. A third important pre-filtering step is to eliminate the duplicates in a library so it will speed up the analysis and reduce redundant analyses.

Example 12

Automated Transcriptions and Encapsidations

Plasmid DNA preparations were subjected to automated transcription reactions in a 96-well plate format using a Tecan Genesis Assay Workstation 200 robotic liquid handling system (Tecan, Inc.; Research Triangle Park, N.C.) according to the manufacturers specifications, operating on the Gemini Software (Tecan, Inc.) program "Automated_Txns.gem. For these reactions, reagents from Ambion, Inc. (Austin, Tex.) were used according to the manufacturers specifications at 0.4.times. reaction volumes. Following the robotic set-up of transcription reactions, 96-well plates were removed from the Tecan, shaken on a platform shaker for 30 sec., centrifuged in an Eppendorf tabletop centrifuge fitted with a 96-well plate rotor at 700 rcf for 1 minute and incubated at 37.degree. C. for 1.5 hours.

During the transcription reaction incubation, encapsidation mixture was prepared according to the following recipe:

TABLE-US-00002 1X Solution Sterile ddi H.sub.2O 100.5 .mu.L 1 M Sodium Phosphate 13.0 .mu.L TMV Coat Protein (20 mg/mL) 6.5 .mu.L 120 .mu.L per well

This mixture was placed in a reservoir of the Tecan and added to the 96-well plates containing transcription reaction following the incubation period using Gemini software program "9_Plates.gem". After adding encapsidation mixture, plates were shaken for 30 sec. on a platform shaker, briefly centrifuged as described above, and incubated at room temperature overnight. Prior to inoculation, encapsidated transcript was sampled and subjected to agarose gel analysis for QC.

Example 13

Infection of N. benthamiana Plants with GENEWARE.TM. Viral Transcripts

Plant Growth

N. benthamiana seeds were sown in 6.5 cm pots filled with Redi-earth medium (Scotts) that had been pre-wetted with fertilizer solution (147 kg Peters Excel 15-5-15 Cal-Mag (The Scotts Company; Marysville, Ohio), 68 kg Peters Excel 15-0-0 Cal-Lite, and 45 kg Peters Excel 10-0-0 MagNitrate in 596L hot tap H.sub.2O, injected (H. E. Anderson; Muskogee, Okla.) into irrigation water at a ratio of 200:1). Seeded pots were placed in the greenhouse for 1 d, transferred to a germination chamber, set to 27.degree. C., for 2 d (Carolina Greenhouses; Kinston, N.C.), and then returned to the greenhouse. Shade curtains (33% transmittance) were used to reduce solar intensity in the greenhouse and artificial lighting, a 1:1 mixture of metal halide and high pressure sodium lamps (Sylvania) that delivered an irradiance of approximately 220 .mu.mol m.sup.2s.sup.-1 was used to extend day length to 16 h and to supplement solar radiation on overcast days. Evaporative cooling and steam heat were used to regulate greenhouse temperature, maintaining a daytime set point of 27.degree. C. and a nighttime set point of 22.degree. C. At approximately 7 days post sowing (dps), seedlings were thinned to one seedling per pot and at 17 to 21 dps, the pots were spaced farther apart to accommodate plant growth. Plants were watered with Hoagland nutrient solution as required. Following inoculation, waste irrigation water was collected and treated with 0.5% sodium hypochlorite for 10 minutes to neutralize any viral contamination before discharging into the municipal sewer.

Inoculation of Plants

For each GENEWARE.TM. clone, 180 .mu.L of inoculum was prepared by combining equal volumes of encapsidated RNA transcript and FES buffer [0.1M glycine, 0.06 M K.sub.2HPO.sub.4, 1% sodium pyrophosphate, 1% diatomaceous earth (Sigma), and either 1% silicon carbide (Aldrich), or 1% Bentonite (Sigma)]. The inoculum was applied to three greenhouse-grown Nicotiana benthamiana plants at 14 or 17 days post sowing (dps) by distributing it onto the upper surface of one pair of leaves of each plant (.about.30 .mu.L per leaf). Either the first pair of leaves or the second pair of leaves above the cotyledons was inoculated on 14 or 17 dps plants, respectively. The inoculum was spread across the leaf surface using one of two different procedures. The first procedure utilized a Cleanfoam swab (Texwipe Co, NJ) to spread the inoculum across the surface of the leaf while the leaf was supported with a plastic pot label (3/4.times.5 2M/RL, White Thermal Pot Label, United Label). The second implemented a 3'' cotton tipped applicator (Calapro Swab, Fisher Scientific) to spread the inoculum and a gloved finger to support the leaf. Following inoculation the plants were misted with deionized water and maintained in a greenhouse.

Infection Scoring

At 13 days post inoculation (dpi), the plants were examined visually and a numerical score was assigned to each plant to indicate the extent of viral infection symptoms based on phenotypic characteristics. 0=no infection, 1=possible infection, 2=infection symptoms limited to leaves <50-75% fully expanded, 3=typical infection, 4=atypically severe infection, often accompanied by moderate to severe wilting and/or necrosis.

Example 14

Metabolic Screens

A. Sample Generation. Individual dwarf tobacco Nicotiana benthamiana, (Nb) plants were manually transfected with an unique DNA sequence at 14 or 17 days post sowing using the GENEWARE.TM. viral vector technology, Example 13. Plants were grown and maintained under greenhouse conditions.

Samples were grouped into sets of up to 96 samples per set for inoculation, harvesting and analysis. Each sample set included 8 negative control (reference samples), up to 80 unknown (test) samples, and 8 quality control samples.

B. Harvesting. At 14 days after infection, infected leaf tissue, excluding stems and petioles, was harvested from plants with an infection score of 3. Infected tissue was placed in a labeled, 50-milliliter (mL), plastic centrifuge tube containing a tungsten carbide ball approximately 1 cm in diameter. The tube was immediately capped, and dipped in liquid nitrogen for approximately 20 seconds to freeze the sample as quickly as possible to minimize degradation of the sample due to biological processes triggered by the harvesting process. Harvested samples were maintained at -80.degree. C. between harvest and analysis. Each sample was assigned a unique identifier, which was used to correlate the plant tissue to the DNA sequence that the plant was transfected with. Each sample set was assigned a unique identifier, referred to as the harvest or meta rack ID.

C. Extraction. Prior to analysis, the frozen sample was homogenized by placing the centrifuge tube on a mechanical shaker. The action of the tungsten carbide ball during approximately 30 seconds of vigorous shaking reduced the frozen whole leaf tissue to a finely homogenized frozen powder. Approximately 1 gram of the frozen powder was extracted with approximately 7.5 mL of a solution of isopropanol (IPA):water 70:30 (v:v), to achieve a 0.133 g/ml ratio, containing pentadecanoic acid ethyl ester (C 15:0 EE) and p-hydroxybenzioc acid as surrogate standards by shaking at room temperature for 30 minutes.

D. Fractionation. A 200 .mu.L aliquot of the IPA:water extract was transferred to a clean glass container and referred to as Fraction 2 (F2). A 1200 microliter (.mu.L) aliquot of the IPA:water extract was partitioned with 1200 .mu.L of hexane. The hexane layer was removed to a clean glass container. This hexane extract is referred to as Fraction 1 (F1). A 90 .mu.L aliquot of the hexane extracted IPA:water extract was removed to a clean glass container. This aliquot is referred to as Fraction 4 (F4). The remaining hexane extracted IPA:water extract is referred to as Fraction 3 (F3). Each fraction for each sample was assigned a unique fraction aliquot ID (sample name).

E. Sample Preparation & Data Generation

Fraction 1: The hexane extract was evaporated to dryness under nitrogen at room temperature. The sample containers were sealed and stored at 4.degree. C. prior to analysis, if storage was required. Immediately prior to capillary gas chromatographic analysis using flame ionization detection (GC/FID), the F1 residue was reconstituted with 180 .mu.L of hexane containing pentacosane and hexatriacontane which were used as internal standards for the F1 analyses. The chromatographic data files generated following GC separation and flame ionization detection were named with the Fraction 1 aliquot ID for each sample and stored in a folder named after the harvest rack (sample set) ID. FIG. 7a summarizes the GC/FID parameters used to analyze Fraction 1 samples.

Fraction 2: The F2 aliquot was evaporated to dryness under nitrogen at room temperature and reconstituted in heptane containing 2 internal standards, undecanoate methyl ester (C11:0) and lignoceroate methyl ester (C24:0). In general, Fraction 2 is designed to analyze esterified fatty acids, such as phospholipids, triacylglycerides, and thioesters. In order to analyze these compounds by GC/FED, they were transmethylated to their respective methyl esters by addition of sodium methoxide in methanol and heat. Excess reagent was quenched by the addition of a small amount of water, which results in phase separation. The fatty acid methyl esters (FAMEs) were contained in the organic phase. FIG. 7b summarizes the GC/FID parameters used to analyze Fraction 2 samples.

Fraction 3: The remaining hexane extracted IPA:water extract (F3) was evaporated under nitrogen at 55.degree. C. The residue was reconstituted with 400 .mu.L of pyridine containing hydroxylamine hydrochloride and the internal standards, n-octyl-.beta.-D-glucopyranoside and tetraclorobenzene (OXIME solution). The derivatization was completed by the addition of 400 .mu.L of the commercially available reagent (N,O-bis[trimethylsily]trifluoroacetamide)+1% trimethylchlorosilane (BSTFA +1% TMCS). The chromatographic data files generated following GC separation and flame ionization detection were named with the Fraction 3 aliquot ID for each sample and stored in a folder named after the harvest rack (sample set) ID. FIG. 7c summarizes the GC/FID parameters used to analyze Fraction 3 samples.

Fraction 4: The F4 aliquot was diluted with 90 .mu.L of distilled water and 20 .mu.L of an 0.1 N hydrochloric acid solution containing norvaline and sarcosine, which are amino acids that are used as internal standards for the amino acids analysis. Immediately prior to high performance liquid chromatographic analysis using fluorescence detection (HPLC/FLD), the amino acids in F4 are mixed in the HPLC injector at room temperature with buffered orthophtaldehyde solution, which derivatizes primary amino acids, followed by fluorenyl methyl chloroformate, which derivatizes secondary amino acids. Following HPLC separation and fluorescence detection, chromatographic data files were generated for each sample, named with a sequential number which can be tracked back to the F4 aliquot ID, and stored in a folder named after the harvest rack (sample set) ID. FIG. 7d summarizes the HPLC/FLD parameters used to analyze Fraction 4 samples.

F. Carbohydrate analysis by digestion of the Dionex extracted plant residue. Approximately 1 gram of frozen homogenized plant tissue was weighed into stainless steel extraction cartridges sandwiched between fiberglass filters. The samples were extracted with approximately 12 mL of 50:50 (v:v) isopropanol:water containing 0.1 N potassium hydroxide at 120.degree. C. and 2000 psi. The extracted sample residue was dried for 2 hours at 100.degree. C. The dry residue, 10 to 20 mg, was transferred from the extraction cartridge into a 13.times.100 mm glass tube containing 0.5 mL of 0.5 N hydrochloric acid in methanol and 0.12 mL of methyl acetate, blanketed with nitrogen, then sealed with TEFLON coated screw cap and heated for 16 hours at 80.degree. C. The liquid phase was then transferred (using an 8 channel pipetor) to a glass insert supported by a 96 well aluminum block plate containing 10 uL of t-butanol, which was then evaporated to dryness under hot flowing nitrogen. The methyl-glycosides and methyl-glycosides methyl ester residues generated by the methanolysis were silylated in 0.1 mL pyridine and 0.1 mL of the commercially available reagent (N,O-bis[Trimethylsily]trifluoroacetamide)+1% Trimethylchlorosilane (BSTFA+1% TMCS) at room temperature for one hour. The derivatized sample is analyzed by GC/FID using a DB1 capillary column (15 meters, 250 microns I.D., 0.25 microns film thickness) with an 11 minute temperature program (160.degree. C. to 190.degree. C. at 5.degree. C./min, then 190.degree. C. to 298.degree. C. at 36.degree. C. /minute followed by a 2 minute hold). The dual injection GC analyzes 192 samples/day from two 96 well plates. Table 1 lists the carbohydrates that have been identified and quantitated in the sample extracts.

G. Data Analysis & Hit Detection. Two complementary methods were used to identify modifications in the metabolic profile of samples. These data analysis methods are called automated data analysis (ADA) and quantitative data analysis (QUANT. ADA was used to identify hits in Fractions 1, 2, 3, and 4; whereas, QUANT was used to identify hits in the Carbohydrate Fraction and Fractions 2 and 4. If either method identified a fraction as a hit, the nucleic acid sequence used to transfect the plant tissue from which that fraction was obtained is classified as a hit.

ADA employs a qualitative pattern recognition approach using ABNORM (The Dow Chemical Company). ABNORM is described in the U.S. Pat. No. 5,592,402 (incorporated herein by reference). ADA was performed on chromatograms from Fractions 1 through 4. The ADA process developed a statistical model from a subset of all the available chromatograms within an analysis group. The chromatograms in this subset are defined to be "normal" and a chromatogram cannot be part of the normal subset if it fails any one of the following three tests. The first test involves alignment. A chromatogram may fail to aligned if, A) correlation is too small in any one of the correlation windows, B) extrapolation is too great as defined by the user, or C) maximum correlation is found too near the edge of the correlation window. The second test involves the internal standard. If a chromatogram has any one of the user selected internal standards with an area count more than 5 median absolute deviations (MADs) from the median area count of the chromatograms in the analysis group it is marked as "bad" and can not be included in the normal subset. The third test involves the total area count. If a chromatogram has a total area count more than 5 MADs below the median total area count of the chromatograms in the analysis

Chromatograms that pass the above three tests are defined as "good" samples. The normal subset was further determined as a certain percent (user defined) of the "good" samples that are nearest in terms of the 2-norm. The model developed on the "normal" subset was then used to test all of the chromatograrns in the analysis group for statistically significant differences from the model.

Updated models for each fraction were generated for each sample set. Detection limits were generated for each of the fractions.

Quantitative data analysis was based on individual peak areas. Quantitative data analysis was applied to specific compounds of interest in Fraction 2 (fatty acids), Fraction 4 (amino acids) and carbohydrate analysis. Peaks for these compounds were identified based on retention time. The peak areas corresponding to these compounds were generated. For Fraction 2 and the Carbohydrate Fraction, the relative percent of the peak areas for the compounds in Table I were calculated. The average ( x) and standard deviation (STD) of the relative % of the peak areas for the individual compounds were calculated from the reference chromatograms generated with the sample set. The average and STD were used to calculate a range for each compound. Depending on the compound, this range was typically x+/-3 or 5 STDs. If the relative percent of the peak area from a sample was outside this range, the compound was considered to be significantly different from the `normal` level and the sample was identified as a hit for F2 or the Carbohydrate Fraction. For Fraction 4, the concentration, in micrograms/gram was calculated for each of the amino acids listed in Table 1 from calibration standards analyzed within the same sample set. The amino acid concentrations from references were used to calculate the `normal` range from the x and STD for each amino acid. If the amino acid concentration for a sample fell outside this range, the amino acid was considered to be different from `normal` and the sample was identified as a hit for F4.

TABLE-US-00003 TABLE 1 Tobacco Metabolites Monitored in Fractions 2, 4 & Carbohydrates by Quantitative Analysis Fraction 4 Fraction 2 (Fatty Acids) (Amino Acids) Carbohydrates undecanoic acid methyl ester* C11:0 Aspartic Acid ASP Arabinose Pentadecanoic acid methyl C15:0 Glutamic Acid GLU Rhamnose ester** Pentadecanoic acid ethyl ester** C15:0 Serine SER Xylose palmitic acid methyl ester C16:0 Histidine HIS Mannose palmitoleic acid methyl ester C16:1 Glycine GLY Galactose iso methylpentadecanoic acid C16:0:Me Threonine THR Galacturonic Acid methyl ester palmitoleic acid methyl ester C16:2 Alanine ALA Glucose palmitolenic acid methyl ester C16:3 Arginine ARG iso methylhexadecanoic acid C17:0Me Tyrosine TYR methyl ester Stearic acid methyl ester C18:0 Cystine CY2 Oleic acid methyl ester C18:1 Valine VAL Linoleic acid methyl ester C18:2 Methionine MET Linolenic acid methyl ester C18:3 Norvaline* NVA Arachidic acid methyl ester C20:0 Tryptohane TRP Lignoceric acid methyl ester* C24:0 Phenylalanine PHE Isoleucine ILE Leucine LEU Lysine LYS Sarcosine* SAR Proline PRO

H. Shipping Hits. Shipping Hits. Any F1, F2, F3, or Carbohydrate Fractions identified as hits by ADA and/or quantitative analysis, and the most typical null for each fraction for each sample set as identified by ADA, were sent to the Function Discovery Laboratory of Analytical Sciences Capability within Corporate R&D (The Dow Chemical Company; Midland, Mich.) for structural characterization and quantification of relative change of the specific biochemical compounds altered (see Example 15). Samples were sealed, packaged on dry ice and shipped for overnight delivery.

Example 15

Identification of Metabolic Changes

This Example describes the identification of the chemical nature of genetic modifications made in tobacco plants using GENEWARE.TM. viral vector technology. The protocols involved the use of gas chromatography/mass spectrometry (GC/MS) for the analyses of three primary fractions obtained from extraction and fractionation processes.

A. Methods. Major instruments and accessories used included bioinformatics computer programs (see the description of the Maxwell program in WO 02/10486, hereby incorporated by reference); mass spectral libraries [includes, Biotech FDL, which is also described in WO 02/10486, and two commercial libraries: NIST Standard Reference Database-NIST98-(National Institute of Standards and Technology) and the Wiley Registry of Mass Spectral Data-WILEY275-(John Wiley and Sons, Inc.)]; biotechnology database (FDL is described in WO 02/10486)-the FDL Biotechnology Database is based on the MICROSOFT ACCESS database program from MICROSOFT (Redmond, Wash.) and utilizes ACCORD FOR ACCESS (available from Accelrys Inc.; San Diego, Calif.) to incorporate chemical structures; BLIMS, a customized LIMS (Nautilus 99; Thermo LabSystems Ltd., Manchester, England) for sample tracking and information transfer; biotechnology database (eBRAD; Dow/DAS/LSBC) an ORACLE (Redwood Shore, Calif.) based rational database that is a depository containing data from various screens and associated sequencing data; HP Model 6890 capillary gas chromatograph (GC; Agilent Technologies); HP Model 5973 Mass Selective Detector (MSD; Agilent Technologies); auto-sampler and sample preparation station (LEAP Technologies); large volume injector system (APEX); Ultra Freezer (Revco); and ChemStation GC/MS Software (Agilent Technologies).

Subject samples (those exhibiting an altered metabolic characteristic) and corresponding References (also referred to as controls or nulls) were shipped via overnight mail from the Metabolic Screening Laboratory (Indianapolis, Ind.) to the Function Discovery Laboratory (Midland, Mich.). Samples were removed from the shipping container, inspected for damage, and then placed in a freezer until analysis by GC/MS.

Samples were received in vials or in titer plates with a titer plate (TP) number, also referred to as a Rack Identification number that was used to track the sample in the BLIMS system. The titer plate number was used by the FDL to extract from BLIMS pertinent information from ADA (Automated Data Analysis chromatographic pattern recognition software) HIT reports and/or QUANT (a quantitative data analysis approach that makes use of individual peak areas of select peaks corresponding to specific compounds of interest in the fatty acid Fraction 2) HIT reports generated by the Metabolic Screening Laboratory. The information in these reports included the well position of the respective HITs (Subject), the corresponding well position of the Reference, and other pertinent information, such as, aliquot identification. This information was used to generate ChemStation and LEAP sequences for FDL analyses.

Samples were sequenced for analysis in the following order:

TABLE-US-00004 TABLE 2 Analysis Order Solvent Blank Instrument Performance Standard Subjects and Associated Reference . . . Performance Standard Solvent Blank

Samples were analyzed on GC/MS systems using the following procedures. Fraction 1 samples were shipped dry and required a hexane reconstitution step. Fraction 2 and Fraction 3 samples were analyzed as received. Internal standards and surrogate standards were added to the samples prior to GC/FD analysis (see Example 14).

B. Fraction 1 Analysis. The name of the GC/MS method used is BIONEUTx (where x is a revision number of the core GC/MS method). The method is retention-time locked to the retention time of pentacosane, an internal standard, using the ChemStation RT Locking algorithm. Internal/Surrogate* Standard(s) Pentacosane Hexatriacontane *Pentadeconoic acid, ethyl ester

TABLE-US-00005 Chromatography Column: J&W DB-5MS 50M .times. 0.320 mm .times. 0.25 .mu.m film Mode: constant flow Flow: 2.0 mL/min Detector: MSD Outlet psi: vacuum Oven: 40.degree. C. for 2.0 min 20.degree. C./min to 350.degree. C., hold 15.0 min Equilibration time: 1 min Inlet: Mode: split Inj Temp: 250.degree. C. Split ratio: 50:1 Gas Type: Helium

TABLE-US-00006 LEAP Injector: Injector: Inj volume: optimized to pentacosane peak intensity (typically 20 .mu.L) Sample pumps: 2 Wash solvent A: Hexane Wash solvent B: Acetone Preinj Solvent A washes: 2 Preinj Solvent B washes: 2 Postinj Solvent A washes: 2 Postinj Solvent B washes: 2

TABLE-US-00007 APEX Injector Method Name: BIONEUTx (where x is a revision number of the core APEX method). Modes: Initial: Standby (GC Split) Splitless: (Purge Off) 0.5 min GC Split: (Standby) 4 min ProSep Split: (Flow Select) 23 min Temps: 50.degree. C. for 0.0 min. 300.degree. C./min to 350.degree. C., hold for 31.5 min

TABLE-US-00008 Mass Spectrometer Scan: 35-800 Da at sampling rate 2 (1.96 scans/sec) Solvent delay: 4.0 min Detector: EM absolute: False EM offset: 0 Temps: Transfer line: 280.degree. C. Ion source: 150.degree. C. MS Source: 230.degree. C.

C. Fraction 2 Analysis: The name of the GC/MS method used is BIOFAMEx (where x is a revision number of the core GC/MS method). The method is retention-time locked to RT of undecanoic acid, methyl ester, an internal standard, using the ChemStation RT Locking algorithm. Internal/Surrogate* Standard(s) Undecanoic acid, methyl ester Tetracosanoic acid, methyl ester *Pentadecanoic acid, methyl ester *Pentadecanoic acid, ethyl ester

TABLE-US-00009 Chromatography Column: J & W DB-23 FAME 60M .times. 0.250 mm .times. 0.15 .mu.m film Mode: constant flow Flow: 2.0 mL/min Detector: MSD Outlet psi: vacuum Oven: 50.degree. C. for 2.0 min 20.degree. C./min to 240.degree. C., hold 10.0 min Equilibration time: 1 min Inlet: Mode: split Inj Temp: 240.degree. C. Split ratio: 50:1 Gas Type: Helium

TABLE-US-00010 LEAP Injector: Injector: Inj volume: optimized to undecanoic acid, methyl ester peak intensity (Typically 10 .mu.L) Sample pumps: 2 Wash solvent A: Methanol Wash solvent B: Methanol Preinj Solvent A washes: 2 Preinj Solvent B washes: 2 Postinj Solvent A washes: 2 Postinj Solvent B washes: 2

TABLE-US-00011 APEX Injector Method Name: BIOFAMEx (where x is a revision number of the core APEX method). Modes: Initial: GC Split Splitless: 0.5 min GC Split: 4 min ProSep Split: 21 min Temps: 60.degree. C. for 0.5 min. 300.degree. C./min to 250.degree. C., hold for 20 min 300.degree. C./min to 260.degree. C., hold for 5 min

TABLE-US-00012 Mass Spectrometer Scan: 35-800 Da at sampling rate 2 (1.96 scans/sec) Solvent delay: 4.5 min Detector: EM absolute: False EM offset: 0 Temps: Transfer line: 200.degree. C. Ion source: 150.degree. C. MS Source: 230.degree. C.

D. Fraction 3 Analysis. The name of the GC/MS method used is BIOAQUAx (where x is a revision number of the core GC/MS method). Method is retention-time locked to the RT of n-octyl-.beta.-D-glucopyranoside, an internal standard, using the ChemStation RT Locking algorithm.

Internal/Surrogate* Standard(s)

n-Octyl-.beta.-D-glucopyranoside *Tetrachlorobenzene *p-Hydroxybenzoic acid

TABLE-US-00013 Chromatography Column: Chrompack 7454 CP-SIL 8 60M .times. 0.320 mm .times. 0.25 .mu.m film Mode: constant flow Flow: 2.0 mL/min Detector: MSD Outlet psi: vacuum Oven: 40.degree. C. for 2.0 min 20.degree. C./min to 350.degree. C., hold 10.0 min Equilibration time: 1 min Inlet: Mode: split Inj Temp: 250.degree. C. Split ratio: 50:1 Gas Type: Helium

TABLE-US-00014 LEAP Injector: Injector: Inj volume: Optimized to n-octyl-.beta.-D-glucopyranoside peak intensity (Typically 2.5 .mu.L) Sample pumps: 2 Wash solvent A: Hexane Wash solvent B: Acetone Preinj Solvent A washes: 2 Preinj Solvent B washes: 2 Postinj Solvent A washes: 2 Postinj Solvent B washes: 2

TABLE-US-00015 APEX Injector Method Name: BIOAQUAx (where x is a revision number of the core APEX method). Modes: Initial: GC Split Splitless: 0.5 min GC Split: 4 min ProSep Split: 20 min Temps: 60.degree. C. for 0.5 min. 300.degree. C./min to 350.degree. C., hold for 21.1 min

TABLE-US-00016 Mass Spectrometer Scan: 35-800 Da at sampling rate 2 (1.96 scans/sec) Solvent delay: 4.0 min Detector: EM absolute: False EM offset: 0 Temps: Transfer line: 280.degree. C. Ion source: 150.degree. C. MS Source: 230.degree. C.

E. Performance Standard: Two mixtures were used as instrument performance standards. One standard was run with Fraction 1 and 3 samples and the second was run with Fraction 2 samples. Below is the composition of the standards as well as approximate retention time values observed when run under the GC/MS conditions previously described. These retention time values are subject to change depending upon specific instrument and chromatographic conditions.

TABLE-US-00017 TABLE 3 Fraction 1 and 3 Performance Standard Time Compound 6.25 dimethyl malonate 7.25 dimethyl succinate 8.15 dimethyl glutarate 8.98 dimethyl adipate 11.06 dimethyl azelate 11.42 hexadecane 11.70 dimethyl sebacate 13.57 eicosane 15.36 tetracosane 16.88 octacosane 18.26 dotriacontane 19.95 hexatriacontane

TABLE-US-00018 TABLE 4 Fraction 2 Performance Standard Time Compound 8.82 undecanoic acid, methyl ester 9.32 dodecanoic acid, methyl ester 10.24 tetradecanoic acid, methyl ester 11.07 hexadecanoic acid, methyl ester 11.84 octadecanoic acid, methyl ester 11.90 oleic acid, methyl ester 12.14 linoleic acid, methyl ester 12.39 linolenic acid, methyl ester 12.60 eicosanoic acid, methyl ester 13.42 docosanoic acid, methyl ester

F. Data Analysis. Subject and Reference data sets were processed using the Bioinformatics computer program Maxwell. The principal elements of the program are 1) Data Reduction, 2) Two-Dimensional Peak Matching, 3) Quantitative Peak Differentiation (Determination of Relative Quantitative Change), 4) Peak Identification, 5) Data Sorting, and 6) Customized Reporting.

The program queries the user for the filenames of the Reference data set and Subject data set(s) to compare against the Reference. A complete listing of user inputs with example input is shown below.

TABLE-US-00019 TABLE 5 Bioinformatics Analysis USER QUERY EXAMPLE USER INPUT Operator Name M. Maxwell Total number of data files to process 5 Which Fraction 3 Reference (Control) File Name AAPR0020.D Process a specific RT Range Y Specific RT range 6.5-23 Internal Standard Retention Time 14.902 +/- variation in Internal Std. RT .004 Variation in peak RI, ChemStation .005 Percent variation in peak RI, Biotech .010 Database Threshold for determining Area % change 60 Spectral Matching Value (Threshold MS- .95 XCR for peaks to be a match) Percent to determine LOP-PM* Value 1 Percent to determine LOP-SRT** Value 3 Quality Level for Library (Library match) 80 Subtract Background Y Time Range for Background 21.5-22.6 SHORT SUMMARY (y/n, y = no Y chromatograms) *LOP-PM - Limit of Processing for Peak Matching **LOP-SRT - Limit of Processing for Sorting

The program integrates the Total Ion Chromatogram (TIC) of the data sets using Agilent Technologies HP ChemStation RTE integrator parameters determined by the analyst. The corresponding raw peak areas are then normalized to the respective Internal Standard peak area. It should be noted that before the normalization is performed, the program chromatographically and spectrally identifies the Internal Standard peak. Should the identification of the Internal Standard not meet established criteria for a given Fraction, then the data set will not be further processed and it will be flagged for analyst intervention.

Peak tables from the Reference and each Subject were generated. The peak tables are comprised of retention time (RT), retention index (RI)--the retention time relative to the Internal Standard RT, raw peak areas, peak areas normalized to the Internal Standard, and other pertinent information.

The first of two filtering criteria, established by the analyst was then invoked and must be met before a peak is further processed. The criterion is based upon a peak's normalized area. All normalized peaks having values below the Limit of Processing for Peak Matching (LOP-PM), were considered to be "background". These "peaks" were not carried forth for any type of mathematical calculation or spectral comparison.

In the initial peak-matching step, the Subject peak table was compared to the Reference peak table and peaks between the two were paired based upon their respective RI values matching one another (within a given variable window). The next step in the peak matching routine utilized mass spectral data. Subject and Reference peaks that have been chromatographically matched were then compared spectrally. The spectral matching was performed using a mass spectral cross-correlation algorithm within the Agilent Technologies HP ChemStation software. The cross-correlation algorithm generates an equivalence value based upon spectral "fit" that was used to determine whether the chromatographically matched peaks are spectrally similar or not. This equivalence value is referred to as the MS-XCR value and must meet or exceed a predetermined value for a pair of peaks to be "MATCHED", which means they appear to be the same compound in both the Reference and the Subject. The MS-XCR value can also be used to judge peak purity. This two-dimensional peak matching process was repeated until all potential peak matches were processed. At the end of the process, peaks are categorized into two categories, MATCHED and UNMATCHED.

A second filtering criterion was next invoked, again based upon the normalized area of the MATCHED or UNMATCHED peak. For a peak to be reported and further processed, its normalized area must meet or exceed the predetermined Limit of Processing for Sorting (LOP-SRT).

Peaks that are UNMATCHED are immediately flagged as different. UNMATCHED peaks are of two types. There are those that are reported in the Reference but appear to be absent in the Subject (based upon criteria for quantitation and reporting). These peaks were designated in the Analyst Report with a percent change of "-100 percent" and the description "UNMATCHED IN SAMPLE". The second types of peaks are those that were not reported in the Reference (again, based upon criteria for quantitation and reporting) but were reported in the Subject, thus appearing to be "new" peaks. These peaks were designated in the Analyst Report with a percent change of "100 percent" and the description "NEW PEAK UNMATCHED IN NULL".

MATCHED peaks were processed further for relative quantitative differentiation. This quantitative differentiation is expressed as a percent change of the Subject peak area relative to the area of the Reference peak. A predetermined threshold for change must be observed for the change to be determined biochemical and statistically significant. The change threshold is based upon previously observed biological and analytical variability factors. Only changes above the threshold for change were reported.

Peaks were then processed through the peak identification process as follows. The mass spectra of the peaks were first searched against mass spectral plant metabolite libraries. The equivalence value assigned to the library match was used as an indication of a proper identification.

To provide additional confirmation to the identity of a peak, or to suggest other possibilities, library hits were searched further against a Biotechnology database (FDL; Dow). The Biotechnology database is based on the MICROSOFT ACCESS database program from MICROSOFT and utilizes ACCORD FOR ACCESS (available from Accelrys Inc.) to incorporate chemical structures into the database.

The Compound Identification Number (CIN) of the compound from the library was searched against those contained in the database. If a match was found, the CIN in the database was then correlated to the data acquisition method for that record. If the method was matched, the program then compared the retention index (RI), in the Peak Table, of the component against the value contained in the database for that given method. Should the RI's match (within a given window of variability) then the peak identity was given a high degree of certainty. Components in the Subject that are not identified by this process were assigned a unique Compound Identification Number based upon Fraction Number and RI (example: F1-U0.555). The unique identifier was used to track unknown components. The program then sorts the data and generates an Analyst Report.

An Analyst Report is an interim report consisting of PBM algorithm match quality value (equivalence value), RT, Normalized Peak Area, RI (Sample), RI (database) Peak Identification status [peak identity of high certainty (peaks were identified by the program based on the pre-established criteria) or criteria not met (program did not positively identify the component)], Component Name, CIN (a unique identifier, which could be a CAS number), Mass Spectral Library (containing spectrum most closely matched to that of the component), Unknown ID (unique identifier used to track unidentified components), MS-XCR value, Relative % Change, Notes (MATCHED/UNMATCHED), and other miscellaneous information. The Analyst Report was reviewed manually by the analyst who determined what further analysis was necessary. The analyst also generated a modified report, for further processing by the program, by editing the Analyst Report accordingly.

For Fractions 2 and 3, derivatization procedures were performed prior to analysis to make the certain components more amenable to gas chromatography. Thus, the compound names in the modified analyst report (MAR) were those of the derivatives. To accurately reflect the true components of these fractions, the MAR was further processed using information contained in an additional database. This database cross-references the observed derivatized compound to that of the original, underivatized "parent" compound by way of their respective compound identification numbers and replaces derivatives with parent names and information for the final report.

The Modified Analyst Report also contains a HIT Score of 0, 1, or 2. The value is assigned by the analyst to the data set of the Subject aliquot based on the following criteria:

TABLE-US-00020 0 No FDL data on Subject 1 FDL data collected; Subject not FDL HIT 2 FDL data collected; Subject is FDL HIT

An FDL HIT is defined as a reportable percent change (modification) observed in a Subject relative to Reference in a component of biochemical significance.

An electronic copy of the final report is entered into the Nautilus LIMS system (BLIMS) and subsequently into eBRAD (Biotech database). The program also generated a hardcopy of the pinpointed TIC and the respective mass spectrum of each component that was reported to have changed.

"NQ" and "NEW" are two terms used in the final report. Both terms refer to UNMATCHED peaks whose percent changes cannot be reported in a numerically quantitative fashion. These terms are defined as follows: "NQ" is used in the case where there was a peak reported in the Reference for which there was no match in the Subject (either because there was no peak in the Subject or, if there was, the area of the peak did not satisfy the Limit of Processing for Peak Matching). The percent change designation of "-100%" used in the Analyst report is replaced with "NQ". "NEW" is used in those situations where a peak was reported in the Subject but for which there was no corresponding match in the Reference (either because there was no peak in the Reference or, if there was, the area of the peak did not satisfy the Limit of Processing for Peak Matching). For these situations, the percent change designation of "100%" used in the Analyst Report is replaced with "NEW". The designation of "NEW" in the final report to a component that is present in the Sample but not in the Reference was necessary to eliminate any ambiguity with the appearance of "100%" for MATCHED peaks. A "100%" designation in the final report exclusively refers to a component with modification that doubled in the Subject relative to the Reference.

G. Results. The results of the identification of metabolic changes are shown in FIG. 6.

Example 16

Bioinformatic Analysis of Hits

A. Phred and Phrap: Phred is a UNIX based program that can read DNA sequencer traces and make nucleotide base calls independent of any software provided by the DNA sequencer manufacturer. Phred also provides a quality score for each base that can be used by the investigator to trim those sequences or preferably by Phrap to help its assembly process.

Phrap is another UNIX based program which takes the output of Phred and tries to assemble the individual sequencing runs into larger contiguous segments on the assumption that they all belong to a single DNA molecule. While this is clearly not the case with collections of Expressed Sequence Tags (ESTs) or with heterogeneous collections of sequencing runs belonging to more than one contiguous segment, the program does a very good job of uniquely assembling these collections with the proper manipulation of its parameters (mainly-penalty and -minscore; settings of 15 and 40 respectively provide contiguous sequences with exact homology approaching 95% over lengths of approximately 50 nucleotide base pairs or more). As with all assemblies it is possible for proper assemblies to be missed and for improper assemblies to be constructed, but the use of the above parameters and judicious-use of input sequences will keep these to a minimum.

Detailed descriptions of the Phred and Phrap software and it's use can be found in the following references which are hereby incorporated herein by reference: Ewing et al., Genome Res. 8:175 [1998]; Ewing & Green, Genome Res. 8:186 [1998]; Ewing et al., Genome Res. 8:195 [1998].

BLAST

The BLAST set of programs may be used to compare a set of sequences against databases composed of large numbers of nucleotide or protein sequences and obtain homologies to sequences with known function or properties. Detailed description of the BLAST software and its uses can be found in the following references which are hereby incorporated herein by reference: Altschul et al., J. Mol. Biol. 215:403 [1990]; Altschul et al., J. Mol. Biol. 219:555 [1991].

Generally, BLAST performs sequence similarity searching and is divided into 5 basic subroutines of which 3 were used: (1) BLASTN compares a nucleotide sequence to a nucleic acid sequence database; (2) BLASTX compares translated protein sequences from a nucleotide sequence done in six frames to a protein sequence database; (3) TBLASTX compares translated protein sequences from a nucleotide sequence done in six frames to the six frame translation of a nucleotide database. BLASTX and TBLASTX are used to identify homologies at the protein level of the nucleotide sequence.

B. Contig Sequence Assembly for Hits. Phred sequence calls and quality data for the individual sequencing runs associated with the above SEQ IDs are stored in a relational database. All the sequence runs stored in the database for the SEQ IDs to be assembled were extracted from the database and the files needed by Phrap recreated with the aid of a Perl script. Perl is an interpreted computer language useful for data manipulation. The same script ran Phrap on the assembled files and then stored the assembled contiguous sequences and singletons in a relational database. The script then assembled two files. One file was a FASTA format file of the sequences of the assembled contigs and singletons. The other file was a record of the assembled sequences and which sequencing runs they contained. FASTA format is a standard DNA sequence format recognized by the BLAST suite of programs as well as by Phrap. Both of these files were then inspected manually to detect incorrect assemblies or to add sequence information not present in the relational database. Any incorrect assemblies found were corrected before this file was used in BLAST searches to identify function and well as other homologous sequences in our databases. Correct assemblies that contained more than one SEQ ID were separated. Although these represent parts of the same sequence, since these are ESTs and contain limited gene sequence data, a one-to-one nucleotide match cannot be predicted at this time for the entire length of a contig representing a single SEQ ID with those containing multiple SEQ IDs. Some full length sequences were obtained and are designated with an FL.

C. Identification of Function. The FASTA formatted file obtained as described above was used to run a BLASTX query against the GenBank non-redundant protein database using a Perl script. The data from this analysis was parsed out by the Perl script such that the following information was extracted:the query sequence name, the level of homology to the hit and the description of the hit sequence (the highest scoring hit from the analysis). The script filtered all hits less than 1.00E-04, to eliminate spurious homologies. The data from this file was used to identify putative functions and properties for the query sequences

D. Identification of Similar Sequences in Derwent.TM.. The FASTA formatted file obtained as described above was used to run a BLASTN query against the Derwent.TM. non-redundant nucleotide database as well as a BLASTX against the Derwent.TM. non-redundant protein database using Perl scripts. These Derwent.TM. non-redundant databases were created by extracting all the sequence information in the Derwent.TM. database. The data from this analysis was parsed out by the Perl script such that the following information was extracted, the query sequence name, the level of homology to the hit and the description of the hit sequence (the highest scoring hit from the analysis). The script filtered all hits less than 1.00E-04, to eliminate spurious homologies.

E. Identification of Homologous Sequences. An internal relational database contains sequences from a large number of SEQ IDs belonging to a diverse group of organisms. In order to identify sequences in the database with high levels of homology to the sequences functionally identified as hits and contained in the FASTA formatted file described above, the following analysis was performed.

All the sequences were extracted in FASTA format from our relational database with standard SQL commands and converted into a searchable BLAST database using tools provided in the BLAST download from the National Center for Biotechnology Information (NCBI). A Perl script then ran a BLASTN search of our query file against our internal nucleotide database containing all relevant sequences. The script then extracted from all hits the following information: the query name, the level of homology and the hit seqID. The script then filtered all homologies less than 1.00E-04 as well as all the redundant seqIDs.

This analysis was repeated again using a TBLASTX query. Both files were then combined and the redundancies eliminated. Since the query sequences are also present in the database, these redundancies were manually eliminated from the results file. Lastly, all hit SEQ IDs with homology scores less than 1.00E-20 were filtered from the results list.

These results were used to extract the sequence and quality score data from the relational database in order to repeat the analysis described in "Contig Sequence Assembly for Hits". The final product consisted of two files. One file (FIG. 1) contains a record of the assembled sequences and the sequencing runs they include. The other file (FIG. 2) lists the search hits with homologies better than 1.00E-20 to the query contigs and singletons.

All publications and patents mentioned in the above specification are herein incorporated by reference. Various modifications and variations of the described compositions and methods of the invention will be apparent to those skilled in the art without departing from the scope and spirit of the invention. Although the invention has been described in connection with particular preferred embodiments, it should be understood that the inventions claimed should not be unduly limited to such specific embodiments. Indeed, various modifications of the described modes for carrying out the invention which are obvious to those skilled in the art and in fields related thereto are intended to be within the scope of the following claims.

SEQUENCE LISTINGS

0 SQTB SEQUENCE LISTING The patent contains a lengthy "Sequence Listing" section. A copy of the "Sequence Listing" is available in electronic form from the USPTO web site (http://seqdata.uspto.gov/?pageRequest=docDetail&DocID=US07635798B2)- . An electronic copy of the "Sequence Listing" will also be available from the USPTO upon request and payment of the fee set forth in 37 CFR 1.19(b)(3).