GeneSeeker: extraction and integration of human disease-related information from web-based genetic databases.

The identification of genes underlying human genetic disorders requires the combination of data related to cytogenetic localization, phenotypes and expression patterns, to generate a list of candidate genes. In the field of human genetics, it is normal to perform this combination analysis by hand. We report on GeneSeeker (http://www.cmbi.ru.nl/GeneSeeker/), a web server that gathers and combines data from a series of databases. All database searches are performed via the web interfaces provided with the original databases, guaranteeing that the most recent data are queried, and obviating data warehousing. GeneSeeker makes the same selection of candidate genes as the human geneticists would have performed, and thus reducing the time-consuming process to a few minutes. GeneSeeker is particularly well suited for syndromes in which the disease gene displays altered expression patterns in the affected tissue(s).

USAGE: a web-based approach towards the analysis of SAGE data. Serial Analysis of Gene Expression.

MOTIVATION: SAGE enables the determination of genome-wide mRNA expression profiles. A comprehensive analysis of SAGE data requires software, which integrates (statistical) data analysis methods with a database system. Furthermore, to facilitate data sharing between users, the application should reside on a central server and be accessed via the internet. Since such an application was not available we developed the USAGE package. RESULTS: USAGE is a web-based application that comprises an integrated set of tools, which offers many functions for analysing and comparing SAGE data. Additionally, USAGE includes a statistical method for the planning of new SAGE experiments. USAGE is available in a multi-user environment giving users the option of sharing data. USAGE is interfaced to a relational database to store data and analysis results. The USAGE query editor allows the composition of queries for searching this database. Several database functions have been included which enable the selection and combination of data. USAGE provides the biologist increased functionality and flexibility for analysing SAGE data. AVAILABILITY: USAGE is freely accessible for academic institutions at http://www.cmbi.kun.nl/usage/. The source code of USAGE is freely available for academic institutions on request from the first author.

Multiple acquisition of methanogenic archaeal symbionts by anaerobic ciliates.

Anaerobic heterotrichous ciliates (Armophoridae and Clevelandellidae) possess hydrogenosomes that generate molecular hydrogen and ATP. This intracellular source of hydrogen provides the basis for a stable endosymbiotic association with methanogenic archaea. We analyzed the SSU rRNA genes of 18 heterotrichous anaerobic ciliates and their methanogenic endosymbionts in order to unravel the evolution of this mutualistic association. Here, we show that the anaerobic heterotrichous ciliates constitute at least three evolutionary lines. One group consists predominantly of gut-dwelling ciliates, and two to three, potentially four, additional clades comprise ciliates that thrive in freshwater sediments. Their methanogenic endosymbionts belong to only two different taxa that are closely related to free-living methanogenic archaea from the particular ecological niches. The close phylogenetic relationships between the endosymbionts and free-living methanogenic archaea argue for multiple acquisitions from environmental sources, notwithstanding the strictly vertical transmission of the endosymbionts. Since phylogenetic analysis of the small-subunit (SSU) rRNA genes of the hydrogenosomes of these ciliates indicates a descent from the mitochondria of aerobic ciliates, it is likely that anaerobic heterotrichous ciliates hosted endosymbiotic methanogens prior to their radiation. Therefore, our data strongly suggest multiple acquisitions and replacements of endosymbiotic methanogenic archaea during their host's adaptation to the various ecological niches.

Oreochromis mossambicus (tilapia) corticotropin-releasing hormone: cDNA sequence and bioactivity.

Although hypothalamic corticotropin-releasing hormone (CRH) is involved in the stress response in all vertebrate groups, only a limited number of studies on this neuroendocrine peptide deals with non-mammalian neuroendocrine systems. We determined the cDNA sequence of the CRH precursor of the teleost Oreochromis mossambicus (tilapia) and studied the biological potency of the CRH peptide in a homologous teleost bioassay. Polymerase chain reaction (PCR) with degenerate and specific primers yielded fragments of tilapia CRH cDNA. Full-length CRH cDNA (988 nucleotides) was obtained by screening a tilapia hypothalamus cDNA library with the tilapia CRH PCR products. The precursor sequence (167 amino acids) contains a signal peptide, the CRH peptide and a motif conserved among all vertebrate CRH precursors. Tilapia CRH (41 aa) displays between 63% and 80% amino acid sequence identity to CRH from other vertebrates, whereas the degree of identity to members of the urotensin I/urocortin lineage is considerably lower. In a phylogenetic tree, based on alignment of all full CRH peptide precursors presently known, the three teleost CRH precursors (tilapia; sockeye salmon, Oncorhynchus nerka; white sucker, Catostomus commersoni) form a monophyletic group distinct from amphibian and mammalian precursors. Despite the differences between the primary structures of tilapia and rat CRH, maximally effective concentrations of tilapia and rat CRH were equally potent in stimulating adrenocorticotropic hormone (ACTH) and alpha-MSH release by tilapia pituitaries in vitro. The tilapia and salmon CRH sequences show that more variation exists between orthologous vertebrate CRH structures, and teleost CRHs in particular than previously recognized. Whether the structural differences reflect different mechanisms of action of this peptide in the stress response remains to be investigated.

Cloning and expression of two proopiomelanocortin mRNAs in the common carp (Cyprinus carpio L.).

Proopiomelanocortin (POMC) is the precursor for a number of biologically active peptides such as adrenocorticotropic hormone (ACTH), alpha-melanocyte-stimulating hormone (alpha-MSH) and beta-endorphin. It is well known that these peptides are involved in the stress response in fish as well as in mammals. We have cloned two different carp POMC cDNAs called, POMC-I and POMC-II. The nucleotide sequences of 955 bp for POMC-I and 959 bp for POMC-II share 93.5% identity in their cDNAs, and the deduced amino acid sequences (both 222 amino acids) are 91.4% identical. In the ACTH and beta-MSH domain, two amino acid substitutions are found, whereas alpha-MSH and beta-endorphin are identical. For beta-MSH, the serine replacement (in POMC-I) by a glycine (in POMC-II) results in a putative amidation site Pro-X-Gly for POMC-II. We used RT-PCR to show that both POMC mRNAs are expressed in the hypophysis, hypothalamus and other parts of the brain of a single fish. Furthermore, in a phylogenetic tree based on POMC sequences the divergence of carp POMC-I and -II from tetraploid animals (salmon, trout and xenopus) is demonstrated.

Definition of natural T cell antigens with mimicry epitopes obtained from dedicated synthetic peptide libraries.

Progress has recently been made in the use of synthetic peptide libraries for the identification of T cell-stimulating ligands. T cell epitopes identified from synthetic libraries are mimics of natural epitopes. Here we show how the mimicry epitopes obtained from synthetic peptide libraries enable unambiguous identification of natural T cell Ags. Synthetic peptide libraries were screened with Mycobacterium tuberculosis-reactive and -autoreactive T cell clones. In two cases, database homology searches with mimicry epitopes isolated from a dedicated synthetic peptide library allowed immediate identification of the natural antigenic protein. In two other cases, an amino acid pattern that reflected the epitope requirements of the T cell was determined by substitution and omission mixture analysis. Subsequently, the natural Ag was identified from databases using this refined pattern. This approach opens new perspectives for rapid and reliable Ag definition, representing a feasible alternative to the biochemical and genetic approaches described thus far.

Distinct frequency-distributions of homopolymeric DNA tracts in different genomes.

The unusual base composition of the genome of the human malaria parasite Plasmodium falciparum prompted us to systematically investigate the occurrence of homopolymeric DNA tracts in the P. falciparum genome and, for comparison, in the genomes of Homo sapiens , Saccharomyces cerevisiae , Caenorhabditis elegans , Arabidopsis thaliana , Escherichia coli and Mycobacterium tuberculosis. Comparison of theobserved frequencies with the frequencies as expected for random DNA revealed that homopolymeric (dA:dT) tracts occur well above chance in the eukaryotic genome. In the majority of these genomes, (dA:dT) tract overrepresentation proved to be an exponential function of the tract length. (dG:dC) tract overrepresentation was absent or less pronounced in both prokaryotic and eukaryotic genomes. On the basis of our results, we propose that homopolymeric (dA:dT) tracts are expanded via replication slippage. This slippage-mediated expansion does not operate on tracts with lengths below a critical threshold of 7-10 bp.

Genealogy of the alpha-crystallin--small heat-shock protein superfamily.

Sequences of 40 very diverse representatives of the alpha-crystallin-small heat-shock protein (alpha-Hsp) superfamily are compared. Their characteristic C-terminal 'alpha-crystallin domain' of 80-100 residues contains short consensus sequences that are highly conserved from prokaryotes to eukaryotes. There are, in addition, some positions that clearly distinguish animal from non-animal alpha-Hsps. The alpha-crystallin domain is predicted to consist of two hydrophobic beta-sheet motifs, separated by a hydrophilic region which is variable in length. Combination of a conserved alpha-crystallin domain with a variable N-terminal domain and C-terminal extension probably modulates the properties of the various alpha-Hsps as stress-protective and structural oligomeric proteins. Phylogeny reconstruction indicates that multiple alpha-Hsps were already present in the last common ancestor of pro- and eukaryotes. It is suggested that during eukaryote evolution, animal and non-animal alpha-Hsps originated from different ancestral gene copies. Repeated gene duplications gave rise to the multiple alpha-Hsps present in most organisms.

Subtilases: the superfamily of subtilisin-like serine proteases.

Subtilases are members of the clan (or superfamily) of subtilisin-like serine proteases. Over 200 subtilases are presently known, more than 170 of which with their complete amino acid sequence. In this update of our previous overview (Siezen RJ, de Vos WM, Leunissen JAM, Dijkstra BW, 1991, Protein Eng 4:719-731), details of more than 100 new subtilases discovered in the past five years are summarized, and amino acid sequences of their catalytic domains are compared in a multiple sequence alignment. Based on sequence homology, a subdivision into six families is proposed. Highly conserved residues of the catalytic domain are identified, as are large or unusual deletions and insertions. Predictions have been updated for Ca(2+)-binding sites, disulfide bonds, and substrate specificity, based on both sequence alignment and three-dimensional homology modeling.

Molecular characterization of the human peroxisomal branched-chain acyl-CoA oxidase: cDNA cloning, chromosomal assignment, tissue distribution, and evidence for the absence of the protein in Zellweger syndrome.

Peroxisomes in human liver contain two distinct acyl-CoA oxidases with different substrate specificities: (i) palmitoyl-CoA oxidase, oxidizing very long straight-chain fatty acids and eicosanoids, and (ii) a branched-chain acyl-CoA oxidase (hBRCACox), involved in the degradation of long branched fatty acids and bile acid intermediates. The accumulation of branched fatty acids and bile acid intermediates leads to severe mental retardation and death of the diseased children. In this study, we report the molecular characterization of the hBRCACox, a prerequisite for studying mutations in patients with a single enzyme deficiency. The composite cDNA sequence of hBRCACox, derived from overlapping clones isolated via immunoscreening and hybridization of human liver cDNA expression libraries, consisted of 2225 bases and contained an open reading frame of 2046 bases, encoding a protein of 681 amino acids with a calculated molecular mass of 76,739 Da. The C-terminal tripeptide of the protein is SKL, a known peroxisome targeting signal. Sequence comparison with the other acyl-CoA oxidases and evolutionary analysis revealed that, despite its broader substrate specificity, the hBRCACox is the human homolog of rat trihydroxycoprostanoyl-CoA oxidase (rTHCCox) and that separate gene duplication events led to the occurrence in mammals of acyl-CoA oxidases with different substrate specificities. Northern blot analysis demonstrated that--in contrast to the rTHCCox gene--the hBRCACox gene is transcribed also in extrahepatic tissues such as heart, kidney, skeletal muscle, and pancreas. The highest levels of the 2.6-kb mRNA were found in heart, followed by liver. The enzyme is encoded by a single-copy gene, which was assigned to chromosome 3p14.3 by fluorescent in situ hybridization. It was absent from livers of Zellweger patients as shown by immunoblot analysis and immunocytochemistry.

Localization of a gene for Möbius syndrome to chromosome 3q by linkage analysis in a Dutch family.

Möbius syndrome (MIM no. 157900) consists of a congenital paresis or paralysis of the VIIth cranial nerve, frequently accompanied by paralysis of other cranial nerves, orofacial and limb malformations, defects of the musculoskeletal system and mental retardation. Although most patients are sporadic cases, familial recurrence is not rare. Different pedigrees suggest different modes of inheritance. We performed linkage analysis in a large family with autosomal dominantly inherited Möbius syndrome, consisting essentially of asymmetric bilateral facial pareses. After exclusion of the candidate region for Möbius syndrome on 13q12.2-q13, we localized the gene to chromosome 3q21-22, indicating genetic heterogeneity of Möbius syndrome. This heterogeneity is further proven by the exclusion of both loci in a second family with Möbius syndrome.

Cloning and sequence analysis of a hypothalamic cDNA encoding a D1c dopamine receptor in tilapia.

Physiological and pharmacological studies have indicated that during acid stress a D1-like dopamine receptor becomes functional on intermediate pituitary melanocyte-stimulating hormone cells of tilapia (Oreochromis mossambicus). As a first step towards physiological expression studies we isolated a D1-like dopamine receptor from a tilapia hypothalamus cDNA library. Construction of a phylogenetic tree of most of the D1-like receptors known in human, rat, Xenopus, goldfish and Drosophila revealed that the here presented clone is most likely the tilapia equivalent of the Xenopus D1c dopamine receptor.

Protein sequences indicate that turtles branched off from the amniote tree after mammals.

The phylogenetic relationships among the major groups of amniote vertebrates remain a matter of controversy. Various alternatives for the position of the turtles have been proposed, branching off either before or after the mammals. To discover the phylogenetic position of turtles in relation to mammals and birds, we have determined cDNA sequences for the eye lens proteins alpha A- and alpha B-crystallin of the red-eared slider turtle (Trachemys scripta elegans). In addition, databases were searched for turtle protein sequences, for which mammalian, avian, and outgroup orthologs were available. All sequences were analyzed by three phylogenetic tree reconstruction methods (neighbor-joining, maximum parsimony, and maximum likelihood). Including the alpha-crystallins, 7 out of 12 proteins support a sister-group relation of turtles and birds with all 3 methods. For each of the other five proteins no topology was consistently preferred by the three approaches. Analyses of the combined amino acid data (1,695 aligned sites) also give extremely strong evidence that turtles are nearer to birds, indicating that mammals branched off before the divergence between turtles and birds occurred.

Evidence for limited within-person evolution of the V3 domain of the HIV-1 envelope in the amsterdam population.

OBJECTIVE: To study the development of the V3 region of the HIV-1 envelope over time, both within subjects and population-wide. METHODS: Direct V3 sequences were obtained from viral RNA from seroconversion samples of 138 individuals [32 intravenous drug users (IVDU), 106 homosexual men], as well as from 5-year follow-up samples of 45 of these individuals (11 IVDU, 34 homosexual men). RESULTS: The population-wide variation of the V3 region in both the seroconversion samples and the 5-year samples steadily increased over consecutive years and were of similar magnitude in each calendar year. The variation in the IVDU group was slightly lower (presumably lagging behind) than in the homosexual group, but also increased over time. The consensus sequence, representing the centre of the swarm of variants, remained almost stationary in 10 years of evolution. The V3 sequences from virions in serum collected 5 years after seroconversion still resembled those from the seroconversion sample, either in overall similarity or in specific (signature) amino acids. Seroconversion and late sequences from a donor-recipient pair were also very similar. CONCLUSIONS: The variation in V3 sequences from seroconversion samples is as large as that in 5-year follow-up samples from the same calendar year, suggesting that there is no strong selection for a particular V3 genotype at transmission. The HIV-1 subtype B quasispecies in a naive population appears to evolve through unbiased expansion around a stationary consensus sequence. Despite its large variability, the V3 region retains many of its individual characteristics after 5 years of infection. This indicates that the sampling moment (relative to the seroconversion data) will not greatly influence the results of phylogenetic analyses.

The expanding small heat-shock protein family, and structure predictions of the conserved "alpha-crystallin domain".

The ever-increasing number of proteins identified as belonging to the family of small heat-shock proteins (shsps) and alpha-crystallins enables us to reassess the phylogeny of this ubiquitous protein family. While the prokaryotic and fungal representatives are not properly resolved, most of the plant and animal shsps and related proteins are clearly grouped in distinct clades, reflecting a history of repeated gene duplications. The members of the shsp family are characterized by the presence of a conserved homologous "alpha-crystallin domain," which sometimes is present in duplicate. Predictions are made of secondary structure and solvent accessibility of this domain, which together with hydropathy profiles and intron positions support the presence of two similar hydrophobic beta-sheet-rich motifs, connected by a hydrophilic alpha-helical region. Together with an overview of the newly characterized members of the shsp family, these data help to define this family as being involved as stable structural proteins and as molecular chaperones during normal development and induced under pathological and stressful conditions.

Parasitic apicomplexans harbor a chlorophyll a-D1 complex, the potential target for therapeutic triazines.

Ultrastructural evidence is presented for the presence of plastid-like organelles in Toxoplasma gondii, Sarcocystis muris, Babesia ovis, and Plasmodium falciparum. In addition, it was shown that merozoites of T. gondii contain protochlorophyllidae a and traces of chlorophyll a bound to the photosynthetic reaction centers I PS I and PS II. A psbA gene was isolated from merozoites of S. muris by the polymerase chain reaction (PCR). Partial sequencing of the PCR product revealed that the herbicide-binding region is highly conserved. Therefore, it is likely that the sensitivity of apicomplexans to the herbicide toltrazuril depends on the interaction of the herbicide with the D1 protein of the photosynthetic reaction center of the parasite's organelles.

Molecular evolution of interleukin-3.

Chimpanzee, tamarin, and marmoset interleukin-3 (IL-3) genes were cloned, sequenced, and expressed. Western blot analysis demonstrated that functional genes were isolated. IL-3 sequences were compared with those of mouse, rat, rhesus monkey, gibbon, and man. Multiple alignment of the IL-3 coding regions showed that only a few regions had been conserved during mammalian evolution, which are likely associated with functional domains of the IL-3 protein. Substitution rates for the various lineages were calculated and the numbers of synonymous and nonsynonymous substitutions were estimated separately. Distance matrices of the IL-3 coding regions were used to construct phylogenetic trees which revealed large differences in IL-3 evolution rate as well as a more rapid substitution rate for rodents and a rate slowdown during hominoid evolution. Extremes were rhesus monkey IL-3, which accumulated few synonymous substitutions, and gibbon IL-3, which had almost exclusively synonymous substitutions. In rhesus monkey IL-3, nonsynonymous substitutions outnumbered synonymous substitutions, which could not be readily explained by a random process of substitutions. We assume that during evolution of IL-3, the majority of the amino acid replacements and the impaired interspecies functional cross-reactivity originate from selection mechanisms with the most likely selective force being the structure of the heterodimeric IL.3 cell-surface receptor. Insight into IL-3 architecture and structural analysis of the IL-3 receptor are needed to analyze the unusually fast evolution of IL-3 in more detail.

Divergence in codon usage of Lactobacillus species.

We have analyzed codon usage patterns of 70 sequenced genes from different Lactobacillus species. Codon usage in lactobacilli is highly biased. Both inter-species and intra-species heterogeneity of codon usage bias was observed. Codon usage in L. acidophilus is similar to that in L. helveticus, but dissimilar to that in L. bulgaricus, L. casei, L. pentosus and L. plantarum. Codon usage in the latter three organisms is not significantly different, but is different from that in L. bulgaricus. Inter-species differences in codon usage can, at least in part, be explained by differences in mutational drift. L. bulgaricus shows GC drift, whereas all other species show AT drift. L. acidophilus and L. helveticus rarely use NNG in family-box (a set of synonymous) codons, in contrast to all other species. This result may be explained by assuming that L. acidophilus and L. helveticus, but not other species examined, use a single tRNA species for translation of family-box codons. Differences in expression level of genes are positively correlated with codon usage bias. Highly expressed genes show highly biased codon usage, whereas weakly expressed genes show much less biased codon usage. Codon usage patterns at the 5'-end of Lactobacillus genes is not significantly different from that of entire genes. The GC content of codons 2-6 is significantly reduced compared with that of the remainder of the gene. The possible implications of a reduced GC content for the control of translation efficiency are discussed.

A molecular dynamics approach for the generation of complete protein structures from limited coordinate data.

Generation of full protein coordinates from limited information, e.g., the C alpha coordinates, is an important step in protein homology modeling and structure determination, and molecular dynamics (MD) simulations may prove to be important in this task. We describe a new method, in which the protein backbone is built quickly in a rather crude way and then refined by minimization techniques. Subsequently, the side chains are positioned using extensive MD calculations. The method is tested on two proteins, and results compared to proteins constructed using two other MD-based methods. In the first method, we supplemented an existing backbone building method with a new procedure to add side chains. The second one largely consists of available methodology. The constructed proteins are compared to the corresponding X-ray structures, which became available during this study, and they are in good agreement (backbone RMS values of 0.5-0.7 A, and all-atom RMS values of 1.5-1.9 A). This comparative study indicates that extensive MD simulations are able, to some extent, to generate details of the native protein structure, and may contribute to the development of a standardized methodology to predict reliably (parts of) protein structures when only partial coordinate data are available.

Evolution of the alpha-crystallin/small heat-shock protein family.

The common characteristic of the alpha-crystallin/small heat-shock protein family is the presence of a conserved homologous sequence of 90-100 residues. Apart from the vertebrate lens proteins--alpha A- and alpha B-crystallin--and the ubiquitous group of 15-30-kDa heat-shock proteins, this family also includes two mycobacterial surface antigens and a major egg antigen of Schistosoma mansoni. Multiple small heat-shock proteins are especially present in higher plants, where they can be distinguished in at least two classes of cytoplasmic proteins and a chloroplast-located class. The alpha-crystallins have recently been found in many tissues outside the lens, and alpha B-crystallin, in particular, behaves in many respects like a small heat-shock protein. The homologous sequences constitute the C-terminal halves of the proteins and probably represent a structural domain with a more variable C-terminal extension. These domains must be responsible for the common structural and functional properties of this protein family. Analysis of the phylogenetic tree and comparison of the biological properties of the various proteins in this family suggest the following scenario for its evolution: The primordial role of the small heat-shock protein family must have been to cope with the destabilizing effects of stressful conditions on cellular integrity. The alpha-crystallin-like domain appears to be very stable, which makes it suitable both as a surface antigen in parasitic organisms and as a long-living lens protein in vertebrates. It has recently been demonstrated that, like the other heat-shock proteins, the alpha-crystallins and small heat-shock proteins function as molecular chaperones, preventing undesired protein-protein interactions and assisting in refolding of denatured proteins. Many of the small heat-shock proteins are differentially expressed during normal development, and there is good evidence that they are involved in cytomorphological reorganizations and in degenerative diseases. In conjunction with the stabilizing, thermoprotective role of alpha-crystallins and small heat-shock proteins, they may also be involved in signal transduction. The reversible phosphorylation of these proteins appears to be important in this respect.