A comparison of the Celera and Ensembl predicted gene sets reveals little overlap in novel genes

TAGADA: a scalable pipeline to improve genome annotations with RNA-seq data

Genome annotation plays a crucial role in providing comprehensive catalog of genes and transcripts for a particular species. As research projects generate new transcriptome data worldwide, integrating this information into existing annotations becomes essential. However, most bioinformatics pipelines are limited in their ability to effectively and consistently update annotations using new RNA-seq data. Here we introduce TAGADA, an RNA-seq pipeline for Transcripts And Genes Assembly, Deconvolution, and Analysis. Given a genomic sequence, a reference annotation and RNA-seq reads, TAGADA enhances existing gene models by generating an improved annotation. It also computes expression values for both the reference and novel annotation, identifies long non-coding transcripts (lncRNAs), and provides a comprehensive quality control report. Developed using Nextflow DSL2, TAGADA offers user-friendly functionalities and ensures reproducibility across different computing platforms through its containerized environment. In this study, we demonstrate the efficacy of TAGADA using RNA-seq data from the GENE-SWiTCH project alongside chicken and pig genome annotations as references. Results indicate that TAGADA can substantially increase the number of annotated transcripts by approximately [Formula: see text] in these species. Furthermore, we illustrate how TAGADA can integrate Illumina NovaSeq short reads with PacBio Iso-Seq long reads, showcasing its versatility. TAGADA is available at github.com/FAANG/analysis-TAGADA.

Rat mir-155 generated from the lncRNA Bic is 'hidden' in the alternate genomic assembly and reveals the existence of novel mammalian miRNAs and clusters

MicroRNAs (miRNAs) are a class of small noncoding RNAs acting as post-transcriptional gene expression regulators in many physiological and pathological conditions. During the last few years, many novel mammalian miRNAs have been predicted experimentally with bioinformatics approaches and validated by next-generation sequencing. Although these strategies have prompted the discovery of several miRNAs, the total number of these genes still seems larger. Here, by exploiting the species conservation of human, mouse, and rat hairpin miRNAs, we discovered a novel rat microRNA, mir-155. We found that mature miR-155 is overexpressed in rat spleen myeloid cells treated with LPS, similarly to humans and mice. Rat mir-155 is annotated only on the alternate genome, suggesting the presence of other "hidden" miRNAs on this assembly. Therefore, we comprehensively extended the homology search also to mice and humans, finally validating 34 novel mammalian miRNAs (two in humans, five in mice, and up to 27 in rats). Surprisingly, 15 of these novel miRNAs (one for mice and 14 for rats) were found only on the alternate and not on the reference genomic assembly. To date, our findings indicate that the choice of genomic assembly, when mapping small RNA reads, is an important option that should be carefully considered, at least for these animal models. Finally, the discovery of these novel mammalian miRNA genes may contribute to a better understanding of already acquired experimental data, thereby paving the way to still unexplored investigations and to unraveling the function of miRNAs in disease models.

ANALYSIS OF THE EVOLVING PROTEOMES: PREDICTIONS OF THE NUMBER OF PROTEIN DOMAINS IN NATURE AND THE NUMBER OF GENES IN EUKARYOTIC ORGANISMS

Motivation: Obtaining accurate estimates of the numbers of protein-coding genes and protein domains in a proteome, and the number of protein domains in nature is a daunting challenge. Computational analysis of the protein domain sets in the proteomes of many species allows us to estimate these numbers and to find their evolution relationships.

Results: We have analyzed the distributions of the number of occurrences of protein domains in sample proteomes of the 70 fully sequenced genome organisms of three major kingdoms of life: Archaea, Bacteria and Eukaryota. We found that a large fraction of the identified distinct protein domains (i.e., unique domains and homologous domain families) in these 70 proteomes (1051 (23%) out of 4493) are found in at least one organism in each of these kingdoms of life and that 43 (1%) of these domains are common to all the 70 organisms. All the observed domain occurrence frequency distributions for these 70 proteomes are well fitted by a family of Pareto-like functions, associated with the steady state distributions of a linear Markov random process. We present explicit formulas that accurately predict the number of distinct protein domains and the number of protein-coding genes for a given organism as functions of the number of non-redundant domain-to-protein links in the proteomes. These functions allows us to predict that there are 42,740, 27,900, and 21,200 protein-coding genes/open reading frames in the human, A. thaliana, and mouse genomes, respectively. We also estimate that there are 5271, 2955, and 4915 distinct protein domains in the human, A. thaliana, and mouse proteomes, respectively, and about 5500 distinct protein domains in the entire "proteome world".

Identifying HIV-1 host cell factors by genome-scale RNAi screening

Advances in the application of RNA interference (RNAi) have facilitated the establishment of systematic cell-based loss-of-function screening platforms. Widespread implementation of this technology has enabled genome-wide genetic analysis of a diverse array of cellular phenotypes, including the identification of host cell factors involved in viral replication. Four recent studies employed whole-genome RNAi technologies to elucidate cellular genes important for the replication of HIV-1. While these four genome-scale screens shared a common objective, they differ in their scope and experimental design. In this review we explore alternative strategies for developing RNAi screens, and discuss potential pitfalls of the technology. Important technical considerations include the choice of silencing reagents, experimental systems, assay readout and analysis methods. We focus on experimental and computational parameters that can impact the outcome of high-throughput genetic screens, and provide guidelines for the development of reliable cell-based RNAi screens.

Using the DFCI gene index databases for biological discovery

The DFCI Gene Index Web pages provide access to analyses of ESTs and gene sequences for nearly 114 species, as well as a number of resources derived from these. Each species-specific database is presented using a common format with a home page. A variety of methods exist that allow users to search each species-specific database. Methods implemented currently include nucleotide or protein sequence queries using WU-BLAST, text-based searches using various sequence identifiers, searches by gene, tissue and library name, and searches using functional classes through Gene Ontology assignments. This protocol provides guidance for using the Gene Index Databases to extract information.

High throughput techniques for characterizing the expression profile of Barrett's esophagus

