NEMBASE4: the nematode transcriptome resource

Nematode parasites are of major importance in human health and agriculture, and free-living species deliver essential ecosystem services. The genomics revolution has resulted in the production of many datasets of expressed sequence tags (ESTs) from a phylogenetically wide range of nematode species, but these are not easily compared. NEMBASE4 presents a single portal into extensively functionally annotated, EST-derived transcriptomes from over 60 species of nematodes, including plant and animal parasites and free-living taxa. Using the PartiGene suite of tools, we have assembled the publicly available ESTs for each species into a high-quality set of putative transcripts. These transcripts have been translated to produce a protein sequence resource and each is annotated with functional information derived from comparison with well-studied nematode species such as Caenorhabditis elegans and other non-nematode resources. By cross-comparing the sequences within NEMBASE4, we have also generated a protein family assignment for each translation. The data are presented in an openly accessible, interactive database. To demonstrate the utility of NEMBASE4, we have used the database to examine the uniqueness of the transcriptomes of major clades of parasitic nematodes, identifying lineage-restricted genes that may underpin particular parasitic phenotypes, possible viral pathogens of nematodes, and nematode-unique protein families that may be developed as drug targets.

DETECT--a density estimation tool for enzyme classification and its application to Plasmodium falciparum

MOTIVATION

A major challenge in genomics is the accurate annotation of component genes. Enzymes are typically predicted using homology-based search methods, where the membership of a protein to an enzyme family is based on single-sequence comparisons. As such, these methods are often error-prone and lack useful measures of reliability for the prediction.

RESULTS

Here, we present DETECT, a probabilistic method for enzyme prediction that accounts for the sequence diversity across enzyme families. By comparing the global alignment scores of an unknown protein to those of all known enzymes, an integrated likelihood score can be readily calculated, ranking the reaction classes relevant for that protein. Comparisons to BLAST reveal significant improvements in enzyme annotation accuracy. Applied to Plasmodium falciparum, we identify potential annotation errors and predict novel enzymes of therapeutic interest.

AVAILABILITY

A standalone application is available from the website: http://www.compsysbio.org/projects/DETECT/

Structural and functional characterization of SporoSAG: a SAG2-related surface antigen from Toxoplasma gondii

Toxoplasma gondii, the etiological agent of toxoplasmosis, utilizes stage-specific expression of antigenically distinct glycosylphosphatidylinositol-tethered surface coat proteins to promote and establish chronic infection. Of the three infective stages of T. gondii, sporozoites are encapsulated in highly infectious oocysts that have been linked to large scale outbreaks of toxoplasmosis. SporoSAG (surface antigen glycoprotein) is the dominant surface coat protein expressed on the surface of sporozoites. Using a bioinformatic approach, we show that SporoSAG clusters with the SAG2 subfamily of the SAG1-related superfamily (SRS) and is non-polymorphic among the 11 haplogroups of T. gondii strains. In contrast to the immunodominant SAG1 protein expressed on tachyzoites, SporoSAG is non-immunogenic during natural infection. We report the 1.60 A resolution crystal structure of SporoSAG solved using cadmium single anomalous dispersion. SporoSAG crystallized as a monomer and displays unique features of the SRS beta-sandwich fold relative to SAG1 and BSR4. Intriguingly, the structural diversity is localized to the upper sheets of the beta-sandwich fold and may have important implications for multimerization and host cell ligand recognition. The structure of SporoSAG also reveals an unexpectedly acidic surface that contrasts with the previously determined SAG1 and BSR4 structures where a basic surface is predicted to play a role in binding negatively charged glycosaminoglycans. Our structural and functional characterization of SporoSAG provides a rationale for the evolutionary divergence of this key SRS family member.

The origins of apicomplexan sequence innovation

