Postnatal changes in the expression of genes located in the callipyge region in sheep skeletal muscle

Identification of genes directly responding to DLK1 signaling in Callipyge sheep

BACKGROUND

In food animal agriculture, there is a need to identify the mechanisms that can improve the efficiency of muscle growth and protein accretion. Callipyge sheep provide excellent machinery since the up-regulation of DLK1 and RTL1 results in extreme postnatal muscle hypertrophy in distinct muscles. The aim of this study is to distinguish the genes that directly respond to DLK1 and RTL1 signaling from the genes that change as the result of muscle specific effects.

RESULTS

The quantitative PCR results indicated that DLK1 expression was significantly increased in hypertrophied muscles but not in non-hypertrophied muscles. However, RTL1 was up-regulated in both hypertrophied and non-hypertrophied muscles. Five genes, including PARK7, DNTTIP1, SLC22A3, METTL21E and PDE4D, were consistently co-expressed with DLK1, and therefore were possible transcriptional target genes responding to DLK1 signaling. Treatment of myoblast and myotubes with DLK1 protein induced an average of 1.6-fold and 1.4-fold increase in Dnttip1 and Pde4d expression respectively. Myh4 expression was significantly elevated in DLK1-treated myotubes, whereas the expression of Mettl21e was significantly increased in the DLK1-treated myoblasts but reduced in DLK1-treated myotubes. DLK1 treatment had no impact on Park7 expression. In addition, Park7 and Dnttip1 increased Myh4 and decreased Myh7 promoter activity, resemble to the effects of Dlk1. In contrast, expression of Mettl21e increased Myh7 and decreased Myh4 luciferase activity.

CONCLUSION

The study provided additional supports that RTL1 alone was insufficient to induce muscle hypertrophy and concluded that DLK1 was likely the primary effector of the hypertrophy phenotype. The results also suggested that DNTTIP1 and PDE4D were secondary effector genes responding to DLK1 signaling resulting in muscle fiber switch and muscular hypertrophy in callipyge lamb.

Expressed alleles of imprinted IGF2, DLK1 and MEG3 colocalize in 3D-preserved nuclei of porcine fetal cells

Background

To explore the relationship between spatial genome organization and gene expression in the interphase nucleus, we used a genomic imprinting model, which offers parental-specific gene expression. Using 3D FISH in porcine fetal liver cells, we compared the nuclear organization of the two parental alleles (expressed or not) of insulin-like growth factor 2 (IGF2), a paternally imprinted gene located on chromosome 2. We investigated whether its nuclear positioning favors specific locus associations. We also tested whether IGF2 is implicated in long-range chromatin trans-associations as previously shown in the mouse model species for its reciprocal imprinted gene H19.

Results

We focused on the 3D position of IGF2 alleles, with respect to their individual chromosome 2 territories. The paternally expressed allele was tagged with nascent RNA. There were no significant differences in the position of the two alleles (p = 0.06). To determine long-range chromatin trans-interactions, we chose 12 genes, some of which are known to be imprinted in mammalian model species and belong to a network of imprinted genes (i.e. SLC38A4, DLK1, MEG3, and ZAC1). We screened them and ABCG2, OSBP2, OSBPL1, RPL32, NF1, ZAR1, SEP15, GPC3 for associations with IGF2 in liver cells. All imprinted genes tested showed an association with IGF2. The DLK1/MEG3 locus showed the highest rate of colocalization. This gene association was confirmed by 3D FISH (in 20 % of the nuclei analyzed), revealing also the close proximity of chromosomes 2 and 7 (in 60 % of nuclei). Furthermore, our observations showed that the expressed paternal IGF2 allele is involved in this association. This IGF2-(DLK1/MEG3) association also occurred in a high percentage of fetal muscle cells (36 % of nuclei). Finally, we showed that nascent IGF2, DLK1 and MEG3 RNAs can associate in pairs or in a three-way combination.

Conclusion

Our results show that trans-associations occur between three imprinted genes IGF2, DLK1 and MEG3 both in fetal liver and muscle cells. All three expressed alleles associated in muscle cells. Our findings suggest that the 3D nuclear organization is linked to the transcriptional state of these genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12860-016-0113-9) contains supplementary material, which is available to authorized users.

Delta-like 1 homolog in Capra hircus: molecular characteristics, expression pattern and phylogeny