Barrett's esophagus (BE) is the metaplastic change of the normal lined squamous epithelium of the distal esophagus to a columnar type of epithelium as a result of chronic long-standing gastroesophageal reflux disease. Patients with BE have a significantly increased risk of developing an esophageal adenocarcinoma, with an estimated annual incidence varying from 0.4 to 1.8%. Over the last 3 decades, the incidence of BE and its associated adenocarcinoma has increased in Western countries at a rate that exceeds that of any other malignancy. Despite all the research performed on BE, there is still an inadequate understanding of the biological basis of this mucosal transformation. With the upcoming modern high throughput technologies, important progression has been made in unraveling the expression profiles and gaining more insight in the biology of BE and esophageal adenocarcinoma. Several studies reported genome, transcriptome, proteome, and kinome investigations using high throughput techniques. These studies were conducted to find biomarkers that can be used to detect BE patients with increased risk for malignant progression or to obtain more insight in the mechanism underlying BE development. In the following review, we first discuss the different techniques that are currently available and summarize findings in this field, including several recent publications of our group.

Using the TIGR gene index databases for biological discovery

The TIGR Gene Index web pages provide access to analyses of ESTs and gene sequences for nearly 60 species, as well as a number of resources derived from these. Each species-specific database is presented using a common format with a homepage. A variety of methods exist that allow users to search each species-specific database. Methods implemented currently include nucleotide or protein sequence queries using WU-BLAST, text-based searches using various sequence identifiers, searches by gene, tissue and library name, and searches using functional classes through Gene Ontology assignments. This protocol provides guidance for using the Gene Index Databases to extract information.

Genome assembly forensics: finding the elusive mis-assembly

A collection of software tools is combined for the first time in an automated pipeline for detecting large-scale genome assembly errors and for validating genome assemblies.

We present the first collection of tools aimed at automated genome assembly validation. This work formalizes several mechanisms for detecting mis-assemblies, and describes their implementation in our automated validation pipeline, called amosvalidate. We demonstrate the application of our pipeline in both bacterial and eukaryotic genome assemblies, and highlight several assembly errors in both draft and finished genomes. The software described is compatible with common assembly formats and is released, open-source, at .

Cloning of novel hepatoma gene in rats by invigorating the spleen to supplement Qi

AIM: To clone the integrated cDNA expression sequence of EST segments in down-regulated genes by invigorating the spleen to supplement Qi.

METHODS: Integrated cDNA expression sequences of EST gene segments were cloned for the 689 down-regulated genes by invigorating the spleen to supplement Qi, with electron cloning in combination with PCR.

RESULTS: cDNA expression sequences were detected in 11 EST segments (G2, G4, G5, G6, G11, G14, G15, G16, G17, G18, G20). BLAST analysis showed that G14, G15 and G20 were novel genes which were submitted to GenBank (their accession number is DQ480745, DQ480746 and DQ480747, respectively).

CONCLUSION: Invigorating the spleen to supplement Qi can clone the cDNA expression sequences of EST segments in the down-regulated genes. Further study is needed to observe the functions of these novel genes and the mechanism of action of different TCM therapies.

Driven polymer transport through a nanopore controlled by a rotating electric field: off-lattice computer simulations

The driven translocation kinetics of a single strand polynucleotide chain through a nanopore is studied using off-lattice Monte Carlo simulations, by which the authors demonstrate a novel method in controlling the driven polymer transport through a nanopore by a rotating electric field. The recorded time series of blockade current from the driven polynucleotide transport are used to determine the sequence of polynucleotides by implementing a modified Monte Carlo algorithm, in which the energy landscape paving technique is incorporated to avoid trapping at deep local minima. It is found that only six-time series of block current are required to completely determine the polynucleotide sequence if the average missing rate (AMR) of current signals in these time series is smaller than 20%. For those time series with AMR greater than 20%, the error rate in sequencing an unknown polynucleotide decreases rapidly by increasing the number of time series. To find the most appropriate experimental conditions, the authors have investigated the dependence of AMR of current signals and qualified rate of measured time series of blockade current on various controllable experimental variables.

Detecting transcriptionally active regions using genomic tiling arrays

We have developed a method for interpreting genomic tiling array data, implemented as the program TranscriptionDetector. Probed loci expressed above background are identified by combining replicates in a way that makes minimal assumptions about the data. We performed medium-resolution Anopheles gambiae tiling array experiments and found extensive transcription of both coding and non-coding regions. Our method also showed improved detection of transcriptional units when applied to high-density tiling array data for ten human chromosomes.

Differential regulation of gene products in newly synthesized Brassica napus allotetraploids is not related to protein function nor subcellular localization

Background

Allopolyploidy is a preeminent process in plant evolution that results from the merger of distinct genomes in a common nucleus via inter-specific hybridization. Allopolyploid formation is usually related to genome-wide structural and functional changes though the underlying mechanisms operating during this "genomic shock" still remain poorly known. The aim of the present study was to investigate the modifications occurring at the proteomic level following an allopolyploidization event and to determine whether these changes are related to functional properties of the proteins. In a previous report, we applied comparative proteomics to synthetic amphiploids of Brassica napus and to its diploid progenitors B. rapa and B. oleracea. Although several hundred polypeptides displayed additivity (i.e. mid-parent values) in the amphiploids, many of them showed non-additivity. Here, we report the in silico functional characterization of the "non-additive" proteins (the ones with a non-additive pattern of regulation) in synthetic B. napus.

Results

The complete set of non-additive proteins (335 in the stem and 205 in the root), as well as a subset of additive polypeptides (200 per organ), was identified by mass spectrometry. Several protein isoforms were found, and most of them (~55%) displayed "different" or "opposite" patterns of regulation in the amphiploids, i.e. isoforms of the same protein showing both up-regulation and down-regulation in the synthetic B. napus compared to the mid-parent value. Components of protein complexes were identified of which ~50% also displayed "different" or "opposite" patterns of regulation in the allotetraploids. In silico functional categorization of the identified proteins was carried out, and showed that neither functional category nor metabolic pathway were systematically affected by non-additivity in the synthetic amphiploids. In addition, no subcellular compartment was found to be over- or under-represented among the proteins displaying non-additive values in the allopolyploids.

Conclusion

Protein identification showed that functionally related polypeptides (isoforms and complex subunits) could be differentially regulated in synthetic B. napus in comparison to its diploid progenitors while such proteins are usually expected to display co-regulation. The genetic redundancy within an allopolyploid could explain why functionally related proteins could display imbalanced levels of expression. No functional category, no metabolic pathway and no subcellular localization was found to be over- or under-represented within non-additive polypeptides, suggesting that the differential regulation of gene products was not related to functional properties of the proteins. Thus, at the protein level, there is no evidence for the "genomic shock" expected in neo-polyploids and the overall topology of protein networks and metabolic pathways is conserved in synthetic allotetraploids of B. napus in comparison to its diploid progenitors B. rapa and B. oleracea.

