Detection of insertions/deletions (InDels) within the goat Runx2 gene and their association with litter size and growth traits

Potential Candidate Genes Associated with Litter Size in Goats: A Review

This review examines genetic markers associated with litter size in goats, a key reproductive trait impacting productivity in small ruminant farming. Goats play a vital socioeconomic role in both low- and high-income regions; however, their productivity remains limited due to low reproductive efficiency. Litter size, influenced by multiple genes and environmental factors, directly affects farm profitability and sustainability by increasing the output per breeding cycle. Recent advancements in genetic research have identified key genes and pathways associated with reproductive traits, including gonadotropin-releasing hormone (GnRH), inhibin (INHAA), Kit ligand (KITLG), protein phosphatase 3 catalytic subunit alpha (PPP3CA), prolactin receptor (PRLR), POU domain class 1 transcription factor 1 (POU1F1), anti-Müllerian hormone (AMH), bone morphogenetic proteins (BMP), growth differentiation factor 9 (GDF9), and KISS1 and suppressor of mothers against decapentaplegic (SMAD) family genes, among others. These genes regulate crucial physiological processes such as folliculogenesis, hormone synthesis, and ovulation. Genome-wide association studies (GWASs) and transcriptomic analyses have pinpointed specific genes linked to increased litter size, highlighting their potential in selective breeding programs. By incorporating genomic data, breeding strategies can achieve higher selection accuracy, accelerate genetic gains, and improve reproductive efficiency. This review emphasizes the importance of genetic markers in optimizing litter size and promoting sustainable productivity in goat farming.

Exploring the distribution of polymorphism across diverse breeds Worldwide in the bovine NR5A2 gene and its correlation with number of mature follicles and corpus albicans

The Nuclear receptor subfamily 5 group A member 2 (NR5A2) gene plays a pivotal role in ovarian development, ovulation, and reproductive traits. There is a lack of studies on its impact on ovarian traits and reproductive traits in cattle. This study aimed to explore NR5A2 gene polymorphisms associations with reproductive traits and investigate the distribution of NR5A2 gene polymorphisms across diverse bovine breeds worldwide. We identified a novel 17-bp deletion within the NR5A2 gene specifically in Chinese Holstein cows (n = 1033) leading to the observation of two genotypes DD and ID. Subsequent association analysis revealed a significant correlation between the 'ID' genotype at this locus and a larger number of corpus albicans (p = 0.042) in diestrus, as well as a higher number of mature follicles (p = 0.038) in estrus. In addition, we also found that the distribution of this deletion exhibits strong regionality across different cattle breeds globally. These findings indicate that the 17-bp deletion mutation within the NR5A2 gene is significantly associated with an increased corpus luteum diameter and a greater number of mature follicles, suggesting its potential utility as a valuable DNA marker for enhancing cow fertility.

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

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

The 12-bp indel in the SMAD family member 2 gene is associated with goat growth traits

SMAD family member 2 (SMAD2) is a member of the TGFβ signaling pathway and functions as an essential regulator in the processes of development, cell proliferation, and bone formation. A previous observation reported that a 12-bp deletion of this gene affected the litter size in goats. However, according to our knowledge, no study has reported an association between this polymorphism and goat body measurement traits. The purpose of this study was to investigate the association of the insertion/deletion (indel) within the SMAD2 gene with the growth traits of goats. The indel polymorphism was found to be significantly associated with chest width and bust (p < 0.05), while cannon circumference was significantly the strongest compared to other traits (p < 0.01) and individuals with the DD genotypes were more dominant genotypes than other genotypes. In summary, we found evidence that the 12-bp indel within the SMAD2 gene could improve goat body measurement traits, paving the way for marker-assisted selection in the field of goat genetics and breeding.

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

InDel and CNV within the AKAP13 Gene Revealing Strong Associations with Growth Traits in Goat

A-kinase-anchoring protein 13 (AKAP13) is a member of the AKAP protein family that has been found to be associated with bone formation. Thus, we investigated the AKAP13 gene as a potential candidate gene for molecular-marker-assisted selection (MAS) in breeding. Our aim was to explore genetic variations (InDel and CNV) within the AKAP13 gene of Shaanbei white cashmere (SBWC) goats and analyze their relationship with growth traits. Ultimately, we identified three InDel loci (16-bp deletion, 15-bp insertion, and 25-bp deletion) and three CNVs, and the 16-bp and 15-bp loci were significantly associated with goat body length (p < 0.05). Both the 16-bp deletion variant and the 15-bp insertion variant facilitated an increase in body length in goats. In addition to this, there was a certain superposition effect between 16-bp and 15-bp loci, although there was no linkage. Additionally, the CNV1 locus was significantly correlated with body height and body length of goats (p < 0.05), and CNV2 was significantly correlated with chest depth, chest circumference, and cannon circumference of goats (p < 0.05). Individuals with gain type showed excellent growth performance. In conclusion, the InDel and CNV loci that we have identified could possibly serve as effective molecular markers in goat breeding, which is very essential for improving efficiency and success of breeding. Moreover, our findings provide a new avenue for further research into the function of the AKAP13 gene.

