Association of MYF5 and KLF15 gene polymorphisms with carcass traits in domestic pigeons (Columba livia)

1. Single nucleotide polymorphisms (SNPs) in the exons of the myogenic factor 5 (MYF5) and Kruppel-like factor 15 (KLF15) genes were identified and analysed by using DNA sequencing methods in 60 female domestic pigeons (Columba livia).

2. Five SNPs (T5067A, C5084T, C5101T, T5127A and C5154G) were detected in exon 3 of MYF5 and 6 SNPs (C1398T, C1464T, G1542A, C1929T, G1965A and A2355G) were found in exon 2 of KLF15, respectively. The analysis revealed three genotypes, in which the AA genotype was dominant and the A allele showed a dominant advantage.

3. For the MYF5 gene, the C5084T and T5127A SNP genotypes were significantly associated with carcass traits of pigeons. Within those two SNPs, the BB genotype showed relatively higher trait association values than those of AA or AB genotypes. No significant association was observed between the KLF15 SNP genotypes and carcass traits.

4. These results indicated that the MYF5 gene is a potential major gene affecting carcass traits in domestic pigeons. The BB genotype of the C5084T and T5127A SNPs could be a potential candidate genetic marker for marker-assisted selection in pigeon.

     

Association of ADSL gene polymorphisms with meat quality and carcass traits in domestic pigeons (Columba livia)

ABSTRACT

1. Adenylosuccinate lyase (ADSL) plays an important role in the synthesis of inosine monophosphate (IMP). In this trial, a total of 200 pigeons were sampled and slaughtered. Seven meat quality traits and 11 carcass traits were measured. DNA sequencing was used to detect nucleotide mutations, and associations between ADSL gene polymorphisms and meat quality and carcass traits were analysed.

2. Sequencing results showed that 9 nucleotide mutations were found in the exons of the ADSL gene. All the mutations were synonymous except C13065G, which caused a change in amino acids (Ser to Arg). In addition, two of the detected single nucleotide polymorphisms (SNPs) had significant associations with meat quality and carcass traits.

3. For the C13065G SNP located in exon11, the IMP content of breast muscle in AA and AB genotype was higher than in the BB genotype (P < 0.01). The SNP G10696A located in exon10 was significantly associated with carcass rate, in which AA and AB genotype were higher than the BB genotype (P < 0.01).

4. The results indicated that the ADSL gene has a close association with meat quality and carcass traits in domestic pigeons, and G10696A and C13065G SNPs could be applied as genetic markers in molecular marker-assisted breeding of pigeons.     

Impact of MYOD family genes on pork traits in Large White and Landrace pigs

Summary Porcine myogenic differentiation genes ( MYOD ) family play a key role in growth and muscle development and are therefore considered as candidate genes for meat production traits. The objective of the study was to investigate the polymorphisms at four loci belonging to the MYOD genes family and analyse their associations with variation in meat production traits in Czech pig breeds. To verify the associations between the polymorphisms and the selected meat traits, altogether 254 pigs, including full- and half-sibs, of Large White and Landrace breeds were tested. The studied meat characteristics were weight of neck, loin, shoulder and ham, lean meat content (LMC), backfat thickness, intramuscular fat (IMF), remission, dry matter content and test daily gain. Statistically significant associations were observed between MYOG gene and fat and neck weight, and between MYF5 gene and IMF and LMC. High significant differences were observed between genotypes AA and AB of MYOD1 in IMF and between genotypes AB and BB of MYF5 in loin weight.     

A genomewide association mapping study using ultrasound‐scanned information identifies potential genomic regions and candidate genes affecting carcass traits in Nellore cattle

The aim of this study was to identify candidate genes and genomic regions associated with ultrasound‐derived measurements of the rib‐eye area (REA), backfat thickness (BFT) and rumpfat thickness (RFT) in Nellore cattle. Data from 640 Nellore steers and young bulls with genotypes for 290 863 single nucleotide polymorphisms (SNPs) were used for genomewide association mapping. Significant SNP associations were explored to find possible candidate genes related to physiological processes. Several of the significant markers detected were mapped onto functional candidate genes including ARFGAP3, CLSTN2 and DPYD for REA; OSBPL3 and SUDS3 for BFT; and RARRES1 and VEPH1 for RFT. The physiological pathway related to lipid metabolism (CLSTN2, OSBPL3, RARRES1 and VEPH1) was identified. The significant markers within previously reported QTLs reinforce the importance of the genomic regions, and the other loci offer candidate genes that have not been related to carcass traits in previous investigations.     