Insulin-like growth factor binding protein 2: gene expression microarrays and the hypothesis-generation paradigm

A major goal of modern medicine is to identify key genes and their products that are altered in the diseased state and to elucidate the molecular mechanisms underlying disease development, progression, and resistance to therapy. This is a daunting task given the exceptionally high complexity of the human genome. The paradigm for research has historically been hypothesis-driven despite the fact that the hypotheses under scrutiny often rest on tenuous subjective grounds or are derived from and dependent on chance observation. The imminent deciphering of the complete human genome, coupled with recent advances in high-throughput bioanalytical technology, has made possible a new paradigm in which data-based hypothesis-generation is the initial step in the investigative process, followed by hypothesis-testing. Genomics technologies are the primary source of the new hypothesis-generating capabilities that are now empowering biomedical researchers. The synergistic interaction between contemporary genomics technologies and the hypothesis-generation paradigm is well-illustrated by the discovery and subsequent ongoing study of the role of insulin-like growth factor binding protein 2 (IGFBP2) in human glioma biology. Using gene expression microarray technology, the IGFBP2 gene was recently found to be highly and differentially overexpressed in the most advanced grade of human glioma, glioblastoma. Based on this discovery, subsequent functional studies were initiated that suggest that IGFBP2 overexpression may contribute to the invasive nature of glioblastoma, and that IGFBP2 may exert its function via a newly identified novel binding protein. The IGFBP2 story is but one example of the power and potential of the new molecular methodologies that are transforming modern diagnostic and investigative neuropathology.

A high throughput messenger RNA differential display screen identifies discrete domains of gene expression and novel patterning processes along the developing neural tube

BACKGROUND

During early development the vertebrate neural tube is broadly organized into the forebrain, midbrain, hindbrain and spinal cord regions. Each of these embryonic zones is patterned by a combination of genetic pathways and the influences of local signaling centres. However, it is clear that much remains to be learned about the complete set of molecular cues that are employed to establish the identity and intrinsic neuronal diversity of these territories. In order to address this, we performed a high-resolution messenger RNA differential display screen to identify molecules whose expression is regionally restricted along the anteroposterior (AP) neuraxis during early chick development, with particular focus on the midbrain and hindbrain vesicles.

RESULTS

This approach identified 44 different genes, with both known and unknown functions, whose transcription is differentially regulated along the AP axis. The identity and ontological classification of these genes is presented. The wide variety of functional classes of transcripts isolated in this screen reflects the diverse spectrum of known influences operating across these embryonic regions. Of these 44 genes, several have been selected for detailed in situ hybridization analysis to validate the screen and accurately define the expression domains. Many of the identified cDNAs showed no identity to the current databases of known or predicted genes or ESTs. Others represent genes whose embryonic expression has not been previously reported. Expression studies confirmed the predictions of the primary differential display data. Moreover, the nature of identified genes, not previously associated with regionalisation of the brain, identifies novel potential mechanisms in that process.

CONCLUSION

This study provides an insight into some of the varied and novel molecular networks that operate during the regionalization of embryonic neural tissue and expands our knowledge of molecular repertoire used during development.

The future of dermatopathology

Of the major issues that dermatopathology will face in the immediate future, two powerful challenges loom large. The first is the application of novel nondestructive imaging technologies to in vivo diagnosis in humans. The second is the application of molecular technologies to a diagnostic arena which formerly belonged exclusively to the light microscopist. The first to be considered in this context is the application of near infrared spectroscopy to the noninvasive in vivo diagnosis of neoplastic skin disease. The second will be a discussion of application, methodology and the current state of the art in microarray technologies as they apply to neoplastic dermatopathology and, in particular, the diagnosis and prognostication of melanoma.

Genome sequence, comparative analysis and haplotype structure of the domestic dog

Here we report a high-quality draft genome sequence of the domestic dog (Canis familiaris), together with a dense map of single nucleotide polymorphisms (SNPs) across breeds. The dog is of particular interest because it provides important evolutionary information and because existing breeds show great phenotypic diversity for morphological, physiological and behavioural traits. We use sequence comparison with the primate and rodent lineages to shed light on the structure and evolution of genomes and genes. Notably, the majority of the most highly conserved non-coding sequences in mammalian genomes are clustered near a small subset of genes with important roles in development. Analysis of SNPs reveals long-range haplotypes across the entire dog genome, and defines the nature of genetic diversity within and across breeds. The current SNP map now makes it possible for genome-wide association studies to identify genes responsible for diseases and traits, with important consequences for human and companion animal health.

Echoes from the past--are we still in an RNP world?

Availability of the human genome sequence and those of other species is unmeasured in their value for a comprehensive understanding of the architecture, function and evolution of genomes and cells. Various mechanisms keep genomes in flux and generate intra- and interspecies variation. The conversion of RNA modules into DNA and their more or less random integration into chromosomes (retroposition) is in many lineages including our own the most pervasive and perhaps the most enigmatic. The proclivity of such events in extant multicellular eukaryotes, even in more recent evolutionary times, gives the impression that the transition period from the RNP (ribonucleoprotein) world to the emergence of modern cells, where DNA became the predominant carrier of genetic information, has lasted billions of years and is an endlessly drawn-out process rather than the punctuated event one might expect. Apart from the impact of such RNA-mediated processes as retroposition, the role of RNA in a wide variety of cellular functions has only recently become more widely appreciated.

DNA microarrays: from structural genomics to functional genomics. The applications of gene chips in dermatology and dermatopathology

