Selection signatures of litter size in Dazu black goats based on a whole genome sequencing mixed pools strategy

Runs of Homozygosity Detection and Selection Signature Analysis for Local Goat Breeds in Yunnan, China

Runs of Homozygosity (ROH) are continuous homozygous DNA segments in diploid genomes, which have been used to estimate the genetic diversity, inbreeding levels, and genes associated with specific traits in livestock. In this study, we analyzed the resequencing data from 10 local goat breeds in Yunnan province of China and five additional goat populations obtained from a public database. The ROH analysis revealed 21,029 ROH segments across the 15 populations, with an average length of 1.27 Mb, a pattern of ROH, and the assessment of the inbreeding coefficient indicating genetic diversity and varying levels of inbreeding. iHS (integrated haplotype score) was used to analyze high-frequency Single-Nucleotide Polymorphisms (SNPs) in ROH regions, specific genes related to economic traits such as coat color and weight variation. These candidate genes include OCA2 (OCA2 melanosomal transmembrane protein) and MLPH (melanophilin) associated with coat color, EPHA6 (EPH receptor A6) involved in litter size, CDKAL1 (CDK5 regulatory subunit associated protein 1 like 1) and POMC (proopiomelanocortin) linked to weight variation and some putative genes associated with high-altitude adaptability and immune. This study uncovers genetic diversity and inbreeding levels within local goat breeds in Yunnan province, China. The identification of specific genes associated with economic traits and adaptability provides actionable insights for utilization and conservation efforts.

IGF2BP2 regulates the proliferation and migration of endometrial stromal cells through the PI3K/AKT/mTOR signaling pathway in Hu sheep

Insulin-like growth factor 2 mRNA binding protein 2 (IGF2BP2), a significant member of the conserved RNA-binding protein family, plays various roles in numerous physiological and pathological processes. However, the specific function of IGF2BP2 in regulating endometrial function in sheep remains largely unknown. In this study, we observed a significant upregulation in IGF2BP2 mRNA abundance in the endometrium during the luteal phase compared to the follicular phase in Hu sheep. The knockdown of IGF2BP2 resulted in accelerated cell proliferation and migration of Hu sheep endometrial stromal cells (ESCs). Moreover, RNA sequencing analysis revealed that genes with significantly altered expression in IGF2BP2 knockdown cells were predominantly enriched in endometrial receptivity-related signaling pathways, such as cytokine-cytokine receptor interaction, NOD-like receptor, PI3K-AKT, and JAK-STAT signaling pathway. Additionally, the knockdown of IGF2BP2 significantly increased the expression of matrix metalloprotein 9 (MMP9), vascular endothelial growth factor, and prolactin (PRL) in ESCs. The knockdown of IGF2BP2 was also observed to stimulate the PI3K/AKT/mTOR pathway by upregulating integrin β4 (ITGB4) expression. Notably, the downregulation of ITGB4 attenuates IGF2BP2 knockdown-induced facilitation of proliferation and migration of Hu sheep ESCs by inhibiting the PI3K/AKT/mTOR pathway. Collectively, these findings highlight the important role of IGF2BP2 in regulating endometrial function, particularly through the modulation of ESC proliferation and migration via the PI3K/AKT/mTOR pathway.

Candidate genes for litter size in Xinjiang sheep identified by Specific Locus Amplified Fragment (SLAF) sequencing

The aim of this study was to investigate the selection signatures at a genome-wide level in 'Pishan' sheep using Specific Locus Amplified Fragment (SLAF)-seq. Blood samples from 126 ewes were sequenced using SLAF tags, and the ovarian tissues from 8 ewes (Bashbay sheep, a single litter size group (SG group); 'Pishan' sheep, double litter size group (DG group)) were collected to detect expression levels by quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR). Selection signature analysis was performed using global fixation index (Fst) and nucleotide diversity (π) ratio. A total of 1,192,168 high-quality SLAFs were identified. Notably, 2380 candidate regions under selection using two approaches were identified. A total of 2069 genes were identified, which were involved in dopaminergic synapses, thyroid hormone synthesis, ovarian steroidogenesis and thyroid hormone signalling pathways. Furthermore, Growth Differentiation Factor 9 (GDF9), Period Circadian Regulator 2 (PER2), Thyroid Stimulating Hormone Receptor (TSHR), and Nuclear Receptor Coactivator 1 (NCOA1) reside within these regions and pathways. The expression levels of GDF9 and PER2 genes in sheep tissue of the DG group were significantly higher than those in the SG group. These genes are interesting candidates for litter size and provide a starting point for further identification of conservation strategies for 'Pishan' sheep.

