Reconstitution of EBV-directed T cell immunity by adoptive transfer of peptide-stimulated T cells in a patient after allogeneic stem cell transplantation for AITL.

Reconstitution of the T cell repertoire after allogeneic stem cell transplantation is a long and often incomplete process. As a result, reactivation of Epstein-Barr virus (EBV) is a frequent complication that may be treated by adoptive transfer of donor-derived EBV-specific T cells. We generated donor-derived EBV-specific T cells by stimulation with peptides representing defined epitopes covering multiple HLA restrictions. T cells were adoptively transferred to a patient who had developed persisting high titers of EBV after allogeneic stem cell transplantation for angioimmunoblastic T-cell lymphoma (AITL). T cell receptor beta (TCRß) deep sequencing showed that the T cell repertoire of the patient early after transplantation (day 60) was strongly reduced and only very low numbers of EBV-specific T cells were detectable. Manufacturing and in vitro expansion of donor-derived EBV-specific T cells resulted in enrichment of EBV epitope-specific, HLA-restricted T cells. Monitoring of T cell clonotypes at a molecular level after adoptive transfer revealed that the dominant TCR sequences from peptide-stimulated T cells persisted long-term and established an EBV-specific TCR clonotype repertoire in the host, with many of the EBV-specific TCRs present in the donor. This reconstituted repertoire was associated with immunological control of EBV and with lack of further AITL relapse.

T-cell repertoires in refractory coeliac disease.

OBJECTIVE: Refractory coeliac disease (RCD) is a potentially hazardous complication of coeliac disease (CD). In contrast to RCD type I, RCD type II is a precursor entity of enteropathy-associated T-cell lymphoma (EATL), which is associated with clonally expanding T-cells that are also found in the sequentially developing EATL. Using high-throughput sequencing (HTS), we aimed to establish the small-intestinal T-cell repertoire (TCR) in CD and RCD to unravel the role of distinct T-cell clonotypes in RCD pathogenesis. DESIGN: DNA extracted from duodenal mucosa specimens of controls (n=9), active coeliacs (n=10), coeliacs on a gluten-free diet (n=9), RCD type I (n=8), RCD type II (n=8) and unclassified Marsh I cases (n=3) collected from 2002 to 2013 was examined by TCRß-complementarity-determining regions 3 (CDR3) multiplex PCR followed by HTS of the amplicons. RESULTS: On average, 106 sequence reads per sample were generated consisting of up to 900 individual TCRß rearrangements. In RCD type II, the most frequent clonotypes (ie, sequence reads with identical CDR3) represent in average 42.6% of all TCRß rearrangements, which was significantly higher than in controls (6.8%; p<0.01) or RCD type I (6.7%; p<0.01). Repeat endoscopies in individual patients revealed stability of clonotypes for up to several years without clinical symptoms of EATL. Dominant clonotypes identified in individual patients with RCD type II were unique and not related between patients. CD-associated, gliadin-dependent CDR3 motifs were only detectable at low frequencies. CONCLUSIONS: TCRß-HTS analysis unravels the TCR in CD and allows detailed analysis of individual TCRß rearrangements. Dominant TCRß sequences identified in patients with RCD type II are unique and not homologous to known gliadin-specific TCR sequences, supporting the assumption that these clonal T-cells expand independent of gluten stimulation.

A new method to prevent carry-over contaminations in two-step PCR NGS library preparations.

Two-step PCR procedures are an efficient and well established way to generate amplicon libraries for NGS sequencing. However, there is a high risk of cross-contamination by carry-over of amplicons from first to second amplification rounds, potentially leading to severe misinterpretation of results. Here we describe a new method able to prevent and/or to identify carry-over contaminations by introducing the K-box, a series of three synergistically acting short sequence elements. Our K-boxes are composed of (i) K1 sequences for suppression of contaminations, (ii) K2 sequences for detection of possible residual contaminations and (iii) S sequences acting as separators to avoid amplification bias. In order to demonstrate the effectiveness of our method we analyzed two-step PCR NGS libraries derived from a multiplex PCR system for detection of T-cell receptor beta gene rearrangements. We used this system since it is of high clinical relevance and may be affected by very low amounts of contaminations. Spike-in contaminations are effectively blocked by the K-box even at high rates as demonstrated by ultra-deep sequencing of the amplicons. Thus, we recommend implementation of the K-box in two-step PCR-based NGS systems for research and diagnostic applications demanding high sensitivity and accuracy.