The human genome project was successful in sequencing the entire human genome and ended earlier than expected. The vast genetic information now available will have far-reaching consequences for medicine in the twenty-first century. The knowledge gained from the mapping and sequencing of human genes on a genome-wide scale--commonly referred to as structural genomics--is prerequisite for studies that focus on the functional aspects of genes. A recently invented technique, known as gene chip, or DNA microarray, technology, allows the study of the function of thousands of genes at once, thereby opening the door to the new field of functional genomics. At its core, the DNA microarray utilizes a unique feature of DNA known as complementary hybridization. As such, it is not different from Southern (DNA) blot or northern (RNA) blot hybridizations, or the polymerase chain reaction, with the exception that it allows expression profiling of the entire human genome in a single hybridization experiment. The article highlights the principles, technology, and applications of DNA microarrays as they pertain to the field of dermatology and dermatopathology. The most important applications are the gene expression profiling of skin cancer, especially of melanoma. Other potential applications include gene expression profiling of inflammatory skin diseases, the mutational analysis of genodermatoses, and polymorphism screening, as well as drug development and chemosensitivity prediction. cDNA microarrays will shape the diagnostic approach of the dermatology and the dermatopathology of the future and may lead to new therapeutic options.

Gene discovery and expression profiling in porcine Peyer's patch

Peyer's patches of the intestinal mucosa are essential for host defense and immune regulation in the enteric system. To better understand molecular mechanisms of Peyer's patch function, we have screened for differentially expressed genes specific to Peyer's patch. cDNA libraries were created from normal Peyer's patch, immune stimulated Peyer's patch, and pooled cDNA subtracted with fibroblast RNA. From the subtracted library, 3687 expressed sequence tags (ESTs), representing 2414 unique nucleotide sequences, were isolated, identified by BLAST searches against public databases, and spotted onto a microarray for gene expression profiling. Approximately 30% of these ESTs BLAST to genes of unknown function and 20% have no known homology in the public databases (novel genes). Of the novel genes, 70% are expressed in normal immune tissues by microarray analysis, suggesting that at least 371 of the unidentified EST sequences from the subtracted library are novel porcine genes and can now be further characterized to determine their function in the porcine Peyer's patch. We surmise that the products of these genes participate in biochemical and cellular functions related to the unique immunological and gastroenterological functions of the small intestine. The BLAST and gene ontology information for each of the subtracted library EST sequences, the normal and immune stimulated libraries, and the microarray are all valuable resources that will facilitate further examination of the biological function of porcine Peyer's patch tissue.

Comparison of the current RefSeq, Ensembl and EST databases for counting genes and gene discovery

Large amounts of refined sequence material in the form of predicted, curated and annotated genes and expressed sequences tags (ESTs) have recently been added to the NCBI databases. We matched the transcript-sequences of RefSeq, Ensembl and dbEST in an attempt to provide an updated overview of how many unique human genes can be found. The results indicate that there are about 25000 unique genes in the union of RefSeq and Ensembl with 12-18% and 8-13% of the genes in each set unique to the other set, respectively. About 20% of all genes had splice variants. There are a considerable number of ESTs (2200000) that do not match the identified genes and we used an in-house pipeline to identify 22 novel genes from Genscan predictions that have considerable EST coverage. The study provides an insight into the current status of human gene catalogues and shows that considerable refinement of methods and datasets is needed to come to a conclusive gene count.

Large-scale RT-PCR recovery of full-length cDNA clones

Pseudogenes, alternative transcripts, noncoding RNA, and polymorphisms each add extensive complexity to the mammalian transcriptome and confound estimation of the total number of genes. Despite advanced algorithms for gene prediction and several large-scale efforts to obtain cDNA clones for all human open reading frames (ORFs), no single collection is complete. To enhance this effort, we have developed a high-throughput pipeline for reverse transcription PCR (RT-PCR) gene recovery. Most importantly, novel molecular strategies for improving RT-PCR yield of transcripts that have been difficult to isolate by other means and computational strategies for clone sequence validation have been developed and optimized. This systematic gene recovery pipeline allows both rescue of predicted human and rat genes and provides insight into the complexity of the transcriptome through comparisons with existing data sets.

A comprehensive transcript index of the human genome generated using microarrays and computational approaches

BACKGROUND

Computational and microarray-based experimental approaches were used to generate a comprehensive transcript index for the human genome. Oligonucleotide probes designed from approximately 50,000 known and predicted transcript sequences from the human genome were used to survey transcription from a diverse set of 60 tissues and cell lines using ink-jet microarrays. Further, expression activity over at least six conditions was more generally assessed using genomic tiling arrays consisting of probes tiled through a repeat-masked version of the genomic sequence making up chromosomes 20 and 22.

RESULTS

The combination of microarray data with extensive genome annotations resulted in a set of 28,456 experimentally supported transcripts. This set of high-confidence transcripts represents the first experimentally driven annotation of the human genome. In addition, the results from genomic tiling suggest that a large amount of transcription exists outside of annotated regions of the genome and serves as an example of how this activity could be measured on a genome-wide scale.

CONCLUSIONS

These data represent one of the most comprehensive assessments of transcriptional activity in the human genome and provide an atlas of human gene expression over a unique set of gene predictions. Before the annotation of the human genome is considered complete, however, the previously unannotated transcriptional activity throughout the genome must be fully characterized.

ORFeome projects: gateway between genomics and omics

The availability of entire genome sequences is expected to revolutionize the way in which biology and medicine are conducted for years to come. However, achieving this promise still requires significant effort in the areas of gene annotation, cloning and expression of thousands of known and heretofore unknown protein-encoding genes. Traditional technologies of manipulating genes are too cumbersome and inefficient when one is dealing with more than a few genes at a time. Entire libraries composed of all protein-encoding open reading frames (ORFs) cloned in highly flexible vectors will be needed to take full advantage of the information found in any genome sequence. The creation of such ORFeome resources using novel technologies for cloning and expressing entire proteomes constitutes an effective gateway from whole genome sequencing efforts to downstream 'omics' applications.

Identification of alternatively spliced mRNA variants related to cancers by genome-wide ESTs alignment

Several databases have been published to predict alternative splicing of mRNAs by analysing the exon linkage relationship by alignment of expressed sequence tags (ESTs) to the genome sequence; however, little effort has been made to investigate the relationship between cancers and alternative splicing. We developed a program, Alternative Splicing Assembler (ASA), to look for splicing variants of human gene transcripts by genome-wide ESTs alignment. Using ASA, we constructed the biosino alternative splicing database (BASD), which predicted splicing variants for reference sequences from the reference sequence database (RefSeq) and presented them in both graph and text formats. EST clusters that differ from the reference sequences in at least one splicing site were counted as splicing variants. Of 4322 genes screened, 3490 (81%) were observed with at least one alternative splicing variants. To discover the variants associated with cancers, tissue sources of EST sequences were extracted from the UniLib database and ESTs from the same tissue type were counted. These were regarded as the indicators for gene expression level. Using Fisher's exact test, alternative splicing variants, of which EST counts were significantly different between cancer tissues and their counterpart normal tissues, were identified. It was predicted that 2149 variants, or 383 variants after Bonferroni correction, of 26 812 variants were likely tumor-associated. By reverse transcription-PCR, 11 of 13 novel alternative splicing variants and eight of nine variants' tissue specificity were confirmed in hepatocellular carcinoma and in lung cancer. The possible involvement of alternative splicing in cancer is discussed.

