Heat stress has a severe effect on animal health and can reduce the productivity and reproductive efficiency; it is therefore necessary to explore the molecular mechanism involved in heat stress response, which is helpful for the cultivation of an animal breed with resistance to heat stress. However, little research about heat stress-responsive molecular analysis has been reported in sheep. Therefore, in this study, RNA sequencing (RNA-Seq) was used to investigate the transcriptome profiling in the liver of Hu sheep with and without heat stress. In total, we detected 520 and 22 differentially expressed mRNAs and lncRNAs, respectively. The differentially expressed mRNAs were mainly associated with metabolic processes, the regulation of biosynthetic processes, and the regulation of glucocorticoid; additionally, they were significantly enriched in the heat stress related pathways, including the carbon metabolism, the PPAR signaling pathway, and vitamin digestion and absorption. The co-located differentially expressed lncRNA Lnc_001782 might positively influence the expression of the corresponding genes APOA4 and APOA5, exerting co-regulative effects on the liver function. Thus, we made the hypothesis that Lnc_001782, APOA4 and APOA5 might function synergistically to regulate the anti-heat stress ability in Hu sheep. This study provides a catalog of Hu sheep liver mRNAs and lncRNAs, and will contribute to a better understanding of the molecular mechanism underlying heat stress responses.

Background: Fertility is an important economic trait in the production of meat goat, and follicular development plays an important role in fertility. Although many mRNAs and microRNAs (miRNAs) have been found to play critical roles in ovarian biological processes, the interaction between mRNAs and miRNAs in follicular development is not yet completely understood. In addition, less attention has been given to the study of single follicle (dominant or atretic follicle) in goats. This study aimed to identify mRNAs, miRNAs, and signaling pathways as well as their interaction networks in the ovarian follicles (large follicles and small follicles) of uniparous and multiple Chuanzhong black goats at estrus phase using RNA-sequencing (RNA-seq) technique. Results: The results showed that there was a significant difference in the number of large follicles between uniparous and multiple goats (P < 0.05), but no difference in the number of small follicles was observed (P > 0.05). For the small follicles of uniparous and multiple goats at estrus phase, 289 differentially expressed mRNAs (DEmRNAs) and 16 DEmiRNAs were identified; and for the large follicles, 195 DEmRNAs and 7 DEmiRNAs were identified. The functional enrichment analysis showed that DE genes in small follicles were significantly enriched in ovarian steroidogenesis and steroid hormone biosynthesis, while in large follicles were significantly enriched in ABC transporters and steroid hormone biosynthesis. The results of quantitative real-time polymerase chain reaction were consistent with those of RNA-seq. Analysis of the mRNA-miRNA interaction network suggested that CD36 (miR-122, miR-200a, miR-141), TNFAIP6 (miR-141, miR-200a, miR-182), CYP11A1 (miR-122), SERPINA5 (miR-1, miR-206, miR-133a-3p, miR-133b), and PTGFR (miR-182, miR-122) might be related to fertility, but requires further research on follicular somatic cells. Conclusions: This study was used for the first time to reveal the DEmRNAs and DEmiRNAs as well as their interaction in the follicles of uniparous and multiple goats at estrus phase using RNA-seq technology. Our findings provide new clues to uncover the molecular mechanisms and signaling networks of goat reproduction that could be potentially used to increase ovulation rate and kidding rate in goat.

