Sex-specific genetic parameter estimates of body weight in Mazandaran indigenous chickens

Body weight is an economically important trait in poultry that shows sexual dimorphism (SD). In the present study, variation in SD in Mazandaran native chickens was investigated in terms of the (Co) variance components and genetic parameters of body weight between males and females. Studied traits were body weights at hatch (BW1), 8 weeks (BW8) and 12 weeks of age (BW12). Also, for weight at sexual maturity (WSM) covariance components were only estimated in females. Cross-sex direct and maternal correlations were also estimated for studied traits except for WSM. For this purpose, a deep 21-generation pedigree and body weight data (57,576 BW1, 72,925 BW8, 62,727 BW12 and, 42,496 WSM) were used. Evaluation of SD of body weight was performed using six bivariate animal models with and without considering the genetic and permanent maternal environmental effects under the restricted maximum likelihood method in WOMBAT software. Model with direct additive genetic effects and maternal genetic effects without covariance between them was identified as the best model for BW1 and BW8. The Model including direct additive genetic effects and permanent maternal environmental effects was the best model for BW12 and WSM. Direct heritability (h2 ) estimates for BW1, BW8 and, BW12 were, respectively, 0.05 ± 0.013, 0.17 ± 0.02 and, 0.25 ± 0.03 in males and, 0.05 ± 0.012, 0.15 ± 0.01 and 0.21 ± 0.01 in females. Also, the direct heritability of WSM based on univariate analysis in females was estimated to be 0.40 ± 0.01. Maternal heritability (h m 2 $$ {h}_m^2 $$ ) varied from 0.39 ± 0.01 (BW1) to 0.04 ± 0.009 (BW8) in males, and 0.36 ± 0.10 (BW1) to 0.04 ± 0.006 (BW8) in females. The correlation between direct genetic effects between males and females for BW1 was not significantly different from one. The direct genetic correlation between the two sexes for BW8 and BW12 was significantly different from 1 concluding that these traits are dimorphic in terms of direct genetic effects and therefore independent selection in both sexes is possible.

Genome-wide study of linkage disequilibrium, population structure, and inbreeding in Iranian indigenous sheep breeds

Knowledge of linkage disequilibrium (LD), genetic structure and genetic diversity are some key parameters to study the breeding history of indigenous small ruminants. In this study, the OvineSNP50 Bead Chip array was used to estimate and compare LD, genetic diversity, effective population size (Ne) and genomic inbreeding in 186 individuals, from three Iranian indigenous sheep breeds consisting of Baluchi (n = 96), Lori-Bakhtiari (n = 47) and Zel (n = 47). The results of principal component analysis (PCA) revealed that all animals were allocated to the groups that they sampled and the admixture analysis revealed that the structure within the populations is best explained when separated into three groups (K  =  3). The average r2 values estimated between adjacent single nucleotide polymorphisms (SNPs) at distances up to 10Kb, were 0.388±0.324, 0.353±0.311, and 0.333±0.309 for Baluchi, Lori-Bakhtiari and Zel, respectively. Estimation of genetic diversity and effective population size (Ne) showed that the Zel breed had the highest heterozygosity and Ne, whereas the lowest value was found in Baluchi breed. Estimation of genomic inbreeding using FROH (based on the long stretches of consecutive homozygous genotypes) showed the highest inbreeding coefficient in Baluchi and the lowest in Zel breed that could be due to higher pressure of artificial selection on Baluchi breed. The results of genomic inbreeding and Ne showed an increase in sharing haplotypes in Baluchi, leading to the enlargement of LD and the consequences of linkage disequilibrium and haplotype blocks confirmed this point. Also, the persistence of the LD phase between Zel and Lori-Bakhtiari was highest indicating that these two breeds would be combined in a multi-breed training population in genomic selection studies.

Autosomal and Z-linked genetic evaluation for body weight in Mazandaran native chicken using different models for dosage compensation on the Z chromosome