Whole‐genome resequencing to identify candidate genes for the QTL for oleic acid percentage in Japanese Black cattle

In our previous study, we detected a QTL for the oleic acid percentage (C18:1) on BTA9 in Japanese Black cattle through a genome‐wide association study (GWAS). In this study, we performed whole‐genome resequencing on eight animals with higher and lower C18:1 to identify candidate polymorphisms for the QTL. A total of 39,658 polymorphisms were detected in the candidate region, which were narrowed to 1993 polymorphisms within 23 genes based on allele differences between the high and low C18:1 groups. We subsequently selected three candidate genes, that is, CYB5R4, MED23, and VNN1, among the 23 genes based on their function in fatty acid metabolism. In each candidate gene, three SNPs, that is, CYB5R4 c.*349G > T, MED23 c.3700G > A, and VNN1 c.197C > T, were selected as candidate SNPs to verify their effect on C18:1 in a Japanese Black cattle population (n = 889). The statistical analysis showed that these SNPs were significantly associated with C18:1 (p < 0.05), suggesting that they were candidates for the QTL. In conclusion, we successfully narrowed the candidates for the QTL by detecting possible polymorphisms located within the candidate region. It is expected that the responsible polymorphism can be identified by demonstrating their effect on the gene's function.     

北京黑猪NR6A1、VSX2、VRTN、LTBP2基因多态性与脊椎数及胴体性状的关联分析

旨在研究NR6A1、VSX2、VRTN、LTBP2基因在北京黑猪群体中基因的多态性与脊椎数和胴体性状(胴体长、胴体直长、胴体斜长及胴体重)的关联。本研究提取410头(210±7)日龄健康北京黑猪耳组织DNA,通过PCR技术和Sanger测序技术对4个基因外显子区进行基因分型,计算变异位点基因型频率以及等位基因频率的分布,并与脊椎数和胴体性状进行关联分析。结果表明,在VSX2、VRTN、LTBP2三个基因中共筛选到15个SNPs,其中包含1个CDS区的同义突变和14个3'UTR突变。使用Duncan's多重检验统计分析发现,VSX2基因中CDS区的1个突变c.536 G＞A与脊椎数关联显著(P＜0.05)。VRTN中3'UTR区共4个突变与性状关联显著(P＜0.05),其中c.2598 A＞G和c.2607 G＞T与脊椎数、胴体长、胴体直长及胴体斜长关联显著,c.3087 G＞C只与胴体斜长关联显著,c.3094 A＞G与脊椎数及胴体斜长关联显著(P＜0.05)。LTBP2上6个3'UTR突变位点(c.7158 A＞G、c.6869 C＞T、c.6319 G＞C、c.6246 G＞T、c.6170 A＞G、c.6079 G＞T)与脊椎数关联显著(P＜0.05),但与胴体性状关联不显著。综上所述,VSX2、VRTN、LTBP2基因中共有10个SNPs与脊椎数性状关联显著(P＜0.05),仅VRTN基因中4个SNPs与胴体性状关联显著(P＜0.05)。这些关联显著的SNPs可以作为北京黑猪脊椎数及胴体性状变异的候选基因功能位点。     

Carcass characteristics of Hanwoo (Korean cattle) from different sex conditions, raising altitudes and slaughter seasons

This study was carried out to investigate the carcass characteristics of Hanwoo (Korean cattle) from different sex conditions, raising altitudes and slaughter seasons. The total number of cattle used in this study was 3608 heads which comprised of 1336 bulls and 1660 steers. The data was analyzed according to sex condition (bull and steer), raising altitude (low: 0–100 m and high: 700–800 m) and slaughter season (spring, summer, autumn and winter). At 24 h post-slaughter, the carcasses were weighed and evaluated for carcass traits according to the Korean carcass grading standard by an official grader. Carcass weight, ribeye area, yield index and grade, and meat color, firmness and maturity scores of carcass from bulls were significantly higher (P < 0.001) than those from steers. Inversely, the backfat thickness, marbling score and quality grade of carcass from steers were significantly higher (P < 0.001) than those from bulls. The maturity score of carcass from a high area was significantly higher (P < 0.01) than that from a low area. The ribeye area of carcass from summer was significantly smaller (P < 0.05) than that from other seasons. The maturity score of carcass from the summer season was highest among that from other seasons and the lowest was that from the winter season. Marbling score and carcass quality grade from the winter season were highest among those from other seasons and the lowest were those from the autumn season. There were no significant interactions between sex condition and raising altitude on carcass traits except the ribeye area. There were significant interactions between sex condition and slaughter season on marbling score and carcass quality grade. There were no significant interactions between raising altitude and slaughter season on all of the carcass traits. There were significant interactions among sex condition, raising altitude and slaughter season on meat color score. It was concluded that sex condition affected muscle and fat depositions on the carcass, raising altitude affected maturity and slaughter season affected ribeye area, maturity and marbling.     