Circular RNAs (circRNAs) play important regulatory roles in a variety of biological processes in mammals. Multiple birth-traits in goats are affected by several factors, but the expression and function of circRNAs in follicular development of goats are not clear. In this study, we aimed to investigate the possible regulatory mechanisms of circRNA and collected five groups of large follicles (Follicle diameter > 6 mm) and small follicles (1 mm < Follicle diameter < 3 mm) from Leizhou goats in estrus for RNA sequencing. RNA sequencing showed that 152 circRNAs were differentially expressed in small and large follicles. Among them, 101 circRNAs were up-regulated in large follicles and 51 circRNAs were up-regulated in small follicles. GO and KEGG enrichment analyses showed that parental genes of the differential circRNAs were significantly enriched in important pathways, such as ovarian steroidogenesis, GnRH signaling pathway, animal autophagy and oxytocin signalling pathway. BioSignal analysis revealed that 152 differentially expressed circRNAs could target 91 differential miRNAs including miR-101 family (chi-miR-101-3p, chi-miR-101-5p), miR-202 family (chi-miR-202-5p, chi-miR-202-3p),60 circRNAs with translation potential. Based on the predicted sequencing results, the ceRNA networks chicirc_008762/chi-miR-338-3p/ARHGAP18 and chicirc_040444/chi-miR-338-3p/STAR were constructed in this study. Importantly, the new gene circCFAP20DC was first discovered in goats. The EDU assay and flow cytometry results indicated that circCFAP20DC enhanced the proliferation of follicular granulosa cells(GCs). Real-time quantitative PCR and western blotting assays showed that circCFAP20DC activated the Retinoblastoma(RB) pathway and promoted the progression of granulosa cells from G1 to S phase. Differential circRNAs in goat size follicles may have important biological functions for follicular development. The novel gene circCFAP20DC activates the RB pathway, promoting the progression of GCs from G1 to S phase. This, in turn, enhances the proliferation of follicular GCs in goats.

The present study aimed to investigate the colonization process of epithelial bacteria attached to the rumen and intestinal tract tissue during the development of goats after birth. However, this process from fetus to adulthood was very limited. In goats, the rumen and cecum are two important fermentation organs, and it is important to study the acquisition and development of the neonatal microbiome, as well as the difference between these two organs. To characterize the microbial establishment and dynamic changes in the rumen and cecum from fetus to adulthood, we performed 16S rRNA gene sequencing for 106 samples from 47 individuals of nine pregnant mother–fetus pairs and 16 kids from birth up to 6 months. The diversity, structure and composition of the microbial communities were distinct between the rumen and cecum after birth, while they were similar in the fetal period. The study showed a rapid loss and influx of microbes at birth, followed by slight selection after drinking colostrum, and then a strong selection after weaning, suggesting that the establishment and dynamic fluctuations of the gut microbiome undergoes three distinct phases of microbiome progression in life: a conserved phase (during late pregnancy in the fetus), a transitional phase (newborn until weaning), and a stable phase (from weaning to adulthood). The results supported the view that microbes exist in the fetus, and revealed the establishment and dynamic fluctuations of the gut microbiome from fetus to adulthood in goats.

Chinese local cattle with a high level of genetic diversity mainly originate from two subspecies; the cattle in northern China are primarilyBos Taurus, and the cattle in southern China are primarilyBos indicus. Cattle from southern China are characterized by a specific phenotype and adapted to the local environment. This study explored the genetic diversity, degree of admixture, and selection signature in eight local cattle breeds in southern China. The lowest level of heterozygosity was found in Hainan and Nandan cattle from Hainan and Guangxi province, respectively, whereas the highest level of heterozygosity was detected in Zhaotong cattle from Yunnan province. A neighbor-joining phylogenetic tree analysis clearly separated Lufeng cattle from other breeds, whereas Leiqiong and Hainan cattle have some crossover. Based on linkage disequilibrium-filtered single nucleotide polymorphisms (SNPs), the admixture analysis revealed two clusters corresponding to the taurine and indicine cattle lineages, and the local cattle breeds from southern China showed a certain degree of admixture. WhenK= 4 and 9, we found a slight separation among Leiqiong, Lufeng, and Hainan cattle. Meanwhile, we performed a selection signature analysis in Hainan, Leiqiong, and Lufeng cattle distributed in the extreme south of China, using the integrated haplotype score (iHS), Rsb statistic, and BayeScan software. Using the iHS approach, we identified 251, 270, and 256 candidate regions in Lufeng, Leiqiong, and Hainan cattle, respectively. Moreover, we identified 184, 174, and 146 candidate regions in pairwise comparisons of Leiqiong vs. Lufeng, Leiqiong vs. Hainan, and Hainan vs. Lufeng cattle using the Rsb approach. In addition, we identified 76 loci with a total of 48 genes under selection, based on the F_STapproach. Several candidate genes under selection were found to be related to meat quality, immunity, and adaptation to the local environment in southern China. Our results provide significant information about the genetic differences among the cattle breeds from southern China and the possible cause of difference in breed-specific characteristics. Selection signature analysis identified a few candidate SNPs and genes related to certain important traits of these cattle. In general, our results provide valuable insights into the genetic basis of specific traits under selection in certain local cattle breeds.