The objectives of this study were to assess the autosomal and sex-linked genetic inheritance of growth traits and identify the effective dosage compensation on the Z chromosome in Mazandaran native chickens. The data included body weights at hatching (BW0), 8 weeks (BW8) and 12 weeks (BW12) of age, related to the first 21 generations of selection, were collected from Mazandaran native chicken breeding centre. The fixed effects included sex of birds in two classes, hatch in five classes and generation in 21 classes. The inverse of the sex-linked additive genetic relationship matrix was constructed using nadiv package in R considering different models for dosage compensation on the Z chromosome. The setup inversed matrix was then supplied externally to WOMBAT using the GIN option. Twelve univariate animal models separating participation of autosomal additive genetic, sex-linked additive genetic and maternal effects (both genetic and permanent environment effects) with considering the five different dosage compensation methods for models with sex-linked effects were analysed by WOMBAT software. BW0 was not affected by sex-linked additive genetic effects. For BW8 and BW12 the model which included autosomal, sex-linked direct additive and maternal effects with no global dosage compensation for the Z chromosome was the most appropriate model. Autosomal heritability estimates were 0.05 ± 0.02, 0.10 ± 0.01 and 0.11 ± 0.01, for BW0, BW8 and BW12, respectively. For BW8 and BW12, sex-linked heritability estimates were 0.07 and 0.27, respectively. Spearman rank correlation coefficient between autosomal and sex-linked breeding values were 0.45 and 0.12 for BW8 and BW12, respectively. Spearman rank correlation coefficient between autosomal and sex-linked breeding values were 0.45 and 0.12 for BW8 and BW12, respectively. The autosomal direct additive genetic correlations between all traits were positive. The estimate of direct sex-linked additive genetic correlation between BW8 and BW12 was high (0.88). Also, maternal genetic correlations were 0.53, 0.54 and 0.91 between BW0-BW8, BW0-BW12 and BW8-BW12, respectively. Given the importance of Z-linked genes for BW8 and BW12, it is recommended that Z-linked effects be separated from autosomal effects in order to increase the accuracy of genetic evaluation of birds for these traits.

Genome-wide identification of copy number variation and association with fat deposition in thin and fat-tailed sheep breeds

Copy number variants (CNVs) are a type of genetic polymorphism which contribute to phenotypic variation in several species, including livestock. In this study, we used genomic data of 192 animals from 3 Iranian sheep breeds including 96 Baluchi sheep and 47 Lori-Bakhtiari sheep as fat-tailed breeds and 47 Zel sheep as thin-tailed sheep breed genotyped with Illumina OvineSNP50K Beadchip arrays. Also, for association test, 70 samples of Valle del Belice sheep were added to the association test as thin-tailed sheep breed. PennCNV and CNVRuler software were, respectively, used to study the copy number variation and genomic association analyses. We detected 573 and 242 CNVs in the fat and thin tailed breeds, respectively. In terms of CNV regions (CNVRs), these represented 328 and 187 CNVRs that were within or overlapping with 790 known Ovine genes. The CNVRs covered approximately 73.85 Mb of the sheep genome with average length 146.88 kb, and corresponded to 2.6% of the autosomal genome sequence. Five CNVRs were randomly chosen for validation, of which 4 were experimentally confirmed using Real time qPCR. Functional enrichment analysis showed that genes harbouring CNVs in thin-tailed sheep were involved in the adaptive immune response, regulation of reactive oxygen species biosynthetic process and response to starvation. In fat-tailed breeds these genes were involved in cellular protein modification process, regulation of heart rate, intestinal absorption, olfactory receptor activity and ATP binding. Association test identified one copy gained CNVR on chromosomes 6 harbouring two protein-coding genes HGFAC and LRPAP1. Our findings provide information about genomic structural changes and their association to the interested traits including fat deposition and environmental compatibility in sheep.

Meta-analysis of genome-wide association studies for litter size in sheep