To research the molecular characteristics, expression pattern and phylogeny of the Delta-like 1 homolog gene (Dlk1) in goats. Dlk1 transcripts were identified in the Jianyang Da'er goats by reverse-transcription polymerase chain reaction (RT-PCR). Phylogenetic trees were constructed by Bayesian inference and neighbor-joining methods. Quantitative real-time PCR (qPCR), western blotting and in situ hybridization were performed to analyze the expression pattern of Dlk1. Five alternatively transcripts were identified in different tissues and designated as Dlk1-AS1, 2, 3, 4 and 5. Compared with the normal transcript Dlk1-AS1, Dlk1-AS4 and Dlk1-AS5 retained the identical open reading frame (ORF) and encoded proteins with truncated epidermal-growth-factor like repeats of 121 and 83 amino acids, respectively. Using the Bayesian inference method, the consensus phylogenetic tree indicated that caprine Dlk1 had a closer relationship with bovine Dlk1 than with Dlk1 from pigs, humans and mice. qPCR revealed high expression levels of Dlk1 in the kidney (P < 0.01). However, mRNA and protein levels presented an inconsistent correlation, possibly because of post-transcriptional regulation. RNA in situ hybridization indicated that Dlk1 mRNA was localized in the interlobular bile duct and alongside the hepatocyte nuclei, in the epithelial cells of proximal and distal convoluted tubules and in the connective region between the mesothelium and myocardium in the heart. The Dlk1 gene in goats produces alternatively spliced transcripts, with specific expression and cellular localization patterns. These findings would lay the foundation for further study.

Regulation of DLK1 by the maternally expressed miR-379/miR-544 cluster may underlie callipyge polar overdominance inheritance

Inheritance of the callipyge phenotype in sheep is an example of polar overdominance inheritance, an unusual mode of inheritance. To investigate the underlying molecular mechanism, we profiled the expression of the genes located in the Delta-like 1 homolog (Dlk1)-type III iodothyronine deiodinase (Dio3) imprinting region in mice. We found that the transcripts of the microRNA (miR) 379/miR-544 cluster were highly expressed in neonatal muscle and paralleled the expression of the Dlk1. We then determined the in vivo role of the miR-379/miR-544 cluster by establishing a mouse line in which the cluster was ablated. The maternal heterozygotes of young mutant mice displayed a hypertrophic tibialis anterior muscle, extensor digitorum longus muscle, gastrocnemius muscle, and gluteus maximus muscle and elevated expression of the DLK1 protein. Reduced expression of DLK1 was mediated by miR-329, a member of this cluster. Our results suggest that maternal expression of the imprinted miR-379/miR-544 cluster regulates paternal expression of the Dlk1 gene in mice. We therefore propose a miR-based molecular working model for polar overdominance inheritance.

Investigation of biochemical changes of the ovine calpain 3 exon-10 polymorphism

Calpain 3 (CAPN3) is a tissue specific calpain, and its mRNA is the most expressed calpain isoform in skeletal muscles. Many mutations and polymorphisms within the human CAPN3 gene have been reported and related to limb-girdle muscular dystrophy. Several reports link CAPN3 polymorphisms and meat quality. An association between three allele variants in exon-10 of ovine CAPN3 and the yield of fat trimmed meat cuts has been reported. This research investigated the biochemical significance of polymorphic variation in CAPN3. CAPN3 mRNA sequences were obtained from muscle samples collected from lambs which were homozygous for each of the three alleles. Four single base substitutions were found besides those in exon-10, but none of them, including the variations within exon-10, caused a change in amino acid sequence. The expression of CAPN3 mRNA and the amounts of CAPN3 protein were also compared among genotypes, and no significant differences were found. These results suggest that the reported association of specific allele variants within CAPN3 exon-10 to phenotype variations were not direct effects of CAPN3 polymorphisms. Interspecies analyses of the CAPN3 sequences indicated that the sequence reported here is more likely to be the correct common ovine CAPN3 sequence than the reference sequence.

New insights into polar overdominance in callipyge sheep

