Topogenesis of microbody enzymes: a sequence comparison of the genes for the glycosomal (microbody) and cytosolic phosphoglycerate kinases of Trypanosoma brucei.

To determine how microbody enzymes enter microbodies, we are studying the genes for cytosolic and glycosomal (microbody) isoenzymes in Trypanosoma brucei. We have found three genes (A, B and C) coding for phosphoglycerate kinase (PGK) in a tandem array in T. brucei. Gene B codes for the cytosolic and gene C for the glycosomal isoenzyme. Genes B and C are 95% homologous, and the predicted protein sequences share approximately 45% amino acid homology with other eukaryote PGKs. The microbody isoenzyme differs from the cytosolic form and other PGKs in two respects: a high positive charge and a carboxy-terminal extension of 20 amino acids. Our results show that few alterations are required to redirect a protein from cytosol to microbody. From a comparison of our results with the unpublished data for three other glycosomal glycolytic enzymes we infer that the high positive charge represents the major topogenic signal for uptake of proteins into glycosomes.     

Cloning and characterization of a variant surface glycoprotein expression site from Trypanosoma equiperdum.

Variant surface glycoprotein (VSG) genes of African trypanosomes are expressed when they are inserted into one of several telomere-linked expression sites. We cloned and characterized an 11-kilobase (kb) DNA fragment located upstream of an expressed VSG gene. A DNA sequence of 1.8 kb that is located immediately upstream of the inserted VSG gene contains sequences homologous to the 76-base-pair repeats described as being upstream of VSG genes in Trypanosoma brucei (D. A. Campbell, M. P. Van Bree, and J. C. Boothroyd, Nucleic Acids Res. 12:2759-2774). There are no such sequences elsewhere in the 11-kb cloned region. Southern blot analysis using probes from the cloned region revealed multiple unlinked copies of the same or very similar regions. At least three of these are located near telomeres, and two have been shown to be used for the expression of known Trypanosoma equiperdum VSG genes. Like VSG genes, the upstream sequences themselves can be duplicated and deleted. The choice of expression site to be used by a duplicated VSG gene is nonrandom; the site used for expression of the parental VSG gene is strongly favored for use in the daughter variant. Furthermore, even when the parental expression site is not used, the VSG gene occupying it is replaced. Thus, an active expression site is a preferential target for gene conversion in the next variation event.     

Early Expression of a Trypanosoma brucei VSG Gene Duplicated from an Incomplete Basic Copy

Intrachromosomal variant surface glycoprotein (VSG) genes in Trypanosoma brucei are expressed by a mechanism involving gene conversion. The 3'boundary of gene conversion is usually within the last 130 bp of the VSG gene, a region of partially conserved sequences. We report here the loss of the predominant telomeric A VSG gene in the cloned variant antigenic type (VAT) 5^A3, leaving only an intrachromosomal A VSG gene (the A-B gene). The nucleotide sequence of the A-B VSG gene reveals that it lacks the normal VSG 3' sequence. Surprisingly, we find cells expressing this A-B VSG gene in relapse populations arising from VAT 5^A3. Since the A VSG mRNAs from these cells have a normal 3' sequence, the incomplete A-B VSG gene must be expressed via a partial gene conversion that supplies the functional 3'end. Although the A-B VSG gene is no longer predominant like the telomeric A VSG gene, it is still expressed more frequently than other intrachromosomal VSG genes, suggesting that factors other than a telomeric location determine whether a VSG gene is expressed early in a serodeme.     

Conservative evolution in duplicated genes of the primate Class I ADH cluster

