Protect Our Kids: a novel program bringing hemorrhage control to schools

Background

Following the shooting at Sandy Hook Elementary School, the Hartford Consensus produced the Stop the Bleed program to train bystanders in hemorrhage control. In our region, the police bureau delivers critical incident training to public schools, offering instruction in responding to violent or dangerous situations. Until now, widespread training in hemorrhage control has been lacking. Our group developed, implemented and evaluated a novel program integrating hemorrhage control into critical incident training for school staff in order to blunt the impact of mass casualty events on children.

Methods

The staff of 25 elementary and middle schools attended a 90-minute course incorporating Stop the Bleed into the critical incident training curriculum, delivered on-site by police officers, nurses and doctors over a three-day period. The joint program was named Protect Our Kids. At the conclusion of the course, hemorrhage control kits and educational materials were provided and a four-question survey to assess the quality of training using a ten-point Likert scale was completed by participants and trainers.

Results

One thousand eighteen educators underwent training. A majority were teachers (78.2%), followed by para-educators (5.8%), counselors (4.4%) and principals (2%). Widely covered by local and state media, the Protect Our Kids program was rated as excellent and effective by a majority of trainees and all trainers rated the program as excellent.

Conclusions

Through collaboration between trauma centers, police and school systems, a large-scale training program for hemorrhage control and critical incident response can be effectively delivered to schools.

Section 7: Bioinformatics: Computational Approaches to Analysis of DNA Microarray Data

To review the current state of the art in computational methods for the analysis of DNA microarray data.

The review considers methods of microarray data collection, transformation and representation, comparisons and predictions of gene expression from the data, their mechanistic analysis, related systems biology, and the application of clustering techniques.

Functional genomics approaches have greatly increased the rate at which data on biological systems is generated, leading to corresponding challenges in analyzing the data through advanced computational techniques . The paper compares and contrasts the application of computational clustering for discovery, comparison, and prediction of gene expression classes, together with their evaluation and relation to mechanistic analyses of biological systems.

Methods for assaying gene expression levels by DNA microarray experiments produce considerably more data than other techniques, and require a wide variety of computational techniques for identifying patterns of expression that may be biologically significant. These will have to be verified and validated by comparison to results from other methods, integrated with other systems data, and provide the feedback for further experimentation for testing mechanistic or other biological hypotheses.

Abstract P2-09-09: Polyligand profiling differentiates cancer patients according to their benefit of treatment

Introduction: Deconvolution of multi-nodal perturbations in cancer network architecture demands highly multiplexed profiling assays. We demonstrate the value of polyligand profiling of tumor systems states using libraries of single stranded oligodeoxynucleotides (ssODN) to distinguish between tumor tissue from breast cancer patients who did or did not derive benefit from treatment regimens containing trastuzumab.

Methods: This study included cases from women with invasive breast cancer who received chemotherapy+ trastuzumab (C+T) or trastuzumab monotherapy with available retrospective data on the time to next treatment (TTNT). A library of 2x1012 unique ssODN was exposed to FFPE tissues from patients who benefited (B) or not (NB) from trastuzumab-based regimens in several rounds of positive and negative selection. Two enriched libraries were screened on independent set of 42 B and 19 NB cases using a modified IHC protocol for detection of bound ssODNs. Poly-Ligand Profiles (PLP) were scored by a blinded pathologist. Two libraries, EL-NB and EL-B, showed significant p-values between groups of responders and non-responders. A Cox-PH model was fitted using either tumors' HER2 status or PLP test results as the independent variable. Median survival time was calculated from the Kaplan-Meier estimate. A separate group of 63 cases with TTNT data from chemotherapy without trastuzumab was used as a control to distinguish prognostic from predictive performance.

Results: The PLP scores of EL-NB and EL-B were assessed by receiver operating characteristic (ROC) curves and resulted in a combined AUC value of 0.81. EL-NB and EL-B were able to effectively classify B and NB patients with either HER2-negative/equivocal (AUC = 0.73) or HER2-positive cancers (AUC = 0.84). In contrast, HER2 status alone yielded an AUC value of 0.47. The combined PLP scores for the independent set of 63 patients treated with C excluding trastuzumab resulted in an AUC value of 0.53, indicating that the assay was predictive and not simply prognostic. Kaplan-Meier curves analysis shows that PLP+ cases have 429 days median TTNT, while PLP- cases have 129 days (HR = 0.38, log-rank p = 0.001). Analysis based on HER2 status showed no significant difference in TTNT between patients that were HER2+ (280 days) or HER2-negative/equivocal (336 days, HR = 1.27, log-rank p =0.45).

