Multisequence comparisons in protein coding genes

A very powerful method for detecting functional constraints operative in biological macromolecules is presented. This method entails performing a base permanence analysis of protein coding genes at each codon position simultaneously in different species. It calculates the degree of permanence of subregions of the gene by dividing it into segments,c codons long, counting how many sites remain unchanged in each segment among all species compared. By comparing the base permanence among several sequences with the expectations based on a stochastic evolutionary process, gene regions showing different degrees of conservation can be selected. This means that wherever the permanence deviates significantly from the expected value generated by the simulation, the corresponding regions are considered “constrained” or “hypervariable”. The constrained regions are of two types: α and β. The α regions result from constraints at the amino acid level, whereas the β regions are those probably involved in “control” processing. The method has been applied to mitochondrial genes coding for subunit 6 of the ATPase and subunit 1 of the cytochrome oxidase in four mammalian species: human, rat, mouse, and cow. In the two mitochondrial genes a few regions that are highly conserved in all codon positions have been identified. Among these regions a sequence, common to both genes, that is complementary to a strongly conserved region of 12S rRNA has been found. This method can also be of great help in studying molecular evolution mechanisms.     

G蛋白亚单位基因家族研究进展

G蛋白由α、β、γ三个亚单位组成异源三聚体。目前已发现16个α、6个β和12个γ基因。G蛋白亚单位基因家族相当保守并且原始,几乎所有G蛋白基因外显子-内含子连接均遵从GT-AG规则,并且各亚单位基因编码区内含子结构和位置显示出很高的保守性。多数G蛋白基因具有持家基因的特点。G蛋白基因在基因组中的分布存在着丛集的倾向,有5对α基因呈二联串连排列。     

Distinct structural features of the alpha and beta subunits of nitrogenase molybdenum-iron protein of Clostridium pasteurianum: an analysis of amino acid sequences

Nitrogenase is composed of two separately purified proteins, a molybdenum-iron (MoFe) protein and an iron (Fe) protein. Structural genes (nifD and nifK) encoding alpha and beta subunits of the MoFe protein of Clostridium pasteurianum (Cp) have been cloned and sequenced. The deduced amino acid sequences were analyzed for structures that could be related to the unique properties of the Cp protein, particularly its low capacity to form an active enzyme with a heterologous Fe protein. Cp nifK is located immediately downstream from Cp nifD, with the start codon of nifK overlapping by one base with the stop codon of nifD. An open reading frame following nifK was identified as nifE. The amino acid sequence deduced from nifK encompasses the partial amino acid sequences previously reported from the isolated beta subunit. Cp nifK encodes a polypeptide of 458 amino acid residues (Mr 50 115) whose amino-terminal region is about 50 residues shorter than the otherwise conserved corresponding polypeptides from four other organisms. In contrast, Cp alpha subunit (nifD product) contains an additional stretch of 50 amino acid residues in the 380-430 region, which is unique to the Cp protein. It therefore appears that the combined size of the alpha and beta subunits could be important to nitrogenase function. An analysis of the predicted secondary structure from the amino acid sequence of each subunit from three species (C. pasteurianum, Azotobacter vinelandii, and Rhizobium japonicum) further revealed structural features, including regions adjacent to some of the conserved cysteine residues, differentiating the Cp MoFe protein from others. These different regions may be further tested for correlation with distinct properties of Cp nitrogenase.     

Cytochrome oxidase subunit II gene in mitochondria of Oenothera has no intron

The cytochrome oxidase subunit II gene has been localized in the mitochondrial genome of Oenothera berteriana and the nucleotide sequence has been determined. The coding sequence contains 777 bp and, unlike the corresponding gene in Zea mays, is not interrupted by an intron. No TGA codon is found within the open reading frame. The codon CGG, as in the maize gene, is used in place of tryptophan codons of corresponding genes in other organisms. At position 742 in the Oenothera sequence the TGG of maize is changed into a CGG codon, where Trp is conserved as the amino acid in other organisms. Homologous sequences occur more than once in the mitochondrial genome as several mitochondrial DNA species hybridize with DNA probes of the cytochrome oxidase subunit II gene.     

