An STS-based map of the human genome

The Field-Dependent Nature of PageRank Values in Citation Networks

The value of scientific research can be easier to assess at the collective level than at the level of individual contributions. Several journal-level and article-level metrics aim to measure the importance of journals or individual manuscripts. However, many are citation-based and citation practices vary between fields. To account for these differences, scientists have devised normalization schemes to make metrics more comparable across fields. We use PageRank as an example metric and examine the extent to which field-specific citation norms drive estimated importance differences. In doing so, we recapitulate differences in journal and article PageRanks between fields. We also find that manuscripts shared between fields have different PageRanks depending on which field's citation network the metric is calculated in. We implement a degree-preserving graph shuffling algorithm to generate a null distribution of similar networks and find differences more likely attributed to field-specific preferences than citation norms. Our results suggest that while differences exist between fields' metric distributions, applying metrics in a field-aware manner rather than using normalized global metrics avoids losing important information about article preferences. They also imply that assigning a single importance value to a manuscript may not be a useful construct, as the importance of each manuscript varies by the reader's field.

Decoding psychosis: from national genome project to national brain project

The mind has puzzled humans for centuries, and its disorders, such as psychoses, have caused tremendous difficulties. However, relatively recent biotechnological breakthroughs, such as DNA technology and neuroimaging, have empowered scientists to explore the more fundamental aspects of psychosis. From searching for psychosis-causing genes to imaging the depths of the brain, scientists worldwide seek novel methods to understand the mind and the causes of its disorders. This article will briefly review the history of understanding and managing psychosis and the main findings of modern genetic research and then attempt to stimulate thought for decoding the biological mechanisms of psychosis in the present era of brain science.

Discovery of genomic variation across a generation

Over the past 30 years (the timespan of a generation), advances in genomics technologies have revealed tremendous and unexpected variation in the human genome and have provided increasingly accurate answers to long-standing questions of how much genetic variation exists in human populations and to what degree the DNA complement changes between parents and offspring. Tracking the characteristics of these inherited and spontaneous (or de novo) variations has been the basis of the study of human genetic disease. From genome-wide microarray and next-generation sequencing scans, we now know that each human genome contains over 3 million single nucleotide variants when compared with the ~ 3 billion base pairs in the human reference genome, along with roughly an order of magnitude more DNA—approximately 30 megabase pairs (Mb)—being ‘structurally variable’, mostly in the form of indels and copy number changes. Additional large-scale variations include balanced inversions (average of 18 Mb) and complex, difficult-to-resolve alterations. Collectively, ~1% of an individual’s genome will differ from the human reference sequence. When comparing across a generation, fewer than 100 new genetic variants are typically detected in the euchromatic portion of a child’s genome. Driven by increasingly higher-resolution and higher-throughput sequencing technologies, newer and more accurate databases of genetic variation (for instance, more comprehensive structural variation data and phasing of combinations of variants along chromosomes) of worldwide populations will emerge to underpin the next era of discovery in human molecular genetics.

Construction and integration of three de novo Japanese human genome assemblies toward a population-specific reference

The complete human genome sequence is used as a reference for next-generation sequencing analyses. However, some ethnic ancestries are under-represented in the reference genome (e.g., GRCh37) due to its bias toward European and African ancestries. Here, we perform de novo assembly of three Japanese male genomes using > 100× Pacific Biosciences long reads and Bionano Genomics optical maps per sample. We integrate the genomes using the major allele for consensus and anchor the scaffolds using genetic and radiation hybrid maps to reconstruct each chromosome. The resulting genome sequence, JG1, is contiguous, accurate, and carries the Japanese major allele at most loci. We adopt JG1 as the reference for confirmatory exome re-analyses of seven rare-disease Japanese families and find that re-analysis using JG1 reduces total candidate variant calls versus GRCh37 while retaining disease-causing variants. These results suggest that integrating multiple genomes from a single population can aid genome analyses of that population.

The Bermuda Triangle: The Pragmatics, Policies, and Principles for Data Sharing in the History of the Human Genome Project

The Bermuda Principles for DNA sequence data sharing are an enduring legacy of the Human Genome Project (HGP). They were adopted by the HGP at a strategy meeting in Bermuda in February of 1996 and implemented in formal policies by early 1998, mandating daily release of HGP-funded DNA sequences into the public domain. The idea of daily sharing, we argue, emanated directly from strategies for large, goal-directed molecular biology projects first tested within the "community" of C. elegans researchers, and were introduced and defended for the HGP by the nematode biologists John Sulston and Robert Waterston. In the C. elegans community, and subsequently in the HGP, daily sharing served the pragmatic goals of quality control and project coordination. Yet in the HGP human genome, we also argue, the Bermuda Principles addressed concerns about gene patents impeding scientific advancement, and were aspirational and flexible in implementation and justification. They endured as an archetype for how rapid data sharing could be realized and rationalized, and permitted adaptation to the needs of various scientific communities. Yet in addition to the support of Sulston and Waterston, their adoption also depended on the clout of administrators at the US National Institutes of Health (NIH) and the UK nonprofit charity the Wellcome Trust, which together funded 90% of the HGP human sequencing effort. The other nations wishing to remain in the HGP consortium had to accommodate to the Bermuda Principles, requiring exceptions from incompatible existing or pending data access policies for publicly funded research in Germany, Japan, and France. We begin this story in 1963, with the biologist Sydney Brenner's proposal for a nematode research program at the Laboratory of Molecular Biology (LMB) at the University of Cambridge. We continue through 2003, with the completion of the HGP human reference genome, and conclude with observations about policy and the historiography of molecular biology.

Sharing Data to Build a Medical Information Commons: From Bermuda to the Global Alliance

The Human Genome Project modeled its open science ethos on nematode biology, most famously through daily release of DNA sequence data based on the 1996 Bermuda Principles. That open science philosophy persists, but daily, unfettered release of data has had to adapt to constraints occasioned by the use of data from individual people, broader use of data not only by scientists but also by clinicians and individuals, the global reach of genomic applications and diverse national privacy and research ethics laws, and the rising prominence of a diverse commercial genomics sector. The Global Alliance for Genomics and Health was established to enable the data sharing that is essential for making meaning of genomic variation. Data-sharing policies and practices will continue to evolve as researchers, health professionals, and individuals strive to construct a global medical and scientific information commons.

Screening for single nucleotide polymorphisms in highly degraded DNA by using the amplified fragment length polymorphism technique

