Exon skipping in human beta-casein

Casein Gene Cluster in Camelids: Comparative Genome Analysis and New Findings on Haplotype Variability and Physical Mapping

The structure of casein genes has been fully understood in llamas, whereas in other camelids, this information is still incomplete. In fact, structure and polymorphisms have been identified in three (CSN1S1, αs1-CN; CSN2, β-CN; CSN3, κ-CN) out of four casein genes, whereas controversial information is available for the CSN1S2 (αs2-CN) in terms of structure and genetic diversity. Data from the genome analysis, whose assembly is available for feral camel, Bactrian, dromedary, and alpaca, can contribute to a better knowledge. However, a majority of the scaffolds available in GenBank are still unplaced, and the comparative annotation is often inaccurate or lacking.Therefore, the aims of this study are 1) to perform a comparative genome analysis and synthesize the literature data on camelids casein cluster; 2) to analyze the casein variability in two dromedary populations (Sudanese and Nigerian) using polymorphisms at CSN1S1 (c.150G > T), CSN2 (g.2126A > G), and CSN3 (g.1029T > C); and 3) to physically map the casein cluster in alpaca. Exon structures, gene and intergenic distances, large insertion/deletion events, SNPs, and microsatellites were annotated. In all camelids, the CSN1S2 consists of 17 exons, confirming the structure of llama CSN1S2 gene. The comparative analysis of the complete casein cluster (∼190kb) shows 12,818 polymorphisms. The most polymorphic gene is the CSN1S1 (99 SNPs in Bactrian vs. 248 in dromedary vs. 626 in alpaca). The less polymorphic is the CSN3 in the Bactrian (22 SNPs) and alpaca (301 SNPs), whereas it is the CSN1S2 in dromedary (79 SNPs). In the two investigated dromedary populations, the allele frequencies for the three markers are slightly different: the allele C at CSN1S1 is very rare in Nigerian (0.054) and Sudanese dromedaries (0.094), whereas the frequency of the allele G at CSN2 is almost inverted. Haplotype analysis evidenced GAC as the most frequent (0.288) and TGC as the rarest (0.005). The analysis of R-banding metaphases hybridized with specific probes mapped the casein genes on chromosome 2q21 in alpaca. These data deepen the information on the structure of the casein cluster in camelids and add knowledge on the cytogenetic map and haplotype variability.

Molecular Characterization of the Llamas (Lama glama) Casein Cluster Genes Transcripts (CSN1S1, CSN2, CSN1S2, CSN3) and Regulatory Regions

In the present paper, we report for the first time the characterization of llama (Lama glama) caseins at transcriptomic and genetic level. A total of 288 casein clones transcripts were analysed from two lactating llamas. The most represented mRNA populations were those correctly assembled (85.07%) and they encoded for mature proteins of 215, 217, 187 and 162 amino acids respectively for the CSN1S1, CSN2, CSN1S2 and CSN3 genes. The exonic subdivision evidenced a structure made of 21, 9, 17 and 6 exons for the αs1-, β-, αs2- and κ-casein genes respectively. Exon skipping and duplication events were evidenced. Two variants A and B were identified in the αs1-casein gene as result of the alternative out-splicing of the exon 18. An additional exon coding for a novel esapeptide was found to be cryptic in the κ-casein gene, whereas one extra exon was found in the αs2-casein gene by the comparison with the Camelus dromedaries sequence. A total of 28 putative phosphorylated motifs highlighted a complex heterogeneity and a potential variable degree of post-translational modifications. Ninety-six polymorphic sites were found through the comparison of the lama casein cDNAs with the homologous camel sequences, whereas the first description and characterization of the 5'- and 3'-regulatory regions allowed to identify the main putative consensus sequences involved in the casein genes expression, thus opening the way to new investigations -so far- never achieved in this species.

Mineralized tissue and vertebrate evolution: the secretory calcium-binding phosphoprotein gene cluster

Gene duplication creates evolutionary novelties by using older tools in new ways. We have identified evidence that the genes for enamel matrix proteins (EMPs), milk caseins, and salivary proteins comprise a family descended from a common ancestor by tandem gene duplication. These genes remain linked, except for one EMP gene, amelogenin. These genes show common structural features and are expressed in ontogenetically similar tissues. Many of these genes encode secretory Ca-binding phosphoproteins, which regulate the Ca-phosphate concentration of the extracellular environment. By exploiting this fundamental property, these genes have subsequently diversified to serve specialized adaptive functions. Casein makes milk supersaturated with Ca-phosphate, which was critical to the successive mammalian divergence. The innovation of enamel led to mineralized feeding apparatus, which enabled active predation of early vertebrates. The EMP genes comprise a subfamily not identified previously. A set of genes for dentine and bone extracellular matrix proteins constitutes an additional cluster distal to the EMP gene cluster, with similar structural features to EMP genes. The duplication and diversification of the primordial genes for enameldentinebone extracellular matrix may have been important in core vertebrate feeding adaptations, the mineralized skeleton, the evolution of saliva, and, eventually, lactation. The order of duplication events may help delineate early events in mineralized skeletal formation, which is a major characteristic of vertebrates.