Summary: Performance of the PLP assay in differentiating patients who did or did not benefit from trastuzumab therapy outperforms the standard IHC assay for HER2 status. These results represent a promising step towards the development of a CDx to identify the 50-70% of HER2+ patients who will not benefit from trastuzumab. In addition, PLP also has the potential to identify the HER2-negative/equivocal patients who may benefit from trastuzumab-containing regimens.

Citation Format: Domenyuk V, Gatalica Z, Santhanam R, Wei X, Stark A, Kennedy P, Toussaint B, Levenberg S, Wang R, Xiao N, Greil R, Rinnerthaler G, Gampenrieder S, Heimberger AB, Berry DJ, Barker A, Demetri GD, Quackenbush J, Marshall JL, Poste G, Vacirca JL, Vidal GA, Schwartzberg LS, Halbert DD, Voss A, Miglarese MR, Famulok M, Mayer G, Spetzler D. Polyligand profiling differentiates cancer patients according to their benefit of treatment [abstract]. In: Proceedings of the 2017 San Antonio Breast Cancer Symposium; 2017 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2018;78(4 Suppl):Abstract nr P2-09-09.

From 'omes to biology

Technologies that have emerged from the genome project have dramatically increased our ability to generate data on the way in which organisms respond to their environments, how they execute their programmes of development and growth, and how these are altered in the development of disease states. However, our ability to analyse these large datasets has not kept pace with our ability to generate them and consequently new strategies must be developed to address the issues associated with their analysis. One approach that we have employed quite successfully is to look at data from microarrays (or proteomics or metabolomics experiments) not as independent datasets, but rather as elements of a much larger body of biological information across various scales that must be integrated with, and interpreted within, the context of such ancillary data. Here we outline the general approach and provide three examples from published studies of the way in which we have applied this strategy.

Computational approaches to analysis of DNA microarray data

OBJECTIVES

To review the current state of the art in computational methods for the analysis of DNA microarray data.

METHODS

The review considers methods of microarray data collection, transformation and representation, comparisons and predictions of gene expression from the data, their mechanistic analysis, related systems biology, and the application of clustering techniques.

RESULTS

Functional genomics approaches have greatly increased the rate at which data on biological systems is generated, leading to corresponding challenges in analyzing the data through advanced computational techniques. The paper compares and contrasts the application of computational clustering for discovery, comparison, and prediction of gene expression classes, together with their evaluation and relation to mechanistic analyses of biological systems.

CONCLUSION

Methods for assaying gene expression levels by DNA microarray experiments produce considerably more data than other techniques, and require a wide variety of computational techniques for identifying patterns of expression that may be biologically significant. These will have to be verified and validated by comparison to results from other methods, integrated with other systems data, and provide the feedback for further experimentation for testing mechanistic or other biological hypotheses.

The transcriptional landscape of the mammalian genome

This study describes comprehensive polling of transcription start and termination sites and analysis of previously unidentified full-length complementary DNAs derived from the mouse genome. We identify the 5' and 3' boundaries of 181,047 transcripts with extensive variation in transcripts arising from alternative promoter usage, splicing, and polyadenylation. There are 16,247 new mouse protein-coding transcripts, including 5154 encoding previously unidentified proteins. Genomic mapping of the transcriptome reveals transcriptional forests, with overlapping transcription on both strands, separated by deserts in which few transcripts are observed. The data provide a comprehensive platform for the comparative analysis of mammalian transcriptional regulation in differentiation and development.

BmiGI: a database of cDNAs expressed in Boophilus microplus, the tropical/southern cattle tick

