Ovulation rate and litter size are important reproductive traits in sheep with high economic value. Recent work has revealed a potential link between DNA methylation and prolificacy. However, a genome-wide study that sought to identify potential DNA methylation sites involved in sheep prolificacy indicated that it is still unknown. Here, we aimed to investigate the genome-wide DNA methylation profiles of Hu sheep ovaries by comparing a high-prolificacy group (HP, litter size of three for at least 2 consecutive lambings) and low prolificacy group (LP, litter size of one for at least 2 consecutive lambings) using deep whole-genome bisulfite sequencing (WGBS). First, our results demonstrated lower expression levels of DNA methyltransferase (DNMT) genes in the ovaries of the HP group than that in the ovaries of the LP group. Both groups showed similar proportions of methylation at CpG sites but different proportions at non-CpG sites. Subsequently, we identified 70,899 differential methylated regions (DMRs) of CG, 16 DMRs of CHG, 356 DMRs of CHH and 12,832 DMR-related genes(DMGs). Gene Ontology (GO) analyses revealed that some DMGs were involved in regulating female gonad development and ovarian follicle development. Finally, we found that 10 DMGs, including BMP7, BMPR1B, CTNNB1, FST, FSHR, LHCGR, TGFB2 and TGFB3, are more likely to be involved in prolificacy of Hu sheep, as assessed by correlation analysis and listed in detail. This study revealed the global DNA methylation pattern of sheep ovaries associated with high and low prolificacy groups, which may contribute to a better understanding of the epigenetic regulation of sheep reproductive capacity.

Myostatin (Mstn) is a conserved negative regulator of skeletal muscle mass in mammals. However, whether precise disruption of Mstn in livestock can be achieved and safely used to improve meat productivity has not been proven. We applied CRISPR/Cas9 system to generate Mstn knock-out (KO) rabbits and goats and then analyzed the changes in their phenotypes to answer this question. We efficiently generated 24 Mstn KO rabbits out of 32 newborn infants after embryo injection with two sgRNAs targeting rabbit Mstn, and found that the Mstn KO rabbits exhibited increased birthweight and a significantly increase in the weight ratios of the quadriceps and biceps muscles to the whole body. Mstn KO also caused high probability of enlarged tongue phenomenon and severe health problems such as stillbirth and early stage death. Using the same method, one out of four goats was generated with edition at Mstn locus. The early stage growth rate of this goat outperformed the control goats. In conclusion, we efficiently generated Mstn KO rabbits and goats using CRISPR/Cas9 technology. However, Mstn KO causes severe health problems and may also have the same effects on other species. This safety issue must be studied further before applied to animal reproduction processes.

Goat’s milk, considered a substitute for cow’s milk, has a high nutritional value. However, goat’s milk contains various allergens, predominantly β-lactoglobulin (BLG). In this study, we employed the CRISPR/Cas9 system to target the BLG locus in goat fibroblasts for sgRNA optimization and generate BLG knock-out goats through co-injection of Cas9 mRNA and small guide RNAs (sgRNAs) into goat embryos at the one-cell stage. We firstly tested sgRNA editing efficiencies in goat fibroblast cells, and approximately 8.00%–9.09% of the cells were modified in single sgRNA-guided targeting experiment. Among the kids, the genome-targeting efficiencies of single sgRNA were 12.5% (10 ng/μL sg1) and 0% (10 ng/μL sg2) and efficiencies of dual sgRNAs were 25.0% (25 ng/μL sg2+sg3 group) and 28.6% (50 ng/μL sg2+sg3 group). Relative expression of BLG in BLG knock-out goat mammary glands significantly (p < 0.01) decreased as well as other milk protein coding genes, such as CSN1S1, CSN1S2, CSN2, CSN3 and LALBA (p < 0.05). As expected, BLG protein had been abolished in the milk of the BLG knock-out goat. In addition, most of the targeted kids were chimeric (3/4), and their various body tissues were edited simultaneously. Our study thus provides a basis for optimizing the quality of goat milk, which can be applied to biomedical and agricultural research.

