Additional file 19: Supplementary material 9A. Co-expression details of DE lncRNA-mRNA after lncRNA targets coincided with DE mRNAs in MM_FP vs MM_LP. Supplementary material 9B. Co-expression details of DE lncRNA-mRNA after lncRNA targets coincided with DE mRNAs in MM_FP vs ww_FP. Supplementary material 9C. Co-expression details of DE lncRNA-mRNA after lncRNA targets coincided with DE mRNAs in MM_LP vs ww_LP. Supplementary material 9D. Co-expression details of DE lncRNA-mRNA after lncRNA targets coincided with DE mRNAs in ww_FP vs ww_LP.
In previous studies, NOX4, PDE11A and GHR genes have been screened as important candidate genes for litter size in sheep by using the GWAS method; however, neither their effects on litter size nor the loci associated with litter size have been identified. In this study, three candidate loci (c.1057-4C > T in NOX4, c.1983C > T in PDE11A and c.1618C > T in GHR) were first screened based on our previous resequencing data of 10 sheep breeds. After the three loci were genotyped using Sequenom MassARRAY technology, we carried out population genetics analysis on the three loci and performed association analysis between the polymorphism of the three loci and the litter size of sheep. The results of population genetics analysis suggested that c.1057-4C > T in NOX4 and c.1983C > T in PDE11A may be subject to natural or artificial selection. The results of association analysis indicated that litter size was significantly associated with c.1057-4C > T in NOX4 and c.1983C > T in PDE11A (p < 0.05) in Small Tail Han sheep, and there was no significant interaction effect between the two loci on the litter size. In summary, c.1057-4C > T in NOX4 and c.1983 C > T in PDE11A can be considered candidate molecular markers for improving litter size in sheep.
A normal estrus cycle is essential for the breeding of goats, and the luteal phase accounts for most of the estrus cycle. The corpus luteum (CL) formed during the luteal phase is a transient endocrine gland that is crucial for the reproductive cycle and pregnancy maintenance, and is controlled by many regulatory factors. However, the molecular mechanism of the hypothalamus effect on the reproductive performance of different litter sizes during the luteal phase of goats has not been elucidated. In this study, RNA-sequencing was used to analyze the mRNA and miRNA expression profiles of the hypothalamic tissues with the high-fecundity goats during the luteal phase (LP-HF) and low-fecundity goats during the luteal phase (LP-LF). The RNA-seq results found that there were 1963 differentially expressed genes (DEGs) (890 up-regulated and 1073 down-regulated). The miRNA-seq identified 57 differentially expressed miRNAs (DEMs), including 11 up-regulated and 46 down-regulated, of which 199 DEGs were predicted to be potential target genes of DEMs. Meanwhile, the functional enrichment analysis identified several mRNA-miRNA pairs involved in the regulation of the hypothalamic activity, such as the common target gene MEA1 of novel-miR-972, novel-miR-125 and novel-miR-403, which can play a certain role as a related gene of the reproductive development in the hypothalamic-pituitary-gonadal (HPG) axis and its regulated network, by regulating the androgen secretion. While another target gene ADIPOR2 of the novel-miR-403, is distributed in the hypothalamus and affects the reproductive system through a central role on the HPG axis and a peripheral role in the gonadal tissue. An annotation analysis of the DE miRNA-mRNA pairs identified targets related to biological processes, such as anion binding (GO:0043168) and small molecule binding (GO: 0036094). Subsequently, the KEGG(Kyoto Encyclopedia of Genes and Genomes) pathways were performed to analyze the miRNA-mRNA pairs with negatively correlated miRNAs. We found that the GnRH signaling pathway (ko04912), the estrogen signaling pathway (ko04915), the Fc gamma R-mediated phagocytosis (ko04666), and the IL-17 signaling pathway (ko04657), etc., were directly and indirectly associated with the reproductive process. These targeting interactions may be closely related to the reproductive performance of goats. The results of this study provide a reference for further research on the molecular regulation mechanism for the high fertility in goats.
Abstract Prolificacy of most local goat breeds in China is low. Jining Grey goat is one of the most prolific goat breeds in China, it is an important goat breed for the rural economy. ASMT (acetylserotonin O‐methyltransferase) and ADAMTS1 (ADAM metallopeptidase with thrombospondin type 1 motif) are essential for animal reproduction. Single nucleotide polymorphisms (SNPs) of ASMT and ADAMTS1 genes in the highly prolific breed (Jining Grey goats), medium prolific breed (Boer goats and Guizhou White goats) and low prolific breeds (Angora goats, Liaoning Cashmere goats and Inner Mongolia Cashmere goats) were detected by polymerase chain reaction‐restriction fragment length polymorphism and sequencing. Two SNPs (g.158122T>C, g.158700G>A) of ASMT gene and two SNPs (g.7979798A>G, g.7979477C>T) of ADAMTS1 gene were identified. For g.158122T>C of ASMT gene, further analysis revealed that genotype TC or CC had 0.66 ( p < 0.05) or 0.75 ( p < 0.05) kids more than those with genotype TT in Jining Grey goats. No significant difference ( p > 0.05) was found in litter size between TC and CC genotypes. The SNP (g.158122T>C) caused a p.Tyr298His change and this SNP mutation resulted in changes in protein binding sites and macromolecule‐binding sites. The improvement in reproductive performance may be due to changes in the structure of ASMT protein. For g.7979477C>T of ADAMTS1 gene, Jining Grey does with genotype CT or TT had 0.82 ( p < 0.05) or 0.86 ( p < 0.05) more kids than those with genotype CC. No significant difference ( p > 0.05) was found in litter size between CT or TT genotypes. These results preliminarily indicated that C allele (g.158122T>C) of ASMT gene and T allele (g.7979477C>T) of ADAMTS1 gene are potential molecular markers which could improve litter size of Jining Grey goats and be used in goat breeding.
With the intensified and large-scale development of sheep husbandry and global warming, sheep heat stress has become an increasingly important issue. However, little is known about the molecular mechanisms related to sheep responses to heat stress. In this study, transcriptomic analysis of liver tissues of sheep in the presence and absence of heat stress was conducted, with the goal of identifying genes and pathways related to regulation when under such stress. After a comparison with the sheep reference genome, 440,226,436 clean reads were obtained from eight libraries. A p-value ≤ 0.05 and fold change ≥ 2 were taken as thresholds for categorizing differentially expressed genes, of which 1137 were identified. The accuracy and reliability of the RNA-Seq results were confirmed by qRT-PCR. The identified differentially expressed genes were significantly associated with 419 GO terms and 51 KEGG pathways, which suggested their participation in biological processes such as response to stress, immunoreaction, and fat metabolism. This study’s results provide a comprehensive overview of sheep heat stress-induced transcriptional expression patterns, laying a foundation for further analysis of the molecular mechanisms of sheep heat stress.
The thyroid gland is a vital endocrine organ involved in the reproduction of animals via the regulation of hormone synthesis and secretion. LncRNAs have been proven to play important roles in reproductive regulation; however, the associated mechanism in the thyroid gland has not been clarified. In this study, we investigated to identify photoperiod-induced lncRNAs and mRNAs in the thyroid gland in Sunite ewes by comparing the expression profiles of short photoperiod (SP) and long photoperiods (LP). A total of 41,088 lncRNAs were identified in the thyroid gland through RNA-Seq. Functional analysis of differentially expressed lncRNAs using the R package revealed that reproductive hormone- and photoperiod response-related pathways, including the prolactin signaling, cAMP signaling, and circadian rhythm pathways, were significantly enriched. An mRNA-lncRNA interaction analysis suggested that the lncRNA LOC1056153S88 trans targets ARG2 and CCNB3, and the lncRNA LOC105607004 trans targets DMXL2, both of these might be involved in seasonal sheep breeding reproduction. Together, these results will provide resources for further studies on seasonal reproduction in sheep.
Alpha 1-antitrypsin (α1-AT) may affect the susceptibility of mastitis in dairy cattle for its possible role in the protection of lactoferrin from proteolytic degradation in the mammary. Milk somatic cell score (SCS) is a logarithmic transformation of the milk somatic cell count widely used as an index to evaluate mastitis. To study the relationships of α1-AT gene and SCS in Chinese Holstein cows, methods of PCR-SSCP, DNA sequencing, PCR-RFLP, and CRS-PCR technologies were used to detect single nucleotide polymorphisms sites in α1-AT gene. Two polymorphic sites at G5503A and G5746C of α1-AT gene were found. AA (0.3633), AB (0.4644) and BB (0.1723) genotypes were detected at G5503A site, CC (0.3483), CD (0.4906) and DD (0.1611) genotypes were found at G5746C in Chinese Holstein cows. Least squares mean of SCS for individuals with BB genotype was significantly lower than that with AA and AB genotype (p < 0.01), and that with AB genotype was significantly lower than that with AA (p < 0.05). There was no significant difference among individuals with CC, CD and DD genotypes (p > 0.05). Least squares mean of SCS for individuals with BBDD genotype combination were significantly lower than those with AACC and AACD (p < 0.05). Statistical analysis indicated that B allele and BBDD genotype combination of α1-AT can improve mastitis resistance in dairy cattle.