Short tandem repeat (STR) analysis is generally used for human identification of forensic samples; however, standard STR analysis sometimes fails to generate full profiles since DNA is frequently degraded by various environmental factors. Recently, single nucleotide polymorphism (SNP) analysis has attracted attention for human identification since the shorter amplicons are better suited for degraded samples. Though various SNP loci are used for analysis of degraded samples, it is unclear which ones are more appropriate. To identify SNPs that were resistant to degradation, we artificially degraded DNA obtained from the buccal swabs of six volunteers and the K562 cell line by heat treatment. Subsequently, the amplified fragment length polymorphism (AFLP) technique was used for SNP screening. We focused on the AFLP bands detected in both the heat-treated and untreated samples, and DNA extracted from these bands was directly sequenced. DNA degradation increased as the duration of heat treatment increased, and no STR profiles could be generated after 6h of heat treatment. When the AFLP band patterns were compared between 6h heat-treated and untreated samples, eight common bands were detected. The sequences of the DNA fragments of these common bands had higher adenine-thymine (A-T) content and included 17 SNPs. The SNPs detected in the heat-treated and untreated samples were considered to be resistant to degradation. Although there was a little information available in databases regarding the nine SNPs identified in this study, this study shows that some of these SNPs might be useful for human identification of extremely degraded DNA.

Unravelling the genetics of inherited retinal dystrophies: Past, present and future

The identification of the genes underlying monogenic diseases has been of interest to clinicians and scientists for many years. Using inherited retinal dystrophies as an example of monogenic disease we describe the history of molecular genetic techniques that have been pivotal in the discovery of disease causing genes. The methods that were developed in the 1970's and 80's are still in use today but have been refined and improved. These techniques enabled the concept of the Human Genome Project to be envisaged and ultimately realised. When the successful conclusion of the project was announced in 2003 many new tools and, as importantly, many collaborations had been developed that facilitated a rapid identification of disease genes. In the post-human genome project era advances in computing power and the clever use of the properties of DNA replication has allowed the development of next-generation sequencing technologies. These methods have revolutionised the identification of disease genes because for the first time there is no need to define the position of the gene in the genome. The use of next generation sequencing in a diagnostic setting has allowed many more patients with an inherited retinal dystrophy to obtain a molecular diagnosis for their disease. The identification of novel genes that have a role in the development or maintenance of retinal function is opening up avenues of research which will lead to the development of new pharmacological and gene therapy approaches. Neither of which can be used unless the defective gene and protein is known. The continued development of sequencing technologies also holds great promise for the advent of truly personalised medicine.

Radiation Hybrid Map of Barley Chromosome 3H

Assembly of the barley (Hordeum vulgare L.) genome is complicated by its large size (5.1 Gb) and proportion of repetitive elements (84%). This process is facilitated by high resolution maps for aligning bacterial artificial chromosome (BAC) contigs along chromosomes. Available genetic maps, however, do not provide accurate information on the physical position of a large portion of the genome located in recombination-poor regions. Radiation hybrid (RH) mapping is an alternative approach, which is based on radiation-induced deletions along the length of chromosomes. In this study, the first RH map for barley chromosome 3H was developed. In total, 373 in vivo RH lines were generated by irradiating wheat (Triticum aestivum L.)-barley chromosome 3H addition lines and crossing them to a normal wheat cultivar. Each RH informative line (containing deletions) had, on average, three deletions. The induced deletion size varied from 36.58 Kb to 576.00 Mb, with an average length of 52.42 Mb. This initial chromosome 3H radiation hybrid (3H-RH) map had a 9.53× higher resolution than an analogous genetic map, reaching a maximum of >262.40× resolution in regions around the centromere. The final RH map was 3066.1 cR in length, with a 0.76 Mb resolution. It was estimated that the map resolution can be improved to an average of 30.34 Kb by saturating the 3H-RH map with molecular markers. The generated RH panel enabled alignment of BAC and sequenced contigs as small as 1.50 Kb in size. The high resolution and the coverage of poor-recombination regions make RH maps an ideal resource for barley genome assembly, as well as other genetic studies.

Genetic differentiation of invasive Aedes albopictus by RAPD-PCR: implications for effective vector control

Aedes albopictus is an invasive mosquito species of great concern to public health as it is responsible for the biological transmission of several pathogens causing dengue, chikungunya, yellow fever, etc. In 2009, this mosquito was detected for the first time in Agra City. This study represents the first genetic analysis of A. albopictus from India. Random amplified polymorphic DNA (RAPD) was used to study the genetic structure of A. albopictus in four populations from different larval habitats. Seven RAPD primers produced 141 loci. The results displayed rich genetic variation among larval populations which is evident from high value of genetic differentiation (G ST), i.e. 0.280, indicating a very great genetic differentiation. Effective migration rates were observed to be 1.28, depicting a limited gene flow. According to analysis of molecular variance (AMOVA), the genetic distance between populations was significant (P < 0.05), showing a very high intrapopulation variation (96%) with only 4% variation among populations. Average genetic distances between populations were also calculated using PopGene software. Nei's average genetic distance between these populations was 0.112 (0.05-0.18). The cluster analysis technique of unweighted pair-group mean analysis (UPGMA) method of arithmetic averages was used to develop the phylogenetic tree which clearly shows two clusters of different larval habitats. The findings highlight high genetic differentiation indicating a slight migration rate confirming the recent introduction of this species in Agra region.

Clinical, molecular- and cytogenetic analysis of a case of severe radio-sensitivity

In radiotherapy the normal tissue reaction is often a limiting factor for radiation treatment. Still there is no screening method, which predicts normal tissue reaction on radiotherapy, especially in comparison to tumor tissue, and therefore allows tailoring of the radiation dose to each patient. Here, we present a case of severe radiation-related side effects. We applied classical cytogenetic techniques (Giemsa-banding and staining of centromeric regions), the comet assay as well as multicolor fluorescence in situ hybridization on peripheral blood lymphocytes of this patient in order to determine the radio-sensitivity on the DNA level and to correlate these findings with the clinical outcome. Our investigations revealed abnormalities on chromosome 9, deficiencies in the DNA-repair capacity after radiation exposure and a high number of radiation induced chromosomal aberrations. A detected high amount of residual damage two or three hours after radiation exposure and repair as well as the high number of chromosomal aberrations (ChAs) suggests a correlation between repair capacity and radiation induced ChAs. We concluded that the detected abnormalities might serve as a genetic basis for the radio-sensitive phenotype of this patient. Taken together this report strengthens the idea that intensive DNA genomic analysis of individual patients can serve as the basis for more favourable treatment of cancer patients.