Uterus proliferative period ceRNA network of Yunshang black goat reveals candidate genes on different kidding number trait

Pregnancy loss that occurs in the uterus is an important and widespread problem in humans and farm animals and is also a key factor affecting the fecundity of livestock. Understanding the differences in the fecundity of goats may be helpful in guiding the breeding of goats with high fecundity. In this study, we performed RNA sequencing and bioinformatics analysis to study the uterus of Yunshang black goats with high and low fecundity in the proliferative period. We identified mRNAs, long non-coding RNAs (lncRNAs), and microRNAs (miRNAs) by analyzing the uterine transcriptomes. The target genes of the identified miRNAs and lncRNAs were predicted, and miRNA-mRNA interaction and competitive endogenous RNA (ceRNA) networks were constructed. By comparisons between low- and high-fecundity groups, we identified 1,674 differentially expressed mRNAs (914 were upregulated, and 760 were downregulated), 288 differentially expressed lncRNAs (149 were upregulated, and 139 were downregulated), and 17 differentially expressed miRNAs (4 were upregulated, and 13 were downregulated). In addition, 49 miRNA-mRNA pairs and 45 miRNA-lncRNA pairs were predicted in the interaction networks. We successfully constructed a ceRNA interaction network with 108 edges that contained 19 miRNAs, 11 mRNAs, and 73 lncRNAs. Five candidate genes (PLEKHA7, FAT2, FN1, SYK, and ITPR2) that were annotated as cell adhesion or calcium membrane channel protein were identified. Our results provide the overall expression profiles of mRNAs, lncRNAs, and miRNAs in the goat uterus during the proliferative period and are a valuable reference for studies into the mechanisms associated with the high fecundity, which may be helpful to guide goat to reduce pregnancy loss.

Whole genome sequencing of simmental cattle for SNP and CNV discovery

BACKGROUD

The single nucleotide polymorphisms (SNPs) and copy number variations (CNVs) are two major genomic variants, which play crucial roles in evolutionary and phenotypic diversity.

RESULTS

In this study, we performed a comprehensive analysis to explore the genetic variations (SNPs and CNVs) of high sperm motility (HSM) and poor sperm motility (PSM) Simmental bulls using the high-coverage (25×) short-read next generation sequencing and single-molecule long reads sequencing data. A total of ~ 15 million SNPs and 2,944 CNV regions (CNVRs) were detected in Simmental bulls, and a set of positive selected genes (PSGs) and CNVRs were found to be overlapped with quantitative trait loci (QTLs) involving immunity, muscle development, reproduction, etc. In addition, we detected two new variants in LEPR, which may be related to the artificial breeding to improve important economic traits. Moreover, a set of genes and pathways functionally related to male fertility were identified. Remarkably, a CNV on SPAG16 (chr2:101,427,468 - 101,429,883) was completely deleted in all poor sperm motility (PSM) bulls and half of the bulls in high sperm motility (HSM), which may play a crucial role in the bull-fertility.

CONCLUSIONS

In conclusion, this study provides a valuable genetic variation resource for the cattle breeding and selection programs.

Integration of transcriptome sequencing and whole genome resequencing reveal candidate genes in egg production of upright and pendulous-comb chickens