Differential splicing of pre-messenger RNA produces multiple forms of mature caprine alpha(s1)-casein

The identity of multiple forms of caprine alpha(s1)-casein in variants A, B, and C has been determined by structural characterisation using mass spectrometry, automated Edman degradation and peptide mapping. Mature goat alpha(s1)-casein exists as a mixture of at least four molecular species which differ in peptide chain length. The main component corresponds to the 199-residues form already described. The other three, in lesser amounts, were shorter forms of alpha(s1)-casein and differed for the deleted peptides 141-148, as shown previously for ovine alpha(s1)-casein, peptide 110-117, or Gln78. Analysis of alpha(s1)-casein mRNA from milk somatic cells demonstrated that these forms originated from skipping events at the level of exon 13 (codifying for peptide 110-117) and 16 (codifying for peptide 141-148) and from the presence of a cryptic splice site within exon 11 (whose first CAG triplet encodes Gln78) during primary transcript processing. The finding of these splicing abnormalities in the three common variants A, B, and C suggests that this is a general feature of alpha(s1)-casein in goat. A further source of heterogeneity of caprine alpha(s1)-casein was identified in the discrete phosphorylation of seryl residues. Eight serine residues (at positions 44, 46, 64 to 68 and 75) are fully phosphorylated (except in variant A because of the replacement Glu77-->Gln which prevents phosphorylation of Ser75). Conversely, Ser115 and Ser41 are phosphorylated only to about 50% and 20%, respectively. Ser12, although located in a consensus triplet, is never phosphorylated, similarly to the ovine alpha(s1)-casein variants. These results confirm that there are stabilised mechanisms of simultaneous synthesis of alpha(s1)-casein at different length and of post-translational modification in both caprine and ovine species.

Structure and expression of the mouse casein gene locus

The analysis of yeast artificial chromosomes (YACs) containing the complete mouse casein gene locus revealed the presence of five casein genes, alpha-, beta-, gamma-, delta-, and kappa-casein, in this order, in the locus. The alpha- and beta-casein genes are only 10 kb apart and have convergent transcriptional orientations. The distance between the beta-casein gene and the alpha s2-like gamma-casein gene is about 70 kb, and these genes have divergent transcriptional orientations. The gamma- and delta-casein genes, both encoding a alpha s2-like casein, are linked within 60 kb and convergently transcribed. The kappa-casein gene is located about 100 kb from the delta-gene. Except for the presence of the delta-casein gene, the organization of the mouse casein locus resembles that of the bovine locus, including the transcriptional orientation of the genes. In contrast to the other casein genes, which are strongly induced at mid-lactation, expression of the delta-casein gene is abruptly induced upon parturition. Comparative analysis of alpha s2-like sequences from various species suggests that the ancestral alpha s2-like gene duplicated around the time of radiation of the rodent and artiodactylid ancestors.

Identification of a novel opioid peptide (Tyr-Val-Pro-Phe-Pro) derived from human alpha S1 casein (alpha S1-casomorphin, and alpha S1-casomorphin amide)

A new casomorphin pentapeptide (alpha S1-casomorphin) has been isolated from the sequence of human alpha S1-casein [alpha S1-casein-(158-162)], with the sequence Tyr-Val-Pro-Phe-Pro. This peptide was found to bind with high affinity to all three subtypes of the kappa-opioid receptor (kappa 1-kappa 2). When amidated at the C-terminus, alpha S1-casomorphin amide binds to the delta- and kappa 3-opioid sites. Both alpha S1-casomorphin and its amide inhibit in a dose-dependent and reversible manner the proliferation of T47D human breast cancer cells. This anti-proliferative activity was greater for alpha S1-casomorphin, which was the most potent opioid in inhibiting T47D cell proliferation. In T47D breast cancer cells, other casomorphins have been found to bind to somatostatin receptors in addition to opioid sites. In contrast, alpha S1-casomorphin and its amide do not interact with somatostatin receptors in our system.

Exon skipping in the ovine alpha s1-casein gene