Impact of microarray technology in clinical oncology

Genomic technologies are rapidly evolving and have demonstrated utility in augmenting oncological pathology or clinical presentation in disease classification and risk of relapse assessment. Numerous malignancies have been subject to microarray examination, and through a variety of analysis methodologies, groups of reporter genes have been identified to generate 'molecular portraits' of these diseases. Once validated, it is likely that assessment of the expression levels of subsets of reporter genes will contribute to personalized genomic medicine through diagnosis and selection of treatment options for patients. The dynamic nature of this field ensures that new developments are missing from this review.

A two-dimensional proteome map of maize endosperm

We have established a proteome reference map for maize (Zea mays L.) endosperm by means of two-dimensional gel electrophoresis and protein identification with LC-MS/MS analysis. This investigation focussed on proteins in major spots in a 4-7 pI range and 10-100 kDa M(r) range. Among the 632 protein spots processed, 496 were identified by matching against the NCBInr and ZMtuc-tus databases (using the SEQUEST software). Forty-two per cent of the proteins were identified against maize sequences, 23% against rice sequences and 21% against Arabidopsis sequences. Identified proteins were not only cytoplasmic but also nuclear, mitochondrial or amyloplastic. Metabolic processes, protein destination, protein synthesis, cell rescue, defense, cell death and ageing are the most abundant functional categories, comprising almost half of the 632 proteins analyzed in our study. This proteome map constitutes a powerful tool for physiological studies and is the first step for investigating the maize endosperm development.

A gene atlas of the mouse and human protein-encoding transcriptomes

The tissue-specific pattern of mRNA expression can indicate important clues about gene function. High-density oligonucleotide arrays offer the opportunity to examine patterns of gene expression on a genome scale. Toward this end, we have designed custom arrays that interrogate the expression of the vast majority of protein-encoding human and mouse genes and have used them to profile a panel of 79 human and 61 mouse tissues. The resulting data set provides the expression patterns for thousands of predicted genes, as well as known and poorly characterized genes, from mice and humans. We have explored this data set for global trends in gene expression, evaluated commonly used lines of evidence in gene prediction methodologies, and investigated patterns indicative of chromosomal organization of transcription. We describe hundreds of regions of correlated transcription and show that some are subject to both tissue and parental allele-specific expression, suggesting a link between spatial expression and imprinting.

Systematic genome-wide approach to positional candidate cloning for identification of novel human disease genes

BACKGROUND

Recent large-scale genome projects afford a unique opportunity to identify many novel disease genes and thereby better understand the genetic basis of human disease. Functional Annotation of Mouse (FANTOM) 2, the largest mouse transcriptome project yet, provides a wealth of data on novel genes, splice variants and non-coding RNA, and provides a unique opportunity to identify novel human disease genes.

AIMS

To demonstrate the power of combining the FANTOM 2 cDNA dataset with a positional candidate approach and bioinformatics analysis to identify genes underlying human genetic disease.

RESULTS

By mapping all FANTOM 2 cDNA to the human genome, we were able to identify mouse clones that co-localised on the human genome with mapped but uncloned human disease loci. By this method we identified mouse and corresponding human genes mapping within the loci of 100 different human genetic diseases (mapped interval of <5 cM). Of particular interest was the elucidation through FANTOM 2 novel mouse gene data of candidate human genes for the following: (i) developmental -disorders: neural tube defect, Meckel syndrome, Wolf--Hirschhorn syndrome and keratosis follicularis spinulosa decalvans cum ophiasi; (ii) neurological disorders: benign familial infantile convulsions 3, early-onset cerebellar ataxia with retained tendon reflexes, infantile-onset spinocerebellar ataxia and vacuolar neuro-myopathy and (iii) cancer-related syndromes: tylosis with oesophageal cancer and low-grade B-cell chronic lymphatic leukaemia.

CONCLUSIONS

The FANTOM 2 data will dramatically accelerate efforts to identify genes underlying human disease. It will also facilitate the creation of transgenic mouse models to help elucidate the function of potential human disease genes.

Numerous novel annotations of the human genome sequence supported by a 5'-end-enriched cDNA collection

A collection of 90,000 human cDNA clones generated to increase the fraction of "full-length" cDNAs available was analyzed by sequence alignment on the human genome assembly. Five hundred fifty-two gene models not found in LocusLink, with coding regions of at least 300 bp, were defined by using this collection. Exon composition proposed for novel genes showed an average of 4.7 exons per gene. In 20% of the cases, at least half of the exons predicted for new genes coincided with evolutionary conserved regions defined by sequence comparisons with the pufferfish Tetraodon nigroviridis. Among this subset, CpG islands were observed at the 5' end of 75%. In-frame stop codons upstream of the initiator ATG were present in 49% of the new genes, and 16% contained a coding region comprising at least 50% of the cDNA sequence. This cDNA resource also provided candidate small protein-coding genes, usually not included in genome annotations. In addition, analysis of a sample from this cDNA collection indicates that approximately 380 gene models described in LocusLink could be extended at their 5' end by at least one new exon. Finally, this cDNA resource provided an experimental support for annotations based exclusively on predictions, thus representing a resource substantially improving the human genome annotation.

Application of DNA Microarrays in Occupational Health Research

The profiling of gene expression patterns with DNA microarrays is recently being widely used not only in basic molecular biological studies but also in the practical fields. In clinical application, for example, this technique is expected to be quite useful in making a correct diagnosis. In the pharmacological area, the microarray analysis can be applied to drug discovery and individualized drug treatment. Although not so popular as these examples, DNA microarrays could also be a powerful tool in studies relevant to occupational health. This review will describe the outline of gene expression profiling with DNA microarrays and prospects in occupational health research.