Egg production performance plays an important role in the poultry industry across the world. Previous studies have shown a great difference in egg production performance between pendulous-comb (PC) and upright-comb (UC) chickens. However, there are no reports to identify potential candidate genes for egg production in PC and UC chickens. In the present study, 1,606 laying chickens were raised, and the egg laid by individual chicken was collected for 100 d. Moreover, the expression level of estrogen and progesterone hormones was measured at the start-laying and peak-laying periods of hens. Besides, 4 PC and 4 UC chickens were selected at 217 d of age to perform transcriptome sequencing (RNA-seq) and whole genome resequencing (WGS) to screen the potential candidate genes of egg production. The results showed that PC chicken demonstrated better egg production performance (P < 0.05) and higher estrogen and progesterone hormone expression levels than UC chicken (P < 0.05). RNA-seq analysis showed that 341 upregulated and 1,036 downregulated differentially expressed genes (DEGs) were identified in the ovary tissues of PC and UC chickens. These DEGs were mainly enriched in protein-related, lipid-related, and nucleic acids-related biological processes including ribosome, peptide biosynthetic process, lipid transport terms, and catalytic activity acting on RNA which can significantly affect egg production in chicken. The enrichment results of WGS analysis were consistent with RNA-seq. Further, joint analysis of WGS and RNA-seq data was utilized to screen 30 genes and CAMK1D, CLSTN2, MAST2, PIK3C2G, TBC1D1, STK3, ADGRB3, and PPARGC1A were identified as potential candidate genes for egg production in PC and UC chickens. In summary, our study provides a wealth of information for a better understanding of the genetic and molecular mechanism for the future breeding of PC and UC chickens for egg production.

Investigation of Copy Number Variations (CNVs) of the Goat PPP3CA Gene and Their Effect on Litter Size and Semen Quality

Simple Summary

PPP3CA is one of the candidate genes for goat reproduction, but no studies have been carried out yet. Therefore, the purpose of this study was to determine the associations between copy number variations in the goat PPP3CA gene and litter size and semen quality in goats, including Shaanbei white cashmere goats (SBWC) (n = 353) and Guizhou Heima (GZHM) goats (n = 64). Based on the association analysis, the results showed that only CNV1 (copy number variation 1) and CNV2 (copy number variation 2) were distinctly related to the first-birth litter size in female goats (p = 7.6802 × 10⁻¹¹; p = 5.0895 × 10⁻⁹), and they were also significantly associated with the semen quality of SBWC goats (p < 0.05). These findings prove that the PPP3CA gene plays an important role in reproduction traits in goats.

Abstract

Copy number variations (CNVs) have many forms of variation structure, and they play an important role in the research of variety diversity, biological evolution and disease correlation. Since CNVs have a greater impact on gene regulation and expression, more studies are being finalized on CNVs in important livestock and poultry species. The protein phosphatase 3 catalytic subunit alpha (PPP3CA) is a key candidate gene involved in the goat fecundity trait, and has important effects on precocious puberty, estrogen signal transduction pathways and oocyte meiosis. Additionally, PPP3CA also has a dephosphorylation effect in the process of spermatogonial stem cell meiosis and spermatogenesis. So far, there is no research on the relationship between the copy number variations of the PPP3CA gene and reproduction traits. Therefore, the purpose of this study was to determine the association between copy number variations in the goat PPP3CA gene and litter size and semen quality in Shaanbei white cashmere goats (SBWC) (n = 353) and Guizhou Heima goats (n = 64). Based on the association analysis, the results showed that only CNV1 and CNV2 within the PPP3CA gene were distinctly related to the first-birth litter size in female goats (p = 7.6802 × 10⁻¹¹; p = 5.0895 × 10⁻⁹, respectively) and they were also significantly associated with the semen quality of SBWC goats (p < 0.05). In addition, individuals with Loss genotypes demonstrated better phenotypic performance compared to those with other types. Therefore, CNV1 and CNV2 of the PPP3CA gene are potentially useful for breeding, as they are linked to important goat reproduction traits.

The SNP-Based Profiling of Montecristo Feral Goat Populations Reveals a History of Isolation, Bottlenecks, and the Effects of Management

The Montecristo wild goat is an endangered feral population that has been on the homonymous island in the Tuscan Archipelago since ancient times. The origins of Montecristo goats are still debated, with authors dating their introduction either back to Neolithic times or between the 6th and 13th century of the Common Era. To investigate the evolutionary history and relationships of this population we assembled a 50K SNP dataset including 55 Mediterranean breeds and two nuclei of Montecristo goats sampled on the island and from an ex situ conservation project. Diversity levels, gene flow, population structure, and genetic relationships were assessed through multiple approaches. The insular population scored the lowest values of both observed and expected heterozygosity, highlighting reduced genetic variation, while the ex situ nucleus highlighted a less severe reduction. Multivariate statistics, network, and population structure analyses clearly separated the insular nucleus from all other breeds, including the population of Montecristo goats from the mainland. Moreover, admixture and gene flow analyses pinpointed possible genetic inputs received by the two Montecristo goat nuclei from different sources, while Runs of Homozygosity (ROHs) indicated an ancient bottleneck/founder effect in the insular population and recent extensive inbreeding in the ex situ one. Overall, our results suggest that Montecristo goats experienced several demographic fluctuations combined with admixture events over time and highlighted a noticeable differentiation between the two nuclei.

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.

