Population Structure Analysis Reveals the Rich Genetic Diversity of Honeybee (Apis mellifera L.) Populations in Kazakhstan Kairat Dossybayev, Aidar Tapelov, Ulzhan Nuraliyeva, Gaukhar Moldakhmetova, Tilek Kapassuly, et al. Insects, 2026 Honeybee (Apis mellifera L.) populations are a vital resource for pollination and honey production, yet their genetic diversity in Central Asia remains poorly understood. This study provides a comprehensive genetic assessment of 16 honeybee populations from Kazakhstan, with comparative samples from Russia, Georgia and Kyrgyzstan, utilizing mitochondrial COI–COII intergenic region and 12 highly polymorphic nuclear STR markers. Mitochondrial DNA analysis revealed the predominance of the Eastern European C lineage (A. m. carnica), while a few populations from East Kazakhstan and Russia attributed the M lineage (A. m. mellifera), indicating local introgression and the persistence of relict lineages. STR analyses showed high levels of polymorphism and genetic diversity, with variation in heterozygosity and inbreeding across populations. Analyses of population genetic structure delineated four principal genetic clusters shaped by regional differentiation, historical gene flow, and sporadic admixture. Concordance between mitochondrial and nuclear markers confirms the robustness of these findings. Overall, this study highlights the rich genetic diversity of honeybees from Kazakhstan and emphasizes the importance of conserving local populations and implementing selective breeding programs to sustain adaptive potential and long-term apiculture.
Genome-Wide SNP Analysis Reveals the Unique Genetic Diversity Represented by Fat-Tailed Coarse-Wooled Sheep Breeds of Kazakhstan Kairat Dossybayev, Makpal Amandykova, Daniya Ualiyeva, Tilek Kapassuly, Altynay Kozhakhmet, et al. Biology, 2025 Background: The fat-tailed coarse-wooled sheep breeds exhibit excellent reproductive performance, exceptional adaptability to pasture conditions, and high precocity, contributing to enhanced meat, fat, and wool productivity in sheep breeding. Despite the significant role of these sheep breeds in Kazakhstan’s livestock production, their genetics remain poorly studied. This raises concerns about the potential loss of unique, breed-specific traits that could be important for the future development and resilience of Kazakh stan’s sheep farming sector. This study aimed to analyze genome-wide genotyping SNP data of local fat-tailed coarse-wooled sheep breeds (Kazakh fat-tailed coarse-wooled, Edilbay, and Gissar) to reveal their genetic diversity, breed characteristics, and phylogenetic relationships with worldwide domestic sheep breeds and wild sheep. Methods: The OvineSNP50 Genotyping BeadChip was used to obtain genome-wide SNP genotyping data from 160 fat-tailed coarse-wooled sheep from Kazakhstan. Population structure analysis, principal component analysis, phylogenetic and the maximum likelihood tree analysis were performed in comparison with foreign domestic sheep breeds and wild sheep populations. Results: Kazakh breeds exhibited high genetic diversity, with Edilbay showing the greatest allelic richness. PCA and Admixture revealed clear differentiation among the three breeds: Edilbay and Gissar formed homogeneous clusters, while Kazakh fat-tailed coarse-wooled sheep displayed admixture and substructure. Evidence of gene flow from Edilbay into other Kazakh populations supports its role as a genetic source for regional breeds. Phylogenetic analysis placed Kazakhstani sheep close to other Central Asian breeds, while clearly distinct from East Asian and European populations. Wild sheep (Argali and Urial) formed separate clades, with Kerman wild sheep clustering closer to Urial. Conclusions: Our results highlight the value of genotyping data for studying genetic diversity and population structure. Developing genetic resources for Kazakhstan’s native sheep breeds will help preserve their unique diversity and ensure it remains available for future use in breeding and adaptation efforts.
Genetic Diversity and Population Structure of Fat-Tailed Coarse-Wooled Sheep Breeds Ovis aries from Kazakhstan Kairat Dossybayev, Daniya Ualiyeva, Tilek Kapassuly, Makpal Amandykova, Altynay Kozhahmet, et al. Veterinary Sciences, 2025 Sheep play a central role in Kazakhstan’s pastoral economy, yet the maternal genetic composition of its traditional breeds remains poorly characterized. We analyzed partial mitochondrial D-loop sequences (848 bp) from 115 individuals of three fat-tailed coarse-wooled breeds (Edilbay, Kazakh fat-tailed coarse-wooled, and Gissar) to assess genetic diversity, population structure, and phylogenetic relationships. Ninety-eight haplotypes were identified, indicating high haplotype diversity (Hd = 0.996 ± 0.002) and moderate nucleotide diversity (π = 0.02624 ± 0.00048). Haplotypes clustered into haplogroups A (57.4%) and B (42.6%), with Edilbay dominating the star-like cluster of haplogroup A, consistent with recent expansion. AMOVA revealed that most variation (92.03%) occurred within populations, with no significant differentiation among breeds. Phylogenetic analyses placed Edilbay close to the most recent common ancestor of fat-tailed domestic sheep and the wild Ovis species, suggesting retention of an ancestral lineage. These findings highlight Kazakhstan as a genetic crossroads in sheep history and underscore the conservation value of its maternal diversity.
Genetic Variability and Population Structure of Camelus from Kazakhstan Inferred from 17 STR Markers Gulfairuz Shaltenbay, Daniya Ualiyeva, Tilek Kapassuly, Altynay Kozhakhmet, Zarina Orazymbetova, et al. Diversity, 2025 Camels have been essential to human survival and development across the arid Central Asian steppes, particularly in Kazakhstan, where the breeding of one-humped and two-humped camels is a longstanding tradition supporting the nomadic lifestyle. This study aimed to assess the genetic diversity and population structure of these camels across their distribution range in Kazakhstan. Blood samples from 100 individuals were collected from five locations, Almaty (ALA), Atyrau (ATR), Shymkent (SHK), Kyzylorda (KZL), and Taraz (TRZ), and genotyped using 17 microsatellite markers. All loci were polymorphic, with a mean observed heterozygosity of 0.707 in C. dromedarius and 0.643 in C. bactrianus. The highest expected heterozygosity (He = 0.939) was observed at VOLP67 in C. bactrianus and at VOLP03 in C. dromedarius. Genetic differentiation was low (FST = 0.021), indicating a weak population structure between the two species with substantial gene flow (Nm = 19.972). The hybrid analysis identified 31% hybrids, including F1, F2, and backcrosses, with the highest frequencies in KZL and TRZ, moderate frequencies in ATR, and lowest frequencies in SHK and ALA. These patterns, consistent with STRUCTURE clustering, reflect widespread but regionally variable hybridization. The phylogenetic analysis revealed three clades, separating Bactrian camels (ALA), dromedaries (SHK), and a hybrid group (ATR, KZL, and TRZ). These findings enhance our understanding of the genetic diversity of Kazakhstan’s camels and support effective conservation, breeding strategies, and genotyping applications in camel husbandry.
Genome-Wide Association Studies Revealed Several Candidate Genes of Meat Productivity in Saryarka Fat-Tailed Coarse-Wool Sheep Breed Kairat Dossybayev, Makpal Amandykova, Ainur Orakbayeva, Sholpan Adylkanova, Altynay Kozhakhmet, et al. Genes, 2024 Background: Saryarka sheep belong to fat-tailed coarse-wool sheep breed. This breed is distinguished by increased meat productivity while being competitive in young lamb production. Live weight and body indices are relevant data for assessing sheep body constitution, which directly affects the breeding characteristics and meat productivity of animals. Objectives: This study aimed to find associations with SNPs and nine phenotypic characteristics of the Saryarka fat-tailed coarse-wool sheep breed including live weight and eight body indices (wither height, rump height, bicoastal diameter, body depth, body length, rump width, heart girth, and cannon bone circumference), and find candidate genes related to these characteristics. Methods: A total of 100 animals from the Karaganda region of Kazakhstan were used in this study. Live weight and eight body indices of sheep were measured using tape and electronic scales. The blood samples of the animals were used for DNA extraction. DNA samples were genotyped with the OvineSNP50 Genotyping BeadChip and analyzed using GWAS. Statistically significant SNPs were identified for each characteristic trait referencing the genome of Ovis aries (Oar_v3.1) using BioMart. Results: The GWAS results demonstrated a substantial chromosomal-level correlation between 32 chromosome-wide significant and suggestively significant SNPs in the studied sheep breed. Overall, seven SNPs located in seven different genes were revealed as candidates for live weight and four body indices: s20793.1 SNP in the IGFBP6 gene for live weight, OAR4_54217431.1 SNP in the ST7 gene for bicoastal diameter, s25229.1 in the SCD5 gene, and s01175.1 SNP in the DTNBP1 gene for rump width, OAR2_175574781.1 SNP in the KYNU gene for heart girth, and OAR1_209022621.1 SNP in the FGF12 gene and s15415.1 SNP in the FTO gene for cannon bone circumference. Some of these genes were previously reported to be involved in body constitution and fat deposit in other sheep breeds. Conclusions: The results of the present study open up new opportunities for targeted sheep breeding for meat and fat productivity.
Detection of genetic variations in the GDF9 and BMP15 genes in Kazakh meat–wool sheep Makpal Amandykova, Zarina Orazymbetova, Tilek Kapassuly, Altynay Kozhakhmet, Saltanat Khamzina, et al. Archives Animal Breeding, 2023 Kazakh meat–wool sheep are of great interest because of the intrabreed multifetal type's high productivity of 140 %–160 %. Genes encoding growth differentiation factor-9 (GDF9) and bone morphogenetic protein 15 (BMP15) are promising candidates for studying sheep productivity, as they affect fertility in mammals, including sheep. Thus, the purpose of this study was to assess the fertility of the Kazakh meat–wool sheep breed based on GDF9 and BMP15 candidate genes of fecundity for the selection of animals with valuable genotypes. We selected 300 heads of the Kazakh meat–wool sheep breed from two populations for PCR-RFLP (polymerase chain reaction–restriction fragment length polymorphism) analysis, 15 of which were subsequently used for sequencing of exon regions of the GDF9 and BMP15 genes. The sheep populations were tested for G1 and G8 mutations of the GDF9 gene and B2 and B4 mutations of the BMP15 gene. The PCR-RFLP analysis revealed that 59 (19.7 %) of the 300 Kazakh meat–wool breed sheep were heterozygous carriers of the G1 mutation (genotype AG) of the GDF9 gene, and sequencing analysis supported these results. The comparative phylogenetic analysis showed a clear separation of Kazakh meat–wool sheep wild types and carriers of the G1 mutation. This mutation was reported to have a relationship with the animals' litter size in other sheep breeds. For this reason, similar relationships should be investigated in Kazakh meat–wool sheep. However, G8, B2, and B4 mutations were not detected among the studied animal populations, showing that these mutations are not characteristic of the Kazakh meat–wool sheep breed.
Distribution of Runs of Homozygosity and Their Relationship with Candidate Genes for Productivity in Kazakh Meat–Wool Sheep Breed Makpal Amandykova, Zhanerke Akhatayeva, Altynay Kozhakhmet, Tilek Kapassuly, Zarina Orazymbetova, et al. Genes, 2023 Increasing the fertility of sheep remains one of the crucial issues of modern sheep breeding. The Kazakh meat–wool sheep is an excellent breed with high meat and wool productivity and well adapted to harsh conditions. Nowadays, runs of homozygosity (ROHs) are considered a suitable approach for studying the genetic characteristics of farm animals. The aims of the study were to analyze the distribution of ROHs, describe autozygosity, and detect genomic regions with high ROH islands. In this study, we genotyped a total of 281 Kazakh meat–wool sheep using the Illumina iScan® system (EquipNet, Canton, MA, USA) via Ovine SNP50 BeadChip array. As a results, a total of 15,069 ROHs were found in the three Kazakh meat–wool sheep populations. The mean number of ROH per animal across populations varied from 40.3 (POP1) to 42.2 (POP2) in the category 1+ Mb. Furthermore, the number of ROH per animal in ROH1–2 Mb were much higher than ROH2–4 Mb and ROH8–16 Mb in the three sheep populations. Most of individuals had small number of ROH>16 Mb. The highest and lowest genomic inbreeding coefficient values were observed in POP2 and POP3, respectively. The estimated FROH presented the impact that recent inbreeding has had in all sheep populations. Furthermore, a set of interesting candidate genes (BMP2, BMPR2, BMPRIB, CLOCK, KDM2B, TIAM1, TASP1, MYBPC1, MYOM1, and CACNA2D1), which are related to the productive traits, were found. Collectively, these findings will contribute to the breeding and conservation strategies of the Kazakh meat–wool sheep breed.
A Study of the Genetic Structure of Hybrid Camels in Kazakhstan Makpal Amandykova, Kairat Dossybayev, Aizhan Mussayeva, Naruya Saitou, Zhazira Zhunusbayeva, et al. Genes, 2023 Camel farming is gaining scientific interest due to its unique agricultural characteristics. Camels are versatile for milk and meat production, wool, racing, transport, and tourism. To use their full potential, it is essential to improve our understanding of the genetic structure of these animals. One-humped and two-humped camels have received detailed genetic descriptions, while there is no such information for their hybrids, which outperform their parent species in several agricultural characteristics. Thus, in this study, for the first time, the whole genome sequencing data (WGS) of five hybrid camels bred in the Almaty region of Kazakhstan are presented in comparison with the WGS data of one-humped, two-humped, and wild camels. A total of 43,552,164 single-nucleotide polymorphisms were found across the studied groups. Further comparison of these SNPs showed the following number of private SNPs among the populations: hybrid camels (3,271,083), wild camels (2,515,591), Bactrians (1,244,694), and dromedaries (531,224). The genetic structure of the studied animals was described, and a phylogenetic tree was built to assess their genetic distance. It was found that the studied hybrids are genetically closer to dromedaries since they were on the close branch of the phylogenetic tree.
Comparative Analysis of the Polymorphism of the Casein Genes in Camels Bred in Kazakhstan Makpal Amandykova, Kairat Dossybayev, Aizhan Mussayeva, Bakytzhan Bekmanov, Naruya Saitou Diversity, 2022 Caseins play an important role in determining the technological properties and quantitative characteristics of camel milk. To date, only a few studies on the genetic polymorphism of casein genes have been reported in the camel populations of Kazakhstan. Therefore, this work aimed to identify the genetic polymorphism level of casein genes among camel populations of the Almaty region of Kazakhstan. The PCR-RFLP method was used for this purpose and the following genotypes were revealed as a result: CSN3 gene—CC, CT, TT, where the T allele predominated in all populations, with a frequency of 0.60; CSN2 gene—AA, AG, GG, with the predomination of A allele (0.64); and CSN1S1 gene—GG and GT, with the predomination of G allele (0.94). Statistical analysis was carried out using the POPGENE and GenAlEx software. The χ2 values were equal to 12.1 (CSN3), 8.6 (CSN2), and 14.5 (CSN1S1). As a result, three out of 53 animals were designated as the “core” of the population—animals with the desired genotypes: CC genotype for the CSN3 gene and AA genotype for the CSN2 gene. Such animals can be selected for further use with an increase in the number of livestock with high productivity rates.
Ecological risk assessment and long-term environmental pollution caused by obsolete undisposed organochlorine pesticides Natalya Mit, Oksana Cherednichenko, Aizhan Mussayeva, Ozada Khamdiyeva, Almira Amirgalieva, et al. Journal of Environmental Science and Health Part B Pesticides Food Contaminants and Agricultural Wastes, 2021 Obsolete organochlorine pesticides (OSPs) are currently prohibited as persistent organic pollutants that contaminate the environment. If undisposed, they continue to pollute soil and water, to accumulate in the food chain and to harm plants, animals and the human body. The aim of the study was to assess water and soil pollution around the storehouses of undisposed, banned OSPs and their possible genotoxic effect. The storehouses in four villages near Almaty, Kazakhstan were investigated. Chemical analysis confirmed contamination of water and soil around storehouses with OSPs. The genotoxic effect of water and soil samples was evaluated using model objects: S.typhymurium, D.melanogaster, sheep lymphocytes cultures and human lymphocytes cultures. It was found that water and soil samples caused mutagenic effect in all model systems. They increased the frequency of revertants in Salmonella, the frequency of lethal mutations in Drosophila chromosomes, and the frequency of chromosome aberrations in cultures of human and sheep lymphocytes. Although a genotoxic effect was demonstrated for each of these models, various models showed different sensitivity to the effects of pesticides and they varied degree of response. The association between the total content of OCPs in soil and the level of mutations for different model systems was discovered.
