Genetic Diversity, Selection Signatures, and Genome-Wide Association Study Identify Candidate Genes Related to Litter Size in Hu Sheep.

Hu sheep is a renowned prolific local sheep breed in China, widely distributed across the country due to its excellent reproductive performance. Deciphering the molecular mechanisms underlying the high fecundity of Hu sheep is crucial for improving the litter size of ewes. In this study, we genotyped 830 female Hu sheep using the Illumina OvineSNP50 BeadChip and performed genetic diversity analysis, selection signature detection, and a genome-wide association study (GWAS) for litter size. Our results revealed that the Hu sheep population exhibits relatively high genetic diversity. A total of 4927 runs of homozygosity (ROH) segments were detected, with the majority (74.73%) being short in length. Different genomic inbreeding coefficients (FROH, FHOM, FGRM, and FUNI) ranged from -0.0060 to 0.0126, showing low levels of inbreeding in this population. Additionally, we identified 91 candidate genomic regions through three complementary selection signature methods, including ROH, composite likelihood ratio (CLR), and integrated haplotype score (iHS), and annotated 189 protein-coding genes. Moreover, we observed two significant SNPs related to the litter size of Hu sheep using GWAS analysis based on a repeatability model. Integrating the selection signatures and the GWAS results, we identified 15 candidate genes associated with litter size, among which BMPR1B and UNC5C were particularly noteworthy. These findings provide valuable insights for improving the reproductive performance and breeding of high-fecundity lines of Hu sheep.

Genotyping both live and dead animals to improve post-weaning survival of pigs in breeding programs.

BACKGROUND: In this study, we tested whether genotyping both live and dead animals (GSD) realises more genetic gain for post-weaning survival (PWS) in pigs compared to genotyping only live animals (GOS). METHODS: Stochastic simulation was used to estimate the rate of genetic gain realised by GSD and GOS at a 0.01 rate of pedigree-based inbreeding in three breeding schemes, which differed in PWS (95%, 90% and 50%) and litter size (6 and 10). Pedigree-based selection was conducted as a point of reference. Variance components were estimated and then estimated breeding values (EBV) were obtained in each breeding scheme using a linear or a threshold model. Selection was for a single trait, i.e. PWS with a heritability of 0.02 on the observed scale. The trait was simulated on the underlying scale and was recorded as binary (0/1). Selection candidates were genotyped and phenotyped before selection, with only live candidates eligible for selection. Genotyping strategies differed in the proportion of live and dead animals genotyped, but the phenotypes of all animals were used for predicting EBV of the selection candidates. RESULTS: Based on a 0.01 rate of pedigree-based inbreeding, GSD realised 14 to 33% more genetic gain than GOS for all breeding schemes depending on PWS and litter size. GSD increased the prediction accuracy of EBV for PWS by at least 14% compared to GOS. The use of a linear versus a threshold model did not have an impact on genetic gain for PWS regardless of the genotyping strategy and the bias of the EBV did not differ significantly among genotyping strategies. CONCLUSIONS: Genotyping both dead and live animals was more informative than genotyping only live animals to predict the EBV for PWS of selection candidates, but with marginal increases in genetic gain when the proportion of dead animals genotyped was 60% or greater. Therefore, it would be worthwhile to use genomic information on both live and more than 20% dead animals to compute EBV for the genetic improvement of PWS under the assumption that dead animals reflect increased liability on the underlying scale.

Genetic regulation of ovulation rate and multiple births.

Ovulation rate in many mammalian species is controlled to regulate the numbers of offspring and maximise reproductive success. Pathways that regulate ovulation rate still respond to genetic and environmental factors and show considerable variation within and between species. Genetic segregation, positional cloning, and association studies have discovered numerous mutations and genetic risk factors that contribute to this variation. Notable among the discoveries has been the role of mutations in bone morphogenetic protein 15 (BMP15 ), growth differentiation factor 9 (GDF9 ) and bone morphogenetic protein receptor type 1B (BMPR1B ) from the intra-ovarian signalling pathway contributing to the evidence that signalling from the oocyte is the key driver in follicle regulation rather than circulating gonadotrophin concentrations. Multiple variants in different domains of BMP15 and GDF9 result in partial or complete loss of function of the proteins providing insights into their functional roles and differential regulation contributing to species differences in ovulation rate. Early success encouraged many more studies in prolific strains of sheep, cattle and goats providing a valuable catalogue of genetic variants of large effect increasing ovulation rate and litter size. More recently, genetic association studies are beginning to identify genetic risk factors with smaller effects. Most genes implicated are from pathways with defined roles in regulation of the ovarian function. However, some genomic regions suggest regulation by novel genes. Continuing genetic and related functional studies will add further to our understanding of the detailed regulation of ovulation rate and litter size with implications for health and animal production systems.