Variations in genes involved in fat metabolism and their association with ultrasonic and carcass traits in Japanese Black steers

Polymorphisms in genes involved in lipid metabolism have been reported to be associated with fatty acid composition of adipose tissue. However, the relationship of these polymorphisms with premortem ultrasonic traits is unknown. The objective of this study, therefore, was to assess the association between polymorphisms in fatty acid synthase, stearoyl‐coenzyme A desaturase (SCD), sterol regulatory element‐binding protein 1 (SREBP1), diacylglycerol acyltransferase 1, and nuclear receptor subfamily 1, group H, number 3 genes with ultrasonic and carcass traits in Japanese Black steers (n = 300) under progeny testing at the Livestock Improvement Association of Miyazaki. To have a comprehensive analysis of the association between the aforementioned genetic polymorphisms and ultrasonic traits, longitudinal measurements of ultrasonic traits were taken. Furthermore, the association of these genetic polymorphisms and carcass traits was evaluated. The polymorphisms in the SCD gene and SREBP1 were associated (P < 0.05) with some ultrasonic traits at multiple stages. To add to that, the polymorphisms were associated (P < 0.05) with some carcass traits. These findings suggest that the polymorphisms in SCD and SREBP1 are functional mutations or could be related to mutations that can aid in selection to improve some ultrasonic and carcass traits.     

Signatures of selection in Charolais beef cattle identified by genome‐wide analysis

Charolais cattle are one of the most important breeds for meat production worldwide; in México, its selection is mainly made by live weight traits. One strategy for mapping important genomic regions that might influence productive traits is the identification of signatures of selection. This type of genomic features contains loci with extended linkage disequilibrium (LD) and homozygosity patterns that are commonly associated with sites of quantitative trait locus (QTL). Therefore, the objective of this study was to identify the signatures of selection in Charolais cattle genotyped with the GeneSeek Genomic Profiler Bovine HD panel consisting of 77 K single nucleotide polymorphisms (SNPs). A total 61,311 SNPs and 819 samples were used for the analysis. Identification of signatures of selection was carried out using the integrated haplotype score (iHS) methodology implemented in the rehh R package. The top ten SNPs with the highest piHS values were located on BTA 4, 5, 6 and 14. By identifying markers in LD with top ten SNPs, the candidate regions defined were mapped to 52.8–59.3 Mb on BTA 4; 67.5–69.3 on BTA 5; 39.5–41.0 Mb on BTA 6; and 26.4–29.6 Mb on BTA 14. The comparison of these candidate regions with the bovine QTLdb effectively confirmed the association (p < 0.05) with QTL related to growth traits and other important productive traits. The genomic regions identified in this study indicated selection for growth traits on the Charolais population via the conservation of haplotypes on various chromosomes. These genomic regions and their associated genes could serve as the basis for haplotype association studies and for the identification of causal genes related to growth traits.     

Distribution of IS900 restriction fragment length polymorphism types among animal Mycobacterium avium subsp. paratuberculosis isolates from Argentina and Europe

