Comparison of the nucleotide sequences of a yeast gene family. I. Distribution and spectrum of spontaneous base substitutions

The nucleotide sequences of closely related members of a gene family can be used to investigate spontaneous mutations. Here we analyse the sequences of different yeast invertase genes which are more than 93% identical in the coding region and share some very similar, but not identical sequences in the noncoding flanking regions. Since all except one of the invertase genes are active, most of the base substitutions are silent. Within the coding region the base substitutions are unevenly distributed, indicating that parts of the genes were homogenized, probably via gene conversion. Transitions occurred more frequently than transversions in both, coding and noncoding regions. In the coding region pyrimidine transitions were the most abundant event due to silent changes mainly in the third codon position. In the noncoding region pyrimidine and purine transitions were found at equal frequencies. Transversions inverting base pairs (A-T and G-C) outnumber transversions changing base pairs (A-C and G-T). While the spectrum of mutations in the coding region is influenced by selective pressure to maintain the amino acid sequence, the spectrum in the noncoding region may be much less affected by selective pressure.     

Evolution of genes for allelic and isotypic forms of immunoglobulin kappa chains and of the genes for T-cell receptor beta chains in rabbits

New insights into the evolution of the families of genes encoding immunoglobulins and T-cell receptors of rabbits (Oryctolagus cuniculus) have come from molecular genetic studies. In contrast to human and mouse, rabbits were shown to have two genes for the constant region of immunoglobulin light chains (C kappa 1 and C kappa 2 isotypes) and complex allelic variants of K1 (allotypes). Although K1 allotype protein sequences differed at up to 41% of the amino acid positions, 3' untranslated, 5', and 3' flanking regions were conserved, and in the coding regions 78-80% of the codons with differences had replacement changes. Proportions of silent changes and changes in noncoding regions were comparable. Thus, in spite of their markedly different protein sequences, the K1b4, b5, and b9 allotypes appeared to be products of allelic genes. Molecular genetic analyses suggested that they may have undergone rapid divergence after an ancestral K2-like gene duplicated. Some rabbits were found to have two similar T-cell receptor C beta genes as do humans and many strains of mice, but others appeared to have three different C beta. In addition, we found allotypic forms of C beta. Some of the C beta allotypic differences occurred at positions where analogous C kappa allotypic differences were found. We also found V beta in mouse and human that were more similar to rabbit V beta than closely linked rabbit genes were to each other. This contrasts with rabbit immunoglobulin VH gene sequences that reflect concerted evolution. The data suggested that T-cell receptor V beta genes duplicated prior to mammalian radiation.     

Rat kappa-chain J-segment genes: two recent gene duplications separate rat and mouse

We have cloned DNA segments containing the Jk genes from LOUVAIN rat liver, and have determined their nucleotide sequence. Seven readily identifiable Jk-coding regions (six expressible) are evident in the rat, compared with five in the mouse (four expressible). The two additional J segments in the rat appear to be the result of two sequential gene duplications occurring since the divergence of rats and mice. The first involved a homologous but unequal crossing-over in a 14 bp region spanning the 3' end of the coding region of J1 and J2. The second involved a crossing-over following unequal pairing of the two newly duplicated regions. We propose that the probability of a second duplication was greatly increased following the first as a result of the increased target for unequal pairing (370 bp of good homology versus 27 bp in the original pairing). Comparisons of rat and mouse J genes show a surprisingly high degree of sequence conservation, both inside and outside the coding regions, similar to the pattern we reported previously for the kappa constant-region gene. This provides additional evidence that constraints exist on the nucleotide sequences of these genes independent of the function of the encoded proteins.     