Humans have seven alcohol dehydrogenase genes (ADH) falling into five classes. Three out of the seven genes (ADH1A, ADH1B and ADH1C) belonging to Class I are expressed primarily in liver and code the main enzymes catalyzing ethanol oxidization. The three genes are tandemly arrayed within the ADH cluster on chromosome 4 and have very high nucleotide similarity to each other (exons: >90%; introns: >70%), suggesting the genes have been generated by duplication event(s). One explanation for maintaining similarity of such clustered genes is homogenization via gene conversion(s). Alternatively, recency of the duplications or some other functional constraints might explain the high similarities among the genes. To test for gene conversion, we sequenced introns 2, 3, and 8 of all three Class I genes (total>15.0 kb) for five non-human primates--four great apes and one Old World Monkey (OWM)--and compared them with those of humans. The phylogenetic analysis shows each intron sequence clusters strongly within each gene, giving no evidence for gene conversion(s). Several lines of evidence indicate that the first split was between ADH1C and the gene that gave rise to ADH1A and ADH1B. We also analyzed cDNA sequences of the three genes that have been previously reported in mouse and Catarrhines (OWMs, chimpanzee, and humans) and found that the synonymous and non-synonymous substitution (dN/dS) ratios in all pairs are less than 1 representing purifying selection. This suggests that purifying selection is more important than gene conversion(s) in maintaining the overall sequence similarity among the Class I genes. We speculate that the highly conserved sequences on the three duplicated genes in primates have been achieved essentially by maintaining stability of the hetero-dimer formation that might have been related to dietary adaptation in primate evolution.     

A single gene encodes multiple neuropeptides mediating a stereotyped behavior

Egg laying in Aplysia is characterized by a stereotyped behavioral array which is mediated by several neuroactive peptides. We have sequenced two genes encoding the A and B peptides thought to initiate the egg-laying process, as well as a gene encoding egg-laying hormone (ELH) which directly mediates the behavioral array. The three genes share 90% sequence homology and are representatives of a small multigene family. Each gene encodes a protein precursor in which the active peptides are flanked by internal cleavage sites providing the potential to generate multiple small peptides. Each of the three genes consists of sequences homologous to A or B peptide as well as ELH. Although these genes share significant nucleotide homology, they have diverged such that different member genes express functionally related but nonoverlapping sets of neuroactive peptides in different tissues.     

The tubulin genes of Trypanosoma cruzi

The organization of the alpha- and beta-tubulin genes in the genome of Trypanosoma cruzi have been analysed by Southern blotting using tubulin probes derived from Trypanosoma brucei. The tubulin array appears to be more complex in this organism than in other members of the same family. Some tubulin genes are tightly clustered in an alternating (alpha-beta)n array with a basic repeat unit length of 4.3 kb. However, other pairs of alternating alpha- and beta-tubulin sequences appear to be physically separated from the basic group. This finding indicates that the tubulin gene cluster present in T. cruzi is less perfectly conserved than in T. brucei. T. (Herpetosoma) rangeli is similar to T. (Schizotrypanum) cruzi in its tubulin gene organization whereas most of these genes are tandemly clustered in the genome of T. (Trypanozoon) evansi, with a basic repeat unit length of 3.6 kb as previously described for T. (Trypanozoon) brucei. Two overlapping recombinant clones containing T. cruzi tubulin sequences have been isolated from a genomic cosmid library of T. cruzi epimastigotes using the T. brucei tubulin probes. Partial sequencing of the T. cruzi beta-tubulin gene has confirmed its identity and shows more than 70% homology with the sea urchin, chicken and T. b. rhodesiense beta-tubulin reported gene sequences. Analysis of tubulin gene organization through the parasite life cycle does not show evidence of major rearrangements within the repeat unit. Several T. cruzi strains and cloned lines whilst sharing the 4.3-kb tubulin repeat unit, exhibited very variable tubulin gene organization with tubulin probes. These striking differences in the organization of this structural gene among T. cruzi strains and cloned lines suggest that the heterogeneity previously reported in parasite populations may be related to a very dynamic, diploid genome.     

Structural diversity and organization of three gene families for Kunitz-type enzyme inhibitors from potato tubers ( Solanum tuberosum L.)

In the potato, Kunitz-type enzyme inhibitors are abundant and highly polymorphic small proteins found in tubers. DNA sequence analysis of 1596 unselected ESTs (expressed sequence tags) from mature tubers of the cultivars Provita and Saturna resulted in the identification of 55 different DNA sequences with high sequence similarity to Kunitz-type enzyme inhibitors. The frequency of Kunitz-type inhibitor ESTs in Provita was four times higher than in Saturna tubers, and none of the Provita ESTs was identical to any of the Saturna ESTs. A phenogram constructed from the deduced amino acid sequences of the inhibitors revealed three major homology groups-A, B and C. Group A inhibitors were all derived from Provita ESTs. Inhibitor groups A and B were more similar to each other than to group C inhibitors, and for most members within-group similarity was at least 90%. Non-conservative amino acid substitutions and insertion/deletion polymorphisms suggest functional differentiation between members of the gene family. A minimum of 21 genes for Kunitz-type enzyme inhibitors (six for group A, nine for group B and six for group C) was estimated to exist in the potato genome. Genetic mapping and the identification of BAC (bacterial artificial chromosome) clones containing more than one member of the gene family indicated that most inhibitor genes of groups A, B and C are organized in a cluster that maps to a single region on potato chromosome III.     