The aryl hydrocarbon receptor (AhR) mediates induction of CYP1A1 and CYP1B1 by 2,3,7,8-tetrachlorodibenzo-ρ-dioxin (dioxin) via binding to xenobiotic-responsive elements (XREs) in their enhancer regions. CYP1A1 and CYPIB1 were both inducible by dioxin in human MCF-7 cells. However, only CYP1A1 was inducible in human HepG2 cells. Further experiments focused on providing an explanation for this last observation. Dioxin induced the recruitment of AHR and the transcriptional coactivators p300 and p300/cAMP response element-binding protein binding protein-associated factor (PCAF) to the CYP1B1 enhancer in HepG2 cells but failed to induce recruitment of RNA polymerase II (polII) or the TATA binding protein (TBP) and acetylations of histones 3 and 4 or methylation of histone 3 at the promoter. Because p300 was required for dioxin induction of the aforementioned histone modifications at the CYP1B1 promoter and for induction of CYP1B1 transcription (in MCF-7 cells), the recruitments of p300 and AhR, although necessary, are not sufficient for eliciting the above responses to dioxin. Cytosine residues within CpG dinucleotides at the enhancer, including those within the XREs, were partially methylated, whereas those at the promoter were fully methylated. Treatment of HepG2 cells with 5-aza-2′-deoxycytidine led to partial demethylation of the promoter, restored polII and TBP binding, and CYP1B1 inducibility. Thus, the deficiency of CYP1B1 induction in HepG2 cells is ascribable to cytosine methylation at the promoter, which prevents recruitment of TBP and polII. It is noteworthy that our data indicate that stable recruitment of p300 and PCAF to the CYP1B1 gene does not require their tethering to the promoter and to the enhancer.

N6-methyladenosine (m⁶A) modification plays a critical role in mammalian development. However, the role of m⁶A in the skeletal muscle development remains largely unknown. Here, we report a global m⁶A modification pattern of goat skeletal muscle at two key development stages and identified that the m⁶A modification regulated the expression of the growth arrest and DNA damage-inducible 45B (GADD45B) gene, which is involved in myogenic differentiation. We showed that GADD45B expression increased during myoblast differentiation, whereas the downregulation of GADD45B inhibits myogenic differentiation and mitochondrial biogenesis. Moreover, the expression of GADD45B regulates the expression of myogenic regulatory factors and peroxisome proliferator-activated receptor gamma coactivator 1 alpha by activating the p38 mitogen-activated protein kinase (MAPK) pathway. Conversely, the inactivation of p38 MAPK abolished the GADD45B-mediated myogenic differentiation. Furthermore, we found that the knockdown of fat mass and obesity-associated protein (FTO) increases GADD45B m⁶A modification and decreases the stability of GADD45B mRNA, which impairs myogenic differentiation. Our results indicate that the FTO-mediated m⁶A modification in GADD45B mRNA drives skeletal muscle differentiation by activating the p38 MAPK pathway, which provides a molecular mechanism for the regulation of myogenesis via RNA methylation.

Long non-coding RNA (LncRNA) have been identified as important regulators in the hypothalamic-pituitary-ovarian axis associated with sheep prolificacy. However, their expression pattern and potential roles in the pituitary are yet unclear. To explore the potential mRNAs and lncRNAs that regulate the expression of the genes involved in sheep prolificacy, we used stranded specific RNA-seq to profile the pituitary transcriptome (lncRNA and mRNA) in high prolificacy (genotype FecB BB, litter size = 3; H) and low prolificacy sheep (genotype FecB B+; litter size = 1; L). Our results showed that 57 differentially expressed (DE) lncRNAs and 298 DE mRNAs were found in the pituitary between the two groups. The qRT-PCR results correlated well with the RNA-seq results. Moreover, functional annotation analysis showed that the target genes of the DE lncRNAs were significantly enriched in pituitary function, hormone-related pathways as well as response to stimulus and some other terms related to reproduction. Furthermore, a co-expression network of lncRNAs and target genes was constructed and reproduction related genes such as SMAD2, NMB and EFNB3 were included. Lastly, the interaction of candidate lncRNA MSTRG.259847.2 and its target gene SMAD2 were validated in vitro of sheep pituitary cells. These differential mRNA and lncRNA expression profiles provide a valuable resource for understanding the molecular mechanisms underlying Hu sheep prolificacy.