Two V kappa germ-line genes related to the GAT idiotypic network (Ab1 and Ab3/Ab1') account for the major subfamilies of the mouse V kappa-1 variability subgroup

V kappa Ig germ-line genes have been isolated from recombinant clones prepared in separate libraries constructed from adult BALB/c liver DNA. Three different clones that strongly hybridized with a V kappa-GAT-specific probe were completely characterized and sequenced. All three genes exhibited common characteristic features in their sequences encompassing the 5' to the 3' noncoding region, with coding sections 95% homologous. A comparison with other V kappa genes shows that the size of the first intron is variability subgroup specific. Moreover, a direct correlation exists between the size of this intron and the entire length of the coding region. Nucleotide sequences of these genes were compared with V kappa chains expressed at the Ab1 and Ab1' levels of the GAT idiotypic network: Ag----Ab1----Ab2----Ab3 (Ab1'). K1A5 and K5.1 genes account for V kappa chains in Ab1 and Ab1' hybridomas, respectively. The high conservation of Ab1' sequences in light chain was also recently reported for the heavy chains, suggesting that immunization with Ab2 (anti-idiotypic) antibodies preferentially stimulates the direct expression of germ-line genes. K5.1 and K1A5 genes belong to the V kappa-1 variability subgroup and encode, without any amino acid substitution, V kappa domain in myeloma TEPC 105 and MOPC 467, which are V kappa-1A and V kappa-1C subgroup prototypes, respectively. These genes are extensively used in different mouse strains and in a number of antibodies of discrete specificities, such as anti-GAT, anti-DNP, anti-flagellin, anti-phosphorylcholine, anti-digoxin, anti-phenyloxazolone, and anti-DNA.     

Importance of both the coding and the segment-specific noncoding regions of the influenza A virus NS segment for its efficient incorporation into virions

The genome of influenza A virus consists of eight single-strand negative-sense RNA segments, each comprised of a coding region and a noncoding region. The noncoding region of the NS segment is thought to provide the signal for packaging; however, we recently showed that the coding regions located at both ends of the hemagglutinin and neuraminidase segments were important for their incorporation into virions. In an effort to improve our understanding of the mechanism of influenza virus genome packaging, we sought to identify the regions of NS viral RNA (vRNA) that are required for its efficient incorporation into virions. Deletion analysis showed that the first 30 nucleotides of the 3' coding region are critical for efficient NS vRNA incorporation and that deletion of the 3' segment-specific noncoding region drastically reduces NS vRNA incorporation into virions. Furthermore, silent mutations in the first 30 nucleotides of the 3' NS coding region reduced the incorporation efficiency of the NS segment and affected virus replication. These results suggested that segment-specific noncoding regions together with adjacent coding regions (especially at the 3' end) form a structure that is required for efficient influenza A virus vRNA packaging.     

Structure and chromosome localization of the human eosinophil-derived neurotoxin and eosinophil cationic protein genes: Evidence for intronless coding sequences in the ribonuclease gene superfamily

Human genomic DNAs for the eosinophil granule proteins, eosinophil-derived neurotoxin (EDN) and eosinophil cationic protein (ECP), were isolated from genomic libraries. Alignment of EDN (RNS2) and ECP (RNS3) gene sequences demonstrated remarkable nucleotide similarities in noncoding sequences, introns, and flanking regions, as well as in the previously known coding regions. Detailed examination of the 5'-noncoding regions yielded putative TATA and CAAT boxes, as well as similarities to promoter motifs from unrelated genes. A single intron of 230 bases was found in the 5' untranslated region and we suggest that a single intron in this region and an intronless coding region are features common to many members of the RNase gene superfamily. The RNS2 and RNS3 genes were localized to the q24-q31 region of human chromosome 14. It is likely that these two genes arose as a consequence of a gene duplication event that took place approximately 25-40 million years ago and that a subset of anthropoid primates possess both of these genes or closely related genes.     

Evolution of nucleic acids coding for ribonucleases: the mRNA sequence of mouse pancreatic ribonuclease

