Genome-Wide Analyses Reveal Genetic Convergence of Prolificacy between Goats and 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.

Whole-genome sequencing identified candidate genes associated with high and low litter size in Chuanzhong black goats

The reproductive performance of goats significantly influences breeding efficiency and economic returns, with litter size serving as a comprehensive indicator. Despite this, research on the genetic control of litter size remains limited. Therefore, we aimed to explore the candidate genes affecting fecundity and compared the whole-genome sequences (WGS) of 15 high-litter (HL) and 15 low-litter (LL) size in Chuanzhong black goats. Then genetic diversity and genomic variation patterns were analyzed by phylogenetic, principal component and population genetic structure analysis, it was found that HL and LL subpopulations diverged. Population evolutionary selection elimination analysis was performed by Fst and θπ resulted in 506 genes were annotated in HL and 528 genes in LL. These genes were mainly related to Hippo signaling pathway, G protein-coupled signaling pathway, G protein-coupled receptor activity, cell surface receptor signaling pathway, gonadal and reproductive structure development. According to the significantly selected genomic regions and important pathways, we found that the g.89172108T > G variant locus in the exon of the AMH gene was significantly associated with litter size (P < 0.05), which could be used as an auxiliary selection gene for the high fertility of Chuanzhong black Goat.

Whole-genome resequencing reveals diversity and selective signals in the Wuxue goat

Animal genetic resources are crucial for ensuring global food security. However, in recent years, a noticeable decline in the genetic diversity of livestock has occurred worldwide. This decline is pronounced in developing countries, where the management of these resources is insufficient. In the current study, we performed whole genome sequencing for 20 Wuxue (WX) and five Guizhou White (GW) goats. Additionally, we utilized the published genomes of 131 samples representing five different goat breeds from various regions in China. We investigated and compared the genetic diversity and selection signatures of WX goats. Whole genome sequencing analysis of the WX and GW populations yielded 120 425 063 SNPs, which resided primarily in intergenic and intron regions. Population genetic structure revealed that WX exhibited genetic resemblance to GW, Chengdu Brown, and Jintang Black and significant differentiation from the other goat breeds. In addition, three methods (nucleotide diversity, linkage disequilibrium decay, and runs of homozygosity) showed moderate genetic diversity in WX goats. We used nucleotide diversity and composite likelihood ratio methods to identify within-breed signatures of positive selection in WX goats. A total of 369 genes were identified using both detection methods, including genes related to reproduction (GRID2, ZNF276, TCF25, and SPIRE2), growth (HMGA2 and GJA3), and immunity (IRF3 and SRSF3). Overall, this study explored the adaptability of WX goats, shedding light on their genetic richness and potential to thrive in challenges posed by climatic changes and diseases. Further investigations are warranted to harness these insights to enhance more efficient and sustainable goat breeding initiatives.

Genomic Regions Associated with Resistance to Gastrointestinal Parasites in Australian Merino Sheep

The objective of this study was to identify genomic regions and genes associated with resistance to gastrointestinal nematodes in Australian Merino sheep in Uruguay, using the single-step GWAS methodology (ssGWAS), which is based on genomic estimated breeding values (GEBVs) obtained from a combination of pedigree, genomic, and phenotypic data. This methodology converts GEBVs into SNP effects. The analysis included 26,638 animals with fecal egg count (FEC) records obtained in two independent parasitic cycles (FEC1 and FEC2) and 1700 50K SNP genotypes. The comparison of genomic regions was based on genetic variances (gVar(%)) explained by non-overlapping regions of 20 SNPs. For FEC1 and FEC2, 18 and 22 genomic windows exceeded the significance threshold (gVar(%) ≥ 0.22%), respectively. The genomic regions with strong associations with FEC1 were located on chromosomes OAR 2, 6, 11, 21, and 25, and for FEC2 on OAR 5, 6, and 11. The proportion of genetic variance attributed to the top windows was 0.83% and 1.9% for FEC1 and FEC2, respectively. The 33 candidate genes shared between the two traits were subjected to enrichment analysis, revealing a marked enrichment in biological processes related to immune system functions. These results contribute to the understanding of the genetics underlying gastrointestinal parasite resistance and its implications for other productive and welfare traits in animal breeding programs.

