Genetic effects of DSCAML1 identified in genome-wide association study revealing strong associations with litter size and semen quality in goat (Capra hircus)

Detection of small sequence variations within the goat GHR gene and its effects on growth traits

Growth hormone receptor (GHR) gene is considered to be an important candidate gene in growth traits. Therefore, the purpose of this study was to detect whether there were potential indel variations in the GHR gene that were related to the growth traits of the Shaanbei white cashmere goats (SBWC). In this study, genomic DNA from 931 healthy SBWC individuals were used to verify the relationship between the indel of the GHR gene and growth traits. Two indel variants, P49-bp indel in intron 1 and P1410-bp indel in 3'-UTR, were confirmed. Association analyses demonstrated that these two indel polymorphism loci were associated with the chest circumference and chest width of SBWC. Additionally, for the P49-bp and P1410-bp indel loci, the ID and II genotypes were dominant genotypes, respectively. Moreover, the genotypic distributions of these two indel loci in SBWC were significantly different from those in three other Chinese indigenous goat breeds (HNBG, GZDG and IMWC) (p < 0.05). Taken together, two indel loci (P49-bp indel and P1410-bp indel) both significantly affected the growth traits of goats. This illustrated that these two indel loci might be the potential DNA marker for use in improving the selection and breeding of goats.

Insertion/deletion (Indel) variant of the goat RORA gene is associated with growth traits

RAR related orphan receptor A (RORA), which encodes the retinoid-acid-related orphan receptor alpha (RORα), is a clock gene found in skeletal muscle. Several studies have shown that RORα plays an important role in bone formation, suggesting that RORA gene may take part in the regulation of growth and development. The purpose of this research is to study the insertion/deletion (indel) variations of the RORA gene and investigate the relationship with the growth traits of Shaanbei white cashmere (SBWC) goats. Herein, the current study identified that the P4-11-bp and P11-28-bp deletion sites are polymorphic among 12 pairs of primers within the RORA gene in the SBWC goats (n = 641). Moreover, the P11-28-bp deletion locus was significantly related to the body height (p = 0.046), height at hip cross (p = 0.012), and body length (p = 0.003). Both of P4-11-bp and P11-28-bp indels showed the moderate genetic diversity (0.25 Collapse

Key Words

• RAR related orphan receptor A (RORA)
• clock gene
• goat
• growth traits
• insertion/deletion

Collapse

MESH Headings

• Pregnancy
• Female
• Animals
• Litter Size/genetics
• Goats/physiology
• INDEL Mutation/genetics
• Phenotype

Collapse

Grants Collapse

Affiliation(s)

Qian Zhou College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
Huina Hu College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
Yuta Yang College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
Yuxin Kang College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
Xianyong Lan College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
Xianfeng Wu Institute of Animal Husbandry and Veterinary, Fujian Academy of Agricultural Sciences, Fuzhou, Fujian, China
Zhengang Guo College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China Animal Husbandry and Veterinary Science Institute of Bijie city, Bijie, Guizhou, China
Chuanying Pan College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China

Collapse

Investigation on mRNA expression and genetic variation within goat SMAD2 gene and its association with litter size

The SMAD family member 2 (SMAD2), a member of the TGF-beta superfamily, executes a significant part in the oogenesis and ovulation process. A genome-wide selective sweep analysis also found SMAD2 was different in the fertility groups of Laoshan dairy goats; whether this gene was linked to litter size was unknown. Therefore, SMAD2 was chosen to study its effects on Shaanbei white cashmere goat reproduction and mRNA expression profile. Herein, the mRNA expression level of SMAD2 was firstly determined in female goat tissues, revealing significant differences in mRNA levels of different tissues (p < 0.05), including ovary tissue, indicating a potential role for SMAD2 in goat prolificacy. Then, using six pairs of primers, only one indel locus (P3-Del-12-bp) was found to be polymorphic in goat SMAD2 (n = 501). ANOVA also revealed that a P3-Del-12-bp deletion was significantly related to first-born litter size (p = 0.037). The Chi-square (χ2) test revealed that the ID genotype was significantly more prevalent in mothers with multiple lambs (p = 0.01), indicating that heterozygous individuals (ID) are more likely to produce multiple lambs. Our findings suggest that the SMAD2 gene's P3-Del-12-bp deletion could be used to improve goat breeds by assisting with litter size selection.

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.

A genome-wide perspective on the diversity and selection signatures in indigenous goats using 53 K single nucleotide polymorphism array

Tibetan goats, Taihang goats, Jining grey goats, and Meigu goats are the representative indigenous goats in China, found in Qinghai-Tibet Plateau, Western pastoral area, Northern and Southern agricultural regions. Very few studies have conducted a comprehensive analysis of the genomic diversity and selection of these breeds. We genotyped 96 unrelated individuals, using goat 53 K Illumina BeadChip array, of the following goat breeds: Tibetan (TG), Taihang (THG), Jining grey (JGG), and Meigu (MGG). A total of 45 951 single nucleotide polymorphisms were filtered to estimate the genetic diversity and selection signatures. All breeds had a high proportion (over 95%) of polymorphic loci. The observed and excepted heterozygosity ranged from 0.338 (MGG) to 0.402 (JGG) and 0.339 (MGG) to 0.395 (JGG), respectively. Clustering analysis displayed a genetically distinct lineage for each breed, and their Fst were greater than 0.25, indicating that they had a higher genetic differentiation between groups. Furthermore, effective population size reduced in all four populations, indicating a loss of genetic diversity. In addition, runs of homozygosity were mainly distributed in 5-10 Mb. Lastly, we identified signature genes, which were closely related to high-altitude adaptation (ADIRF) and prolificity (CNTROB, SMC3, and PTEN). This study provides a valuable resource for future studies on genome-wide perspectives on the diversity and selection signatures of Chinese indigenous goats.

Whole-genome selective scans detect genes associated with important phenotypic traits in goat (Capra hircus)

Goats with diverse economic phenotypic traits play an important role in animal husbandry. However, the genetic mechanisms underlying complex phenotypic traits are unclear in goats. Genomic studies of variations provided a lens to identify functional genes. In this study, we focused on the worldwide goat breeds with outstanding traits and used whole-genome resequencing data in 361 samples from 68 breeds to detect genomic selection sweep regions. We identified 210-531 genomic regions with six phenotypic traits, respectively. Further gene annotation analysis revealed 332, 203, 164, 300, 205, and 145 candidate genes corresponding with dairy, wool, high prolificacy, poll, big ear, and white coat color traits. Some of these genes have been reported previously (e.g., KIT, KITLG, NBEA, RELL1, AHCY, and EDNRA), while we also discovered novel genes, such as STIM1, NRXN1, LEP, that may be associated with agronomic traits like poll and big ear morphology. Our study found a set of new genetic markers for genetic improvement in goats and provided novel insights into the genetic mechanisms of complex traits.

A 24-bp indel within the sheep AHR gene is associated with litter size

Aryl Hydrocarbon Receptor (AHR) is a member of the PER-ARNT-SIM (PAS) family, which could mediate various biological processes, for instance, the balance of the immune system, cell proliferation, differentiation, vascular tissue remodeling and reproduction ability regulation. A previous research showed that the AHR gene exerted important functions on the pig reproduction, implying that it could serve as a candidate gene related to animal reproductive traits. Here, the aim of this work was to identify potential insertion/deletion (indel) mutations of the AHR gene in three sheep breeds and analyze the associations between these mutations and reproductive traits. Results showed that a 24-bp indel was uncovered three genotypes (II, ID and DD) in the Australian White sheep (AuW) and Lanzhou fat-tail sheep (LZFT) population, while there were only two genotypes (ID and DD) in Luxi black-headed sheep (LXBH). Moreover, the Fisher's exact test showed that the 24-bp indel mutation was significantly associated with litter size and live litter size in AuW sheep (Fisher's p < 0.05). Therefore, the 24-bp indel of sheep AHR gene can contribute to sheep marker-assisted selection breeding and further improve the sheep reproductive performance.

