The rat casein multigene family. Fine structure and evolution of the beta-casein gene

Eight overlapping phage clones, spanning 34.4 kilobase pairs of genomic DNA, containing the 7.2-kilobase pair rat beta-casein gene have been isolated and characterized. The first 510 base pairs (bp) of 5' flanking, 110 bp of 3' flanking, and all the exon/intron junctions have been sequenced. The beta-casein gene contains 9 exons ranging in size from 21 to 525 bp. We have attempted to identify potential regulatory elements by searching for regions of sequence homology shared between milk protein genes which respond similarly to lactogenic hormones and by searching for previously reported hormone receptor-binding sites. Within the conserved first 200 bp of 5' flanking sequences 3 regions of greater than 70% homology were observed between the rat beta- and gamma-casein genes. One of these contains a region 90% homologous to the chicken progesterone receptor-binding site. The conserved 5' noncoding region, the highly conserved signal peptide, and the hydrophobic carboxyl-terminal region of the protein are each encoded by a separate exon. In contrast the evolutionarily conserved phosphorylation site of beta-casein is formed by an RNA-splicing event. The exons which encode the phosphorylation sites of beta-casein appear to have resulted from an intragenic duplication. Based upon the exon structure of the casein genes, an evolutionary model of intragenic and intergenic exon duplications for this gene family is proposed.     

Evolution of the casein multigene family: conserved sequences in the 5'' flanking and exon regions.

The rat alpha- and bovine alpha s1-casein genes have been isolated and their 5' sequences determined. The rat alpha-, beta-, gamma- and bovine alpha s1-casein genes contain similar 5' exon arrangements in which the 5' noncoding, signal peptide and casein kinase phosphorylation sequences are each encoded by separate exons. These findings support the hypothesis that during evolution, the family of casein genes arose by a process involving exon recruitment followed by intragenic and intergenic duplication of a primordial gene. Several highly conserved regions in the first 200 base pairs of the 5' flanking DNA have been identified. Additional sequence homology extending up to 550 base pairs upstream of the CAP site has been found between the rat alpha- and bovine alpha s1-casein sequences. Unexpectedly, the 5' flanking promoter regions are conserved to a greater extent than both the entire mature coding and intron regions of these genes. These conserved 5' flanking sequences may contain potential cis regulatory elements which are responsible for the coordinate expression of the functionally-related casein genes during mammary gland development.     

The genes for the frog skin peptides GLa, xenopsin, levitide and caerulein contain a homologous export exon encoding a signal sequence and part of an amphiphilic peptide

From genomic libraries of Xenopus laevis, parts of the genes coding for the precursors of the skin peptides GLa (peptide with amino-terminal glycine and carboxy-terminal leucinamide), xenopsin and levitide have been isolated and sequenced. The gene for prepropeptide GLa comprises four exons, separated by relatively small introns. The gene for preproxenopsin is composed of five exons, of which all but the last one have been analyzed. This is a large gene encompassing at least 25,000 base pairs. In addition, two exons of the gene for preprolevitide have been isolated. A comparison of these genes reveals the presence of a homologous exon. This exon contains 161 bp, starts one base pair prior to the initiation codon and encodes a signal peptide and part of a pro region with processing sites. In addition, the two genes for preprocaerulein analyzed previously [Vlasak et al. (1987) Eur. J. Biochem. 169, 53-58] also contain a similar exon. This demonstrates the existence of a homologous export exon in genes encoding the precursors of different skin peptides.     

Conserved exon-intron organization in two different caerulein precursor genes of Xenopus laevis. Additional detection of an exon potentially coding for a new peptide