Sixty-one Mycobacterium avium subsp. paratuberculosis isolates from cattle and deer from the Buenos Aires province, an important livestock region in Argentina, were typed by restriction fragment length polymorphisms (RFLP) analysis based on IS900. Four different RFLP patterns (designated ‘A’, ‘B’, ‘C’ and ‘E’) were identified in BstEII digests of genomic DNA. The most frequently observed type, pattern ‘A’, was found in 46 isolates (75%). The second, pattern ‘E’, included 8 isolates (13%), while the third, pattern ‘B’, included 6 isolates (10%). Pattern ‘C’ was found for only one isolate. All of the deer isolates were classified as pattern ‘A’, while cattle isolates represented all four RFLP patterns. Twenty-one isolates representing the four different BstEII-RFLP patterns were digested with PstI. Twenty isolates showed identical PstI-RFLP pattern. BstEII-RFLP patterns from Argentine cattle and deer were compared with patterns found in cattle, goat, deer, rabbit, and human isolates from Europe. The most common pattern in Argentina, pattern ‘A’, was identical to a less frequently occurring pattern R9 (C17) from Europe. The other Argentine patterns ‘B’, ‘C’ and ‘E’, were not found in the Europe. These results indicate that the distribution of M. avium subsp. paratuberculosis genotypes in the Buenos Aires province of Argentina is different from that found in Europe.     

Polymorphisms in the promoter region of myostatin gene are associated with carcass traits in pigs

Higher average daily gain, more lean meat yield and less fat yield of porcine carcass increase selling profits for animal producers. Myostatin (MSTN), previously called GDF8, is a member of transforming growth factor‐β (TGF‐β) superfamily. It is a negative regulator for both embryonic development and adult homeostasis of skeletal muscle. In this study, the genotypes of the previously described SNPs MSTN g.435G>A and g.447A>G SNPs in 66 Duroc pigs, 33 Landrace pigs, 180 Duroc × Landrace (DL) pigs and 155 Duroc × Yorkshire × Landrace (DYL) pigs were determined by Taqman SNP Genotyping Assays. For Duroc and Landrace pigs, MSTN g.435GG/g.447AA individual had greater backfat thickness (p < 0.05) than g.435AA/g.447GG individual, whereas MSTN g.435AA/g.447GG had greater meat (p < 0.05) and meat percentage (p < 0.05) than g.435GA/g.447AG individual. For DL and DYL pigs, the MSTN g.435GG/g.447AA animals were greater in backfat at ultrasound 10th rib (p < 0.05) and carcass 10th rib (p < 0.01) than g.435AA/g.447GG individual. The MSTN g.435AA/g.447GG individual also had higher values than g.435GG/g.447AA for anterior‐end meat (p < 0.05), posterior‐end meat (p < 0.01), total meat weight (p < 0.01) and meat percentage (p < 0.01). This study confirmed evidence that MSTN g.435G>A and g.447A>G affected carcass traits in pigs. The effects of the mutated alleles were additive with the maximal effects resulting from two copies of the mutated allele. Selection for MSTN g.435A/g.447G allele is expected to increase muscle of limb and total meat production and decrease backfat thickness.     

Single nucleotide polymorphisms associated with thermoregulation in lactating dairy cows exposed to heat stress

Dairy cows with increased rectal temperature experience lower milk yield and fertility. Rectal temperature during heat stress is heritable, so genetic selection for body temperature regulation could reduce effects of heat stress on production. One aim of the study was to validate the relationship between genotype and heat tolerance for single nucleotide polymorphisms (SNPs) previously associated with resistance to heat stress. A second aim was to identify new SNPs associated with heat stress resistance. Thermotolerance was assessed in lactating Holsteins during the summer by measuring rectal temperature (a direct measurement of body temperature regulation; n = 435), respiration rate (an indirect measurement of body temperature regulation, n = 450) and sweating rate (the major evaporative cooling mechanism in cattle, n = 455). The association between genotype and thermotolerance was evaluated for 19 SNPs previously associated with rectal temperature from a genomewide analysis study (GWAS), four SNPs previously associated with change in milk yield during heat stress from GWAS, 2 candidate gene SNPs previously associated with rectal temperature and respiration rate during heat stress (ATPA1A and HSP70A) and 66 SNPs in genes previously shown to be associated with reproduction, production or health traits in Holsteins. For SNPs previously associated with heat tolerance, regions of BTA4, BTA6 and BTA24 were associated with rectal temperature; regions of BTA6 and BTA24 were associated with respiration rate; and regions of BTA5, BTA26 and BTA29 were associated with sweating rate. New SNPs were identified for rectal temperature (n = 12), respiration rate (n = 8) and sweating rate (n = 3) from among those previously associated with production, reproduction or health traits. The SNP that explained the most variation were PGR and ASL for rectal temperature, ACAT2 and HSD17B7 for respiration rate, and ARL6IP1 and SERPINE2 for sweating rate. ARL6IP1 was associated with all three thermotolerance traits. In conclusion, specific genetic markers responsible for genetic variation in thermoregulation during heat stress in Holsteins were identified. These markers may prove useful in genetic selection for heat tolerance in Holstein cattle.     