Molecular Mechanisms by Which a Fucus vesiculosus Extract Mediates Cell Cycle Inhibition and Cell Death in Pancreatic Cancer Cells.

Pancreatic cancer is one of the most aggressive cancer entities, with an extremely poor 5-year survival rate. Therefore, novel therapeutic agents with specific modes of action are urgently needed. Marine organisms represent a promising source to identify new pharmacologically active substances. Secondary metabolites derived from marine algae are of particular interest. The present work describes cellular and molecular mechanisms induced by an HPLC-fractionated, hydrophilic extract derived from the Baltic brown seaweed Fucus vesiculosus (Fv1). Treatment with Fv1 resulted in a strong inhibition of viability in various pancreatic cancer cell lines. This extract inhibited the cell cycle of proliferating cells due to the up-regulation of cell cycle inhibitors, shown on the mRNA (microarray data) and protein level. As a result, cells were dying in a caspase-independent manner. Experiments with non-dividing cells showed that proliferation is a prerequisite for the effectiveness of Fv1. Importantly, Fv1 showed low cytotoxic activity against non-malignant resting T cells and terminally differentiated cells like erythrocytes. Interestingly, accelerated killing effects were observed in combination with inhibitors of autophagy. Our in vitro data suggest that Fv1 may represent a promising new agent that deserves further development towards clinical application.

Mutations in DNAH5 cause primary ciliary dyskinesia and randomization of left-right asymmetry.

Primary ciliary dyskinesia (PCD, MIM 242650) is characterized by recurrent infections of the respiratory tract due to reduced mucociliary clearance and by sperm immobility. Half of the affected offspring have situs inversus (reversed organs), which results from randomization of left-right (LR) asymmetry. We previously localized to chromosome 5p a PCD locus containing DNAH5, which encodes a protein highly similar to the Chlamydomonas gamma-dynein heavy chain. Here we characterize the full-length 14-kb transcript of DNAH5. Sequence analysis in individuals with PCD with randomization of LR asymmetry identified mutations resulting in non-functional DNAH5 proteins.

Treatment with ribociclib shows favourable immunomodulatory effects in patients with hormone receptor-positive breast cancer-findings from the RIBECCA trial.

BACKGROUND: Inhibitors of the cyclin-dependent kinases 4 and 6 (CDK4/6i) have significantly improved clinical outcomes in patients with advanced hormone receptor-positive (HR+) breast cancer and have demonstrated favourable antitumour immune responses in preclinical studies. METHODS: Here, we investigated peripheral immune responses to ribociclib in patients with metastatic HR+ breast cancer as a preplanned exploratory subanalysis of the RIBECCA trial (NCT03096847). Peripheral blood mononuclear cells were subjected to immune cell profiling, gene expression analysis of immune-related signatures, and deep T cell receptor profiling before treatment started and after 12 weeks of treatment with ribociclib. RESULTS: Gene expression analysis revealed an upregulation of signatures associated with an activated adaptive immune system and a decrease in immunosuppressive cytokine signalling during treatment with ribociclib. Profiling of peripheral immune cell subpopulations showed a decrease in Treg cell frequencies, which was associated with treatment response. Furthermore, induction of CD4+ naive T cells could be seen, whereas effector and memory T cell populations remained largely unchanged. Correspondingly, T cell repertoire diversity remained mostly unchanged during treatment, although an increase in clonality could be observed in single patients. CONCLUSIONS: We show that treatment with ribociclib has significant effects on the peripheral innate and adaptive immune response in patients with HR+ breast cancer. Our data suggest that these effects lead to an activation of an already existing immune response rather than a de novo induction and make a strong case for future combination strategies of CDK4/6i with immunotherapies to enhance the adaptive immune response in HR+ breast cancer.

Transcript profiles in sugar beet genotypes uncover timing and strength of defense reactions to Cercospora beticola infection.

Cercospora leaf spot disease, caused by the fungus Cercospora beticola, is the most destructive foliar disease of sugar beet (Beta vulgaris) worldwide. Despite the great agronomical importance of this disease, little is known about its underlying molecular processes. Technical resources are scarce for analyzing this important crop species. We developed a sugar beet microarray with 44,000 oligonucleotides that represent 17,277 cDNAs. During the four stages of C. beticola-B. vulgaris interactions, we profiled the transcriptional responses of three genotypes: susceptible, polygenic partial resistance, and monogenic resistant. Similar genes were induced in all three genotypes during infection but with striking differences in timing. The monogenic resistant genotype displayed strong defense responses at 1 day postinoculation (dpi). The other genotypes displayed defense responses in a later phase (15 dpi) of the infection cycle. The partially resistant genotype displayed a strong defense response in the late phase of the infection cycle. Furthermore, the partially resistant genotype expressed pathogen-related transcripts that the susceptible genotype lacked. These results indicate that resistance was achieved by the ability to mount an early defense response, and partial resistance was determined by additional defense and signaling transcripts that allowed effective defense in the late phase of the infection cycle.