The cDNA of mouse pancreatic mRNA has been cloned. After the library was screened with a rat ribonuclease cDNA probe, the positive clones were isolated and sequenced. There were no differences from the previously determined protein sequence. The mRNA codes for a preribonuclease of 149 amino acid residues including a signal peptide of 25 amino acids. The 3' noncoding region has a length of 260 bp, and the total mRNA length is approximately 940 bp. Comparison with the rat pancreatic ribonuclease sequence showed a high rate of nucleotide substitution. Within the coding region, nonsynonymous and synonymous substitution rates are 4.3 X 10(-9) and 15 X 10(-9) nucleotide substitutions/site/year, respectively. The latter value is one of the highest rates observed in the molecular evolution of mammalian nuclear genes. In the signal sequences the synonymous substitution rate is much lower and about the same as the nonsynonymous rate. Signal sequences of other mouse and rat proteins also exhibit little difference between synonymous and nonsynonymous rates. The sequences of rat and mouse pancreatic ribonuclease messengers were compared with those of bovine pancreatic, seminal, and brain ribonuclease. While the 3' noncoding regions of rat and mouse are very similar, as are those of the three bovine messengers, there is no significant similarity between both rodent and the three bovine messengers for the greater part of these regions. There is a duplication of approximately 50 nucleotides in the 3' noncoding region of the bovine messengers, with a region rich in A and C in between. The presence of this structural feature may be correlated with recent gene duplications that have occurred in the bovine genome.     

The Z family, a group of transposed human immunoglobulin V kappa genes

A group of highly homologous transposed human V kappa I genes, which we call the Z family, was characterized. To date four members, ZI-ZIV, comprising about 230 kb, have been analyzed on cosmid clones. The largest region (ZI) has a length of 85 kb. The Z regions show extensive homology to each other according to restriction maps and hybridization data. In each Z region a solitary V kappa I gene was found. No V kappa genes of other subgroups were detected by hybridization. The nucleotide sequence of the ZI gene revealed a non-processed V kappa I pseudogene. Hybridization experiments with DNAs from rodent/human cell hybrids and other experimental data indicate that some and possibly all members of the Z family lie outside of the kappa locus which is located on chromosome 2; they have been transposed to other chromosomes. Because of their separation from the J kappa C kappa gene segment, the Z genes can be classified as pseudogenes independent of their sequences. We postulate that the Z family arose by amplification event(s). The Z regions can also be regarded as a small family of very long repetitive sequences.     

The nucleotide sequence of a 5.5-kilobase DNA segment containing the mouse kappa immunoglobulin J and C region genes

We have determined the nucleotide sequence of a 5.5-kilobase segment of cloned mouse DNA which includes regions encoding two parts of the mouse kappa immunoglobulin gene: the J regions (amino acids 96-108 of the kappa chain) and the C region (residues 109-214). This sequence allows us to rule out interruptions in the germline constant region coding segment as well as the presence of additional functional J genes in the sequenced DNA segment, although two weak homologies to J regions have been found. The complete sequence also allows us to identify a single occurrence of the heptanucleotide palindrome thought to play a role in V/J joining. This palindrome, midway between J and C regions, is the site of aberrant joining in the plasmacytoma MPC11 and may be a target for such aberrant recombination of other kappa genes. In addition, computer analysis of the J sequences suggests that those closest to the C region arose by the most recent duplication event.     

The nucleotide sequence of a human immunoglobulin C gamma1 gene

We report the nucleotide sequence of a gene encoding the constant region of a human immunoglobulin gamma 1 heavy chain (C gamma 1). A comparison of this sequence with those of the C gamma 2 and C gamma 4 genes reveals that these three human C gamma genes share considerable homology in both coding and noncoding regions. The nucleotide sequence differences indicate that these genes diverged from one another approximately 608 million years ago. An examination of hinge exons shows that these coding regions have evolved more rapidly than any other areas of the C gamma genes in terms of both base substitution and deletion/insertion events. Coding sequence diversity also is observed in areas of CH domains which border the hinge.     

Molecular comparison of the negative-acting nitrogen control gene,nmr, inNeurospora crassa and otherNeurospora and fungal species