Structural variant landscapes reveal convergent signatures of evolution in sheep and goats

BACKGROUND

Sheep and goats have undergone domestication and improvement to produce similar phenotypes, which have been greatly impacted by structural variants (SVs). Here, we report a high-quality chromosome-level reference genome of Asiatic mouflon, and implement a comprehensive analysis of SVs in 897 genomes of worldwide wild and domestic populations of sheep and goats to reveal genetic signatures underlying convergent evolution.

RESULTS

We characterize the SV landscapes in terms of genetic diversity, chromosomal distribution and their links with genes, QTLs and transposable elements, and examine their impacts on regulatory elements. We identify several novel SVs and annotate corresponding genes (e.g., BMPR1B, BMPR2, RALYL, COL21A1, and LRP1B) associated with important production traits such as fertility, meat and milk production, and wool/hair fineness. We detect signatures of selection involving the parallel evolution of orthologous SV-associated genes during domestication, local environmental adaptation, and improvement. In particular, we find that fecundity traits experienced convergent selection targeting the gene BMPR1B, with the DEL00067921 deletion explaining ~10.4% of the phenotypic variation observed in goats.

CONCLUSIONS

Our results provide new insights into the convergent evolution of SVs and serve as a rich resource for the future improvement of sheep, goats, and related livestock.

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.

Genomic predictions and GWAS for heat tolerance in pigs based on reaction norm models with performance records and data from public weather stations considering alternative temperature thresholds

Genetic improvement of livestock productivity has resulted in greater production of metabolic heat and potentially greater susceptibility to heat stress. Various studies have demonstrated that there is genetic variability for heat tolerance and genetic selection for more heat tolerant individuals is possible. The rate of genetic progress tends to be greater when genomic information is incorporated into the analyses as more accurate breeding values can be obtained for young individuals. Therefore, this study aimed (1) to evaluate the predictive ability of genomic breeding values for heat tolerance based on routinely recorded traits, and (2) to investigate the genetic background of heat tolerance based on single-step genome-wide association studies for economically important traits related to body composition, growth and reproduction in Large White pigs. Pedigree information was available for 265,943 animals and genotypes for 8686 animals. The studied traits included ultrasound backfat thickness (BFT), ultrasound muscle depth (MDP), piglet weaning weight (WW), off-test weight (OTW), interval between farrowing (IBF), total number of piglets born (TNB), number of piglets born alive (NBA), number of piglets born dead (NBD), number of piglets weaned (WN) and weaning-to-estrus interval (IWE). The number of phenotypic records ranged from 6059 (WN) to 172,984 (TNB). Single-step genomic reaction norm predictions were used to calculate the genomic estimated breeding values for each individual. Predictions of breeding values for the validation population individuals were compared between datasets containing phenotypic records measured in the whole range of temperatures (WR) and datasets containing only phenotypic records measured when the weather station temperature was above 10°C (10C) or 15°C (15C), to evaluate the usefulness of these datasets that may better reflect the within-barn temperature. The use of homogeneous or heterogeneous residual variance was found to be trait-dependent, where homogeneous variance presented the best fit for MDP, BFT, OTW, TNB, NBA, WN and IBF, while the other traits (WW and IWE) had better fit with heterogeneous variance. The average prediction accuracy, dispersion and bias values considering all traits for WR were 0.36 ± 0.05, -0.07 ± 0.13 and 0.76 ± 0.10, respectively; for 10C were 0.39 ± 0.02, -0.05 ± 0.07 and 0.81 ± 0.05, respectively; and for 15C were 0.32 ± 0.05, -0.05 ± 0.11 and 0.84 ± 0.10, respectively. Based on the studied traits, using phenotypic records collected when the outside temperature (from public weather stations) was above 10°C provided better predictions for most of the traits. Forty-three and 62 candidate genomic regions were associated with the intercept (overall performance level) and slope term (specific biological mechanisms related to environmental sensitivity), respectively. Our results contribute to improve genomic predictions using existing datasets and better understand the genetic background of heat tolerance in pigs. Furthermore, the genomic regions and candidate genes identified will contribute to future genomic studies and breeding applications.