A 7-nt nucleotide sequence variant within the sheep KDM3B gene affects female reproduction traits

Lysine demethylase 3B (KDM3B) gene is a histone demethylase, demonstrating specific demethylation of the histone H3 lysine 9. It was detected as a sheep reproductive candidate gene by genome-wide scans, and related studies also showed its significance in female reproductive process. However, rare study researched its polymorphism. Herein, we hypothesized that the polymorphisms of KDM3B gene were associated with sheep reproduction traits. A 7-nt nucleotide sequence variant (rs1088697156) within KDM3B gene was identified in a total of 888 individuals, including the Australian White (AUW) sheep and Lanzhou Fat-tailed (LFT) sheep. II (insertion/insertion) and ID (insertion/deletion) genotypes of 7-nt variant were detected, which were at Hardy-Weinberg equilibrium (HWE) in detected breeds. Association analysis illustrated the 7-nt variant was significantly associated with the litter size, duration of pregnancy, live lamb number, live lamb rate, stillbirth number, stillbirth rate of average and different parity (P < 0.05) in AUW sheep. Moreover, 'ID' was the dominant genotype with excellent consistency in reproductive traits. It is instrumental to select individuals with ID genotype for improving the sheep reproduction traits. These findings suggest that the 7-nt variant within KDM3B gene can be used as a candidate marker of reproduction traits for sheep breeding improvement by marker-assisted selection.

Goat SMAD family member 1 (SMAD1): mRNA expression, genetic variants, and their associations with litter size

SMAD family member 1 (SMAD1) is phosphorylated and activated by the BMP receptors, which help regulate ovulation rate, cell growth, apoptosis, and development. Previously, the genome-wide association study revealed that it has been associated with fecundity in sheep. However, its effect on litter size has not been investigated in goats. Therefore, this study aimed to determine the level of SMAD1 mRNA expression in various tissues and to identify its polymorphisms and their association with litter size in Shaanbei white cashmere goat (SBWC). As a result, RT-qPCR analysis showed that SMAD1 was expressed in various tissues in female SBWC goats, including the ovary (P < 0.05). Importantly, the mRNA expression level in the ovaries of mothers of multi-lambs had a higher level than the mothers of single lambs (P < 0.05). Moreover, two InDels (18-bp and 7-bp) in intron 1 of SMAD1 were polymorphic among ten potential loci. Both 18-bp and 7-bp InDels were significantly correlated with litter size (P = 0.014) and (P = 0.0001), respectively. As shown by the chi-squared test, genotypic distributions of 18-bp and 7-bp were significantly distinct between single-lamb (P = 0.02) and multi-lamb mothers (P = 0.002). Our findings confirm that two InDels in SMAD1 were significantly associated with litter size and suggest that they could be used to improve fertility traits in goat breeding strategies.

Short variation of the sheep PDGFD gene is correlated with litter size

Platelet-derived growth factor (PDGF) family, exert plays a key role in embryonic development, cell proliferation, cell migration, angiogenesis and reproduction. Related studies about GWAS analyses have found that PDGFD significantly affected deposition of tail fat in sheep, but there are no studies on reproduction in animals. In this study, three breed of sheep were used to find insertion/deletion (indel) fragment polymorphism of PDGFD which including Australian white (AUW) sheep (Meat type, n = 932), Guiqian semi-fine wool (GSFW) sheep (wool type, n = 60) and East Friensian milk (EFM) sheep (dairy type, n = 60). Only a 18-bp variation was polymorphic in the study AUW sheep population and the genotypes of different sheep breed are also specific. Moreover, the association analysis indicated that this variant was associated with litter size of AUW sheep in the first parity (p < 0.05). The litter size of II genotype was significantly lower than other genotypes in the first parity (p < 0.05). We also revealed that the PDGFD gene was relatively conservative in eight species, PDGFD mRNA expression in 832 sheep samples implying this gene was related to reproduction traits. Hence, these finding demonstrated the one-cause multipotency of PDGFD gene. Collectively, these results suggest that this indel can be used as an effective marker for sheep breeding.

Goat SNX29: mRNA expression, InDel and CNV detection, and their associations with litter size

The sorting nexin 29 (SNX29) gene, a member of the SNX family, is associated with material transport and lipid metabolism. Previous studies have shown that lipid metabolism affects reproductive function in animals. Thus, we hypothesized there is a correlation between the SNX29 gene and reproductive trait. To date, studies on the relationship between the SNX29 gene and reproductive traits are limited. Therefore, the purpose of this study was to examine the polymorphism in the SNX29 gene and its correlation with litter size. Herein, the mRNA expression levels of SNX29 were assayed in various goat tissue. Surprisingly, we found that SNX29 was highly expressed in the corpus luteum, large and small follicles. This result led us to suggest that the SNX29 gene has a critical role in reproduction. We further detected potential polymorphisms in Shaanbei white cashmere (SBWC) goats, including insertion/deletion (InDel, n = 2,057) and copy number variation (CNV, n = 1,402), which were related to fertility. The 17 bp deletion (n = 1004) and the 20 bp deletion (n = 1,053) within the SNX29 gene were discovered to be significantly associated with litter size (P < 0.05), and individuals the ID genotype of P1-Del-17 bp and the DD genotype of P2-Del-20bp had larger litter size. Additionally, the four CNV loci had significant correlations with litter size (P < 0.01) in our detected population. In CNV5, individuals with the median genotype were superior compared to those with loss or gain genotype in term of litter size, and in other three CNVs showed better reproductive trait in the gain genotype. Briefly, these findings suggest that SNX29 could be used as a candidate gene for litter size in goat breeding through marker-assisted selection (MAS).

Genetic Variations and mRNA Expression of Goat DNAH1 and Their Associations with Litter Size

Dynein Axonemal Heavy Chain 1 (DNAH1) encodes proteins which provide structural support for the physiological function and motor structure of spermatozoa (hereafter referred to as sperm) and ova. This study found that three single nucleotide polymorphisms (SNPs), the 27-bp insertion/deletion (InDel) mutations and three exonic copy number variations (CNVs) within DNAH1 were significantly associated with litter size of Shaanbei white cashmere goats (n = 1101). Goats with the wildtypes of these three SNPs had higher litter sizes than other carriers (p < 0.05). II genotype of the 27-bp InDel had the highest litter size compared with ID carriers (p = 0.000022). The gain genotype had the largest litter sizes compared with the loss or medium carriers for the three CNV mutations (p < 0.01). Individuals with the AA-TT-CC-II-M1-M2-M3 and AA-TT-CC-II-G1-G2-M3 combination genotypes had larger litter sizes compared with the other genotypes. This study also showed the DNAH1 expression in mothers of multiple kids was higher than mothers of single kids. These three SNPs, the 27-bp InDel and three CNVs in DNAH1 could be used as molecular markers for the selection of goat reproductive traits.