In Neurospora crassa, the expression of unlinked structural genes which encode nitrogen catabolic enzymes is subject to genetic and metabolic regulation. The negative-acting nmr regulatory gene appears to play a role in nitrogen catabolite repression. Using the N. crassa nmr gene as a probe, homologous sequences were identified in a variety of other filamentous fungi. The polymerase chain reaction was used to isolate the nmr-like gene from the exotic Mauriceville strain of N. crassa and from the two related species, N. intermedia and N. sitophila. Sequence comparisons were carried out with a 1.7-kb DNA segment which includes the entire coding region of nmr plus 5' and 3' noncoding sequences. The size of the nmr coding region was identical in all three Neurospora species. Approximately 30 nucleotide base substitutions were found in the coding region of the nmr gene of each of the sister species when compared to the standard N. crassa sequence. However, most of the base changes occurred in third codon positions and were silent. The NMR proteins of N. sitophila and of N. intermedia display only three and four amino acid substitutions, respectively, from the N. crassa protein. Two regions of high variability, which include deletions and insertions of bases, were found in the 5' and 3' noncoding regions of the gene.     

DNA variation at the rp49 gene region of Drosophila simulans: evolutionary inferences from an unusual haplotype structure.

An approximately 1.3-kb region including the rp49 gene plus its 5' and 3' flanking regions was sequenced in 24 lines of Drosophila simulans (10 from Spain and 14 from Mozambique). Fifty-four nucleotide and 8 length polymorphisms were detected. All nucleotide polymorphisms were silent: 52 in noncoding regions and 2 at synonymous sites in the coding region. Estimated silent nucleotide diversity was similar in both populations (pi = 0.016, for the total sample). Nucleotide variation revealed an unusual haplotype structure showing a subset of 11 sequences with a single polymorphism. This haplotype was present at intermediate frequencies in both the European and the African samples. The presence of such a major haplotype in a highly recombining region is incompatible with the neutral equilibrium model. This haplotype structure in both a derived and a putatively ancestral population can be most parsimoniously explained by positive selection. As the rate of recombination in the rp49 region is high, the target of selection should be close to or within the region studied.     

Silent nucleotide substitutions and G+C content of some mitochondrial and bacterial genes

Summary The G+C content of DNA varies widely in different organisms, especially microorganisms. This variation is accompanied by changes in the nucleotide composition of silent positions in codons. (Silent positions are defined and explained in the text.) These changes are mostly neutral or near neutral, and appear to result from mutation pressure in the direction of increasing either A+T (AT pressure) or G+C(GC pressure) content. Variations in G+C content are also accompanied by substitutions at replacement positions in codons. These substituions produce changes in the amino acid content of homologous proteins. The examples studied were genes for 13 mitochondrial proteins in five species, and A and B genes for bacterial tryptophan synthase in four species.In microorganisms, varying AT and GC mutational pressures, presumably resulting from shifts in the DNA polymerase system, exert strong effects on molecular evolution by changing the G+C content of DNA. These effects may be greater than those of random drift. The effects of GC pressure on silent substitutions in the systems examined are several times as great as the effects on replacement substitutions.GC pressure is exerted on noncoding as well as coding regions in mitochondrial DNA. This is shown by the close correlation (correlation coefficient, 0.99) of the G+C content of the noncoding D loop of mitochondria with the G+C content of silent positions in the corresponding mitochondrial genes.     

Molecular evolution of the 5'-flanking regions of the duplicated Amy genes in Drosophila melanogaster species subgroup

The nucleotide sequences of the 5'-flanking regions of the duplicated Amy genes in eight sibling species belonging to the melanogaster species subgroup are analyzed. In Drosophila melanogaster, a region of about 450 bp immediately upstream of the translation initiation site of the two paralogous genes (the proximal and distal genes) has sequence similarities. However, we could not detect any significant sequence similarity in the region more upstream than -450. This result indicates that the coding regions of the ancestral Amy gene were duplicated together with 450 bp of the 5'-flanking region as one unit. Multiple alignment of these 450-bp sequences in the proximal and distal genes of all eight species revealed a mosaic pattern of highly conserved and divergent regions. The conserved regions included almost all the putative regulatory elements identified in previous analyses of the sequences. A phylogenetic analysis of the aligned sequences shows that these 450-bp sequences are clustered into the proximal and the distal groups. As a whole, the divergence between groups in this region is very large in contrast to that in the coding regions. Based on the divergence between groups, the 450-bp region is divided into two subregions. We found that the ratios of the divergence between groups to that within groups differ in the two subregions. From these observations, we discuss a possibility of positive selection acting on the subregion immediately upstream of the Amy coding region to cause divergence of regulatory elements of the paralogous genes.      