A Genome-Wide Association Study of the Chest Circumference Trait in Xinjiang Donkeys Based on Whole-Genome Sequencing Technology

Animal genotyping by means of genome-wide association studies is important for connecting phenotypes of interest with their underlying genetics in livestock. However, the use of whole genome sequencing to investigate chest circumference (CC) in donkeys has rarely been reported. We aimed to use the genome-wide association study approach to detect significant single nucleotide polymorphisms (SNPs) and key genes associated with chest circumference traits in Xinjiang donkeys. We assessed 112 Xinjiang donkeys in this study. The chest circumference of each was measured 2 h before milking. We re-sequenced blood samples from the Xinjiang donkeys, and genome-wide association study analyses were performed using a mixed model with the PLINK, GEMMA, and REGENIE programs. We tested 38 donkeys for candidate SNPs for genome-wide association study using three software programs. Additionally, 18 SNP markers reached genome-wide significance (p < 1.61 × 10^-9). On the basis of these, 41 genes were identified. Previously proposed candidate genes for CC traits were supported by this study, including NFATC2 (Nuclear Factor of Activated T Cells 2), PROP1 (PROP Paired-Like Homeobox 1), UBB (Ubiquitin B), and HAND2 (Heart and Neural Crest Derivatives Expressed 2). These promising candidates provide a valuable resource for validating potential meat production genes and will facilitate the development of high-yielding Xinjiang donkey breeds through marker-assisted selection or gene editing.

Screening Indels from the whole genome to identify the candidates and their association with economic traits in several goat breeds

In the present study, the re-sequencing of our previous whole-genome sequencing (WGS) for selected individuals of Dazu-black goat (DBG) and Inner-Mongolia Cashmere goat (IMCG) breeds was used to detect and compare the differentiation in Indels depending on the reference genome of goat. Then, three selected candidate Indels rs668795676, rs657996810, and rs669452874 of the three genes SUFU, SYCP2L and GLIPR1L1, respectively, have been chosen, based on the results of prior GWAS across the genome, and examined for their association with body weight and dimensions (body height, body length, heart girth, chest width, cannon circumference, and chest depth) by kompetitive allele specific PCR assay for 342 goats from the three studied goat breeds (DBG, n = 203; ♂99, ♀104), IMCG (n = 65; 15♂, 50♀), and Hechuan white goat (HWG, n = 74; 34♂, 40♀) breeds. The analysis of 192.747 Gb WGS revealed an average 334,151 Indels in the whole genome of DBG and IMCG breeds. Chromosome 1 had a maximum number of mutations (Indels) of 58,497 and 55,527 for IMCG and DBG, respectively, while chromosome 25 had the least number of mutations of 15,680 and 16,103 for IMCG and DBG, respectively. The majority of Indels were either Ins or Del of short fragments of 1-5 bp, which covered 79.06 and 71.78% of the total number of Indels mutations in IMCG and DBG, respectively. Comparing the differences of Indels between the studied goat breeds revealed 100 and 110 unique Indels for IMCG and DBG, respectively. The Indels loci in the intron region were unique for both studied goat breeds which were related to 30 and 38 candidate genes in IMCG and DBG, respectively, including SUFU, SYCP2L, and GLIPR1L1 genes. Concerning rs669452874 locus, body height and body length of Del/Del genotype in DBG were significantly higher (P < 0.05) than that of Ins/Del genotype, while body height and body length of Del/Del genotype in IMCG were significantly higher (P < 0.01) than those of Ins/Ins genotype, whereas body height and body length and heart girth of Del/Del genotype in HWG were significantly higher (P < 0.01) than those of the Ins/Del and Ins/Ins genotypes. Thus, Del/Del genotype of rs669452874 locus can be used as a candidate molecular marker related to the body dimensions in the studied goat breeds.

Are Copy Number Variations within the FecB Gene Significantly Associated with Morphometric Traits in Goats?