Large-scale identification of disease genes involved in acute myeloid leukemia

Acute myeloid leukemia (AML) is a heterogeneous group of diseases in which chromosomal aberrations, small insertions or deletions, or point mutations in certain genes have profound consequences for prognosis. However, the majority of AML patients present without currently known genetic defects. Retroviral insertion mutagenesis in mice has become a powerful tool for identifying new disease genes involved in the pathogenesis of leukemia and lymphoma. Here we have used the Graffi-1.4 strain of murine leukemia virus, which causes predominantly AML, in a screen to identify novel genes involved in the pathogenesis of this disease. We report 79 candidate disease genes in common integration sites (CISs) and 15 genes whose family members previously were found to be affected in other studies. The majority of the identified sequences (60%) were not found in lymphomas and monocytic leukemias in previous screens, suggesting a specific involvement in AML. Although most of the virus integrations occurred in or near the 5' or 3' ends of the genes, suggesting deregulation of gene expression as a consequence of virus integration, 18 CISs were located exclusively within the genes, conceivably causing gene disruption.

Altered gene expression in acute systemic inflammation detected by complete coverage of the human liver transcriptome

The goal of the current study was to provide complete coverage of the liver transcriptome with human probes corresponding to every gene expressed in embryonic, adult, and/or cancerous liver. We developed dedicated tools, namely, the Liverpool nylon array of complementary DNA (cDNA) probes for approximately 10,000 nonredundant genes and the LiverTools database. Inflammation-induced transcriptome changes were studied in liver tissue samples from patients with an acute systemic inflammation and from control subjects. One hundred and fifty-four messenger RNAs (mRNA) correlated statistically with the extent of inflammation. Of these, 134 mRNA samples were not associated previously with an acute-phase (AP) response. The hepatocyte origin and proinflammatory cytokine responsiveness of these mRNAs were confirmed by quantitative reverse-transcription polymerase chain reaction (Q-RT-PCR) in cytokine-challenged hepatoma cells. The corresponding gene promoters were enriched in potential binding sites for inflammation-driven transcription factors in the liver. Some of the corresponding proteins may provide novel blood markers of clinical relevance. The mRNAs whose level is most correlated with the AP extent (P <.05) were enriched in intracellular signaling molecules, transcription factors, glycosylation enzymes, and up-regulated plasma proteins. In conclusion, the hepatocyte responded to the AP extent by fine tuning some mRNA levels, controlling most, if not all, intracellular events from early signaling to the final secretion of proteins involved in innate immunity. Supplementary material for this article can be found on the HEPATOLOGY website (http://interscience.wiley.com/jpages/0270-9139/suppmat/index.html).

Evidence for existence of thirty hypothetical proteins in rat brain

Background

The rapid completion of genome sequences has created an infrastructure of biological information and provided essential information to link genes to gene products, proteins, the building blocks for cellular functions. In addition, genome/cDNA sequences make it possible to predict proteins for which there is no experimental evidence. Clues for function of hypothetical proteins are provided by sequence similarity with proteins of known function in model organisms.

Results

We constructed a two-dimensional protein map and searched for expression of hypothetical proteins in rat brain. Two-dimensional electrophoresis (2-DE) with subsequent in-gel digestion of spots and matrix-assisted laser desorption/ionization (MALDI) spectrometric identification were applied. In total about 3700 spots were analysed, which resulted in the identification of about 1700 polypeptides, that were the products of 190 different genes. A number of hypothetical gene products were detected (30 of 190, 15.8%) and are considered brain proteins.

Conclusions

A major finding of this study is the demonstration of the existence of putative proteins that were so far only deduced from their nucleic acid structure by a protein chemical method independent of antibody availability and specificity and unambiguously identifying proteins.

Serial analysis of gene expression (SAGE) during porcine embryo development

Functional genomics provides a powerful means for delving into the molecular mechanisms involved in pre-implantation development of porcine embryos. High rates of embryonic mortality (30%), following either natural mating or artificial insemination, emphasise the need to improve the efficiency of reproduction in the pig. The poor success rate of live offspring from in vitro-manipulated pig embryos also hampers efforts to generate transgenic animals for biotechnology applications. Previous analysis of differential gene expression has demonstrated stage-specific gene expression for in vivo-derived embryos and altered gene expression for in vitro-derived embryos. However, the methods used to date examine relatively few genes simultaneously and, thus, provide an incomplete glimpse of the physiological role of these genes during embryogenesis. The present review will focus on two aspects of applying functional genomics research strategies for analysing the expression of genes during elongation of pig embryos between gestational day (D) 11 and D12. First, we compare and contrast current methodologies that are being used for gene discovery and expression analysis during pig embryo development. Second, we establish a paradigm for applying serial analysis of gene expression as a functional genomics tool to obtain preliminary information essential for discovering the physiological mechanisms by which distinct embryonic phenotypes are derived.

Autosomal similarity revealed by eukaryotic genomic comparison

To describe eukaryotic autosomes quantitatively and determine differences between them in terms of amino acid sequences of genes, functional classification of proteins, and complete DNA sequences, we applied two theoretical methods, the Proteome-vector method and the function of degree of disagreement (FDOD) method, that are based on function and sequence similarity respectively, to autosomes from nine eukaryotes. No matter what aspect of the autosome is considered, the autosomal differences within each organism were less than that between species. Our results show that eukaryotic autosomes resemble each other within a species while those from different organisms differ. We propose a hypothesis (named intra-species autosomal random shuffling) as an explanation for our results and suggest that lateral gene transfer (LGT) did not occur frequently during the evolution of eukarya.

An integrated gene annotation and transcriptional profiling approach towards the full gene content of the Drosophila genome

A novel Drosophila microarray constructed on the basis of an integrated in silico/wet biology approach provides evidence for the transcription of approximately 2,600 additional genes. Validation indicates a lower limit of 2,000 novel annotations, thus raising the number of genes that make a fly.

Background

While the genome sequences for a variety of organisms are now available, the precise number of the genes encoded is still a matter of debate. For the human genome several stringent annotation approaches have resulted in the same number of potential genes, but a careful comparison revealed only limited overlap. This indicates that only the combination of different computational prediction methods and experimental evaluation of such in silico data will provide more complete genome annotations. In order to get a more complete gene content of the Drosophila melanogaster genome, we based our new D. melanogaster whole-transcriptome microarray, the Heidelberg FlyArray, on the combination of the Berkeley Drosophila Genome Project (BDGP) annotation and a novel ab initio gene prediction of lower stringency using the Fgenesh software.

Results

Here we provide evidence for the transcription of approximately 2,600 additional genes predicted by Fgenesh. Validation of the developmental profiling data by RT-PCR and in situ hybridization indicates a lower limit of 2,000 novel annotations, thus substantially raising the number of genes that make a fly.

Conclusions

The successful design and application of this novel Drosophila microarray on the basis of our integrated in silico/wet biology approach confirms our expectation that in silico approaches alone will always tend to be incomplete. The identification of at least 2,000 novel genes highlights the importance of gathering experimental evidence to discover all genes within a genome. Moreover, as such an approach is independent of homology criteria, it will allow the discovery of novel genes unrelated to known protein families or those that have not been strictly conserved between species.

Whole-genome screening indicates a possible burst of formation of processed pseudogenes and Alu repeats by particular L1 subfamilies in ancestral primates

BACKGROUND

Abundant pseudogenes are a feature of mammalian genomes. Processed pseudogenes (PPs) are reverse transcribed from mRNAs. Recent molecular biological studies show that mammalian long interspersed element 1 (L1)-encoded proteins may have been involved in PP reverse transcription. Here, we present the first comprehensive analysis of human PPs using all known human genes as queries.

RESULTS

The human genome was queried and 3,664 candidate PPs were identified. The most abundant were copies of genes encoding keratin 18, glyceraldehyde-3-phosphate dehydrogenase and ribosomal protein L21. A simple method was developed to estimate the level of nucleotide substitutions (and therefore the age) of PPs. A Poisson-like age distribution was obtained with a mean age close to that of the Alu repeats, the predominant human short interspersed elements. These data suggest a nearly simultaneous burst of PP and Alu formation in the genomes of ancestral primates. The peak period of amplification of these two distinct retrotransposons was estimated to be 40-50 million years ago. Concordant amplification of certain L1 subfamilies with PPs and Alus was observed.

CONCLUSIONS

We suggest that a burst of formation of PPs and Alus occurred in the genome of ancestral primates. One possible mechanism is that proteins encoded by members of particular L1 subfamilies acquired an enhanced ability to recognize cytosolic RNAs in trans.

Conserved noncoding sequences in the grasses

As orthologous genes from related species diverge over time, some sequences are conserved in noncoding regions. In mammals, large phylogenetic footprints, or conserved noncoding sequences (CNSs), are known to be common features of genes. Here we present the first large-scale analysis of plant genes for CNSs. We used maize and rice, maximally diverged members of the grass family of monocots. Using a local sequence alignment set to deliver only significant alignments, we found one or more CNSs in the noncoding regions of the majority of genes studied. Grass genes have dramatically fewer and much smaller CNSs than mammalian genes. Twenty-seven percent of grass gene comparisons revealed no CNSs. Genes functioning in upstream regulatory roles, such as transcription factors, are greatly enriched for CNSs relative to genes encoding enzymes or structural proteins. Further, we show that a CNS cluster in an intron of the knotted1 homeobox gene serves as a site of negative regulation. We showthat CNSs in the adh1 gene do not correlate with known cis-acting sites. We discuss the potential meanings of CNSs and their value as analytical tools and evolutionary characters. We advance the idea that many CNSs function to lock-in gene regulatory decisions.

GANESH: software for customized annotation of genome regions

GANESH is a software package designed to support the genetic analysis of regions of human and other genomes. It provides a set of components that may be assembled to construct a self-updating database of DNA sequence, mapping data, and annotations of possible genome features. Once one or more remote sources of data for the target region have been identified, all sequences for that region are downloaded, assimilated, and subjected to a (configurable) set of standard database-searching and genome-analysis packages. The results are stored in compressed form in a relational database, and are updated automatically on a regular schedule so that they are always immediately available in their most up-to-date versions. A Java front-end, executed as a stand alone application or web applet, provides a graphical interface for navigating the database and for viewing the annotations. There are facilities for importing and exporting data in the format of the Distributed Annotation System (DAS), enabling a GANESH database to be used as a component of a DAS configuration. The system has been used to construct databases for about a dozen regions of human chromosomes and for three regions of mouse chromosomes.

The T cell surface--how well do we know it?

The overall degree of complexity of the T cell surface has been unclear, constraining our understanding of its biology. Using global gene expression analysis, we show that 111 of 374 genes encoding well-characterized leukocyte surface antigens are expressed by a resting cytotoxic T cell. Unexpectedly, of 97 stringently defined, T cell-specific transcripts with unknown functions that we identify, none encode proteins with the modular architecture characteristic of 80% of leukocyte surface antigens. Only two encode proteins with membrane topologies found exclusively in cell surface molecules. Our analysis indicates that the cell type-specific composition of the resting CD8+ T cell surface is now largely defined, providing an insight into the overall compositional complexity of the mammalian cell surface and a framework for formulating systematic models of T cell surface-dependent processes, such as T cell receptor triggering.

Diversity, taxonomy and evolution of medium-chain dehydrogenase/reductase superfamily

A comprehensive, structural and functional, in silico analysis of the medium-chain dehydrogenase/reductase (MDR) superfamily, including 583 proteins, was carried out by use of extensive database mining and the blastp program in an iterative manner to identify all known members of the superfamily. Based on phylogenetic, sequence, and functional similarities, the protein members of the MDR superfamily were classified into three different taxonomic categories: (a) subfamilies, consisting of a closed group containing a set of ideally orthologous proteins that perform the same function; (b) families, each comprising a cluster of monophyletic subfamilies that possess significant sequence identity among them and might share or not common substrates or mechanisms of reaction; and (c) macrofamilies, each comprising a cluster of monophyletic protein families with protein members from the three domains of life, which includes at least one subfamily member that displays activity related to a very ancient metabolic pathway. In this context, a superfamily is a group of homologous protein families (and/or macrofamilies) with monophyletic origin that shares at least a barely detectable sequence similarity, but showing the same 3D fold. The MDR superfamily encloses three macrofamilies, with eight families and 49 subfamilies. These subfamilies exhibit great functional diversity including noncatalytic members with different subcellular, phylogenetic, and species distributions. This results from constant enzymogenesis and proteinogenesis within each kingdom, and highlights the huge plasticity that MDR superfamily members possess. Thus, through evolution a great number of taxa-specific new functions were acquired by MDRs. The generation of new functions fulfilled by proteins, can be considered as the essence of protein evolution. The mechanisms of protein evolution inside MDR are not constrained to conserve substrate specificity and/or chemistry of catalysis. In consequence, MDR functional diversity is more complex than sequence diversity. MDR is a very ancient protein superfamily that existed in the last universal common ancestor. It had at least two (and probably three) different ancestral activities related to formaldehyde metabolism and alcoholic fermentation. Eukaryotic members of this superfamily are more related to bacterial than to archaeal members; horizontal gene transfer among the domains of life appears to be a rare event in modern organisms.

ORFeome-based search of airway epithelial cell-specific novel human [beta]-defensin genes

beta-Defensin is one of the major host defense shields produced by various tissues and organs against microbial infection. To date, four human beta-defensins (DEFBs) gene products that share a consensus six-cysteine motif have been discovered. The hidden Markov model (HMM) profile was constructed from the common features of those known beta-defensin peptides to search for additional novel DEFB genes. A genome-wide search of the profile against ORFeome-based peptide databases (e.g., Ensembl project) led to the identification of six new DEFB members that also shared the conserved six-cysteine motif. Phylogenetic analysis supported a close relationship of these six new members with existing DEFB genes. Polymerase Chain Reaction studies of human tissue cDNA panels confirmed the expression of all six novel DEFB genes in various tissues. Two of them, DEFB106 and DEFB109, were expressed in the lung. A pilot study with cRNA probes for in situ hybridization and a synthetic propeptide for the functional characterization demonstrated the tissue-/cell-specific expression and the strong antimicrobial activity of DEFB106. These results support the utility of ORFeome-based HMM search in gene discovery for members of a specific gene family. The novel DEFB genes identified in this study may significantly contribute to overall antimicrobial host defenses.

TMC and EVER genes belong to a larger novel family, the TMC gene family encoding transmembrane proteins

BACKGROUND

Mutations in the transmembrane cochlear expressed gene 1 (TMC1) cause deafness in human and mouse. Mutations in two homologous genes, EVER1 and EVER2 increase the susceptibility to infection with certain human papillomaviruses resulting in high risk of skin carcinoma. Here we report that TMC1, EVER1 and EVER2 (now TMC6 and TMC8) belong to a larger novel gene family, which is named TMC for trans membrane channel-like gene family.

RESULTS

Using a combination of iterative database searches and reverse transcriptase-polymerase chain reaction (RT-PCR) experiments we assembled contigs for cDNA encoding human, murine, puffer fish, and invertebrate TMC proteins. TMC proteins of individual species can be grouped into three subfamilies A, B, and C. Vertebrates have eight TMC genes. The majority of murine TMC transcripts are expressed in most organs; some transcripts, however, in particular the three subfamily A members are rare and more restrictively expressed.

CONCLUSION

The eight vertebrate TMC genes are evolutionary conserved and encode proteins that form three subfamilies. Invertebrate TMC proteins can also be categorized into these three subfamilies. All TMC genes encode transmembrane proteins with intracellular amino- and carboxyl-termini and at least eight membrane-spanning domains. We speculate that the TMC proteins constitute a novel group of ion channels, transporters, or modifiers of such.

Human-mouse gene identification by comparative evidence integration and evolutionary analysis

The identification of genes in the human genome remains a challenge, as the actual predictions appear to disagree tremendously and vary dramatically on the basis of the specific gene-finding methodology used. Because the pattern of conservation in coding regions is expected to be different from intronic or intergenic regions, a comparative computational analysis can lead, in principle, to an improved computational identification of genes in the human genome by using a reference, such as mouse genome. However, this comparative methodology critically depends on three important factors: (1) the selection of the most appropriate reference genome. In particular, it is not clear whether the mouse is at the correct evolutionary distance from the human to provide sufficiently distinctive conservation levels in different genomic regions, (2) the selection of comparative features that provide the most benefit to gene recognition, and (3) the selection of evidence integration architecture that effectively interprets the comparative features. We address the first question by a novel evolutionary analysis that allows us to explicitly correlate the performance of the gene recognition system with the evolutionary distance (time) between the two genomes. Our simulation results indicate that there is a wide range of reference genomes at different evolutionary time points that appear to deliver reasonable comparative prediction of human genes. In particular, the evolutionary time between human and mouse generally falls in the region of good performance; however, better accuracy might be achieved with a reference genome further than mouse. To address the second question, we propose several natural comparative measures of conservation for identifying exons and exon boundaries. Finally, we experiment with Bayesian networks for the integration of comparative and compositional evidence.

Trapping and sequence analysis of 1138 putative exons from human chromosome 18

In a search for novel genes on chromosome 18 (HC18), on which several regions have been linked to bipolar disorder, we applied exon trapping to HC18-specific cosmids. Among the 1138 exons trapped, 1052 of them have been mapped to HC18, and the remaining 86 have not been localized. No exons were localized to genomic regions other than HC18. BLAST database search revealed that 190 exons were identical to 98 Unigenes on HC18; 98 identical to additional 82 clusters of ESTs not present in the HC18 Unigene set; 39 homologous to genes from human and other species (e<10(-3)); and the remaining 811 exons had no significant homology to transcripts in public databases. The mapped exons were compared to the 867 annotated genes on HC18 in the Celera databases; 216 exons were identical to 104 Celera 'genes' and the remaining 836 exons were not found in the Celera databases. On average, there were two exons for a matched transcript (known genes and ESTs). Therefore, the 850 novel exons may represent hundreds of novel genes on chromosome 18.

The Dynamics of Chromosome Organization and Gene Regulation

With the sequence of the human genome now complete, studies must focus on how the genome is functionally organized within the confines of the cell nucleus and the dynamic interplay between the genome and its regulatory factors to effectively control gene expression and silencing. In this review I describe our current state of knowledge with regard to the organization of chromosomes within the nucleus and the positioning of active versus inactive genes. In addition, I discuss studies on the dynamics of chromosomes and specific genetic loci within living cells and its relationship to gene activity and the cell cycle. Furthermore, our current understanding of the distribution and dynamics of RNA polymerase II transcription factors is discussed in relation to chromosomal loci and other nuclear domains.