OverGeneDB: a database of 5' end protein coding overlapping genes in human and mouse genomes

Gene overlap plays various regulatory functions on transcriptional and post-transcriptional levels. Most current studies focus on protein-coding genes overlapping with non-protein-coding counterparts, the so called natural antisense transcripts. Considerably less is known about the role of gene overlap in the case of two protein-coding genes. Here, we provide OverGeneDB, a database of human and mouse 5′ end protein-coding overlapping genes. The database contains 582 human and 113 mouse gene pairs that are transcribed using overlapping promoters in at least one analyzed library. Gene pairs were identified based on the analysis of the transcription start site (TSS) coordinates in 73 human and 10 mouse organs, tissues and cell lines. Beside TSS data, resources for 26 human lung adenocarcinoma cell lines also contain RNA-Seq and ChIP-Seq data for seven histone modifications and RNA Polymerase II activity. The collected data revealed that the overlap region is rarely conserved between the studied species and tissues. In ∼50% of the overlapping genes, transcription started explicitly in the overlap regions. In the remaining half of overlapping genes, transcription was initiated both from overlapping and non-overlapping TSSs. OverGeneDB is accessible at http://overgenedb.amu.edu.pl.

Overactive BRCA1 Affects Presenilin 1 in Induced Pluripotent Stem Cell-Derived Neurons in Alzheimer's Disease

The BRCA1 protein, one of the major players responsible for DNA damage response has recently been linked to Alzheimer's disease (AD). Using primary fibroblasts and neurons reprogrammed from induced pluripotent stem cells (iPSC) derived from familial AD (FAD) patients, we studied the role of the BRCA1 protein underlying molecular neurodegeneration. By whole-transcriptome approach, we have found wide range of disturbances in cell cycle and DNA damage response in FAD fibroblasts. This was manifested by significantly increased content of BRCA1 phosphorylated on Ser1524 and abnormal ubiquitination and subcellular distribution of presenilin 1 (PS1). Accordingly, the iPSC-derived FAD neurons showed increased content of BRCA1(Ser1524) colocalized with degraded PS1, accompanied by an enhanced immunostaining pattern of amyloid-β. Finally, overactivation of BRCA1 was followed by an increased content of Cdc25C phosphorylated on Ser216, likely triggering cell cycle re-entry in FAD neurons. This study suggests that overactivated BRCA1 could both influence PS1 turnover leading to amyloid-β pathology and promote cell cycle re-entry-driven cell death of postmitotic neurons in AD.

Comparison of Highly and Weakly Virulent Dickeya solani Strains, With a View on the Pangenome and Panregulon of This Species

Bacteria belonging to the genera Dickeya and Pectobacterium are responsible for significant economic losses in a wide variety of crops and ornamentals. During last years, increasing losses in potato production have been attributed to the appearance of Dickeya solani. The D. solani strains investigated so far share genetic homogeneity, although different virulence levels were observed among strains of various origins. The purpose of this study was to investigate the genetic traits possibly related to the diverse virulence levels by means of comparative genomics. First, we developed a new genome assembly pipeline which allowed us to complete the D. solani genomes. Four de novo sequenced and ten publicly available genomes were used to identify the structure of the D. solani pangenome, in which 74.8 and 25.2% of genes were grouped into the core and dispensable genome, respectively. For D. solani panregulon analysis, we performed a binding site prediction for four transcription factors, namely CRP, KdgR, PecS and Fur, to detect the regulons of these virulence regulators. Most of the D. solani potential virulence factors were predicted to belong to the accessory regulons of CRP, KdgR, and PecS. Thus, some differences in gene expression could exist between D. solani strains. The comparison between a highly and a low virulent strain, IFB0099 and IFB0223, respectively, disclosed only small differences between their genomes but significant differences in the production of virulence factors like pectinases, cellulases and proteases, and in their mobility. The D. solani strains also diverge in the number and size of prophages present in their genomes. Another relevant difference is the disruption of the adhesin gene fhaB2 in the highly virulent strain. Strain IFB0223, which has a complete adhesin gene, is less mobile and less aggressive than IFB0099. This suggests that in this case, mobility rather than adherence is needed in order to trigger disease symptoms. This study highlights the utility of comparative genomics in predicting D. solani traits involved in the aggressiveness of this emerging plant pathogen.

Transcriptional interference by small transcripts in proximal promoter regions

Proximal promoter regions (PPR) are heavily transcribed yielding different types of small RNAs. The act of transcription within PPRs might regulate downstream gene expression via transcriptional interference (TI). For analysis, we investigated capped and polyadenylated small RNA transcripts within PPRs of human RefSeq genes in eight different cell lines. Transcripts of our datasets overlapped with experimentally determined transcription factor binding sites (TFBS). For TFBSs intersected by these small RNA transcripts, we established negative correlation of sRNA expression levels and transcription factor (TF) DNA binding affinities; suggesting that the transcripts acted via TI. Accordingly, datasets were designated as TFbiTrs (TF-binding interfering transcripts). Expression of most TFbiTrs was restricted to certain cell lines. This facilitated the analysis of effects related to TFbiTr expression for the same RefSeq genes across cell lines. We consistently uncovered higher relative TF/DNA binding affinities and concomitantly higher expression levels for RefSeq genes in the absence of TFbiTrs. Analysis of corresponding chromatin landscapes supported these results. ChIA-PET revealed the participation of distal enhancers in TFbiTr transcription. Enhancers regulating TFbiTrs, in effect, act as repressors for corresponding downstream RefSeq genes. We demonstrate the significant impact of TI on gene expression using selected small RNA datasets.