Bioinformatic perspectives in the neuronal ceroid lipofuscinoses

The neuronal ceroid lipofuscinoses (NCLs) are a group of rare genetic diseases characterised clinically by the progressive deterioration of mental, motor and visual functions and histopathologically by the intracellular accumulation of autofluorescent lipopigment - ceroid - in affected tissues. The NCLs are clinically and genetically heterogeneous and more than 14 genetically distinct NCL subtypes have been described to date (CLN1-CLN14) (Haltia and Goebel, 2012 [1]). In this review we will chronologically summarise work which has led over the years to identification of NCL genes, and outline the potential of novel genomic techniques and related bioinformatic approaches for further genetic dissection and diagnosis of NCLs. This article is part of a Special Issue entitled: The Neuronal Ceroid Lipofuscinoses or Batten Disease.

Loss of TIMP-3 promotes tumor invasion via elevated IL-6 production and predicts poor survival and relapse in HPV-infected non-small cell lung cancer

Human papillomavirus (HPV) 16/18 E6 oncoprotein is expressed in lung tumors and is associated with p53 inactivation. The tissue inhibitor of metalloproteinase 3 (TIMP-3) is essential for limiting inflammation; therefore, we expected that TIMP-3 loss might induce chronic inflammation, thereby promoting tumor malignancy as well as poor survival and relapse in patients with HPV-infected non-small cell lung cancer. In this study, the loss of TIMP-3 by loss of heterozygosity and/or promoter hypermethylation was more frequent in HPV16/18 E6-positive tumors than in E6-negative tumors. To explore the possible underlying mechanism, E6-negative TL4 and CL1-0 cells were transfected with an E6 cDNA plasmid. A marked decrease in TIMP-3 expression was caused by promoter hypermethylation via increased DNA (cytosine-5-)-methyltransferase 1 (DNMT1) expression. Mechanistic studies indicated that TIMP-3 loss promoted interleukin-6 (IL-6) production, which led to cell invasion and anchorage-independent growth on soft agar plates. Kaplan-Meier and Cox regression models showed that patients with low-TIMP-3/high-IL-6 tumors had shorter overall survival and relapse-free survival periods when compared with patients with high-TIMP-3/low-IL-6 tumors. In summary, loss of TIMP-3 may increase IL-6 production via the tumor necrosis factor α/nuclear factor κB axis, thereby promoting tumor malignancy and subsequent relapse and poor survival in patients with HPV-infected non-small cell lung cancer.

A high-resolution 15,000(Rad) radiation hybrid panel for the domestic cat

The current genetic and recombination maps of the cat have fewer than 3,000 markers and a resolution limit greater than 1 Mb. To complement the first-generation domestic cat maps, support higher resolution mapping studies, and aid genome assembly in specific areas as well as in the whole genome, a 15,000(Rad) radiation hybrid (RH) panel for the domestic cat was generated. Fibroblasts from the female Abyssinian cat that was used to generate the cat genomic sequence were fused to a Chinese hamster cell line (A23), producing 150 hybrid lines. The clones were initially characterized using 39 short tandem repeats (STRs) and 1,536 SNP markers. The utility of whole-genome amplification in preserving and extending RH panel DNA was also tested using 10 STR markers; no significant difference in retention was observed. The resolution of the 15,000(Rad) RH panel was established by constructing framework maps across 10 different 1-Mb regions on different feline chromosomes. In these regions, 2-point analysis was used to estimate RH distances, which compared favorably with the estimation of physical distances. The study demonstrates that the 15,000(Rad) RH panel constitutes a powerful tool for constructing high-resolution maps, having an average resolution of 40.1 kb per marker across the ten 1-Mb regions. In addition, the RH panel will complement existing genomic resources for the domestic cat, aid in the accurate re-assemblies of the forthcoming cat genomic sequence, and support cross-species genomic comparisons.

Prevalence of the DNA repair enzyme-NEIL1 gene mutation in patients with type 2 diabetes in the Turkish population

In spite of the fact that an indirect relationship between NEIL1 gene and Type 2 diabetes has been demonstrated in animal model studies, there have been no human studies showing this relationship. In our study, we aimed to show the relationship between NEIL1 mutation and Type 2 diabetes in humans. The study group consisted 70 patients with Type 2 diabetes and the control group consisted of 50 healthy individuals. The mean age was 53±11 yr and 49±11 yr, respectively. Two NEIL1 mutations (2.9%) were detected in the patient group. There was A→G change (133A→G) at the 133. position of the 8th exon with 257 bp length in base sequencing. There was no mutation in the control group. We searched NEIL1 gene mutation for the first time in patients with Type 2 diabetes. This mutation was "silent" as it did not cause any amino acid change. The effects of these mutations on the etiopathogenesis of disease are not known. Although the lysine encoded by AAG was identical to the lysine encoded by AAA, it is not clear if they have functional differences due to the changing environmental conditions. NEIL1 gene mutation may have causative role in the development of Type 2 diabetes.

Progress and promise of genome-wide association studies for human complex trait genetics

Enormous progress in mapping complex traits in humans has been made in the last 5 yr. There has been early success for prevalent diseases with complex phenotypes. These studies have demonstrated clearly that, while complex traits differ in their underlying genetic architectures, for many common disorders the predominant pattern is that of many loci, individually with small effects on phenotype. For some traits, loci of large effect have been identified. For almost all complex traits studied in humans, the sum of the identified genetic effects comprises only a portion, generally less than half, of the estimated trait heritability. A variety of hypotheses have been proposed to explain why this might be the case, including untested rare variants, and gene-gene and gene-environment interaction. Effort is currently being directed toward implementation of novel analytic approaches and testing rare variants for association with complex traits using imputed variants from the publicly available 1000 Genomes Project resequencing data and from direct resequencing of clinical samples. Through integration with annotations and functional genomic data as well as by in vitro and in vivo experimentation, mapping studies continue to characterize functional variants associated with complex traits and address fundamental issues such as epistasis and pleiotropy. This review focuses primarily on the ways in which genome-wide association studies (GWASs) have revolutionized the field of human quantitative genetics.

Derivation, characterization, differentiation, and registration of seven human embryonic stem cell lines (VAL-3, -4, -5, -6M, -7, -8, and -9) on human feeder