Candidate Genes in Bull Semen Production Traits: An Information Approach Review

Semen quality plays a crucial role in the successful implementation of breeding programs, especially where artificial insemination (AI) is practiced. Bulls with good semen traits have good fertility and can produce a volume of high semen per ejaculation. The aim of this review is to use an information approach to highlight candidate genes and their relation to bull semen production traits. The use of genome-wide association studies (GWAS) has been demonstrated to be successful in identifying genomic regions and individual variations associated with production traits. Studies have reported over 40 genes associated with semen traits using Illumina BeadChip single-nucleotide polymorphism (SNPs).

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.

Screen of small fragment mutations within the sheep thyroid stimulating hormone receptor gene associated with litter size

The thyroid stimulating hormone receptor (TSHR), a glycoprotein hormone receptor, plays an important role in metabolic regulation and photoperiod control in the time of reproduction in birds and mammals. Previous genome-wide association studies revealed that the TSHR gene was related to reproduction and its function was identified in female reproduction, but rare studies reported the polymorphism of TSHR gene. However, the molecular mutations of the TSHR gene in sheep have not been reported so far. Herein, we explored potential polymorphisms of the sheep TSHR gene, and a 29 bp nucleotide sequence variant (rs1089565492) was identified in the AUW sheep. There were three genotypes of the 29 bp variant locus detected which named 'II' 'DD' and 'ID' been identified. Association analysis results showed the 29 bp variant was significantly associated with the litter size of the AUW sheep (p < 0.05). This finding suggests that the 29 bp nucleotide sequence variant within TSHR gene could be a candidate marker of reproduction traits for sheep breeding improving through the marker-assisted selection (MAS).

Novel InDel variations of the Cry2 gene are associated with litter size in Australian White sheep

Cryptochrome 2 (Cry2) gene regulates circadian rhythm and affects reproduction and pregnancy. Therefore, this study aimed to explore polymorphisms of the Cry2 gene and their associations with litter size at different parity in Australian White (AuW) ewes. Five putative insertion or deletion mutations within the Cry2 gene were selected to study their association with litter size. Two novel deletion mutations were identified in intronic region of Cry2 gene and were genotyped by agarose gel electrophoresis and DNA sequencing. The polymorphism information content (PIC) indicated that both mutations were low polymorphism in tested groups. Statistical analysis revealed that the P1-Del-6-bp was significantly correlated with litter size at third parity (P = 0.010), in which individuals with insertion/deletion (ID) genotype had larger litter size than insertion/insertion (II) genotype (P < 0.05). Whereas, the P2-Del-6-bp was significantly correlated with litter size at first parity (P = 0.036), in which individuals with insertion/insertion (II) genotype had larger litter size than insertion/deletion (ID) genotype (P < 0.05). Collectively, these findings may provide new insights to expedite molecular breeding in sheep through marker-assisted selection strategies (MAS).

Palliative effects of metformin on testicular damage induced by triptolide in male rats

As a widely existing traditional Chinese medicine component, TP (triptolide) has serious reproductive toxicity which causes severe damage to the reproductive system and limits its application prospect. TP and MET (metformin) have shown great potential in combined with each other in anticancer and anti-inflammatory. Whether metformin can resist the reproductive toxicity caused by triptolide, the effects of MET on TP-induced reproductive capacity has not been reported. In this study, metformin was used to investigate the therapeutic effect on reproductive toxicity induced by TP in rat. The results showed that metformin had significant therapeutic effects on oxidative stress damage, destruction of the blood-testosterone barrier and apoptosis. And it proved that its therapeutic effect is mainly to restore the structural and functional stability of testis through antioxidant stress. It will provide guidance for the treatment of reproductive toxicity caused by TP and the adjuvant detoxification of TP application.

Genetic variations of bovine PCOS-related DENND1A gene identified in GWAS significantly affect female reproductive traits

Genome-wide association studies (GWAS) have identified DENND1A as a potential candidate gene linked to the fertility-related phenotypes in dairy cows. However, to date, no studies have examined the association of the DENND1A insertion/deletions (indels) to bovine fertility on a large scale. Herein, two indel sites, including P4-del-26-bp and P8-ins-15-bp were identified in 1064 Holstein cows. The values of the minor allelic frequency (MAF) ranged between 0.471 (deletion) and 0.230 (deletion), respectively, and combined four different haplotypes by analyzing the haplotype combination. It is noteworthy that P4-del-26-bp is associated with the ovarian width (P = 0.0004) and corpus luteum diameter (P = 0.004). Meanwhile, P8-ins-15-bp was found to have a significant association with the ovarian width (P = 0.020), ovarian weight (P = 0.004), the number of mature follicles (P = 0.020), and diameter of the mature follicles (P = 0.016). Furthermore, the combinatorial analysis showed that the two indel combined-genotypes were significantly related to several reproductive traits (ovarian width, ovarian weight, etc.). Collectively, our findings indicated that these two novel indels and their combinations are correlated with the reproductive traits, and hence, they can serve in the marker-assisted selection (MAS) in cattle breeding. Nevertheless, further functional experiments are needed for understanding the mechanisms of these indels in cattle reproduction in a better way.

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

A-kinase anchoring protein 12 (AKAP12) plays key roles in male germ cells and female ovarian granulosa cells, whereas its influence on livestock litter size remains unclear. Herein we detected the genetic variants of AKAP12 gene and their effects on litter size as well as alternative splicing variants expression in Shaanbei white cashmere (SBWC) goats, aiming at exploring theoretical basis for goat molecular breeding. We identified two Insertion/deletions (Indels) (7- and 13-bp) within the AKAP12 gene. Statistical analyses demonstrated that the 13-bp indel mutation in the 3′ UTR was significantly associated with litter size (n = 1,019), and the carriers with DD genotypes presented lower litter sizes compared with other carriers (P < 0.01). Bioinformatics analysis predicted that this 13-bp deletion sequence could bind to the seed region of miR-181, which has been documented to suppress porcine reproductive and respiratory syndrome virus (PRRSV) infection by targeting PRRSV receptor CD163 and affect the pig litter size. Therefore, luciferase assay for this 13-bp indel binding with miRNA-181 was performed, and the luciferase activity of pcDNA-miR-181-13bp-Deletion-allele vector was significantly lower than that of the pcDNA-miR-181-13bp-Insertion-allele vector (P < 0.05), suggesting the reduced binding capability with miR-181 in DD genotype. Given that alternative spliced variants and their expression considerably account for the Indel genetic effects on phenotypic traits, we therefore detected the expression of the alternative spliced variants in different tissues and identified that AKAP12-AS2 exhibited the highest expression levels in testis tissues. Interestingly, the AKAP12-AS2 expression levels of homozygote DD carriers were significantly lower than that of individuals with heterozygote ID, in both testis and ovarian tissues (P < 0.05), which is consistent with the effect of the 13-bp deletion on the reduced litter size. Taken together, our results here suggest that this 13-bp indel mutation within goat AKAP12 might be utilized as a novel molecular marker for improving litter size in goat breeding.

A novel (ATC)n microsatellite locus is associated with litter size in an indigenous Chinese pig

