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.     

VH genes encoding the immune response to beta-(1,6)-galactan: somatic mutation in IgM molecules

The immune response to beta-(1,6)-galactan in the BALB/c mouse has been well characterized and includes the amino acid sequence determination of 13 monoclonal antibodies. The genetic potential encoding the VH regions of these antibodies has been determined by isolation and sequencing of homologous germline genes. The germline repertoire encoding these proteins was found to consist of two closely related genes. One of these directly encodes the VH segments of seven Gal-binding proteins, and the second directly encodes one additional protein sequence. Sequence variations found in the VH regions of five other Gal-binding proteins can be explained by somatic mutations leading to single base substitutions in the more frequently used gene. Since four of the hybridoma proteins exhibiting somatic mutations are of the IgM class, these results indicate that somatic mutation, in this system, is not associated with class switching and can apparently be initiated early in B-cell development. The two Gal genes are the only members of a very restricted multigene family and probably result from a gene duplication estimated to occur 1.4-2.8 million years ago. Three other genes hybridizing at moderate stringency to a VHGal probe were also sequenced and were found to be members of two additional VHIII families. Studies of the silent to replacement substitution ratios of these and other VH genes indicate that the number of silent substitutions found in immunoglobulin VH genes is lower than expected when compared with proteins such as preproinsulin and globin. Analysis of base composition reflected in these sequences indicates a marked increase of A-T% in the first and second codon positions of complementarity determining regions (CDR) which may be important in facilitating point mutations.     

Intra- and interbreed genetic variations of mitochondrial DNA major non-coding regions in Japanese native dog breeds [Canis familiaris)

Mitochondrial DNA (mtDNA) major non-coding regions were amplified from 73 dogs of eight Japanese native dog breeds and from 21 dogs of 16 non-Japanese dog breeds by the polymerase chain reaction and their DNA sequences were determined. A total of 51 nucleotide positions within the non-coding region (969–972 base pairs) showed nucleotide variations of which 48 were caused by transition. These nucleotide substitutions were abundant in the region proximate to tRNA^Pro. In addition to the nucleotide substitutions, the dog mtDNA D-loop sequences had a heteroplasmic repetitive sequence (TACACGTÀCG) involving size variation. The DNA sequences of the non-coding region were classified into four different groups by phylogenetic analysis and the deepest branchpoints of this dog phylogeny was calculated to about 100 000 years before the present. Phylogenetic analysis showed that Japanese native dog breeds could not be clearly delimited as distinct breeds. Many haplotypes found in members of some clustering groups were seen in each dog breed, and interbreed nucleotide differences between Japanese dog breeds were almost the same as the intrabreed nucleotide diversities.     

The silkmoth late chorion locus. I. Variation within two paired multigene families

The 140 X 10(3) base late chorion locus of Bombyx mori contains two 15-member multigene families arranged in tightly linked pairs, which are divergently transcribed (the high-cysteine A (HcA) and the high-cysteine B (HcB) families). Previous DNA hybridization experiments have indicated that all members of these gene families contain a complex pattern of shared sequence variation. The sequence analysis in this paper involving all 15 gene pairs allows a comprehensive examination of the nature of this variation. Average sequence homology between gene pairs is: 95% for the protein-encoding regions; 93% for the common 272 base-pair 5' flanking region; 87% for the introns; and 88% for the 3' untranslated regions. Considering the great degree of sequence homology in the coding regions, an unexpectedly high level of variation is found in the deduced protein sequences. Over 50% of the nucleotide substitutions in the protein-encoding regions lead to amino acid replacements, most of which involve a change in charge or effect the secondary structure of the protein. In addition, significant differences in length between the proteins occur in the carboxyl-terminal arm. In both families, the major portion of this arm is composed of Cys-Gly-Gly and Cys-Gly subrepeats forming a (Cys-Gly-Gly)2-(Cys-Gly)2 major repeat. Differences in the number of complete and partial repeats results in deduced protein sequences that contain arms varying from 32 to 54 amino acid residues for members of the HcA family and 14 to 88 residues for the HcB family. The high level of variation in protein composition indicates a lack of strong selective pressure. We suggest the high level of DNA sequence homology maintained by these genes in the coding as well as in the non-coding regions is the result of sequence exchange between family members.     