Identification of overlay differentially expressed genes in both rats and goats with blast lung injury through comparative transcriptomics

PURPOSE

To identify the potential target genes of blast lung injury (BLI) for the diagnosis and treatment.

METHODS

This is an experimental study. The BLI models in rats and goats were established by conducting a fuel-air explosive power test in an unobstructed environment, which was subsequently validated through hematoxylin-eosin staining. Transcriptome sequencing was performed on lung tissues from both goats and rats. Differentially expressed genes were identified using the criteria of q ≤ 0.05 and |log2 fold change| ≥ 1. Following that, enrichment analyses were conducted for gene ontology and the Kyoto Encyclopedia of Genes and Genomes pathways. The potential target genes were further confirmed through quantitative real-time polymerase chain reaction and enzyme linked immunosorbent assay.

RESULTS

Observations through microscopy unveiled the presence of reddish edema fluid, erythrocytes, and instances of focal or patchy bleeding within the alveolar cavity. Transcriptome sequencing analysis identified a total of 83 differentially expressed genes in both rats and goats. Notably, 49 genes exhibited a consistent expression pattern, with 38 genes displaying up-regulation and 11 genes demonstrating down-regulation. Enrichment analysis highlighted the potential involvement of the interleukin-17 signaling pathway and vascular smooth muscle contraction pathway in the underlying mechanism of BLI. Furthermore, the experimental findings in both goats and rats demonstrated a strong association between BLI and several key genes, including anterior gradient 2, ankyrin repeat domain 65, bactericidal/permeability-increasing fold containing family A member 1, bactericidal/permeability-increasing fold containing family B member 1, and keratin 4, which exhibited up-regulation.

CONCLUSIONS

Anterior gradient 2, ankyrin repeat domain 65, bactericidal/permeability-increasing fold containing family A member 1, bactericidal/permeability-increasing fold containing family B member 1, and keratin 4 hold potential as target genes for the prognosis, diagnosis, and treatment of BLI.

Identifying Candidate Genes for Litter Size and Three Morphological Traits in Youzhou Dark Goats Based on Genome-Wide SNP Markers

This study aimed to reveal the potential genetic basis for litter size, coat colour, black middorsal stripe and skin colour by combining genome-wide association analysis (GWAS) and selection signature analysis and ROH detection within the Youzhou dark (YZD) goat population (n = 206) using the Illumina GoatSNP54 BeadChip. In the GWAS, we identified one SNP (snp54094-scaffold824-899720) on chromosome 11 for litter size, two SNPs on chromosome 26 (snp11508-scaffold142-1990450, SORCS3) and chromosome 12 (snp55048-scaffold842-324525, LOC102187779) for coat colour and one SNP on chromosome 18 (snp56013-scaffold873-22716, TCF25) for the black middorsal stripe. In contrast, no SNPs were identified for skin colour. In selection signature analysis, 295 significant iHS genomic regions with a mean |iHS| score > 2.66, containing selection signatures encompassing 232 candidate genes were detected. In particular, 43 GO terms and one KEGG pathway were significantly enriched in the selected genes, which may contribute to the excellent environmental adaptability and characteristic trait formation during the domestication of YZD goats. In ROH detection, we identified 4446 ROH segments and 282 consensus ROH regions, among which nine common genes overlapped with those detected using the iHS method. Some known candidate genes for economic traits such as reproduction (TSHR, ANGPT4, CENPF, PIBF1, DACH1, DIS3, CHST1, COL4A1, PRKD1 and DNMT3B) and development and growth (TNPO2, IFT80, UCP2, UCP3, GHRHR, SIM1, CCM2L, CTNNA3 and CTNNA1) were revealed by iHS and ROH detection. Overall, this study is limited by the small population size, which affects the results of GWAS to a certain extent. Nevertheless, our findings could provide the first overview of the genetic mechanism underlying these important traits and provide novel insights into the future conservation and utilisation of Chinese goat germplasm resources.