The Apicomplexa are a group of phylogenetically related parasitic protists that include Plasmodium, Cryptosporidium, and Toxoplasma. Together they are a major global burden on human health and economics. To meet this challenge, several international consortia have generated vast amounts of sequence data for many of these parasites. Here, we exploit these data to perform a systematic analysis of protein family and domain incidence across the phylum. A total of 87,736 protein sequences were collected from 15 apicomplexan species. These were compared with three protein databases, including the partial genome database, PartiGeneDB, which increases the breadth of taxonomic coverage. From these searches we constructed taxonomic profiles that reveal the extent of apicomplexan sequence diversity. Sequences without a significant match outside the phylum were denoted as apicomplexan specialized. These were collated into 9134 discrete protein families and placed in the context of the apicomplexan phylogeny, identifying the putative origin of each family. Most apicomplexan families were associated with an individual genus or species. Interestingly, many genera-specific innovations were associated with specialized host cell invasion and/or parasite survival processes. Contrastingly, those families reflecting more ancestral relationships were enriched in generalized housekeeping functions such as translation and transcription, which have diverged within the apicomplexan lineage. Protein domain searches revealed 192 domains not previously reported in apicomplexans together with a number of novel domain combinations. We highlight domains that may be important to parasite survival.

On the extent and origins of genic novelty in the phylum Nematoda

BACKGROUND

The phylum Nematoda is biologically diverse, including parasites of plants and animals as well as free-living taxa. Underpinning this diversity will be commensurate diversity in expressed genes, including gene sets associated specifically with evolution of parasitism.

METHODS AND FINDINGS

Here we have analyzed the extensive expressed sequence tag data (available for 37 nematode species, most of which are parasites) and define over 120,000 distinct putative genes from which we have derived robust protein translations. Combined with the complete proteomes of Caenorhabditis elegans and Caenorhabditis briggsae, these proteins have been grouped into 65,000 protein families that in turn contain 40,000 distinct protein domains. We have mapped the occurrence of domains and families across the Nematoda and compared the nematode data to that available for other phyla. Gene loss is common, and in particular we identify nearly 5,000 genes that may have been lost from the lineage leading to the model nematode C. elegans. We find a preponderance of novelty, including 56,000 nematode-restricted protein families and 26,000 nematode-restricted domains. Mapping of the latest time-of-origin of these new families and domains across the nematode phylogeny revealed ongoing evolution of novelty. A number of genes from parasitic species had signatures of horizontal transfer from their host organisms, and parasitic species had a greater proportion of novel, secreted proteins than did free-living ones.

CONCLUSIONS

These classes of genes may underpin parasitic phenotypes, and thus may be targets for development of effective control measures.

Ubiquitination screen using protein microarrays for comprehensive identification of Rsp5 substrates in yeast

Ubiquitin-protein ligases (E3s) are responsible for target recognition and regulate stability, localization or function of their substrates. However, the substrates of most E3 enzymes remain unknown. Here, we describe the development of a novel proteomic in vitro ubiquitination screen using a protein microarray platform that can be utilized for the discovery of substrates for E3 ligases on a global scale. Using the yeast E3 Rsp5 as a test system to identify its substrates on a yeast protein microarray that covers most of the yeast (Saccharomyces cerevisiae) proteome, we identified numerous known and novel ubiquitinated substrates of this E3 ligase. Our enzymatic approach was complemented by a parallel protein microarray protein interaction study. Examination of the substrates identified in the analysis combined with phage display screening allowed exploration of binding mechanisms and substrate specificity of Rsp5. The development of a platform for global discovery of E3 substrates is invaluable for understanding the cellular pathways in which they participate, and could be utilized for the identification of drug targets.

PartiGene--constructing partial genomes