The reported cDNA sequences for the bovine (Bos taurus) and ovine (Ovis aries) alpha s1-caseins display a high degree of identity with the exception that a 24 bp region, corresponding to bovine exon 16, is absent in the ovine sequence. Here we show that the ovine gene for alpha s1-casein contains a sequence block displaying 23/24 identity to bovine exon 16, indicating that the absence of this block from ovine mRNA is due not to genomic deletion but to exon skipping. Analysis of the products obtained by reverse transcription of ovine alpha s1-casein mRNA followed by amplification, demonstrated the presence of mRNA species containing the exon 16 sequence as well as the species in which it had been spliced out. It was estimated that the latter constitutes 20% of the total ovine alpha s1-casein mRNA. We propose that a substitution within the donor splice site is responsible for the partial skipping of exon 16, possibly through the formation of an inhibitory RNA secondary structure.

Complete sequence of the ovine beta-casein-encoding gene and interspecies comparison

The 9149-bp transcription unit encoding ovine beta-casein (Cas) and 4636 bp of 5' flanking region were completely sequenced. The gene is composed of nine exons and its overall organization is similar to that of its counterparts from other species. Intron 4, the largest, shares three similar stretches (sizes ranging from 0.1 to 0.3 kb) with the region upstream from the transcription unit. These common sequences are part of short interspersed nuclear elements (SINE) specific to Bovidae (Bov). Intron 4 contains two 274-bp Bov-A2 SINE in opposite orientation, as well as a full-length 569-bp Bov-B SINE. This latter SINE, also present in caprine intron 4, is missing in cattle. This suggests that the amplification of Bov-SINE has continued after the divergence of cattle from sheep and goats, assuming that the presently known sequences are representative of these species.

Structure of the human beta-casein encoding gene

The entire human beta-casein-encoding gene, Bca, was cloned and sequenced. The gene consists of eight exons ranging from 21 to 531 nucleotides (nt) in length and extending over 10,466 nt. Exon-2 contains the translational start, the entire signal sequence and the codons for the two first amino acids of the mature protein. This corresponds to the organization found in other species. The translational stop is localized to exon-7. Exon/intron boundaries are in accordance with the AG/GT rule and conform to suggested consensus sequences. Splice junctions are located between coding triplets. In all other species analyzed, Bca has been found to consist of nine exons; however, within intron-2 of the human gene, a sequence omitted from human mRNA, but corresponding to exon-3 of other known Bca genes, was revealed.

Structure and function of milk protein genes

Interspecies comparisons of cDNA and mosaic milk protein genes have confirmed their high rate of evolution, but the overall gene organization has been conserved. The three Ca-sensitive casein genes, which share common motifs in the promoter region and contain similar sequences that encode signal peptide and multiple phosphorylation sites, probably derived from a common ancestor. alpha s1- and alpha s2-casein genes, divided into many small exons, undergo complex splicing, and the deleted caseins arise from exon skipping. The four bovine casein genes are clustered on 200 kb of chromosome 6. alpha-Lactalbumin and beta-lactoglobulin pseudogenes occur in ruminants. Study of the expression of native and modified milk protein genes in mammary cell lines and transgenic animals and DNA footprinting have shown the occurrence of important regulatory motifs in the proximal 5' flanking region, including one recognized by a specific mammary nuclear factor. Good stage- and tissue-specific expression has been obtained in transgenic animals with milk protein genes having less than a 3-kb 5' flanking region. Better knowledge of both the structure and function of milk protein genes, which has already allowed the use of powerful techniques for the rapid identification of alleles, offers the potential for the genetic modification of milk composition.

Bovine alpha s2-casein D is generated by exon VIII skipping

Bovine alpha s2-casein D (CasD) differs from the common type A by the deletion of a stretch of 9 amino acids (aa) starting at a position not precisely known, either at aa 50, 51, or 52. The sequence of cloned PCR-amplified genomic DNA from three homozygous cows, two unrelated females carrying the CasD allele and one carrying the CasA allele, did not reveal any deletion and showed two identical nucleotide (nt) substitutions in the 1.7-kb region of both CasD alleles encompassing codons 43-75 in the cDNA encoding alpha s2-CasA. This strongly suggests that the deleted bovine alpha s2-CasD arises from skipping the 27-nt exon, now identified as exon VIII, which encodes aa 51-59 of alpha s2-CasA. The G-->T transversion (allele A-->D) affecting the last nt of exon VIII, i.e., the 5' consensus splicing site, might be responsible for the altered splicing of the primary transcript of alpha s2-CasD.