Genetic assessment of litter size, body weight, carcass traits and gene expression profiles in exotic and indigenous rabbit breeds: a study on New Zealand White, Californian, and Gabali rabbits in Egypt.

Rabbits are essential for commercial meat production due to their efficient growth and productivity, breeds like New Zealand White (NZW), Californian (CAL), and Gabali (GAB) rabbits offer unique genetic traits in litter, growth, and carcass traits. This study aimed to evaluate heritability (h2), genetic and phenotypic correlations (rg and rp) for litter size, body weight and carcass traits across California (CAL), New Zealand white (NZW) and Gabali (GA) rabbits. Along with exploring gene expression profiles of TBC1D1, NPY, AGRP, POMC, Leptin, GH, GHR, IGF-1, CAA, GPR, ACC, CPT1, FAS, and CART in the brain, liver, and meat tissues of different rabbit breeds. The breed genotype had a significant impact on litter size (LS), litter weight (LW), body weight at 12 weeks (BW12), and daily weight gain (DWG) traits. NZW rabbits displayed superior performance in terms of litter size and litter weight, while CAL rabbits recorded the highest values for BW12 and DWG. Heritability estimates (h2) were generally low for litter size (ranging from 0.05 to 0.12) and medium for body weight (ranging from 0.16 to 0.31). Both genetic (rg) and phenotypic (rp) correlations for litter size were positive and moderate (ranging from 0.08 to 0.48), while correlations for body weight ranged from 0.21 to 0.58. Additionally, CAL rabbits exhibited higher carcass traits compared to NZW and GA rabbits. In terms of breed-specific gene expression patterns, New Zealand White (NZW) rabbits displayed the highest expression levels of key genes related to energy metabolism (TBC1D1), appetite regulation (NPY, AGRP, POMC), nutrient transport (CAA), and G protein-coupled receptors (GPR) in both brain and liver tissues. Californian (CAL) rabbits exhibited superior gene expression of the ACC gene in brain tissue and GH, GHR, and IGF-1 genes in brain and meat tissues. Gabali (GAB) rabbits demonstrated the highest expression levels of TBC1D1, NPY, AGRP, GPR, and ACC genes in meat tissues. These breed-specific gene expression differences, combined with genetic evaluation efforts, have the potential to enhance reproductive and productive performance in rabbits, offering valuable insights for rabbit breeding programs and genetic selection.

Association Analysis of METTL23 Gene Polymorphisms with Reproductive Traits in Kele Pigs.

Methyltransferase-like 23 (METTL23) is a kind of RNA methyltransferase that catalyzes the methylation transfer to the N6-adenosine of RNA, serving as one of the key mediators in this process. However, the METTL23 gene has been poorly researched in pigs. In this study, we investigated the genetic effects of METTL23 single-nucleotide polymorphism(SNPs) on reproductive traits in Kele pigs. The DNA was extracted from 228 healthy multiparous Kele sows, and Sanger sequencing revealed three SNPs, g.4804958 G > T (intron 2), g.4805082 C > T (exon 2), and g.4806821 A > G (exon 3). The polymorphism information content (PIC) for each SNP was 0.264, 0.25, and 0.354, indicating moderate polymorphism (0.25 < PIC < 0.5) and providing genetic information. Linkage disequilibrium analysis showed no strong linkage disequilibrium between the three SNPs. The association analysis revealed that in the SNP g.4804958 G > T individuals with the GG genotype had a significantly higher number of piglets born alive, litter birth weight, number of weaned piglets, and weaning litter weight compared to those with the TT genotype (p < 0.05). Individuals with the GG genotype in the SNP g.4806821 A > G group had significantly higher litter birth weight and average birth weight than those with the AA genotype (p < 0.05). The H4H4 diplotype showed significant effects on the number of piglets born alive, litter birth weight, number of weaned piglets, weaning litter weight, and weaning weight (p < 0.05). Together, the METTL23 gene could be used as a candidate gene for the selection of reproductive traits in Kele pigs.