RetrogeneDB-a database of plant and animal retrocopies

For a long time, retrocopies were considered ‘junk DNA’, but numerous studies have shown that retrocopies may gain functionality and become so-called retrogenes. Retrogenes may code fully functional proteins that coexist with parental gene products or may even replace them. Retrocopies may also function as regulatory RNAs and, for example, become a source of small interfering RNAs, act as trans natural antisense transcripts or as alternative targets for miRNAs. Numerous researchers have emphasized that retrogenes play a crucial role in various organisms’ developmental stages and diseases. Despite the ever-growing evidence of the importance of retrocopies, resources dedicated to retroposition are very limited. Here, we report an update of the RetrogeneDB, which, to the best of our knowledge, is the largest database dedicated to retrocopies. It provides annotations of 86 458 retrocopies in 62 animal and 37 plant species. The database contains information about the retrocopies’ localization, open reading frame conservation, expression, RNA Polymerase II activity and the alternative transcription start site studies. Orthologous relationships between retrogenes were also determined, which made retrocopy conservation studies much more valuable. Additionally, based on the RNA-Seq data from the Geuvadis project, the expression levels of retrocopies were estimated in a total of 50 individuals from 5 human populations. The information is now presented in a new, more user-friendly web interface, with easy access to the source data, which may be used for the downstream analysis. RetrogeneDB is freely available at http://yeti.amu.edu.pl/retrogenedb.

Database URL:http://yeti.amu.edu.pl/retrogenedb

Secondary database URL:http://rhesus.amu.edu.pl/retrogenedb

Rootstock-regulated gene expression patterns associated with fire blight resistance in apple

BACKGROUND

Desirable apple varieties are clonally propagated by grafting vegetative scions onto rootstocks. Rootstocks influence many phenotypic traits of the scion, including resistance to pathogens such as Erwinia amylovora, which causes fire blight, the most serious bacterial disease of apple. The purpose of the present study was to quantify rootstock-mediated differences in scion fire blight susceptibility and to identify transcripts in the scion whose expression levels correlated with this response.

RESULTS

Rootstock influence on scion fire blight resistance was quantified by inoculating three-year old, orchard-grown apple trees, consisting of 'Gala' scions grafted to a range of rootstocks, with E. amylovora. Disease severity was measured by the extent of shoot necrosis over time. 'Gala' scions grafted to G.30 or MM.111 rootstocks showed the lowest rates of necrosis, while 'Gala' on M.27 and B.9 showed the highest rates of necrosis. 'Gala' scions on M.7, S.4 or M.9F56 had intermediate necrosis rates. Using an apple DNA microarray representing 55,230 unique transcripts, gene expression patterns were compared in healthy, un-inoculated, greenhouse-grown 'Gala' scions on the same series of rootstocks. We identified 690 transcripts whose steady-state expression levels correlated with the degree of fire blight susceptibility of the scion/rootstock combinations. Transcripts known to be differentially expressed during E. amylovora infection were disproportionately represented among these transcripts. A second-generation apple microarray representing 26,000 transcripts was developed and was used to test these correlations in an orchard-grown population of trees segregating for fire blight resistance. Of the 690 transcripts originally identified using the first-generation array, 39 had expression levels that correlated with fire blight resistance in the breeding population.

CONCLUSIONS

Rootstocks had significant effects on the fire blight susceptibility of 'Gala' scions, and rootstock-regulated gene expression patterns could be correlated with differences in susceptibility. The results suggest a relationship between rootstock-regulated fire blight susceptibility and sorbitol dehydrogenase, phenylpropanoid metabolism, protein processing in the endoplasmic reticulum, and endocytosis, among others. This study illustrates the utility of our rootstock-regulated gene expression data sets for candidate trait-associated gene data mining.

miRNEST database: an integrative approach in microRNA search and annotation