Primers for the amplification of nuclear introns in sea stars of the family Asteriidae

We describe polymerase chain reaction primers that consistently amplify three intron regions of approximately 2 kb in total length for two nuclear protein‐coding genes (ATP synthase beta subunit and elongation factor‐1 alpha subunit) in sea stars of the family Asteriidae. The introns are moderately polymorphic at the species level (average within‐species percentage of site differences = 0.42%, range 0–1.44%), are evolving at about 29% as fast as mitochondrial sequences in the same species and are alignable at the genus or family level, making them suitable for phylogenetic and population genetic analyses.     

Nucleotide sequence of the pntA and pntB genes encoding the pyridine nucleotide transhydrogenase of Escherichia coli

A 3240-base-pair DNA fragment spanning the pyridine nucleotide transhydrogenase (pnt) genes of Escherichia coli has been sequenced. The sequence contains two open-reading frames, pntA and pntB of 1506 and 1386 base pairs, coding for the transhydrogenase alpha and beta subunits, respectively. The coding sequences are preceded by a promoter-like structure and are most likely co-transcribed. Each coding sequence is preceded by a Shine-Dalgarno sequence. The amino-terminal amino acid sequences were determined from the purified alpha and beta subunits of the transhydrogenase. These sequences agree with those predicted from the nucleotide sequences of the pntA and pntB genes. The predicted relative molecular masses of 53906 (alpha) and 48667 (beta) are close to the values obtained by analysis of the subunits by sodium dodecyl sulfate/polyacrylamide gel electrophoresis. Several hydrophobic regions large enough to span the cytoplasmic membrane were observed in each subunit. These results indicate that transhydrogenase is an intrinsic membrane protein.     

A discontinuous small subunit ribosomal RNA in Tetrahymena pyriformis mitochondria

We show here that in the mitochondria of Tetrahymena pyriformis, the small subunit (SSU) rRNA is discontinuous, being comprised of two separate components which we term "alpha" (a novel low molecular weight RNA, approximately equal to 200 nucleotides long) and "beta" (a previously described 14 S RNA). The SSU alpha rRNA has been sequenced in its entirety; it represents the immediate 5'-terminal domain of conventional SSU rRNA. The sequences at the ends of the SSU beta rRNA have also been determined; they show that this molecule corresponds to the 3'-terminal 7/8 of conventional SSU rRNA. A 2.5-kilobase pair XbaI restriction fragment of T. pyriformis mitochondrial DNA which contains the SSU alpha and SSU beta rRNA genes was cloned and its complete nucleotide sequence was determined. This revealed that the genes encoding the two segments of SSU rRNA are separated by a 54-base pair (A + T)-rich spacer. The alpha and beta sequences can be fitted to a generalized secondary structure model for eubacterial 16 S rRNA, with the two RNA species associating through long range interactions to form base-paired regions characteristic of SSU rRNA. In this model, the spacer is situated in a region of pronounced primary and secondary structural variation among SSU rRNAs. The significance of these findings with respect to rRNA biosynthesis and processing and the possible evolutionary relationship between spacers and variable regions in rRNA genes is discussed.     

Codon usage and base composition inRickettsia prowazekii

Codon usage and base composition in sequences from the A + T-rich genome ofRickettsia prowazekii, a member of the alpha Proteobacteria, have been investigated. Synonymous codon usage patterns are roughly similar among genes, even though the data set includes genes expected to be expressed at very different levels, indicating that translational selection has been ineffective in this species. However, multivariate statistical analysis differentiates genes according to their G + C contents at the first two codon positions. To study this variation, we have compared the amino acid composition patterns of 21R. prowazekii proteins with that of a homologous set of proteins fromEscherichia coli. The analysis shows that individual genes have been affected by biased mutation rates to very different extents: genes encoding proteins highly conserved among other species being the least affected. Overall, protein coding and intergenic spacer regions have G + C content values of 32.5% and 21.4%, respectively. Extrapolation from these values suggests thatR. prowazekii has around 800 genes and that 60–70% of the genome may be coding. Correspondence to: S.G.E. Andersson     

