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.     

Characterization of a mouse interferon gene locus II. Differential expression of alpha-interferon genes.

A cluster of four MuIFN-alpha genes was recently isolated and characterized (1); one of the genes in this cluster had, in the coding region, an internal deletion of 5 amino acids. Bacterial expression plasmids were constructed to examine the effect of this deletion on the antiviral activity of the MuIFN-alpha 4 peptide and it was found that the alpha 4 interferon peptide had a 100-fold lower antiviral activity than full length alpha-interferon proteins when expressed in E. coli. Three of the four MuIFN-alpha genes identified were expressed coordinately in L-cells infected with NDV. The relative levels of alpha 4 mRNA were substantially higher than the levels of the other alpha mRNAs. Comparison of the 5' end flanking sequences of these four alpha interferon genes revealed that the promoter sequences of alpha 1, alpha 5 and alpha 6 are more homologous to each other than to the alpha 4 promoter which also contains a G rich cluster not seen in the other three promoters.     

Identification and mapping on human chromosome 19 of chromosome-specific repeated element preferentially binding with nuclear matrix]

From a library of sequences binding preferentially to nuclear matrix (matrix attachment regions, MARs), a fragment of about 300 bp in length (CEA (carcinoembryonic antigen)-MAR) was isolated and characterized. The CEA-MAR sequence was found in more than ten loci of chromosome 19 containing elements similar to genes of the CEA family. No sequences of this group were found on other human chromosomes. Two CEA-MAR-containing loci were sequenced, and sequences for another seven loci were found in GenBank. A comparative analysis of CEA-MARs and the flanking sequences is reported. Based on the sequence of the CEA-containing chromosome 19 loci, a hypothetical model of the domain structure of a 2-Mb chromosome region was constructed and the mutual arrangement of CEA-MARs and genes of CEA family was elucidated. The CEA-MARs were located 5-20 kb downstream of the CEA genes. These results suggest that the duplication unit of the CEA family may coincide with chromatin domains containing these genes.     

A unique structure of a mouse gamma-actin processed pseudogene

We have isolated several gamma-actin-related genes from a mouse genomic library. One of these has been shown to be a gamma-actin processed pseudogene (Tokunga, K., Yoda, K. and Sakiyama, S. (1985) Nucleic Acids Res. 13, 3031-3042). Here, we report the structure of another pseudogene (pMA131). pMA131 contained the sequences corresponding to the carboxyl half of a cytoskeletal actin in which random point mutations as well as insertion and deletion events took place. This region was flanked at its 5' end by the sequences related to mouse repetitive sequences, including the MIF-1 family, and was interrupted by the sequence homologous to the R family which is also a mouse repetitive sequence. The coding region was followed by the sequence corresponding to 3' untranslated region of gamma-actin mRNA.     

One of the copia genes is adjacent to satellite DNA in Drosophila melanogaster.

A method for purifying sequences adjacent to satellite DNA in the heterochromatin of D. melanogaster is described. A cloned DNA segment containing part of a copia gene adjacent to 1.688 g/cm3 satellite DNA has been isolated. The copia genes compose a repeated gene family which codes for abundant cytoplasmic poly(a)-containing RNA (Young and Hogness, 1977; Finnegan et al., 1978). We have identified two major poly (A)-containing RNA species [5.2 and 2.1 kilobases (kb)] produced by the copia gene family. The cloned segment contains copia sequences homologous to the 5' end of RNA within 0.65 kb of the 1.688 satellite DNA sequences. Seven different cloned copia genes from elsewhere in the genome have also been isolated, and a 5.2 kb region present in five of the clones was identified as copia by heteroduplex analysis. In addition, three ususual copies of copia were found: a "partial" copy of the gene (3.7 kb) which has one endpoint in common with the 5.2 kb unit; a copia gene flanked on one side by a 1.6 kb sequence and on the other by the same 1.6 kb sequence in the inverted orientation; and a copia gene flanked only on one side by the same sequence.     

PCR-based isolation and identification of full-length low-molecular-weight glutenin subunit genes in bread wheat (<Emphasis Type="Italic">Triticum aestivum</Emphasis> L.)

