Study of Population Structure and Genetic Prediction of Buffalo from Different Provinces of Iran using Machine Learning Method

Considering breeding livestock programs to milk production and type traits based on existence two different ecotypes of Iranian’s buffalo, a study carried out to investigate the population structure of Iranian buffalo and validate its classification accuracy according to different ecotypes from Iran (Azerbaijan and North) using data SNP chip 90K by means Support vector Machine (SVM), Random Forest (RF) and Discriminant Analysis Principal Component (DAPC) methods. A total of 258 buffalo were sampled and genotyped. The results of admixture, multidimensional scaling (MDS), and DAPC showed a close relationship between the animals of different provinces. Two ecotypes indicated higher accuracy of 96% that the Area Under Curve (AUC) confirmed the obtained result of the SVM approach while the DAPC and RF approach demonstrated lower accuracy of 88% and 80 %, respectively. SVM method proved high accuracy compared with DAPC and RF methods and assigned animals to their herds with more accuracy. According to these results, buffaloes distributed in two different ecotypes are one breed, and therefore the same breeding program should be used in the future. The water buffalo ecotype of the northern provinces of Iran and Azerbaijan seem to belong to the same population

Genome Diversity and the Origin of the Arabian Horse

The Arabian horse, one of the world\u27s oldest breeds of any domesticated animal, is characterized by natural beauty, graceful movement, athletic endurance, and, as a result of its development in the arid Middle East, the ability to thrive in a hot, dry environment. Here we studied 378 Arabian horses from 12 countries using equine single nucleotide polymorphism (SNP) arrays and whole-genome re-sequencing to examine hypotheses about genomic diversity, population structure, and the relationship of the Arabian to other horse breeds. We identified a high degree of genetic variation and complex ancestry in Arabian horses from the Middle East region. Also, contrary to popular belief, we could detect no significant genomic contribution of the Arabian breed to the Thoroughbred racehorse, including Y chromosome ancestry. However, we found strong evidence for recent interbreeding of Thoroughbreds with Arabians used for flat-racing competitions. Genetic signatures suggestive of selective sweeps across the Arabian breed contain candidate genes for combating oxidative damage during exercise, and within the Straight Egyptian subgroup, for facial morphology. Overall, our data support an origin of the Arabian horse in the Middle East, no evidence for reduced global genetic diversity across the breed, and unique genetic adaptations for both physiology and conformation

Integrated analysis of inflammatory mRNAs, miRNAs, and lncRNAs elucidates the molecular interactome behind bovine mastitis

Abstract Mastitis is known as intramammary inflammation, which has a multifactorial complex phenotype. However, the underlying molecular pathogenesis of mastitis remains poorly understood. In this study, we utilized a combination of RNA-seq and miRNA-seq techniques, along with computational systems biology approaches, to gain a deeper understanding of the molecular interactome involved in mastitis. We retrieved and processed one hundred transcriptomic libraries, consisting of 50 RNA-seq and 50 matched miRNA-seq data, obtained from milk-isolated monocytes of Holstein–Friesian cows, both infected with Streptococcus uberis and non-infected controls. Using the weighted gene co-expression network analysis (WGCNA) approach, we constructed co-expressed RNA-seq-based and miRNA-seq-based modules separately. Module-trait relationship analysis was then performed on the RNA-seq-based modules to identify highly-correlated modules associated with clinical traits of mastitis. Functional enrichment analysis was conducted to understand the functional behavior of these modules. Additionally, we assigned the RNA-seq-based modules to the miRNA-seq-based modules and constructed an integrated regulatory network based on the modules of interest. To enhance the reliability of our findings, we conducted further analyses, including hub RNA detection, protein–protein interaction (PPI) network construction, screening of hub-hub RNAs, and target prediction analysis on the detected modules. We identified a total of 17 RNA-seq-based modules and 3 miRNA-seq-based modules. Among the significant highly-correlated RNA-seq-based modules, six modules showed strong associations with clinical characteristics of mastitis. Functional enrichment analysis revealed that the turquoise module was directly related to inflammation persistence and mastitis development. Furthermore, module assignment analysis demonstrated that the blue miRNA-seq-based module post-transcriptionally regulates the turquoise RNA-seq-based module. We also identified a set of different RNAs, including hub-hub genes, hub-hub TFs (transcription factors), hub-hub lncRNAs (long non-coding RNAs), and hub miRNAs within the modules of interest, indicating their central role in the molecular interactome underlying the pathogenic mechanisms of S. uberis infection. This study provides a comprehensive insight into the molecular crosstalk between immunoregulatory mRNAs, miRNAs, and lncRNAs during S. uberis infection. These findings offer valuable directions for the development of molecular diagnosis and biological therapies for mastitis