Simple sequence repeats (SSRs) are an important part of the genome and have become powerful auxiliary DNA markers in animal breeding using marker-assisted selection (MAS). Based on previous sequencing data of Qinghai Bamei pigs, a total of three novel candidate SSR loci were analysed in this study. Time-of-flight mass spectrometry (TOF-MS) was used for SSR genotyping, and association analyses between SSRs and the litter size of Qinghai Bamei sows was also performed. The results of genotyping showed that the (ATC)_n -P1, (AC)_n -P2 and (AC)_n -P3 loci had 2, 3 and 18 genotypes, respectively; 2, 3 and 8 alleles were also identified at these loci. Except for the (AC)_n -P2 locus, the polymorphism information content (PIC) values of other loci were greater than 0.25. Association analyses indicated that only the (ATC)_n -P1 locus was significantly associated with the litter size of Qinghai Bamei sows (p = .047). Compared to 189-/189- genotype, individuals with the 189-/195- genotype had the senior litter size, which was 9.04 ± 0.21. Our results enrich the data on SSRs in Qinghai Bamei pigs and indicate that (ATC)_n -P1 is a candidate locus for MAS in the pig industry.

Indel mutations within the bovine HSD17B3 gene are significantly associated with ovary morphological traits and mature follicle number

Given the intensive selection for increased milk production, it is imperative that the problem of declining fertility in dairy cows be reversed. In female mammals their reproductive traits and functioning is controlled by a finely tuned process balancing estrogens and androgens, in which androgens (e.g., testosterone) as a precursor substance can participate in estrogen synthesis by activating 17β-hydroxysteroid dehydrogenase (17betaHSD). Being a key catalyst for testosterone synthesis, we hypothesized HSD17B3 gene is involved in the ovary's development and thereby capable of influencing cows' fecundity. Herein, to investigate the relationship between polymorphisms of the HSD17B3 gene and cow fertility, we characterized three insertion/deletion (indels) polymorphisms of this gene in 1110 healthy bovine ovaries. Their respective minimum allelic frequency (MAF) ranged from 0.180 to 0.482. For the ovary morphological traits, correlations revealed that both P1-D_15-bp and P4-D_19-bp demonstrated significant associations with ovarian height (P = 0.007 and 0.004, respectively), while P5-I_5-bp was found to be significantly associated with the ovarian weight (P = 0.024). For ovarian volume, a significant correlation was uncovered between it and both polymorphisms of P4-D_19-bp (P = 0.036) and P5-I_5-bp (P = 0.045). Cows with either the DD genotype of P4-D_19-bp or P5-I_5-bp tended to have greater ovarian volume, a result consistent with their relationship to ovarian weight (P5-I_5-bp) or height (P4-D_19-bp). For the mature follicle traits, polymorphisms of P4-D_19-bp were found significantly associated with the number of mature follicles (P = 0.045). Furthermore, expression levels of HSD17B3 differed significantly between the maximal and minimum groups of ovarian weight or volume, and the transcription factors GATA-1 and USF were predicted to bind P1-D_15-bp and P4-D_19-bp, respectively. This suggested the detected intron mutations could affect HSD17B3's transcription by regulating the binding of transcription factors, thereby affecting ovarian weight and other reproductive traits. As a potential effective molecular marker loci significantly related to traits of ovary and follicle, these three indels could be used in practical molecular marker-assisted selection (MAS) breeding programs, to optimize female fertility and enhance economic efficiency in the dairy cow industry.

The mRNA expression profile of the goat prion protein testis-specific (PRNT) gene and its associations with litter size

The goat PRNT gene was initially identified as a testis-specific gene with a role in spermatogenesis. In this study, we used quantitative real-time PCR (qPCR) to first determine the mRNA expression profile of this gene in different goat tissues. Surprisingly, we found that PRNT was expressed not only in the testis but also in nine other tissues in goats. Moreover, PRNT was weakly expressed in the testis, while its expression was strongest in the ovary. These results, combined with those of other studies, led us to hypothesize that the goat PRNT gene has a role in both male and female reproduction. We further used direct DNA sequencing to detect potential SNPs within this gene in Shaanbei whit cashmere (SBWC) rams and ewes, and identified three SNPs within the PRNT gene, namely, c.-58C > T, c.71A > G (p.Alanine24Valine), and c.102C > T (synonymous). In rams, c.-58C > T and c.102C > T were strongly linked with each other (D' = 1.000, r² = 0.504), whereas no significant association (P > 0.05) was found between the three SNPs and semen quality, which was consistent with the low expression of the PRNT gene in the testis. Interestingly, in ewes (n = 502), c.-58C > T and c.71A > G were also strongly linked with each other (D' = 0.973, r² = 0.537). Additionally, the c.71A > G locus, especially the AA genotype, had a significant influence on litter size (P = 0.006), consistent with the high PRNT expression in the ovary. Combined, the results of the expression profiling and analysis of the association between the SNPs and reproductive traits showed that two strongly linked nucleotide sequence variants within PRNT were significantly associated with goat litter size. These findings provide potential DNA markers for use in the marker-assisted selection (MAS) of goats with high-fertility traits.

Identification of three new microsatellites and their effects on body measurement traits in pigs using time of flight-mass spectrometry (TOF-MS)

The body status of livestock affects their physiological function and productive performances. Microsatellites, one of the most used DNA markers, have been found to be associated with pig productive traits. However, their identifications and effects on body measurement traits of the Chinese Qinghai Bamei pig still uncovered. According to our previous sequencing data, in this study, three novel microsatellites were found in this breed. Using time of flight-mass spectrometry (TOF-MS) method, these microsatellites were further identified in a large Bamei pig population. TOF-MS spectra showed that there are three microsatellites loci, named P1, P2 and P3. These microsatellites were linkage equilibrium based on the values of D' and r² tests. Association results demonstrated that P1 locus was associated with the body length, body height and chest width and the beneficial genotype was 150-/150-bp (p < 0.05); and P2 locus was associated with the body height (p < 0.05), and the 145-/145-bp, 145-/147-bp and 145-/149-bp were claimed as favorable genotypes and 145-bp allele was considered as the favorable allele. These findings suggested that P1 and P2 microsatellites might be considered as the candidate genetic markers to select pigs with superior body sizes, especially in local breed.

The KMT2A gene: mRNA differential expression in the ovary and a novel 13-nt nucleotide sequence variant associated with litter size in cashmere goats

A genome-wide association study had shown that lysine methyltransferase 2A (KMT2A), which encodes the histone 3 lysine 4 methyltransferase and reportedly can regulate gametogenesis, steroidogenesis, and development as well as other biological processes, is a potential candidate gene influencing litter size in the dairy goat, suggesting its key function in animal reproduction. Here, we aimed to explore the genetic effects of the KMT2A gene on litter size in females of the Chinese indigenous cashmere goat, using a large sample size (n > 1,000), based on their levels of RNA transcription and DNA variation. First, mRNA expression levels of this gene in ovarian tissues between the low-prolific group (first-born litter size = 1) and high-prolific group (first-born litter size ≥2) were significantly different, revealing the potential functioning of KMT2A in goat prolific. Moreover, a novel 13-nt nucleotide sequence variant was identified in Shaanbei white cashmere goats (n = 1,616). In accordance with the independent chi-square (χ²) analysis, the distribution of genotypes (P = 2.57 × 10^-9) and allelotypes (P = 3.00 × 10^-7) between the low- and high-prolific groups differed significantly, indicating the 13-nt mutation was associated with litter size. Further analysis showed that the insertion/insertion (II) genotype was significantly different with insertion/deletion (ID) (P = 1.76 × 10^-9) and deletion/deletion (DD) (P = 7.00 × 10^-6), with goats having the DD genotype producing an average litter size larger than the other genotypes. Taken together, these findings suggest KMT2A can serve as a candidate gene for breeding goats, which may have implications for improving the future development of the goat industry.

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.

