Meta-analysis of genome-wide association studies for litter size in sheep

Detection of polymorphisms in six genes and their association analysis with litter size in sheep

Litter size in sheep is a complex trait controlled by micro-effective polygenes. APAF1, CLSTN2, CTH, PLCB1, PLCB4, and CHST11 are all involved in mammalian reproduction. However, the effects of these genes on litter size in sheep are still unclear. Therefore, in this study, we used Sequenom MassARRAY® SNP assay technology to type the single nucleotide polymorphisms (SNPs) loci of six genes in five sheep breeds. The results showed that most sheep breeds contain three genotypes at each locus. Then, we conducted population genetic analysis on the SNPs of six genes and found that the polymorphic information content in all sheep breeds ranged from 0 to 0.37, and most sheep breeds were in Hardy-Weinberg equilibrium (p > 0.05). In addition, association analysis in Small Tail Han sheep indicated that the rs399534524 locus in CLSTN2 was highly associated with first parity litter size, and litter size in ewes with CT genotype was higher than that in ewes with CC genotype or TT genotype. Furthermore, the rs407142552 locus in CTH was highly associated with second parity litter size in Small Tail Han sheep, and litter size in ewes with CT genotype was higher than that in ewes with TT genotype. Finally, we predicted the CTH and CLSTN2 protein interaction network and found that HTR1E, NOM1, CCDC174 and ALPK3 interact with CLSTN2 and have been reported as candidate genes related to litter size in sheep. These results suggest that they may be useful genetic markers for increasing litter size in sheep.

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.

Multiomics Analyses Provide New Insight into Genetic Variation of Reproductive Adaptability in Tibetan Sheep

Domestication and artificial selection during production-oriented breeding have greatly shaped the level of genomic variability in sheep. However, the genetic variation associated with increased reproduction remains elusive. Here, two groups of samples from consecutively monotocous and polytocous sheep were collected for genome-wide association, transcriptomic, proteomic, and metabolomic analyses to explore the genetic variation in fecundity in Tibetan sheep. Genome-wide association study revealed strong associations between BMPR1B (p.Q249R) and litter size, as well as between PAPPA and lambing interval; these findings were validated in 1,130 individuals. Furthermore, we constructed the first single-cell atlas of Tibetan sheep ovary tissues and identified a specific mural granulosa cell subtype with PAPPA-specific expression and differential expression of BMPR1B between the two groups. Bulk RNA-seq indicated that BMPR1B and PAPPA expressions were similar between the two groups of sheep. 3D protein structure prediction and coimmunoprecipitation analysis indicated that mutation and mutually exclusive exons of BMPR1B are the main mechanisms for prolific Tibetan sheep. We propose that PAPPA is a key gene for stimulating ovarian follicular growth and development, and steroidogenesis. Our work reveals the genetic variation in reproductive performance in Tibetan sheep, providing insights and valuable genetic resources for the discovery of genes and regulatory mechanisms that improve reproductive success.

Investigation of SNP markers for the melatonin production trait in the Hu sheep with bulked segregant analysis

BACKGROUND

As an important reproductive hormone, melatonin plays an important role in regulating the reproductive activities of sheep and other mammals. Hu sheep is a breed favoring for meat, with prolific traits. In order to explore the relationship between melatonin and reproductive function of Hu sheep, 7,694,759 SNPs were screened out through the whole genome sequencing analysis from high and low melatonin production Hu sheep.

RESULTS

A total of 68,673 SNPs, involving in 1126 genes, were identified by ED association analysis. Correlation analysis of SNPs of AANAT/ASMT gene and MTNR1A/MTNR1B gene were carried out. The melatonin level of CG genotype 7,981,372 of AANAT, GA genotype 7,981,866 of ASMT and GG genotype 17,355,171 of MTNR1A were higher than the average melatonin level of 1.64 ng/mL. High melatonin Hu sheep appear to have better multiple reproductive performance.

CONCLUSIONS