中国西门塔尔牛血管生成素样蛋白3基因多态性及其与肉质性状的关联性分析

本研究旨在分析血管生成素样蛋白3(angiopoietin-like protein 3,ANGPTL3)基因的单核苷酸多态性(SNPs)及其与中国西门塔尔牛肉质性状的关联性,寻找可用于辅助选择的分子标记。随机选取饲养条件相同的98头中国西门塔尔牛为研究对象,采用PCR-SSCP技术结合DNA测序技术检测ANGPTL3基因的SNP。运用SPSS 19.0软件对ANGPTL3基因的不同基因型与肉质性状进行关联性分析。经测序发现,ANGPTL3基因外显子4中存在1个突变位点G7358C,存在2种基因型:CD、DD。关联性分析结果表明,G7358C突变位点的不同基因型与胴体脂肪覆盖率、眼肌面积、肌内脂肪和大理石花纹差异显著相关(P0.05),CC基因型个体胴体脂肪覆盖率、眼肌面积、肌内脂肪和大理石花纹显著高于CD基因型个体(P0.05)。试验结果表明,G7358C位点的CC基因型为优势基因型,与胴体脂肪覆盖率、眼肌面积、肌内脂肪和大理石花纹等肉用性状有相关性,ANGPTL3基因有望作为分子标记辅助选择的候选基因。     

Identification of genomic regions that exhibit sexual dimorphism for size and muscularity in cattle

Sexual dimorphism, the phenomenon whereby males and females of the same species are distinctive in some aspect of appearance or size, has previously been documented in cattle for traits such as growth rate and carcass merit using a quantitative genetics approach. No previous study in cattle has attempted to document sexual dimorphism at a genome level; therefore, the objective of the present study was to determine whether genomic regions associated with size and muscularity in cattle exhibited signs of sexual dimorphism. Analyses were undertaken on 10 linear-type traits that describe the muscular and skeletal characteristics of both males and females of five beef cattle breeds: 1,444 Angus (AA), 6,433 Charolais (CH), 1,129 Hereford, 8,745 Limousin (LM), and 1,698 Simmental. Genome-wide association analyses were undertaken using imputed whole-genome sequence data for each sex separately by breed. For each single-nucleotide polymorphism (SNP) that was segregating in both sexes, the difference between the allele substitution effect sizes for each sex, in each breed separately, was calculated. Suggestively (P ≤ 1 × 10⁻⁵) sexually dimorphic SNPs that were segregating in both males and females were detected for all traits in all breeds, although the location of these SNPs differed by both trait and breed. Significantly (P ≤ 1 × 10⁻⁸) dimorphic SNPs were detected in just three traits in the AA, seven traits in the CH, and three traits in the LM. The vast majority of all segregating autosomal SNPs (86% in AA to 94% in LM) had the same minor allele in both males and females. Differences (P ≤ 0.05) in allele frequencies between the sexes were observed for between 36% (LM) and 66% (AA) of the total autosomal SNPs that were segregating in both sexes. Dimorphic SNPs were located within a number of genes related to muscularity and/or size including the NAB1, COL5A2, and IWS1 genes on BTA2 that are located close to, and thought to be co-inherited with, the MSTN gene. Overall, sexual dimorphism exists in cattle at the genome level, but it is not consistent by either trait or breed.     

Polymorphisms of the PNPLA3 gene and their associations with chicken growth and carcass traits

1. An F₂ resource population of Gushi chickens crossed with Anka broilers was used to investigate the genetic effects of the chicken PNPLA3 gene on growth and adipose accumulation.

2. Associations between three SNPs (g.40006G?>?T, g.42344T?>?C and g.42404A?>?T) and broiler traits were determined using linkage disequilibrium, haplotype construction and association analysis.