Genome-wide association studies for sperm traits in Assaf sheep breed

Sperm quality traits routinely collected by artificial insemination (AI) center for rams progeny test are related with the capacity to produce sperm doses for AI and, in more or less grade, with males' fertility. Low-quality ejaculates are unuseful to perform AI sperm doses, which suppose high economic loses for the AI center. Moreover, sperm quality traits have low heritability values which make traditional genetic selection little efficient to its improvement. In this work, a genome-wide association study (GWAS) was conducted by using sperm quality traits data and 50 K Affymetrix custom chip genotypes of 429 rams of Assaf breed from OVIGEN AI centre. Furthermore, 47 of these rams were also genotyped with the Illumina HD Ovine BeadChip, and therefore HD genotypes were imputed for all rams with phenotype data. Previous to the GWAS, a linear regression model was fitted including sperm traits as dependent variables; the flock of origin, date of sperm collection, and jump number as fixed effects; rams age at collection in months as covariate; and ram permanent effect as random. Pseudo-phenotypes obtained from this model were used as input for GWAS. Associations at the chromosome-wise level (FDR 10%) of 76 single-nucleotide polymorphisms (SNPs) in 4 chromosomes for ejaculate concentration (CON), 20 SNPs in 3 chromosomes for ejaculate volume (VOL), 32 SNPs in 1 chromosome for ejaculate number of spermatozoa (SPZ), and 23 SNPs for spermatozoa mass motility (MOT) in 17 chromosomes were found. Only SNPs associated with MOT overcame the genome-wide significance level. Some candidate genes for sperm traits variability were SLC9C1 (OAR1), TSN (OAR2), and FUT10 (OAR26) for MOT;. DOCK2, CPLANE1, SPEF2, and RAI14 (OAR16) for CON; SCAPER and PSMA4 (OAR18) for VOL; and PARM1 and LOC101110593 (OAR6) for SPZ. SNPs associated with sperm traits were not found to be correlated with milk production genetic variation; however, the high frequencies of some SNPs with negative effect over sperm traits found in animals at the top milk yield estimated breeding values (EBVs) ranking would allow to exert some selective presure to improve rams sperm performances. Effects and frequencies of some of the SNPs detected over sperm quality traits make these variants good candidates to be used in marker-assisted selection to improve sperm characteristics of Assaf rams and AI center efficiency to produce sperm doses.

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.

Identification and Characterization of Hypothalamic Alternative Splicing Events and Variants in Ovine Fecundity-Related Genes

Simple Summary

Previous studies revealed that alternative splicing (AS) events and gene variants played key roles in reproduction. However, their location and distribution in hypothalamic fecundity-related genes in sheep without the FecB mutation remain largely unknown. In this study, we performed a correlation analysis of transcriptomics and proteomics, and the results suggested several differentially expressed genes (DEGs)/differentially expressed proteins (DEPs), including galectin 3 (LGALS3), aspartoacylase (ASPA) and transthyretin (TTR), could be candidate genes influencing ovine litter size. Further analysis suggested that AS events, single nucleotide polymorphisms (SNPs) and microRNA (miRNA)-binding sites existed in key DEGs/DEPs, such as ASPA and TTR. This study provides a new insight into ovine and even other mammalian reproduction.

Abstract

Previous studies revealed that alternative splicing (AS) events and gene variants played key roles in reproduction; however, their location and distribution in hypothalamic fecundity-related genes in sheep without the FecB mutation remain largely unknown. Therefore, in this study, we described the hypothalamic AS events and variants in differentially expressed genes (DEGs) in Small Tail Han sheep without the FecB mutation at polytocous sheep in the follicular phase vs. monotocous sheep in the follicular phase (PF vs. MF) and polytocous sheep in the luteal phase vs. monotocous sheep in the luteal phase (PL vs. ML) via an RNA-seq study for the first time. We found 39 DEGs with AS events (AS DEGs) in PF vs. MF, while 42 AS DEGs were identified in PL vs. ML. No DEGs with single nucleotide polymorphisms (SNPs) were observed in PF vs. MF, but five were identified in PL vs. ML. We also performed a correlation analysis of transcriptomics and proteomics, and the results suggested several key DEGs/differentially expressed proteins (DEPs), such as galectin 3 (LGALS3) in PF vs. MF and aspartoacylase (ASPA) and transthyretin (TTR) in PL vs. ML, could be candidate genes influencing ovine litter size. In addition, further analyses suggested that AS events, SNPs and miRNA-binding sites existed in key DEGs/DEPs, such as ASPA and TTR. All in all, this study provides a new insight into ovine and even other mammalian reproduction.

Relationships between novel nucleotide variants within the colony-stimulating factor 1 receptor (CSF1R) gene and mastitis indicators in sheep

Colony-stimulating factor 1 receptor (CSF1R) plays an important role in the process of innate immunity and inflammation, thus it was hypothesized that the CSF1R gene might affect the occurrence of mammalian mastitis. The purpose of this study was to investigate the association between nucleotide variations of CSF1R gene and mastitis in Australian white sheep (AUWs). Two indel variants (Intron5-27 bp and Intron5-22 bp) within the CSF1R gene have been found in AUWs. The Chi-square test for different mastitis symptoms demonstrated that individuals without symptoms of mastitis had higher 'I' allele frequencies and 'II' genotype frequencies (p < 0.01). We found strong correlation between mastitis and lactation score through Pearson correlation analysis. Therefore, we also analyzed the relationship between the two indel loci and lactation, we found that the lactation ability of individuals with type II was stronger than that of DD genotype at the Intron5-22 bp (p < 0.05). Additionally, we found that the combined genotype of the two loci was significantly associated with mastitis (p < 0.01). These findings indicated that CSF1R mutations were significantly associated with mastitis, and could affect lactation performance, suggesting that two deletion sites could be used as the effective molecular markers against mastitis in sheep breeding.

A 17-bp InDel (rs668420586) within goat CHCHD7 gene located in growth-related QTL affecting body measurement traits

The Coiled-Coil-Helix-Coiled-Coil-Helix Domain Containing 7 (CHCHD7) gene was located in a growth-related major QTL that participated in the process of bone cells metabolism in animals by regulating mitochondrial copper homeostasis and cytochrome C oxidase assembly. Therefore, we speculated that CHCHD7 gene might be involved in animal growth and body size. Herein, we discovered a 17-bp insertion/deletion (indel) within the goat CHCHD7 gene. Then, we detected this variation in Shaanbei White Cashmere (SBWC, n = 1055) goats and Inner Mongolia White Cashmere (IMWC, n = 743) goats (Alathai type) using the mathematical expectation (ME) method. We then analyzed the correlation between these genotypes and goat body measurement traits. The results showed that the minor allelic frequency (MAF) was 0.011 in SBWC, and 0.048 in IMWC. In SBWC and IMWC, the reaction time by ME method was reduced by 36.78% and 27.59%, respectively, compared to the traditional method of screening samples one by one. Moreover, in SBWC goats, the 17-bp indel was significantly associated with body measurement traits (e.g. body height, and body length) in adults. In IMWC goats, the 17-bp indel was correlated with body measurement traits (e.g. body height) in weaners. In SBWC and IMWC goat populations, the body measurement traits of the individuals homozygous for 17-bp indel were higher than those in heterozygous individuals, except for the case of cannon circumference in IMWC weaners. These findings showed that the 17-bp insertion mutation within the goat CHCHD7 gene significantly affected body morphometric traits, and could provide a basis for marker-assisted selection (MAS) breeding of cashmere goats.

