Goat AKAP12: Indel Mutation Detection, Association Analysis With Litter Size and Alternative Splicing Variant Expression

Whole-genome resequencing reveals new mutations in candidate genes for Beichuan-white goat prolificacya

In this study, we evaluated the copy number variation in the genomes of two groups of Beichuan-white goat populations with large differences in litter size by FST method, and identified 1739 genes and 485 missense mutations in the genes subject to positive selection. Through functional enrichment, ITGAV, LRP4, CDH23, TPRN, RYR2 and CELSR1 genes, involved in embryonic morphogenesis, were essential for litter size trait, which received intensive attention. In addition, some mutation sites of these genes have been proposed (ITGAV: c.38C > T; TPRN: c.133A > T, c.1192A > G, c.1250A > C; CELSR1: c.7640T > C), whose allele frequencies were significantly changed in the high fecundity goat group. Besides, we found that new mutations at these sites altered the hydrophilicity and 3D structure of the protein. Candidate genes related to litter size in this study and their missense mutation sites were identified. These candidate genes are helpful to understand the genetic mechanism of fecundity in Beichuan white goat, and have important significance for future goat breeding.

Goat CLSTN2 gene: tissue expression profile, genetic variation, and its associations with litter size

Calsyntenin-2 (CLSTN2) is involved in cell proliferation, differentiation, cell death, tumorigenesis, and follicular expression. Although CLSTN2 has been identified as a potential candidate gene for sheep prolificacy, no studies have been done on its effect on goat prolificacy. The purpose of this study was to identify mRNA expression and genetic variation within goat CLSTN2, and its association with prolificacy. Herein, we uncovered significant differences in mRNA levels of the CLSTN2 gene in different tissues in female goats (p < 0.01), including ovary tissue. Nine putative indels were designed to investigate their correlation to litter size, but only one 16-bp deletion was discovered in female Shaanbei white cashmere goats (n = 902). We discovered that a 16-bp deletion within the CLSTN2 gene was significantly correlated with first-born litter size (p = 0.0001). As shown by the chi-squared test, the genotypic II of single-lambs and multi-lambs was dramatically higher than with genotype ID (p = 0.005). Our findings suggest that indel within the CLSTN2 gene is a candidate gene affecting prolificacy in goats and may be used for Marker Assisted Selection (MAS) in goats.

Polymorphism within the GATA binding protein 4 gene is significantly associated with goat litter size

GATA binding protein 4 (GATA4) is a typical transcription binding factor, and its main functions include regulating the proliferation, differentiation and apoptosis of ovarian granulosa cells, promoting spermatogenesis and sex differentiation, implying that this gene have possibly roles in animal reproduction. This study aims to detect five potential insertion/deletions (indels) of the GATA4 gene in 606 healthy unrelated Shaanbei white cashmere (SBWC) goats and analyze its association with the litter size. The electrophoresis and DNA sequencing identified two polymorphic indels (e.g., P4-Del-8bp and P5-Ins-9bp indel). Then T-test analysis showed that P4-Del-8bp was significantly correlated with litter size (p = 0.022) because of two different genotypes detected, e.g., insertion-deletion (ID) and deletion-deletion (DD), and the average litter size of individuals with DD genotype goats was higher than that of others. However, there was no correlation between P5-Ins-9bp and lambing of goats. Chi-square (X2) test found that the distribution of and P4-Del-8bp genotypes (X2 = 6.475, p = 0.011) was significantly different between single and multiple-lamb groups, while P5-Ins-9bp (X2 = 0.030, p = 0.862) was not. Therefore, these findings revealed that P4-Del-8bp polymorphism of goat GATA4 gene was a potential molecular marker significantly associated with litter size, which can be used for the marker-assisted selection (MAS) breeding to improve goat industry.

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

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

Identification of WNT4 alternative splicing patterns and effects on proliferation of granulosa cells in goat

Increasing ovulation numbers is one of the most important ways to promote reproduction in mammals, and follicular granulosa cells (GCs) provide the necessary nutrients and microenvironment for oocytes to ovulate. WNT4 has been shown to be a key factor in regulating the proliferation of GCs in mammalian ovarian tissues. Our previous transcriptome sequencing (RNA-seq) results have identified two alternatively spliced products of WNT4;however, little is known about the splicing mechanism and its effect on GC proliferation. In this study, two alternatively spliced products of WNT4, designated WNT4-α and WNT4-β, were identified by cloning and analyzed for their function by bioinformatics. The RT-qPCR and Western blot results showed that the expression of WNT4-α was significantly higher than that of WNT4-β in the ovary tissues and GCs of Yunshang black goats. We therefore hypothesized that WNT4-α was the main isoform affecting the proliferation of goat GCs. Subsequently, goat GC proliferation assays showed that overexpression of WNT4-α significantly promoted GC proliferation, and the opposite was true after WNT4-α inhibition. The expression of marker genes of the Wnt signaling pathway was also examined and WNT4-α was found to affect the proliferation and hormone secretion of goat GCs by regulating the Wnt signaling pathway. In addition, a series of splicing factors were involved in in the alternative splicing; in this study, SRSF6 was found to be involved as a splicing factor in the generation of WNT4 alternative splicing. In summary, WNT4 alternative splicing was mediated by the splicing factor SRSF6, and WNT4-α alternative splicing played an important role in follicle development and had a significant effect on the proliferation of goat GCs. The results of this study provide a theoretical foundation for further understanding the molecular regulatory mechanisms of the WNT4 in follicle development in goats.