A second generation radiation hybrid map to aid the assembly of the bovine genome sequence.

BACKGROUND: Several approaches can be used to determine the order of loci on chromosomes and hence develop maps of the genome. However, all mapping approaches are prone to errors either arising from technical deficiencies or lack of statistical support to distinguish between alternative orders of loci. The accuracy of the genome maps could be improved, in principle, if information from different sources was combined to produce integrated maps. The publicly available bovine genomic sequence assembly with 6x coverage (Btau_2.0) is based on whole genome shotgun sequence data and limited mapping data however, it is recognised that this assembly is a draft that contains errors. Correcting the sequence assembly requires extensive additional mapping information to improve the reliability of the ordering of sequence scaffolds on chromosomes. The radiation hybrid (RH) map described here has been contributed to the international sequencing project to aid this process. RESULTS: An RH map for the 30 bovine chromosomes is presented. The map was built using the Roslin 3000-rad RH panel (BovGen RH map) and contains 3966 markers including 2473 new loci in addition to 262 amplified fragment-length polymorphisms (AFLP) and 1231 markers previously published with the first generation RH map. Sequences of the mapped loci were aligned with published bovine genome maps to identify inconsistencies. In addition to differences in the order of loci, several cases were observed where the chromosomal assignment of loci differed between maps. All the chromosome maps were aligned with the current 6x bovine assembly (Btau_2.0) and 2898 loci were unambiguously located in the bovine sequence. The order of loci on the RH map for BTA 5, 7, 16, 22, 25 and 29 differed substantially from the assembled bovine sequence. From the 2898 loci unambiguously identified in the bovine sequence assembly, 131 mapped to different chromosomes in the BovGen RH map. CONCLUSION: Alignment of the BovGen RH map with other published RH and genetic maps showed higher consistency in marker order and chromosome assignment than with the current 6x sequence assembly. This suggests that the bovine sequence assembly could be significantly improved by incorporating additional independent mapping information.

GOblet: a platform for Gene Ontology annotation of anonymous sequence data.

GOblet is a comprehensive web server application providing the annotation of anonymous sequence data with Gene Ontology (GO) terms. It uses a variety of different protein databases (human, murines, invertebrates, plants, sp-trembl) and their respective GO mappings. The user selects the appropriate database and alignment threshold and thereafter submits single or multiple nucleotide or protein sequences. Results are shown in different ways, e.g. as survey statistics for the main GO categories for all sequences or as detailed results for each single sequence that has been submitted. In its newest version, GOblet allows the batch submission of sequences and provides an improved display of results with the aid of Java applets. All output data, together with the Java applet, are packed to a downloadable archive for local installation and analysis. GOblet can be accessed freely at http://goblet.molgen.mpg.de.

Automated Gene Ontology annotation for anonymous sequence data.

Gene Ontology (GO) is the most widely accepted attempt to construct a unified and structured vocabulary for the description of genes and their products in any organism. Annotation by GO terms is performed in most of the current genome projects, which besides generality has the advantage of being very convenient for computer based classification methods. However, direct use of GO in small sequencing projects is not easy, especially for species not commonly represented in public databases. We present a software package (GOblet), which performs annotation based on GO terms for anonymous cDNA or protein sequences. It uses the species independent GO structure and vocabulary together with a series of protein databases collected from various sites, to perform a detailed GO annotation by sequence similarity searches. The sensitivity and the reference protein sets can be selected by the user. GOblet runs automatically and is available as a public service on our web server. The paper also addresses the reliability of automated GO annotations by using a reference set of more than 6000 human proteins. The GOblet server is accessible at http://goblet.molgen.mpg.de.

Mutational analysis of the NPHP4 gene in 250 patients with nephronophthisis.