Study of whole genome linkage disequilibrium patterns of Iranian water buffalo breeds using the Axiom Buffalo Genotyping 90K Array.

Accuracy of genome-wide association studies, and the successful implementation of genomic selection depends on the level of linkage disequilibrium (LD) across the genome and also the persistence of LD phase between populations. In the present study LD between adjacent SNPs and LD decay between SNPs was calculated in three Iranian water buffalo populations. Persistence of LD phase was evaluated across these populations and effective population size (Ne) was estimated from corrected r2 information. A set of 404 individuals from three Iranian buffalo populations were genotyped with the Axiom Buffalo Genotyping 90K Array. Average r2 and |D'| between adjacent SNP pairs across all chromosomes was 0.27 and 0.66 for AZI, 0.29 and 0.68 for KHU, and 0.32 and 0.72 for MAZ. The LD between the SNPs decreased with increasing physical distance from 100Kb to 1Mb between markers, from 0.234 to 0.018 for AZI, 0.254 to 0.034 for KHU, and 0.297 to 0.119 for MAZ, respectively. These results indicate that a density of 90K SNP is sufficient for genomic analyses relying on long range LD (e.g. GWAS and genomic selection). The persistence of LD phase decreased with increasing marker distances across all the populations, but remained above 0.8 for AZI and KHU for marker distances up to 100Kb. For multi-breed genomic evaluation, the 90K SNP panel is suitable for AZI and KHU buffalo breeds. Estimated effective population sizes for AZI, KHU and MAZ were 477, 212 and 32, respectively, for recent generations. The estimated effective population sizes indicate that the MAZ is at risk and requires careful management

A genome-wide scan for signatures of selection in Azeri and Khuzestani buffalo breeds

Abstract Background Identification of genomic regions that have been targets of selection may shed light on the genetic history of livestock populations and help to identify variation controlling commercially important phenotypes. The Azeri and Kuzestani buffalos are the most common indigenous Iranian breeds which have been subjected to divergent selection and are well adapted to completely different regions. Examining the genetic structure of these populations may identify genomic regions associated with adaptation to the different environments and production goals. Results A set of 385 water buffalo samples from Azeri (N = 262) and Khuzestani (N = 123) breeds were genotyped using the Axiom® Buffalo Genotyping 90 K Array. The unbiased fixation index method (FST) was used to detect signatures of selection. In total, 13 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 regions. Putative candidate genes identified include FBXO9, NDFIP1, ACTR3, ARHGAP26, SERPINF2, BOLA-DRB3, BOLA-DQB, CLN8, and MYOM2. Conclusions Candidate genes identified in regions potentially under selection were associated with physiological pathways including milk production, cytoskeleton organization, growth, metabolic function, apoptosis and domestication-related changes include immune and nervous system development. The QTL identified are involved in economically important traits in buffalo related to milk composition, udder structure, somatic cell count, meat quality, and carcass and body weight

Estimation of economic values in three breeding perspectives for longevity and milk production traits in Holstein dairy cattle in Iran

The objectives of this study were to estimate economic values (EVs) for three production traits (milk, fat and protein yields) and longevity and to develop a national selection index. The proposed Iranian selection index was compared with selection indices of three other countries in the world. A simple and appropriate model was used to describe the Holstein dairy cattle industry under an Iranian production system. Production parameters and economic data were collected from two Holstein dairy farms in Tehran province. The EVs were estimated at farm level for three breeding perspectives (maximized profit, minimized costs, and economic efficiency) and two restrictions in production system (fixed herd size and fixed total input). The average absolute EVs on profit perspective and herd size restriction for milk, fat, and protein yields (based on $/kg) and longevity ($/month) were 0.11, 0.89, -0.20, and 6.20, respectively. The average absolute EVs under minimized costs per unit of product interest for milk, fat, protein yields and longevity were -0.30, -3.43, 0.88 and -20.40, respectively. The average absolute EVs under maximized economic efficiency for milk, fat and protein yields and longevity were 0.34, 2.73, -0.99 and 36.33, respectively. Relative emphasis for three production traits and longevity were 59.7, 14.3, -3.0 and 23.1, respectively. The comparison of the proposed Iranian index with those countries where most of the semen and embryos are imported points out that developing a national selection index to improve cow profitability and optimum generic trends is necessary. Sensitivity analysis indicated that the influence of milk payment changes on EVs was the greatest as its influence on fat and protein EVs is substantial. EVs for milk and fat yields, with respect to price changes (milk, feed and non-feed), were the least sensitive and most sensitive, respectively

Additional file 4: of A genome-wide scan for signatures of selection in Azeri and Khuzestani buffalo breeds

Haploview LD graph for selected regions. (PDF 904 kb

Additional file 12: of A genome-wide scan for signatures of selection in Azeri and Khuzestani buffalo breeds

FST output file. (XLS 9242 kb