Low-molecular-weight glutenin subunits (LMW-GSs) are encoded by a multi-gene family and are essential for determining the quality of wheat flour products, such as bread and noodles. However, the exact role or contribution of individual LMW-GS genes to wheat quality remains unclear. This is, at least in part, due to the difficulty in characterizing complete sequences of all LMW-GS gene family members in bread wheat. To identify full-length LMW-GS genes, a polymerase chain reaction (PCR)-based method was established, consisting of newly designed conserved primers and the previously developed LMW-GS gene molecular marker system. Using the PCR-based method, 17 LMW-GS genes were identified and characterized in Xiaoyan 54, of which 12 contained full-length sequences. Sequence alignments showed that 13 LMW-GS genes were identical to those found in Xiaoyan 54 using the genomic DNA library screening, and the other four full-length LMW-GS genes were first isolated from Xiaoyan 54. In Chinese Spring, 16 unique LMW-GS genes were isolated, and 13 of them contained full-length coding sequences. Additionally, 16 and 17 LMW-GS genes in Dongnong 101 and Lvhan 328 (chosen from the micro-core collections of Chinese germplasm), respectively, were also identified. Sequence alignments revealed that at least 15 LMW-GS genes were common in the four wheat varieties, and allelic variants of each gene shared high sequence identities (>95%) but exhibited length polymorphism in repetitive regions. This study provides a PCR-based method for efficiently identifying LMW-GS genes in bread wheat, which will improve the characterization of complex members of the LMW-GS gene family and facilitate the understanding of their contributions to wheat quality.     

Nucleotide and deduced amino acid sequences of mutanase-like genes from Paenibacillus isolates: proposal of a new family of glycoside hydrolases

Three mutanase (alpha-1,3-glucanase)-producing microorganisms isolated from soil samples were identified as a relatives of Paenibacillus. A mutanase was purified to homogeneity from cultures of each, and the molecular masses of the purified enzymes were approximately 132, 141, and 141kDa, respectively. The corresponding three genes for mutanases were cloned by PCR using primers designed from each N-terminal amino acid sequence. Another mutanase-like gene from one strain was also cloned by PCR using primers designed from conserved amino acid sequences among known mutanases. Consequently, four mutanase-like genes were sequenced. The genes contained long open reading frames of 3411 to 3915bp encoding 1136 to 1304 amino acids. The deduced amino acid sequences of the mutanases showed relatively high similarity to those of a mutanase (E16590) from Bacillus sp. RM1 with 46.9% to 73.2% identity and an alpha-1,3-glucanase (AB248056) from Bacillus circulans KA-304 with 46.7% to 70.4% identity. Phylogenetic analysis based on the amino acid sequences of the enzymes showed bacterial mutanases form a new family between fungal mutanases (GH family 71) and Streptomycetes mycodextranases (GH family 87).     

Sequence homologies in the mouse protamine 1 and 2 genes

To identify candidates for cis-acting sequences that regulate the stage and cell-specific expression of the two coordinately regulated protamine genes in the mouse, genomic clones were isolated and the nucleotide sequences of the 5′ flanking regions and coding regions were compared. Unlike most histone genes and the multigene family of trout protamine genes which are intronless, each mouse protamine gene has a single, short intervening sequence. Although the coding regions do not share significant nucleotide homology, the 5′ flanking regions contain several short homologous sequences that may be involved in gene regulation. An additional shared sequence is present in the 3′ untranslated region surrounding the poly(A) addition signal in both genes.     

Rat ribosomal protein L35a multigene family: molecular structure and characterization of three L35a-related pseudogenes

