Accepted mutations in a gene family: Evolutionary diversification of duplicated DNA

Summary We report and compare the DNA sequences of 14 silkmoth (Antheraea polyphemus) chorion genes, derived from either cDNA or chromosomal DNA clones. Seven of these genes are members of the A multigene family, and seven are members of the B family. Where available, the previously reported (Jones and Kafatos 1980) intronic and extragenic flanking DNA sequences are also considered. Closely related sequences are compared, revealing the types of spontaneous mutations that were fixed during paralogous evolution. Segmental mutations (i.e. mutations other than substitutions) are nearly always interpretable as small duplications or deletions. related to small direct repeats. Segmental mutations are strongly constrained in the coding regions, although they do occur. Nucleotide substitutions also appear to be under selective constraints: relatively few substitutions leading to amino acid replacements are accepted, silent substitutions leading to some codons (especially purine-terminated ones) are disfavored, and different compositional biases are maintained in different parts of the sequences. Other sequence differences can be interpreted as indicative of neutral drift, including most differences in non-coding regions and most T/C transitions in third-base positions. In the non-coding regions, which are thought to be only loosely constrained by selection, transitions are observed more frequently than might be expected: they account for 52% of all substitutions, and they appear to be favored two to threefold over transversions when allowance is made for the skewed base composition of these regions.     

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.     

Comparison of the nucleotide sequences of a yeast gene family. II. Analysis of spontaneous deletions and insertions

We compared the nucleotide sequences of 3 yeast invertase genes in regions where the homology is better than 90%. In the noncoding region 40 gaps of 1-61 bases were found. This is about half as much as the nucleotide substitutions in the same sequences. We grouped the gaps into 5 categories by their length and the characteristics of their sequences. Group I gaps are about 20 nucleotides long and are flanked by repeated sequence of 6 bases which may trigger the deletion of one of the repeats and the sequence between the repeats. Group II gaps are characterized by a small repeated sequence which is missing in one of the invertase genes. Gaps which occur in sequences exclusively made up of one of the 4 bases are summarized in group III. The 4 gaps in group IV do not show any of these sequence characteristics and they are all just one base long. A 61 nucleotide sequence found in only one of the invertase genes seems to be of complex origin. We conclude that small repeated sequences or monotonous sequences are prone to deletion or insertion mutations.     

Excision repair influences the site and strand specificity of sunlight mutagenesis in yeast.

A collection of 384 mutations recovered in a tRNA gene (SUP4-o) following exposure of isogenic excision-repair-proficient (RAD1) or deficient (rad1) strains of the yeast Saccharomyces cerevisiae to sunlight was characterized by DNA sequencing. In each case, greater than 90% of the mutations were single base-pair substitutions with events at G.C pairs constituting most of the changes. However, more than half of these substitutions were transversions in the RAD1 strain whereas transitions predominated in the rad1 strain. Tandem double substitutions were recovered in both strains and the individual changes were exclusively G.C----A.T transitions. The majority of single substitutions, and all tandem double changes, were at base-pairs where the pyrimidine(s) was part of a dipyrimidine sequence and the site specificities were consistent with cyclobutane dimers and/or pyrimidine (6-4) pyrimidone photoproducts contributing to sunlight mutagenesis. Yet, the data also pointed to an important role for lesions that form at G.C pairs and give rise to transversions. Analysis of the strand specificity of sunlight mutagenesis indicated that transitions or transversions at G.C pairs occurred preferentially in SUP4-o at sites where a dipyrimidine or a guanine, respectively, was on the transcribed strand. These biases required a functional excision-repair system.     

Patterns of ribosomal RNA evolution in salamanders

Sequence comparisons are presented for four segments of the large subunit of ribosomal RNA, including divergent domains D7a and D7b, portions of the large divergent domains D2, D3, and D8, and evolutionarily conservative sequences flanking divergent domains. These results resolve phylogenetic relationships among exemplars of seven families of salamanders and the three amphibian orders. Phylogenetic analysis confirms the prediction that divergent domains feature the highest relative rates of base substitution and length variation within the ribosome, but the divergent domains evolve more slowly than nuclear noncoding DNA and the silent sites of structural genes. Base substitutions demonstrate approximately twice as many transitions as transversions and an uneven distribution among sites within the divergent domains but no apparent bias in base composition. Length mutations are primarily small insertions and deletions, with deletions predominating. The divergent domains appear to be a good source of phylogenetic information for evolutionary events occurring approximately 100-200 million years ago.      

