Ethics statement
Collection of biological samples and experimental procedures involved in this study were approved by the Institutional Animal Care and Use Committee in Sichuan Animal Science Academy, China.
Animals and body size traits
The experiment was carried out from 2019-05 to 2020-4 at the Sichuan Animal Science Academy. A total of 169 healthy female Shuxuan cattle with more than 12 months of age were involved in the present study. According to pedigree records, the genetic relationships within three generations were avoided for all animals. They were subjected to a collection of jugular venous blood samples and measurements of body size traits. We totally measured and collected 11 body size traits, including the withers height (WH, from the highest point of shoulder blade to the ground), body length (BL, from shoulder joint to tuber ischiadicum), head length (HL, from poll to nose), forehead width (FW, maximum distance between the temples), chest width (HW, height immediately posterior to the front leg), chest depth (CD, width immediately posterior to the front leg), chest girth (CG, body circumference immediately posterior to the front leg), shin circumference (SC, circumference of left metacarpus at its narrowest), rump width (RW, distance between tuber ischiadicum), hip width (PW, distance between the lateral surfaces of the tuber coxae) and hip height (HH, from hip to the ground). All body size traits were measured when the animal was standing naturally.
Genotyping of MTNR1A gene
The candidate SNP (c.A455G with GenBank accession number of rs208053155) was selected from the former publication (Yang
et al., 2015) and subjected to genotyping in the present study. Genomic DNA was extracted using AxyPrep Genomic DNA Miniprep Kit (Axygen Bioscience, USA) and then stored at -20°C for later analysis. This candidate SNP was genotyped using a custom-by-design SNPscan
TM Kit (Genesky Biotechnologies Inc., Shanghai, China). Briefly, 100-200 ng genomic DNA was used for producing 20 μL ligation mixture containing 1X Ligase buffer, 1X probe and 0.5 μL ligase. The ligation reaction was performed on a ABI 2720 Thermal Cycler (Applied Biosystems, USA) under the following program: 2 min at 98°C, 4 cycles of (94°C 1 min, 58°C 4 hr), 94°C 2 min, hold at 4°C and immediately stopped by adding 20μl of 2xStop Buffer. Multiplex fluorescence PCR reactions were performed for each ligation product. Each PCR reaction was prepared in a 20 μL mixture containing 1xPCR Master Mix, 1 μL Primer Mix SetA or SetB and 1 μL ligation product. PCR products were separated and detected by capillary electrophoresis on ABI3730XL sequencer (Applied Biosystems, USA).
Data analysis
We first analyzed the overall distribution and pairwise relationships among the 11 body size traits using the R package of Performance Analytics, in which the spearman correlation coefficient was derived. After both the genotype and allele frequencies of SNP were calculated, the Hardy-Weinberg equilibrium (HWE) was measured through the chi-square test (c
2 test). The polymorphism information content (PIC) and heterozygosity (H
e) were calculated. The associations between genotypes and body size traits were investigated under the linear mixed model (LMM) as:
Where
Y
ijk is measurement values of body size traits, μ is overall of mean for each trait, G
i is the effect of genotypes (AA, AG and GG); Agej is the fixed effect of animal age (having three levels of 12 ~ 16 months, 17 ~ 23 months and > 24 months); e
ijk is random residual effects. The model fitting and statistical tests were conducted using R package of lmerTest
(Kuznetsova et al., 2017). The Bonferroni method was used for adjustment for multiple comparisons.