Pattern and rate of indel evolution inferred from whole chloroplast intergenic regions in sugarcane, maize and rice.

Microstructural changes such as insertions and deletions (=indels) are a major driving force in the evolution of non-coding DNA sequences. To better understand the mechanisms by which indel mutations arise, as well as the molecular evolution of non-coding regions, the number and pattern of indels and nucleotide substitutions were compared in the whole chloroplast genomes. Comparisons were made for a total of over 38 kb non-coding DNA sequences from 126 intergenic regions in two data sets representing species with different divergence times: sugarcane and maize and Oryza sativa var. indica and japonica. The main findings of this study are: (i) Approximately half of all indels are single nucleotide indels. This observation agrees with previous studies in various organisms. (ii) The distribution and number of indels was different between two data sets, and different patterns were observed for tandem repeat and non-repeat indels. (iii) Distribution pattern of tandem repeat indels showed statistically significant bias towards A/T-rich. (iv) The rate of indel mutation was estimated to be approximately 0.8 +/- 0.04 x 10(-9) per site per year, which was similar to previous estimates in other organisms. (v) The frequencies of nucleotide substitutions and indels were significantly lower in inverted repeat (IR).     

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.     

Identification of 142 single nucleotide polymorphisms in 41 candidate genes for rheumatoid arthritis in the Japanese population

