Mahdi JavanNikkhah; Hossein Moradi Shahrbabak; Mohammad Hossein Moradi; Ali Sadeghi-Sefidmazgi
Abstract
Bovine viral leukemia (BLV) is one of the deadliest viral diseases that is associated with many economic losses in the dairy industry, such as reduced production capacity and reproductive performance of infected animals and their eventual culling. The aim of this study was to identify the selection signatures ...
Read More
Bovine viral leukemia (BLV) is one of the deadliest viral diseases that is associated with many economic losses in the dairy industry, such as reduced production capacity and reproductive performance of infected animals and their eventual culling. The aim of this study was to identify the selection signatures associated with resistance to BLV in Iranian Holstein cows. For this purpose, a total of 152 animals were genotyped for 30,105 SNP markers using GGP Bovine LD v4 chips. After quality control of the initial data, 23,513 SNP markers in 140 animals of cattle were finally entered for further analysis. The animals used were classified into two groups consisting resistant to disease or healthy (68 animals) and sick (77 animals) animals, and then the regions of the genome that were divergently selected in these animals were evaluated using the unbiased Theta method. The results of this study showed that four genomic regions on chromosomes 1, 13, 20 and 22 were divergently selected in these groups. Study of the genes reported in these regions revealed that some genes such as the STGA1 on chromosome 1, the STK35, EBF4, and PDYN on chromosome 13 and the SLC38A3, RASSF1, and RBM6 on chromosome 22 were previously reported within or adjacent to these genomic regions. Study the function of these genes showed that the genes are involved in the immune system, the regulation of mitotic and meiotic cycles and cancer suppression. Overall, the results of this study can provide a valuable source of information to identify candidate genomic regions or causal genes associated with this disease.
Mahdi Mokhber; Mohammad Moradi Shahre Babak; Mostafa Sadeghi; Hossein Moradi Shahrbabak; Javad Rahmani-Nia
Abstract
In order to estimate the effective population size (Ne) in Iranian water buffalo blood and hair samples of 407 individual from Azari (N=260), Khuzestani (N=120) and Mazandrani (N=27) buffalo populations were gathered. After DNA extaraction, the samples were genotyped using Axiom® Buffalo Genotyping ...
Read More
In order to estimate the effective population size (Ne) in Iranian water buffalo blood and hair samples of 407 individual from Azari (N=260), Khuzestani (N=120) and Mazandrani (N=27) buffalo populations were gathered. After DNA extaraction, the samples were genotyped using Axiom® Buffalo Genotyping 90K Array. The Ne was estimated from 700 to 4 generations ago and also for the present generation by linkage disequiblirum data and based on heterozygote-excess method using NeEstimator (V2), respectively. Estimated Ne for Azari, Khuzestani and Mazandarani were calculated 1530, 1375 and 1141, respectively, for 700 generations ago. Ne for the present generation in Azeri, Khuzestani and Mazandarani were estimated 447, 226 and 35, respectively. The Ne for Azeri and Khuzestani were relatively high and these two populations were not endanger to extinction, but their Ne has been declined in the resent generations massively and it is necessary to care about the maintenance of Ne and relatively high diversity for these populations. However, the Mazandarani population is endangered because of low Ne and so it is necessary to carefully monitor their effective population size, improve the profitability of production and planning a suitable mating scheme to control inbreeding and genetically conserve this population.
Mahdi Mokhber; Mohammad Moradi Shahre Babak; Mostafa Sadeghi; Hossein Moradi Shahrbabak; Javad Rahmani-Nia
Abstract
In order to detect signature of selection on buffalo genome, a set of 287 water buffalo samples from 260 Azari and 27 Mazandarani buffalo breeds were genotyped using the Axiom® Buffalo Genotyping 90K Array. The unbiased fixation index method (FST) was used to detect signatures of selection. In total, ...
Read More
In order to detect signature of selection on buffalo genome, a set of 287 water buffalo samples from 260 Azari and 27 Mazandarani buffalo breeds were genotyped using the Axiom® Buffalo Genotyping 90K Array. The unbiased fixation index method (FST) was used to detect signatures of selection. In total, 14 regions with outlier FST values (0.1%) were identified. Annotation of these regions using the UMD3.1 Bos taurus Genome Assembly was performed to find putative candidate genes and QTLs within the selected and 105 genes and 28 QTLs with selection signatures were detected. A high proportion of identified genes (N=27) in regions under selection were involved in olfactory receptor, also some of the detected genes were associated with growth and body development, metabolicand apoptosis possesses, immune system development, and mammary gland development. Some of the identified QTLs in regions under selection were associated with growth traits such as body weight at birth, weaning and mature, subcutaneous fat, meat yield and carcass weight. The detected QTL for milk traits were only associated with milk contents and somatic cell count. However, it is recommended to carry out association studies to show the actual function of these genes.
Mohamadreza Attari; Hossein Moradi Shahrbabak; Gholamali Nehzati Paghale; Mohamad Hossein Banabazi; Mohamadreza Hashemi
Abstract
The aim of this study was to investigate the gene expression profile and find the important genes in the differentiation and evolution of the queen, worker and drone honey bee at Stage 5 Larvae. Therefore, transcripts (total mRNA sequence) of 15 samples of Italian honey bee (A. m. ligustica) including ...
Read More
The aim of this study was to investigate the gene expression profile and find the important genes in the differentiation and evolution of the queen, worker and drone honey bee at Stage 5 Larvae. Therefore, transcripts (total mRNA sequence) of 15 samples of Italian honey bee (A. m. ligustica) including 5 drone, 5 worker and 5 queen bees were aligned to the reference genome of the honey bee. In gene expression analysis of RNA-Seq data, 15962 genes and 31297 isoforms were identified. In our differential gene expression analyses, 465 genes between drone and queen bees, 495 genes between worker and queen bees and 764 genes between drone and worker bees, were expressed differently (P <0.000005). The largest difference in expression of genes was observed between drone and workers were for GB45614 and GB42053, with log2 fold change that was -10 and 11.5, respectively. In drone and queen bees comparison, GB45614 and GB48020 genes with log2 fold change 7.11 and 8.11, and in queen and worker bees comparison, GB43508 and GB42053 with log2 fold change 6.6 and -9, had the largest difference in gene expression, respectively. The analysis of the gene ontology (GO) and the pathways involved showed that the function of many of these genes has yet to be found. However, a large number of the genes expressed defiantly in drone and queen bees were related to integral component of membrane, calcium ion binding, carboxypeptidase, cholecystokinin receptor,chitin metabolic process, chymotrypsin inhibitor, haemolymph juvenile hormone binding and pupal cuticle protein, while differentially expressed genes in queen and worker bees comparison were related to metabolic pathways, enzymes metabolism, Pyrimidine metabolism, lipid metabolism, protein kinase, ATP binding and nucleic acid, intracellular cholesterol transport, chitin metabolic process, nitrogen compound metabolic process, hydrolase activity and chitin binding. The genes expressed at different levels in worker bees and drones were related to structural elements, the metabolism membrane and transfer of amino acids, calcium-ion-bound, ion-binder, intracellular cholesterol transport and chitin metabolism.
Ali Jalil Sarghale; Mohammad Moradi Shahre Babak; Hossein Moradi Shahrbabak; Ardeshir Nejati Javaremi; Mahdi Saatchi
Abstract
The methane production from ruminant production system was estimated to reach 250-500 L per animal per day which has been reported to contribute up to 8-10 % of global warming during the next 50-100 years. The aim of this study was to investigate the correlation among methane emission (predicted by volatile ...
Read More
The methane production from ruminant production system was estimated to reach 250-500 L per animal per day which has been reported to contribute up to 8-10 % of global warming during the next 50-100 years. The aim of this study was to investigate the correlation among methane emission (predicted by volatile fatty acids) with milk production traits, its components and breeding values (BV) of these traits in Iranian Holstein cattle. The rumen digesta was obtained from 150 cattle through stomach tubing and this population divided into 2 groups with 75 cattle in each (the groups have different milk production BV). Data were analyzed by R.3.3.0. The results showed that methane emission per unit of milk and fat were different in the two groups (P<0.0001). Also, the BVs of milk production, fat and protein traits and daily production of milk, fat and protein had weak to moderate negative correlation with methane emission per unit(P<0.05). The highest correlation was observed between daily production of fat with methane emission per unit of fat (-0.79) as well as daily milk production with methane emission per unit of milk (-0.62). These results showed that methane emission may be reduced by indirect selection per generation for the traits had a high correlation with the gas (daily production of milk and fat).
Mohammad Hossein Fallahi; Hossein Moradi Shahrbabak; Mohammad Moradi Shahrbabak; Rostam Abdollahi Arpanahi
Abstract
The aim of this study was to determine genomewide linkage disequilibrium (LD), Haplotype block and effective population size using the information obtained from 243 Azarbaijani breed buffalo using a high density SNP panel (Axiom® Buffalo Genotyping 90K). After quality control of SNP markers data, ...
Read More
The aim of this study was to determine genomewide linkage disequilibrium (LD), Haplotype block and effective population size using the information obtained from 243 Azarbaijani breed buffalo using a high density SNP panel (Axiom® Buffalo Genotyping 90K). After quality control of SNP markers data, 62,141 SNP markers remained for identification of linkage disequiliberum, haplotype blocks and effective population size. LD was measured by the square of the correlation coefficient (r2) between alleles. The maximum LD measured by r2 varied from 0.25 to 0.29 at a distance of < 2.5 kb, and the minimum average values of r2 varied from 0.012 to 0.014 at distances ranging from 900 to 1000 kb, clearly showing that the average r2 reduced with the increase in SNP pair distances. Overall, 1693 blocks were observed through the genome. Eleven percent of all SNPs were clustered into haplotype blocks, covering 202 ±3.4 Mb of the total autosomal genome size. Effective population size (Ne) was estimated based on expected linkage disequilibrium. Ne was estimated to be 422 in our population.
Hamidreza Abdollahi; Hassan Mehrabani Yeganheh; Hossein Moradi Shahrbabak
Abstract
The aim of this research was to utilize a simple, rapid, cheap and safe method for sex determination of newly born canaries through feather to help accurately breeding this bird. The external morphology of males and females of many birds such as canary are fairly similar to each other, Therefore, determining ...
Read More
The aim of this research was to utilize a simple, rapid, cheap and safe method for sex determination of newly born canaries through feather to help accurately breeding this bird. The external morphology of males and females of many birds such as canary are fairly similar to each other, Therefore, determining their sexes is very difficult. The different length of CHD (chromo-helicase-DNA-binding) gene in Z and W chromosomes can recignized by using molecular technique based and PCR method, for the sex determination of young birds. In this study, blood and feathers samples of 140 canaries from 10 different breeds were used and their DNA were extracted. It should be mentioned that for PCR, at first 5 primers designed and finally 1 specific primer were confirmed and used. The results showed that the amplification of a fragment with the length of 252 bp for males (ZZ) and two fragments with the lengths of 252 and 309 bp for female (ZW). The bands had a length difference of 57 bp. The DNA sequencing of these two fragments was done for the first time, the sequence of these two fragments of the CHD gene in canary species was submitted and accepted in the NCBI site under the names of CHDZ-1 and CHDW-1 with the accession numbers of MG679825.1 and MG679826.1.
Hossein Mohammadi; Seyed Abbas Rafat; Hossein Moradi Shahrbabak; Jalil Shoja; mohammad hosein moradi
Abstract
Identification of selection targeted genomic regions is one of the main aims of biological research.The objective of this study was a genome-wide scan to identify the genomic regions that have been under artificial and natural selection in Zandi sheep breed. For this purpose, 96 animal of Zandi breed ...
Read More
Identification of selection targeted genomic regions is one of the main aims of biological research.The objective of this study was a genome-wide scan to identify the genomic regions that have been under artificial and natural selection in Zandi sheep breed. For this purpose, 96 animal of Zandi breed have been genotyped using the Illumina ovine SNP50 BeadChip. The intergrated hapl-type score (iHS) test was used to detect the selection sweep, due to linkage disequilibrium, associated with these signatures. The results revealed eleven genomic regions on 1 (two areas), 3, 6, 7,8, 9, 10, 17, 22 and 26 chromosomes. Bioinformatics analysis demonstrated that some of these genomic regions overlapped with reported genes that directly and indirectly influenced traits for adaptation to hot arid environments (DNAJB4, HSPA4L, MSRB3), immune response (IL23A, STAT6, LY96), body size and development (STAC3, LAP3), development of the skeletal system (SPP1, MEPE, IBSP) and energy and digestive metabolism (ATP5B, GLS2, CS). Finally, study of the reported QTL in these regions of the sheep genome showed that they overlapped with QTL of economically important traits such as carcass yield, growth and wool traits. However, it will be necessary to carry out more association and functional studies to demonstrate the implication of these genes.
Zahra Azizi; Hossein Moradi Shahrbabak; Mohammad Moradi Shahrbabak
Abstract
Understanding of population genetic structure is valuable for better implementation of breeding programs and most importantly, preservation of genetic resources. Genomic data provide an opportunity to consider complex evolutionary history of populations and reconstruct rare historical events. In this ...
Read More
Understanding of population genetic structure is valuable for better implementation of breeding programs and most importantly, preservation of genetic resources. Genomic data provide an opportunity to consider complex evolutionary history of populations and reconstruct rare historical events. In this research, the structure of Iranian buffalo populations was studied by using principal component analysis and discriminant analysis principal component methods. For this purpose, the number of 404 buffalos from three breeds including North, Azari and Khozestani were sampled and genotyped by SNPChip 90k from Padano Company in Italy. The results of principal component analysis and discriminant analysis principal component showed a clear picture of the genetic structure of the studied populations. Assessing the optimal number of clusters with criteria BIC, K = 3 by the DAPC method showed the best results. The result of cross-validation for retaining principal components was optimized to 50 first components that showed the lowest MSE. In this study, DAPC predicted assignment of individuals to clusters and membership probabilities with 100% accuracy. PCA method was not able to provide a group assessment and DAPC method outperformed than PCA in achieving a clear variance difference between populations. DAPC method can be applied in quality control and stratification population correction of GWAS as an alternative to the PCA because of summarizing the genetic differentiation between groups and overlooking within-group variation and providing better population structure.
Zahra Azizi; Hossein Moradi Shahrbabak; Mohammad Moradi Shahrbabak; Abbas Rafat; Jalil Shodja
Abstract
The purpose of this research was to classify buffaloes from different areas of the two Azari (West and East Azarbayjan and Ardabil provinces) and North (Guilan province) ecotypes using support vector machine method. A total of 258 buffalo were sampled and genotyped using the Axiom Buffalo 90K Genotyping ...
Read More
The purpose of this research was to classify buffaloes from different areas of the two Azari (West and East Azarbayjan and Ardabil provinces) and North (Guilan province) ecotypes using support vector machine method. A total of 258 buffalo were sampled and genotyped using the Axiom Buffalo 90K Genotyping Array at the Parco Technologic Padano lab in Italy. Two metric methods of cross validation and the area under the receiver operating characteristic (AUC) were used to determine the predictive performance of support vector machine (SVM) to classify individuals. The results of cross validation and methods for classifying different regions of the two ecotypes (4 provinces) were 92% and 96%, respectively that showed despite the difficulty of identifying individuals from provinces close to each other, support vector machine (SVM) method shows higher accuracy in assigning animals to their herds. Result of two ecotypes showed accuracy about 96% and 98% which represents the better ability to separate the two ecotypes. Machine learning method provides predictions for classification of each individual which can be efficient in quality control and genetic studies.
Homa Aarabi; Mohammad Moradi Shahrbabak; Abbas Pakdel; Hossein Moradi Shahrbabak; Ali Esmailizadeh koshkoiyeh
Abstract
Insulin-like growth factor-1 gene (IGF1) is a biological candidate gene for investigation traits such as growth, body weight and muscles growth in different species. In this research eight pairs of white and wild quails were crossed reciprocally as a base population. A total of 34 quails were produced ...
Read More
Insulin-like growth factor-1 gene (IGF1) is a biological candidate gene for investigation traits such as growth, body weight and muscles growth in different species. In this research eight pairs of white and wild quails were crossed reciprocally as a base population. A total of 34 quails were produced in first generation (F1) and 422 quails were generated by intercrossing the F1 population. Body weights at the time of hatching and different weeks were recorded in the second generation. The association between IGF1 promoter region polymorphism and body weight in different ages was investigated in F2 generation. Using PCR-SSCP assay and sequencing, a novel single nucleotide polymorphism (SNP) was identified in IGF1 promoter region in 472 birds from three Generations. Genotypic frequencies of AA, AG and GG genotypes were 0.05, 0.62 and 0.33, respectively for all generations. The frequency of A and G alleles were 0.36 and 0.64, respectively. The AA genotype was lower than AG and GG genotypes for body weights in ages of one to five weeks and slaughter time and the average of daily gain was also lower from week first to four but these differences weren’t significant. The average of daily gain was higher in females compared to males for AG genotype (P<0.001). Comparison between detected allele in the present study with reported allele by other research groups in chickens revealed a novel SNP in promoter of IGF1 gene in Japanese quail, but association between this polymorphism with body weight and growth rate were not significant.
Hossein Moradi Shahrbabak; Rasoul Khodabakhshzadeh; Mostafa Sadeghi; Amir TaheriYeganeh
Abstract
The growth hormone (GH) gene is a candidate for growth trait in farm animals and plays an important role in growth metabolism. Growth hormone can contribute in growth metabolism when it\'s receptor (GHR) are on cells target. The single nucleotide polymorphism was occurred in GH and GHR genes and that ...
Read More
The growth hormone (GH) gene is a candidate for growth trait in farm animals and plays an important role in growth metabolism. Growth hormone can contribute in growth metabolism when it\'s receptor (GHR) are on cells target. The single nucleotide polymorphism was occurred in GH and GHR genes and that associated with carcass traits in Lori-Bakhtiari sheep breed with using PCR-SSCP. In this research, blood samples were collected from the left jugular vein from 152 Lori-Bakhtiari sheep breed. Genomic DNA was extracted from blood samples by using the salting-out procedure and PCR were used to amplify the regions are located in exon 4 (214 bp segments) of the ovine GH gene and exon 10 (218 bp segments) of the ovine GHR Gene. The single stranded conformation polymorphism (SSCP) patterns of PCR products were studied using acrylamide gel electrophoresis and silver-nitrate staining method. The result showed polymorphism in GH gene but any difference between banding patterns for GHR gene was not observed. The sequencing results showed the presence of 5 Single nucleotide polymorphism for GH gene in the studied population. No significant associations of the available genotypes in the exon 4 of the ovine GH Gene with carcass traits.
Javad Rahmaninia; Seyed Reza Miraei-Ashtiani; Hossein Moradi Shahrbabak
Abstract
High through put sequencing of single nucleotide polymorphisms (SNP) has revolutionized the fine scale analysis of the population structure in different species. Various methods have been proposed and used for the study of population structure using whole-genome marker data that each has advantages and ...
Read More
High through put sequencing of single nucleotide polymorphisms (SNP) has revolutionized the fine scale analysis of the population structure in different species. Various methods have been proposed and used for the study of population structure using whole-genome marker data that each has advantages and disadvantages with respect to their characteristics. Super Paramagnetic Clustering (SPC) which is based on data mining was used in this study in order to investigate the population and sub-population structures in simulated populations. The purpose of applying this method was to achieve population structure without using any information from ancestral population. After editing the data, 29209 autosomal markers from 159 animals were analyzed. The results showed that animals are placed properly in their respective population and sub-populations based on their similarities and dissimilarities. The main advantages of this method are the computational efficiency and not requiring any prior assumptions. Therefore, it might be used to analyze the data from thousands of animals without any pedigree and ancestry information to reveal their population structure.
Mahdi Mokhber; Mohammad Moradi Shahrbabak; Mostafa Sadeghi; Hossein Moradi Shahrbabak; John Williams
Abstract
Identification of selection targeted genomic regions is one of the most challenging areas of research in animal genetics, particularly in livestock. We carried out a genome-wide scan for signatures of positive selection to identify genomic regions that had been under selection in Iranian Khuzestani and ...
Read More
Identification of selection targeted genomic regions is one of the most challenging areas of research in animal genetics, particularly in livestock. We carried out a genome-wide scan for signatures of positive selection to identify genomic regions that had been under selection in Iranian Khuzestani and Mazandrani buffalo breeds. A total of 148 water buffalo from Khuzestani (N=121) and Mazandrani (N=27) buffalo breeds were genotyped using Axiom® Buffalo Genotyping 90K Array. Unbiased method of population differentiation index (FST) was applied to detect signatures of selection. In total, 23 regions exceeding the 0.1 percent threshold of the empirical posterior distribution were identified as extremely differentiated. These selected genomic regions were surveyed to find encoding putative candidate genes and 64 genes and 27 QTL were extracted from the corresponding areas in UMD3.1 Bos Taurus Genome Assembly. Some of these genes have previously reported as signature of positive selection in the last studies. Some of these genes were also found to be involved in milk production traits and domestication-related changes include sensory perceptions, brain and neural system development, pigmentation, and geographic adaptation. Also, survey on extracted QTLs was shown that these QTLs involved in some economicl important traits in buffalo such as feed conversion ratio, subcutaneous fat, body weight, average daily gain, type, Meat tenderness, milk production conentent, udder attachment, calf size and calving ease traits. However, it will be necessary to carry out more association and functional studies to demonstrate the implication of these genes.
Hamed Rezaei; Hosein Moradi Shahrbabak; Mohammad Moradi Shahrbabak; Mahdi Bayeriyar
Abstract
The data used in this research were according to survival rate records of 9700 lambs from 735 rams and 8082 ewes, which collected at Moghani Flock Breeding Research Station during the year 1365 to 1390. The frequency distribution causes of culling and survival distribution function of lambs to yearling ...
Read More
The data used in this research were according to survival rate records of 9700 lambs from 735 rams and 8082 ewes, which collected at Moghani Flock Breeding Research Station during the year 1365 to 1390. The frequency distribution causes of culling and survival distribution function of lambs to yearling age were estimated by SAS (2000). Results revealed that 49.81 percent of total lambs were culled from the flock until the first year of life because of reasons including illness, lack of breed purity, etc. The mortality rate over this period was estimated 20.81 percent which occurred 8.03, 7.5, 4.43 and 0.85 percent in first three months, second three months, third three months and fourth three months of lambs life, respectively. According to percentage of mortality rate, the accumulative survival of lambs from birth to yearling age was calculated 79.19 percent. The regression coefficient of lambs survival from birth up to yearling age showed a monthly 1.77 percent decrease in lambs survival of this breed. The most mortality rate occurred in first, fifth, second, fourth, sixth and ninth months of lambs life, respectively. In comparison with third and fourth three months of lambs life, the fatality in first and second three months were higher because of heavy rainfall and frigid weather in winter quarters.
Mahdi Mokhber; Hosein Moradi Shahr Babak; Amirhosein Khelt Abadi Farahani
Abstract
The objective followed in the present study was to survey the relationship between 18 body trait measurements (live weight, height at wither, paunch girth, neck diameter, body length, girth around the body, width of fat tail at above, below and midpoint of fat tail, fat tail length lowers right and left ...
Read More
The objective followed in the present study was to survey the relationship between 18 body trait measurements (live weight, height at wither, paunch girth, neck diameter, body length, girth around the body, width of fat tail at above, below and midpoint of fat tail, fat tail length lowers right and left sides, fat tail gap length, fat tail depth at the above, below, and midpoint, and girth around fat tail at the above, mid and down point ) and the traits of fat tail weight and carcass (weight with and without fat tail) in Makooei Sheep and to Predict these traits, performance though multivariate Linear Regression Method Based on Principal Component Analysis. Sex showed significant effects on all the measured traits expect fat tail, depth at above, midpoint and down point, fat tail length towards right side plus gap length, width of fat tail at the above and below point of tail. Means and standard error for each trait were evaluated by sex. Multicollinearity was detected through a survey of the relationship among these traits, variance inflation factor and tolerance value. Principal component analysis was employed to resolve multicollinearity problem among independent variables and for a clearer explanation of the results. R2 range for different regression models varied between 0.973 for carcass weight without tail in male sex and 0.561 for carcass weight with fat tail in female sex of Makooei lambs.