From a genomic library of Xenopus laevis, two genes coding for different preprocaeruleins have been isolated and sequenced. These correspond to the type I and type III precursors analyzed previously at the cDNA level [Richter, K., Egger, R. and Kreil, G. (1986) J. Biol. Chem. 261, 3676-3680]. The type III gene comprises eight exons; the type I apparently contains eight exons as well, of which six have been sequenced. The genetic information for the dekapeptide caerulein is present on small exons of 45 base pairs. The two genes are highly homologous in their 5'-flanking region, the exon/intron boundaries, and long stretches of intron sequences. A possible scheme for the evolution of this small family of genes through exon and gene duplications is presented. In the type I gene, in place of one of the caerulein exons, a potential exon with conserved splice sites was discovered. If expressed in some frog cells, this exon would code for a new peptide 60% homologous to caerulein.     

Predicting internal exons by oligonucleotide composition and discriminant analysis of spliceable open reading frames.

A new method which predicts internal exon sequences in human DNA has been developed. The method is based on a splice site prediction algorithm that uses the linear discriminant function to combine information about significant triplet frequencies of various functional parts of splice site regions and preferences of oligonucleotides in protein coding and intron regions. The accuracy of our splice site recognition function is 97% for donor splice sites and 96% for acceptor splice sites. For exon prediction, we combine in a discriminant function the characteristics describing the 5'-intron region, donor splice site, coding region, acceptor splice site and 3'-intron region for each open reading frame flanked by GT and AG base pairs. The accuracy of precise internal exon recognition on a test set of 451 exon and 246693 pseudoexon sequences is 77% with a specificity of 79%. The recognition quality computed at the level of individual nucleotides is 89% for exon sequences and 98% for intron sequences. This corresponds to a correlation coefficient for exon prediction of 0.87. The precision of this approach is better than other methods and has been tested on a larger data set. We have also developed a means for predicting exon-exon junctions in cDNA sequences, which can be useful for selecting optimal PCR primers.     

Complete nucleotide sequence of the human corticotropin-beta-lipotropin precursor gene.

The nucleotide sequence of an 8658-base-pair human genomic DNA segment containing the entire corticotropin-beta-lipotropin precursor gene has been determined, and some sequence features of the gene and its flanking regions have been analysed. The gene is composed of 7665 base pairs including two introns of 3708 and 2886 base pairs. Comparison of the 5'-flanking sequences of the human, bovine and mouse corticotropin-beta-lipotropin precursor genes reveals the presence of a highly conserved region, which contains sequences of 14-15 base pairs homologous with sequences located upstream of the mRNA start site of other glucocorticoid-regulated genes.     

The isolation and characterization of the murine T cell antigen receptor zeta chain gene

The T cell antigen receptor (TCR) is a multisubunit complex which has a dual function of antigen recognition and signal transduction. One of its invariant subunits, the zeta chain, has been shown to have a significant role in the expression and function of the TCR on the cell surface. The mouse and human zeta cDNAs share significant homologies to each other but are distinct from all of the previously characterized TCR components. We now report the isolation and structural analysis of the complete murine zeta gene. This gene spans at least 31 kilobases and divides into eight exons. The first exon, which is located at least 20 kilobases upstream from the second exon, codes for the 5'-untranslated region and most of the signal peptide. The second exon codes for the remainder of the signal peptide, the extracellular domain, the transmembrane domain, and the first three amino acids of the intracytoplasmic domain. Exons 3-7 encode the majority of the intracytoplasmic domain. The eight exon encodes the carboxyl-terminal 21 amino acids and the 3'-untranslated region. Four groups of mRNA initiation sites have been identified at approximately 140 base pairs upstream to the AUG codon. No TATA-like box has been detected. The gene is localized to the distal part of chromosome 1 in a linkage group highly conserved between man and mouse.     

Structural organization of a 17 KB segment of the alpha 2 collagen gene: evaluation by R loop mapping.

A recombinant phage, SpC3, containing a 17 kb genomic DNA insert representing approximately 60% of the 3' portion of the sheep collagen alpha 2 gene, was evaluated by electron microscopic R loop analysis. A minimum of 17 intervening sequences (introns) and 18 alpha 2 coding sequences (exons) were mapped. With the exception of the 850 base pair exon located at the extreme 3' end of the insert, all exons contained 250 base pairs or less. The total length of all the exons in SpC3 was 3,014 base pairs. The length distribution of the 17 introns ranged from 300 to 1600 base pairs; together, all of the introns comprised 14,070 base pairs of SpC3 DNA. Thus, the DNA region required for coding the interspersed 3 kb of alpha 2 collagen genetic information was 5.6 fold longer than the corresponding alpha 2 mRNA coding sequences.     