Litter size and ovulation rate are important reproduction traits in sheep and have important impacts on the profitability of farm animals. To investigate the genetic architecture of litter size, we report the first meta-analysis of genome-wide association studies (GWAS) using 522 ewes and 564,377 SNPs from six sheep breeds. We identified 29 significant associations for litter size which 27 of which have not been reported in individual GWAS for each population. However, we could confirm the role of BMPR1B in prolificacy. Our gene set analysis discovered biological pathways related to cell signaling, communication, and adhesion. Functional clustering and enrichment using protein databases identified epidermal growth factor-like domain affecting litter size. Through analyzing protein-protein interaction data, we could identify hub genes like CASK, PLCB4, RPTOR, GRIA2, and PLCB1 that were enriched in most of the significant pathways. These genes have a role in cell proliferation, cell adhesion, cell growth and survival, and autophagy. Notably, identified SNPs were scattered on several different chromosomes implying different genetic mechanisms underlying variation of prolificacy in each breed. Given the different layers that make up the follicles and the need for communication and transfer of hormones and nutrients through these layers to the oocyte, the significance of pathways related to cell signaling and communication seems logical. Our results provide genetic insights into the litter size variation in different sheep breeds.

Genes and Pathways Affecting Sheep Productivity Traits: Genetic Parameters, Genome-Wide Association Mapping, and Pathway Enrichment Analysis

Ewe productivity is a composite and maternal trait that is considered the most important economic trait in sheep meat production. The objective of this study was the application of alternative genome-wide association study (GWAS) approaches followed by gene set enrichment analysis (GSEA) on the ewes’ genome to identify genes affecting pregnancy outcomes and lamb growth after parturition in Iranian Baluchi sheep. Three maternal composite traits at birth and weaning were considered. The traits were progeny birth weight, litter mean weight at birth, total litter weight at birth, progeny weaning weight, litter mean weight at weaning, and total litter weight at weaning. GWASs were performed on original phenotypes as well as on estimated breeding values. The significant SNPs associated with composite traits at birth were located within or near genes RDX, FDX1, ARHGAP20, ZC3H12C, THBS1, and EPG5. Identified genes and pathways have functions related to pregnancy, such as autophagy in the placenta, progesterone production by the placenta, placental formation, calcium ion transport, and maternal immune response. For composite traits at weaning, genes (NR2C1, VEZT, HSD17B4, RSU1, CUBN, VIM, PRLR, and FTH1) and pathways affecting feed intake and food conservation, development of mammary glands cytoskeleton structure, and production of milk components like fatty acids, proteins, and vitamin B-12, were identified. The results show that calcium ion transport during pregnancy and feeding lambs by milk after parturition can have the greatest impact on weight gain as compared to other effects of maternal origin.

Performance evaluation of various evapotranspiration modeling scenarios based on METRIC method and climatic indexes

Evapotranspiration (ET) is one of the most important factors controlling hydrologic, agricultural, and weather cycles. It also converts a large portion of rainfall into vapor, being known as the largest water flux from the earth into the atmosphere. Since ET is affected by many factors, such as land surface characteristics and climatic conditions, it undergoes considerable spatiotemporal variations, particularly at the watershed scale. Hence, to obtain a more accurate estimation of ET, it is required to identify homogenous and uniform regions, each represented by a meteorological station. In this study, three scenarios were proposed in order to identify homogenous regions to estimate ET based on METRIC method, and the scenarios were tested in Sefidrood Watershed in the north of Iran. The first scenario included only vegetation factor with one representative station for the entire case study watershed and ignored diverse conditions affecting ET across the watershed. The second scenario incorporated not only the vegetation factor but also the altitudinal variations of the watershed. In the second scenario, the watershed was divided into two distinct altitudinal sections, each with a representative station with a specific influenced area, with ET being estimated separately for each section. Finally, the third scenario incorporated the altitudinal and climatic variations. The results indicated that the second scenario performed better than two other scenarios in ET estimation. In other words, altitude and vegetation strongly influenced spatial and temporal distributions of ET, leading to considerable variations of it in the watershed.

Genome-wide association study and pathway analysis identify NTRK2 as a novel candidate gene for litter size in sheep