Nephronophthisis (NPH), a recessive cystic kidney disease, is the most frequent genetic cause for end-stage renal disease in the first two decades of life. Mutations in three genes (NPHP1, 2, and 3) were identified as causative. Extrarenal manifestations are known, such as retinitis pigmentosa (Senior-Loken syndrome, SLS) and ocular motor apraxia type Cogan. Recently, we identified a novel gene (NPHP4) as mutated in NPH. To date, a total of only 13 different NPHP4 mutations have been described. To determine the frequency of NPHP4 mutations, we performed mutational analysis by direct sequencing of all 30 NPHP4 exons in 250 different patients with isolated NPH, SLS, or Cogan syndrome ascertained worldwide over 14 years. We identified 23 novel NPHP4 sequence variants in 26/250 different patients (10%). Interestingly, we detected homozygous or compound heterozygous mutations of NPHP4 in only 6/250 different patients (2.4%), but only one heterozygous NPHP4 sequence variant in 20/250 different patients (8%). In the six patients with two NPHP4 mutations, 5/8 mutations (63%) were likely loss-of-function mutations, whereas in the 20 patients with only one sequence variant, only 1/20 (5%) was a likely loss-of-function (i.e., truncating) mutation. We conclude that: i) two recessive mutations in NPHP4 are a rare cause of nephronophthisis; ii) single heterozygous NPHP4 sequence variants are three times more prevalent than two recessive mutations; iii) there is no genotype/phenotype correlation; iv) there must exist further genes causing nephronophthisis, since in 224/250 (90%) patients, no sequence variants in either of the four NPH genes were detected.

Genome-wide array analysis of normal and malformed human hearts.

BACKGROUND: We present the first genome-wide cDNA array analysis of human congenitally malformed hearts and attempted to partially elucidate these complex phenotypes. Most congenital heart defects, which account for the largest number of birth defects in humans, represent complex genetic disorders. As a consequence of the malformation, abnormal hemodynamic features occur and cause an adaptation process of the heart. METHODS AND RESULTS: The statistical analysis of our data suggests distinct gene expression profiles associated with tetralogy of Fallot, ventricular septal defect, and right ventricular hypertrophy. Applying correspondence analysis, we could associate specific gene functions to specific phenotypes. Furthermore, our study design allows the suggestion that alterations associated with primary genetic abnormalities can be distinguished from those associated with the adaptive response of the heart to the malformation (right ventricular pressure overload hypertrophy). We provide evidence for the molecular transition of the hypertrophic right ventricle to normal left ventricular characteristics. Furthermore, we present data on chamber-specific gene expression. CONCLUSIONS: Our findings propose that array analysis of malformed human hearts opens a new window to understand the complex genetic network of cardiac development and adaptation. For detailed access, see the online-only Data Supplement.

A first generation physical map of the medaka genome in BACs essential for positional cloning and clone-by-clone based genomic sequencing.

In order to realize the full potential of the medaka as a model system for developmental biology and genetics, characterized genomic resources need to be established, culminating in the sequence of the medaka genome. To facilitate the map-based cloning of genes underlying induced mutations and to provide templates for clone-based genomic sequencing, we have created a first-generation physical map of the medaka genome in bacterial artificial chromosome (BAC) clones. In particular, we exploited the synteny to the closely related genome of the pufferfish, Takifugu rubripes, by marker content mapping. As a first step, we clustered 103,144 public medaka EST sequences to obtain a set of 21,121 non-redundant sequence entities. Avoiding oversampling of gene-dense regions, 11,254 of EST clusters were successfully matched against the draft sequence of the fugu genome, and 2363 genes were selected for the BAC map project. We designed 35mer oligonucleotide probes from the selected genes and hybridized them against 64,500 BAC clones of strains Cab and Hd-rR, representing 14-fold coverage of the medaka genome. Our data set is further supplemented with 437 results generated from PCR-amplified inserts of medaka cDNA clones and BAC end-fragment markers. Our current, edited, first generation medaka BAC map consists of 902 map segments that cover about 74% of the medaka genome. The map contains 2721 markers. Of these, 2534 are from expressed sequences, equivalent to a non-redundant set of 2328 loci. The 934 markers (724 different) are anchored to the medaka genetic map. Thus, genetic map assignments provide immediate access to underlying clones and contigs, simplifying molecular access to candidate gene regions and their characterization.

Expression profiling of human idiopathic dilated cardiomyopathy.