Sheep litter size heredity basis using genome-wide selective analysis.

Sheep are important herbivorous domestic animal globally, and the Chinese indigenous sheep breed has a multitude of economically significant variations due to the diverse geographical and ecological conditions. In particular, certain native breeds exhibit a visible high litter size phenotype due to the selection pressure of natural and artificial for thousands of years, offering an ideal animal model for investigating sheep's fecundity. In this study, selective signal analysis was performed on public whole-genome sequencing data from 60 sheep across eight breeds to identify candidate genes related to litter size. Results revealed that a total of 34,065,017 single-nucleotide polymorphisms (SNPs) were identified from all sheep, and 65 candidate genes (CDGs) were pinpointed from the top 1% of interacted windows and SNPs between the pairwise fixation index (FST, >0.149543) and cross-population extended haplotype homozygosity (XP-EHH, >0.701551). A total of 41 CDGs (e.g. VRTN, EYA2 and MCPH1) were annotated to 576 GO terms, of which seven terms were directly linked to follicular and embryonic development (e.g. TBXT, BMPR1B, and BMP2). In addition, 73 KEGG pathways were enriched by 21 CDGs (e.g. ENTPD5, ABCD4 and RXFP2), mainly related to Hippo (TCF4, BMPR1B and BMP2), TGF-ß (BMPR1B and BMP2), PI3K-Akt (ITGB4, IL4R and PPP2R5A) and Jak-STAT signalling pathways (IL20RA and IL4R). Notably, a series of CDGs was under strong selection in sheep with high litter size traits. These findings result could improve the comprehension of the genetic underpinnings of sheep litter size. Furthermore, it provides valuable CDGS for future molecular breeding.

Assessing genetic parameters for first parity reproductive traits for early selection in Saanen × Beetal goats of India.

The present study aimed to evaluate the genetic parameters of first parity reproductive traits. Information on 762 reproductive records on Saanen × Beetal (S × B) goats reared for approximately five decades was collected from ICAR-National Dairy Research Institute, Karnal, Haryana (1973-2020). For genetic analysis, single-trait and multiple-trait animal models were used. Gibbs sampler for animal model (GSAM) approach was used for estimating (co)variance components of reproductive traits. Six different single-trait animal models (with or without maternal and environmental effects) were used and the deviance information criterion (DIC) determined the best model. The least squares mean for age at first service (AFS), age at first kidding (AFK), service period (SP), dry period (DP), gestation length (GL), kidding interval (KI), litter weight (LW), number of kids born (NKB) and number of female kids born (NFKB) in first parity were 526.99 ± 4.86, 662.96 ± 5.03, 219.11 ± 6.25, 109.38 ± 6.00, 150.48 ± 0.27, 356.63 ± 4.80 days, 3.87 ± 0.05 kg, 1.27 ± 0.02 and 0.67 ± 0.03, respectively. Lower heritability estimates for these reproductive traits revealed a sparse scope for genetic improvement. Multivariate analysis using Model 1 was carried out to evaluate the genetic and phenotypic correlation of these nine reproductive traits. The genetic correlation of DP and SP was negatively with LW, NKB and NFKB, which is favourable as reduction in SP and DP can improve these economically important traits through indirect selection. Consistent efforts towards genetic improvement of these goat flock poses a promising future for meat industry owing to high prolificacy and good reproductive potential in this flock.

Genetic Links between Reproductive Traits and Amino Acid Pairwise Distances of Swine Leukocyte Antigen Alleles among Mating Partners in Microminipigs.

Previously, we found that a greater dissimilarity in swine leukocyte antigen (SLA) class I and class II alleles between mating partners resulted in increased farrowing rates in a highly inbred population of Microminipigs (MMPs). In this follow-up study, we have analyzed the effects of dissimilarity in SLA alleles between mating partners for seven different reproductive traits, including litter size and the number of stillborn and live or dead weaned piglets. We determined the relationships among reproductive traits within each mating event and the amino acid distances of SLA alleles as markers of diversity between mating partners. Our results indicate that mating partners with greater amino acid pairwise genetic distances in the SLA-1 class I gene or DQB1 class II gene alleles were associated with significantly larger litter sizes and higher numbers of live piglets at birth and weaning. Also, partners with greater pairwise distances in the SLA-2 class I gene alleles exhibited fewer pre-weaning deaths. These findings suggest that the dissimilarity in SLA class I and class II alleles between mating partners may affect not only farrowing rates but also other key reproductive traits such as litter size and improved piglet survival rates. Consequently, SLA alleles could serve as valuable genetic markers for selecting mating partners in breeding programs and for conducting epistatic studies on various reproductive traits in MMPs.