Litter size is one of the most important economic traits in sheep. Identification of gene variants that are associated with the prolificacy rate is an important step in breeding program success and profitability of the farm. So, to identify genetic mechanisms underlying the variation in litter size in Iranian Baluchi sheep, a two-step genome-wide association study (GWAS) was performed. GWAS was conducted using genotype data from 91 Baluchi sheep. Estimated breeding values (EBVs) for litter size calculated for 3848 ewes and then used as the response variable. Besides, a pathway analysis using GO and KEGG databases were applied as a complementary approach. A total of three single nucleotide polymorphisms (SNPs) associated with litter size were identified, one each on OAR2, OAR10, and OAR25. The SNP on OAR2 is located within a novel putative candidate gene, Neurotrophic receptor tyrosine kinase 2. This gene product works as a receptor which is essential for follicular assembly, early follicular growth, and oocyte survival. The SNP on OAR25 is located within RAB4A which is involved in blood vessel formation and proliferation through angiogenesis. The SNP on OAR10 was not associated with any gene in the 1Mb span. Moreover, gene-set analysis using the KEGG database identified several pathways, such as Ovarian steroidogenesis, Steroid hormone biosynthesis, Calcium signaling pathway, and Chemokine signaling. Also, pathway analysis using the GO database revealed several functional terms, such as cellular carbohydrate metabolic, biological adhesion, cell adhesion, cell junction, and cell-cell adherens junction, among others. This is the first study that reports the NTRK2 gene affecting litter size in sheep and our study of this gene functions showed that this gene could be a good candidate for further analysis.

A genome-wide association study revealed five SNPs affecting 8-month weight in sheep

Lamb weight at 8 months of age is an important trait in the sheep industry in terms of the onset of puberty around this age; however, knowledge of its effective genetic factors is limited. Therefore, a GWAS using the 50K SNP-Chip was performed on 96 Baluchi sheep to identify the genomic regions associated with 8-month weight. The results of the present study revealed five SNPs on chromosomes 4, 14 and 16 at 5% chromosome-wide significance level, jointly accounting for 0.95% of total genetic variance. Four genes - MTPN, HYDIN, LRGUK and ZFP90 - were found in 50 kb intervals around the significant SNPs, of which MTPN is involved in regulation of skeletal muscle growth. Our results may provide a new vision to identify the genomic regions affecting growth traits in sheep.

Quantifying the effect of inbreeding on average daily gain and Kleiber ratio in Mazandaran native chickens