Kappa-chain constant-region gene sequences in genus Rattus: coding regions are diverging more rapidly than noncoding regions

We have determined the nucleotide sequence of a 1,200-base pair (bp) genomic fragment that includes the kappa-chain constant-region gene (C kappa) from two species of native Australian rodents, Rattus leucopus cooktownensis and Rattus colletti. Comparison of these sequences with each other and with other rodent C kappa genes shows three surprising features. First, the coding regions are diverging at a rate severalfold higher than that of the nearby noncoding regions. Second, replacement changes within the coding region are accumulating at a rate at least as great as that of silent changes. Third, most of the amino acid replacements are localized in one region of the C kappa domain--namely, the carboxy-terminal "bends" in the alpha-carbon backbone. These three features have previously been described from comparisons of the two allelic forms of C kappa genes in R. norvegicus. These data imply the existence of considerable evolutionary constraints on the noncoding regions (based on as yet undetermined functions) or powerful positive selection to diversify a portion of the constant-region domain (whose physiological significance is not known). These surprising features of C kappa evolution appear to be characteristic only of closely related C kappa genes, since comparison of rodent with human sequences shows the expected greater conservation of coding regions, as well as a predominance of silent nucleotide substitutions within the coding regions.      

The complete pattern of mutagenesis arising from the repair of site-specific psoralen crosslinks: analysis by oligonucleotide hybridization

Psoralen crosslinks were site-specifically placed in plasmid pBR322 near the BamHI site in the tet gene by enzymatically inserting mercurated nucleotides and reacting at the target site with a sulfhydryl-containing psoralen. The damaged plasmid was repaired in SOS-induced E. coli cells. Mutants were detected by colony hybridization to oligonucleotides in the target region, and their sequences were determined. The mutations are all base substitutions, 80% transitions and 20% transversions, similar to the mutations previously identified by the loss of tetracycline resistance. However, the mutation sites detected by a physical method, unconstrained by phenotypic changes, follow a broader distribution than those identified genetically. They occur primarily at favored psoralen crosslinking sites, where T-T and T-C interstrand crosslinks can be formed. A majority of these mutations are silent.     

Molecular Analysis of the Hot Spot Region Related to Length Mutations in Wheat Chloroplast Dnas. I. Nucleotide Divergence of Genes and Intergenic Spacer Regions Located in the Hot Spot Region

The nucleotide divergence of chloroplast DNAs around the hot spot region related to length mutation in Triticum (wheat) and Aegilops was analyzed. DNA sequences (ca. 4.5 kbp) of three chloroplast genome types of wheat complex were compared with one another and with the corresponding region of other grasses. The sequences region contained rbcL and psaI, two open reading frames, and a pseudogene, rpl23' (pseudogene for ribosomal protein L23) disrupted by AT-rich intergic spacer regions. The evolution of these genes in the closely related wheat complex is characterized by nonbiased nucleotide substitutions in terms of being synonymous/nonsynonymous, having A-T pressure transitions over transversions, and frequent changes at the third codon position, in contrast with the gene evolution among more distant plant groups where biased nucleotide substitutions have frequently occurred. The sequences of these genes had diverged almost in proportion to taxonomic distance. The sequence of the pseudogene rpl23' changed approximately two times faster than that of the coding region. Sequence comparison between the pseudogene and its protein-coding counterpart revealed different degrees of nucleotide homology in wheat, rice and maize, suggesting that the transposition timing of the pseudogene differed and/or that different rates of gene conversion operated on the pseudogene in the cpDNA of the three plant groups in Gramineae. The intergenic spacer regions diverged approximately ten times faster than the genes. The divergence of wheat from barley, and that from rice are estimated based on the nucleotide similarity to be 1.5, 10 and 40 million years, respectively.     