Derivation of human embryonic stem cell lines has been a remarkable scientific achievement during the last decade. Human embryonic stem cells are regarded as an unlimited cell source for replacement therapy in regenerative medicine. Clearly, the scientific community requires proper derivation, characterization, and registration with the purpose of making them available for research and future medical applications worldwide. In this paper, we report our derivation work as the Valencian Node of the Spanish Stem Cell Bank in the generation, characterization, and registration of VAL-3, -4, -5, -6M, -7, -8, and 9 (www.isciii/htdocs/terapia/terapia_bancocelular.jsp). The derivation process was performed on microbiologically tested and irradiated human foreskin fibroblasts and designed to minimize contact with xeno-components in knockout Dulbecco's modified Eagle's medium supplemented with knockout serum replacement and basic fibroblast growth factor. Fingerprinting of the cell lines was performed to allow their identification and traceability. All lines were expressed at the mRNA and specific protein markers for undifferentiation and were found to be negative for classical differentiation markers such as neurofilament heavy chain (ectoderm), renin (mesoderm), and amylase (endoderm). All lines displayed high levels of telomerase activity and were shown to successfully overcome cryopreservation and thawing. Finally, we demonstrated the potential to differentiate in vitro (embryoid body formation) and in vivo (teratoma formation) into cell types from all three germ layers. Teratoma derived from all human embryonic stem cell lines present similar morphological features except VAL-8 that display more aggressive tumor behavior with a larger proportion of solid tissues, as opposed to cyst formation in the other cell lines.

Sequence assembly

Despite the rapidly increasing number of sequenced and re-sequenced genomes, many issues regarding the computational assembly of large-scale sequencing data have remain unresolved. Computational assembly is crucial in large genome projects as well for the evolving high-throughput technologies and plays an important role in processing the information generated by these methods. Here, we provide a comprehensive overview of the current publicly available sequence assembly programs. We describe the basic principles of computational assembly along with the main concerns, such as repetitive sequences in genomic DNA, highly expressed genes and alternative transcripts in EST sequences. We summarize existing comparisons of different assemblers and provide a detailed descriptions and directions for download of assembly programs at: http://genome.ku.dk/resources/assembly/methods.html.

Genetic mapping in human disease

Genetic mapping provides a powerful approach to identify genes and biological processes underlying any trait influenced by inheritance, including human diseases. We discuss the intellectual foundations of genetic mapping of Mendelian and complex traits in humans, examine lessons emerging from linkage analysis of Mendelian diseases and genome-wide association studies of common diseases, and discuss questions and challenges that lie ahead.

Application of radiation hybrid in gene mapping

Radiation hybrid (RH) mapping technique was exploited to determine chromosome locations of 26 human novel full length cDNAs recently cloned. All these cDNA clones were isolated from human cord blood CD (+) (34) cells and may be related to regulation of hematopoiesis. 23 genes were successfully mapped to chromosomal positions, while RH analyses were not possible in the remaining 3 cases. RH technique is indeed a powerful tool for mapping novel cDNA sequences due to its rapidity, precision, convenience and reproducibility.

Radiation hybrid (RH) and HAPPY mapping in plants

Radiation hybrid (RH) and HAPPY mapping are two technologies used in animal systems that have attracted the attention of the plant genetics community because they bridge the resolution gap between meiotic and BAC-based physical mapping that would facilitate the analysis of plant species lacking substantial genomics resources. Research has shown that the essence of these approaches can be applied and that a variety of strategies can be used to produce mapping panels. Mapping panels composed of live plants, protoplast fusion cultures, and sub-genomic DNA samples have been described. The resolution achievable by RH mapping panels involving live-plant derivatives of a monosomic maize (Zea mays) chromosome 9 addition in allohexaploid oat (Avena sativa), a monosomic chromosome 1D addition in allotetraploid durum wheat (Triticum turgidum), and interspecific hybrids between two tetraploid cotton species (G. hirsutum and G. barbadense), has been estimated to range from 0.6 to 6 Mb. On the other hand, a more comprehensive evaluation of one panel from durum wheat suggests that a higher mapping resolution (approximately 200 kb) is possible. In cases involving RH mapping panels based on barley (Hordeum vulgare)-tobacco (Nicotiana tabacum) protoplast fusions or a HAPPY mapping panel based on genomic DNA from Arabidopsis thaliana, the potential mapping resolution appears to be higher (50 to 200 kb). Despite these encouraging results, the application of either RH or HAPPY mapping in plants is still in the experimental phase and additional work is clearly needed before these methods are more routinely utilized.

Genomic DNA sequence comparison between two inbred soybean cyst nematode biotypes facilitated by massively parallel 454 micro-bead sequencing

Heterodera glycines, the soybean cyst nematode (SCN), is a damaging agricultural pest that could be effectively managed if critical phenotypes, such as virulence and host range could be understood. While SCN is amenable to genetic analysis, lack of DNA sequence data prevents the use of such methods to study this pathogen. Fortunately, new methods of DNA sequencing that produced large amounts of data and permit whole genome comparative analyses have become available. In this study, 400 million bases of genomic DNA sequence were collected from two inbred biotypes of SCN using 454 micro-bead DNA sequencing. Comparisons to a BAC, sequenced by Sanger sequencing, showed that the micro-bead sequences could identify low and high copy number regions within the BAC. Potential single nucleotide polymorphisms (SNPs) between the two SCN biotypes were identified by comparing the two sets of sequences. Selected resequencing revealed that up to 84% of the SNPs were correct. We conclude that the quality of the micro-bead sequence data was sufficient for de novo SNP identification and should be applicable to organisms with similar genome sizes and complexities. The SNPs identified will be an important starting point in associating phenotypes with specific regions of the SCN genome.

Constructing contigs from large-insert clones

This unit describes three approaches that are widely used to define alignments between overlapping clones bearing large-insert genomic DNA and to generate extensive contiguous overlapping sets of clones (contigs). The three approaches are sequence-tagged site (STS) content mapping, repetitive-element hybridization fingerprinting, and Alu-PCR fingerprinting. Methods for isolating the necessary BAC DNA suitable for automated fluorescent sequencing and generating new STS markers are discussed in support protocols. An alternate protocol presents repetitive-element hybridization fingerprinting to detect overlaps and build contigs with full-genomic YAC libraries.

Navigating public physical mapping databases

This unit provides concise overviews of the many physical mapping resources available and relates them to the genetic and transcript maps. Useful information on resolution of the maps, how to access them, and how to interpret them is compiled and presented in a clear fashion. Especially useful is a set of detailed protocols describing how to construct an STS marker and how to map it by means of available yeast artificial chromosomes (YACs). An additional protocol describes accessing EST marker maps.

Radiation hybrid maps of Medaka chromosomes LG 12, 17, and 22