The aims of this study was to evaluate inbreeding effects for growth rate (average daily gain from hatch to 8 weeks of age (ADG1-8), average daily gain from hatch to 12 weeks of age (ADG1-12), average daily gain from 8 weeks of age to 12 weeks (ADG8-12), average daily gain from hatch to sex maturity (ADG1-SM), average daily gain from 8 weeks of age to sex maturity (ADG8-SM), average daily gain from 12 weeks to sex maturity (ADG12-SM), and Kleiber ratios (KR1-8, KR1-12, KR8-12) in Mazandaran native chickens. The data set of 82,446 pedigreed individuals over 21 generations was analyzed. All individuals were grouped into three classes based on the inbreeding coefficient: First class contained non-inbred birds (F = 0) and second and third classes contained inbred individuals (0 < F ≤ 5% and F > 5%, respectively). The effects of inbreeding coefficient on the studied traits were estimated by the restricted maximum likelihood method (REML) applying the WOMBAT software and fitting individual increase in inbreeding coefficient (ΔFi) of birds as linear covariates under 6 different animal models. ADG8-SM and ADG12-SM increased non-significantly by 0.003 and 0.005 g, respectively, due to 1% increase in ΔFi, while ADG1-8, ADG1-12, ADG8-12, ADG1-SM, KR1-8, and KR1-12 decreased (P < 0.001), respectively, by 0.012 g, 0.011 g, 0.010 g, 0.014 g, 0.003 g, and 0.002 g. Also, KR8-12 non-significantly decreased by 0.001 g. Estimates of direct heritability of traits based on the most appropriate models ranged from 0.05 (KR8-12) to 0.26 (ADG12-SM). Results showed that it is important to regularly analyze inbreeding level in the herd in an attempt to prevent the reduction of performance.

Comparison of parametric, semiparametric and nonparametric methods in genomic evaluation

Access to dense panels of molecular markers has facilitated genomic selection in animal breeding. The purpose of this study was to compare the nonparametric (random forest and support vector machine), semiparametric reproducing kernel Hilbert spaces (RKHS), and parametric methods (ridge regression and Bayes A) in prediction of genomic breeding values for traits with different genetic architecture. The predictive performance of different methods was compared in different combinations of distribution of QTL effects (normal and uniform), two levels of QTL numbers (50 and 200), three levels of heritability (0.1, 0.3 and 0.5), and two levels of training set individuals (1000 and 2000). To do this, a genome containing four chromosomes each 100-cM long was simulated on which 500, 1000 and 2000 evenly spaced single-nucleotide markers were distributed. With an increase in heritability and the number of markers, all the methods showed an increase in prediction accuracy (P<0.05). By increasing the number of QTLs from 50 to 200, we found a significant decrease in the prediction accuracy of breeding value in all methods (P<0.05). Also, with the increase in the number of training set individuals, the prediction accuracy increased significantly in all statistical methods (P<0.05). In all the various simulation scenarios, parametric methods showed higher prediction accuracy than semiparametric and nonparametric methods. This superior mean value of prediction accuracy for parametric methods was not statistically significant compared to the semiparametric method, but it was statistically significant compared to the nonparametric method. Bayes A had the highest accuracy of prediction among all the tested methods and, is therefore, recommended for genomic evaluation.

Detecting selection signatures in three Iranian sheep breeds

The objective of genome mapping is to achieve valuable insight into the connection between gene variants (genotype) and observed traits (phenotype). Part of that objective is to understand the selective forces that have operated on a population. Finding links between genotype-phenotype changes makes it possible to identify selective sweeps by patterns of genetic variation and linkage disequilibrium. Based on Illumina 50KSNP chip data, two approaches, XP-EHH (cross-population extend haplotype homozygosity) and F_ST (fixation index), were carried out in this research to identify selective sweeps in the genome of three Iranian local sheep breeds: Baluchi (n = 86), Lori-Bakhtiari (n = 45) and Zel (n = 45). Using both methods, 93 candidate genomic regions were identified as harboring putative selective sweeps. Bioinformatics analysis of the genomic regions showed that signatures of selection related to multiple candidate genes, such as HOXB9, HOXB13, ACAN, NPR2, TRIL, AOX1, CSF2, GHR, TNS2, SPAG8, HINT2, ALS2, AAAS, RARG, SYCP2, CAV1, PPP1R3D, PLA2G7, TTLL7 and C20orf10, that play a role in skeletal system and tail, sugar and energy metabolisms, growth, reproduction, immune and nervous system traits. Our findings indicated diverse genomic selection during the domestication of Iranian sheep breeds.

A genome scan for quantitative trait loci affecting average daily gain and Kleiber ratio in Baluchi Sheep

Genomewide association study (GWAS) is an efficient tool for the detection of SNPs and candidate genes in quantitative traits. Growth rate is an important trait for increasing the meat production in sheep. A total of 96 Baluchi sheep were genotyped using Illumina Ovine SNP50 BeadChip to run a GWAS for an average daily gain (ADG) and Kleiber ratio (KR) traits in different periods of age in sheep. Traits included were average daily gain from birth to three months (ADG_0-3), from three months to six months (ADG_3-6), from six months to nine months (ADG_6-9), from nine months to yearling (ADG_9-12), from birth to six months (ADG_0-6), from three months to nine months (ADG_3-9), from three months to yearling (ADG_3-12) and corresponding Kleiber ratios (KR_0-3, KR_3-6, KR_6-9, KR_9-12, KR_0-6, KR_3-9 and KR_3-12, respectively). A total of 42,243 SNPs passed the quality-control filters and were analysed by PLINK software in a linear mixed model. Two SNPs were identified on two chromosomes at the 5% genomewide significance level for KR(3-9) and KR(6-9). Two candidate genes, namely MAGI1 and ZNF770, were identified correspondingly harbouring and close to these QTL. Also, a total of 21 SNPs were found on chromosomes 2, 3, 5, 6, 7, 10, 17, 19, 20 and 25 at the 5% chromosomewide significance level for ADG and KR traits. Thus, we suggest more studies to discover the causative variants for growth traits in sheep.

Detection of QTL for greasy fleece weight in sheep using a 50 K single nucleotide polymorphism chip

Genome-wide association studies (GWAS) have introduced an influential tool in the search for quantitative trait loci (QTL) influencing economically important traits in sheep. To identify QTL associated with greasy fleece weight, a GWAS with 50 K single nucleotide polymorphisms (SNPs) was performed in a Baluchi sheep population. Association with greasy fleece weights was tested using the software Plink. The results of our GWAS provided three novel SNP markers and candidate genes associated with greasy fleece weight. A total of three chromosome-wide significant associations were detected for SNP on chromosomes 17 and 20 affecting greasy fleece weight across the four shearing. One of the significant SNP markers was located within ovine known genes namely FAM101A. Further investigation of these identified regions in validation studies will facilitate the identification of strong candidate genes for wool production in sheep.

Population structure of Mazandaran native fowls using pedigree analysis

The objective of this study was to use pedigree analysis to evaluate the population structure and genetic variability of the Mazandaran native fowls in Iran by quantifying the pedigree completeness index, effective population size, genetic diversity, inbreeding level, and individual increase in inbreeding. The pedigree completeness analysis showed 3.31 full, 10.19 maximum, and 6.30 equivalent generations. The effective number of founders (f _e) was 131, representing 5% of the potential number of founders. The effective number of ancestors (f _a) was 81, and the genetic contribution of the 37 most influent ancestors explained 50% of the genetic variability in the population. The ratio f _e/f _a (effective number of founders/effective number of ancestors), which expresses the effect of population bottlenecks, was 1.62. The inbreeding coefficient increased over generations and the average was 1.93%. The average relatedness coefficient between individuals of the population was estimated to be 2.59%. The effective population size, based on the number of full generations, was 56. Family size analysis showed that fewer males than females were used, resulting in the observed levels of inbreeding. Average inbreeding coefficient in the Mazandaran native fowls can be regarded to be below critical levels. However, considering the relationship coefficients of individuals is recommended to aid maintaining genetic diversity of Mazandaran native fowls.

A mechanistic insight into the effect of piperidine as an organocatalyst on the [3 + 2] cycloaddition reaction of benzalacetone with phenyl azide from a computational study

Several transition structures (TSs) for catalyst-free [3 + 2] cycloaddition and two plausible mechanistic pathways for the organocatalyzed [3 + 2] cycloaddition (32CA) between benzalacetone and phenyl azide were located by quantum chemistry methods. Calculations were carried out with B3LYP, MPWB1K and M06-2X functionals using 6-31G(d) and 6-311G(d,p) basis sets in gas and solvent phases. The calculated activation barriers imply that the lowest barrier pathway is the catalyzed process producing 3-regioisomers through the iminium intermediate and not through the dienamine route. Electronic displacements along the reaction path have been examined using a topological analysis of the electron-localization function (ELF). ELF topological analyses along the intrinsic reaction coordinates (IRC) of both catalyzed and uncatalyzed 32CA reactions indicated that while the first C1-N1 single bond is formed as a dative bond, the formation of the second C2-N3 bond takes place via a C-to-N coupling between the interacting centers of the reagents. Moreover, the ELF analyses imply that the reaction mechanism is a two-stage one-step process in the presence of a piperidine organocatalyst, while bond formation in an uncatalyzed process is almost synchronous.

The origin of regio- and stereoselectivity in the 1,3-dipolar cycloaddition of nitrile oxides with C1-substituted 7-oxabenzonorbornadienes, a DFT study

The origin ofexo/antiselectivity in the 1,3-dipolar cycloaddition of nitrile oxides with C₁-substituted 7-oxabenzonorbornadienes has been investigated theoretically.

Measurement of within and between genetic variability in duck breeds by RAPD markers

A total of 200 genomic DNAs were isolated from the four breeds of duck: Muscovy, Pekin, Khaki Campbell and Native, through a modified salting out procedure. The samples were used in a Polymerase Chain Reaction (PCR) with 27 RAPD markers. Amplified PCR-products with the markers were separated on a 2% agarose gel and stained with ethidium bromide. To evaluate the bands, polymorphic and monomorphic bands were described. Genetic similarity and genetic distance were calculated. The RAPD analysis data from 7 primers were utilized in estimating genetic similarity between four breeds, which ranged from 0.46 to 0.77 between the breeds. The maximum genetic distance was observed between Muscovy and Native duck breeds (0.76).