The 5S RNA genes of Schizosaccharomyces pombe.

The genomic arrangement and sequences of S. pombe 5S RNA genes are reported here. The 5S gene sequences appear to be dispersed within the genome, and are found independently of other rRNA genes. The sequences of two 5S genes examined show identical coding regions of 119 base pairs but have widely varying flanking sequences. A tRNAAsp gene is found in the 3' flanking region of one of the 5S genes. The tRNAAsp gene is faithfully transcribed in an X. laevis in vitro system, while the 5S genes are not transcribed in this system. The phylogenetic position of S. pombe is examined through comparison of 5S RNA sequences.     

Spontaneous Mutation in the Escherichia Coli Laci Gene

To gain more detailed insight into the nature and mechanisms of spontaneous mutations, we undertook a DNA sequence analysis of a large collection of spontaneous mutations in the N-terminal region of the Escherichia coli lacI gene. This region of circa 210 base pairs is the target for dominant lacI mutations (i-d) and is suitable for studies of mutational specificity since it contains a relatively high density of detectable mutable sites. Among 414 independent i-d mutants, 70.8% were base substitutions, 17.2% deletions, 7.7% additions and 4.3% single-base frameshifts. The base substitutions were both transitions (60%) and transversions (40%), the largest single group being G.C----A.T (47% of base substitutions). All four transversions were observed. Among the 71 deletions, a hotspot (37 mutants) was present: an 87-bp deletion presumably directed by an 8-bp repeated sequence at its endpoints. The remaining 34 deletions were distributed among 29 different mutations, either flanked (13/34) or not flanked (21/34) by repeated sequences. The 32 additions comprised 29 different events, with only two containing a direct repeat at the endpoints. The single-base frameshifts were the loss of a single base from either repeated (67%) or nonrepeated (33%) bases. A comparison with the spectrum obtained previously in strains defective in DNA mismatch correction (mutH, mutL, mutS strains) yielded information about the apparent efficiency of mismatch repair. The overall effect was 260-fold but varied substantially among different classes of mutations. An interesting asymmetry was uncovered for the two types of transitions, A.T----G.C and G.C----A.T being reduced by mismatch repair 1340- and 190-fold, respectively. Explanations for this asymmetry and its possible implications for the origins of spontaneous mutations are discussed.     

Purifying and directional selection in overlapping prokaryotic genes

In overlapping genes, the same DNA sequence codes for two proteins using different reading frames. Analysis of overlapping genes can help in understanding the mode of evolution of a coding region from noncoding DNA. We identified 71 pairs of convergent genes, with overlapping 3' ends longer than 15 nucleotides, that are conserved in at least two prokaryotic genomes. Among the overlap regions, we observed a statistically significant bias towards the 123:132 phase (i.e. the second codon base in one gene facing the degenerate third position in the second gene). This phase ensures the least mutual constraint on nonconservative amino acid replacements in both overlapping coding sequences. The excess of this phase is compatible with directional (positive) selection acting on the overlapping coding regions. This could be a general evolutionary mode for genes emerging from noncoding sequences, in which the protein sequence has not been subject to selection.     

Mutations in the NAM2 genes of Saccharomyces cerevisiae and S. douglasii are clustered non-randomly as a result of constraints on the nucleic acid sequence and not on the protein

Summary We present a statistical study of the nature and distribution of mutations along the NAM2 gene coding for the mitochondrial leucyl tRNA synthetase in Saccharomyces cerevisiae and S. douglasii (Herbert et al. 1988). Two important facts are observed: (1) the relative frequency of transitions and transversions is the same among silent substitutions and replacements. (2) The two kinds of mutations (silent substitutions and replacements) are distributed in the same way along the gene. This distribution is not random; the mutations are clustered and the clusters are regularly spaced along the gene. The NAM2 gene offers an example spaced along the gene. The NAM2 gene offers an example of recent divergence. We show that, in this case, the fixation of mutations is the result of genetic drift and of constraints on the nucleic acid sequence and not on that of the protein.     

Transitions,transversions, and the molecular evolutionary clock