The Medaka is an excellent genetic system for studies of vertebrate development and disease and environmental and evolutionary biology studies. To facilitate the mapping of markers or the cloning of affected genes in Medaka mutants identified by forward-genetic screens, we have established a panel of whole-genome radiation hybrids (RHs) and RH maps for three Medaka chromosomes. RH mapping is useful, since markers to be mapped need not be polymorphic and one can establish the order of markers that are difficult to resolve by genetic mapping owing to low genetic recombination rates. RHs were generated by fusing the irradiated donor, OLF-136 Medaka cell line, with the host B78 mouse melanoma cells. Of 290 initial RH clones, we selected 93 on the basis of high retention of fragments of the Medaka genome to establish a panel that allows genotyping in the 96-well format. RH maps for linkage groups 12, 17, and 22 were generated using 159 markers. The average retention for the three chromosomes was 19% and the average break point frequency was ∼33 kb/cR. We estimate the potential resolution of the RH panel to be ∼186 kb, which is high enough for integrating RH data with bacterial artificial chromosome clones. Thus, this first RH panel will be a useful tool for mapping mutated genes in Medaka.

Derivation and characterization of three new Spanish human embryonic stem cell lines (VAL -3 -4 -5) on human feeder and in serum-free conditions

A total of 184 human embryos, frozen for >5 years, were donated; informed consent was obtained according to Spanish law 45/2003. Survival rate was 40% and three out of 24 blastocysts (12.5%) developed into putative hESC lines, named VAL-3, VAL-4, and VAL-5. The derivation process was performed on microbiologically tested and irradiated human foreskin fibroblasts and designed to minimize contact with xeno-components in knockout DMEM supplemented with knockout serum replacement, and basic fibroblast growth factor. Fingerprinting and HLA typing of the cell lines allowed their identification and traceability. Karyotype was normal for VAL-3 (46XY), VAL-4 (46XX) and VAL-5 (46XX). All three hESC lines expressed specific markers for non-differentiation (Nanog, stage-specific embryonic antigen-4 [SSEA-4], tumour-related antigen [TRA]-1-60, and TRA-1-81) and were negative for SSEA-1. RT-PCR further demonstrated the expression of Oct-4, Sox2, Rex-1, Nanog, Cripto, Thy-1, and Lefty-A. Furthermore, they were found to be negative for classical differentiation markers such as neurofilament heavy chain (ectoderm), renin (mesoderm), and amylase (endoderm). All three cell lines displayed high levels of telomerase activity, and were shown to successfully overcome cryopreservation and thawing. Finally, these three new hESC lines have demonstrated the potential to differentiate in vitro and in vivo (teratoma formation) into cell types originating from all three germ layers.

Multiple upstream modules regulate zebrafish myf5 expression

Background

Myf5 is one member of the basic helix-loop-helix family of transcription factors, and it functions as a myogenic factor that is important for the specification and differentiation of muscle cells. The expression of myf5 is somite- and stage-dependent during embryogenesis through a delicate regulation. However, this complex regulatory mechanism of myf5 is not clearly understood.

Results

We isolated a 156-kb bacterial artificial chromosome clone that includes an upstream 80-kb region and a downstream 70-kb region of zebrafish myf5 and generated a transgenic line carrying this 156-kb segment fused to a green fluorescent protein (GFP) reporter gene. We find strong GFP expression in the most rostral somite and in the presomitic mesoderm during segmentation stages, similar to endogenous myf5 expression. Later, the GFP signals persist in caudal somites near the tail bud but are down-regulated in the older, rostral somites. During the pharyngula period, we detect GFP signals in pectoral fin buds, dorsal rostral myotomes, hypaxial myotomes, and inferior oblique and superior oblique muscles, a pattern that also corresponds well with endogenous myf5 transcripts. To characterize the specific upstream cis-elements that regulate this complex and dynamic expression pattern, we also generated several transgenic lines that harbor various lengths within the upstream 80-kb segment. We find that (1) the -80 kb/-9977 segment contains a fin and cranial muscle element and a notochord repressor; (2) the -9977/-6213 segment contains a strong repressive element that does not include the notochord-specific repressor; (3) the -6212/-2938 segment contains tissue-specific elements for bone and spinal cord; (4) the -2937/-291 segment contains an eye enhancer, and the -2937/-2457 segment is required for notochord and myocyte expression; and (5) the -290/-1 segment is responsible for basal transcription in somites and the presomitic mesoderm.

Conclusion

We suggest that the cell lineage-specific expression of myf5 is delicately orchestrated by multiple modules within the distal upstream region. This study provides an insight to understand the molecular control of myf5 and myogenesis in the zebrafish.

PeanutMap: an online genome database for comparative molecular maps of peanut

BACKGROUND

Molecular maps have been developed for many species, and are of particular importance for varietal development and comparative genomics. However, despite the existence of multiple sets of linkage maps, databases of these data are lacking for many species, including peanut.

DESCRIPTION

PeanutMap http://peanutgenetics.tamu.edu/cmap provides a web-based interface for viewing specific linkage groups of a map set. PeanutMap can display and compare multiple maps of a set based upon marker or trait correspondences, which is particularly important as cultivated peanut is a disomic tetraploid. The database can also compare linkage groups among multiple map sets, allowing identification of corresponding linkage groups from results of different research projects. Data from the two published peanut genome map sets, and also from three maps sets of phenotypic traits are present in the database. Data from PeanutMap have been incorporated into the Legume Information System website http://www.comparative-legumes.org to allow peanut map data to be used for cross-species comparisons.

CONCLUSION

The utility of the database is expected to increase as several SSR-based maps are being developed currently, and expanded efforts for comparative mapping of legumes are underway. Optimal use of these data will benefit from the development of tools to facilitate comparative analysis.

Genetic characterization of asymmetric somatic hybrids between Bupleurum scorzonerifolium Willd and Triticum aestivum L.: potential application to the study of the wheat genome

In this paper, we describe how Bupleurum scorzonerifolium/Triticum aestivum asymmetric somatic hybrids can be exploited to study the wheat genome. Protoplasts of B. scorzonerifolium Willd were irradiated with ultraviolet light (UV) and fused with protoplasts of common wheat (T. aestivum L.). All cell clones were similar in appearance to those of B. scorzonerifolium, while the regenerated plantlets were either intermediate or B. scorzonerifolium-like. Genotypic screening using isozymes showed that 39.3% of cell clones formed were hybrid. Some of the hybrid cell clones grew vigorously, and differentiated green leaves, shoots or plantlets. DNA marker analysis of the hybrids demonstrated that wheat DNA was integrated into the nuclear genomes of B. scorzonerifolium and in situ karyotyping cells revealed that a few wheat chromosome fragments had been introgressed into B. scorzonerifolium. The average wheat SSR retention frequency of the RH panel was 20.50%, but was only 6.67% in fusions with a non-irradiated donor. B. scorzonerifolium chromosomes and wheat SSR fragments in most asymmetric hybrid cell lines remained stable over a period of 2.5-3.5 years. We suggest the UV-induced asymmetric somatic hybrids between B. scorzonerifolium Willd and T. aestivum L. have the potential for use in the construction of an RH map of the wheat genome.