Nucleotide sequences of switch regions of immunoglobulin C epsilon and C gamma genes and their comparison

Immunoglobulin class switch involves a unique recombination event that takes place at the switch (S) region which is located 5' to each constant region (C) gene of the heavy (H) chain. For example, differentiation of the B lymphocyte from a mu-chain producer to an epsilon-chain producer is mediated by the switch recombination between the S mu and S epsilon regions. In order to elucidate the molecular mechanism for the switch recombination, we have determined nucleotide sequences surrounding the class switch recombination sites of the C epsilon and C gamma 3 genes and those in the 5' flanking regions of the C gamma 2a and C delta genes. The results indicate that the 5' flanking regions of all the CH genes except for the C delta gene contain the S regions which comprise tandem repetition of short unit sequences in agreement with the previous analyses of the S gamma 1, S gamma 2b, S mu, and S alpha regions. Comparison of the nucleotide sequences of all the S regions revealed that length as well as nucleotide sequences of the S regions vary among different classes of the CH gene, but they share short common sequences, (G)AGCT and TGGG(G). The nucleotide sequence of the S mu region is homologous to those of the other S regions in the decreasing order of the S epsilon, S alpha, S gamma 3, and (S gamma 1, S gamma 2b, s gamma 2a) regions. We have compared the nucleotide sequences immediately adjacent to the recombination sites of seven rearranged genes and have always fund tetranucleotides TGAG and/or TGGG, except for one case. Such tetranucleotides may constitute a part of the recognition sequence of a putative recombinase. These results provide further support for our previous proposal that the switch recombination may be facilitated by short common sequences dispersed in all the S regions.     

The cDNA and protein sequences of mouse lactate dehydrogenase B. Molecular evolution of vertebrate lactate dehydrogenase genes A (muscle), B (heart) and C (testis)

Mouse lactate dehydrogenase-B cDNAs were isolated from cDNA libraries of macrophage (ICR strain) and thymus (F1 hybrid of C57BL/6 and CBA strains), and their nucleotide sequences determined. The lactate dehydrogenase-B cDNA insert of thymus clone mB188 consists of the protein-coding sequence (1002 nucleotides), the 5' (46 nucleotides) and 3' (190 nucleotides) non-coding regions, and poly(A) tail (19 nucleotides), while macrophage clone mB168 contains a partial lactate dehydrogenase cDNA insert from codon no. 55 to the poly(A) tail. Seven silent nucleotide substitutions at codon no. 142, 143, 186, 187, 241, 285 and 292, as well as a single nucleotide change in the 3' non-coding region, were found between these different strains of mice. The predicted sequence of 333 amino acids, excluding initiation methionine, was confirmed by sequencing and/or compositional analyses of a total of 103 (31%) amino acids from tryptic peptides of mouse lactate dehydrogenase-B protein. The nucleotide sequence of the mouse coding region for lactate dehydrogenase B shows 86% identity with that of the human isoenzyme, and only eight of the 139 nucleotide differences resulted in amino acid substitutions at residues 10, 13, 14, 17, 52, 132, 236 and 317. The rates of nucleotide substitutions at synonymous and nonsynonymous sites in the mammalian lactate dehydrogenase genes are calculated. The rates of synonymous substitutions for lactate dehydrogenase genes A (muscle) and B (heart) are considerably higher than the average rate computed from human and rodent genes. The rates of nonsynonymous substitutions for lactate dehydrogenase genes A (muscle) and B (heart), particularly the latter, are highly conservative. The rates of synonymous and nonsynonymous substitutions for the lactate dehydrogenase-C gene are about the same as the average rates for mammalian genes. A phylogenetic tree of vertebrate lactate dehydrogenase protein sequences is constructed. In agreement with the previous results, this analysis further indicates that lactate dehydrogenase-C gene branched off earlier than did lactate dehydrogenase-A and lactate dehydrogenase-B genes.