Complete structure of the mouse mast cell receptor for IgE (Fc epsilon RI) and surface expression of chimeric receptors (rat-mouse-human) on transfected cells

The high affinity receptor for IgE (Fc epsilon RI) found on mast cells and basophils is a tetrameric complex of a single alpha subunit, a single beta subunit, and two identical gamma subunits. The genes for the three subunits of mouse Fc epsilon RI have now been cloned from the mast cell line, PT18. When compared at the DNA level, the rat and mouse subunits are similarly conserved. However, at the protein level the homology between mouse and rat alpha is surprisingly low (71% identities) especially in the cytoplasmic regions (57% identities) which are of different length (25 and 20 residues, respectively). By contrast the beta and gamma are homogeneously conserved between mouse and rat (83 and 93% identities, respectively). The consensus amino acid sequence of the alpha subunit derived from three species (rat, mouse, and human) shows that the cytoplasmic tail diverges to the same extent as the leader peptide. Conversely, the transmembrane domain of the alpha is highly conserved and contains 10 consecutive residues that are identical. Comparisons between mouse Fc epsilon RI and other mouse proteins reveal regions of high homology between the alpha subunit and Fc gamma RIIa and between the gamma subunit and the zeta chain of the T cell receptor. Cells transfected with the alpha gene express the alpha subunit on their surface very inefficiently. Efficient expression is only achieved after co-transfection of the three rodent genes or of the human alpha gene together with the rodent gamma without apparent need for beta. The subunits are completely interchangeable upon transfection so that various chimeric mouse-rat-human receptors can be expressed.     

Mechanism of expression of the overlapping genes of Bacillus subtilis aspartokinase II

The mechanism of expression of the overlapping genes that encode the alpha and beta subunits of aspartokinase II of Bacillus subtilis was studied by specific mutagenesis of the cloned coding sequence. Escherichia coli or B. subtilis VB31 (aspartokinase II-deficient), transformed with plasmids carrying either a deletion of the translation start site and about one-half of the coding region for the larger alpha subunit or a frameshift mutation early in the alpha subunit coding region, produced the smaller beta subunit in the absence of alpha subunit synthesis, indicating that beta subunit is not derived from alpha subunit and that its synthesis does not depend on the alpha subunit translation initiation site. The beta subunit translation start site was identified by oligonucleotide-directed mutagenesis of the putative translation start codon. Modification of the nucleotide sequence encoding methionine residue 247 of the alpha subunit from ATG to either TTA or AAT (but not GTG) abolished beta subunit synthesis but had no effect on the production of alpha subunit. This observation is consistent with peptide chain initiation by N-formylmethionine, which specifically requires an ATG or GTG sequence, and indicates that translation of the beta subunit starts at a site corresponding to Met247 of the alpha subunit. Initial studies on the function of the aspartokinase II subunits, using E. coli as a heterologous host, showed that beta subunit was not essential for the expression of the catalytic function of aspartokinase, measured in vitro and in vivo, nor for its allosteric regulation by L-lysine. Whether the beta subunit has a function specific to B. subtilis needs to be explored in a homologous expression system.     

DNA sequence of the tryptophan synthase genes of Pseudomonas putida