Expressed sequence tags (ESTs) offer a low-cost approach to gene discovery and are being used by an increasing number of laboratories to obtain sequence information for a wide variety of organisms. The challenge lies in processing and organizing this data within a genomic context to facilitate large scale analyses. Here we present PartiGene, an integrated sequence analysis suite that uses freely available public domain software to (1) process raw trace chromatograms into sequence objects suitable for submission to dbEST; (2) place these sequences within a genomic context; (3) perform customizable first-pass annotation of the data; and (4) present the data as HTML tables and an SQL database resource. PartiGene has been used to create a number of non-model organism database resources including NEMBASE (http://www.nematodes.org) and LumbriBase (http://www.earthworms.org/). The packages are readily portable, freely available and can be run on simple Linux-based workstations.

AVAILABILITY

PartiGene is available from http://www.nematodes.org/PartiGene and also forms part of the EST analysis software, associated with the Natural Environmental Research Council (UK) Bio-Linux project (http://envgen.nox.ac.uk/biolinux.html).

Identification and characterization of two novel tetratricopeptide repeat-containing genes

The tetratricopeptide repeat (TPR) is a degenerate, repeating amino acid motif of 34 residues that has been identified in a variety of proteins; however, no biochemical function has been established for it. In a two-hybrid screen for interaction with the GAP-related domain of neurofibromin, the product of the NF1 gene, we have identified two novel human genes encoding proteins with TPR motifs. The first, represented by cDNA tpr1, is located in chromosome 5q32-33.2. It is ubiquitously expressed as a 1.6-kb transcript that encodes three tandem TPR motifs, but is not related to any other known gene outside this domain. The second, defined by cDNA tpr2, maps to human chromosome 17q11.2-23. It is ubiquitously expressed as a 2.2-kb transcript encoding seven TPR units. Interestingly, a separate region of the tpr2 cDNA has striking similarity to the "J region" of the DnaJ family. The products of the tpr1 and tpr2 cDNAs interact preferentially with a truncated form of the NF1 GAP-related domain via their TPR regions, suggesting that they may be targeted to an abnormality of protein folding.

Clinical evidence of genetic anticipation in adult-onset idiopathic dystonia

Idiopathic dystonia occurs in both hereditary and sporadic forms. In this report, we studied the age of onset and family history of 260 patients (probands) with idiopathic adult-onset dystonia (IAD), cranial or cervical. The mean age at onset of these patients was (45.71 +/- 15.85) years. Forty-nine probands had a positive family history of dystonia or tremor in first- and second-degree relatives, and 7 had affected siblings only. The significance of tremor as a part of clinical manifestation of dystonia was evidenced by a high frequency of postural or action tremor in patients and relatives. Retrospectively, we examined the age of onset of dystonia (cervical or cranial) on successive generations in 49 families. Age of onset of clinical symptoms was earlier, by an average of 21.25 years, in the second generation than in the first generation. The mean age at onset of affected family members differed significantly between successive generations in these 49 families (p = 1.11 x 10(-8)). Our results suggest a tendency for earlier onset of dystonia and worsening of disease phenotype in succeeding generations in the same family. These findings are most compatible with genetic anticipation and suggest that an unstable trinucleotide repeat is most likely involved in adult-onset primary cranial or cervical dystonia. In addition, tremor as an integral part of dystonia needs further evaluation by molecular genetic studies.

A contig of non-chimaeric YACs containing the spinal muscular atrophy gene in 5q13

We have constructed a contig of non-chimaeric yeast artificial chromosomes (YACs) across the candidate region for childhood autosomal recessive spinal muscular atrophy (SMA) in 5q13. A novel microsatellite reduces the candidate region to approximately 400kb of DNA distal to D5S435. The candidate region contains blocks of chromosome 5 specific repeats which have copies on 5p as well as elsewhere on 5q. Restriction mapping of the YACs reveals at least one CpG island in the SMA gene region. The YAC maps indicate that the contig contains minimal rearrangements or deletions. The data show the value of screening several YAC libraries simultaneously in order to construct a set of overlapping sequences suitable for candidate gene searches and direct genomic sequencing.

Assay by polymerase chain reaction (PCR) of multi-allele polymorphisms in the Huntington's disease region of chromosome 4

The Huntington's disease-linked D4S115 marker has been converted from a DNA blot assay to a more sensitive and rapid polymerase chain reaction (PCR) assay. PCR amplification of a tandem repeat at D4S115 revealed 7 allelic fragments, ranging in size from approximately 610 to 915 bp, differing in their apparent copy number of a approximately 55 bp core repeat. This repeat unit differs strikingly in sequence from the repeat units of other multi-allele markers from chromosome region 4p 16.3, arguing that the VNTR (Variable Number of Tandem Repeats) loci clustered in this region did not arise from a common ancestral sequence. The D4S115 marker can be assayed simultaneously with PCR products from D4S125, D4S95 and D4S43 on a single agarose gel, providing a rapid scan for successful amplification of these difficult-to-assay VNTRs, and for inheritance of the entire candidate Huntington's disease region. This approach should help to increase the speed, informativeness and accuracy of presymptomatic and prenatal linkage testing in this devastating disorder.

FLT4, a novel class III receptor tyrosine kinase in chromosome 5q33-qter

The receptors for at least two hematopoietic growth factors, namely the stem cell factor and colony-stimulating factor 1, belong to class III receptor tyrosine kinases. Here we describe cloning of a partial complementary DNA for FLT4, an additional member of this gene family from human leukemia cells. The FLT4 tyrosine kinase domain is 79% homologous with the previously cloned FLT1 (M. Shibuya et al., Oncogene, 5: 519-524, 1990) tyrosine kinase and maps to the chromosomal region 5q33-qter. We have found FLT4 expression in human placenta, lung, heart, and kidney, whereas the pancreas and brain appeared to contain very little if any FLT4 RNA. The results suggest that FLT4 functions in multiple adult tissues.

Complex patterns of linkage disequilibrium in the Huntington disease region

The genetic defect causing Huntington disease (HD) has been mapped to 4p16.3 by linkage analysis using DNA markers. Two apparently contradictory classes of recombination events in HD kindreds preclude precise targeting of efforts to clone the disease gene. Here, we report a new recombination event that increases support for an internal candidate region of 2.5 Mb between D4S10 and D4S168. Analysis of 23 DNA polymorphisms in 4p16.3 revealed a complex pattern of association with the disease gene that failed to narrow the size of the candidate region. The degree of linkage disequilibrium did not show a continuous increase across the physical map, nor was a region of extreme disequilibrium identified. Markers displaying no association with the disorder were interspersed with and, in many cases, close to markers displaying significant disequilibrium. Comparison of closely spaced marker pairs on normal and HD chromosomes, as well as analysis of haplotypes across the HD region, suggest that simple recombination subsequent to a single original HD mutation cannot easily explain the pool of HD chromosomes seen today. A number of different mechanisms could contribute to the diversity of haplotypes observed on HD chromosomes, but it is likely that there has been more than one and possibly several independent origins of the HD mutation.

A genetic linkage map of human chromosome 5 with 60 RFLP loci

A genetic map of human chromosome 5 that contains 60 restriction fragment length polymorphism (RFLP) loci in one linkage group has been constructed. Segregation data using these markers and 40 large multigenerational families supplied by the Centre d'Etude du Polymorphisme Humain have been collected. Linkage analyses were performed with the program package CRI-MAP; using odds greater than 1000:1, 30 RFLP loci could be placed on the map. This genetic map spans 289 cM sex-equal, 353 cM in females, and 244 cM in males. While the relative rate of recombination for female meioses is nearly twice that of males over much of the chromosome, several instances of statistically significant excess male recombination were observed. The order of probes on the genetic map has been confirmed by their physical order as determined by somatic cell hybrid lines containing deletions of normal chromosome 5. There is concordance between the physical positions of markers and their genetic positions. Our most distal probes on the genetic map are cytologically localized to the most distal portions of the chromosome. This suggests that our genetic map spans most of chromosome 5.

Evidence from family studies that the gene causing Huntington disease is telomeric to D4S95 and D4S90

A DNA probe (D4S95) that detects a variable number of tandem repeats and a single-site-variation polymorphism after digestion with a single restriction enzyme, AccI, has previously been described. The order of this probe relative to the gene for Huntington disease (HD) and other previously described markers has not been established. Analysis of 24 affected families with HD has shown that D4S95 is in tight linkage with the gene causing HD, with a maximal Lod score of 12.489 at a theta of .03. D4S90 is a probe which maps to 4p16.3, telomeric to D4S95, and detects polymorphisms with HincII and other enzymes. In one affected person, recombination has occurred between D4S10 and HD, between D4S95 and HD, and in all likelihood also between D4S90 and HD, which strongly suggests that the gene for HD is telomeric to all these DNA probes. This suggests that the gene causing HD is located in the most distal region of the short arm of chromosome 4, flanked by D4S90 and the telomere, and supports the locus order D4S10-D4S95-D4S90-HD-telomere. D4S95 is a most useful DNA marker for predictive testing programs, while D4S90 will serve as a useful starting point for identifying DNA fragments closer to the gene for HD.

Localization of a susceptibility locus for schizophrenia on chromosome 5

Schizophrenia is a common disorder with a life time prevalence of approximately 1 per cent. The illness often develops in young adults, who were previously normal, and is characterized by a constellation of symptoms including hallucinations and delusions (psychotic symptoms) and symptoms such as severely inappropriate emotional responses, a disorder of thinking and concentration, erratic behaviour as well as social and occupational deterioration. A considerable proportion of the variance in the liability to develop schizophrenia may be genetic, but segregation analysis, to establish a mode of transmission, has not produced a consistent result. One of these studies was carried out in Iceland and made use of the large family size and extensive geneaological information present in that country. Here we demonstrate genetic linkage of two DNA polymorphisms on the long arm of human chromosome 5 to schizophrenia in seven British and Icelandic families with multiple affected members. The results indicate the existence of a gene locus with a dominant schizophrenia-susceptibility allele. Inheritance of the allele in the families studied suggests that it may also predispose to psychiatric conditions such as schizophrenia spectrum disorders and a variety of other disorders. This report provides the first strong evidence for the involvement of a single gene in the causation of schizophrenia.

Localization of the genetic defect in familial adenomatous polyposis within a small region of chromosome 5

Familial adenomatous polyposis (FAP), a Mendelian disorder that includes familial polyposis coli (FPC) and Gardner syndrome (GS), has an autosomal dominant mode of inheritance. It is characterized by hundreds to thousands of adenomatous polyps that can progress to carcinoma of the colon, suggesting that the gene that harbors the FAP germ-line mutation may play an important role in the somatic genetic pathway to colon cancer. The defect responsible for FAP was recently mapped to the long arm of chromosome 5 by linkage between the FPC phenotype and a locus defined by DNA probe pC11p11 (D5S71), located at 5q21-22. Because an important next step in the paradigm for identification of a disease gene is to obtain a more precise localization, we isolated and mapped by linkage six additional polymorphic DNA markers in the FAP region. Subsequent linkage analysis in six pedigrees, three having the FPC phenotype and three segregating GS, placed the FAP locus very close to a new marker, YN5.48 (D5S81), that is approximately 17 centimorgans distal to C11p11 on the genetic map. The analysis revealed no evidence of genetic heterogeneity between the two phenotypes, a question that had not been clearly resolved by the earlier studies. The new set of markers in the near vicinity of the FAP locus represents a further step toward isolation of the genetic defect and provides the opportunity for preclinical diagnosis of risk status for colon cancer among individuals in families that are segregating adenomatous polyposis.

A DNA probe, D5 [D4S90] mapping to human chromosome 4p16.3

Images

Ribosomal protein gene sequences map to human chromosomes 5, 8, and 17

DNA sequences complementary to six mammalian ribosomal protein (r-protein) cDNAs are assigned to human chromosomal linkage groups in human-Chinese hamster hybrid cell clones. Ten r-protein DNA fragments map to chromosomes 5, 8 and 17, indicating that these important, housekeeping genes are distributed to multiple sites in the human genome. Each of the chromosome assignments, determined initially by surveying Chinese hamster-human hybrid cell clones with complex karyotypes using Chinese hamster and human cDNA probes, were confirmed in critical minipanels of highly reduced or monochromosomal hybrid cells. As all 10 fragments mapped to only three human chromosomes, r-protein sequences appear to be distributed nonrandomly within human DNA. The r-protein S14 sequence assigned to human chromosome 5 (5q23-5q33) rescues Chinese hamster emetine-resistance mutations (emt b) in interspecific hybrids. Therefore, this sequence corresponds to the transcriptionally active human RPS14 gene. In contrast, other r-protein DNA sequences examined likely are a mixture of functional genes and inactive pseudogenes.

Assignment of human 3-hydroxy-3-methylglutaryl coenzyme A reductase gene to q13----q23 region of chromosome 5

We have used hamster cDNA probes for 3-hydroxy-3-methylglutaryl coenzyme A (HMG CoA) reductase (HMGCR) to determine the chromosomal location of the human gene for HMG CoA reductase. Southern blot analysis of genomic DNA from 16 independent mouse-human somatic cell hybrids showed that the human gene for HMG CoA reductase resides on chromosome 5. Analysis of Chinese hamster-human somatic cell hybrids selectively retaining human 5 or a portion of it showed that the gene locus for HMG CoA reductase can be assigned to the q13----q23 region of chromosome 5.

Assignment of human dihydrofolate reductase gene to band q23 of chromosome 5 and of related pseudogene psi HD1 to chromosome 3

The chromosomal location of the human dihydrofolate reductase (DHFR; EC 1.5.1.3) gene that is amplified in a methotrexate-resistant human cell line has been investigated by screening a number of human-Chinese hamster ovary cell hybrids containing terminal and interstitial deletions in human chromosome 5. A correlation of genomic blotting data with the chromosome 5 constitution of the individual hybrids has allowed the assignment of the human DHFR gene to 5q23. The present work also establishes the location of the related intronless pseudogene psi HD1 in chromosome 3.