Novel indel variations of the sheep FecB gene and their effects on litter size

Various studies had shown that the FecB gene (also known as BMPRIB, BMPR1B or ALK-6) was the major gene influencing sheep litter size, for which its SNPs' variations were reportedly linked. Yet, surprisingly, there was no published information on the insertion/deletion (indel) variation of this gene. Herein, using a population of Chinese Australian White sheep (n = 932), we identified five novel indels in the different introns of the FecB gene. Among them, the 12-bp indel was distinguished as a splicing region variation that was completely linked to the 17-bp indel. Analysis of variance revealed that only the 10-bp indel was significantly associated with sheep litter size (P = 0.010), for which the deletion/deletion (DD) genotype was the harmful one for fecundity. Nevertheless, the combined genotypes of these five indels were significantly (P = 0.033) correlated with the litter size. These findings would provide fresh insight into developing a sounder basis to accelerate molecular breeding in sheep via DNA markers in a marker-assisted selection strategy.

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.

Polymorphic variants of bovine ADCY5 gene identified in GWAS analysis were significantly associated with ovarian morphological related traits

The reproductive performance (e.g. fertility) of dairy cows, which declined over past few decades due to the intense and intensive selection, needs to be improved. Previous genome-wide association study (GWAS) of female Holstein screened the Adenylate cyclase 5 (ADCY5) as the candidate gene for cow fertility. As a member of the adenylyl cyclases family, adenylate cyclase 5 (ADCY5) is famous for regulating extrapyramidal motor system related various neuropsychiatric diseases, and its genetic variant is reported to associate with lower birth and placenta weight which leads to asymmetric fetal growth restriction. It was hypothesized that ADCY5 may affect the fertility of cows by regulating the processes of ovarian development. Herein, genomic DNA from 768 ovaries samples of healthy unrelated Holstein cow were used to screen potential insertion/deletion (indel) mutations using eight pairs of primers, and we found three novel polymorphic indel variants, namely, rs385624978 (P3-D_11-bp), rs433028962 (P5-I_19-bp) and rs382393457 (P8-D_19-bp). The minor allelic frequencies (MAF) of P3-D_11-bp, P5-I_19-bp and P8-D_19-bp loci were 0.188, 0.365 and 0.06, respectively, and there were 7 different haplotypes. Additionally, linkage disequilibrium analysis demonstrated no linkage among them. Importantly, P3-D_11-bp locus was significantly related to both ovarian width (P = 1.0E-6) and corpus luteum diameter (P = 0.015); P5-I_19-bp locus had a significant relation with corpus albicans diameter (P = 0.030) and ovaries with mutational homozygous genotype produced a superior corpus albicans diameter than those with other genotypes. Briefly, three novel indel mutations of bovine ADCY5 gene were identified and two of them were uncovered to be significantly correlated with ovarian phenotypic traits or corpus luteum or albicans traits. These findings contributed to the application of molecular marker-assisted selection (MAS) in improving female fertility in cattle, which could accelerate the development of the cattle industry.

Goat sperm associated antigen 17 protein gene (SPAG17): Small and large fragment genetic variation detection, association analysis, and mRNA expression in gonads

Transcriptome sequencing analyses have suggested that sperm associated antigen 17 protein gene (SPAG17) may play important regulating roles in litter size. In this study, the expression profiles and genetic variations of the SPAG17 were studied in Shaanbei White Cashmere (SBWC) goats (n=1567). SPAG17 was highly expressed in testis and ovary of SBWC goats. At different developmental stages, it also continued to be highly expressed in testis. In addition, two variations of SPAG17, one indel locus and one copy number variation locus, were significantly associated with first-born litter size. Joint analysis results suggested that two polymorphic loci of the SPAG17 gene may regulate host gene expression in goat ovary and testis. Overall, the results indicated the important role of SPAG17 in the reproductive process of goats.

Determination of genetic variation within the DYRK2 gene and its associations with milk traits in cattle

To speed up the progress of marker-assisted selection (MAS) in cattle breeding, the dual-specificity tyrosine phosphorylation-regulated kinase 2 (DYRK2), cadherin 2 (CDH2), and kinesin family member 1A (KIF1A) genes were chosen based on our pervious genome-wide association study (GWAS) analysis results. DYRK2 is a kinase that may participate in cell growth and/or development; it shows phosphorylation activity toward serine, threonine, and tyrosine fragments of proteins, and it is different from other protein kinases. The CDH2 gene encodes a classic cadherin, which is a member of the cadherin superfamily. The protein encoded by KIF1A is a member of the kinesin family and plays a role in the transportation of membrane organelles along axon microtubules. We detected insertion/deletion (InDel) variation in these three candidate genes in 438 individual cattle (Xinjiang Brown cattle and Wagyu × Luxi crossbreed cattle). Only DYRK2-P3-11 bp was polymorphic and genotyped. The polymorphism information content of DYRK2-P3-11 bp was 0.336. Correlation analyses showed that InDel polymorphism was significantly associated with six different milk traits. These findings may aid future analyses of InDel genotypes in cattle breeds, and speed up the progress of MAS in cattle breeding.

Indel variants within the PRL and GHR genes associated with sheep litter size

Growth hormone and prolactin belong to the class of peptide hormones that have a wide range of regulatory functions. In this study, polymorphisms of growth hormone receptor (GHR) and prolactin (PRL) genes were analysed as candidate genes, which are responsible for the litter size in Australian White (AUW) sheep. According to the statistical analyses results, the polymorphism information content (PIC) values of the PRL-P1-ins-23 bp, GHR-P2-del-23 bp and GHR-P8-del-23 bp were 0.371, 0.366 and 0.375, respectively, which indicates the high genetic polymorphism in AUW sheep. Moreover, all indel loci are not conformed to the HWE (p < .05). Further, our findings revealed that the PRL-P1-ins-23 bp polymorphism in the ovine PRL gene was significantly related to the first parity litter size (p = .001) and the DD genotype displaying the highest genotypic mean. Meanwhile, the GHR-P2-del-23 bp and GHR-P8-23 bp indels in the ovine GHR gene were significantly correlated with first parity litter size (p < .05), and the individuals with the genotype II showed significantly higher litter size than others. Collectively, these results demonstrated that our findings could be useful for future sheep breeding strategies based on the molecular-assisted selection (MAS).

A 31-bp indel in the 5' UTR region of GNB1L is significantly associated with chicken body weight and carcass traits

BACKGROUND

G-protein subunit beta 1 like (GNB1L) encodes a G-protein beta-subunit-like polypeptide. Chicken GNB1L is upregulated in the breast muscle of high feed efficiency chickens, and its expression is 1.52-fold that in low feed efficiency chickens. However, no report has described the effects of GNB1L indels on the chicken carcass and growth traits.

RESULTS

