Screen of small fragment mutations within the sheep thyroid stimulating hormone receptor gene associated with litter size

Comparative whole-genome resequencing to uncover selection signatures linked to litter size in Hu Sheep and five other breeds

Hu sheep (HS), a breed of sheep carrying the FecB mutation gene, is known for its "year-round estrus and multiple births" and is an ideal model for studying the high fecundity mechanisms of livestock. Through analyzing and comparing the genomic selection features of Hu sheep and other sheep breeds, we identified a series of candidate genes that may play a role in Hu sheep's high fecundity mechanisms. In this study, we conducted whole-genome resequencing on six breeds and screened key mutations significantly correlated with high reproductive traits in sheep. Notably, the CC2D1B gene was selected by the fixation index (FST) and the cross-population composite likelihood ratio (XP-CLR) methods in HS and other five breeds. It was worth noting that the CC2D1B gene in HS was different from that in other sheep breeds, and seven missense mutations have been identified. Furthermore, the linkage disequilibrium (LD) analysis revealed a strong linkage disequilibrium in this specific gene region. Subsequently, by performing different grouping based on FecB genotypes in Hu sheep, genome-wide selective signal analysis screened several genes related to reproduction, such as BMPR1B and PPM1K. Besides, FST analysis identified functional genes related to reproductive traits, including RHEB, HSPA2, PPP1CC, HVCN1, and CCDC63. Additionally, a missense mutation was found in the CCDC63 gene and the haplotype was different between the high reproduction (HR) group and low reproduction (LR) group in HS. In summary, we discovered genetic differentiation among six distinct breeding sheep breeds at the whole genome level. Additionally, we identified a set of genes which were associated with reproductive performance in Hu sheep and visualized how these genes differed in different breeds. These findings laid a theoretical foundation for understanding genetic mechanisms behind high prolific traits in sheep.

Genetic diversity, population structure, and selective signature of sheep in the northeastern Tarim Basin

Local sheep in the northeastern Tarim Basin can adapt to dry and low-rainfall regional environments. In this study, three local sheep breeds in the northeastern Tarim Basin, LOP (LOP) sheep, Bayinbuluke (BYK) sheep, and Kunlun (KUN, also known as the Qiemo sheep) sheep, and three introduced sheep breeds, Suffolk (SUF) sheep, Dorset (APD) sheep, and Texel (TEX) sheep, were analyzed for genetic diversity, population structure, and selective signature using the Illumina OvineSNP50K BeadChip. We found that LOP, BYK, and KUN had lower observed heterozygosity and expected heterozygosity than TEX, SUF, and ADP, which were differentiated based on geographic distribution. We performed fixation index (FST) analysis on three local sheep breeds in the northeastern Tarim Basin (LOP, BYK, and KUN) and introduced sheep breeds (TEX, SUF, and ADP) to measure genetic differentiation. Nucleotide diversity (PI) analysis was performed on single-nucleotide polymorphism (SNP) data of LOP, BYK, and KUN. A total of 493 candidate genes were obtained by taking the intersection at a threshold of 5%. Among them, SMAD2, ESR2, and HAS2 were related to reproductive traits. PCDH15, TLE4, and TFAP2B were related to growth traits. SOD1, TSHR, and DNAJB5 were related to desert environmental adaptation. Analyzing the genetic patterns of local sheep in the northeastern Tarim Basin can protect the germplasm resources of local sheep and promote the development and utilization of sheep genetic resources.

Whole-genome resequencing reveals genetic diversity and selection characteristics of dairy goat

The dairy goat is one of the earliest dairy livestock species, which plays an important role in the economic development, especially for developing countries. With the development of agricultural civilization, dairy goats have been widely distributed across the world. However, few studies have been conducted on the specific characteristics of dairy goat. In this study, we collected the whole-genome data of 89 goat individuals by sequencing 48 goats and employing 41 publicly available goats, including five dairy goat breeds (Saanen, Nubian, Alpine, Toggenburg, and Guanzhong dairy goat; n = 24, 15, 11, 6, 6), and three goat breeds (Guishan goat, Longlin goat, Yunshang Black goat; n = 6, 15, 6). Through compared the genomes of dairy goat and non-dairy goat to analyze genetic diversity and selection characteristics of dairy goat. The results show that the eight goats could be divided into three subgroups of European, African, and Chinese indigenous goat populations, and we also found that Australian Nubian, Toggenburg, and Australian Alpine had the highest linkage disequilibrium, the lowest level of nucleotide diversity, and a higher inbreeding coefficient, indicating that they were strongly artificially selected. In addition, we identified several candidate genes related to the specificity of dairy goat, particularly genes associated with milk production traits (GHR, DGAT2, ELF5, GLYCAM1, ACSBG2, ACSS2), reproduction traits (TSHR, TSHB, PTGS2, ESR2), immunity traits (JAK1, POU2F2, LRRC66). Our results provide not only insights into the evolutionary history and breed characteristics of dairy goat, but also valuable information for the implementation and improvement of dairy goat cross breeding program.

Photoperiods induced the circRNA differential expression in the thyroid gland of OVX+E2 ewes

Circular RNAs (circRNAs) are non-coding RNAs newly identified and play important roles in RNA regulation. However, little is known regarding photoperiods induced circRNAs in the thyroid gland. In this study, we performed a comprehensive analysis of circRNA profiles in the thyroid gland of OVX+E₂ ewes at different photoperiods by whole transcriptome sequencing. A total of 37,470 novel circRNAs were detected in different photoperiods (42 days of short photoperiod treatment, SP42; 42 days of long photoperiod treatment, LP42; SP42 transfer to LP42, SPLP42), with a total of 817 circRNAs for SP42-LP42 (down: 132; up: 114), LP42-SPLP42 (down: 136; up: 112) and SP42-SPLP42 (down: 182; up: 141) having differentially expressed. Functional enrichment annotation analysis of DE-circRNAs for GO and KEGG by R package, features that influence photoperiod response in Sunite ewes through the Inositol phosphate metabolism, cGMP-PKG signaling pathway, Calcium signaling pathway, MAPK signaling pathway, and Oocyte meiosis. In addition, competitive endogenous RNA (ceRNA) network analysis revealed target binding sites for identified miRNAs in DE-cirRNAs such as oar-miR-10b, oar-miR-200c, oar-miR-21, oar-miR-370-3p, oar-miR-377-3p, oar-miR-181a, oar-miR-432, and oar-miR-495-3p. These results of this study will provide some new information for understanding circRNA function as well as the changes in the sheep thyroid gland under different photoperiods.