By using different methods, three SNPs which are associated with high melatonin production trait have been identified in Hu sheep. These 3 SNPs are located in melatonin synthetase AANAT/ASMT and receptor MTNR1A, respectively. Considering the positive association between melatonin production and reproductive performance in ruminants, these three SNPs can be served as the potential molecular markers for breading Hu sheep with the desirable reproductive traits.

Comparative Transcriptomics Identify Key Pituitary Circular RNAs That Participate in Sheep (Ovis aries) Reproduction

CircRNAs have been found to play key roles in many biological processes and have diverse biological functions. There have been studies on circRNAs in sheep pituitary, and some important circRNAs have been found. But there are still few studies on circRNAs in sheep pituitary with different fecundity. In this study, we obtained the circRNAs expression profiles in the pituitary of FecB ++ genotype Small Tail Han sheep with different fecundity and estrous phases. A total of 34,878 circRNAs were identified in 12 pituitary samples, 300 differentially expressed circRNAs (DE circRNAs) (down: 104; up: 196) were identified in polytocous sheep in the follicular phase (PF) and monotocous sheep in the follicular phase (MF) (PF vs. MF), and 347 DE circRNAs (down: 162; up: 185) were identified in polytocous sheep in the luteal phase (PL) and monotocous sheep in the luteal phase (ML) (PL vs. ML). Cortisol synthesis and secretion pathway (follicular phase) and estrogen signaling pathway (luteal phase) were obtained by functional enrichment analysis of circRNAs source genes. Competing endogenous RNA (ceRNA) network analysis of key DE circRNAs revealed that oar-circ-0022776 (source gene ITPR2, follicular phase) targeted oar-miR-432, oar-circ-0009003 (source gene ITPR1, luteal phase) and oar-circ-0003113 (source gene PLCB1, luteal phase) targeted oar-miR-370-3p. We also explored the coding ability of DE circRNAs. In conclusion, our study shows that changes in the pituitary circRNAs may be related to the response of the pituitary to steroid hormones and regulate the reproductive process of sheep by affecting the pituitary function. Results of this study provide some new information for understanding the functions of circRNAs and the fecundity of FecB ++ genotype sheep.

Identification and characterization of circRNAs in peri-implantation endometrium between Yorkshire and Erhualian pigs

BACKGROUND

One of the most critical periods for the loss of pig embryos is the 12th day of gestation when implantation begins. Recent studies have shown that non-coding RNAs (ncRNAs) play important regulatory roles during pregnancy. Circular RNAs (circRNAs) are a kind of ubiquitously expressed ncRNAs that can directly regulate the binding proteins or regulate the expression of target genes by adsorbing micro RNAs (miRNA).

RESULTS

We used the Illumina Novaseq6,000 technology to analyze the circRNA expression profile in the endometrium of three Erhualian (EH12) and three Yorkshire (YK12) pigs on day 12 of gestation. Overall, a total of 22,108 circRNAs were identified. Of these, 4051 circRNAs were specific to EH12 and 5889 circRNAs were specific to YK12, indicating a high level of breed specificity. Further analysis showed that there were 641 significant differentially expressed circRNAs (SDEcircRNAs) in EH12 compared with YK12 (FDR < 0.05). Functional enrichment of differential circRNA host genes revealed many pathways and genes associated with reproduction and regulation of embryo development. Network analysis of circRNA-miRNA interactions further supported the idea that circRNAs act as sponges for miRNAs to regulate gene expression. The prediction of differential circRNA binding proteins further explored the potential regulatory pathways of circRNAs. Analysis of SDEcircRNAs suggested a possible reason for the difference in embryo survival between the two breeds at the peri-implantation stage.

CONCLUSIONS

Together, these data suggest that circRNAs are abundantly expressed in the endometrium during the peri-implantation period in pigs and are important regulators of related genes. The results of this study will help to further understand the differences in molecular pathways between the two breeds during the critical implantation period of pregnancy, and will help to provide insight into the molecular mechanisms that contribute to the establishment of pregnancy and embryo loss in pigs.

Association of Novel C319T Variant of PITX2 Gene 3'UTR Region With Reproductive Performance in Awassi Sheep