OBJECTIVE: To investigate the global changes accompanying human dilated cardiomyopathy (DCM) we performed a large-scale expression screen using myocardial biopsies from a group of DCM patients with moderate heart failure. By hierarchical clustering and functional annotation of the deregulated genes we examined extensive changes in the cellular and molecular processes associated to DCM. METHODS: The expression profiles were obtained using a whole genome covering library (UniGene RZPD1) comprising 30336 cDNA clones and amplified RNA from myocardiac biopsies from 10 DCM patients in comparison to tissue samples from four non-failing, healthy donors. RESULTS: By setting stringent selection criteria 364 differentially expressed, sequence-verified non-redundant transcripts were identified with a false discovery rate of <0.001. Numerous genes and ESTs were identified representing previously recognised, as well as novel DCM-associated transcripts. Many of them were found to be upregulated and involved in cardiomyocyte energetics, muscle contraction or signalling. Two hundred and twenty-two deregulated transcripts were functionally annotated and hierarchically clustered providing an insight into the pathophysiology of DCM. Data was validated using the MLP-deficient mouse, in which several differentially expressed transcripts identified in the human DCM biopsies could be confirmed. CONCLUSIONS: We report the first genome-wide expression profile analysis using cardiac biopsies from DCM patients at various stages of the disease. Although there is a diversity of links between the cytoskeleton and the initiation of DCM, we speculate that genes implicated in intracellular signalling and in muscle contraction are associated with early stages of the disease. Altogether this study represents the most comprehensive and inclusive molecular portrait of human cardiomyopathy to date.

Influence of maternal low protein diet during pregnancy on hepatic gene expression signature in juvenile female porcine offspring.

SCOPE: Epidemiological and experimental evidence indicates that maternal nutrition status contributes to long-term changes in the metabolic phenotype of the offspring, a process known as fetal programming. METHODS AND RESULTS: We have used a swine model (Sus scrofa) to analyze consequences of a maternal low protein diet (about 50% of control) during pregnancy on hepatic lipid metabolism and genome-wide hepatic gene expression profile of juvenile female offspring (mean age 85 days). We found 318 S. scrofa genes to be differentially expressed in the liver at age 85 days. In the low protein offspring group key genes of fatty acid de novo synthesis were downregulated whereas several genes of lipolysis and phospholipid biosynthesis were upregulated. qRT-PCR analysis of selected genes verified microarray data and revealed linear correlations between gene expression levels and slaughter weight. Hepatic cholesterol 7α hydroxylase protein expression tended to be lower in the low protein group. Total lipid and triglyceride content and fatty acid composition of total lipids were not different between groups. CONCLUSION: A maternal low protein diet during pregnancy induces a distinct hepatic gene expression signature in juvenile female pigs which was not translated into phenotypical changes of liver lipid metabolism.

Effects on transcriptional regulation and lipid droplet characteristics in the liver of female juvenile pigs after early postnatal feed restriction and refeeding are dependent on birth weight.

Epidemiological and experimental data indicate that caloric restriction in early postnatal life may improve liver lipid metabolism in low birth weight individuals. The present study investigated transcriptional and metabolic responses to low (U) and normal (N) birth weight (d 75, T1) and postnatal feed restriction (R, 60% of controls, d 98, T2) followed by subsequent refeeding until d 131 of age (T3). Liver tissue studies were performed with a total of 42 female pigs which were born by multiparous German landrace sows. Overall, 194 genes were differentially expressed in the liver of U vs. N (T1) animals with roles in lipid metabolism. The total mean area and number of lipid droplets (LD) was about 4.6- and 3.7 times higher in U compared to N. In U, the mean LD size (µm(2)) was 24.9% higher. 3-week feed restriction reduced total mean area of LDs by 58.3 and 72.7% in U and N, respectively. A functional role of the affected genes in amino acid metabolism was additionally indicated. This was reflected by a 17.0% higher arginine concentration in the liver of UR animals (vs. NR). To evaluate persistency of effects, analyses were also done after refeeding period at T3. Overall, 4 and 22 genes show persistent regulation in U and N animals after 5 weeks of refeeding, respectively. These genes are involved in e.g. processes of lipid and protein metabolism and glucose homeostasis. Moreover, the recovery of total mean LD area in U and N animals back to the previous T1 level was observed. However, when compared to controls, the mean LD size was still reduced by 23.3% in UR, whereas it was increased in NR (+24.7%). The present results suggest that short-term postnatal feed restriction period programmed juvenile U animals for an increased rate of hepatic lipolysis in later life.

GOblet: annotation of anonymous sequence data with gene ontology and pathway terms.