Single-nucleotide polymorphisms (SNPs) can make an important contribution to our understanding of genetic backgrounds that may influence medical conditions and ethnic diversity. We undertook a systematic survey of genomic DNA for SNPs located not only in coding sequences but also in non-coding regions (e.g., introns and 5' flanking regions) of selected genes. Using DNA samples from 48 Japanese patients with rheumatoid arthritis (RA) as templates, we surveyed 41 genes that represent candidates for RA, screening a total of 104 kb of DNA (30 kb of coding sequences and 74 kb of non-coding DNA). Within this 104 kb of genomic sequences we identified 163 polymorphisms (1 per 638 bases on average), of which 142 were single-nucleotide substitutions and the remainder, insertions or deletions. Of the coding SNPs, 52% were non-synonymous substitutions, and non-conservative amino acid changes were observed in a quarter of those. Sixty-nine polymorphisms showed high frequencies for minor alleles (more than 15%) and 20 revealed low frequencies (<5%). Our results indicated a greater average distance between SNPs than others have reported, but this disparity may reflect the type of genes surveyed and/or the relative ethnic homogeneity of our test population.     

Genetic relationship amongst the major non-coding regions of mitochondrial DNAs in wild boars and several breeds of domesticated pigs

We completed phylogenetic analysis of the major non-coding region of the mitochondrial DNA (mtDNA) from 159 animals of eight Euro-American and six East Asian domesticated pig breeds and 164 Japanese and five European wild boars. A total of 62 mtDNA haplotypes were detected. Alignment of these regions revealed nucleotide variations (including gaps) at 73 positions, including 58 sites with transition nucleotide substitutions, and two transversion substitutions. Phylogenetic analysis of the sequences could not organize domestic pig breeds into discrete clusters. In addition, many of the haplotypes found in members of diverged clustering groups were found primarily in Euro-American pig breeds, indicating extensive introgression of Asian domestic pigs into European breeds. Furthermore, phylogenetic analysis allocated the DNA sequences of non-coding regions into two different groups, and the deepest branchpoint of this porcine phylogeny corresponded to 86 000-136 000 years before present. This time of divergence would predate the historical period when the pig is thought to have been domesticated from the wild boar.     

Unusual mutation clusters provide insight into class I gene conversion mechanisms.

Genetic diversity among the K and D alleles of the mouse major histocompatibility complex is generated by gene conversion among members of the class I multigene family. The majority of known class I mutants contain clusters of nucleotide changes that can be traced to linked family members. However, the details of the gene conversion mechanism are not known. The bm3 and bm23 mutations represent exceptions to the usual pattern and provide insight into intermediates generated during the gene conversion process. Both of these variants contain clusters of five nucleotide substitutions, but they differ from the classic conversion mutants in the important respect that no donor gene for either mutation could be identified in the parental genome. Nevertheless, both mutation clusters are composed of individual mutations that do exist within the parent. Therefore, they are not random and appear to be templated. Significantly, the bm3 and bm23 mutation clusters are divided into overlapping regions that match class I genes which have functioned as donor genes in other characterized gene conversion events. The unusual structure of the mutation clusters indicates an underlying gene conversion mechanism that can generate mutation clusters as a result of the interaction of three genes in a single genetic event. The unusual mutation clusters are consistent with a hypothetical gene conversion model involving extrachromosomal intermediates.     

Characterization of ω-secalin genes from rye, triticale, and a wheat 1BL/1RS translocation line

Sixty-two DNA sequences for the coding regions of omega-secalin (ω-secalin) genes have been characterized from rye (Secale cereale L.), hexaploid and octoploid triticale (× Triticosecale Wittmack), and wheat (Triticum aestivum L.) 1BL/1RS translocation line. Only 19 out of the 62 ω-secalin gene sequences were full-length open reading frames (ORFs), which can be expressed into functional proteins. The other 43 DNA sequences were pseudogenes, as their ORFs were interrupted by one or a few stop codons or frameshift mutations. The 19 ω-secalin genes have a typical primary structure, which is different from wheat gliadins. There was no cysteine residue in ω-secalin proteins, and the potential celiac disease (CD) toxic epitope (PQQP) was identified to appear frequently in the repetitive domains. The ω-secalin genes from various cereal species shared high homology in their gene sequences. The ω-secalin gene family has involved fewer variations after the integration of the rye R chromosome or whole genome into the wheat or triticale genome. The higher Ka/Ks ratio (i.e. non-synonymous to synonymous substitutions per site) in ω-secalin pseudogenes than in ω-secalin ORFs indicate that the pseudogenes may be subject to a reduced selection pressure. Based on the conserved sequences of ω-secalin genes, it will be possible to manipulate the expression of this gene family in rye, triticale, or wheat 1BL/1RS translocation lines, to reduce its negative effects on grain quality.     

Patterns of Somatic Mutations in Immunoglobulin Variable Genes

The mechanism responsible for somatic mutation in the variable genes of antibodies is unknown and may differ from previously described mechanisms that produce mutation in DNA. We have analyzed 421 somatic mutations from the rearranged immunoglobulin variable genes of mice to determine if the nucleotide substitutions differ from those generated during meiosis and if the presence of nearby direct and inverted repeated sequences could template mutations around the variable gene. The results reveal a difference in the pattern of substitutions obtained from somatic mutations vs. meiotic mutations. An increased frequency of T:A to C:G transitions and a decreased frequency of mutations involving a G in the somatic mutants compared to the meiotic mutants is indicated. This suggests that the mutational processes responsible for somatic mutations in antibody genes differs from that responsible for mutation during meiosis. An analysis of the local DNA sequences revealed many direct repeats and palindromic sequences that were capable of templating some of the known mutations. Although additional factors may be involved in targeting mutations to the variable gene, mistemplating by nearby repeats may provide a mechanism for the enhancement of somatic mutation.     

CG dinucleotide transitions in the factor IX gene account for about half of the point mutations in hemophilia B patients: a Seattle series

Summary Hemophilia B is due to multiple molecular defects in the factor IX gene. Over 80% of mutations are single base substitutions. By amplification and direct sequencing, 51 single base substitutions were found in the transcribed sequence of the factor IX genes of patients from 50 distinct families with hemophilia B. These include 30 mutations in 29 families not previously reported by us; of these, 12 are novel, i.e., not previously published in other series. Of the 51 substitutions in our overall series 23 (45%) occurred as C-to-T or G-to-A transitions at 11 sites within CG dinucleotides. It is estimated that CG transitions occur from one to two orders of magnitude more frequently than mutations in nucleotides that are not within a CG pair. More than one family had identical defects for 6 of the CG mutations. At 4 of these sites, most patients had different haplotypes compatible with distinct mutations. Non-CG-type mutations occurred thoughout the coding regions with only one mutation in more than one family. The latter included 7 families with a 397 Ile-to-Thr defect that all share a rare haplotype, suggesting a common ancestor.     

Intra- and interspecies analyses of the carcinoembryonic antigen (CEA) gene family reveal independent evolution in primates and rodents

Summary Various rodent and primate DNAs exhibit a stronger intra- than interspecies cross-hybridization with probes derived from the N-terminal domain exons of human and rat carcinoembryonic antigen (CEA)-like genes. Southern analyses also reveal that the human and rat CEA gene families are of similar complexity. We counted at least 10 different genes per human haploid genome. In the rat, approximately seven to nine different N-terminal domain exons that presumably represent different genes appear to be present. We were able to assign the corresponding genomic restriction endonuclease fragments to already isolated CEA gene family members of both human and rat. Highly similar subgroups, as found within the human CEA gene family, seem to be absent from the rat genome. Hybridization with an intron probe from the human nonspecific cross-reacting antigen (NCA) gene and analysis of DNA sequence data indicate the conservation of noncoding regions among CEA-like genes within primates, implicating that whole gene units may have been duplicated. With the help of a computer program and by calculating the rate of synonymous substitutions, evolutionary trees have been derived. From this, we propose that an independent parallel evolution, leading to different CEA gene families, must have taken place in, at least, the primate and rodent orders.     

Discrete size classes of monkey extrachromosomal circular DNA containing the L1 family of long interspersed nucleotide sequences are produced by a general non-sequence specific mechanism

The L1 family of long interspersed nucleotide sequences (LINES) has recently been identified and characterized in the small polydisperse circular DNA (spc-DNA) populations of monkey (1), human (2) and mouse (3) cells. In monkey spc-DNA, the L1 (also known as Kpn I) family is present in discrete size classes (ranging from 300 to 6000 base pairs (bp)) which appear to be generated by non homologous recombination events within chromosomal elements. In this communication it is shown that different regions of the consensus L1 family are present at different frequencies in monkey spc-DNA (as they are in chromosomal DNA), that all regions of the family are present in extrachromosomal DNA, and that each region appears to be represented in an identical discrete spc-DNA size distribution. This size distribution reflects a non-sequence specific mechanism that generates spc-DNA size classes by chromosomal DNA recombination events that are in some way constrained to occur between sites separated by relatively defined lengths.     

A new family of interspersed repetitive DNA sequences in the mouse genome

Repetitive DNA sequences near immunoglobulin genes in the mouse genome (Steinmetz et al., 1980a,b) were characterized by restriction mapping and hybridization. Six sequences were determined that turned out to belong to a new family of dispersed repetitive DNA. From the sequences, which are called R1 to R6, a 475 base-pair consensus sequence was derived. The R family is clearly distinct from the mouse B1 family (Krayev et al., 1980). According to saturation hybridization experiments, there are about 100,000 R sequences per haploid genome, and they are probably distributed throughout the genome. The individual R sequences have an average divergence from the consensus sequence of 12.5%, which is largely due to point mutations and, among those, to transitions. Some R sequences are severly truncated. The R sequences extend into A-rich sequences and are flanked by short direct repeats. Also, two large insertions in the R2 sequence are flanked by direct repeats. In the neighbourhood of and within R sequences, stretches of DNA have been identified that are homologous to parts of small nuclear RNA sequences. Mouse satellite DNA-like sequences and members of the B1 family were also found in close proximity to the R sequences. The dispersion of R sequences within the mouse genome may be a consequence of transposition events. The possible role of the R sequences in recombination and/or gene conversion processes is discussed.     

The loss of dinucleotides CpG from DNA. IV. Methylation and divergence of genes and pseudogenes of small nuclear RNA

The frequency of neighboring base pairs in nucleotide sequences of over 80 genes and pseudogenes of low molecular weight RNAs U1-U8, 4.5S and 7S in different eukaryotes was determined. The probable frequency of CpG----TpG + CpA substitutions, caused as a result of the deamination of 5-methylcytosine residues in DNA, was determined. It was found that the genes of small RNAs do not reveal a single level of CpG methylation for all the species studied. In most cases CpG in the genes of U1, 4.5S and 7S RNA are methylated, whereas in the genes of U2-U6-RNA these sites must have never been subjected to methylation. Nearly all the investigated pseudogenes of different small RNAs are strongly methylated due to a considerable lack of CpG. It was established that CpG----TpG + CpA transitions may amount to as much third of all the mutations accumulated in the genes of the same RNAs in different species. Such transitions in pseudogenes may account for 40% of all the nucleotide substitutions. This disproportionately high level of mutations in CpG dinucleotides (3-5-fold higher than in other DNA dupletes) must be the direct result of the methylation of these sites. Consequently, CpG methylation causes a dramatic acceleration of the divergence rate of DNA sequences. It has been concluded that protection of most vital genes against methylation is one of the essential conditions for sustaining the high level of stability of the macromolecular structure and for the reliability of macromolecular functioning in a cell.     

Using Triplet Periodicity of Nucleotide Sequences for Finding Potential Reading Frame Shifts in Genes

We introduce a novel approach for the detection of possible mutations leading to a reading frame (RF) shift in a gene. Deletions and insertions of DNA coding regions are considerable events for genes because an RF shift results in modifications of the extensive region of amino acid sequence coded by a gene. The suggested method is based on the phenomenon of triplet periodicity (TP) in coding regions of genes and its relative resistance to substitutions in DNA sequence. We attempted to extend 326 933 regions of continuous TP found in genes from the KEGG databank by considering possible insertions and deletions. We revealed totally 824 genes where such extension was possible and statistically significant. Then we generated amino acid sequences according to active (KEGG''s) and hypothetically ancient RFs in order to find confirmation of a shift at a protein level. Consequently, 64 sequences have protein similarities only for ancient RF, 176 only for active RF, 3 for both and 581 have no protein similarity at all. We aimed to have revealed lower bound for the number of genes in which a shift between RF and TP is possible. Further ways to increase the number of revealed RF shifts are discussed.     

Evidence for local DNA influences on patterns of substitutions in the human alpha-interferon gene family

The evolutionary history of genes can be used to examine patterns of spontaneous mutation if the sequences are sufficiently extensive to provide reliable data. Many human alpha-interferon genes have been sequenced and they form a large multigene family including several pseudogenes. A phylogenetic history for 15 human interferon sequences was reconstructed and their ancestral sequences inferred using a maximum parsimony method. This evolutionary history provided a record of more than 738 spontaneous mutations that have occurred in man's recent evolution. Of these mutations, more than 267 base substitution and deletion-insertion events were analyzed to determine the possible effects of nearby DNA sequences. Many substitutions occur at the end of long runs of identical bases and some dinucleotide pairs may mutate more often than others. Because templating by local DNA sequences has been implicated in prokaryotic mutation, the sequences were also examined for nearby repeats that include the substituted nucleotide and hence are potentially capable of templating the substitution. The majority of sequence alterations examined have either a similar direct repeat or palindrome nearby. Often such templates can account for simultaneous multiple mutations. These results suggest that sequence-directed events may occur occasionally in eukaryotes and that neighbouring DNA sequences can influence both the occurrence and types of mutations in several different ways.