Samples and traits
DNA samples were isolated from blood that has collected from hundred and three standing female animals of Indian origin maintained at Livestock Research Complex, ICAR-National Dairy Research Institute, India under the same management regimen which completed their first lactation. Around 10 ml of venous blood was collected aseptically from the jugular vein of the animals in a 15 ml polypropylene centrifuge tube under sterile condition using 0.5 ml of EDTA as an anticoagulant. The samples were stored at 20°C for further utilization. The different traits considered for association analysis included
viz. body weight at different age
i.e. Birth, 6 Months (6M), 12 Months (12M), 18 Months (18M), 24 Months (24M), 30 Months (30M), Weight at First Service (WFS), Weight at First Calving (WFC), reproduction
i.e. Age at First Service (AFS), Age at First Calving (AFC), First Service Period (FSP), First Calving Interval (FCI) and first lactation production traits
i.e. First Lactation 305 Day or less Milk Yield (FL305DMY), First Lactation Total Milk Yield (FLTMY), First Lactation Length (FLL), First Lactation 305 day or less Fat Yield (FL305DFY), First Lactation 305 day or less SNF Yield (FL305DSNFY).
DNA extraction and genotyping of animal
The DNA was extracted from blood samples by commercially available kit [Genomic DNA Purification Kit (Cat.#A1620, Promega)]. Three regions
i.e. partial intron 1, exon 2 and partial intron 2 (454 bp), partial intron 2, exon 3 and partial intron 3 (437 bp) and partial intron 3 and exon 4 (277 bp) of NPY were amplified. Amplification of the desired NPY gene fragments was performed with designed primer pairs (Table 1).
PCR was performed in 25-μl reactions with a PCR mix containing 15 mM MgCl
2, 200 μM dNTP, 30 pM of each primer, 100 ng of genomic DNA and 1 U of Taq polymerase. Thermal cycling conditions were included an initial denaturation for 5 min at 95°C followed by 32 cycles of 95°C (30 s), annealing temperature 64°C, 59°C and 63°C (30 s) and ending with a final extension for 10 min at 72°C. The PCR products were confirmed by agarose gel electro- phoresis methods using 1.5% agarose gel with ethidium bromide for 20-30 min and visualized under Gel Documentation System.
The sequencing of the targeted regions was done to screen the DNA polymorphisms of the NPY gene. Sequence data were analyzed using multiple sequence alignment software, freely available online. The coding regions of NPY gene were subjected to basic local alignment search tool at NCBI (www.ncbi.nlm.org) to identify the sequence homology with the corresponding region of other species available with the data bank. Sequence data from different variants obtained were analyzed to identify the SNP variation present in the Sahiwal cattle maintained at ICAR-NDRI, Karnal.
Statistical analysis
Genotypic and allelic frequencies
Genotypic frequencies for variant genotypes were calculated using the formula:
Allelic frequencies were calculated as follows:
Allelic frequency of A = AA + ½ AB
Allelic frequency of B = BB + ½ AB
Where
AA and BB = Genotypic frequency of homozygote
AB = Genotypic frequency of heterozygote
A and B = Allelic frequencies
Least squares analysis of variance was carried out for unequal and non-orthogonal data using the technique described by Harvey (1990) to study effect of non-genetic factors like season and period. The model was used with assumptions that different components being fitted into the model were linear, independent and additive.
Association analysis
Association of the animal genotypes with growth, reproduction and milk production traits were determined by analyses of quantitative traits. Statistical analysis was carried out using SAS Enterprise Guide 4.2 software (SAS Institute Inc. 2009) and significance of differences based on genotypes effect of growth, reproduction and production traits were tested by following general linear models:
Yij= µ + Gi + eij
Where
Y
ij is the trait value of j
th animal of i
th genotype, μ is the overall
mean, G
i is the effect of i
th genotypes and e
ij is the residual error NID (0, σ
2e).