The rat ribosomal protein L35a gene comprises a multigene family which contains 15-20 members as shown by the Southern blot analysis using L35a cDNA as a probe. We isolated 15 independent clones which contained distinct genes from a rat genomic library. Analysis of the restriction sites showed that all of them lacked the intervening sequences. Thermal stability of the hybrid molecules between these genes and the cDNA indicated that the similarity of the genes to the cDNA sequence varied. The nucleotide sequences of three genes gRL35a-A, gRL35a-B and gRL35a-G were determined. They shared some characteristics; namely: they lacked the intervening sequences, they contained (A)-rich tracts, and they were flanked by direct repeats. Two genes, gRL35a-A and gRL35a-B, contained a sequence completely identical to that of the cDNA. The nucleotide sequence of the 5' flanking region of gRL35a-B showed a significant homology with that of the same region of mouse ribosomal protein L32-related unmutated processed genes. Although this region of gRL35a-B contained the sequences homologous to the TATA box and the CCAAT box, gRL35a-B was not transcribed in an in vitro assay system. Thus, the L35a gene family comprises mostly processed pseudogenes. Further, Southern blot analysis in various animals indicated that the multigene construction of this ribosomal protein gene was a feature of mammalian genes. The origin and the evolutionary aspect of processed pseudogenes are discussed.     

Human histone genes are interspersed with members of the Alu family and with other transcribed sequences

We have isolated a series of recombinant λCh4A phages containing human histone genes. Histone H2A, H2B, H3 and H4 genes have been found to be clustered, but are not present in any simple repeat pattern. Hybridization of a blot containing phage DNA with S phase polysomal cDNA indicates the presence of additional sequences complementary to HeLa polysomal RNA sequences. Northern blot analysis using these clones as probes has also shown the presence of sequences complementary to non-histone-coding RNAs, some of which accumulate differentially in different stages of the cell cycle. We have also found, by hybridization with appropriate probes, that histone genes are interspersed with several copies of the Alu DNA family; however, not all of the histone genes are associated with an Alu DNA sequence.     

Retropseudogenes for human chromosomal protein HMG-17

The human genome contains multiple copies of sequences homologous to the cDNA coding for non-histone chromosomal protein HMG-17. To study the mechanism of generation and dispersion of the HMG-17 multigene family a human genomic library was screened and 70 clones isolated and studied by Southern transfer and restriction site analysis. The results suggest that most of the clones contain unique sequences. Sequence analysis of two genomic clones indicates that they contain elements typical of processed retropseudogenes. Even though both sequences contained open reading frames the sequences lacked introns, were flanked by short, direct repeats and lacked elements associated with functional genes. The sequences of the two pseudogenes were 85% homologous to each other and each was 90% homologous to the human cDNA. Based on the sequence difference in the open reading frame between the pseudogenes and the cDNA it can be estimated that the sequences arose approximately ten million years ago from a common precursor. The present paper, which is the first study on genes coding for this nucleosomal binding protein, indicates that the HMG-17 multigene family is the largest known human retropseudogene family.     

Evolution of human cytoplasmic actin gene sequences: chromosome mapping and structural characterizations of three cytoplasmic actin-like pseudogenes including one on the Y chromosome

Eight recombinant phage clones containing cytoplasmic actin-like gene sequences have been isolated from a human genomic library for structural characterization. Kpn I family repeat sequences flank six of these actin genes isolated, and Alu family repeats are scattered throughout the DNA inserts of all eight phage clones. Three of these genes are γ actin-like, and the other five are β actin-like. The complete nucleotide sequence analysis of one β and one γ actin-like genes and their flanking regions demonstrates that they both are processed pseudogenes. Using unique DNA sequences flanking these two pseudogenes as hybridization probes for human-mouse somatic cell hybrid DNAs, we have mapped the two actin pseudogenes on human chromosomes 8 and 3, respectively. We have also determined the DNA sequence of a human Y chromosome-linked, processed actin pseudogene. The different values of sequence divergence of these processed pseudogenes and their functional counterparts allow us to estimate the time of generation of the pseudogenes. The results suggest that the cDNA insertion events generating the human cytoplasmic actin-like pseudogenes have occurred at significantly different times during the evolution of primates, after their separation from other mammalian species.     

Isolation and characterization of six different chicken actin genes.