Genomic Regions Associated with Wool, Growth and Reproduction Traits in Uruguayan Merino Sheep

The aim of this study was to identify genomic regions and genes associated with the fiber diameter (FD), clean fleece weight (CFW), live weight (LW), body condition score (BCS), pregnancy rate (PR) and lambing potential (LP) of Uruguayan Merino sheep. Phenotypic records of approximately 2000 mixed-age ewes were obtained from a Merino nucleus flock. Genome-wide association studies were performed utilizing single-step Bayesian analysis. For wool traits, a total of 35 genomic windows surpassed the significance threshold (PVE ≥ 0.25%). The proportion of the total additive genetic variance explained by those windows was 4.85 and 9.06% for FD and CFW, respectively. There were 42 windows significantly associated with LWM, which collectively explained 43.2% of the additive genetic variance. For BCS, 22 relevant windows accounted for more than 40% of the additive genetic variance, whereas for the reproduction traits, 53 genomic windows (24 and 29 for PR and LP, respectively) reached the suggestive threshold of 0.25% of the PVE. Within the top 10 windows for each trait, we identified several genes showing potential associations with the wool (e.g., IGF-1, TGFB2R, PRKCA), live weight (e.g., CAST, LAP3, MED28, HERC6), body condition score (e.g., CDH10, TMC2, SIRPA, CPXM1) or reproduction traits (e.g., ADCY1, LEPR, GHR, LPAR2) of the mixed-age ewes.

Identification of selection signatures in Capra hircus and Capra aegagrus in Iran

Identification of selection signatures may provide a better understanding of domestication process and candidate genes contributing to this process. In this study, two populations of domestic and wild goats from Iran were analyzed to identify selection signatures. RSB, iHS, and XP-EHH statistics were used in order to identify robust selection signatures in the goat genome. Genotype data of domestic and wild goats from the NextGen project was used. The data was related to 18 Capra aegagrus (wild goat) and 20 Capra hircus (domestic goat) from Iran. The iHS method indicated 675 and 441 selection signatures in C. aegagrus and C. hircus, respectively. RSB and XP-EHH methods showed about 370 and 447 selection signatures in C. aegagrus and C. hircus, respectively. These selection signatures were mainly associated with milk production, fleece trait, mammary epithelial cells, reproduction, and immune system.

Trajectory of livestock genomics in South Asia: A comprehensive review

Livestock plays a central role in sustaining human livelihood in South Asia. There are numerous and distinct livestock species in South Asian countries. Several of them have experienced genetic development in recent years due to the application of genomic technologies and effective breeding programs. This review discusses genomic studies on cattle, buffalo, sheep, goat, pig, horse, camel, yak, mithun, and poultry. The frontiers covered in this review are genetic diversity, admixture studies, selection signature research, QTL discovery, genome-wide association studies (GWAS), and genomic selection. The review concludes with recommendations for South Asian livestock systems to increasingly leverage genomic technologies, based on the lessons learned from the numerous case studies. This paper aims to present a comprehensive analysis of the dichotomy in the South Asian livestock sector and argues that a realistic approach to genomics in livestock can ensure long-term genetic advancements.

A Novel SNP in the Promoter Region of IGF1 Associated With Yunshang Black Goat Kidding Number via Promoting Transcription Activity by SP1