Characterization of the DNA duplication-transposition that controls the expression of two genes for variant surface glycoproteins in Trypanosoma brucei

The genome of Trypanosoma brucei carries over a hundred genes coding for different variants of the major surface glycoprotein. Activation of some of these genes is accompanied by a duplication and transposition of the gene (the basic copy) to another region in the genome where it is transcribed. We present here physical maps of the basic and transposition-activated genes for two surface glycoproteins of Trypanosoma brucei, stock 427. In both cases the transposed segment starts 1-2 kb in front of the coding region and ends within the 3'-terminal region of the gene. The DNA segments flanking both transposed genes are indistinguishable and share a 6-kb stretch upstream and a 8-kb stretch downstream of the transposed segment not cut by several restriction endonucleases. The 5' borders of the two transposed segments are homologous and contain sequences present in many copies in the genome. A different repeated sequence has previously been found at the 3' edge of the transposed segment. The replicative transposition may, therefore, involve a unidirectional gene conversion initiated by base pairing between the edges of the transposed sequence and a single expression site elsewhere in the genome.     

African trypanosome glucose transporter genes: organization and evolution of a multigene family

Trypanosoma brucei brucei (EATRO-164) contains a tandem array of six genes encoding a glucose transporter, THT1 (trypanosome hexose transporter), followed by five genes encoding a second isoform, THT2. Two distinct clusters containing THT1 and THT2 genes have been identified in the EATRO-164 clone and in most other African trypanosome clones analyzed. Analysis of progeny from crosses between clones of T. b. brucei displaying polymorphism in THT1 copy number per cluster suggests that the two clusters of THT genes are present on homologous chromosomes. In addition, analysis of 30 African trypanosome clones revealed a high degree of polymorphism in THT1 copy number per cluster. Sequence comparison of five THT1 and two and one-half THT2 unit repeats, present within a 20-kb region, provided information about the genesis and evolution of the THT multigene family. The most divergent regions between THT1 and THT2 unit repeats probably arose from insertion of DNA fragments into an ancestral THT region. Genes of each of the different families are almost identical, and there are large regions of identity shared between THT1 and THT2 members. A mosaic copy containing most of a THT1 gene with the 3' extremity of a THT2 gene is found within the cluster. These results suggest that THT1 and THT2 arose by modification (insertion, mutation, or conversion) of duplicated ancestral genes. Functional constraints and homologous recombination may be evoked to explain the maintenance of the conserved sequences of THT1 and THT2.      

Comparison of the maxicircle (mitochondrial) genomes of Leishmania tarentolae and Trypanosoma brucei at the level of nucleotide sequence

The entire 16.7-kilobase (kb) transcribed region of the Leishmania tarentolae maxicircle was compared to the entire 15-kb transcribed region of the Trypanosoma brucei maxicircle at the nucleotide sequence level by dot matrix analysis and by alignments of individual genes. The L. tarentolae NADH dehydrogenase subunit 1 (ND1) gene was identified in a newly obtained 2.9-kb sequence. All but two regions which flank the cytochrome b gene are highly conserved in both species. One 3.1-kb region in L. tarentolae that contains the cytochrome oxidase subunit III (COIII) gene and several open reading frames corresponds to a 2-kb sequence in T. brucei with limited sequence homology that lacks the COIII gene. Another 0.6-kb region that comprises an unidentified open reading frame (open reading frame 12) in L. tarentolae is substituted by a nonhomologous 0.4-kb open reading frame in T. brucei. A short intergenic region between the ND1 gene and the maxicircle unidentified reading frame 1 gene shows limited sequence homology, and the regions between the ND4 and ND5 genes and between the COI and ND4 genes are not conserved. All of the intergenic regions share G + C richness and a similar pattern of G versus C strand bias. 1.8 kb of the L. tarentolae divergent region (DV) and around 3 kb of the T. brucei DV were also obtained. The T. brucei DV sequences were not homologous to the L. tarentolae DV sequence but were organized in a similar fashion with tandem repeats of varying complexity.     

