Pig fetal skeletal muscle development is associated with genome-wide DNA hypomethylation and corresponding alterations in transcript and microRNA expression.

Fetal myogenesis represents a critical period of porcine skeletal muscle development and requires coordinated expression of thousands of genes. Epigenetic mechanisms, including DNA methylation, drive transcriptional regulation during development; however, these processes are understudied in developing porcine tissues. We performed bisulfite sequencing to assess DNA methylation in pig longissimus dorsi muscle at 41- and 70-days gestation (dg), as well as RNA- and small RNA-sequencing to identify coordinated changes in methylation and expression between myogenic stages. We identified 45 739 differentially methylated regions (DMRs) between stages, and the majority (N = 34 232) were hypomethylated at 70 versus 41 dg. Integration of methylation and transcriptomic data revealed strong associations between differential gene methylation and expression. Differential miRNA methylation was significantly negatively correlated with abundance, and dynamic expression of assayed miRNAs persisted postnatally. Motif analysis revealed significant enrichment of myogenic regulatory factor motifs among hypomethylated regions, suggesting that DNA hypomethylation may function to increase accessibility of muscle-specific transcription factors. We show that developmental DMRs are enriched for GWAS SNPs for muscle- and meat-related traits, demonstrating the potential for epigenetic processes to influence phenotypic diversity. Our results enhance understanding of DNA methylation dynamics of porcine myogenesis and reveal putative cis-regulatory elements governed by epigenetic processes.

Subfamilies of CR1 non-LTR retrotransposons have different 5'UTR sequences but are otherwise conserved.

CR1 elements and CR1-related (CR1-like) elements are a novel family of non-LTR retrotransposons that are found in all vertebrates (reptilia, amphibia, fish, and mammals), whereas more distantly related elements are found in several invertebrate species. CR1 elements have several features that distinguish them from other non-LTR retrotransposons. Most notably, their 3' termini lack a polyadenylic acid (poly A) tail and instead contain 2-4 copies of a unique 8 bp repeat. CR1 elements are present at approximately 100,000 copies in the chicken genome. The vast majority of these elements are severely 5' truncated and mutated; however, six subfamilies (CR1-A through CR1-F) are resolved by sequence comparisons. One of these subfamilies (i.e. CR1-B) previously was analyzed in detail. In the present study, we identified several full-length elements from the CR1-F subfamily. Although regions within the open reading frames and 3' untranslated regions of CR1-F and CR1-B elements are well conserved, their respective 5' untranslated regions are unrelated. Thus, our results suggest that new CR1 subfamilies form when elements with intact open reading frames acquire new 5' UTRs, which could, in principle, function as promoters.

Chicken repeat 1 (CR1) elements, which define an ancient family of vertebrate non-LTR retrotransposons, contain two closely spaced open reading frames.

Chicken repeat 1 (CR1) elements comprise a family of non-long terminal repeat (LTR) retrotransposons that have several noteworthy features. For example, whereas most other non-LTR elements have poly(A) tracts or other simple A-rich repeats at their 3' ends, the 3' ends of CR1 elements conform to the consensus [(CATTCTRT)(GATTCTRT)1-3]. CR1 elements also display an unusual bias for severe 5' truncations: only approx. 30 (out of a total of approx. 30 000) CR1 elements in the chicken genome include significant portions of the pol-like open reading frame (ORF) that we previously identified and partially sequenced [Burch et al. (1993) Proc. Natl. Acad. Sci. USA 90, 8199-8203]. In the present study we derived a consensus sequence for this entire ORF (ORF2) as well as an upstream ORF (ORF1) and part of a 5' untranslated region (UTR). The conceptual translation product of ORF2 is predicted to contain an endonuclease domain in addition to a reverse transcriptase domain. These results suggest that CR1 elements retrotranspose using a "nick and prime" mechanism similar (but not identical) to other families of non-LTR elements.