The Booroola fecundity (FecB) gene is a major fertility-related gene first identified in Booroola sheep. Numerous studies have investigated whether the FecB gene is a major fecundity gene in goats or whether there are other genes that play a critical role in goat fertility. Nevertheless, little attention has been paid to the role of the FecB gene in the body morphometric traits of goats, despite the positive relationship discerned between litter size and growth. We identified five copy number variations (CNVs) within the FecB gene in 641 goats, including 318 Shaanbei white cashmere (SBWC) goats, 203 Guizhou Heima (GZHM) goats, and 120 Nubian goats, which exhibited different distributions among these populations. Our results revealed that these five CNVs were significantly associated with goat morphometric traits (p < 0.05). The normal type of CNV3 was the dominant type and displayed superior phenotypes in both litter size and morphometric traits, making it an effective marker for goat breeding. Consequently, LD blocks in the region of 10 Mb upstream and downstream from FecB and potential transcription factors (TFs) that could bind with the CNVs were analyzed via bioinformatics. Although no significant LD block was detected, our results illustrated that these CNVs could bind to growth-related TFs and indirectly affect the growth development of the goats. We identified potential markers to promote litter size and growth, and we offer a theoretical foundation for further breeding work.

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.

Expression of different genotypes of bovine TRDMT1 gene and its polymorphisms association with body measures in Qinchuan cattle (Bos Taurus)

DNA methyltransferase 2 (DNMT2) was renamed as tRNA aspartic acid methyltransferase 1 (TRDMT1) by catalyzing the methylation of tRNA^Asp anti-codon loop C38. The development of sequencing of nucleic acids and protein detection techniques have prompted the demonstration that TRDMT1 mediated tRNA modification affects protein synthesis efficiency. This process affects the growth and development of animals. The DNA of 224 Qinchuan cattles aged 2-4 years old was collected in this experiment. The genetic variations of TRDMT1 exon and some intron regions were detected by mixed pool sequencing technology. qRT-PCR and Western Blot were used to detect the expression levels of mRNA and protein produced with the combination of different genetic variant loci. Three haplotypes were detected and the distribution ratios were different. Muscle tissue mRNA and protein testing showed that there were differences in mRNA expression levels among different genotypes (P < 0.05) and the protein expression levels between different genotypes show the same trend as mRNA. This study provides potential molecular materials for the improvement of Qinchuan cattle reproductivity and provides theoretical support for studying the effects of livestock TRDMT1 on animal growth and development.

Detection of selection signatures in Limousin cattle using whole-genome resequencing

Limousin, a renowned beef breed originating from central France, has been selectively bred over the last 100 years to improve economically important traits. We used whole-genome sequencing data from 10 unrelated Limousin bull calves to detect polymorphisms and identify regions under selection. A total of 13 943 766 variants were identified. Moreover, 311 852 bi-allelic SNPs and 92 229 indels located on autosomes were fixed for the alternative allele in all sequenced animals, including the previously reported missense deleterious F94L mutation in MSTN. We performed a whole-genome screen to discover genomic regions with excess homozygosity, using the pooled heterozygosity score and identified 171 different candidate selective sweeps. In total, 68 candidate genes were found in only 57 of these regions, indicating that a large fraction of the genome under selection might lie in non-coding regions and suggesting that a majority of adaptive mutations might be regulatory in nature. Many QTL were found within candidate selective sweep regions, including QTL associated with shear force or carcass weight. Among the putative selective sweeps, we located genes (MSTN, NCKAP5, RUNX2) that potentially contribute to important phenotypes in Limousin. Several candidate regions and genes under selection were also found in previous genome-wide selection scans performed in Limousin. In addition, we were able to pinpoint candidate causative regulatory polymorphisms in GRIK3 and RUNX2 that might have been under selection. Our results will contribute to improved understanding of the mechanisms and targets of artificial selection and will facilitate the interpretation of GWASs performed in Limousin.

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.

Genetic Polymorphisms and Forensic Efficiencies of a Set of Novel Autosomal InDel Markers in a Chinese Mongolian Group

Insertion/deletion (InDel) markers have been treated as a prospective and helpful aid in the fields of forensic human identifications and biogeography origin researches for the past few years. In this study, we analyzed genetic polymorphisms and forensic efficiencies of 35 InDels in a novel multiplex PCR-InDel panel in a Chinese Mongolian group. All these 35 InDel loci were observed to conform to Hardy–Weinberg equilibrium and linkage equilibrium. The mean values of expected heterozygosity and observed heterozygosity were 0.4788 and 0.4852, respectively. Besides, the interpopulation differentiations and genetic distributions based on 35 InDels found that the Chinese Mongolian group might have closer genetic relationships and similar population genetic structures with East Asian populations.