Antigenic variation among group A streptococcal M proteins. Nucleotide sequence of the serotype 5 M protein gene and its relationship with genes encoding types 6 and 24 M proteins

The 1479-base pair (bp) nucleotide sequence of the serotype 5 M protein gene (smp5) from Streptococcus pyogenes contains three distinct types of tandemly repeated sequences, designated A, B, and C. Repeat A (21 bp x 6, in the 5'-half of smp5), shares no homology with the types 6 or 24 M protein genes (Hollingshead, S. K., Fischetti, V. A., and Scott, J. R. (1986) J. Biol. Chem. 261, 1677-1686; Mouw, A. R., Beachey, E. H., and Burdett, V. (1988) J. Bacteriol., in press). Repeat B (75 bp x 3.6, in the center of smp5) is also present in the M6, but not in the M24 gene. Repeat C (105 bp x 2.7, just distal to the B repeats) shares homology with repeats in both the M6 and M24 genes. All three genes share extensive homology in their 3'-halves and in 5' sequences encoding the N-terminal signal peptides, but between these two regions there are highly variable sequences that are responsible for antigenic diversity. These relationships suggest that both intergenic and intragenic recombination has occurred during the evolution of distinct M protein serotypes. All three M proteins contain conserved hydrophobic and proline-rich sequences at their C-terminal ends, suggestive of a membrane anchor and a peptidoglycan spanning region.     

山茶属CHS基因家族的组成和分子进化初探

用PCR方法从4种山茶属(Camellia)(山茶科)(Tlaeaceae)植物的总DNA中分别扩增到CHS基因外显子2的部分序列,经克隆、测序得到16个该基因的序列,这些序列与来自GenBank的该属另一种植物的3个序列及作为外类群的大豆(Glycine max (L)Merr)的2个序列一起进行分析。研究表明,山茶属CHS基因家族在进化过程中已分化为A、B、c三个亚家族,包括A1、A2、A3、B1、B2、C等6类不同的基因成员;其中只有A2类成员为全部被研究的5种植物所共有,而其他5类成员只在部分被研究的植物中发现;所有这些CHS成员具有很高的同源性：在核苷酸水平上同一亚家族内基本上高于90％,不同亚家族间也在78％以上。从推测的氨基酸组成看,山茶属内CHS基因的功能已发生了分化,各类成员的碱基替代率有较大差异;从分子系统发育树和可能的氨基酸组成分析,山茶属具有新功能的基因成员是在经过基因重复后,或是由少数几个位点的突变而成,或是由逐渐积累的突变而形成的。进一步分析认为,该属CHS基因的分化直到近期还在活跃地进行,并且不同种的进化式样有一定的差别,这种不同的进化式样可能是物种形成后受不同环境因素影响而形成的。     

山茶属CHS基因家族的组成和分子进化初探

用PCR方法从4种山茶属(Camellia)(山茶科)(Theaceae)植物的总DNA中分别扩增到CHS基因外显子2的部分序列,经克隆、测序得到16个该基因的序列,这些序列与来自GenBank的该属另一种植物的3个序列及作为外类群的大豆(Glycine max (L.) Merr.)的2个序列一起进行分析.研究表明,山茶属CHS基因家族在进化过程中已分化为A、B、C三个亚家族,包括A1、A2、A3、B1、B2、C 等6类不同的基因成员;其中只有A2类成员为全部被研究的5种植物所共有,而其他5类成员只在部分被研究的植物中发现;所有这些CHS成员具有很高的同源性:在核苷酸水平上同一亚家族内基本上高于90%,不同亚家族间也在78%以上.从推测的氨基酸组成看,山茶属内CHS基因的功能已发生了分化,各类成员的碱基替代率有较大差异; 从分子系统发育树和可能的氨基酸组成分析,山茶属具有新功能的基因成员是在经过基因重复后,或是由少数几个位点的突变而成,或是由逐渐积累的突变而形成的.进一步分析认为,该属CHS基因的分化直到近期还在活跃地进行,并且不同种的进化式样有一定的差别,这种不同的进化式样可能是物种形成后受不同环境因素影响而形成的.     

Heterogeneity of the gamma-globin gene sequences in Japanese individuals: implication of gene conversion in generation of polymorphisms

The linked fetal globin genes (the G gamma- and A gamma-globin genes) were cloned from Japanese individuals with three different haplotypes of the HindIII polymorphisms within the gamma-globin genes. Determination of nucleotide sequences of the segment spanning from IVS2 to the 3' flanking region of each gamma-globin gene revealed that nucleotide differences are located at 43 positions and a stretch of simple GT or GC sequences. Almost half of the nucleotide changes could be accounted for by gene conversion between the G gamma- and A gamma-globin genes. We found that gene conversion had created the SacI polymorphic site just downstream of the A gamma-globin coding region. Association of the SacI polymorphic site with the HindIII polymorphic site suggests that the region containing these two sites was derived from that of the linked G gamma-globin gene through a gene conversion event. The nucleotide sequences obtained here are identical to those of the Caucasoid fetal globin genes of the same haplotypes, with the exception of some sequence changes in the hot spots of mutations. These results indicate that the sequence heterogeneity of the gamma-globin genes can be classified into three major categories according to HindIII haplotypes. The possible mechanisms of generation of the heterogeneity of the gamma-globin gene sequences are discussed.     

Bombyxin: An Insect Brain Peptide that Belongs to the Insulin Family

Bombyxin is a 5 kDa secretory brain peptide that belongs to the insulin family. Bombyxin of the silkmoth Bombyx mori can induce adult development when injected into brain-removed dormant pupae of the saturniid moth Samia cynthia ricini by activating the prothoracic glands to synthesize and release ecdysone. Bombyx bombyxin has been shown to lower the concentration of the major haemolymph sugar, trehalose, and to elevate the trehalase activity in the midgut and muscles in Bombyx, but the doses required to be effective are higher than the amounts in the feeding larvae. The exact physiological function of bombyxin in Bombyx itself is therefore still obscure, but its insulin-like structure suggests it has important roles. Bombyxin comprises a mixture of highly heterogeneous molecular forms whose amino acid sequences have 40% identity with human insulin. The Bombyx bombyxin gene encodes a precursor consisting of the signal peptide, B chain, C peptide, and A chain, in that order from the N terminus. So far, 32 bombyxin genes have been identified in Bombyx, and they are classified into 7 families, A to G, according to their sequence similarity. The bombyxin genes have no introns and cluster in unique distribution patterns. The gene arrangement in the cluster has been classified into three categories: gene pairs, gene triplets, and single genes. Nucleotide sequence analysis indicates that equal and unequal crossings-over and duplications may have generated these unique distribution patterns. The Bombyx bombyxin genes are expressed predominantly in the brain and at low levels in a number of other tissues. Genes of all 7 families are expressed in four pairs of the medial neurosecretory cells of the brain. Detailed examination indicated that only a limited number of genes in the A, B and C family members are expressed and that their expression shows a gene-arrangement-dependent pattern.     

The African trypanosome genome

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

The role of gene conversion in determining sequence variation and divergence in the Est-5 gene family in Drosophila pseudoobscura.

Nucleotide sequences of eight Est-5A and Est-5C genes corresponding to previously sequenced Est-5B genes in Drosophila pseudoobscura were determined to compare patterns of polymorphism and divergence among members of this small gene family. The three esterase genes were also sequenced from D. persimilis and D. miranda for interspecific comparisons. The data provide evidence that gene conversion between loci contributes to polymorphism and to the homogenization of the Est5 genes. For Est-5B, which encodes one of the most highly polymorphic proteins in Drosophila, 12% of the segregating amino acid variants appear to have been introduced via gene conversion from other members of the gene family. Interlocus gene conversion can also explain high sequence similarity, especially at synonymous sites, between Est-5B and Est-5A. Tests of neutrality using interspecific comparisons show that levels of polymorphism conform to neutral expectations at each Est-5 locus. However, McDonald-Kreitman tests based on intraspecific gene comparisons indicate that positive selection on amino acids has accompanied Est-5 gene duplication and divergence in D. pseudoobscura.