The callipyge phenotype in sheep involves substantial postnatal muscle hypertrophy and other changes to carcass composition. A single nucleotide polymorphism in the DLK1-DIO3 imprinted gene cluster alters gene expression of the paternal allele-specific protein-coding genes and several maternal allele-specific long noncoding RNA and microRNA when the mutation is inherited in cis. The inheritance pattern of the callipyge phenotype is polar overdominant because muscle hypertrophy only occurs in heterozygous animals that inherit a normal maternal allele and the callipyge SNP on the paternal allele (+/C). We examined the changes of gene expression of four major transcripts from the DLK1-DIO3 cluster and four myosin isoforms during the development of muscle hypertrophy in the semimembranosus as well as in the supraspinatus that does not undergo hypertrophy. The homozygous (C/C) animals had an intermediate gene expression pattern for the paternal allele-specific genes and two myosin isoforms, indicating a biological activity that was insufficient to change muscle mass. Transcriptome analysis was conducted by RNA sequencing in the four callipyge genotypes. The data show that homozygous animals (C/C) have lower levels of gene expression at many loci relative to the other three genotypes. A number of the downregulated genes are putative targets of the maternal allele-specific microRNA with gene ontology, indicating regulatory and cell signaling functions. These results suggest that the trans-effect of the maternal noncoding RNA and associated miRNA is to stabilize the expression of a number of regulatory genes at a functional, but low level to make the myofibers of homozygous (C/C) lambs less responsive to hypertrophic stimuli of the paternal allele-specific genes.

Park7 expression influences myotube size and myosin expression in muscle

Callipyge sheep exhibit postnatal muscle hypertrophy due to the up-regulation of DLK1 and/or RTL1. The up-regulation of PARK7 was identified in hypertrophied muscles by microarray analysis and further validated by quantitative PCR. The expression of PARK7 in hypertrophied muscle of callipyge lambs was confirmed to be up-regulated at the protein level. PARK7 was previously identified to positively regulate PI3K/AKT pathway by suppressing the phosphatase activity of PTEN in mouse fibroblasts. The purpose of this study was to investigate the effects of PARK7 in muscle growth and protein accretion in response to IGF1. Primary myoblasts isolated from Park7 (+/+) and Park7 (−/−) mice were used to examine the effect of differential expression of Park7. The Park7 (+/+) myotubes had significantly larger diameters and more total sarcomeric myosin expression than Park7 (−/−) myotubes. IGF1 treatment increased the mRNA abundance of Myh4, Myh7 and Myh8 between 20-40% in Park7 (+/+) myotubes relative to Park7 (−/−). The level of AKT phosphorylation was increased in Park7 (+/+) myotubes at all levels of IGF1 supplementation. After removal of IGF1, the Park7 (+/+) myotubes maintained higher AKT phosphorylation through 3 hours. PARK7 positively regulates the PI3K/AKT pathway by inhibition of PTEN phosphatase activity in skeletal muscle. The increased PARK7 expression can increase protein synthesis and result in myotube hypertrophy. These results support the hypothesis that elevated expression of PARK7 in callipyge muscle would increase levels of AKT activity to cause hypertrophy in response to the normal IGF1 signaling in rapidly growing lambs. Increasing expression of PARK7 could be a novel mechanism to increase protein accretion and muscle growth in livestock or help improve muscle mass with disease or aging.

Genes contributing to genetic variation of muscling in sheep

Selective breeding programs aiming to increase the productivity and profitability of the sheep meat industry use elite, progeny tested sires. The broad genetic traits of primary interest in the progeny of these sires include skeletal muscle yield, fat content, eating quality, and reproductive efficiency. Natural mutations in sheep that enhance muscling have been identified, while a number of genome scans have identified and confirmed quantitative trait loci (QTL) for skeletal muscle traits. The detailed phenotypic characteristics of sheep carrying these mutations or QTL affecting skeletal muscle show a number of common biological themes, particularly changes in developmental growth trajectories, alterations of whole animal morphology, and a shift toward fast twitch glycolytic fibers. The genetic, developmental, and biochemical mechanisms underpinning the actions of some of these genetic variants are described. This review critically assesses this research area, identifies gaps in knowledge, and highlights mechanistic linkages between genetic polymorphisms and skeletal muscle phenotypic changes. This knowledge may aid the discovery of new causal genetic variants and in some cases lead to the development of biochemical and immunological strategies aimed at enhancing skeletal muscle.

Frequency of DLK1 c.639C>T polymorphism and the analysis of MEG3/DLK1/PEG11 cluster expression in muscle of swine raised in Poland