Wrapping of genomic polydA.polydT tracts around nucleosome core particles.

Five human clones containing genomic regions of polydA have been isolated by their ability to form intermolecular triple helices with agarose cross-linked polyU. All of these clones contain Alu repetitive DNA sequences. End-labelled DNA fragments containing these sequences have been successfully reconstituted onto nucleosome core particles by salt exchange. The structure of these has been examined by digesting with DNase I, hydroxyl radicals or diethylpyrocarbonate. DNase I cleavage of the polydA tracts is poor in the free DNA but is markedly enhanced at certain positions when complexed with nucleosome cores. Phased digestion patterns are observed which continue through the (A)n blocks and reveal an average helical periodicity of about 10 base pairs. The distance between adjacent maxima varies between 8-12 base pairs, suggesting that the exact helical repeat is not necessarily constant. One fragment containing the sequence (TA)11T34 reveals a 12 base pair repeat within the (AT)n region. A pUC19 polylinker fragment containing a block of A69.T69 cloned into the Smal site could also be reconstituted onto nucleosome cores and reveals the same phased DNaseI digestion pattern. The DNase I cleavage pattern is not identical at each of the maxima, suggesting that the structural distortions imposed by the core particles are not constant along the DNA.     

Differential splicing in mouse thymus generates two forms of terminal deoxynucleotidyl transferase.

A new form of TdT mRNA has been identified by screening a mouse thymus cDNA library. It contains an open reading frame of 1527 base pairs corresponding to a protein containing 509 aminoacids, whereas the previously identified mouse TdT mRNA is composed of 1587 base pairs and encodes a protein of 529 aminoacids. Analysis of a mouse genomic clone containing the 3' portion of the TdT gene shows that these twenty additional aminoacids are encoded by an additional exon located between exons X and XI. Both forms of TdT mRNA are present in the thymus and could be generated by alternative splicing. The cDNA reported here corresponds to the major form of TdT mRNA in Balb/c mice and closely resembles human and bovine TdT cDNA. Expression of this cDNA in mammalian cells shows that it encodes a functional protein capable of catalysing N region insertions at the recombination junction of an episomic recombination substrate.     

Identification of genetic variation in the bovine major histocompatibility complex DRß-like genes using sequenced bovine genomic probes

Two bovine genomic clones that crosshybridize with HLA-DR beta cDNA have been isolated. Nucleotide sequence analysis of the beta 1, beta 2 and transmembrane (TM) exon regions for one of these clones revealed 70, 89 and 86% identity with the corresponding HLA-DR beta exons. Stop codons are present in the beta 1 and TM exons and a single base deletion toward the 3' end of the TM exon negates the consensus sequence for exon/intron splicing. Therefore, we conclude this is a bovine DR beta-like pseudogene, BoDR beta I. Exon-containing regions have been used as probes in Southern blot analyses of bovine genomic DNA digested with EcoRI. The beta 2 exon of BoDR beta I results in prominent bands of 18.9, 7.8, 7.2, 6.4, 5.6, 3.6, 3.0 and 2.7 kb. Polymorphisms were observed for all but the 18.9 kb band and at least three of these bands were identified in each of the 185 animals sampled. A probe containing the TM exon of BoDR beta I hybridizes only to the 5.6- and 3.6-kb bands, suggesting that these are allelic bands corresponding to this pseudogene. Results from hybridizations of a TM exon-containing probe of the second bovine DR beta-like clone (BoDR beta II) suggest that the 6.4- and 2.7-kb bands correspond to this second bovine gene. A nonpolymorphic 8.1-kb band results from a probe containing the BoDR beta I beta 1 exon. Major differences in frequency for the 6.4/2.7 alleles were found for the four breeds sampled.