Superior kidding rate is an important economic trait in production of meat goat, and ovulation rate is the precondition of kidding rate. MicroRNAs (miRNAs) play critical roles in almost all ovarian biological processes, including folliculogenesis, follicle development, follicle atresia, luteal development and regression. To find out the different ovarian activity and follicle recruitment with miRNA-mediated posttranscriptional regulation, the small RNAs expressed pattern in the ovarian tissues of multiple and uniparous Anhui White goats during follicular phase was analyzed using Solexa sequencing data. 1008 miRNAs co-expressed, 309 and 433 miRNAs specifically expressed in the ovaries of multiple and uniparous goats during follicular phase were identified. The 10 most highly expressed miRNAs in the multiple library were also the highest expressed in the uniparous library, and there were no significantly different between each other. The highest specific expressed miRNA in the multiple library was miR-29c, and the one in the uniparous library was miR-6406. 35 novel miRNAs were predicted in total. GO annotation and KEGG Pathway analyses were implemented on target genes of all miRNA in two libraries. RT-PCR was applied to detect the expression level of 5 randomly selected miRNAs in multiple and uniparous hircine ovaries, and the results were consistent with the Solexa sequencing data. In the present study, the different expression of miRNAs in the ovaries of multiple and uniparous goats during follicular phase were characterized and investigated using deep sequencing technology. The result will help to further understand the role of miRNAs in kidding rate regulation and also may help to identify miRNAs which could be potentially used to increase hircine ovulation rate and kidding rate in the future.

Asian corn borer (ACB), Ostrinia furnacalis (Guenée), is the major insect pest of maize in China and countries of East and Southeast Asia, the Pacific and Australasia. ACB can develop strong resistance to the transgenic Bt maize expressing Cry1Ab, the most widely commercialized Bt maize worldwide. However, the molecular basis for the resistance mechanisms of ACB to Cry1Ab remained unclear. Two biological replicates of the transcriptome of Bt susceptible (ACB-BtS) and Cry1Ab resistant (ACB-AbR) strains of ACB were sequenced using Solexa/Illumina RNA-Seq technology to identify Cry1Ab resistance-relevant genes. The numbers of unigenes for two biological replications were 63,032 and 53,710 for ACB-BtS and 57,770 and 54,468 for ACB-AbR. There were 35,723 annotated unigenes from ACB reads found by BLAST searching NCBI non-redundant, NCBI non-redundant nucleotide, Swiss-prot protein, Kyoto Encyclopedia of Genes and Genomes, Cluster of Orthologous Groups of proteins, and Gene Ontology databases. Based on the NOISeq method, 3,793 unigenes were judged to be differentially expressed between ACB-BtS and ACB-AbR. Cry1Ab resistance appeared to be associated with change in the transcription level of enzymes involved in growth regulation, detoxification and metabolic/catabolic process. Among previously described Bt toxin receptors, the differentially expressed unigenes associated with aminopeptidase N and chymotrypsin/trypsin were up-regulated in ACB-AbR. Whereas, other putative Cry receptors, cadherin-like protein, alkaline phosphatase, glycolipid, actin, V-type proton ATPase vatalytic, heat shock protein, were under-transcripted. Finally, GPI-anchor biosynthesis was found to be involved in the significantly enriched pathway, and all genes mapped to the pathway were substantially down-regulated in ACB-AbR. To our knowledge, this is the first comparative transcriptome study to discover candidate genes involved in ACB Bt resistance. This study identified differentially expressed unigenes related to general Bt resistance in ACB. The assembled, annotated transcriptomes provides a valuable genomic resource for further understanding of the molecular basis of ACB Bt resistance mechanisms.

The goat is an important farm animal. Reproduction is an important process of goat farming. The ovary is the most important reproductive organ for goats. In recent years, an increasing number of long non-coding RNAs (lncRNAs) have been implicated in the regulation of mammal reproduction. However, there are few studies on the function of lncRNAs in reproduction, particularly lncRNAs in the ovary. The sequencing of goat ovaries generated 1,122,014,112 clean reads, and 4926 lncRNAs and 1454 TUCPs (transcripts of uncertain coding potential) were identified for further analysis by using the coding potential analysis software, CNCI, CPC and Pfam-sca. There were 115 /22 differential lncRNAs /TUCPs transcripts between the ovaries of the luteal phase and the follicular phase. We predicted the related genes of lncRNA /TUCP based on co-expression and co-localization methods. In total, 2584 /904 genes were predicted by co-expression, and 326/73 genes were predicted by co-localization. The functions of these genes were further analyzed with GO and KEGG analysis. The results showed that lncRNAs /TUCPs, which are highly expressed in goat ovaries in the luteal phase, are mainly associated with the synthesis of progesterone, and we filtered the lncRNAs /TUCPs, such as XR_001918177.1 and TUCP_001362, which may regulate the synthesis of progesterone; lncRNAs /TUCPs, which are highly expressed in goat ovaries in the follicular phase, are mainly associated with oogenesis and the maturation of oocytes, and we filtered the lncRNAs /TUCPs that may regulate the oogenesis and maturation of oocyte, such as XR_001917388.1 and TUCP_000849. The present study provided the genome expression profile of lncRNAs /TUCPs in goat ovaries at different estrus periods and filtered the potential lncRNAs /TUCPs associated with goat reproduction. These results are helpful to further study the molecular mechanisms of goat reproduction.