DLK1--(Drosophila like element 1) is a paternally expressed gene, associated with the callipyge phenotype in sheep. In a present study we designed a new real-time PCR alleleic discrimination assay for genotyping of a silent C/T mutation (c.639C>T) in DLK1 gene in swine. The DLK1 c.639C>T mutation was highly polymorphic in all breeds analyzed and C allele was predominant in Landrace and Duroc while T allele was more frequent in Pietrain and Puławska breed. Moreover, we analyzed mRNA expression of DLK1 and adjacent genes--MEG3 and PEG11 in muscles of swines of different breeds raised in Poland. We did not observe significantly different expression of DLK1, MEG3 or PEG11 mRNA in any of analyzed breeds. We also attempted to assess the effect of DLK1 (c.639C>T) on the expression of genes in callipyge locus but did not find significant differences between animals with alternate genotypes (C/C and T/T homozygotes).

Dlk1 is necessary for proper skeletal muscle development and regeneration

Delta-like 1homolog (Dlk1) is an imprinted gene encoding a transmembrane protein whose increased expression has been associated with muscle hypertrophy in animal models. However, the mechanisms by which Dlk1 regulates skeletal muscle plasticity remain unknown. Here we combine conditional gene knockout and over-expression analyses to investigate the role of Dlk1 in mouse muscle development, regeneration and myogenic stem cells (satellite cells). Genetic ablation of Dlk1 in the myogenic lineage resulted in reduced body weight and skeletal muscle mass due to reductions in myofiber numbers and myosin heavy chain IIB gene expression. In addition, muscle-specific Dlk1 ablation led to postnatal growth retardation and impaired muscle regeneration, associated with augmented myogenic inhibitory signaling mediated by NF-κB and inflammatory cytokines. To examine the role of Dlk1 in satellite cells, we analyzed the proliferation, self-renewal and differentiation of satellite cells cultured on their native host myofibers. We showed that ablation of Dlk1 inhibits the expression of the myogenic regulatory transcription factor MyoD, and facilitated the self-renewal of activated satellite cells. Conversely, Dlk1 over-expression inhibited the proliferation and enhanced differentiation of cultured myoblasts. As Dlk1 is expressed at low levels in satellite cells but its expression rapidly increases upon myogenic differentiation in vitro and in regenerating muscles in vivo, our results suggest a model in which Dlk1 expressed by nascent or regenerating myofibers non-cell autonomously promotes the differentiation of their neighbor satellite cells and therefore leads to muscle hypertrophy.

The imprinted retrotransposon-like gene PEG11 (RTL1) is expressed as a full-length protein in skeletal muscle from Callipyge sheep

Members of the Ty3-Gypsy retrotransposon family are rare in mammalian genomes despite their abundance in invertebrates and some vertebrates. These elements contain a gag-pol-like structure characteristic of retroviruses but have lost their ability to retrotranspose into the mammalian genome and are thought to be inactive relics of ancient retrotransposition events. One of these retrotransposon-like elements, PEG11 (also called RTL1) is located at the distal end of ovine chromosome 18 within an imprinted gene cluster that is highly conserved in placental mammals. The region contains several conserved imprinted genes including BEGAIN, DLK1, DAT, GTL2 (MEG3), PEG11 (RTL1), PEG11as, MEG8, MIRG and DIO3. An intergenic point mutation between DLK1 and GTL2 causes muscle hypertrophy in callipyge sheep and is associated with large changes in expression of the genes linked in cis between DLK1 and MEG8. It has been suggested that over-expression of DLK1 is the effector of the callipyge phenotype; however, PEG11 gene expression is also strongly correlated with the emergence of the muscling phenotype as a function of genotype, muscle type and developmental stage. To date, there has been no direct evidence that PEG11 encodes a protein, especially as its anti-sense transcript (PEG11as) contains six miRNA that cause cleavage of the PEG11 transcript. Using immunological and mass spectrometry approaches we have directly identified the full-length PEG11 protein from postnatal nuclear preparations of callipyge skeletal muscle and conclude that its over-expression may be involved in inducing muscle hypertrophy. The developmental expression pattern of the PEG11 gene is consistent with the callipyge mutation causing recapitulation of the normal fetal-like gene expression program during postnatal development. Analysis of the PEG11 sequence indicates strong conservation of the regions encoding the antisense microRNA and in at least two cases these correspond with structural or functional domains of the protein suggesting co-evolution of the sense and antisense genes.