3. The g.40006G?>?T mutation was associated with body weights at 4, 6, 8, 10 and 12 weeks of age, carcass weight, evisceration weight and semi-evisceration weight (P? 4. Haplotypes of the g.42344T?>?C and g.42404A?>?T mutations were associated with body weight at 12 weeks, carcass weight, evisceration weight, and semi-evisceration weight (P? 5. The results suggest that the PNPLA3 gene may be in linkage with the causative mutation or a QTL controlling growth traits in chickens. In contrast to human studies, the polymorphisms were not associated with fat related traits.     

Genetic parameter estimation and gene association analyses for meat quality traits in open-air free-range Iberian pigs

Meat quality of Iberian pigs is defined by the combination of their genetic characteristics and the particular production system. To carry out a genetic analysis of the main meat quality traits, we estimated their heritabilities, genetic correlations and the association effects of 32 selected SNPs of 12 candidate genes. A total of ten traits were measured in longissimus dorsi samples from 1,199 Iberian pigs fattened in the traditional free-range system: water holding capacity (thawing, cooking and centrifuge force water losses), instrumental colour (lightness L*, redness a* and yellowness b*), myoglobin content, shear force on cooked meat, and shear force and maximum compression force on dry-cured loin. Estimated heritability values were low to moderate (0.01 to 0.43) being the lowest for L* and the highest for cooking loss. Strong genetic correlations between water holding capacity traits (0.93 to 0.96) and between myoglobin content and a* (0.94) were observed. The association analyses revealed 19 SNPs significantly associated with different traits. Consistent and strong effects were observed between PRKAG3 SNPs (rs319678464G > C and rs330427832C > T), MYH3_rs81437544T > C, CASP3_rs319658214G > T and CTSL_rs332171512A > G and water losses. Also for CAPN1_rs81358667G > A and CASP3_rs319658214G > T and shear force. The SNPs mapping on PRKAG3 showed the highest effects on Minolta colour traits. Genotyping of these SNPs could be useful for the selection of Iberian young boars with similar estimated breeding values for productive traits.     

Genome wide association study on beef production traits in Marchigiana cattle breed

A genomewide association study was carried out on a sample of Marchigiana breed cattle to detect markers significantly associated with carcass and meat traits. Four hundred and nine young bulls from 117 commercial herds were genotyped by Illumina 50K BeadChip assay. Eight growth and carcass traits (average daily gain, carcass weight, dressing percentage, body weight, skin weight, shank circumference, head weight and carcass conformation) and two meat quality traits (pH at slaughter and pH 24 h after slaughter) were measured. Data were analysed with a linear mixed model that included fixed effects of herd, slaughter date, fixed covariables of age at slaughter and SNP genotype, and random effects of herd and animal. A permutation test was performed to correct SNP genotype significance level for multiple testing. A total of 96 SNPs were significantly associated at genomewide level with one or more of the considered traits. Gene search was performed on genomic regions identified on the basis of significant SNP position and level of linkage disequilibrium. Interesting loci affecting lipid metabolism (SOAT1), bone (BMP4) and muscle (MYOF) biology were highlighted. These results may be useful to better understand the genetic architecture of growth and body composition in cattle.     

Three SNPs within exons of INHA and ACVR2B genes are significantly associated with litter size in Dazu black goats

The aim of this study was to analyse the association between single-nucleotide polymorphisms within INHA and ACVR2B and litter size in Dazu black goats. In total, twenty-two SNPs were genotyped in 190 individuals by SNaPshot and resequencing. The results showed that three SNPs (SNP_1, SNP_12 and SNP_13 in this study) were detected to have significant additive genetic effect on the recorded goat litter size (p < .05). The SNP_1 (NC_030809.1), a non-synonymous substitution of G for T at chr2-g. 28314990 in the exon 2 of INHA gene (NM_001285606.1), resulted in homozygote 2 (HOM2) contributed 0.25 and heterozygote (HET) contributed 0.12 larger litter than homozygote 1 (HOM1). Meanwhile, SNP_12 (Chr22-g. 11721225 A > T) and SNP_13 (Chr22-g. 11721227 A > C) (NC_030829.1) simultaneously mutated at the first and third position of a triplet AAA (lysine, K) in the exon 4 of ACVR2B gene (XM_018066623.1) had estimated genetic effects of HOM1 (0.00) and HOM2 (0.03) larger than HET (−0.12). In conclusion, one SNPs (chr2-g. 28314990 T > G) within the exon 2 of INHA and two SNPs (Chr22-g. 11721225 A > T and Chr22-g. 11721227 A > C) i n the exon 4 of ACVR2B gene were highly recommended as candidate markers of litter size in Dazu black goats. A large-scale association study to assess the impact of these variants on litter size is still necessary.