Genes representing six different actin isoforms were isolated from a chicken genomic library. Cloned actin cDNAs as well as tissue-specific mRNAs enriched in different actin species were used as hybridization probes to group individual actin genomic clones by their relative thermal stability. Restriction maps showed that these actin genes were derived from separate and nonoverlapping regions of genomic DNA. Of the six isolated genes, five included sequences from both the 5' and 3' ends of the actin-coding area. Amino acid sequence analysis from both the NH2- and COOH-terminal regions provided for the unequivocal identification of these genes. The striated isoforms were represented by the isolated alpha-skeletal, alpha-cardiac, and alpha-smooth muscle actin genes. The nonmuscle isoforms included the beta-cytoplasmic actin gene and an actin gene fragment which lacked the 5' coding and flanking sequence; presumably, this region of DNA was removed from this gene during construction of the genomic library. Unexpectedly, a third nonmuscle chicken actin gene was found which resembled the amphibian type 5 actin isoform (J. Vandekerckhove, W. W. Franke, and K. Weber, J. Mol. Biol., 152:413-426). This nonmuscle actin type has not been previously detected in warm-blooded vertebrates. We showed that interspersed, repeated DNA sequences closely flanked the alpha-skeletal, alpha-cardiac, beta-, and type 5-like actin genes. The repeated DNA sequences which surround the alpha-skeletal actin-coding regions were not related to repetitious DNA located on the other actin genes. Analysis of genomic DNA blots showed that the chicken actin multigene family was represented by 8 to 10 separate coding loci. The six isolated actin genes corresponded to 7 of 11 genomic EcoRI fragments. Only the alpha-smooth muscle actin gene was shown to be split by an EcoRI site. Thus, in the chicken genome each actin isoform appeared to be encoded by a single gene.     

Human ventricular/slow twitch myosin alkali light chain gene characterization, sequence, and chromosomal location

The gene coding for the human ventricular/slow twitch myosin alkali light chain isoform was isolated and sequenced. It was found to contain a total of seven exons, the last of which is completely 3'-untranslated sequence. Comparison of this gene sequence with that of the various fast twitch skeletal isoform gene sequences revealed that the exon-intron arrangement is conserved within the myosin alkali light chain gene family. In fact the introns are in exactly the same positions within analogous codons. Comparison of the derived amino acid sequence from the human ventricular/slow twitch isoform gene with that of other isoform protein sequences indicated that the protein encoded by this gene is more homologous to the chicken cardiac isoform protein sequence than to any of the other protein sequences. These results indicate that the gene duplication which gave rise to the ventricular/slow twitch and fast twitch isoform genes must have occurred prior to the divergence of mammals and avians. We have also localized the human ventricular/slow twitch isoform gene to the short arm of human chromosome 3. Interestingly the corresponding mouse gene has been mapped to the distal region of mouse chromosome 9 which contains a conserved syntenic group of genes that map to the short arm of human chromosome 3.     

Nucleotide sequence analysis of the human salivary protein genes HIS1 and HIS2, and evolution of the STATH/HIS gene family

Human histatins are a family of low-M(r), neutral to very basic, histidine-rich salivary polypeptides. They probably function as part of the nonimmune host defense system in the oral cavity. A 39-kb region of DNA containing the HIS1 and HIS2 genes was isolated from two human genomic phage libraries as a series of overlapping clones. The nucleotide sequences of the HIS1 gene and part of the HIS2(1) gene were determined. The transcribed region of HIS1 spans 8.5 kb and contains six exons and five introns. The HIS1 and HIS2(1) genes exhibit 89% overall sequence identity, with exon sequences exhibiting 95% identity. The two loci probably arose by a gene duplication event approximately 15-30 Mya. The HIS1 sequence data were also compared with that of STATH. Human statherin is a low-M(r) acidic phosphoprotein that acts as an inhibitor of precipitation of calcium phosphate salts in the oral cavity. The HIS1 and STATH genes show nearly identical overall gene structures. The HIS1 and STATH loci exhibit 77%-81% sequence identity in intron DNA and 80%-88% sequence identity in noncoding exons but only 38%-43% sequence identity in the protein-coding regions of exons 4 and 5. These unusual data suggest that HIS1, HIS2, and STATH belong to a single gene family exhibiting accelerated evolution between the HIS and STATH coding sequences.