Intron-flanking EST-PCR markers: from genetic marker development to gene structure analysis in Rhododendron

With a long-term goal of constructing a linkage map of Rhododendron enriched with gene-specific markers, we utilized Rhododendron catawbiense ESTs for the development of high-efficiency (in terms of generating polymorphism frequency) PCR-based markers. Using the gene-sequence alignment between Rhododendron ESTs and the genomic sequences of Arabidopsis homologs, we developed 'intron-flanking' EST-PCR-based primers that would anneal in conserved exon regions and amplify across the more highly diverged introns. These primers resulted in increased efficiency (61% vs. 13%; 4.7-fold) of polymorphism-detection compared with conventional EST-PCR methods, supporting the assumption that intron regions are more diverged than exons. Significantly, this study demonstrates that Arabidopsis genome database can be useful in developing gene-specific PCR-based markers for other non-model plant species for which the EST data are available but genomic sequences are not. The comparative analysis of intron sizes between Rhododendron and Arabidopsis (made possible in this study by aligning of Rhododendron ESTs with Arabidopsis genomic sequences and the sequencing of Rhododendron genomic PCR products) provides the first insight into the gene structure of Rhododendron.

Holoprosencephaly and cleidocranial dysplasia in a patient due to two position-effect mutations: case report and review of the literature

Holoprosencephaly (HPE) is a genetically heterogeneous developmental field defect in which midline cleavage of the forebrain and craniofacial structures is impaired. Based on the analysis of HPE patients with chromosome rearrangements, at least six loci for the disorder have been assigned. The sonic hedgehog gene (SHH) at 7q36 has been identified as the HPE3 locus. Cleidocranial dysplasia (CCD) is an autosomal dominant skeletal disorder characterized by clavicular, pelvic and dental anomalies. It is caused by mutations in the osteoblast-specific transcription factor CBFA1/RUNX2, which maps to 6p21. We report a 20-year-old female with premaxillary agenesis (part of the HPE spectrum), as well as skeletal abnormalities and impacted teeth reminiscent of CCD. She carries a de novo 6;7 reciprocal translocation, with breakpoints at 6p21.1 and 7q36. We have shown previously that the 7q36 breakpoint maps 15 kb telomeric to the 5' end of SHH, which explains the patient's HPE phenotype. Now, using fluorescence in situ hybridization, we have identified a P1 artificial chromosome clone 800 kb upstream of CBFA1/RUNX2 that spans the 6p breakpoint. We propose that the proband's complex phenotype is due to two position-effect (PE) mutations, one at each translocation breakpoint, which have altered the expression of the SHH and CBFA1/RUNX2 genes. The role of PE mutations in human disease is also reviewed.

Radiation hybrid mapping of the species cytoplasm-specific (scsae) gene in wheat

Radiation hybrid (RH) mapping is based on radiation-induced chromosome breakage and analysis of chromosome segment retention or loss using molecular markers. In durum wheat (Triticum turgidum L., AABB), an alloplasmic durum line [(lo) durum] has been identified with chromosome 1D of T. aestivum L. (AABBDD) carrying the species cytoplasm-specific (scsae) gene. The chromosome 1D of this line segregates as a whole without recombination, precluding the use of conventional genome mapping. A radiation hybrid mapping population was developed from a hemizygous (lo) scsae--line using 35 krad gamma rays. The analysis of 87 individuals of this population with 39 molecular markers mapped on chromosome 1D revealed 88 radiation-induced breaks in this chromosome. This number of chromosome 1D breaks is eight times higher than the number of previously identified breaks and should result in a 10-fold increase in mapping resolution compared to what was previously possible. The analysis of molecular marker retention in our radiation hybrid mapping panel allowed the localization of scsae and 8 linked markers on the long arm of chromosome 1D. This constitutes the first report of using RH mapping to localize a gene in wheat and illustrates that this approach is feasible in a species with a large complex genome.

Identification of transcribed sequences in Arabidopsis thaliana by using high-resolution genome tiling arrays

Using a maskless photolithography method, we produced DNA oligonucleotide microarrays with probe sequences tiled throughout the genome of the plant Arabidopsis thaliana. RNA expression was determined for the complete nuclear, mitochondrial, and chloroplast genomes by tiling 5 million 36-mer probes. These probes were hybridized to labeled mRNA isolated from liquid grown T87 cells, an undifferentiated Arabidopsis cell culture line. Transcripts were detected from at least 60% of the nearly 26,330 annotated genes, which included 151 predicted genes that were not identified previously by a similar genome-wide hybridization study on four different cell lines. In comparison with previously published results with 25-mer tiling arrays produced by chromium masking-based photolithography technique, 36-mer oligonucleotide probes were found to be more useful in identifying intron-exon boundaries. Using two-dimensional HPLC tandem mass spectrometry, a small-scale proteomic analysis was performed with the same cells. A large amount of strongly hybridizing RNA was found in regions "antisense" to known genes. Similarity of antisense activities between the 25-mer and 36-mer data sets suggests that it is a reproducible and inherent property of the experiments. Transcription activities were also detected for many of the intergenic regions and the small RNAs, including tRNA, small nuclear RNA, small nucleolar RNA, and microRNA. Expression of tRNAs correlates with genome-wide amino acid usage.

Polymorphism profile of nine short tandem repeat Loci in the Han chinese

Nine short tandem repeat (STR) markers (D3S1358, VWA, FGA, THO1, TPOX, CSFIPO, D5S818, D13S317, and D7S820) and a sex-identification marker (Amelogenin locus) were amplified with multiplex PCR and were genotyped with a four-color fluorescence method in samples from 174 unrelated Han individuals in North China. The allele frequencies, genotype frequencies, heterozygosity, probability of discrimination powers, probability of paternity exclusion and Hardy-Weinberg equilibrium expectations were determined. The results demonstrated that the genotypes at all these STR loci in Han population conform to Hardy-Weinberg equilibrium expectations. The combined discrimination power (DP) was 1.05×10⁻¹⁰ within nine STR loci analyzed and the probability of paternity exclusion (EPP) was 0.9998. The results indicate that these nine STR loci and the Amelogenin locus are useful markers for human identification, paternity and maternity testing and sex determination in forensic sciences.