We used an expressed sequence tag approach to initiate a study of the genome of the southern cattle tick, Boophilus microplus. A normalized cDNA library was synthesized from pooled RNA purified from tick larvae which had been subjected to different treatments, including acaricide exposure, heat shock, cold shock, host odor, and infection with Babesia bovis. For the acaricide exposure experiments, we used several strains of ticks, which varied in their levels of susceptibility to pyrethroid, organophosphate and amitraz. We also included RNA purified from samples of eggs, nymphs and adult ticks and dissected tick organs. Plasmid DNA was prepared from 11,520 cDNA clones and both 5' and 3' sequencing performed on each clone. The sequence data was used to search public protein databases and a B. microplus gene index was constructed, consisting of 8270 unique sequences whose associated putative functional assignments, when available, can be viewed at the TIGR website (http://www.tigr.org/tdb/tgi). A number of novel sequences were identified which possessed significant sequence similarity to genes, which might be involved in resistance to acaricides.

Grape-seed derived procyanidins interfere with adipogenesis of 3T3-L1 cells at the onset of differentiation

OBJECTIVE

Our group's previous results on the effects of a grape seed procyanidin extract (GSPE) on adipose metabolism showed that peroxisome proliferator-activated receptor-gamma (PPARgamma) plays a central role in the lipolytic effects of GSPE on adipocytes. Since PPARgamma2 is a main regulator of the differentiation process of adipocytes, we investigated whether GSPE affects the adipogenesis of 3T3-L1 cells.

DESIGN

We performed a time point screening by treating 3T3-L1 cells with GSPE during the differentiation process for 24 h.

MEASUREMENTS

Differentiation markers and differential gene expression due to GSPE treatment (using the microarray technique).

RESULTS

Twenty four hour-GSPE treatment at the onset of differentiation reduces adipose-specific markers and maintains the expression of preadipocyte marker preadipocyte factor-1 (Pref-1) significantly elevated. These effects were not found in other time points. Microarray analysis of gene expression after GSPE treatment at the early stage of differentiation showed a modified gene expression profile in which cell cycle and growth-related genes were downregulated by GSPE.

CONCLUSION

These results suggest that GSPE affects adipogenesis, mainly at the induction of differentiation, and that procyanidins may have a new role in which they impede the formation of adipose cells.

The TIGR Gene Indices: clustering and assembling EST and known genes and integration with eukaryotic genomes

Although the list of completed genome sequencing projects has expanded rapidly, sequencing and analysis of expressed sequence tags (ESTs) remain a primary tool for discovery of novel genes in many eukaryotes and a key element in genome annotation. The TIGR Gene Indices (http://www.tigr.org/tdb/tgi) are a collection of 77 species-specific databases that use a highly refined protocol to analyze gene and EST sequences in an attempt to identify and characterize expressed transcripts and to present them on the Web in a user-friendly, consistent fashion. A Gene Index database is constructed for each selected organism by first clustering, then assembling EST and annotated cDNA and gene sequences from GenBank. This process produces a set of unique, high-fidelity virtual transcripts, or tentative consensus (TC) sequences. The TC sequences can be used to provide putative genes with functional annotation, to link the transcripts to genetic and physical maps, to provide links to orthologous and paralogous genes, and as a resource for comparative and functional genomic analysis.

Sequence analysis of a rainbow trout cDNA library and creation of a gene index

Expressed sequence tag (EST) projects have produced extremely valuable resources for identifying genes affecting phenotypes of interest. A large-scale EST sequencing project for rainbow trout was initiated to identify and functionally annotate as many unique transcripts as possible. Over 45,000 5' ESTs were obtained by sequencing clones from a single normalized library constructed using mRNA from six tissues. The production of this sequence data and creation of a rainbow trout Gene Index eliminating redundancy and providing annotation for these sequences will facilitate research in this species.

Enrichment of gene-coding sequences in maize by genome filtration

Approximately 80% of the maize genome comprises highly repetitive sequences interspersed with single-copy, gene-rich sequences, and standard genome sequencing strategies are not readily adaptable to this type of genome. Methodologies that enrich for genic sequences might more rapidly generate useful results from complex genomes. Equivalent numbers of clones from maize selected by techniques called methylation filtering and High C0t selection were sequenced to generate approximately 200,000 reads (approximately 132 megabases), which were assembled into contigs. Combination of the two techniques resulted in a sixfold reduction in the effective genome size and a fourfold increase in the gene identification rate in comparison to a nonenriched library.

Analysis of the mouse transcriptome based on functional annotation of 60,770 full-length cDNAs

Only a small proportion of the mouse genome is transcribed into mature messenger RNA transcripts. There is an international collaborative effort to identify all full-length mRNA transcripts from the mouse, and to ensure that each is represented in a physical collection of clones. Here we report the manual annotation of 60,770 full-length mouse complementary DNA sequences. These are clustered into 33,409 'transcriptional units', contributing 90.1% of a newly established mouse transcriptome database. Of these transcriptional units, 4,258 are new protein-coding and 11,665 are new non-coding messages, indicating that non-coding RNA is a major component of the transcriptome. 41% of all transcriptional units showed evidence of alternative splicing. In protein-coding transcripts, 79% of splice variations altered the protein product. Whole-transcriptome analyses resulted in the identification of 2,431 sense-antisense pairs. The present work, completely supported by physical clones, provides the most comprehensive survey of a mammalian transcriptome so far, and is a valuable resource for functional genomics.

Use of RNA and genomic DNA references for inferred comparisons in DNA microarray analyses

In most microarray assays, labeled cDNA molecules derived from reference and query RNA samples are co-hybridized to probes arrayed on a glass surface. Gene expression profiles are then calculated for each gene based on the relative hybridization intensities measured between the two samples. The most commonly used reference samples are typically isolates from a single representative RNA source (RNA-0) or pooled mixtures of RNA derived from a plurality of sources (RNA-p). Genomic DNA offers an alternative reference nucleic acid with a number of potential advantages, including stability, reproducibility, and a potentially uniform representation of all genes, as each unique gene should have equal representation in a haploid genome. Using hydrogen peroxide-treated Arabidopsis thaliana plants as a model, we evaluated genomic DNA and RNA-p as reference samples and compared expression levels inferred through the reference relative to unexposed plants with expression levels measured directly using an RNA-0 reference. Our analysis demonstrates that while genomic DNA can serve as a reasonable reference source for microarray assays, a much greater correlation with direct measurements can be achieved using an RNA-based reference sample.

Monitoring of representational difference analysis subtraction procedures by global microarrays

Various approaches to the study of differential gene expression are applied to compare cell lines and tissue samples in a wide range of biological contexts. The compromise between focusing on only the important genes in certain cellular processes and achieving a complete picture is critical for the selection of strategy. We demonstrate how global microarray technology can be used for the exploration of the differentially expressed genes extracted through representational difference analysis (RDA). The subtraction of ubiquitous gene fragments from the two samples was demonstrated using cDNA microarrays including more than 32 000 spotted, PCR-amplified human clones. Hybridizations indicated the expression of 9100 of the microarray elements in a macrophage/foam cell atherosclerosis model system, of which many were removed during the RDA process. The stepwise subtraction procedure was demonstrated to yield an efficient enrichment of gene fragments overrepresented in either sample (18% in the representations, 86% after the first subtraction, and 88% after the second subtraction), many of which were impossible to detect in the starting material. Interestingly, the method allowed for the observation of the differential expression of several members of the low-abundant nuclear receptor gene family. We also observed a certain background level in the difference products of nondifferentially expressed gene fragments, warranting a verification strategy for selected candidate genes. The differential expression of several genes was verified by real-time PCR.

Characterization of a Vibrio vulnificus LysR homologue, HupR, which regulates expression of the haem uptake outer membrane protein, HupA

In Vibrio vulnificus, the ability to acquire iron from the host has been shown to correlate with virulence. Here, we show that the DNA upstream of hupA (haem uptake receptor) in V. vulnificus encodes a protein in the inverse orientation to hupA (named hupR). HupR shares homology with the LysR family of positive transcriptional activators. A hupA-lacZ fusion contained on a plasmid was transformed into Fur(-), Fur(+)and HupR(-)strains of V. vulnificus. The beta-galactosidase assays and Northern blot analysis showed that transcription of hupA is negatively regulated by iron and the Fur repressor in V. vulnificus. Under low-iron conditions with added haemin, the expression of hupA in the hupR mutant was significantly lower than in the wild-type. This diminished response to haem was detected by both Northern blot and hupA-lacZ fusion analysis. The haem response of hupA in the hupR mutant was restored to wild-type levels when complemented with hupR in trans. These studies suggest that HupR may act as a positive regulator of hupA transcription under low-iron conditions in the presence of haemin.

Minimum information about a microarray experiment (MIAME)-toward standards for microarray data

Microarray analysis has become a widely used tool for the generation of gene expression data on a genomic scale. Although many significant results have been derived from microarray studies, one limitation has been the lack of standards for presenting and exchanging such data. Here we present a proposal, the Minimum Information About a Microarray Experiment (MIAME), that describes the minimum information required to ensure that microarray data can be easily interpreted and that results derived from its analysis can be independently verified. The ultimate goal of this work is to establish a standard for recording and reporting microarray-based gene expression data, which will in turn facilitate the establishment of databases and public repositories and enable the development of data analysis tools. With respect to MIAME, we concentrate on defining the content and structure of the necessary information rather than the technical format for capturing it.

Identification of tumor markers in models of human colorectal cancer using a 19,200-element complementary DNA microarray

Metastasis represents a crucial transition in disease development and progression and has a profound impact on survival for a wide variety of cancers. Cell line models of metastasis have played an important role in developing our understanding of the metastatic process. We used a 19,200-element human cDNA microarray to profile transcription in three paired cell-line models of colorectal tumor metastasis. By correlating expression patterns across these cell lines, we have identified 176 genes that appear to be differentially expressed (greater than 2-fold) in all highly metastatic cell lines relative to their reference. An analysis of these genes reiterates much of our understanding of the metastatic process and suggests additional genes, many of previously uncharacterized function, that may be causatively involved in, or at least prognostic of, metastasis. Northern analysis of a limited number of these genes validates the observed pattern of expression and suggests that further investigation and functional characterization of the identified genes is warranted.

The power of public access: the human genome project and the scientific process

The scientific process, and scientific progress, require a critical examination of all published reports. Recent publications detailing errors in the draft human genome sequence are an integral part of our quest to better understand nature and demonstrate the value of free access to scientific data.

Exploring the transcriptome of the malaria sporozoite stage

Most studies of gene expression in Plasmodium have been concerned with asexual and/or sexual erythrocytic stages. Identification and cloning of genes expressed in the preerythrocytic stages lag far behind. We have constructed a high quality cDNA library of the Plasmodium sporozoite stage by using the rodent malaria parasite P. yoelii, an important model for malaria vaccine development. The technical obstacles associated with limited amounts of RNA material were overcome by PCR-amplifying the transcriptome before cloning. Contamination with mosquito RNA was negligible. Generation of 1,972 expressed sequence tags (EST) resulted in a total of 1,547 unique sequences, allowing insight into sporozoite gene expression. The circumsporozoite protein (CS) and the sporozoite surface protein 2 (SSP2) are well represented in the data set. A BLASTX search with all tags of the nonredundant protein database gave only 161 unique significant matches (P(N) < or = 10(-4)), whereas 1,386 of the unique sequences represented novel sporozoite-expressed genes. We identified ESTs for three proteins that may be involved in host cell invasion and documented their expression in sporozoites. These data should facilitate our understanding of the preerythrocytic Plasmodium life cycle stages and the development of preerythrocytic vaccines.

Effects of ischemia on gene expression

Microarray gene expression technology has recently made it feasible to characterize the RNA expression of thousands of genes across numerous tissue samples. We hypothesized that the warm ischemia commonly associated with the surgical extirpation of human tissue would have significant effects on gene expression profiles. To quantitate the effects of warm ischemia on human tissue, we rapidly dissected normal mucosa from a human colon cancer specimen. The specimen was divided and maintained at room temperature until snap-frozen in liquid nitrogen. Aliquots of tissue were frozen at times 5, 10, 15, 20, 40, and 60 min after extirpation. Spotted microarrays composed of 2400 distinct elements were used to assay mRNA derived from each time point in triplicate. Eisen's hierarchical clustering methodology and Bayesean statistical methods were then used to assay the effects of warm ischemia on gene expression. Application of time-course statistical models suggest that three patterns were induced by ischemia, accounting for 68.2, 17.8, and 13.4% of the evaluable genes, respectively. Pattern I corresponds to an average change of 27% over 60 min from 5 min baseline level of expression and 63.8% of the genes with at least 80% probability of membership in this pattern show average increases in expression over 60 min. The remainder decrease on average. Pattern II genes show the least ischemia-related effects, demonstrating an average change of only 12% over 60 min. In contrast to pattern I, we find that 67.5% of the genes with at least 80% probability of membership in this pattern are decreasing in expression on average over time. The remaining 32.5% in this pattern increase an average of 12% over 60 min. Finally, pattern III genes (13.4% of the sample) show the greatest sensitivity to ischemia, changing an average of 50% over 60 min, with about the same number increasing as are decreasing. Fold changes in RNA over- or under-expression were observed up to greater than 20-fold. Warm ischemia associated with the surgical extirpation of human tissues has significant effects on gene expression. These data support the careful monitoring of ischemic time for tissues harvested for the purpose of gene profiling.

Microarray analysis of the in vivo effects of hypophysectomy and growth hormone treatment on gene expression in the rat

Complementary DNA microarrays containing 3000 different rat genes were used to study the consequences of severe hormonal deficiency (hypophysectomy) on the gene expression patterns in heart, liver, and kidney. Hybridization signals were seen from a majority of the arrayed complementary DNAs; nonetheless, tissue-specific expression patterns could be delineated. Hypophysectomy affected the expression of genes involved in a variety of cellular functions. Between 16-29% of the detected transcripts from each tissue changed expression level as a reaction to this condition. Chronic treatment of hypophysectomized animals with human GH also caused significant changes in gene expression patterns. The study confirms previous knowledge concerning certain gene expression changes in the above-mentioned situations and provides new information regarding hypophysectomy and chronic human GH effects in the rat. Furthermore, we have identified several new genes that respond to GH treatment. Our results represent a first step toward a more global understanding of gene expression changes in states of hormonal deficiency.

Computational analysis of microarray data

Microarray experiments are providing unprecedented quantities of genome-wide data on gene-expression patterns. Although this technique has been enthusiastically developed and applied in many biological contexts, the management and analysis of the millions of data points that result from these experiments has received less attention. Sophisticated computational tools are available, but the methods that are used to analyse the data can have a profound influence on the interpretation of the results. A basic understanding of these computational tools is therefore required for optimal experimental design and meaningful data analysis.

Sequence evaluation of four pooled-tissue normalized bovine cDNA libraries and construction of a gene index for cattle

An essential component of functional genomics studies is the sequence of DNA expressed in tissues of interest. To provide a resource of bovine-specific expressed sequence data and facilitate this powerful approach in cattle research, four normalized cDNA libraries were produced and arrayed for high-throughput sequencing. The libraries were made with RNA pooled from multiple tissues to increase efficiency of normalization and maximize the number of independent genes for which sequence data were obtained. Target tissues included those with highest likelihood to have impact on production parameters of animal health, growth, reproductive efficiency, and carcass merit. Success of normalization and inter- and intralibrary redundancy were assessed by collecting 6000-23,000 sequences from each of the libraries (68,520 total sequences deposited in GenBank). Sequence comparison and assembly of these sequences was performed in combination with 56,500 other bovine EST sequences present in the GenBank dbEST database to construct a cattle Gene Index (available from The Institute for Genomic Research at http://www.tigr.org/tdb/tgi.shtml). The 124,381 bovine ESTs present in GenBank at the time of the analysis form 16,740 assemblies that are listed and annotated on the Web site. Analysis of individual library sequence data indicates that the pooled-tissue approach was highly effective in preparing libraries for efficient deep sequencing.

Rice bioinformatics. analysis of rice sequence data and leveraging the data to other plant species

Rice (Oryza sativa) is a model species for monocotyledonous plants, especially for members in the grass family. Several attributes such as small genome size, diploid nature, transformability, and establishment of genetic and molecular resources make it a tractable organism for plant biologists. With an estimated genome size of 430 Mb (Arumuganathan and Earle, 1991), it is feasible to obtain the complete genome sequence of rice using current technologies. An international effort has been established and is in the process of sequencing O. sativa spp. japonica var "Nipponbare" using a bacterial artificial chromosome/P1 artificial chromosome shotgun sequencing strategy. Annotation of the rice genome is performed using prediction-based and homology-based searches to identify genes. Annotation tools such as optimized gene prediction programs are being developed for rice to improve the quality of annotation. Resources are also being developed to leverage the rice genome sequence to partial genome projects such as expressed sequence tag projects, thereby maximizing the output from the rice genome project. To provide a low level of annotation for rice genomic sequences, we have aligned all rice bacterial artificial chromosome/P1 artificial chromosome sequences with The Institute of Genomic Research Gene Indices that are a set of nonredundant transcripts that are generated from nine public plant expressed sequence tag projects (rice, wheat, sorghum, maize, barley, Arabidopsis, tomato, potato, and barrel medic). In addition, we have used data from The Institute of Genomic Research Gene Indices and the Arabidopsis and Rice Genome Projects to identify putative orthologues and paralogues among these nine genomes.

Functional annotation of a full-length mouse cDNA collection

The RIKEN Mouse Gene Encyclopaedia Project, a systematic approach to determining the full coding potential of the mouse genome, involves collection and sequencing of full-length complementary DNAs and physical mapping of the corresponding genes to the mouse genome. We organized an international functional annotation meeting (FANTOM) to annotate the first 21,076 cDNAs to be analysed in this project. Here we describe the first RIKEN clone collection, which is one of the largest described for any organism. Analysis of these cDNAs extends known gene families and identifies new ones.

The TIGR Gene Indices: analysis of gene transcript sequences in highly sampled eukaryotic species

While genome sequencing projects are advancing rapidly, EST sequencing and analysis remains a primary research tool for the identification and categorization of gene sequences in a wide variety of species and an important resource for annotation of genomic sequence. The TIGR Gene Indices (http://www.tigr.org/tdb/tgi. shtml) are a collection of species-specific databases that use a highly refined protocol to analyze EST sequences in an attempt to identify the genes represented by that data and to provide additional information regarding those genes. Gene Indices are constructed by first clustering, then assembling EST and annotated gene sequences from GenBank for the targeted species. This process produces a set of unique, high-fidelity virtual transcripts, or Tentative Consensus (TC) sequences. The TC sequences can be used to provide putative genes with functional annotation, to link the transcripts to mapping and genomic sequence data, to provide links between orthologous and paralogous genes and as a resource for comparative sequence analysis.

Anchoring of rice BAC clones to the rice genetic map in silico

A wealth of molecular resources have been developed for rice genomics, including dense genetic maps, expressed sequence tags (ESTs), yeast artificial chromosome maps, bacterial artificial chromosome (BAC) libraries and BAC end sequence databases. Integration of genetic and physical maps involves labor-intensive empirical experiments. To accelerate the integration of the bacterial clone resources with the genetic map for the International Rice Genome Sequencing Project, we cleaned and filtered the available EST and BAC end sequences for repetitive sequences and then searched all available rice genetic markers with our filtered databases. We identified 418 genetic markers that aligned with at least one BAC end sequence with >95% sequence identity, providing a set of large insert clones with an average separation of 1 Mb that can serve as nucleation points for the sequencing phase of the International Rice Genome Sequencing Project.

An optimized protocol for analysis of EST sequences

The vast body of Expressed Sequence Tag (EST) data in the public databases provide an important resource for comparative and functional genomics studies and an invaluable tool for the annotation of genomic sequences. We have developed a rigorous protocol for reconstructing the sequences of transcribed genes from EST and gene sequence fragments. A key element in developing this protocol has been the evaluation of a number of sequence assembly programs to determine which most faithfully reproduce transcript sequences from EST data. The TIGR Gene Indices constructed using this protocol for human, mouse, rat and a variety of other plant and animal models have demonstrated their utility in a variety of applications and are freely available to the scientific research community.

A concise guide to cDNA microarray analysis

Microarray expression analysis has become one of the most widely used functional genomics tools. Efficient application of this technique requires the development of robust and reproducible protocols. We have optimized all aspects of the process, including PCR amplification of target cDNA clones, microarray printing, probe labeling and hybridization, and have developed strategies for data normalization and analysis.

Gene index analysis of the human genome estimates approximately 120,000 genes

Although sequencing of the human genome will soon be completed, gene identification and annotation remains a challenge. Early estimates suggested that there might be 60,000-100,000 (ref. 1) human genes, but recent analyses of the available data from EST sequencing projects have estimated as few as 45,000 (ref. 2) or as many as 140, 000 (ref. 3) distinct genes. The Chromosome 22 Sequencing Consortium estimated a minimum of 45,000 genes based on their annotation of the complete chromosome, although their data suggests there may be additional genes. The nearly 2,000,000 human ESTs in dbEST provide an important resource for gene identification and genome annotation, but these single-pass sequences must be carefully analysed to remove contaminating sequences, including those from genomic DNA, spurious transcription, and vector and bacterial sequences. We have developed a highly refined and rigorously tested protocol for cleaning, clustering and assembling EST sequences to produce high-fidelity consensus sequences for the represented genes (F.L. et al., manuscript submitted) and used this to create the TIGR Gene Indices-databases of expressed genes for human, mouse, rat and other species (http://www.tigr.org/tdb/tgi.html). Using highly refined and tested algorithms for EST analysis, we have arrived at two independent estimates indicating the human genome contains approximately 120,000 genes.

The African trypanosome genome

The haploid nuclear genome of the African trypanosome, Trypanosoma brucei, is about 35 Mb and varies in size among different trypanosome isolates by as much as 25%. The nuclear DNA of this diploid organism is distributed among three size classes of chromosomes: the megabase chromosomes of which there are at least 11 pairs ranging from 1 Mb to more than 6 Mb (numbered I-XI from smallest to largest); several intermediate chromosomes of 200-900 kb and uncertain ploidy; and about 100 linear minichromosomes of 50-150 kb. Size differences of as much as four-fold can occur, both between the two homologues of a megabase chromosome pair in a specific trypanosome isolate and among chromosome pairs in different isolates. The genomic DNA sequences determined to date indicated that about 50% of the genome is coding sequence. The chromosomal telomeres possess TTAGGG repeats and many, if not all, of the telomeres of the megabase and intermediate chromosomes are linked to expression sites for genes encoding variant surface glycoproteins (VSGs). The minichromosomes serve as repositories for VSG genes since some but not all of their telomeres are linked to unexpressed VSG genes. A gene discovery program, based on sequencing the ends of cloned genomic DNA fragments, has generated more than 20 Mb of discontinuous single-pass genomic sequence data during the past year, and the complete sequences of chromosomes I and II (about 1 Mb each) in T. brucei GUTat 10.1 are currently being determined. It is anticipated that the entire genomic sequence of this organism will be known in a few years. Analysis of a test microarray of 400 cDNAs and small random genomic DNA fragments probed with RNAs from two developmental stages of T. brucei demonstrates that the microarray technology can be used to identify batteries of genes differentially expressed during the various life cycle stages of this parasite.

The TIGR gene indices: reconstruction and representation of expressed gene sequences

Expressed sequence tags (ESTs) have provided a first glimpse of the collection of transcribed sequences in a variety of organisms. However, a careful analysis of this sequence data can provide significant additional functional, structural and evolutionary information. Our analysis of the public EST sequences, available through the TIGR Gene Indices (TGI; http://www.tigr.org/tdb/tdb.html ), is an attempt to identify the genes represented by that data and to provide additional information regarding those genes. Gene Indices are constructed for selected organisms by first clustering, then assembling EST and annotated gene sequences from GenBank. This process produces a set of unique, high-fidelity virtual transcripts, or tentative consensus (TC) sequences. The TC sequences can be used to provide putative genes with functional annotation, to link the transcripts to mapping and genomic sequence data, and to provide links between orthologous and paralogous genes.

Genomic sequence of a Lyme disease spirochaete, Borrelia burgdorferi

The genome of the bacterium Borrelia burgdorferi B31, the aetiologic agent of Lyme disease, contains a linear chromosome of 910,725 base pairs and at least 17 linear and circular plasmids with a combined size of more than 533,000 base pairs. The chromosome contains 853 genes encoding a basic set of proteins for DNA replication, transcription, translation, solute transport and energy metabolism, but, like Mycoplasma genitalium, it contains no genes for cellular biosynthetic reactions. Because B. burgdorferi and M. genitalium are distantly related eubacteria, we suggest that their limited metabolic capacities reflect convergent evolution by gene loss from more metabolically competent progenitors. Of 430 genes on 11 plasmids, most have no known biological function; 39% of plasmid genes are paralogues that form 47 gene families. The biological significance of the multiple plasmid-encoded genes is not clear, although they may be involved in antigenic variation or immune evasion.

The complete genome sequence of the hyperthermophilic, sulphate-reducing archaeon Archaeoglobus fulgidus

Archaeoglobus fulgidus is the first sulphur-metabolizing organism to have its genome sequence determined. Its genome of 2,178,400 base pairs contains 2,436 open reading frames (ORFs). The information processing systems and the biosynthetic pathways for essential components (nucleotides, amino acids and cofactors) have extensive correlation with their counterparts in the archaeon Methanococcus jannaschii. The genomes of these two Archaea indicate dramatic differences in the way these organisms sense their environment, perform regulatory and transport functions, and gain energy. In contrast to M. jannaschii, A. fulgidus has fewer restriction-modification systems, and none of its genes appears to contain inteins. A quarter (651 ORFs) of the A. fulgidus genome encodes functionally uncharacterized yet conserved proteins, two-thirds of which are shared with M. jannaschii (428 ORFs). Another quarter of the genome encodes new proteins indicating substantial archaeal gene diversity.

The complete genome sequence of the gastric pathogen Helicobacter pylori

Helicobacter pylori, strain 26695, has a circular genome of 1,667,867 base pairs and 1,590 predicted coding sequences. Sequence analysis indicates that H. pylori has well-developed systems for motility, for scavenging iron, and for DNA restriction and modification. Many putative adhesins, lipoproteins and other outer membrane proteins were identified, underscoring the potential complexity of host-pathogen interaction. Based on the large number of sequence-related genes encoding outer membrane proteins and the presence of homopolymeric tracts and dinucleotide repeats in coding sequences, H. pylori, like several other mucosal pathogens, probably uses recombination and slipped-strand mispairing within repeats as mechanisms for antigenic variation and adaptive evolution. Consistent with its restricted niche, H. pylori has a few regulatory networks, and a limited metabolic repertoire and biosynthetic capacity. Its survival in acid conditions depends, in part, on its ability to establish a positive inside-membrane potential in low pH.

An STS-based radiation hybrid map of the human genome

We have constructed a physical map of the human genome by using a panel of 83 whole genome radiation hybrids (the Stanford G3 panel) in conjunction with 10,478 sequence-tagged sites (STSs) derived from random genomic DNA sequences, previously mapped genetic markers, and expressed sequences. Of these STSs, 5049 are framework markers that fall into 1766 high-confidence bins. An additional 945 STSs are indistinguishable in their map location from one or more of the framework markers. These 5994 mapped STSs have an average spacing of 500 kb. An additional 4484 STSs are positioned with respect to the framework markers. Comparison of the orders of markers on this map with orders derived from independent meiotic and YAC STS-content maps indicates that the error rate in defining high-confidence bins is < 5%. Analysis of 322 random cDNAs indicates that the map covers the vast majority of the human genome. This STS-based radiation hybrid map of the human genome brings us one step closer to the goal of a physical map containing 30,000 unique ordered landmarks with an average marker spacing of 100 kb.

A nucleotide substitution in the promoter of human angiotensinogen is associated with essential hypertension and affects basal transcription in vitro

In earlier studies, we provided statistical evidence that individual differences in the angiotensinogen gene, the precursor of the vasoactive hormone angiotensin II, constitute inherited predispositions to essential hypertension in humans. We have now identified a common variant in the proximal promoter, the presence of an adenine, instead of a guanine, 6 bp upstream from the initiation site of transcription, in significant association with the disorder. Tests of promoter activity and DNA binding studies with nuclear proteins suggest that this nucleotide substitution affects the basal transcription rate of the gene. These observations provide some biological insight about the possible mechanism of a genetic predisposition to essential hypertension; they may also have important evolutionary implications.