Estimation of genotype by environmental interaction for litter traits by reaction norm model in Taiwan Landrace sows.

The negative effects of heat stress on swine reproduction have been well documented and the recent global warming trend caused by climate change is leading to more days with high temperatures every year. This has caused a reduction in litter trait performance of Landrace sows in Taiwan, a country extending across tropical and subtropical oceanic zones. Therefore, this study developed a modified model to determine which stages of pregnancy, before, early, middle, and late, had the largest impacts of heat stress on litter traits. A reaction norm model (RNM) was used to identify sows with high resilience to heat stress for litter traits followed by analysis of the modified model. Data from Landrace sows were collected from 2 farms in Taiwan between 2008 and 2021. A total of 11,059 records were collected for total number born (TNB), number born alive (NBA), and stillborn rate (STBR). The results showed that the heritabilities of TNB, NBA, and STBR were 0.170, 0.115, and 0.077, respectively. These results were similar between the conventional model and the modified model. In the modified model, the before and early stages of sow pregnancy were the significant periods for TNB and NBA (P < 0.05), while the early and middle stages were significant for STBR (P < 0.05). According to the RNM results, the heritability estimates for TNB, NBA, and STBR were 0.23 to 0.11, 0.18 to 0.08, and 0.10 to 0.04, respectively, showing a decrease from low temperature-humidity index (THI) to high THI. The minimum genetic correlations between the highest and the lowest THI for TNB, NBA, and STBR were 0.85, 0.64, and 0.80, respectively. The results of the RNM for breeding value showed re-ranking across THI values. In conclusion, similar results were obtained for heritability when the model was modified for heat stress estimation. Yet re-ranking of breeding values across THI could help farmers to select not only for improved litter trait performance but also for heat stress resilience of Landrace sows in Taiwan.

Heat stress caused by climate change is a challenge for the pig industry, especially in countries located in tropical and subtropical zones, such as Taiwan. It can adversely affect litter traits, leading to less pork production and higher economic losses to farms. Therefore, identifying sows with the potential to tolerate high heat with high humid condition is an important task for the pig industry. This article proposes a reaction norm model to determine the trend in breeding value across temperature and humidity index values and its implications for litter traits of sows. Our results indicate that litter traits can be used to select sows with high heat tolerance.

Analysis of goat PPP6C mRNA profile, detection of genetic variations, and their associations with litter size.

The Protein Phosphatase 6 Catalytic Subunit (PPP6C) is evolutionarily a conserved gene in eukaryotes known to play a significant role in mammalian reproduction. This study aimed to investigate expression patterns of PPP6C and explore its association with litter size in Shaanbei white cashmere (SBWC) goats. Initially, we determined the mRNA expression levels of PPP6C in both male and female goats across multiple tissues. The results showed that PPP6C mRNA was expressed in multiple tissues, with higher levels in the testis and fallopian tubes, suggesting its involvement in goat reproduction. Additionally, we identified a novel 19 bp InDel within the PPP6C gene in a population of 1030 SBWC goats, which exhibited polymorphism. Statistical analysis revealed a significant association between the19 bp InDel mutation and litter size (P < 0.05). Subsequent, bioinformatics analysis, including linkage disequilibrium (LD) block and selective scanning, highlighted the linkage tendency among most InDel loci did not stand out within B-8 block, there were still some InDel loci linked to the 19 bp within a relatively narrow region. Furthermore, comparative analysis with Bezoars, these selective signals all indicated that this gene was under higher selection pressure, implying that the 19 bp InDel locus within the PPP6C is potentially associated with domesticated traits, particularly in relation to litter size. The results of the present study suggest that the PPP6C is a vital candidate gene affecting prolificacy in goats, with implications for selective breeding programs for goat breeds.

Placental Transcriptome Analysis in Connection with Low Litter Birth Weight Phenotype (LBWP) Sows.