Cloneless genomic DNA analysis: an efficient and simple methods for de novo genomic sequencing projects and gap filling

The utility of using genomic DNA directly in agarose, i.e. cloneless libraries, in place of large clone libraries, radiation hybrid panels, or chromosome dissection was demonstrated. The advantage of the cloneless library approach is that, in principle, a targeted genomic resource can be developed rapidly for any genomic region using any genomic DNA sample. Here, a human chromosome 20 Not I fragment library was generated by slicing a pulsed field gel lane containing fractionating Not I cleaved DNA from a monosomic hybrid cell line into 2 mm pieces. A reliable PCR method using agarose embedded DNA was developed. InterAlu PCR generated unique patterns of products from adjacent slices (e.g. fractions). Further, the specificity of the interAlu products was demonstrated by FISH analysis and in other hybridization experiments to arrayed interAlu products. STS content mapping was used to order the fractions and also demonstrate the unique content of the library fractions.

High-resolution mapping of genotype-phenotype relationships in cri du chat syndrome using array comparative genomic hybridization

We have used array comparative genomic hybridization to map DNA copy-number changes in 94 patients with cri du chat syndrome who had been carefully evaluated for the presence of the characteristic cry, speech delay, facial dysmorphology, and level of mental retardation (MR). Most subjects had simple deletions involving 5p (67 terminal and 12 interstitial). Genotype-phenotype correlations localized the region associated with the cry to 1.5 Mb in distal 5p15.31, between bacterial artificial chromosomes (BACs) containing markers D5S2054 and D5S676; speech delay to 3.2 Mb in 5p15.32-15.33, between BACs containing D5S417 and D5S635; and the region associated with facial dysmorphology to 2.4 Mb in 5p15.2-15.31, between BACs containing D5S208 and D5S2887. These results overlap and refine those reported in previous publications. MR depended approximately on the 5p deletion size and location, but there were many cases in which the retardation was disproportionately severe, given the 5p deletion. All 15 of these cases, approximately two-thirds of the severely retarded patients, were found to have copy-number aberrations in addition to the 5p deletion. Restriction of consideration to patients with only 5p deletions clarified the effect of such deletions and suggested the presence of three regions, MRI-III, with differing effect on retardation. Deletions including MRI, a 1.2-Mb region overlapping the previously defined cri du chat critical region but not including MRII and MRIII, produced a moderate level of retardation. Deletions restricted to MRII, located just proximal to MRI, produced a milder level of retardation, whereas deletions restricted to the still-more proximal MRIII produced no discernible phenotype. However, MR increased as deletions that included MRI extended progressively into MRII and MRIII, and MR became profound when all three regions were deleted.

International Human Genome Sequencing Consortium. Finishing the euchromatic sequence of the human genome. Nature. 2004;431:931-945. [PMID: 15496913 DOI: 10.1038/nature03001]

The sequence of the human genome encodes the genetic instructions for human physiology, as well as rich information about human evolution. In 2001, the International Human Genome Sequencing Consortium reported a draft sequence of the euchromatic portion of the human genome. Since then, the international collaboration has worked to convert this draft into a genome sequence with high accuracy and nearly complete coverage. Here, we report the result of this finishing process. The current genome sequence (Build 35) contains 2.85 billion nucleotides interrupted by only 341 gaps. It covers approximately 99% of the euchromatic genome and is accurate to an error rate of approximately 1 event per 100,000 bases. Many of the remaining euchromatic gaps are associated with segmental duplications and will require focused work with new methods. The near-complete sequence, the first for a vertebrate, greatly improves the precision of biological analyses of the human genome including studies of gene number, birth and death. Notably, the human genome seems to encode only 20,000-25,000 protein-coding genes. The genome sequence reported here should serve as a firm foundation for biomedical research in the decades ahead.

Application of AFLP technology to radiation hybrid mapping

We have investigated the use of AFLP technology as a tool for the high throughput enrichment of Radiation Hybrid (RH) maps. The 3000 rad TM112 bovine RH panel was assayed with 37 EcoRI/TaqI AFLP primer combinations. The number of selective nucleotides used during PCR was increased to seven, to reduce the complexity of the AFLP profile and minimise the overlap between hamster and bovine bands co-amplified from hybrid cell clones. Seven-hundred-forty-seven bovine AFLP bands were amplified that could be distinguished following electrophoresis. Repeatability was tested within and between laboratories on independent template preparations and an error rate of 1.3% found. Two-point linkage analysis clustered 428 AFLP fragments in 39 linkage groups of at least 4 markers. Multi-point maps were constructed for 5 sample linkage groups. The study demonstrated that the AFLP approach could be used to rapidly screen for the most informative clones during panel construction and to increase the number of markers on RH maps, which could be useful for joining linkage groups formed by other markers. The use of AFLP markers as anchor points between existing RH maps and other physical maps, such as BAC contigs, is also discussed.

Analysis of myostatin gene structure, expression and function in zebrafish

Myostatin is a member of the TGF-beta family that functions as a negative regulator of skeletal muscle development and growth in mammals. Recently, Myostatin has also been identified in fish; however, its role in fish muscle development and growth remains unknown. We have reported here the isolation and characterization of myostatin genomic gene from zebrafish and analysis of its expression in zebrafish embryos, larvae and adult skeletal muscles. Our data showed that myostatin was weakly expressed in early stage zebrafish embryos, and strongly expressed in swimming larvae, juvenile and skeletal muscles of adult zebrafish. Transient expression analysis revealed that the 1.2 kb zebrafish myostatin 5' flanking sequence could direct green fluorescent protein (GFP) expression predominantly in muscle cells, suggesting that the myostatin 5' flanking sequence contained regulatory elements required for muscle expression. To determine the biological function of Myostatin in fish, we generated a transgenic line that overexpresses the Myostatin prodomain in zebrafish skeletal muscles using a muscle-specific promoter. The Myostatin prodomain could act as a dominant negative and inhibit Myostatin function in skeletal muscles. Transgenic zebrafish expressing the Myostatin prodomain exhibited no significant change in myogenic gene expression and differentiation of slow and fast muscle cells at their embryonic stage. The transgenic fish, however, exhibited an increased number of myofibers in skeletal muscles, but no significant difference in fiber size. Together, these data demonstrate that Myostatin plays an inhibitory role in hyperplastic muscle growth in zebrafish.