Effect of DLK1 and RTL1 but not MEG3 or MEG8 on muscle gene expression in Callipyge lambs

Callipyge sheep exhibit extreme postnatal muscle hypertrophy in the loin and hindquarters as a result of a single nucleotide polymorphism (SNP) in the imprinted DLK1-DIO3 domain on ovine chromosome 18. The callipyge SNP up-regulates the expression of surrounding transcripts when inherited in cis without altering their allele-specific imprinting status. The callipyge phenotype exhibits polar overdominant inheritance since only paternal heterozygous animals have muscle hypertrophy. Two studies were conducted profiling gene expression in lamb muscles to determine the down-stream effects of over-expression of paternal allele-specific DLK1 and RTL1 as well as maternal allele-specific MEG3, RTL1AS and MEG8, using Affymetrix bovine expression arrays. A total of 375 transcripts were differentially expressed in callipyge muscle and 25 transcripts were subsequently validated by quantitative PCR. The muscle-specific expression patterns of most genes were similar to DLK1 and included genes that are transcriptional repressors or affect feedback mechanisms in beta-adrenergic and growth factor signaling pathways. One gene, phosphodiesterase 7A had an expression pattern similar to RTL1 expression indicating a biological activity for RTL1 in muscle. Only transcripts that localize to the DLK1-DIO3 domain were affected by inheritance of a maternal callipyge allele. Callipyge sheep are a unique model to study over expression of both paternal allele-specific genes and maternal allele-specific non-coding RNA with an accessible and nonlethal phenotype. This study has identified a number of genes that are regulated by DLK1 and RTL1 expression and exert control on postnatal skeletal muscle growth. The genes identified in this model are primary candidates for naturally regulating postnatal muscle growth in all meat animal species, and may serve as targets to ameliorate muscle atrophy conditions including myopathic diseases and age-related sarcopenia.

Analysis on cDNA sequence, mRNA expression and imprinting status of Dlk1 gene in goats

Dlk1 (Delta-like homolog 1) is a cell surface transmembrane glycoprotein belonging to the epidermal growth factor like family of homeotic proteins and plays an important role in regulating fetal and postnatal development. Increased expression of Dlk1 is the primary cause of muscle hypertrophy in the callipyge sheep exhibiting overgrowth of fast-twitch muscles and reduced adiposity. However, the function of Dlk1 in goats remains unknown. In this study, a fragment of 864 bp of goat Dlk1 cDNA, encoding 287 amino acids, was amplified, which has a high homology both in nucleotide sequence and amino acid sequence with the corresponding region of pig, cattle and sheep Dlk1. The Dlk1 was found to be expressed in most tissues of goat fetuses, and in the adrenal gland, pancreas and thymus of adult goats. Two alternative transcripts of Dlk1-C and Dlk1-C2 were expressed in both fetuses and adult goats. One C/T transition in the coding region of goat Dlk1 was identified and by genotyping one segregating goat family and the expressed allele in the tissues of the offspring, Dlk1 was found to be paternally expressed.

Molecular characteristics of the porcine DLK1 and MEG3 genes

Imprinted genes play important roles in embryo survival and postnatal growth regulation. The DLK1 and MEG3 (previously GTL2) genes are linked and reciprocally imprinted in several mammals, but their imprinting status is still unknown in pigs. In this study, we report polymorphisms, imprinting status and QTL analyses of the porcine DLK1 and MEG3 genes. Muscle and adipose DNA and RNA samples from 30-day-old animals generated with reciprocal crosses between the Korean native pig (KNP) and Yorkshire breeds were used to analyse DLK1 and MEG3 variation and expression. The samples exhibited paternal expression of DLK1 and maternal expression of MEG3 in pigs. These results indicated that the imprinting status of the DLK1 and MEG3 genes is conserved across mammalian species. By linkage analyses, we assigned the DLK1 and MEG3 genes to the telomeric region of SSC7. By QTL analyses, we confirmed a significant polar overdominance (POD) effect in DLK1, which was previously detected for several growth traits in pigs. However, no significant POD effect was found with the MEG3 locus.

Myoblasts isolated from hypertrophy-responsive callipyge muscles show altered growth rates and increased resistance to serum deprivation-induced apoptosis