A series of complex processes regulate muscle development, and lncRNAs play essential roles in the regulation of skeletal myogenesis. Using RNA sequencing, we profiled the lncRNA expression during goat (Capra hircus) skeletal muscle development, which included seven stages across fetal 45 (F45), 65 (F65), 90 (F90), 120 (F120), 135 (F135) days, born for 24 h (B1) and 90 (B90) days. A total of 15,079 lncRNAs were identified in the seven stages, and they were less conservative with other species (human, cow, and mouse). Among them, 547 were differentially expressed, and they divided the seven stages into three functional transition periods. Following weighted gene co-expression network analysis (WGCNA), five lncRNA modules specific for developmental stages were defined as three types: ‘Early modules’, ‘late modules’, and ‘individual-stage-specific modules’. The enrichment content showed that ‘early modules’ were related to muscle structure formation, ‘late modules’ participated in the ‘p53 signaling pathway’ and other pathways, the F90-highly related module was involved in the ‘MAPK signaling pathway’, and other pathways. Furthermore, we identified hub-lncRNA in three types of modules, and LNC_011371, LNC_ 007561, and LNC_001728 may play important roles in goat skeletal muscle. These data will facilitate further exploration of skeletal muscle lncRNA functions at different developmental stages in goats.

The desmoglein 4 (DSG4) gene is a potential candidate in the search for genes that may affect wool traits, because of its function. This study aimed to screen for polymorphisms in partial exon 16 and 3ꞌUTR of the sheep desmoglein 4 DSG4 gene, and to test its possible association with wool length and crimp associated with fur. Overall, 326 sheep were scanned via single-strand conformational polymorphism assay, through three pairs of primers. The breeds included Tan, Han, and TanxHan from China, Polled Dorset from Australia, and Suffolk from Britain genotypes AA, BB, and AB for primer2 and genotypes DD, EE, and DE for primer3 were detected in native breeds. Six SNPs and 3-bp insertion/deletions were found in exon 16, of which 4 lead to amino acid substitutions. In addition, 1 SNP was found in 3ꞌUTR. The DSG4 genotype was found to be strongly associated with all wool traits that were considered in this study (P < 0.01). Sheep with the genotype MM had a higher least square mean compared to sheep with the genotype WW or WM with respect to birth scapular wool length (P < 0.01), crimp number of birth scapular wool crimp (P < 0.01), crimp number of weaning scapular wool crimp (P < 0.01), and crimp number of weaning rump wool crimp (P < 0.01, P < 0.05). In conclusion, our study is the first to demonstrate that the DSG4 gene may be a candidate, or major gene, which influences important wool traits.

Background: Circular RNA (circRNA) is produced during the splicing of mRNA (in addition to linear splicing) and is part of the gene regulatory network. The temporal expression patterns the different developmental stages were inseparable from these molecules’ function. Results: Skeletal muscles of Anhui white goat (AWG) across seven fetal to postnatal development stages were sequenced and 21 RNA sequencing libraries were constructed. We thereby identified 9090 circRNAs and analyzed their molecular properties, temporal expression patterns, and potential functions at the different stages. CircRNAs showed complexities and diversity of formation as the same host gene produces multiple isoforms of these nucleic acids with different expression profiles. The differential expression of 2881 circRNAs (DECs, P < 0.05) was identified and four were randomly selected and validated by qPCR. Moreover, 1118 DECs under strict selected (SDECs, |log2^FC| > 2 and P-adj value < 0.01) showed 4 expression trends (Clusters 0, 19, 16 and 18). Cluster 0 molecules had increasing expression at all stages with effects on muscle through metabolism, regulation of enzyme activity, and biosynthesis. Cluster 16 circRNAs had high expression in the early and late stages and are involved in “Wnt signaling pathway”, “AMPK signaling pathway” and others. Cluster 18 molecules were mainly expressed at F120 and participate in “cytoskeletal protein binding”, “Notch signaling pathway” and so on. Cluster 19 circRNAs were down-regulated at all stages and related to muscle structure and development. Lastly, the SDECs divided the period of skeletal muscle development into three transitional stages: stage 1 (F45 to F90), which related to muscle satellite cell proliferation and muscle fiber structure; stage 2 (F90 to B1), in which the attachment of the cytoplasmic surface to the actin cytoskeleton initiates; and stage 3, which involved the “cGMP-PKG signaling pathway”. Moreover, the paraffin sections messages also validated that there are three transitional stages of skeletal muscle development. Conclusion: Our current study provides a catalog of goat muscle-related circRNAs that can stratify skeletal muscle development fetus 45 days to newborn 90 days into three developmental stages. These findings better our understanding of functional transitions during mammalian muscle development.