Integrated and sequence-ordered BAC- and YAC-based physical maps for the rat genome

As part of the effort to sequence the genome of Rattus norvegicus, we constructed a physical map comprised of fingerprinted bacterial artificial chromosome (BAC) clones from the CHORI-230 BAC library. These BAC clones provide approximately 13-fold redundant coverage of the genome and have been assembled into 376 fingerprint contigs. A yeast artificial chromosome (YAC) map was also constructed and aligned with the BAC map via fingerprinted BAC and P1 artificial chromosome clones (PACs) sharing interspersed repetitive sequence markers with the YAC-based physical map. We have annotated 95% of the fingerprint map clones in contigs with coordinates on the version 3.1 rat genome sequence assembly, using BAC-end sequences and in silico mapping methods. These coordinates have allowed anchoring 358 of the 376 fingerprint map contigs onto the sequence assembly. Of these, 324 contigs are anchored to rat genome sequences localized to chromosomes, and 34 contigs are anchored to unlocalized portions of the rat sequence assembly. The remaining 18 contigs, containing 54 clones, still require placement. The fingerprint map is a high-resolution integrative data resource that provides genome-ordered associations among BAC, YAC, and PAC clones and the assembled sequence of the rat genome.

Reading maps of the genes: interpreting the spatiality of genetic knowledge

Genetics has become the pre-eminent interpretation of the body and health and illness. This paper engages with a central technique and metaphor of the new genetics-gene mapping. Through an exploration of the process of gene mapping, the paper argues that the genetic material of the body is spatialised and transformed into a knowable and manipulable entity. Three interpretations of this spatial transformation of the body's materiality are discussed, in turn drawing on Foucault's notion of the construction of medical knowledge, the deconstruction of geographical maps and Haraway's 'fetishised' conception of the gene map. The paper concludes by considering contestations to this dominant discourse, and begins to construct an alternative spatialisation of the body that attempts to 'place' the gene more appropriately in a socially-embedded body and health.

Genetic and cytogenetic analyses of breast cancer yield different perspectives of a complex disease

Genomic instability in breast cancer results in low-level changes in DNA copy number, a significant but poorly understood mechanism underlying the genetic heterogeneity of this disorder. Two different approaches, loss of heterozygosity (LOH) and comparative genomic hybridization (CGH), have been used to probe the genetics of breast cancer evolution. LOH is a locus specific method that detects the variation in the parental origin of DNA, but is not quantitative. CGH provides a genome-wide accounting of the magnitude of DNA copy number changes, but not parental origin. Both methods have identified complex and heterogeneous patterns of DNA losses, duplications, and amplifications during breast cancer evolution. LOH and CGH technologies interrogate very distinct mechanisms driving breast tumor evolution, yet are seldom used in parallel to profile specimens. Thus, the relative significance of genetic versus numerical variations of DNA in breast cancer evolution remains undefined. This review will attempt to summarize some of the successes of these investigations, highlight some complex and confounding observations emerging from these studies, and discuss the potential of these studies to improve our understanding of breast cancer biology and treatment.

GALA, a database for genomic sequence alignments and annotations

We have developed a relational database to contain whole genome sequence alignments between human and mouse with extensive annotations of the human sequence. Complex queries are supported on recorded features, both directly and on proximity among them. Searches can reveal a wide variety of relationships, such as finding all genes expressed in a designated tissue that have a highly conserved noncoding sequence 5' to the start site. Other examples are finding single nucleotide polymorphisms that occur in conserved noncoding regions upstream of genes and identifying CpG islands that overlap the 5' ends of divergently transcribed genes. The database is available online at http://globin.cse.psu.edu/ and http://bio.cse.psu.edu/.

A preliminary linkage map of the hard tick, Ixodes scapularis

A linkage map of the Ixodes scapularis genome was constructed, based upon segregation amongst 127 loci. These included 84 random amplified polymorphic DNA (RAPD) markers, 32 Sequence-Tagged RAPD (STAR) markers, 5 cDNAs, and 5 microsatellites in 232 F1 intercross progeny from a single, field-collected P1 female. A preliminary linkage map of 616 cM was generated across 14 linkage groups with one marker every 10.8 cM. Assuming a genome size of approximately 109 bp, the relationship of physical to genetic distance was found to be approximately 300 kb/cM in the I. scapularis genome.

A novel algorithm for computational identification of contaminated EST libraries

A key goal of the Human Genome Project was to understand the complete set of human proteins, the proteome. Since the genome sequence by itself is not sufficient for predicting new genes and alternative splicing events that lead to new proteins, expressed sequence tags (ESTs) are used as the primary tool for these purposes. The high prevalence of artifacts in dbEST, however, often leads to invalid predictions. Here we describe a novel method for recognizing genomic DNA contamination and other artifacts that cannot be identified using current EST cleaning techniques. Our method uses the alignment of the entire set of ESTs to the human genome to identify highly contaminated EST libraries. We discovered 53 highly contaminated libraries and a subset of 24 766 ESTs from these libraries that probably represent contamination with genomic DNA, pre-mRNA, and ESTs that span non-canonical introns. Although this is only a small fraction of the entire EST dataset, each contaminating sequence could create a spurious transcript prediction. Indeed, in the clustering and assembly tool that we used, these sequences would have caused incorrect inference of 9575 new splice variants and 6370 new genes. Conclusions based on EST analysis, including prediction of alternative splicing, should be re-evaluated in light of these results. Our method, along with the identified set of contaminated sequences, will be essential for applications that depend on large EST datasets.

GenBank

GenBank (R) is a comprehensive sequence database that contains publicly available DNA sequences for more than 119 000 different organisms, obtained primarily through the submission of sequence data from individual laboratories and batch submissions from large-scale sequencing projects. Most submissions are made using the BankIt (web) or Sequin programs and accession numbers are assigned by GenBank staff upon receipt. Daily data exchange with the EMBL Data Library in the UK and the DNA Data Bank of Japan helps ensure worldwide coverage. GenBank is accessible through NCBI's retrieval system, Entrez, which integrates data from the major DNA and protein sequence databases along with taxonomy, genome, mapping, protein structure and domain information, and the biomedical journal literature via PubMed. BLAST provides sequence similarity searches of GenBank and other sequence databases. Complete bimonthly releases and daily updates of the GenBank database are available by FTP. To access GenBank and its related retrieval and analysis services, go to the NCBI home page at: http://www.ncbi.nlm.nih.gov.