Back and hind limb muscles of sheep paternally heterozygous for the callipyge single nucleotide polymorphism undergo extensive hypertrophy shortly after birth. We have established cell cultures from foetal semitendinosus and longissimus dorsi muscles of normal and callipyge animals. Cultures were assessed for rates of proliferation, cell death, myogenicity and DLK1 expression. Myoblasts from callipyge semitendinosus, but not longissimus dorsi muscles, proliferated faster than myoblasts isolated from normal semitendinosus muscle, and cells isolated from either callipyge muscle were more resistant to serum deprivation-induced apoptosis than equivalent cells isolated from normal individuals. These observations indicate that there are intrinsic differences in the behaviour of isolated myoblasts, which are associated with their muscle and genotype of origin. As myoblasts are the cells responsible for hypertrophy of muscle fibres, the observed differences in cell growth may play a role in the hypertrophy of certain muscles in callipyge animals.

Analysis of the callipyge phenotype through skeletal muscle development; association of Dlk1 with muscle precursor cells

The callipyge mutation in sheep in the form of the paternal heterozygote results in skeletal muscle hypertrophy, which is most pronounced in the hindquarters. Overexpression of one of the genes in the region of the causative single-nucleotide polymorphism, Dlk1, is postulated to be a primary cause of the muscle hypertrophy although the mechanism is not clear. This study examined the expression of Dlk1 mRNA and its encoded protein in skeletal muscles of callipyge and wild-type sheep. The muscles examined included those that demonstrate hypertrophy in callipyge sheep as well as an unaffected muscle. The expression pattern of Dlk1 protein in these muscles was also measured over a developmental time course ranging from 80 days of gestation to 12 weeks after birth. Quantitative reverse transcription-polymerase chain reaction demonstrated that Dlk1 mRNA was significantly increased in affected, but not unaffected, muscles from callipyge sheep at 120 days of gestation through to 12 weeks of age. Immuno-localization of Dlk1 was pronounced in the interstitial connective tissue of fetal muscle but was less intense at later ages. No clear difference in Dlk1 immuno-localization was noted between genotypes in the fetal samples. Strong myofiber-specific Dlk1 immuno-localization was observed in hypertrophied callipyge muscles at 12 weeks of age. This staining was exclusively associated with fast type II myofibers and these had a significantly larger mean cross-sectional area, compared with fast type II myofibers in control sheep that did not overexpress Dlk1. In addition, Dlk1 immuno-localization was associated with a sub-population of Pax7-positive mononucleated cells in all skeletal muscles examined during fetal development and at birth, but this was not apparent at 12 weeks. There were no genotype-dependent alterations in the mRNA expression patterns of a number of promyogenic transcription factors indicating that the callipyge mutation was not affecting muscle cell differentiation per se. We postulate that Dlk1 is implicated in the commitment and/or proliferation of fetal myoblasts as well as in the maintenance of hypertrophy in fully differentiated myofibers.

Analysis of gene expression during the onset of muscle hypertrophy in callipyge lambs

The callipyge mutation causes postnatal muscle hypertrophy in heterozygous lambs that inherit a paternal callipyge allele (+/CLPG). Our hypothesis was that the up-regulation of one or both of the affected paternally expressed genes (DLK1 or PEG11) initiates changes in biochemical and physiological pathways in skeletal muscle to induce hypertrophy. The goal of this study was to identify changes in gene expression during the onset of muscle hypertrophy to identify the pathways that are involved in the expression of the callipyge phenotype. Gene expression was analysed in longissimus dorsi total RNA from lambs at 10, 20, and 30 days of age using the Affymetrix Bovine Expression Array. An average of 40.6% of probe sets on the array was detected in sheep muscle. Data were normalized and analysed using a two-way anova for genotype and age effects with a false discovery rate of 0.10. From the anova, 13 genes were significant for the effect of genotype and 13 were significant for effect of age (P < 0.10). No significant age-by-genotype interactions were detected (P > 0.10). Of the 13 genes indicating an effect of genotype, quantitative PCR assays were developed for all of them and tested on a larger group of animals from 10 to 200 days of age. Nine genes had significantly elevated transcript levels in callipyge lambs. These genes included phosphofructokinase, a putative methyltransferase protein, a cAMP phosphodiesterase, and the transcription factor DNTTIP1.