This study identified a 31-bp indel in the 5' untranslated region (UTR) of GNB1L and elucidated the effect of this gene mutation on the carcass and growth traits in chickens. The 31-bp indel showed a highly significant association with the body weight at 8 different stages and was significantly correlated with daily gains at 0 to 4 weeks and 4 to 8 weeks. Similarly, the mutation was significantly associated with small intestine length, breast width, breast depth and breast muscle weight. Moreover, DD and ID were superior genotypes for chicken growth and carcass traits.

CONCLUSIONS

These results show that the 31-bp indel of GNB1L significantly affects chicken body weight and carcass traits and can serve as a candidate molecular marker for chicken genetics and breeding programs.

Molecular characterization and a duplicated 31-bp indel within the LDB2 gene and its associations with production performance in chickens

Many studies have shown that the LDB2 gene plays a regulatory role in retinal development and the cell cycle, but its biological role remains unclear. In this study, a 31-bp indel in the LDB2 gene was found for the first time on the basis of 2797 individuals from 10 different breeds, which led to different genotypes among individuals (II, ID and DD). Among these genotypes, DD was the most dominant. Association analysis of an F₂ resource population crossed with the Gushi (GS) chicken and Anka chicken showed that the DD genotype conferred a significantly greater semi-evisceration weight (SEW, 1108.665 g ± 6.263), evisceration weight (EW, 927.455 g ± 5.424), carcass weight (CW, 1197.306 g ± 6.443), breast muscle weight (BMW, 71.05 g ± 0.574), and leg muscle weight (LMW, 100.303 g ± 0.677) than the ID genotype (SEW, 1059.079 g ± 16.86; EW, 879.459 g ± 14.446; CW, 1141.821 g ± 17.176; BMW, 67.164 g ± 1.523; and LMW, 96.163 g ± 1.823). In addition, LDB2 gene expression in different breeds was significantly higher in the breast muscles and leg muscles than in other tissues. The expression level in the breast muscle differed significantly among stages of GS chicken development, with the highest expression observed at 6 weeks. The expression levels in the pectoral muscles differed significantly among Ross 308 genotypes. In summary, we studied the relationships between a 31-bp indel in the LDB2 gene and economic traits in chickens. The indel was significantly correlated with multiple growth and carcass traits in the F₂ resource population and affected the expression of the LDB2 gene in muscle tissue. In short, our study revealed that the LDB2 gene 31-bp indel can be used as a potential genetic marker for molecular breeding.

Two Novel Rare Strongly Linked Missense SNPs (P27R and A85G) Within the GDF9 Gene Were Significantly Associated With Litter Size in Shaanbei White Cashmere (SBWC) Goats

Growth differentiation factor 9 (GDF9) is a high-fertility candidate gene that plays a crucial role in early folliculogenesis in female mammals. In this study, direct sequencing was used to screen possible SNP loci in the goat GDF9 gene. Three SNP loci, p.proline27alanine (P27R), p.leucine61leucine (L61L), and p.alanine85glycine (A85G), were identified in Shaanbei white cashmere (SBWC) goats. Among the three SNPs, two rare missense SNP loci (P27R and A85G) were discovered to be strongly linked with each other (D′ value = 0.926, r² value = 0.703). Both P27R and A85G loci had two genotypes: wild type and heterozygous type. A85G exerted a significant effect on litter size (P = 0.029) in SBWC goats, and the heterozygous genotype was superior in comparison with the wild type. The heterozygous genotype was also superior in P27R but no significant association was found. However, the combination genotypes of P27R and A85G were identified to have superior effects on litter size (P = 3.8E−15). This information suggested that these two SNPs influenced litter size in goats synergistically. Combining this information with our previous studies, we propose that the GDF9 gene is the principal high-fertility candidate gene and that the A85G locus is a promising SNP that affects litter size in goats. These results may fill a research gap regarding rare mutations as well as provide crucial molecular markers that could be useful in marker-assisted selection (MAS) goat rearing when selecting superior individuals.

Detection of the 23-bp nucleotide sequence mutation in retinoid acid receptor related orphan receptor alpha (RORA) gene and its effect on sheep litter size

Retinoid acid receptor related orphan receptor alpha (RORA) transcribes steroid-related genes to regulate estrogen synthesis. As an important reproductive trait, litter size relates to estrogen synthesis. Therefore, it is important to investigate the association between RORA gene and sheep litter size. In this study, one 23-bp nucleotide sequence mutation was identified in intron 1 of RORA gene in 532 female Australian White Sheep. Moreover, the polymorphic information content (PIC) values of this locus was 0.219. The litter size of ID genotype was significantly better than II genotype and DD genotype in the second born litter size (p < 0.05). This loci was related to third born litter size and the ID is the dominant genotype (p < 0.05). The association between combined genotypes and average litter size showed that sheep with heterozygote (ID) genotypes had larger lamb than homozygous (DD and II) genotypes. To sum, this study provided theoretical references for the comprehensively research of the function of RORA gene and the breeding of Australian White Sheep. The 23-bp indel variants could be considered as molecular markers for the second and third born litter size of sheep for MAS breeding.

Detecting novel Indel variants within the GHR gene and their associations with growth traits in Luxi Blackhead sheep

The growth hormone is important in the regulation of metabolism and energy homeostasis and acts through a growth hormone receptor (GHR). In this work, genetic variations within the ovine GHR gene were identified and tested for associations with body morphometric traits in Chinese Luxi Blackhead (LXBH) sheep. Novel deletion loci in the LXBH GHR gene included P2-del-23 bp and P8-del-23 bp indel variants. The polymorphic information content (PIC) was 0.329 in P2-del-23 bp and 0.257 in P8-del-23 bp. Moreover, both indel polymorphisms were not at Hardy-Weinberg equilibrium (p < 0.05) in the LXBH population. Statistical analyses revealed that the P2-del-23 bp and P8-del-23 bp indels were significantly associated (p < 0.05) with several growth traits in rams and ewes, including body weight, body height, chest depth, chest width, chest circumference, cannon circumference, paunch girth and hip width. Among the tested sheep, the body traits of those with genotype DD were superior to those with II and ID genotypes, suggesting that the 'D' allele was responsible for the positive effects on growth traits. Thus, these results indicate that the P2-del-23 bp and P8-del-23 bp indel sites and the DD genotype can be useful in marker-assisted selection in sheep.

Identification of novel alternative splicing of bovine lncRNA lncFAM200B and its effects on preadipocyte proliferation

Adipogenesis is closely related to human health, livestock growth, and meat quality. A previous study identified that bovine lncFAM200B promoter has high activity in 3T3-L1 mice preadipocytes. Thus, lncFAM200B was a candidate gene for regulating adipogenesis. This study aimed to uncover the role of lncFAM200B in bovine adipogenesis and identify novel genetic variations within the bovine lncFAM200B gene. An expression analysis found that lncFAM200B was expressed higher in fat than that in muscle, but the difference was not related to the total methylation level of the promoter active region. Moreover, the expression of lncFAM200B exhibited a significant positive correlation with the expression of C/EBPa during bovine adipocyte differentiation. To uncover the function of lncFAM200B, the full-length lncFAM200B was cloned, and four kinds of transcript variants were found. Protein-coding potential prediction and prokaryotic expression system analysis showed that these four transcript variants were noncoding RNAs. The quantitative reverse-transcription polymerase chain reaction and 5-ethynyl-2'-deoxyuridine assay showed that the transcript variants decreased the messenger RNA expression of Cyclin D1 and inhibited the proliferation of bovine preadipocytes. Considering the important role of lncFAM200B in adipogenesis, we identified genetic variations in lncFAM200B. Three single-nucleotide polymorphisms (SNPs) were revealed, and two of them (SNP1 and SNP3) were associated with Nanyang cattle body measurement traits. In conclusion, this study found that bovine lncFAM200B inhibited preadipocyte proliferation, and two genetic variations of lncFAM200B could be used in cattle breeding.