Genes encoding the 2 subunits of tryptophan synthase in Pseudomonas putida have been identified and cloned by their similarity to the corresponding genes in Pseudomonas aeruginosa. The deduced amino acid sequences were confirmed by comparison with regions ascertained earlier by protein sequencing. The Pseudomonas amino acid sequences are 85% identical for the beta subunit and 70% identical for the alpha subunit. These sequences are compared to those of Salmonella typhimurium, where the structure is known from X-ray crystallography. Although amino acid conservation drops to 54% and 36% for the beta and alpha subunits, only 3 single residue gaps are required to maintain alignment throughout and most of the residues identified as important for catalysis or cofactor binding are conserved. The 23 residues surrounding the beta chain lysine that enters into a Schiff base linkage with the pyridoxal phosphate cofactor are compared in 13 species, including representatives from the eukaryotic and both prokaryotic kingdoms; appreciable conservation is apparent. The approximately 100 base pairs separating the trpB gene from its divergently transcribed activator gene are similar in the 2 pseudomonads, but do not resemble those of any other bacterium or fungus studied to date.     

Cloning, expression and characterization of phycoerythrin gene from Ceramium boydenn.

Phycobiliproteins function as a major light harvesting protein-pigment complex in the cyanobacteria and the eukaryotic algae. Phycoerythrin (PE) is a kind of phycobiliproteins, widely located in all rhodophytes, some species of cyanobacteria and cryptophytes, and different ecotypes of Prochlorococcus populations. PeBA encoding beta and alpha subunits of PE from Ceramium boydenn was cloned and sequenced in this research. A peBA specific PCR primer was synthesized, based on the peBA gene conserved sequences. The beta subunit encoding gene (peB) contained an open reading frame of 534 bp, while the alpha subunit (peA) was 495 bp. Recombinant expression plasmid pET-peAB was constructed and expressed in Escherichia coli BL21. The molecular weight of expressive product of peB and peA was about 23.3 and 18.2 KD, respectively. Results of codon usage analysis show that G + C content is heterogeneous among different groups of PE and spacers have dramatically lower G + C contents than coding regions. Also there is a high variance in G + C content among sequences at the third position sites. It is also found in this paper that several sequence regions, which might reflect functional or structural requirements of the PE organization, and several residues known for their functional importance are conserved in almost all the sequences.     

The mitochondrial genome of Ruspolia dubia (Orthoptera: Conocephalidae) contains a short A+T-rich region of 70 bp in length.

The complete sequence (14 971 bp) of the Ruspolia dubia mitochondrial genome was determined and annotated. The genome contains the gene content, base composition, and codon usage typical of metazoan mitochondrial genomes. All 37 genes are conserved in the positions observed most frequently in insect mitochondrial genome structures. The secondary structures of both small subunit and large subunit rRNA were predicted. The most unusual features found were the initiation codon (TTA) of COI and a short A+T-rich region of 70 bp in length. In addition, a short, highly conserved polythymidine stretch that was previously described in Orthoptera and Diptera was also present in the A+T-rich region.     

Intergenic DNA sequences flanking the pseudo alpha globin genes of human and chimpanzee.

We have determined the sequence of 2400 base pairs upstream from the human pseudo alpha globin (psi alpha) gene, and for comparison, 1100 base pairs of DNA within and upstream from the chimpanzee psi alpha gene. The region upstream from the promoter of the psi alpha gene shows no significant homology to the intergenic regions of the adult alpha 2 and alpha 1 globin genes. The chimpanzee gene has a coding defect in common with the human psi alpha gene, showing that the product of this gene, if any, was inactivated before the divergence of human and chimpanzee. However the chimpanzee gene contains a normal ATG initiation codon in contrast to the human gene which has GTG as the initiation codon. The psi alpha genes of both human and chimpanzee are flanked by the same Alu family member. The structure and position of this repeat have not been altered since the divergence of human and chimpanzee, and it is at least as well conserved as its immediate flanking sequence. Comparing human and chimpanzee, the 300 bp Alu repeat has accumulated only two base substitutions and one length mutation; the adjacent 300 bp flanking region has accumulated five base substitutions and twelve length mutations.