It is possible to identify sub-populations of sows in every pig herd that consistently give birth to low birth weight (BW) piglets, irrespective of the litter size. A previous study from our group demonstrated that placental development is a main factor affecting the litter birth weight phenotype (LBWP) in sows, thereby impacting the BW of entire litters, but the biological and molecular pathways behind this phenomenon are largely unknown. The aim of this study was to investigate the differential gene expression in placental tissues at day 30 of gestation between low LBWP (LLBWP) vs. high LBWP (HLBWP) sows from a purebred Large White maternal line. Using mRNA sequencing, we found 45 differentially expressed genes (DEGs) in placental tissues of LLBWP and HLBWP sows. Furthermore, (GO) enrichment of upregulated DEGs predicted that there were two biological processes significantly related to cornification and regulation of cell population proliferation. To better understand the molecular interaction between cell proliferation and cornification, we conducted transcriptional factor binding site (TFBS) prediction analysis. The results indicated that a highly significant TFBS was located at the 5' upstream of all four upregulated genes (CDSN, DSG3, KLK14, KRT17), recognized by transcription factors EGR4 and FOSL1. Our findings provide novel insight into how transcriptional regulation of two different biological processes interact in placental tissues of LLBWP sows.

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.

Association of polymorphism in the promotor area of the caprine BMPR1B gene with litter size and body measurement traits in Damani goats.

This study aimed to explore polymorphisms in the promoter region of the caprine BMPR1B (Bone morphogenetic protein receptor 1 beta) gene and its association with body measurement and litter size traits in Damani does. A total of 53 blood samples were collected to analyze the association between the BMPR1B gene polymorphism and 11 phenotypic traits in Damani female goats. The results revealed that three novel SNPs were identified in the promoter region of the caprine BMPR1B gene, including g.67 A > C (SNP1), g.170 G > A(SNP2), and g.501A > T (SNP3), among which the SNP1 and SNP2 were significantly (p < 0.05) associated with litter size and body measurement traits in Damani goats. In SNP1 the AC genotype could be used as a marker for litter size, and the CC genotype for body weight in Damani goats. In SNP2, the genotype GG was significantly (p < 0.05) associated with ear and head length. Therefore, we can conclude from the present study, that genetic variants AC and CC of the caprine BMPR1B gene could be used as genetic markers for economic traits through marker-assisted selection for the breed improvement program of the Damani goat.

Genetic variant of the sheep E2F8 gene and its associations with litter size.

The economic efficiency of sheep breeding, aiming to enhance productivity, is a focal point for improvement of sheep breeding. Recent studies highlight the involvement of the Early Region 2 Binding Factor transcription factor 8 (E2F8) gene in female reproduction. Our group's recent genome-wide association study (GWAS) emphasizes the potential impact of the E2F8 gene on prolificacy traits in Australian White sheep (AUW). Herein, the purpose of this study was to assess the correlation of the E2F8 gene with litter size in AUW sheep breed. This work encompassed 659 AUW sheep, subject to genotyping through PCR-based genotyping technology. Furthermore, the results of PCR-based genotyping showed significant associations between the P1-del-32bp bp InDel and the fourth and fifth parities litter size in AUW sheep; the litter size of those with genotype ID were superior compared to those with DD and II genotypes. Thus, these results indicate that the P1-del-32bp InDel within the E2F8 gene can be useful in marker-assisted selection (MAS) in sheep.

Microbiome composition as a potential predictor of longevity in rabbits.

BACKGROUND: Longevity and resilience are two fundamental traits for more sustainable livestock production. These traits are closely related, as resilient animals tend to have longer lifespans. An interesting criterion for increasing longevity in rabbit could be based on the information provided by its gut microbiome. The gut microbiome is essential for regulating health and plays crucial roles in the development of the immune system. The aim of this research was to investigate if animals with different longevities have different microbial profiles. We sequenced the 16S rRNA gene from soft faeces from 95 does. First, we compared two maternal rabbit lines with different longevities; a standard longevity maternal line (A) and a maternal line (LP) that was founded based on longevity criteria: females with a minimum of 25 parities with an average prolificacy per parity of 9 or more. Second, we compared the gut microbiota of two groups of animals from line LP with different longevities: females that died/were culled with two parities or less (LLP) and females with more than 15 parities (HLP). RESULTS: Differences in alpha and beta diversity were observed between lines A and LP, and a partial least square discriminant analysis (PLS-DA) showed a high prediction accuracy (> 91%) of classification of animals to line A versus LP (146 amplicon sequence variants (ASV)). The PLS-DA also showed a high prediction accuracy (> 94%) to classify animals to the LLP and HLP groups (53 ASV). Interestingly, some of the most important taxa identified in the PLS-DA were common to both comparisons (Akkermansia, Christensenellaceae R-7, Uncultured Eubacteriaceae, among others) and have been reported to be related to resilience and longevity. CONCLUSIONS: Our results indicate that the first parity gut microbiome profile differs between the two rabbit maternal lines (A and LP) and, to a lesser extent, between animals of line LP with different longevities (LLP and HLP). Several genera were able to discriminate animals from the two lines and animals with different longevities, which shows that the gut microbiome could be used as a predictive factor for longevity, or as a selection criterion for these traits.