Goat DNMT3B: An indel mutation detection, association analysis with litter size and mRNA expression in gonads

DNA methyltransferase 3β (DNMT3B) is a gene encoding a de novo methylation enzyme that is required for DNA methylation during mammalian embryo development. Previous genome-wide association analysis suggested DNMT3B is a candidate gene for goat fertility, but there is no study on the effect of DNMT3B on litter size in goat. The aim of this study was to identify possible insertion/deletion (indel) mutations associated with litter size. Seven putative indels were designed to study their association with litter size, but just one 11-bp insertion variant of intron 22 (the last intron) was found in healthy female Shaanbei white cashmere goats (SBWC goats) (n = 1534). Statistical analysis showed that the 11-bp insertion was related to the first-born litter size (P < 0.01) and the goats with the deletion/deletion genotype had a higher average first-born litter size (P < 0.01). In addition, the expression profile of the DNMT3B mRNA in goat was detected, which revealed significant differences in DNMT3B mRNA expression in the gonads. Additionally, the results of western blotting revealed that the ovaries of mothers of multi-lamb (MML) had a higher level of DNMT3B protein than the ovaries of mothers of single-lamb (MSL). Furthermore, the mRNA of DNMT3B was widely expressed in male goats. Differences in mRNA expression levels were observed in the ovaries of MSL and MML. These findings indicated that the 11-bp indel in DNMT3B was significantly associated with first-born litter size, which can be used for marker-assisted selection (MAS) of goats for breeding.

Detection of a 4 bp Mutation in the 3'UTR Region of Goat Sox9 Gene and Its Effect on the Growth Traits

Simple Summary

The sex determining region Y (SRY)-type high mobility group (HMG) box 9 (Sox9) gene is critically important in the formation and development of cartilage and is considered the “main regulator” of chondrogenesis. Additionally, a large number of studies have shown that mutations in a single allele of human Sox9 can lead to campomelic dysplasia syndrome. Therefore, the mutations of Sox9 have been the subject of increasing interest among researchers. However, no studies to date have examined the association between Sox9 gene variants and growth traits in goats. Here, we detected a 4 bp indel in the 3′Untranslated Regions (3′UTR) region of Sox9 in Shaanbei white cashmere (SBWC) goats (n = 1109) and studied the association between this indel and growth traits. The 4 bp indel of Sox9 was significantly associated with body length, heart girth, hip width, and all body measurement indexes (p < 0.05) in SBWC goats. Thus, this deletion could be used as an effective molecular marker for maximizing the growth traits of goats in breeding programs.

Abstract

The SRY-type HMG box 9 (Sox9) gene plays an important role in chondrocyte development as well as changes in hypertrophic chondrocytes, indicating that Sox9 can regulate growth in animals. However, no studies to date have examined the correlation between variations in Sox9 and growth traits in goats. Here, we found a 4 bp indel in the 3′UTR of Sox9 and verified its association with growth traits in Shaanbei white cashmere goats (n = 1109). The frequencies of two genotypes (ID and II) were 0.397 and 0.603, respectively, and polymorphic information content (PIC) values showed that the indel had a medium PIC (PIC > 0.25). The 4 bp indel was significantly correlated with body length (p = 0.006), heart girth (p = 0.001), and hip width (p = 4.37 × 10 ⁻⁴). Notably, individuals with the ID genotype had significantly superior phenotypic traits compared with individuals bearing the II genotype. Hence, we speculated that the 4 bp indel is an important mutation affecting growth traits in goat, and may serve as an effective DNA molecular marker for marker-assisted selection in goat breeding programs.

Goat CMTM2: mRNA expression profiles of different alternative spliced variants and associations analyses with growth traits

CKLF like MARVEL transmembrane domain containing 2 (CMTM2) plays crucial roles in spermiogenesis, skeletogenous, growth, and development through PI3K/Akt and other pathways. The purpose of this study was to explore the expression profile and variation of different spliced CMTM2 gene in Shaanbei white cashmere goats, as well as to find the relationships between a CMTM2 promoter region 14 bp genetic variant and growth traits in 1366 Shaanbei white cashmere goats. In this study, we identified alternative CMTM2 splicing and detected the effects of the spliced variants on mRNA expression levels in tissues. Meanwhile, an unreported spliced variant of CMTM2 in goat was identified using in CDS cloning and RT-PCR, namely, CMTM2-AS2. Compared with the normal transcript (CMTM2-AS1), the novel variant had the higher expression level in muscle and liver tissues, indicating that it plays an effective role in growth traits. Furthermore, a 14 bp deletion was detected within CMTM2 promoter region, and the different genotypes were significantly associated with growth traits (e.g., body length, circumference of cannon bone) in the large group of 1366 individuals in Shaanbei white cashmere goats. We found that the body length of the individuals with II (n = 571) genotype had better phenotypes than those with DD (n = 118) and ID (n = 650) genotypes. These results have direct guiding significance for goat breeding in the future and provide a new idea for studying the characteristics and functions of CMTM2 gene in goats.

Detection of rs665862918 (15-bp Indel) of the HIAT1 Gene and its Strong Genetic Effects on Growth Traits in Goats

Simple Summary

Growth traits are important in goats and can affect their body size and meat production. In this study, the hippocampus abundant transcript 1 (HIAT1) gene, which has been reported as a meat-associated trait in elite goat breeds and a dairy-associated trait in water buffalo, was chosen to detect its correlation with growth traits in goats. The results show that the rs665862918 polymorphism (a 15 bp insertion) in HIAT1 is associated with body length, chest width, chest depth, height at hip cross and heart girth in Shaanbei white cashmere goats (SBWC, n = 1013). Our results reveal that rs665862918 in HIAT1 is relevant to the growth traits of goats and could be used for marker-assisted selection (MAS) as a molecular marker in goat populations.

Abstract

The hippocampus abundant transcript 1 (HIAT1) gene, which was detected by the genome-wide identification of selective sweeps among elite goat breeds and water buffalo, is proposed to play an important role in meat characteristics. Four indels of the HIAT1 gene selected from the NCBI and Ensembl databases were detected via a pooling and sequencing strategy. A 15 bp insertion (rs665862918) in the first intron of HIAT1 was selected and classified on an electrophoresis platform in the Shaanbei white cashmere goat (SBWC) population. The correlation analysis revealed that rs665862918 is significantly highly associated with chest width (p = 1.57 × 10⁻⁵), chest depth (p = 8.85 × 10⁻⁵), heart girth (p = 1.05 × 10⁻⁷), body length (p = 0.022), and height at hip cross (p = 0.023) in the SBWC population (n = 1013). Further analysis revealed that individuals with a genotype insertion/insertion (II) of the rs665862918 locus exhibited better growth trait performance than individuals with an insertion/deletion (ID) or deletion/deletion (DD). These findings verify that HIAT1 affects the body size of goats and that rs665862918 could be a potential molecular marker for growth traits in goat breeding.