IGF1, a member of the insulin-like growth factor (IGF) superfamily, is also known as the growth-promoting factor (somatomedin C). IGF1 is involved in vertebrate growth and development, immunity, cell metabolism, reproduction, and breeding. However, there are relatively few studies on the relationship between IGF1 and goat reproduction. In this study, a new transcription factor SP1 bound to the IGF1 g. 64943050T>C promoted granulosa cell (GC) proliferation. A mutation g.64943050T>C located in the promoter region of IGF1 was identified. Association analysis revealed that the kidding number in the first and second litters and the average number of first three litters of the CC genotype (2.206 ± 0.044, 2.254 ± 0.056, and 2.251 ± 0.031) were significantly higher than those in the TC genotype (1.832 ± 0.049, 1.982 ± 0.06, and 1.921 ± 0.034) and TT genotype (1.860 ± 0.090, 1.968 ± 0.117, and 1.924 ± 0.062) (p < 0.05). The kidding number in the third litter of the CC genotype (2.355 ± 0.057) was significantly higher than that in the TT genotype (2.000 ± 0.107) (p < 0.05). Then, the function of this mutation was validated by the dual-luciferase reporter assay and EMSA. The results showed that the luciferase activity of IGF1-mutant-C was significantly higher than that of IGF1-Wild-T (p < 0.05). The EMSA also showed that the binding ability of IGF1-mutant-C was higher than that of IGF1-Wild-T (p < 0.05). Subsequently, the transcription factor SP1 was predicted to bind to the mutation of IGF1 (g.64943050T>C). Overexpression of SP1 promotes the expression of IGF1 in the primary granulosa cells (GCs). The results of the CCK-8 assay and the expression of GC proliferation factors (CDK4, cyclin D1, and cyclin D2) demonstrated that SP1 promoted GC proliferation by regulating IGF1 expression. Our results suggested that the IGF1 g.64943050T>C was significantly associated with the kidding number of Yunshang black goats, and SP1 as a transcription factor of IGF1 binding to the mutation T>C regulated the expression of IGF1. Furthermore, SP1 promoted goat GC proliferation by regulating the expression of IGF1, which provides a new insight for the goat fertility trait.

Genome-wide detection of selective signals for fecundity traits in goats (Capra hircus)

Fecundity in livestock is an economically important complex quantitative trait that is influenced by both genetics and the environment. However, the underlying genetic mechanism of reproductive performance in goats has not been well investigated. To investigate the genomic basis of fecundity in goats, genomic sequencing data of the Jining grey goat (a high prolificacy breed in China) were collected, as well as data for other commonly available goat breeds, and a mass of genomic variants were generated after variation calling. We screened the Jining grey goat (20 individuals) using a selective sweep with the Asian wild goat population (5 individuals), and potential candidate genes were proposed, such as STIM1, ESR1, LRRC14B and SLC9A3. Among, STIM1 is a most promising one associated with high reproductive capacity. When compared to Chinese domestic goats with low fecundity (17 individuals), the genes including MLLT10, SPIRE2, TCF25, ZNF276 and FANCA were screened, and the SPIRE2 gene was thought to be associated with fecundity traits. Meanwhile, the functional enrichment of these candidate genes revealed that they were involved in biological processes of mammary gland morphogenesis, uterus development, gastrulation, mesoderm morphogenesis and formation, and blood vessel development, which might undergo natural or artificial selection during reproductive trait formation in goats. Thus, our findings could enrich the genetic basis of reproductive trait selection during goat domestication, which may serve to improve goat breeding practices.

Expression Profile Analysis to Identify Circular RNA Expression Signatures in the Prolificacy Trait of Yunshang Black Goat Pituitary in the Estrus Cycle

The pituitary gland is an important organ. It is a complex area of the brain involved in endocrine function and reproductive regulation. However, the function of the pituitary in goat reproduction is still unclear. Herein, RNA sequencing was used to explore the expression patterns of circle RNAs (circRNAs) in the pituitary of Yunshang black goats during the various estrus phases. Then the host genes of the circRNAs were predicted, and a competing endogenous RNA (ceRNA) network was constructed. The results showed a total of 6,705 circRNAs in the pituitary of Yunshang black goats, among which 388 differentially expressed (DE) circRNAs (214 were upregulated, while 174 were downregulated) were identified between high- and low-yield Yunshang black goats in the follicular phase (HF vs. LF); moreover, 361 DE circRNAs (136 were upregulated, while 225 were downregulated) were identified between high- and low-yield Yunshang black goats in the luteal phase (HL vs. LL). There were 65 DE circRNAs targeting 40 miRNAs in the HF vs. LF comparison and 46 DE circRNAs targeting 31 miRNAs in the HL vs. LL comparison. We identified chi_circ_0030920, chi_circ_0043017, chi_circ_0008353, chi_circ_0041580, and chi_circ_0016478 as the key circRNAs through functional enrichment analysis. The ceRNA network analysis showed that chi_circ_0031209 and chi_circ_0019448 might play an important role in reproduction by influencing the expression of prolactin receptor (PRLR) in high- and low-yielding goats during the luteal phase, whereas chi_circ_0014542 regulates the expression of WNT5A during the follicular phase. Our study provided the overall expression profiles of circRNAs in the goat pituitary during the estrus phase, which provides new insight into the mechanism of high-yield goats, which can be helpful to guide goat breeding.