Several genes influence sheep's reproductive performance, among them the paired-like homeodomain transcription factor 2 (PITX2) gene. Thus, this study aimed to examine whether the variability within the PITX2 gene is associated with the reproductive performance of Awassi ewes. A total of 123 single-progeny ewes and 109 twin ewes were used to extract genomic DNA. An amplicon of 4 sequence fragments from exons 2, 4, 5 (upstream portion), and 5 (downstream portion) of the PITX2 gene was generated by polymerase chain reaction (PCR), 228, 304, 381, and 382 bp, respectively. Three genotypes of 382 bp amplicons were identified: CC, CT, and TT. Sequence analysis revealed a novel mutation in the CT genotype 319C > T. Statistical analysis revealed that single-nucleotide polymorphism (SNP) 319C > T was associated with reproductive performance. Single-nucleotide polymorphism 319C > T-carrying ewes had significantly (P ⩽ .01) lower litter sizes, twinning rates, lambing rates, and more days to lambing than those carrying CT and CC genotypes. Based on a logistic regression analysis, it was confirmed that the 319C > T SNP decreased litter size. Ewes with TT genotype produced fewer lambs than ewes with CT and CC genotypes. According to these results, the variant 319C> T SNP negatively affects the reproductive performance of Awassi sheep. Ewes carrying the 319C > T SNP have a lower litter size and are less prolific than those without the SNP.

A Genome-Wide Association Study of the Chest Circumference Trait in Xinjiang Donkeys Based on Whole-Genome Sequencing Technology

Animal genotyping by means of genome-wide association studies is important for connecting phenotypes of interest with their underlying genetics in livestock. However, the use of whole genome sequencing to investigate chest circumference (CC) in donkeys has rarely been reported. We aimed to use the genome-wide association study approach to detect significant single nucleotide polymorphisms (SNPs) and key genes associated with chest circumference traits in Xinjiang donkeys. We assessed 112 Xinjiang donkeys in this study. The chest circumference of each was measured 2 h before milking. We re-sequenced blood samples from the Xinjiang donkeys, and genome-wide association study analyses were performed using a mixed model with the PLINK, GEMMA, and REGENIE programs. We tested 38 donkeys for candidate SNPs for genome-wide association study using three software programs. Additionally, 18 SNP markers reached genome-wide significance (p < 1.61 × 10^-9). On the basis of these, 41 genes were identified. Previously proposed candidate genes for CC traits were supported by this study, including NFATC2 (Nuclear Factor of Activated T Cells 2), PROP1 (PROP Paired-Like Homeobox 1), UBB (Ubiquitin B), and HAND2 (Heart and Neural Crest Derivatives Expressed 2). These promising candidates provide a valuable resource for validating potential meat production genes and will facilitate the development of high-yielding Xinjiang donkey breeds through marker-assisted selection or gene editing.

Detection of Candidate Genes Associated with Fecundity through Genome-Wide Selection Signatures of Katahdin Ewes

One of the strategies to genetically improve reproductive traits, despite their low inheritability, has been the identification of candidate genes. Therefore, the objective of this study was to detect candidate genes associated with fecundity through the fixation index (F_ST) and runs of homozygosity (ROH) of selection signatures in Katahdin ewes. Productive and reproductive records from three years were used and the genotypes (OvineSNP50K) of 48 Katahdin ewes. Two groups of ewes were identified to carry out the genetic comparison: with high fecundity (1.3 ± 0.03) and with low fecundity (1.1 ± 0.06). This study shows for the first time evidence of the influence of the CNOT11, GLUD1, GRID1, MAPK8, and CCL28 genes in the fecundity of Katahdin ewes; in addition, new candidate genes were detected for fecundity that were not reported previously in ewes but that were detected for other species: ANK2 (sow), ARHGAP22 (cow and buffalo cow), GHITM (cow), HERC6 (cow), DPF2 (cow), and TRNAC-GCA (buffalo cow, bull). These new candidate genes in ewes seem to have a high expression in reproduction. Therefore, future studies are needed focused on describing the physiological basis of changes in the reproductive behavior influenced by these genes.

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.