Summary Nucleotide substitutions in the form of transitions (purine-purine or pyrimidine-pyrimidine interchanges) and transversions (purine-pyrimidine interchanges) occur during evolution and may be complied by aligning the sequences of homologous genes. Referring to the genetic code tables, silent transitions take place in third positions of codons in family boxes and two-codon sets. Silent transversions in third positions occur only in family boxes, except for AC transversions between AGR and CGR arginine codons (R=A or G). Comparisons of several protein genes have been made, and various subclasses of transitional and transversional nucleotide substitutions have been compiled. Considerable variations occur among the relative proportions of transitions and transversions. Such variations could possibly be caused by mutator genes, favoring either transitions or, conversely, transversions, during DNA replication. At earlier stages of evolutionary divergence, transitions are usually more frequent, but there are exceptions. No indication was found that transversions usually originate from multiple substitutions in transitions.     

Mutations induced by ultraviolet light

The different ultraviolet (UV) wavelength components, UVA (320-400 nm), UVB (280-320 nm), and UVC (200-280 nm), have distinct mutagenic properties. A hallmark of UVC and UVB mutagenesis is the high frequency of transition mutations at dipyrimidine sequences containing cytosine. In human skin cancers, about 35% of all mutations in the p53 gene are transitions at dipyrimidines within the sequence 5'-TCG and 5'-CCG, and these are localized at several mutational hotspots. Since 5'-CG sequences are methylated along the p53 coding sequence in human cells, these mutations may be derived from sunlight-induced pyrimidine dimers forming at sequences that contain 5-methylcytosine. Cyclobutane pyrimidine dimers (CPDs) form preferentially at dipyrimidines containing 5-methylcytosine when cells are irradiated with UVB or sunlight. In order to define the contribution of 5-methylcytosine to sunlight-induced mutations, the lacI and cII transgenes in mouse fibroblasts were used as mutational targets. After 254 nm UVC irradiation, only 6-9% of the base substitutions were at dipyrimidines containing 5-methylcytosine. However, 24-32% of the solar light-induced mutations were at dipyrimidines that contain 5-methylcytosine and most of these mutations were transitions. Thus, CPDs forming preferentially at dipyrimidines with 5-methylcytosine are responsible for a considerable fraction of the mutations induced by sunlight in mammalian cells. Using mouse cell lines harboring photoproduct-specific photolyases and mutational reporter genes, we showed that CPDs (rather than 6-4 photoproducts or other lesions) are responsible for the great majority of UVB-induced mutations. An important component of UVB mutagenesis is the deamination of cytosine and 5-methylcytosine within CPDs. The mutational specificity of long-wave UVA (340-400 nm) is distinct from that of the shorter wavelength UV and is characterized mainly by G to T transversions presumably arising through mechanisms involving oxidized DNA bases. We also discuss the role of DNA damage-tolerant DNA polymerases in UV lesion bypass and mutagenesis.     

SNP-based markers for discriminating olive (Olea europaea L.) cultivars.

A set of 11 polymorphic markers (1 cleaved amplified polymorphic sequence (CAPS), 2 sequence-characterized amplified regions (SCARs), and 8 single-nucleotide polymorphism (SNP)-derived markers) was obtained for olive cultivar identification by comparing DNA sequences from different accessions. Marker development was more efficient, using sequences from the database rather than cloning arbitrary DNA fragments. Analyses of the sequences of 3 genes from 11 diverse cultivars revealed an SNP frequency of 1 per 190 base pairs in exons and 1 per 149 base pairs in introns. Most mutations were silent or had little perceptible effect on the polypeptide encoded. The higher incidence of transversions (55%) suggests that methylation is not the major driving force for DNA base changes. Evidence of linkage disequilibrium in 2 pairs of markers has been detected. The set of predominantly SNP-based markers was used to genotype 65 olive samples obtained from Europe and Australia, and was able clearly to discriminate 77% of the cultivars. Samples, putatively of the same cultivar but derived from different sources, were revealed as identical, demonstrating the utility of these markers as tools for resolving nomenclature issues. Genotyping data were used for constructing a dendrogram by UPGMA cluster analysis using the simple matching similarity coefficient. Relationships between cultivars are discussed in relation to the route of olive's spread.