The functional annotation of genomic data has become a major task for the ever-growing number of sequencing projects. In order to address this challenge, we recently developed GOblet, a free web service for the annotation of anonymous sequences with Gene Ontology (GO) terms. However, to overcome limitations of the GO terminology, and to aid in understanding not only single components but as well systemic interactions between the individual components, we have now extended the GOblet web service to integrate also pathway annotations. Furthermore, we extended and upgraded the data analysis pipeline with improved summaries, and added term enrichment and clustering algorithms. Finally, we are now making GOblet available as a stand-alone application for high-throughput processing on local machines. The advantages of this frequently requested feature is that a) the user can avoid restrictions of our web service for uploading and processing large amounts of data, and that b) confidential data can be analysed without insecure transfer to a public web server. The stand-alone version of the web service has been implemented using platform independent Tcl-scripts, which can be run with just a single runtime file utilizing the Starkit technology. The GOblet web service and the stand-alone application are freely available at http://goblet.molgen.mpg.de.

Comparative and evolutionary analysis of the rhesus macaque extended MHC class II region.

The sequence-based map of a part of the rhesus macaque major histocompatibility complex (MHC) extended class II region is presented. The sequenced region encompasses 67,401 bp and contains the SACM2L, RING1, FABGL and KE4 genes, as well as the HTATSF1-like and ZNF-like pseudogenes. Similar to human, but different from rat and mouse, no class I genes are found in the SACM2L- RING1 interval. The rhesus macaque extended MHC class II region shows a high degree of conservation of exonic as well as intronic and intergenic sequences compared with the respective human region. It is concluded that this particular genomic organization of the extended class II region-i.e., the absence of class I genes and the presence of the HTATSF1-like and ZNF-like pseudogenes-can be traced back to a common ancestor of humans and rhesus macaques about 23 million years ago.

New evidence for genome-wide duplications at the origin of vertebrates using an amphioxus gene set and completed animal genomes.

The 2R hypothesis predicting two genome duplications at the origin of vertebrates is highly controversial. Studies published so far include limited sequence data from organisms close to the hypothesized genome duplications. Through the comparison of a gene catalog from amphioxus, the closest living invertebrate relative of vertebrates, to 3453 single-copy genes orthologous between Caenorhabditis elegans (C), Drosophila melanogaster (D), and Saccharomyces cerevisiae (Y), and to Ciona intestinalis ESTs, mouse, and human genes, we show with a large number of genes that the gene duplication activity is significantly higher after the separation of amphioxus and the vertebrate lineages, which we estimate at 650 million years (Myr). The majority of human orthologs of 195 CDY groups that could be dated by the molecular clock appear to be duplicated between 300 and 680 Myr with a mean at 488 million years ago (Mya). We detected 485 duplicated chromosomal segments in the human genome containing CDY orthologs, 331 of which are found duplicated in the mouse genome and within regions syntenic between human and mouse, indicating that these were generated earlier than the human-mouse split. Model based calculations of the codon substitution rate of the human genes included in these segments agree with the molecular clock duplication time-scale prediction. Our results favor at least one large duplication event at the origin of vertebrates, followed by smaller scale duplication closer to the bird-mammalian split.

Construction of a 'unigene' cDNA clone set by oligonucleotide fingerprinting allows access to 25 000 potential sugar beet genes.

Access to the complete gene inventory of an organism is crucial to understanding physiological processes like development, differentiation, pathogenesis, or adaptation to the environment. Transcripts from many active genes are present at low copy numbers. Therefore, procedures that rely on random EST sequencing or on normalisation and subtraction methods have to produce massively redundant data to get access to low-abundance genes. Here, we present an improved oligonucleotide fingerprinting (ofp) approach to the genome of sugar beet (Beta vulgaris), a plant for which practically no molecular information has been available. To identify distinct genes and to provide a representative 'unigene' cDNA set for sugar beet, 159 936 cDNA clones were processed utilizing large-scale, high-throughput data generation and analysis methods. Data analysis yielded 30 444 ofp clusters reflecting the number of different genes in the original cDNA sample. A sample of 10 961 cDNA clones, each representing a different cluster, were selected for sequencing. Standard sequence analysis confirmed that 89% of these EST sequences did represent different genes. These results indicate that the full set of 30 444 ofp clusters represent up to 25 000 genes. We conclude that the ofp analysis pipeline is an accurate and effective way to construct large representative 'unigene' sets for any plant of interest with no requirement for prior molecular sequence data.