A unique 15-bp InDel in the first intron of BMPR1B regulates its expression in Taihu pigs

BACKGROUND

BMPR1B (Bone morphogenetic protein receptor type-1B) is a receptor in the bone morphogenetic protein (BMP) family and has been identified as a candidate gene for reproductive traits in pigs. Our previous study in Taihu pigs found a specific estrogen response element (ERE) in the first intron of the BMPR1B gene that is associated with the number born alive trait. However, little is known about the mechanism by which the ERE regulates the expression of BMPR1B in the endometrium.

RESULTS

Here, a 15-bp InDel (insertion/deletion) (AGCCAGAAAGGAGGA) was identified as a unique variation in Taihu pigs, and was shown to be responsible for the binding of the type I receptor of estrogen (ESR1) to the ERE using dual-luciferase assays. Four BMPR1B transcripts (T1, T2, T3, and T4) were identified by 5' RACE in endometrial tissue. Expression of T3 and T4 in the endometrium of Meishan pigs was significantly higher than in Duroc pigs during pregnancy. Luciferase assays showed that three distinct BMPR1B promoters may drive expression of T1, T3, and T4. Interestingly, ERE-mediated enhancement of T4 promoter activity significantly increased expression of Transcript T4 in the endometrium of Taihu pigs (P < 0.05). In contrast, the ERE inhibited activity of the T3 promoter and decreased expression of the T3 transcript in the Duroc background (P < 0.05). In summary, we identified a 15-bp InDel in the Taihu ERE that can be used as a molecular marker for the number born alive trait, characterized the 5' untranslated regions (UTRs) of BMPR1B transcripts in the endometrium, and determined how the transcripts are processed by alternative splicing events.

CONCLUSIONS

Our results provide a foundation for understanding the transcriptional regulation of BMPR1B and its contributions to the unique breeding prolificacy characteristics of Taihu pigs.

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.

Exploration of Genetic Variants within the Goat A-Kinase Anchoring Protein 12 (AKAP12) Gene and Their Effects on Growth Traits

Simple Summary

AKAP12, the family of A-kinase anchoring proteins (AKAPs), plays an important role in the regulation of growth and development. There have been no corresponding studies of the effect of the AKAP12 gene on growth traits in goats. In our previous study, 7 bp (intron 3) and 13 bp (3′UTR) indels within the AKAP12 gene significantly influenced AKAP12 gene expression. This study expected to identify the association between these two genetic variations and growth-related traits in 1405 Shaanbei white cashmere (SBWC) goats. The P1–7 bp indel locus was significantly correlated with height at hip cross (HHC; p < 0.05) and the P2–13 bp indel locus was associated with body weight, body length, chest depth, chest width, hip width, chest circumference and cannon (bone) circumference in SBWC goats (p < 0.05). These results prove that the AKAP12 gene plays an important role in the growth and development of goats.

Abstract

The A-kinase anchoring protein 12 gene (AKAP12) is a scaffold protein, which can target multiple signal transduction effectors, can promote mitosis and cytokinesis and plays an important role in the regulation of growth and development. In our previous study, P1–7 bp (intron 3) and P2–13 bp (3′UTR) indels within the AKAP12 gene significantly influenced AKAP12 gene expression. Therefore, this study aimed to identify the association between these two genetic variations and growth-related traits in Shaanbei white cashmere goats (SBWC) (n = 1405). Herein, we identified two non-linkage insertions/deletions (indels). Notably, we found that the P1–7 bp indel mutation was related to the height at hip cross (HHC; p < 0.05) and the P2–13 bp indel was associated with body weight, body length, chest depth, chest width, hip width, chest circumference and cannon (bone) circumference in SBWC goats (p < 0.05). Overall, the two indels’ mutations of AKAP12 affected growth traits in goats. Compared to the P1–7 bp indel, the P2–13 bp indel is more suitable for the breeding of goat growth traits.