Transcriptome Analysis Reveals Key miRNA-mRNA Pathways in Ovarian Tissues of Yunshang Black Goats With Different Kidding Numbers

The granulosa cell (GC) is the basic functional unit of follicles, and it is important for promoting follicle growth and sex hormones, as well as growth factor secretion in the process of reproduction. A variety of factors influence granulocyte proliferation, yet there are still many gaps to be filled in target and non-coding RNA regulation. In our study, the differentially expressed (DE) mRNAs and miRNAs were detected by using RNA-seq, and we constructed a mRNA-miRNA network related to goat prolificacy. Then, the goat primary GCs were isolated from the follicle for the function validation of candidate genes and their regulator miRNAs. A total of 2,968 DE mRNAs and 99 DE miRNAs were identified in the high- and low-prolificacy goat by RNA-seq, of which there were 1,553 upregulated and 1,415 downregulated mRNAs, and 80 upregulated and 19 downregulated miRNAs, respectively. JAK3 was identified as highly expressed in the low-prolificacy goats (3 times higher than high-prolificacy goats), and the integrated analysis showed that chi-miR-493-3p was a potential regulator of JAK3. The analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) showed that JAK3 was involved in the PI3K-Akt signaling pathway, the Jak-STAT signaling pathway, and signaling pathways regulating pluripotency of stem cells. In particular, the PI3K-Akt signaling pathway was a typical pathway for cell proliferation, differentiation, apoptosis, and migration. We found that the chi-miR-493-3p targets JAK3 directly via RT-qPCR, dual fluorescence assays, and Western blot. Furthermore, the expression of JAK3 was significantly decreased by the chi-miR-493-3p mimic and increased by the chi-miR-493-3p inhibitor. The CCK-8 assay showed that overexpression of JAK3 promoted cell proliferation, while inhibiting JAK3 had the opposite effect. The expression of cell proliferation markers CDK4 and cyclin D2 also showed the same results. Moreover, the enzyme-linked immunosorbent assay showed that steroid hormones E₂ and PROG were increased by overexpressing JAK3 and decreased by inhibiting JAK3. Therefore, our results identified a chi-miR-439-3p-JAK3 regulatory pathway, which provided a new insight into the GC proliferation and prolificacy of goat.

Combined approaches identify known and novel genes associated with sheep litter size and non-seasonal breeding

Improvement of ewe reproduction is considerable by appropriately increasing litter size and sustaining non-seasonal breeding. However, their genetic makeups have not been entirely elucidated. Genome-wide analyses of 821 individuals were performed by combining three genomic approaches (genome-wide association study, XP-nSL, and runs of homozygosity). Consequently, 35 candidate genes including three domestication genes (TSHR, GTF2A1, and KITLG) were identified. Other than the FecB mutation at BMPR1B, we described a significant association of a missense mutation rs406686139 at seasonal lambing-associated TSHR gene with litter size. Some promising novel genes may be relevant for sheep reproduction by multitude biological processes, such as FETUB functioning in fertilization, HNRNPA1 in oogenesis, DCUN1D1 in spermatogenesis, and HRG in fertility outcome. The present study suggests that improvement of ewe reproduction is attributed to selective breeding, and casts light on the genetic basis and improvement of sheep reproduction.