Insertion/deletion variants within the IGF2BP2 gene identified in reported genome-wide selective sweep analysis reveal a correlation with goat litter size

Insulin-like growth factor 2 mRNA-binding protein 2 (IGF2BP2, also called IMP2) plays an essential role in the development and maturation of germ cells and embryos and is a candidate gene for goat litter size, based on a previous genome-wide selective sweep analysis. In this study, the mRNA expression level of IGF2BP2 was found to be significantly higher in a single-lamb group than in a multi-lamb group. Insertions/deletions (indels) within the goat IGF2BP2 gene, including P4-Ins-13bp and P5-Del-12bp, were verified in 918 Shaanbei White Cashmere (SBWC) female goats. The minor allelic frequencies (MAFs) of P4-Ins-13bp and P5-Del-12bp loci were 0.349 and 0.295, respectively. Analysis using the Chi-square (χ²) test showed that the genotype (χ²=14.479, P=0.006) distribution of P4-Ins-13bp was significantly different between the single-lamb and multi-lamb groups. Correlation analysis demonstrated that P4-Ins-13bp was significantly associated with goat litter size (P=0.022), and individual goats with the homozygous deletion/deletion (DD) genotype produced more litters than other goats. Therefore, considered as a potential molecular marker significantly related to lambing traits, the P4-Ins-13bp mutation of the goat IGF2BP2 gene can be used in goat breeding with practical molecular marker-assisted selection (MAS) to optimize female reproduction and improve economic efficiency in the goat industry.

Detecting the footprint of selection on the genomes of Murciano-Granadina goats

Artificial selection is one of the major forces modifying the genetic composition of livestock populations. Identifying genes under selection could be useful to elucidate their impact on phenotypic variation. We aimed to identify genomic regions targeted by selection for dairy and pigmentation traits in Murciano-Granadina goats. Performance of a selection scan based on the integrated haplotype score test in a population of 1183 Murciano-Granadina goats resulted in the identification of 77 candidate genomic regions/SNPs. The most significant selective sweeps mapped to chromosomes 1 (69.86 Mb), 4 (41.80-49.95 Mb), 11 (65.74 Mb), 12 (31.24 and 52.51 Mb), 17 (34.76-37.67 Mb), 22 (31.75 Mb), and 26 (26.69-31.05 Mb). By using previously generated RNA-Seq data, we built a catalogue of 6414 genes that are differentially expressed across goat lactation (i.e. 78 days post-partum, early lactation; 216 days post-partum, late lactation; 285 days post-partum, dry period). Interestingly, 183 of these genes mapped to selective sweeps and several of them display functions related with lipid, protein, and carbohydrate metabolism, insulin signaling, cell proliferation, as well as mammary development and involution. Of particular interest are the CSN3 and CSN1S2 genes, which encode two major milk proteins. Additionally, we found three pigmentation genes (GLI3, MC1R, and MITF) co-localizing with selective sweeps. Performance of a genome-wide association study and Sanger sequencing and TaqMan genotyping experiments revealed that the c.801C>G (p.Cys267Trp) polymorphism in the melanocortin 1 receptor (MC1R) gene is the main determinant of the black (GG or GC genotypes) and brown (CC genotypes) colorations of Murciano-Granadina goats.

Deletion mutation within the goat PPP3CA gene identified by GWAS significantly affects litter size

The protein phosphatase 3 catalytic subunit α (PPP3CA) gene is a high reproduction traits candidate gene for goats as revealed by a genome-wide association study. The aim of this work was to explore the genetic variations of the goat PPP3CA as well as to evaluate the genetic effects on litter size. Three novel insertions/deletions (indels) within the goat PPP3CA were found and their minor allelic frequencies (MAF) were 0.105, 0.066, and 0.042, respectively. The results showed that only the 20bp indel polymorphism was significantly associated with litter size in Shaanbei white cashmere goats (P<0.05) and individuals with deletion/deletion (DD) genotypes demonstrated the junior phenotypes when compared with those with other genotypes. These findings suggested that the 20bp indel is a potential DNA marker for selecting superior individuals in marker-assisted selection for breeding concerning fecundity in goats.

