Circular RNAs are a class of noncoding RNA with a widespread occurrence in eukaryote tissues, and with some having been demonstrated to have clear biological function. In sheep, little is known about the role of circular RNAs in mammary gland tissue, and therefore an RNA sequencing approach was used to compare mammary gland tissue expression of circular RNAs in 9 Small Tail Han sheep at peak lactation, and subsequently when they were not lactating. These 9 sheep had their RNA pooled for analysis into 3 libraries from peak lactation and 3 from the nonlactating period. A total of 3,278 and 1,756 circular RNAs were identified in the peak lactation and nonlactating mammary gland tissues, respectively, and the expression and identity of 9 of them was confirmed using reverse transcriptase-polymerase chain reaction analysis and DNA sequencing. The type, chromosomal location and length of the circular RNAs identified were ascertained. Forty up-regulated and one downregulated circular RNAs were characterized in the mammary gland tissue at peak lactation compared with the nonlactating mammary gland tissue. Gene ontology enrichment analysis revealed that the parental genes of these differentially expressed circular RNAs were related to molecular function, binding, protein binding, ATP binding, and ion binding. Five differentially expression circular RNAs were selected for further analysis to predict their target microRNAs, and some microRNAs reportedly associated with the development of the mammary gland were found in the constructed circular RNA-microRNA network. This study reveals the expression profiles and characterization of circular RNAs at 2 key stages of mammary gland activity, thereby providing an improved understanding of the roles of circular RNAs in the mammary gland of sheep.

MicroRNAs (miRNAs) have been found to be involved in lipid deposition and metabolism. However, there have been no reports on the roles of miR-148a in the proliferation and adipogenesis of preadipocytes in sheep. In this study, the expression of miR-148a was profiled in the eight tissues of Tibetan ewes and differentiated preadipocytes, and the role of miR-148a in differentiation and proliferation of ovine preadipocytes was investigated using Oil Red O staining, CCK-8, EdU staining, cell cycle detection, and RT-qPCR. The effect of PTEN on the differentiation of ovine preadipocytes was also investigated. The miR-148a was widely expressed in the eight tissues investigated and had significantly increased expression in liver, spleen and subcutaneous adipose tissues, and the heart. The expression of miR-148a continued to increase with the differentiation of ovine preadipocytes. The over-expression of miR-148a significantly promoted differentiation but inhibited the proliferation of ovine preadipocytes. The inhibition of miR-148a had the opposite effect on the differentiation and proliferation of ovine preadipocytes with over-expressed miR-148a. The results from the dual luciferase reporter assays showed that miR-148a mimic significantly decreased the luciferase activity of PTEN-3′UTR dual luciferase reporter vector, suggesting that PTEN is a target gene of miR-148a. In over-expressed-PTEN preadipocytes, the number of lipid droplets remarkably decreased, and the expression levels of adipogenesis marker genes PPARγ, FASN, FATP4, GLUT4, C/EBPβ and LPL were also significantly down-regulated. These results suggest that miR-148a accelerated the adipogenic differentiation of ovine preadipocytes by inhibiting PTEN expression, and also inhibited the proliferation of ovine preadipocytes.

The miR-432 inhibits the proliferation of ovine mammary epithelial cells and milk fat synthesis by targeting SCD and LPL in sheep.

Variation in some caprine keratin-associated protein (KAP) genes has been associated with cashmere fiber traits, but many KAP genes remain unidentified in goats. In this study, we confirm the identification of a KAP27-1 gene (KRTAP27-1) and describe its effect on cashmere traits in 248 Longdong cashmere goats. A polymerase chain reaction–single strand conformation polymorphism (PCR-SSCP) analysis was used to screen for sequence variation in this gene, and three sequence variants (named A to C) were found. These sequences have the highest similarity (77% identity) to a human KRTAP27-1 sequence, while sharing some homology with a predicted caprine KRTAP27-1 sequence ENSCHIG00000023347 in the goat genome construct (ARS1:CM004562.1) at chromosome 1 position 3,966,193–3,973,677 in the forward strand. There were two single nucleotide polymorphisms (SNPs) detected in the coding sequence, including one nonsynonymous SNP (c.413C/T; p.Ala138Val) and one synonymous SNP (c.495C/T). The C variant differed from A and B at c.413C/T, having cytosine in its nucleotide sequence, while the B variant differed from A and C at c.495C/T, having thymine in its nucleotide sequence. Goats of the genotypes AB and BB produced cashmere fibers of higher mean fiber diameter (MFD) than goats of genotype AA, but no difference in MFD was detected between the AB and BB goats. These results suggest that B is associated with increased MFD. Expression of the caprine KRTAP27-1 sequence was predominantly detected in the skin tissue of goats but not or only weakly detected in other tissues, including longissimus dorsi muscle, heart, kidney, liver, lung and spleen.

The mammary gland is a crucial tissue for milk synthesis and plays a critical role in the feeding and growth of mammalian offspring. The aim of this study was to use RNA-sequencing (RNA-Seq) technology to provide a transcriptome profile of the ovine mammary gland at the peak of lactation. Small-Tailed Han (STH) sheep (n = 9) and Gansu Alpine Merino (GAM) sheep (n = 9), breeds with phenotypic differences in milk production traits, were selected for the RNA-Seq analysis. This revealed 74 genes that were more highly expressed in the STHs than in the GAMs. Similarly, 143 genes that were expressed at lower levels in the STHs than in the GAMs, were identified. Gene ontogeny (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses revealed that these differentially expressed genes (DEGs) were associated with binding and catalytic activities, hematopoietic cell lineages, oxytocin signaling pathway and neuroactive ligand–receptor interaction. This is the first study of the transcriptome profile of the ovine mammary gland in these Chinese breeds at peak lactation. The results provide for a better understanding of the genetic mechanisms involved in ovine lactation.

The keratin-associated proteins (KAPs) are constituents of cashmere fibers and variation in many KAP genes (KRTAPs) has been found to be associated with fiber traits. The gene encoding the high-sulphur KAP28-1 has been described in sheep, but it has not been identified in the goat genome. In this study, a 255-bp open reading frame on goat chromosome 1 was identified using a search of similar sequence to ovine KRTAP28-1, and that would if transcribed and translated encode a high sulphur KAP. Based on the analysis of polymerase chain reaction amplicons for the goat nucleotide sequences in 385 Longdong cashmere goats in China, five unique banding patterns were detected using single strand conformation polymorphism analysis. These represented five DNA sequences (named variants A to E) and they had the highest resemblance to KRTAP28-1 sequences from sheep, suggesting A–E are variants of caprine KRTAP28-1. DNA sequencing revealed a 2 or 4-bp deletion and eleven nucleotide sequence differences, including four non-synonymous substitutions. Of the four common variants (A, B, C and D) found in these goats, the presence of variant A was associated with decreased mean fiber diameter and this effect appeared to be additive. These results indicate that caprine KRTAP28-1 variation might have value as a molecular marker for reducing cashmere mean fiber diameter.