Cashmere goats are valuable genetic resources which are famous worldwide for their high-quality fiber. Runs of homozygosity (ROHs) have been identified as an efficient tool to assess inbreeding level and identify related genes under selection. However, there is limited research on ROHs in cashmere goats. Therefore, we investigated the ROH pattern, assessed genomic inbreeding levels and examined the candidate genes associated with the cashmere trait using whole-genome resequencing data from 123 goats. Herein, the Inner Mongolia cashmere goat presented the lowest inbreeding coefficient of 0.0263. In total, we identified 57,224 ROHs. Seventy-four ROH islands containing 50 genes were detected. Certain identified genes were related to meat, fiber and milk production (FGF1, PTPRM, RERE, GRID2, RARA); fertility (BIRC6, ECE2, CDH23, PAK1); disease or cold resistance and adaptability (PDCD1LG2, SVIL, PRDM16, RFX4, SH3BP2); and body size and growth (TMEM63C, SYN3, SDC1, STRBP, SMG6). 135 consensus ROHs were identified, and we found candidate genes (FGF5, DVL3, NRAS, KIT) were associated with fiber length or color. These findings enhance our comprehension of inbreeding levels in cashmere goats and the genetic foundations of traits influenced by selective breeding. This research contributes significantly to the future breeding, reservation and use of cashmere goats and other goat breeds.

Chinese Zhongwei goat is a rare and precious fur breed as its lamb fur is a well-known fur product. Wool bending of lamb fur of the Zhongwei goat is its most striking feature. However, the curvature of the wool decreases gradually with growth, which significantly affects its quality and economic value. The mechanism regulating the phenotypic changes of hair bending is still unclear. In the present study, the skin tissues of Zhongwei goats at 45 days (curving wool) and 108 days (slight-curving wool) after birth were taken as the research objects, and the expression profiling of long non-coding RNAs (lncRNAs) and mRNAs were analyzed based on the Ribo Zero RNA sequencing (RNA-seq) method. In total, 46,013 mRNAs and 13,549 lncRNAs were identified, of which 352 were differentially expressed mRNAs and 60 were. lncRNAs. Functional enrichment analysis of the target genes of lncRNAs were mainly enriched in PI3K-Akt, Arachidonic acid metabolic, cAMP, Wnt, and other signaling pathways. The qRT-PCR results of eight selected lncRNAs and target genes were consistent with the sequencing result, which indicated our data were reliable. Through the analysis of the weighted gene co-expression network, 13 co-expression modules were identified. The turquoise module contained a large number of differential expressed lncRNAs, which were mainly enriched in the PI3K-Akt signaling pathway and cAMP signaling pathway. The predicted LOC102172600 and LOC102191729 might affect the development of hair follicles and the curvature of wool by regulating the target genes. Our study provides novel insights into the potential roles of lncRNAs in the regulation of wool bending. In addition, the study offers a theoretical basis for further study of goat wool growth, so as to be a guidance and reference for breeding and improvement in the future.

Skeletal muscle satellite cells (SMSCs), which are highly multifunctional muscle-derived stem cells, play an essential role in myogenesis and regeneration. Here, the transcriptional profile of SMSCs during proliferation and differentiation were constructed using the RNA-Seq method. A total of 1954 differentially expressed genes (DEGs) and 1092 differentially alternative splicing genes (DAGs) were identified including 1288 upregulated genes as well as 666 downregulated genes. GO and KEGG analyses showed that the DEGs and DAGs were enriched in the MAPK (mitogen-activated protein kinase) signaling pathway, the PI3K-Akt (phosphatidylinositol-tris-phosphate kinase 3/protein kinase B) signaling pathway, the Wnt signaling pathway, and the Ras signaling pathway. In total, 1479 alternative splice events (AS) were also identified during SMSC proliferation and differentiation. Among them, a unique AS event was the major per-mRNA splicing type, and SE was the predominant splicing pattern. Furthermore, transcription factors with AS were scanned during SMSC differentiation such as myocyte enhancer factor-2C (MEF2C) and the nuclear receptor subfamily 4 group A member 2 (NR4A2). Our results imply that MEF2C and NR4A2 can interact, and we speculate that NR4A2 and MEF2C might regulate the myogenesis of ovine SMSCs through interaction. Together, our study provides useful information on the transcriptional regulation of SMSCs during proliferation and differentiation at the transcriptional level, and provides a valuable resource for understanding the molecular mechanism of myogenesis and muscle development.

Background: Transposable elements (TEs, or transposons) are repetitive genomic sequences, accounting for half of a mammal genome. Most TEs are transcriptionally silenced, whereas some TEs, especially endogenous retroviruses (ERVs, long terminal repeat retrotransposons), are physiologically expressed in certain conditions. However, the expression pattern of TEs in those less studied species, like goat (Capra hircus), remains unclear. To obtain an overview of the genomic and transcriptomic features of TEs and ERVs in goat, an important farm species, we herein analyzed transcriptomes of ten C. hircus tissues and cells under various physiological and pathological conditions.Method: Distribution of classes, families, and subfamilies of TEs in the C. hircus genome were systematically annotated. The expression patterns of TE-derived transcripts in multiple tissues were investigated at subfamily and location levels. Differential expression of ERV-derived reads was measured under various physiological and pathological conditions, such as embryo development and virus infection challenges. Co-expression between ERV-reads and their proximal genes was also explored to decipher the expression regulation of ERV-derived transcripts.Results: There are around 800 TE subfamilies in the goat genome, accounting for 49.1% of the goat genome sequence. TE-derived reads account for 10% of the transcriptome and their abundance are comparable in various goat tissues, while expression of ERVs are variable among tissues. We further characterized expression pattern of ERV reads in various tissues. Differential expression analysis showed that ERVs are highly active in 16-cell embryos, when the genome of the zygote begins to transcribe its own genes. We also recognized numerous activated ERV reads in response to RNA virus infection in lung, spleen, caecum, and immune cells. CapAeg_1.233:ERVK in chromosome 1 and 17 are dysregulated under endometrium development and infection conditions. They showed strong co-expression with their proximal gene OAS1 and TMPRSS2, indicating the impact of activated proximal gene expression on nearby ERVs.Conclusion: We generated ERV transcriptomes across goat tissues, and identified ERVs activated in response to different physiological and pathological conditions.