An Identification of Functional Genetic Variants in B4GALNT2 Gene and Their Association with Growth Traits in Goats.

ß-1,4-N-acetylgalactosamine transferase 2 (B4GALNT2) is a vital candidate gene that affects the growth traits in sheep. However, whether it has the same function in goats remains to be investigated further. This study selected 348 Nanjiang Yellow goats, screened all exons, and conserved non-coding regions of the B4GALNT2 gene for single-nucleotide polymorphisms (SNPs). Our results revealed the presence of a synonymous mutation, rs672215506, within the exon of the B4GALNT2 gene in the Nanjiang Yellow goat population. The mutation resulted in a decrease in the mRNA stability of the B4GALNT2 gene. The results of SNP detection of the conserved non-coding region of the B4GALNT2 gene showed five potential regulatory SNPs in the Nanjiang Yellow goat population. Except for rs66095343, the ~500 bp fragments of the other four SNPs (rs649127714, rs649573228, rs652899012, and rs639183528) significantly increased the luciferase activity both in goat skeletal muscle satellite cells (MuSCs) and 293T cells. The genetic diversity indexes indicated low or intermediate levels for all six SNPs analyzed, and the genotype frequencies were in Hardy-Weinberg equilibrium. Association analysis showed that rs660965343, rs649127714, and rs649573228 significantly correlate with growth traits in the later stage of growth and development of Nanjiang Yellow goats. The haplotype combinations of H2H3 and H2H2 had higher body weight and greater body size. Moreover, H2H2 haplotype combinations significantly correlated with the litter size of the Nanjiang Yellow goats. The results of our study demonstrate the potential role of the B4GALNT2 gene as a functional genetic marker in the breeding programs of Nanjiang Yellow goats.

Genome-Wide DNA Methylation Analysis and Functional Validation of Litter Size Traits in Jining Grey Goats.

DNA methylation (DNAm) is associated with the reproductive system. However, the genetic mechanism through which DNAm regulates gene expression and thus affects litter size in goats is unclear. Therefore, in the present work, genome-wide DNAm profiles of HP and LP Jining Grey goat ovary tissues were comprehensively analyzed via WGBS, and RNA-Seq data were combined to identify candidate genes associated with litter size traits in the Jining Grey goat. Finally, BSP and RT-qPCR were used to verify the sequencing results of the key genes. Notably, the DNMT genes were downregulated at the expression level in the HP group. Both groups exhibited comparable levels of methylation. A total of 976 differentially methylated regions (DMRs) (973 DMRs for CG and 3 DMRs for CHG) and 310 differentially methylated genes (DMGs) were identified in this study. Through integration of WGBS and RNA-Seq data, we identified 59 differentially methylated and differentially expressed genes (DEGs) and ultimately screened 8 key DMGs (9 DMRS) associated with litter size traits in Jining Grey goats (SERPINB2: chr24_62258801_62259000, NDRG4: chr18_27599201_27599400, CFAP43: chr26_27046601_27046800, LRP1B. chr2_79720201_79720400, EPHA6: chr1_40088601_40088800, TTC29: chr17_59385801_59386000, PDE11A: chr2_117418601_117418800 and PGF: chr10_ 16913801_16914000 and chr10_16916401_16916600). In summary, our research comprehensively analyzed the genome-wide DNAm profiles of HP and LP Jining Grey goat ovary tissues. The data findings suggest that DNAm in goat ovaries may play an important role in determining litter size.

FecB Was Associated with Litter Size and Follows Mendel's Laws of Inheritance When It Transited to Next Generation in Suhu Meat Sheep Breeding Population.