Despite accumulating data on animal and plant microRNAs and their functions, existing public miRNA resources usually collect miRNAs from a very limited number of species. A lot of microRNAs, including those from model organisms, remain undiscovered. As a result there is a continuous need to search for new microRNAs. We present miRNEST (http://mirnest.amu.edu.pl), a comprehensive database of animal, plant and virus microRNAs. The core part of the database is built from our miRNA predictions conducted on Expressed Sequence Tags of 225 animal and 202 plant species. The miRNA search was performed based on sequence similarity and as many as 10 004 miRNA candidates in 221 animal and 199 plant species were discovered. Out of them only 299 have already been deposited in miRBase. Additionally, miRNEST has been integrated with external miRNA data from literature and 13 databases, which includes miRNA sequences, small RNA sequencing data, expression, polymorphisms and targets data as well as links to external miRNA resources, whenever applicable. All this makes miRNEST a considerable miRNA resource in a sense of number of species (544) that integrates a scattered miRNA data into a uniform format with a user-friendly web interface.

Phytophthora Database: A Forensic Database Supporting the Identification and Monitoring of Phytophthora

Phytophthora spp. represent a serious threat to agricultural and ecological systems. Many novel Phytophthora spp. have been reported in recent years, which is indicative of our limited understanding of the ecology and diversity of Phytophthora spp. in nature. Systematic cataloging of genotypic and phenotypic information on isolates of previously described species serves as a baseline for identification, classification, and risk assessment of new Phytophthora isolates. The Phytophthora Database (PD) was established to catalog such data in a web-accessible and searchable format. To support the identification of new Phytophthora isolates via comparison of their sequences at one or more loci with the corresponding sequences derived from the isolates archived in the PD, we generated and deposited sequence data from more than 1,500 isolates representing the known diversity in the genus. Data search and analysis tools in the PD include BLAST, Phyloviewer (a program for building phylogenetic trees using sequences of selected isolates), and Virtual Gel (a program for generating expected restriction patterns for given sequences). The PD also provides a customized means of storing and sharing data via the web. The PD serves as a model that easily can be adopted to develop databases for other important pathogen groups.

The H-Invitational Database (H-InvDB), a comprehensive annotation resource for human genes and transcripts

Here we report the new features and improvements in our latest release of the H-Invitational Database (H-InvDB; http://www.h-invitational.jp/), a comprehensive annotation resource for human genes and transcripts. H-InvDB, originally developed as an integrated database of the human transcriptome based on extensive annotation of large sets of full-length cDNA (FLcDNA) clones, now provides annotation for 120 558 human mRNAs extracted from the International Nucleotide Sequence Databases (INSD), in addition to 54 978 human FLcDNAs, in the latest release H-InvDB_4.6. We mapped those human transcripts onto the human genome sequences (NCBI build 36.1) and determined 34 699 human gene clusters, which could define 34 057 (98.1%) protein-coding and 642 (1.9%) non-protein-coding loci; 858 (2.5%) transcribed loci overlapped with predicted pseudogenes. For all these transcripts and genes, we provide comprehensive annotation including gene structures, gene functions, alternative splicing variants, functional non-protein-coding RNAs, functional domains, predicted sub cellular localizations, metabolic pathways, predictions of protein 3D structure, mapping of SNPs and microsatellite repeat motifs, co-localization with orphan diseases, gene expression profiles, orthologous genes, protein-protein interactions (PPI) and annotation for gene families. The current H-InvDB annotation resources consist of two main views: Transcript view and Locus view and eight sub-databases: the DiseaseInfo Viewer, H-ANGEL, the Clustering Viewer, G-integra, the TOPO Viewer, Evola, the PPI view and the Gene family/group.

Plant pathogen culture collections: it takes a village to preserve these resources vital to the advancement of agricultural security and plant pathology

ABSTRACT Plant pathogen culture collections are essential resources in our fight against plant disease and for connecting discoveries of the present with established knowledge of the past. However, available infrastructure in support of culture collections is in serious need of improvement, and we continually face the risk of losing many of these collections. As novel and reemerging plant pathogens threaten agriculture, their timely identification and monitoring depends on rapid access to cultures representing the known diversity of plant pathogens along with genotypic, phenotypic, and epidemiological data associated with them. Archiving such data in a format that can be easily accessed and searched is essential for rapid assessment of potential risk and can help track the change and movement of pathogens. The underexplored pathogen diversity in nature further underscores the importance of cataloguing pathogen cultures. Realizing the potential of pathogen genomics as a foundation for developing effective disease control also hinges on how effectively we use the sequenced isolate as a reference to understand the genetic and phenotypic diversity within a pathogen species. In this letter, we propose a number of measures for improving pathogen culture collections.

ChloroplastDB: the Chloroplast Genome Database

The Chloroplast Genome Database (ChloroplastDB) is an interactive, web-based database for fully sequenced plastid genomes, containing genomic, protein, DNA and RNA sequences, gene locations, RNA-editing sites, putative protein families and alignments (http://chloroplast.cbio.psu.edu/). With recent technical advances, the rate of generating new organelle genomes has increased dramatically. However, the established ontology for chloroplast genes and gene features has not been uniformly applied to all chloroplast genomes available in the sequence databases. For example, annotations for some published genome sequences have not evolved with gene naming conventions. ChloroplastDB provides unified annotations, gene name search, BLAST and download functions for chloroplast encoded genes and genomic sequences. A user can retrieve all orthologous sequences with one search regardless of gene names in GenBank. This feature alone greatly facilitates comparative research on sequence evolution including changes in gene content, codon usage, gene structure and post-transcriptional modifications such as RNA editing. Orthologous protein sets are classified by TribeMCL and each set is assigned a standard gene name. Over the next few years, as the number of sequenced chloroplast genomes increases rapidly, the tools available in ChloroplastDB will allow researchers to easily identify and compile target data for comparative analysis of chloroplast genes and genomes.

SLC24A5, a putative cation exchanger, affects pigmentation in zebrafish and humans

Lighter variations of pigmentation in humans are associated with diminished number, size, and density of melanosomes, the pigmented organelles of melanocytes. Here we show that zebrafish golden mutants share these melanosomal changes and that golden encodes a putative cation exchanger slc24a5 (nckx5) that localizes to an intracellular membrane, likely the melanosome or its precursor. The human ortholog is highly similar in sequence and functional in zebrafish. The evolutionarily conserved ancestral allele of a human coding polymorphism predominates in African and East Asian populations. In contrast, the variant allele is nearly fixed in European populations, is associated with a substantial reduction in regional heterozygosity, and correlates with lighter skin pigmentation in admixed populations, suggesting a key role for the SLC24A5 gene in human pigmentation.

Overlapping genes in vertebrate genomes

Overlapping genes in mammalian genomes are unexpected phenomena even though hundreds of pairs of protein coding overlapping genes have been reported so far. Overlapping genes can be divided into different categories based on direction of transcription as well as on sequence segments being shared between overlapping coding regions. The biologic functions of natural antisense transcripts, their involvement in physiological processes and gene regulation in living organisms are not fully understood. Number of documented examples indicates that they may exert control at various levels of gene expression, such as transcription, mRNA processing, splicing, stability, transport, and translation. Similarly, evolutionary origin of such genes is not known, existing hypotheses can explain only selected cases of mammalian gene overlaps which could originate as result of rearrangements, overprinting and/or adoption of signals in the neighboring gene locus.

Origin and evolution of the chicken leukocyte receptor complex

In mammals, the cell surface receptors encoded by the leukocyte receptor complex (LRC) regulate the activity of T lymphocytes and B lymphocytes, as well as that of natural killer cells, and thus provide protection against pathogens and parasites. The chicken genome encodes many Ig-like receptors that are homologous to the LRC receptors. The chicken Ig-like receptor (CHIR) genes are members of a large monophyletic gene family and are organized into genomic clusters, which are in conserved synteny with the mammalian LRC. One-third of CHIR genes encode polypeptide molecules that contain both activating and inhibitory motifs. These genes are present in different phylogenetic groups, suggesting that the primordial CHIR gene could have encoded both types of motifs in a single molecule. In contrast to the mammalian LRC genes, the CHIR genes with similar function (inhibition or activation) are evolutionarily closely related. We propose that, in addition to recombination, single nucleotide substitutions played an important role in the generation of receptors with different functions. Structural models and amino acid analyses of the CHIR proteins reveal the presence of different types of Ig-like domains in the same phylogenetic groups, as well as sharing of conserved residues and conserved changes of residues between different CHIR groups and between CHIRs and LRCs. Our data support the notion that the CHIR gene clusters are regions homologous to the mammalian LRC gene cluster and favor a model of evolution by repeated processes of birth and death (expansion-contraction) of the Ig-like receptor genes.

Integrative annotation of 21,037 human genes validated by full-length cDNA clones

The human genome sequence defines our inherent biological potential; the realization of the biology encoded therein requires knowledge of the function of each gene. Currently, our knowledge in this area is still limited. Several lines of investigation have been used to elucidate the structure and function of the genes in the human genome. Even so, gene prediction remains a difficult task, as the varieties of transcripts of a gene may vary to a great extent. We thus performed an exhaustive integrative characterization of 41,118 full-length cDNAs that capture the gene transcripts as complete functional cassettes, providing an unequivocal report of structural and functional diversity at the gene level. Our international collaboration has validated 21,037 human gene candidates by analysis of high-quality full-length cDNA clones through curation using unified criteria. This led to the identification of 5,155 new gene candidates. It also manifested the most reliable way to control the quality of the cDNA clones. We have developed a human gene database, called the H-Invitational Database (H-InvDB; http://www.h-invitational.jp/). It provides the following: integrative annotation of human genes, description of gene structures, details of novel alternative splicing isoforms, non-protein-coding RNAs, functional domains, subcellular localizations, metabolic pathways, predictions of protein three-dimensional structure, mapping of known single nucleotide polymorphisms (SNPs), identification of polymorphic microsatellite repeats within human genes, and comparative results with mouse full-length cDNAs. The H-InvDB analysis has shown that up to 4% of the human genome sequence (National Center for Biotechnology Information build 34 assembly) may contain misassembled or missing regions. We found that 6.5% of the human gene candidates (1,377 loci) did not have a good protein-coding open reading frame, of which 296 loci are strong candidates for non-protein-coding RNA genes. In addition, among 72,027 uniquely mapped SNPs and insertions/deletions localized within human genes, 13,215 nonsynonymous SNPs, 315 nonsense SNPs, and 452 indels occurred in coding regions. Together with 25 polymorphic microsatellite repeats present in coding regions, they may alter protein structure, causing phenotypic effects or resulting in disease. The H-InvDB platform represents a substantial contribution to resources needed for the exploration of human biology and pathology.

Identification and characterization of mouse Rab32 by mRNA and protein expression analysis

Rab proteins, a subfamily of the ras superfamily, are low molecular weight GTPases involved in the regulation of intracellular vesicular transport. Cloning of human RAB32 was recently described. Presently, we report the cloning and characterization of the mouse homologue of Rab32. We show that murine Rab32 exhibits a ubiquitous expression pattern, with tissue-specific variation in expression level. Three cell types with highly specialized organelles, melanocytes, platelets and mast cells, exhibit relatively high level of Rab32. We show that in murine amelanotic in vitro transformed melanocytes as well as in human amelanotic metastatic melanoma cell lines, the expression of Rab32 is markedly reduced or absent, in parallel with the loss of expression of two key enzymes for the production of melanin, tyrosinase and Tyrp1. Therefore, in both mouse and human systems, the expression of Rab32 correlates with the expression of genes involved in pigment production. However, in melanoma samples, amelanotic due to a mutation in the tyrosinase gene, the expression of Rab32 remains at levels comparable to those observed in pigmented melanoma samples. Finally, we observed co-localization of Rab32 and the melanosomal proteins, Tyrp1 and Dct, indicating an association of Rab32 with melanosomes. Based on these data, we propose the inclusion of Rab32 to the so-called melanocyte/platelet family of Rab proteins.

Cloning and characterization of a novel gene, SHPRH, encoding a conserved putative protein with SNF2/helicase and PHD-finger domains from the 6q24 region

Here we report the identification of a novel transcript containing SNF2, PHD-finger, RING-finger, helicase, and linker histone domains mapping to the q24 band region of human chromosome 6. These domains are characteristic of several DNA repair proteins, transcription factors, and helicases. We have cloned both human and mouse homologs of this novel gene using interexon PCR and RACE technologies. The human cDNA, termed SHPRH, is 6018 bp and codes for a putative protein of 1683 amino acids. The mouse cDNA, termed Shprh, is 7225 bp and codes for a putative protein of 1616 amino acids. The deduced amino acid sequences of the two proteins share 86% identity. Both genes are expressed ubiquitously, with a transcript size of approximately 7.5 kb. Mapping of this gene to 6q24, a region reported to contain a tumor suppressor locus, prompted us to evaluate SHPRH by mutation analysis in tumor cell lines. We have identified one truncating and three missense mutations, thus suggesting SHPRH as a possible candidate for the tumor suppressor gene.

Melanoma mouse model implicates metabotropic glutamate signaling in melanocytic neoplasia

To gain insight into melanoma pathogenesis, we characterized an insertional mouse mutant, TG3, that is predisposed to develop multiple melanomas. Physical mapping identified multiple tandem insertions of the transgene into intron 3 of Grm1 (encoding metabotropic glutamate receptor 1) with concomitant deletion of 70 kb of intronic sequence. To assess whether this insertional mutagenesis event results in alteration of transcriptional regulation, we analyzed Grm1 and two flanking genes for aberrant expression in melanomas from TG3 mice. We observed aberrant expression of only Grm1. Although we did not detect its expression in normal mouse melanocytes, Grm1 was ectopically expressed in the melanomas from TG3 mice. To confirm the involvement of Grm1 in melanocytic neoplasia, we created an additional transgenic line with Grm1 expression driven by the dopachrome tautomerase promoter. Similar to the original TG3, the Tg(Grm1)EPv line was susceptible to melanoma. In contrast to human melanoma, these transgenic mice had a generalized hyperproliferation of melanocytes with limited transformation to fully malignant metastasis. We detected expression of GRM1 in a number of human melanoma biopsies and cell lines but not in benign nevi and melanocytes. This study provides compelling evidence for the importance of metabotropic glutamate signaling in melanocytic neoplasia.

GALA, a database for genomic sequence alignments and annotations

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

Mutant deoxynucleotide carrier is associated with congenital microcephaly

The disorder Amish microcephaly (MCPHA) is characterized by severe congenital microcephaly, elevated levels of alpha-ketoglutarate in the urine and premature death. The disorder is inherited in an autosomal recessive pattern and has been observed only in Old Order Amish families whose ancestors lived in Lancaster County, Pennsylvania. Here we show, by using a genealogy database and automated pedigree software, that 23 nuclear families affected with MCPHA are connected to a single ancestral couple. Through a whole-genome scan, fine mapping and haplotype analysis, we localized the gene affected in MCPHA to a region of 3 cM, or 2 Mb, on chromosome 17q25. We constructed a map of contiguous genomic clones spanning this region. One of the genes in this region, SLC25A19, which encodes a nuclear mitochondrial deoxynucleotide carrier (DNC), contains a substitution that segregates with the disease in affected individuals and alters an amino acid that is highly conserved in similar proteins. Functional analysis shows that the mutant DNC protein lacks the normal transport activity, implying that failed deoxynucleotide transport across the inner mitochondrial membrane causes MCPHA. Our data indicate that mitochondrial deoxynucleotide transport may be essential for prenatal brain growth.

Identification of six novel genes by experimental validation of GeneMachine predicted genes

In silico gene identification from finished and unfinished human genome sequence has become critically important in many projects seeking to gain insights into the gene content of genomic regions implicated in diseases. To establish limitations and criteria for in silico gene identification, and to identify novel genes of potential relevance to human prostate cancer and melanoma, 3 Mb of chromosome 1 sequence have been analyzed using GeneMachine. This program is a software suite comprising of sequence similarity programs and four gene identification programs. A total of 49 potential transcripts were selected and 37 of them were selected for experimental validation. We verified 16 of the predicted genes by experimental analysis. The comparison of the predicted transcripts with their cloned forms helped to refine predicted gene models as well as to identify splice variants for several of them. Although sequences matching with ten of our verified genes have been recently deposited in the GenBank, six of them remain novel. Our studies support the feasibility of identifying novel genes from regions of interest using draft human genome sequence.

Germline mutations in the ribonuclease L gene in families showing linkage with HPC1

Although prostate cancer is the most common non-cutaneous malignancy diagnosed in men in the United States, little is known about inherited factors that influence its genetic predisposition. Here we report that germline mutations in the gene encoding 2'-5'-oligoadenylate(2-5A)-dependent RNase L (RNASEL) segregate in prostate cancer families that show linkage to the HPC1 (hereditary prostate cancer 1) region at 1q24-25 (ref. 9). We identified RNASEL by a positional cloning/candidate gene method, and show that a nonsense mutation and a mutation in an initiation codon of RNASEL segregate independently in two HPC1-linked families. Inactive RNASEL alleles are present at a low frequency in the general population. RNASEL regulates cell proliferation and apoptosis through the interferon-regulated 2-5A pathway and has been suggested to be a candidate tumor suppressor gene. We found that microdissected tumors with a germline mutation showed loss of heterozygosity and loss of RNase L protein, and that RNASEL activity was reduced in lymphoblasts from heterozyogous individuals compared with family members who were homozygous with respect to the wildtype allele. Thus, germline mutations in RNASEL may be of diagnostic value, and the 2-5A pathway might provide opportunities for developing therapies for those with prostate cancer.

Isolation and characterization of the human homeobox gene HOX D1

Homeobox genes, first identified in Drosophila, encode transcription factors that regulate embryonic development along the anteroposterior axis of an organism. Vertebrate homeobox genes are described on the basis of their homology to the genes found within the Drosophila Antennapedia and Bithorax homeotic gene complexes. Mammals possess four paralogous homeobox (HOX) gene clusters, HOX A, HOX B, HOX C and HOX D, each located on different chromosomes, consisting of 9 to 11 genes arranged in tandem. We report the characterization of the human HOX D1 gene. This gene consists of two exons, encoding a 328 amino acid protein, separated by an intron of 354 bp. The human HOX D1 protein is one amino acid longer (328 amino acids) than the mouse protein (327 amino acids) and is 82% identical to the mouse HOX D1 homolog. The DNA binding homeodomain region of the human protein exhibits a 97% and 80% identity between mouse Hoxd1 and Drosophila labial homeodomains, respectively. The exon/intron and intron/exon splice junctions are conserved in position between human and mouse genes. Determination of the human HOX D1 gene structure permits the use of PCR based analysis of this gene for the assessment of mutations, for diseases that link to the HOXD cluster (such as Duanes Retraction Syndrome (DRS)), or polymorphisms associated with human variation. Molecular characterization of the HOXD1 gene may also permit analysis of the functional role of this gene in human neurogenisis.

Physical and transcript map of the hereditary prostate cancer region at xq27

We have recently mapped a locus for hereditary prostate cancer (termed HPCX) to the long arm of the X chromosome (Xq25-q27) through a genome-wide linkage study. Here we report the construction of an approximately 9-Mb sequence-ready bacterial clone contig map of Xq26.3-q27.3. The contig was constructed by screening BAC/PAC libraries with markers spaced at approximately 85-kb intervals. We identified overlapping clones by end-sequencing framework clones to generate 407 new sequence-tagged sites, followed by PCR verification of overlaps. Contig assembly was based on clone restriction fingerprinting and the landmark information. We identified a minimal overlap contig for genomic sequencing, which has yielded 7.7 Mb of finished sequence and 1.5 Mb of draft sequence. The transcriptional mapping effort localized 57 known and predicted genes by database searching, STS content mapping, and sequencing, followed by sequence annotation. These transcriptional units represent candidate genes for HPCX and multiple other hereditary diseases at Xq26.3-q27.3.

The Histone Database

Histone proteins are often noted for their high degree of sequence conservation. It is less often recognized that the histones are a heterogeneous protein family. Furthermore, several classes of non-histone proteins containing the histone fold motif exist. Novel histone and histone fold protein sequences continue to be added to public databases every year. The Histone Database (http://genome.nhgri.nih.gov/histones/) is a searchable, periodically updated collection of histone fold-containing sequences derived from sequence-similarity searches of public databases. Sequence sets are presented in redundant and non-redundant FASTA form, hotlinked to GenBank sequence files. Partial sequences are also now included in the database, which has considerably augmented its taxonomic coverage. Annotated alignments of full-length non-redundant sets of sequences are now available in both web-viewable (HTML) and downloadable (PDF) formats. The database also provides summaries of current information on solved histone fold structures, post-translational modifications of histones, and the human histone gene complement.

Web alert. Genetics of disease

A selection of World Wide Web sites relevant to papers published in this issue of Current Opinion in Genetics & Development.

GeneMachine: gene prediction and sequence annotation

MOTIVATION

A number of free-standing programs have been developed in order to help researchers find potential coding regions and deduce gene structure for long stretches of what is essentially 'anonymous DNA'. As these programs apply inherently different criteria to the question of what is and is not a coding region, multiple algorithms should be used in the course of positional cloning and positional candidate projects to assure that all potential coding regions within a previously-identified critical region are identified.

RESULTS

We have developed a gene identification tool called GeneMachine which allows users to query multiple exon and gene prediction programs in an automated fashion. BLAST searches are also performed in order to see whether a previously-characterized coding region corresponds to a region in the query sequence. A suite of Perl programs and modules are used to run MZEF, GENSCAN, GRAIL 2, FGENES, RepeatMasker, Sputnik, and BLAST. The results of these runs are then parsed and written into ASN.1 format. Output files can be opened using NCBI Sequin, in essence using Sequin as both a workbench and as a graphical viewer. The main feature of GeneMachine is that the process is fully automated; the user is only required to launch GeneMachine and then open the resulting file with Sequin. Annotations can then be made to these results prior to submission to GenBank, thereby increasing the intrinsic value of these data.

AVAILABILITY

GeneMachine is freely-available for download at http://genome.nhgri.nih.gov/genemachine. A public Web interface to the GeneMachine server for academic and not-for-profit users is available at http://genemachine.nhgri.nih.gov. The Web supplement to this paper may be found at http://genome.nhgri.nih.gov/genemachine/supplement/.

Cloning and characterization of 13 novel transcripts and the human RGS8 gene from the 1q25 region encompassing the hereditary prostate cancer (HPC1) locus

The aim of this study was to develop a saturated transcript map of the region encompassing the HPC1 locus to identify the susceptibility genes involved in hereditary prostate cancer (OMIM 176807) and hyperparathyroidism-jaw tumor syndrome (OMIM 145001). We previously reported the generation of a 6-Mb BAC/PAC contig of the candidate region and employed various strategies, such as database searching, exon-trapping, direct cDNA hybridization, and sample sequencing of BACs, to identify all potential transcripts. These efforts led to the identification and precise localization on the BAC contig of 59 transcripts representing 22 known genes and 37 potential transcripts represented by ESTs and exon traps. Here we report the detailed characterization of these ESTs into full-length transcript sequences, their expression pattern in various tissues, their genomic organization, and their homology to known genes. We have also identified an Alu insertion polymorphism in the intron of one of the transcripts. Overall, data on 13 novel transcripts and the human RGS8 gene (homologue of the rat RGS8 gene) are presented in this paper. Ten of the 13 novel transcripts are expressed in prostate tissue and represent positional candidates for HPC1.

Tissue-specific expression of a splicing mutation in the IKBKAP gene causes familial dysautonomia

Familial dysautonomia (FD; also known as "Riley-Day syndrome"), an Ashkenazi Jewish disorder, is the best known and most frequent of a group of congenital sensory neuropathies and is characterized by widespread sensory and variable autonomic dysfunction. Previously, we had mapped the FD gene, DYS, to a 0.5-cM region on chromosome 9q31 and had shown that the ethnic bias is due to a founder effect, with >99.5% of disease alleles sharing a common ancestral haplotype. To investigate the molecular basis of FD, we sequenced the minimal candidate region and cloned and characterized its five genes. One of these, IKBKAP, harbors two mutations that can cause FD. The major haplotype mutation is located in the donor splice site of intron 20. This mutation can result in skipping of exon 20 in the mRNA of patients with FD, although they continue to express varying levels of wild-type message in a tissue-specific manner. RNA isolated from lymphoblasts of patients is primarily wild-type, whereas only the deleted message is seen in RNA isolated from brain. The mutation associated with the minor haplotype in four patients is a missense (R696P) mutation in exon 19, which is predicted to disrupt a potential phosphorylation site. Our findings indicate that almost all cases of FD are caused by an unusual splice defect that displays tissue-specific expression; and they also provide the basis for rapid carrier screening in the Ashkenazi Jewish population.

Web alert. Oncogenes and cell proliferation

A selection of World Wide Web sites relevant to papers published in this issue of Current Opinion in Genetics & Development.

The human RGL (RalGDS-like) gene: cloning, expression analysis and genomic organization

Ral GDP dissociation stimulator (RalGDS) and its family members RGL, RLF and RGL2 are involved in Ras and Ral signaling pathways as downstream effector proteins. Here we report the precise localization and cloning of two forms of human RGL gene differing at the amino terminus. Transcript A, cloned from liver cDNA libraries has the same amino terminus as the mouse RGL, whereas transcript B cloned from brain has a substitution of 45 amino acids for the first nine amino acids. At the genomic level, exon 1 of transcript A is replaced by two alternative exons (1B1 and 1B2) in transcript B. Both forms share exons 2 through 18. The human RGL protein shares 94% amino acid identity with the mouse protein. Northern blot analysis shows that human RGL is expressed in a wide variety of tissues with strong expression being seen in the heart, brain, kidney, spleen and testis.

A 6-Mb high-resolution physical and transcription map encompassing the hereditary prostate cancer 1 (HPC1) region

Several hereditary disease loci have been genetically mapped to the chromosome 1q24-q31 interval, including the hereditary prostate cancer 1 (HPC1) locus. Here, we report the construction of a 20-Mb yeast artificial chromosome contig and a high-resolution 6-Mb sequence-ready bacterial artificial chromosome (BAC)/P1-derived artificial chromosome (PAC) contig of 1q25 by sequence and computational analysis, STS content mapping, and chromosome walking. One hundred thirty-six new STSs, including 10 novel simple sequence repeat polymorphisms that are being used for genetic refinement of multiple disease loci, have been generated from this contig and are shown to map to the 1q25 interval. The integrity of the 6-Mb BAC/PAC contig has been confirmed by restriction fingerprinting, and this contig is being used as a template for human chromosome 1 genome sequencing. A transcription mapping effort has resulted in the precise localization of 18 known genes and 31 ESTs by database searching, exon trapping, direct cDNA hybridization, and sample sequencing of BACs from the 1q25 contig. An additional 11 known genes and ESTs have been placed within the larger 1q24-q31 interval. These transcription units represent candidate genes for multiple hereditary diseases, including HPC1.

The histone database: a comprehensive WWW resource for histones and histone fold-containing proteins

The Histone Database (HDB) is an annotated and searchable collection of all full-length sequences and structures of histone and non-histone proteins containing the histone fold motif. These sequences are both eukaryotic and archaeal in origin. Several new histone fold-containing proteins have been identified, including Spt7p, and a few false positives have been removed from the earlier version of HDB. Database contents include compilations of post-translational modifications for each of the core and linker histones, as well as genomic information in the form of map loci for the human histone gene complement, with the genetic loci linked to Online Mendelian Inheritance in Man (OMIM). Conflicts between similar sequence entries from a number of source databases are also documented. Newly added to the HDB are multiple sequence alignments in which predicted functions of histone fold amino acid residues are annotated. The database is freely accessible through the WWW at http://genome.nhgri.nih.gov/histones/

Cloning, genomic organization and expression of a putative human transmembrane protein related to the Caenorhabditis elegans M01F1.4 gene

A novel human transcript CG-2 (C9ORF5), was isolated from the familial dysautonomia candidate region on 9q31 using a combination of cDNA selection and exon trapping. CG-2 was detected as a relatively abundant 8kb transcript in all adult and fetal tissues with the exception of adult thymus. Genomic analysis of CG-2 identified 18 exons that span more than 110kb. The gene encodes a 911-amino-acid protein with a predicted molecular weight of 101kDa and a hypothetical pI of 9.03. Sequence analysis of CG-2 indicates that it is likely to encode a transmembrane protein. Here, we assess CG-2 as a candidate for familial dysautonomia.

CACP, encoding a secreted proteoglycan, is mutated in camptodactyly-arthropathy-coxa vara-pericarditis syndrome

Altered growth and function of synoviocytes, the intimal cells which line joint cavities and tendon sheaths, occur in a number of skeletal diseases. Hyperplasia of synoviocytes is found in both rheumatoid arthritis and osteoarthritis, despite differences in the underlying aetiologies of the two disorders. We have studied the autosomal recessive disorder camptodactyly-arthropathy-coxa vara-pericarditis syndrome (CACP; MIM 208250) to identify biological pathways that lead to synoviocyte hyperplasia, the principal pathological feature of this syndrome. Using a positional-candidate approach, we identified mutations in a gene (CACP) encoding a secreted proteoglycan as the cause of CACP. The CACP protein, which has previously been identified as both 'megakaryocyte stimulating factor precursor' and 'superficial zone protein', contains domains that have homology to somatomedin B, heparin-binding proteins, mucins and haemopexins. In addition to expression in joint synovium and cartilage, CACP is expressed in non-skeletal tissues including liver and pericardium. The similarity of CACP sequence to that of other protein families and the expression of CACP in non-skeletal tissues suggest it may have diverse biological activities.