Whole-genome sequencing to identify candidate genes for litter size and to uncover the variant function in goats (Capra hircus)

To select candidate genes for goat prolificacy, we managed six multi- and six single-kid female goats at the same feeding level and in the same management mode over a 4-year period. These goats showed stable differences in litter size over five continuous parturition records. Whole-genome re-sequencing was used in all 12 to select candidate genes, namely, AURKA, ENDOG, SOX2, RORA, GJA10, RXFP2, CDC25C, and NANOS3, by the strength of their differentiation signals. Most of the selected genes were enriched in the coiled coil process and ovarian development, which suggests that the coiled coil process has a potential regulatory effect on fecundity. Detection of the distribution of variants and association analyses with litter size in 400 goats showed that NANOS3 exon mutations may lead to a transformation of the protein structure. The variation in CDC25C, ENDOG, and NANOS3 showed a significant association with litter size. These results can contribute to the improvement of reproduction traits in the artificial breeding of goats.

Combined approaches to reveal genes associated with litter size in Yunshang black goats

Intensive artificial selection has been imposed in Yunshang black goats, the first black specialist mutton goat breed in China, with a breeding object of improving reproductive performance, which has contributed to reshaping of the genome including the characterization of SNP, ROH and haplotype. However, variation in reproductive ability exists in the present population. A WGS was implemented in two subpopulations (polytocous group, PG, and monotocous group, MG) with evident differences of litter size. Following the mapping to reference genome, and SNP calling and pruning, three approaches - GWAS, ROH analysis and detection of signatures of selection - were employed to unveil candidate genes responsible for litter size. Consequently, 12 candidate genes containing OSBPL8 with the minimum P-value were uncovered by GWAS. Differences were observed in the pattern of ROH between two subpopulations that shared similar low inbreeding coefficients. Two ROH hotspots and 12 corresponding genes emerged from ROH pool association analysis. Based on the nS_L statistic, 15 and 61 promising genes were disclosed under selection for MG and PG respectively. Of them, some promising genes participate in ovarian function (PPP2R5C, CDC25A, ESR1, RPS26 and SERPINBs), seasonal reproduction (DIO3, BTG1 and CRYM) and metabolism (OSBPL8, SLC39A5 and SERPINBs). Our study pinpointed some novel promising genes influencing litter size, provided a comprehensive insight into genetic makeup of litter size and might facilitate selective breeding in goats.

Genome-wide association analysis reveals the genetic locus for high reproduction trait in Chinese Arbas Cashmere goat

BACKGROUND

Litter size is the most important reproductive trait which plays a crucial role in goat production. Therefore, improvement of litter size trait has been of increasing interest in goat industry as small improvement in litter size may lead to large profit. The recent Cashmere goat breeding program produced a high-reproductive genetic line of Arbas Cashmere goat. But the genetic mechanism of high reproduction rate remains largely unknown in this Chinese native goat breed. To address this question, we performed a genome-wide association studies (GWAS) using two groups of goats varying in fecundity.

OBJECTIVES

Our study was aimed to investigate the significant SNPs and genes associated with high reproduction trait in Inner Mongolia Arbas Cashmere Goat.

METHODS

We used logistic model association to perform GWAS using 47 goats from high fecundity group (~ 190%) and 314 goats from low fecundity group (~ 130%) of the Arbas Cashmere goat breed.

RESULTS

We identified 66 genomic regions associated with genome wide significant level wherein six loci were found to be associated with reproduction traits. Further analysis showed that five key candidate genes including KISS1, KHDRBS2, WNT10B, SETDB2 and PPP3CA genes are involved in goat fecundity trait. Gene ontology enrichment analysis revealed that several biological pathways could be involved in the variation of fecundity in female goats.

CONCLUSIONS

The identified significant SNPs or genes provide useful information about the underlying genetic control of fecundity trait which will be helpful to use them in goat breeding programs for improving the reproductive efficiency of goats.