In order to investigate the effect of FecB on litter size and growth and development traits of Suhu meat sheep and the inheritance patterns of FecB between parents and offspring in the population. In this experiment, 2241 sheep from the Suhu meat sheep population were tested for FecB using capillary electrophoresis. We combined the lambing records of 473 ewes, the growth trait records of 881 sheep at both the birth and weaning (2-month-old) stages, and the complete genealogical records of 643 lambs to analysis the distribution of FecB in the Suhu meat sheep breeding population, its effect on litter size of ewes, growth and development of lambs, and the inheritance patterns of FecB. The results showed that there were three genotypes of FecB in the Suhu meat sheep population, namely the AA genotype, AG genotype, and GG genotype. FecB in this population has a moderate polymorphism (0.25 < PIC < 0.5), and deviates from Hardy-Weinberg disequilibrium (p < 0.05). The litter size of GG genotype ewes was significantly higher than that with the AG and AA genotypes (p < 0.01). A Chi-square test showed that the inheritance patterns of FecB follows Mendel's Laws of Inheritance (p > 0.05). An association analysis of different genotypes of FecB with body weight and body size of Suhu meat sheep at birth and weaning revealed that FecB adversely affects the early growth and development of Suhu meat sheep. In summary, FecB can improve the litter size of ewes but it has negative effects on the early growth and survival rate of lambs in sheep. Therefore, FecB test results and feeding management measures should be comprehensively applied to improve the reproductive performance of ewes, the survival rate and production performance of lambs in sheep production, and thus improve the economic benefits of sheep farms.

4D label-free quantitative proteomic analysis identifies CRABP1 as a novel candidate gene for litter size in rabbits†.

In commercial rabbit breeding, litter size is a crucial reproductive trait. This trait directly determines the reproductive ability of female rabbits and is crucial for evaluating the production efficiency. We here compared differentially expressed proteins of in the ovary tissue from New Zealand female rabbits with high (H) and low (L) litter sizes by using 4D label-free quantitative proteomic technology and identified 92 differential proteins. The biological functions of these proteins were revealed through gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Most distributions of GO and KEGG were related to reproduction, growth development, and metabolism. Furthermore, a novel candidate gene cellular retinoic acid binding protein-1 (CRABP1), which was highly expressed in the L group, was selected for further biological function verification. The Cell Counting Kit-8 (CCK-8) assay and flow cytometry analysis revealed that CRABP1 can promote granulosa cell (GC) apoptosis and inhibit GC proliferation. Furthermore, qRT-PCR and western blotting analysis revealed that CRABP1 regulates the genes (HSD17B1, Wnt-10b, FSHR, TAF4B, BMP15, and BMP6) and protein (Wnt-10b) associated with steroid hormone synthesis and follicle development. The PCR product direct sequencing method revealed single nucleotide polymorphisms in the core promoter region of CRABP1. Luciferase activity assays revealed that the transcriptional activity of the GG genotype was significantly higher than that of the TT or TG genotype. Different genotypes are accompanied by changes in transcription factors, which indicates that T-359G polymorphism can regulate CRABP1 expression. In general, we identified litter size-related genes and revealed the mechanism underlying the effect of CRABP1 on litter size. CRABP1 serves as a key factor in the reproductive capacity of rabbits and can act as a molecular biomarker for the breeding of New Zealand rabbits.

Estimation of genetic parameters for reproductive indices in sheep.

This study aimed to estimate two reproductive efficiency indices in sheep based on the ratio between litter weight (at birth and weaning) and dam weight, as well as their genetic parameters. Phenotypic and pedigree data comprising the period from 1990 to 2018 were obtained from the Santa Inês sheep database of Embrapa Tabuleiros Costeiros. For estimation of the genetic parameters of the indices, a repeatability model was applied in single- and two-trait analyses by a Bayesian approach. The mean reproductive efficiency index was 0.069 ± 0.0163 and 0.43 ± 0.0955 at birth and weaning, respectively. These values indicate that, on average, ewes give birth to 69 g of lamb per kg body weight and wean 430 g of lamb per kg body weight. Described here for the first time, the heritability estimate obtained in single- and two-trait analyses was 0.24 for the index based on birth weights and ranged from 0.13 to 0.15 for the index based on weaning weights. The estimates indicate the possibility of genetic gain by selection and are similar to those reported for reproductive traits in sheep, representing an option for selection criterion. The genetic correlation between indices was positive and moderate (0.26). The repeatability estimates were high (0.49 for the birth weight index and 0.71 for the weaning weight index). These values indicate good prediction of future performance with few observations. The weaning weight index might be a good culling criterion of females.