Search

Abstract Background To enhance and extend the knowledge about the global historical and phylogenetic relationships between Merino and Merino-derived breeds, 19 populations were genotyped with the OvineSNP50 BeadChip specifically for this study, while an additional 23 populations from the publicly available genotypes were retrieved. Three complementary statistical tests, Rsb (extended haplotype homozygosity between-populations), XP-EHH (cross-population extended haplotype homozygosity), and runs of homozygosity (ROH) islands were applied to identify genomic variants with potential impact on the adaptability of Merino genetic type in two contrasting climate zones. Results The results indicate that a large part of the Merino’s genetic relatedness and admixture patterns are explained by their genetic background and/or geographic origin, followed by local admixture. Multi-dimensional scaling, Neighbor-Net, Admixture, and TREEMIX analyses consistently provided evidence of the role of Australian, Rambouillet and German strains in the extensive gene introgression into the other Merino and Merino-derived breeds. The close relationship between Iberian Merinos and other South-western European breeds is consistent with the Iberian origin of the Merino genetic type, with traces from previous contributions of other Mediterranean stocks. Using Rsb and XP-EHH approaches, signatures of selection were detected spanning four genomic regions located on Ovis aries chromosomes (OAR) 1, 6 and 16, whereas two genomic regions on OAR6, that partially overlapped with the previous ones, were highlighted by ROH islands. Overall, the three approaches identified 106 candidate genes putatively under selection. Among them, genes related to immune response were identified via the gene interaction network. In addition, several candidate genes were found, such as LEKR1, LCORL, GHR, RBPJ, BMPR1B, PPARGC1A, and PRKAA1, related to morphological, growth and reproductive traits, adaptive thermogenesis, and hypoxia responses. Conclusions To the best of our knowledge, this is the first comprehensive dataset that includes most of the Merino and Merino-derived sheep breeds raised in different regions of the world. The results provide an in-depth picture of the genetic makeup of the current Merino and Merino-derived breeds, highlighting the possible selection pressures associated with the combined effect of anthropic and environmental factors. The study underlines the importance of Merino genetic types as invaluable resources of possible adaptive diversity in the context of the occurring climate changes.

The pig industry is usually considered an intensive livestock industry, mainly supported by hybrid breeding between commercial pig breeds. However, people’s pursuit of a more natural environment and higher meat quality has led to an increasing demand for eco-friendly and diverse pig feeding systems. Therefore, the importance of rearing and conserving local pig breeds is increasing. The Livni pig is a local breed with good adaptability to the environmental and fodder conditions in central Russia. In this study, we aimed to analyze the genetic diversity and population structure of Livni pigs using whole-genome single nucleotide polymorphism (SNP) data. We utilized the Porcine GGP HD BeadChip on genotype samples from old (n = 32, 2004) and modern (n = 32, 2019) populations of Livni pigs. For the museum samples of Livni pigs (n = 3), we extracted DNA from their teeth, performed genomic sequencing, and obtained SNP genotypes from the whole-genome sequences. SNP genotypes of Landrace (n = 32) and Large White (n = 32) pigs were included for comparative analysis. We observed that the allelic richness of Livni pigs was higher than those of Landrace and Large White pigs (AR = 1.775–1.798 vs. 1.703 and 1.668, respectively). The effective population size estimates (NE5 = 108 for Livni pigs, NE5 = 59 for Landrace and Large White pigs) confirmed their genetic diversity tendency. This was further supported by the length and number of runs of homozygosity, as well as the genomic inbreeding coefficient (almost twofold lower in Livni pigs compared to Landrace and Large White pigs). These findings suggest that the Livni pig population exhibits higher genetic diversity and experiences lower selection pressure compared to commercial pig populations. Furthermore, both principal component and network tree analyses demonstrated a clear differentiation between Livni pigs and transboundary commercial pigs. The TreeMix results indicated gene flow from Landrace ancestors to Livni pigs (2019) and from Large White ancestors to Livni pigs (2004), which was consistent with their respective historical breeding backgrounds. The comparative analysis of museum, old, and modern Livni pigs indicated that the modern Livni pig populations have preserved their historical genomic components, suggesting their potential suitability for future design selection programs.

Abstract Background The genomes of worldwide poultry breeds divergently selected for performance and other phenotypic traits may also be affected by, and formed due to, past and current admixture events. Adaptation to diverse environments, including acclimation to harsh climatic conditions, has also left selection footprints in breed genomes. Results Using the Chicken 50K_CobbCons SNP chip, we genotyped four divergently selected breeds: two aboriginal, cold tolerant Ushanka and Orloff Mille Fleur, one egg-type Russian White subjected to artificial selection for cold tolerance, and one meat-type White Cornish. Signals of selective sweeps were determined in the studied breeds using three methods: (1) assessment of runs of homozygosity islands, (2) F ST based population differential analysis, and (3) haplotype differentiation analysis. Genomic regions of true selection signatures were identified by two or more methods or in two or more breeds. In these regions, we detected 540 prioritized candidate genes supplemented them with those that occurred in one breed using one statistic and were suggested in other studies. Amongst them, SOX5, ME3, ZNF536, WWP1, RIPK2, OSGIN2, DECR1, TPO, PPARGC1A, BDNF, MSTN, and beta-keratin genes can be especially mentioned as candidates for cold adaptation. Epigenetic factors may be involved in regulating some of these important genes (e.g., TPO and BDNF). Conclusion Based on a genome-wide scan, our findings can help dissect the genetic architecture underlying various phenotypic traits in chicken breeds. These include genes representing the sine qua non for adaptation to harsh environments. Cold tolerance in acclimated chicken breeds may be developed following one of few specific gene expression mechanisms or more than one overlapping response known in cold-exposed individuals, and this warrants further investigation.

The aim of the current study was to assess the female metabolic rate and test the hypothesis that there is a relationship between the egg productivity of Japanese quails from eight breeds and their morphometric, or growth, parameters. Parameters measured were body weight (B), volume (V), and surface area (S), as well as the metabolism level expressed by the ratio S/V. The collected egg performance traits were as follows: the number of eggs produced (N), the average egg weight (W), and the total egg mass (M) (i.e., N multiplied by W). To measure the S and V values, a novel technique was developed that takes into account the similarity of the quail’s body to an ellipsoid. An analysis of the relationships between productivity indicators allowed us to introduce a new index called the metabolic index, B·S/V, based on all three main growth parameters in quails. Using the values of this index, we were then able to judge indirectly the level of quails’ egg productivity. We went on to assess the N, W, and M values, not only depending on the size of the bird’s growth parameters but also according to the degree of their changes during quail growth. These changes were expressed as the slope angles of trend lines describing the growth process data. This approach produced more accurate results for predicting the egg productivity in terms of W and M.

Traces of long-term artificial selection can be detected in genomes of domesticated birds via whole-genome screening using single-nucleotide polymorphism (SNP) markers. This study thus examined putative genomic regions under selection that are relevant to the development history, divergence and phylogeny among Japanese quails of various breeds and utility types. We sampled 99 birds from eight breeds (11% of the global gene pool) of egg (Japanese, English White, English Black, Tuxedo and Manchurian Golden), meat (Texas White and Pharaoh) and dual-purpose (Estonian) types. The genotyping-by-sequencing analysis was performed for the first time in domestic quails, providing 62,935 SNPs. Using principal component analysis, Neighbor-Net and Admixture algorithms, the studied breeds were characterized according to their genomic architecture, ancestry and direction of selective breeding. Japanese and Pharaoh breeds had the smallest number and length of homozygous segments indicating a lower selective pressure. Tuxedo and Texas White breeds showed the highest values of these indicators and genomic inbreeding suggesting a greater homozygosity. We revealed evidence for the integration of genomic and performance data, and our findings are applicable for elucidating the history of creation and genomic variability in quail breeds that, in turn, will be useful for future breeding improvement strategies.

A study for genomic variation that may reflect putative selective signaling and be associated with economically important traits is instrumental for obtaining information about demographic and selection history in domestic animal species and populations. A rich variety of the Russian chicken gene pool breeds warrants a further detailed study. Specifically, their genomic features can derive implications from their genome architecture and selective footprints for their subsequent breeding and practical efficient exploitation. In the present work, whole genome genotyping of 19 chicken breeds (20 populations with up to 71 samples each) was performed using the Chicken 50 K BeadChip DNA chip. The studied breed sample included six native Russian breeds of chickens developed in the 17th–19th centuries, as well as eight Russian chicken breeds, including the Russian White (RW), created in the 20th century on the basis of improving local chickens using breeds of foreign selection. Five specialized foreign breeds of chickens, including the White Leghorn (WL), were used along with other breeds representing the Russian gene pool. The characteristics of the genetic diversity and phylogenetic relationships of the native breeds of chickens were represented in comparison with foreign breeds. It was established that the studied native breeds demonstrate their own genetic structure that distinguishes them from foreign breeds, and from each other. For example, we previously made an assumption on what could cause the differences between two RW populations, RW1 and RW2. From the data obtained here, it was verified that WL was additionally crossed to RW2, unlike RW1. Thus, inherently, RW1 is a purer population of this improved Russian breed. A significant contribution of the gene pool of native breeds to the global genetic diversity of chickens was shown. In general, based on the results of a multilateral survey of this sample of breeds, it can be concluded that phylogenetic relationships based on their genetic structure and variability robustly reflect the known, previously postulated and newly discovered patterns of evolution of native chickens. The results herein presented will aid selection and breeding work using this gene pool.

Negative heterosis can occur on different economically important traits, but the exact biological mechanisms of this phenomenon are still unknown. The present study focuses on determining the genetic factors associated with negative heterosis in interspecific hybrids between domestic sheep (Ovis aries) and argali (Ovis ammon). One locus (rs417431015) associated with viability and two loci (rs413302370, rs402808951) associated with meat productivity were identified. One gene (ARAP2) was prioritized for viability and three for meat productivity (PDE2A, ARAP1, and PCDH15). The loci associated with meat productivity were demonstrated to fit the overdominant inheritance model and could potentially be involved int negative heterosis mechanisms.

Operating under the assumption that decline in population size increases the threat of loss of genetic diversity, herein, our first aim was to monitor and document the genetic diversity and population structure of the endemic Tuva–Tofalar reindeer, the number of which has seen a dramatic decrease. Second, we were interested in understanding of Tuva reindeer genetic relationship with four officially recognized reindeer breeds, whose wild populations inhabit the Far North of Russia, as well as with the domestic reindeer making up the Mongolian population. Our results, based on the genome-wide SNP genotypes generated by the BovineHD BeadChip, revealed a low level of genetic variability of Tuva population in terms of the values of both allelic and genetic diversity. The ADMIXTURE analysis, the PCA plot, and the Neighbor Net network results, showed that Tuva population shared genetic background with reindeer inhabiting Mongolia, an implication of their common origin. Moreover, all statistic approaches used in our study showed a distinctive shared genetic structure revealed in independent clusters found in the composition of: an admixed cluster of Tuvan and Mongolian populations, a cluster of domestic reindeer breeds, and a cluster of the wild populations—all of which indirectly points to the possibility of the independent origins (with regard to the domestication) of the reindeer groups studied herein. We believe that our findings will contribute to the formation of a rational basis for solving problems related to the conservation of domestic Tuva-Tofalar reindeer in order to increase the number of this substantial reindeer population which will contribute to the sustainable functioning of ecosystems and the lives and culture of indigenous peoples inhabiting the southern sections of Eastern Siberia.

Indigenous pig breeds are mainly associated with the adaptive capacity that is necessary to respond adequately to climate change, food security, and livelihood needs, and natural resources conservation. Livni pigs are an indigenous fat-type breed farmed in a single farm in the Orel region and located in the Central European part of the Russian Federation. To determine the genomic regions and genes that are affected by artificial selection, we conducted the comparative study of two pig breeds with different breeding histories and breeding objectives, i.e., the native fat-type Livni and meat-type Duroc breeds using the Porcine GGP HD BeadChip, which contains ~80,000 SNPs. To check the Livni pigs for possible admixture, the Landrace and the Large White breeds were included into the study of genetic diversity as these breeds participated in the formation of the Livni pigs. We observed the highest level of genetic diversity in Livni pigs compared to commercial breeds (UHE = 0.409 vs. 0.319–0.359, p < 0.001; AR = 1.995 vs. 1.894–1.964, p < 0.001). A slight excess of heterozygotes was found in all of the breeds. We identified 291 candidate genes, which were localized within the regions under putative selection, including 22 and 228 genes, which were specific for Livni and Duroc breeds, respectively, and 41 genes common for both breeds. A detailed analysis of the molecular functions identified the genes, which were related to the formation of meat and fat traits, and adaptation to environmental stress, including extreme temperatures, which were different between breeds. Our research results are useful for conservation and sustainable breeding of Livni breed, which shows a high level of genetic diversity. This makes Livni one of the valuable national pig genetic resources.

Specific local environmental and sociocultural conditions have led to the creation of various goat populations in Russia. National goat diversity includes breeds that have been selected for down and mohair production traits as well as versatile local breeds for which pastoralism is the main management system. Effective preservation and breeding programs for local goat breeds are missing due to the lack of DNA-based data. In this work, we analyzed the genetic diversity and population structure of Russian local goats, including Altai Mountain, Altai White Downy, Dagestan Downy, Dagestan Local, Karachaev, Orenburg, and Soviet Mohair goats, which were genotyped with the Illumina Goat SNP50 BeadChip. In addition, we addressed genetic relationships between local and global goat populations obtained from the AdaptMap project. Russian goats showed a high level of genetic diversity. Although a decrease in historical effective population sizes was revealed, the recent effective population sizes estimated for three generations ago were larger than 100 in all studied populations. The mean runs of homozygosity (ROH) lengths ranged from 79.42 to 183.94 Mb, and the average ROH number varied from 18 to 41. Short ROH segments (16 Mb) was the least frequent. Principal component analysis, Neighbor-Net graph, and Admixture clustering revealed several patterns in Russian local goats. First, a separation of the Karachaev breed from other populations was observed. Moreover, genetic connections between the Orenburg and Altai Mountain breeds were suggested and the Dagestan breeds were found to be admixed with the Soviet Mohair breed. Neighbor-Net analysis and clustering of local and global breeds demonstrated the close genetic relations between Russian local and Turkish breeds that probably resulted from past admixture events through postdomestication routes. Our findings contribute to the understanding of the genetic relationships of goats originating in West Asia and Eurasia and may be used to design breeding programs for local goats to ensure their effective conservation and proper management.

Abstract Background Domestication and centuries of selective breeding have changed genomes of sheep breeds to respond to environmental challenges and human needs. The genomes of local breeds, therefore, are valuable sources of genomic variants to be used to understand mechanisms of response to adaptation and artificial selection. As a step toward this we performed a high-density genotyping and comprehensive scans for signatures of selection in the genomes from 15 local sheep breeds reared across Russia. Results Results demonstrated that the genomes of Russian sheep breeds contain multiple regions under putative selection. More than 50% of these regions matched with intervals identified in previous scans for selective sweeps in sheep genomes. These regions contain well-known candidate genes related to morphology, adaptation, and domestication (e.g., KITLG, KIT, MITF, and MC1R), wool quality and quantity (e.g., DSG@, DSC@, and KRT@), growth and feed intake (e.g., HOXA@, HOXC@, LCORL, NCAPG, LAP3, and CCSER1), reproduction (e.g., CMTM6, HTRA1, GNAQ, UBQLN1, and IFT88), and milk-related traits (e.g., ABCG2, SPP1, ACSS1, and ACSS2). In addition, multiple genes that are putatively related to environmental adaptations were top-ranked in selected intervals (e.g., EGFR, HSPH1, NMUR1, EDNRB, PRL, TSHR, and ADAMTS5). Moreover, we observed that multiple key genes involved in human hereditary sensory and autonomic neuropathies, and genetic disorders accompanied with an inability to feel pain and environmental temperatures, were top-ranked in multiple or individual sheep breeds from Russia pointing to a possible mechanism of adaptation to harsh climatic conditions. Conclusions Our work represents the first comprehensive scan for signatures of selection in genomes of local sheep breeds from the Russian Federation of both European and Asian origins. We confirmed that the genomes of Russian sheep contain previously identified signatures of selection, demonstrating the robustness of our integrative approach. Multiple novel signatures of selection were found near genes which could be related to adaptation to the harsh environments of Russia. Our study forms a basis for future work on using Russian sheep genomes to spot specific genetic variants or haplotypes to be used in efforts on developing next-generation highly productive breeds, better suited to diverse Eurasian environments.

Cattle breeds used in industrial production tend to be larger than local cattle, as increased dairy and beef productivity is closely related to stature. This study aimed to identify single nucleotide polymorphisms (SNPs) significantly associated with stature in diverse cattle breeds. Thirteen local and transboundary cattle breeds (n = 670) were divided into two groups according to their stature. The high-stature group included the Angus, Ayrshire, Black-and-White, Holstein, Kholmogor, Yaroslavl, Tagil, and Istoben breeds. The low-stature group included Jersey, Kalmyk, Kyrgyz, Mongolian, and Yakut. The average height at withers was 136.3 ± 1.6 and 121.6 ± 2.8 cm in the high- and low-stature groups, respectively. The samples of the 11 breeds were genotyped using high-density DNA arrays. Genotypes of the two remaining breeds were downloaded from the publicly available WIDDE database. Genome-wide association studies revealed seven SNPs strongly associated with stature, including three at 77.3–77.8 cM on BTA4 and four at 24.5–25.2 cM on BTA14. Functional annotation showed the localization within identified regions of genes responsible for growth, exterior characteristics, protein and lipid metabolism, and feed intake. The identified SNPs can be considered useful DNA markers for marker-assisted cattle breeding aimed at increasing stature.

The Dagestan Mountain (DM) cattle breed was established to provide cheese and milk products to mountain dwellers in the specific conditions of the Republic of Dagestan in Southern Russia. Only 650 head of DM cattle were registered in 2020, and their “breed purity” is questionable. We aimed to assess the genetic diversity and population structure of modern DM cattle using short tandem repeat (STR) markers. The sample included 150 animals collected from private owners in Dagestan during a scientific expedition (n = 32) and provided by a gene pool farm (n = 118). An additional 166 samples from other cattle breeds distributed in the breeding zone of DM cattle were used as the comparison groups. The genotypes for the 11 STR loci recommended by ISAG were obtained using a genetic analyzer. We identified a high level of genetic diversity as revealed by allelic richness calculations (Ar = 6.82 vs. 4.38–5.82 in other cattle breeds) and observed heterozygosity indices (uHe = 0.76 vs. 0.65–0.72). Based on the STRUCTURE results, animals with low levels of admixture with other breeds were found within the DM cattle, which can be considered as candidates for use in germplasm conservation programs.

Native cattle breeds can carry specific signatures of selection reflecting their adaptation to the local environmental conditions and response to the breeding strategy used. In this study, we comprehensively analysed high-density single nucleotide polymorphism (SNP) genotypes to characterise the population structure and detect the selection signatures in Russian native Yaroslavl and Kholmogor dairy cattle breeds, which have been little influenced by introgression with transboundary breeds. Fifty-six samples of pedigree-recorded purebred animals, originating from different breeding farms and representing different sire lines, of the two studied breeds were genotyped using a genome-wide bovine genotyping array (Bovine HD BeadChip). Three statistical analyses-calculation of fixation index (FST) for each SNP for the comparison of the pairs of breeds, hapFLK analysis, and estimation of the runs of homozygosity (ROH) islands shared in more than 50% of animals-were combined for detecting the selection signatures in the genome of the studied cattle breeds. We confirmed nine and six known regions under putative selection in the genomes of Yaroslavl and Kholmogor cattle, respectively; the flanking positions of most of these regions were elucidated. Only two of the selected regions (localised on BTA 14 at 24.4-25.1 Mbp and on BTA 16 at 42.5-43.5 Mb) overlapped in Yaroslavl, Kholmogor and Holstein breeds. In addition, we detected three novel selection sweeps in the genome of Yaroslavl (BTA 4 at 4.74-5.36 Mbp, BTA 15 at 17.80-18.77 Mbp, and BTA 17 at 45.59-45.61 Mbp) and Kholmogor breeds (BTA 12 at 82.40-81.69 Mbp, BTA 15 at 16.04-16.62 Mbp, and BTA 18 at 0.19-1.46 Mbp) by using at least two of the above-mentioned methods. We expanded the list of candidate genes associated with the selected genomic regions and performed their functional annotation. We discussed the possible involvement of the identified candidate genes in artificial selection in connection with the origin and development of the breeds. Our findings on the Yaroslavl and Kholmogor breeds obtained using high-density SNP genotyping and three different statistical methods allowed the detection of novel putative genomic regions and candidate genes that might be under selection. These results might be useful for the sustainable development and conservation of these two oldest Russian native cattle breeds.

Abstract Insights into the genetic characteristics of a species provide important information for wildlife conservation programs. Here, we used the OvineSNP50 BeadChip developed for domestic sheep to examine population structure and evaluate genetic diversity of snow sheep (Ovis nivicola) inhabiting Verkhoyansk Range and Momsky Ridge. A total of 1,121 polymorphic SNPs were used to test 80 specimens representing five populations, including four populations of the Verkhoyansk Mountain chain: Kharaulakh Ridge–Tiksi Bay (TIK, n = 22), Orulgan Ridge (ORU, n = 22), the central part of Verkhoyansk Range (VER, n = 15), Suntar‐Khayata Ridge (SKH, n = 13), and Momsky Ridge (MOM, n = 8). We showed that the studied populations were genetically structured according to a geographic pattern. Pairwise FST values ranged from 0.044 to 0.205. Admixture analysis identified K = 2 as the most likely number of ancestral populations. A Neighbor‐Net tree showed that TIK was an isolated group related to the main network through ORU. TreeMix analysis revealed that TIK and MOM originated from two different ancestral populations and detected gene flow from MOM to ORU. This was supported by the f3 statistic, which showed that ORU is an admixed population with TIK and MOM/SKH heritage. Genetic diversity in the studied groups was increasing southward. Minimum values of observed (Ho) and expected (He) heterozygosity and allelic richness (Ar) were observed in the most northern population—TIK, and maximum values were observed in the most southern population—SKH. Thus, our results revealed clear genetic structure in the studied populations of snow sheep and showed that TIK has a different origin from MOM, SKH, and VER even though they are conventionally considered a single subspecies known as Yakut snow sheep (Ovis nivicola lydekkeri). Most likely, TIK was an isolated group during the Late Pleistocene glaciations of Verkhoyansk Range.

Abstract Background The origin of native and locally developed Russian cattle breeds is linked to the historical, social, cultural, and climatic features of the diverse geographical regions of Russia. In the present study, we investigated the population structure of nine Russian cattle breeds and their relations to the cattle breeds from around the world to elucidate their origin. Genotyping of single nucleotide polymorphisms (SNPs) in Bestuzhev (n = 26), Russian Black-and-White (n = 21), Kalmyk (n = 14), Kholmogor (n = 25), Kostromsky (n = 20), Red Gorbatov (n = 23), Suksun (n = 20), Yakut (n = 25), and Yaroslavl cattle breeds (n = 21) was done using the Bovine SNP50 BeadChip. SNP profiles from an additional 70 breeds were included in the analysis as references. Results The observed heterozygosity levels were quite similar in eight of the nine studied breeds (HO = 0.337–0.363) except for Yakut (Ho = 0.279). The inbreeding coefficients F IS ranged from -0.028 for Kalmyk to 0.036 for Russian Black-and-White and were comparable to those of the European breeds. The nine studied Russian breeds exhibited taurine ancestry along the C1 axis of the multidimensional scaling (MDS)-plot, but Yakut was clearly separated from the European taurine breeds on the C2 axis. Neighbor-Net and admixture analyses, discriminated three groups among the studied Russian breeds. Yakut and Kalmyk were assigned to a separate group because of their Turano-Mongolian origin. Russian Black-and-White, Kostromsky and Suksun showed transboundary European ancestry, which originated from the Holstein, Brown Swiss, and Danish Red breeds, respectively. The lowest level of introgression of transboundary breeds was recorded for the Kholmogor, Yaroslavl, Red Gorbatov and Bestuzhev breeds, which can be considered as an authentic genetic resource. Conclusions Whole-genome SNP analysis revealed that Russian native and locally developed breeds have conserved authentic genetic patterns in spite of the considerable influence of Eurasian taurine cattle. In this paper, we provide fundamental genomic information that will contribute to the development of more accurate breed conservation programs and genetic improvement strategies.

Abstract Background Russia has a diverse variety of native and locally developed sheep breeds with coarse, fine, and semi-fine wool, which inhabit different climate zones and landscapes that range from hot deserts to harsh northern areas. To date, no genome-wide information has been used to investigate the history and genetic characteristics of the extant local Russian sheep populations. To infer the population structure and genome-wide diversity of Russian sheep, 25 local breeds were genotyped with the OvineSNP50 BeadChip. Furthermore, to evaluate admixture contributions from foreign breeds in Russian sheep, a set of 58 worldwide breeds from publicly available genotypes was added to our data. Results We recorded similar observed heterozygosity (0.354–0.395) and allelic richness (1.890–1.955) levels across the analyzed breeds and they are comparable with those observed in the worldwide breeds. Recent effective population sizes estimated from linkage disequilibrium five generations ago ranged from 65 to 543. Multi-dimensional scaling, admixture, and neighbor-net analyses consistently identified a two-step subdivision of the Russian local sheep breeds. A first split clustered the Russian sheep populations according to their wool type (fine wool, semi-fine wool and coarse wool). The Dagestan Mountain and Baikal fine-fleeced breeds differ from the other Merino-derived local breeds. The semi-fine wool cluster combined a breed of Romanian origin, Tsigai, with its derivative Altai Mountain, the two Romney-introgressed breeds Kuibyshev and North Caucasian, and the Lincoln-introgressed Russian longhaired breed. The coarse-wool group comprised the Nordic short-tailed Romanov, the long-fat-tailed outlier Kuchugur and two clusters of fat-tailed sheep: the Caucasian Mountain breeds and the Buubei, Karakul, Edilbai, Kalmyk and Tuva breeds. The Russian fat-tailed breeds shared co-ancestry with sheep from China and Southwestern Asia (Iran). Conclusions In this study, we derived the genetic characteristics of the major Russian local sheep breeds, which are moderately diverse and have a strong population structure. Pooling our data with a worldwide genotyping set gave deeper insight into the history and origin of the Russian sheep populations.

Comparison of genomic footprints in chicken breeds with different selection history is a powerful tool in elucidating genomic regions that have been targeted by recent and more ancient selection. In the present work, we aimed at examining and comparing the trajectories of artificial selection in the genomes of the native egg-type Russian White (RW) and meat-type White Cornish (WC) breeds. Combining three different statistics (top 0.1% SNP by FST value at pairwise breed comparison, hapFLK analysis, and identification of ROH island shared by more than 50% of individuals), we detected 45 genomic regions under putative selection including 11 selective sweep regions, which were detected by at least two different methods. Four of such regions were breed-specific for each of RW breed (on GGA1, GGA5, GGA8, and GGA9) and WC breed (on GGA1, GGA5, GGA8, and GGA28), while three remaining regions on GGA2 (two sweeps) and GGA3 were common for both breeds. Most of identified genomic regions overlapped with known QTLs and/or candidate genes including those for body temperatures, egg productivity, and feed intake in RW chickens and those for growth, meat and carcass traits, and feed efficiency in WC chickens. These findings were concordant with the breed origin and history of their artificial selection. We determined a set of 188 prioritized candidate genes retrieved from the 11 overlapped regions of putative selection and reviewed their functions relative to phenotypic traits of interest in the two breeds. One of the RW-specific sweep regions harbored the known domestication gene, TSHR. Gene ontology and functional annotation analysis provided additional insight into a functional coherence of genes in the sweep regions. We also showed a greater candidate gene richness on microchromosomes relative to macrochromosomes in these genomic areas. Our results on the selection history of RW and WC chickens and their key candidate genes under selection serve as a profound information for further conservation of their genomic diversity and efficient breeding.

Caucasian tur (Capra caucasica) is native to Greater Caucasus Mountain Chain from Azerbaijan and Georgia in the East to Krasnodar region of Russia in the West. This species is divided into two subspecies (by some authors into species)—East-Caucasian tur and West-Caucasian tur and a subpopulation referred to as Mid-Caucasian tur. Up to date most of the genetic studies of Caucasian tur are based on mitochondrial DNA sequences and comprehensive investigation based on nuclear DNA is required for clarification of its genetic diversity and population structure. In our work, we assessed the applicability of Illumina Goat SNP50 BeadChip for genetic studies of Caucasian tur. Total of 15 specimens of Capra caucasica including East-Caucasian tur from Dagestan (E_TUR, n = 5), West-Caucasian tur from Karachay-Cherkessia (W_TUR, n = 5), and Mid-Caucasian tur from Kabardino-Balkaria (M_TUR, n = 5) were genotyped. After quality control, 5544 polymorphic loci, which were distributed all over 29 autosomes, were detected. The lowest number of SNPs was found on the 25th chromosome—68, and the highest on the 1st chromosome—348. It was shown that all the three groups of Caucasian tur clustered separately. A total of 2061 SNPs were common for all the populations, 594 were found only in W_TUR, 689 in E_TUR, and 530 in M_TUR. Individual heterozygosity ranged from 0.273 to 0.282 in W_TUR, from 0.217 to 0.253 in E_TUR, and from 0.255 to 0.283 in M_TUR. A clinal pattern of genetic variation was revealed. It was suggested to consider Caucasian tur a single species with several ecotypes. Thus, in our study we demonstrated that the Illumina Goat SNP50 BeadChip developed for domestic goats can be used as a useful tool for genetic studies of Caucasian tur.

The comparative molecular genetic study of museum and modern representatives of cattle breeds can help to elucidate the origin and maintenance of historical genetic components in modern populations. We generated the consensus genotypes for 11 microsatellite loci for 24 museum samples of Kalmyk, Kyrgyz, and Kazakh cattle, dated from the first quarter of the 20th century, and compared them with those of modern Kalmyk, Kyrgyz, and Kazakh white-headed breeds. The level of genetic diversity of the modern Kalmyk and Kyrgyz cattle (uHe = 0.771–0.778) was similar to those observed in the museum samples (uHe = 0.772–0.776), while a visible decrease in genetic variability in the modern Kazakh white-headed breed compared to museum Kazakh cattle was detected (uHe = 0.726 and 0.767, respectively). The PCA plot, FST- and Jost’s D-based networks, and STRUCTURE clustering provided strong evidence of the maintenance of the historical genetic background in modern populations of Kalmyk and Kyrgyz cattle. In spite of the allele pool of Kazakh white-headed cattle having undergone great changes compared to the museum Kazakh cattle, several animals still carry the visible aspect of the historical genetic components. Our results can be used for the selection of individuals for the creation of gene banks and may significantly improve the efficiency of conservation programs aimed at preserving genetic diversity in the national genetic resources of cattle.

Karachaev goat (Capra hircus) is a local breed from North-Caucasus region, Russia. Here we present complete mitochondrial genome of Karachaev goat from the republic of Karachaevo-Cherkessia, Russia. The length of the studied sequence was 16,624 bp in size. It was shown that the studied specimen belonged to haplogroup A.

Sheep farming has been an important sector of the UK’s economy and rural life for many centuries. It is the favored source of wool, meat and milk products. In the era of exponential progress in genomic technologies, we can now address the questions of what is special about UK sheep breed genotypes and how they differ genetically form one another and from other countries. We can reflect how their natural history has been determined at the level of their genetic code and what traces have been left in their genomes because of selection for phenotypic traits. These include adaptability to certain environmental conditions and management, as well as resistance to disease. Application of these advancements in genetics and genomics to study sheep breeds of British domestic selection has begun and will continue in order to facilitate conservation solutions and production improvement.

Currently, the intraspecific taxonomy of snow sheep (Ovis nivicola) is controversial and needs to be specified using DNA molecular genetic markers. In our previous work using whole-genome single nucleotide polymorphism (SNP) analysis, we found that the population inhabiting Kharaulakh Ridge was genetically different from the other populations of Yakut subspecies to which it was usually referred. Here, our study was aimed at the clarification of taxonomic status of Kharaulakh snow sheep using mitochondrial cytochrome b gene. A total of 87 specimens from five different geographic locations of Yakut snow sheep as well as 20 specimens of other recognized subspecies were included in this study. We identified 19 haplotypes, two of which belonged to the population from Kharaulakh Ridge. Median-joining network and Bayesian tree analyses revealed that Kharaulakh population clustered separately from all the other Yakut snow sheep. The divergence time between Kharaulakh population and Yakut snow sheep was estimated as 0.48 ± 0.19 MYA. Thus, the study of the mtDNA cytb sequences confirmed the results of genome-wide SNP analysis. Taking into account the high degree of divergence of Kharaulakh snow sheep from other groups, identified by both nuclear and mitochondrial DNA markers, we propose to classify the Kharaulakh population as a separate subspecies.

Three semi-domesticated reindeer breeds inhabiting the Republic of Sakha – Yakutia have been characterized using nine microsatellite markers. Genomic DNA was isolated from tissue samples of 123 individuals of the Chukotka (Khargin) (CHU, n=47), the Evenk (EVK, n=32) and the Even (EVN, n=44) breeds, collected from different regions of Yakutia. Fragment analysis and sizing were run on ABI 3131xl genetic analyzer. Allele frequencies were calculated and used for the characterization of reindeer breeds and the evaluation of their genetic biodiversity. Nei’s standard genetic distance was calculated and used for the construction of a neighbor-joining tree. Statistical analysis was conducted with GenAIEx 6.5.1, PAST2.15 and STRUCTURE2.3.4 software. The highest number of alleles, such as informative (with a frequency more than 5%), effective (Ne) and private (Pr), was detected in the CHU breed: Na≥5%=5.333±0.441, Ne=4.517±0393 and Pr =1.111±0.389, while the EVN breed had the lowest number: 4.778±0.324, 4.315±0.488 and 0.444±0.242, respectively. The EVN breed occupied an intermediate position (5.000±0.373, 4.408±0.315 and 0.889±0.261). Among reindeer breeds, observed heterozygosity ranged from 0.729 ± 0.026 to 0.608±0.050 with the lowest value found in CHU reindeer and the highest in EVK reindeer. A heterozygotes’ deficiency was observed in all reindeer breeds. At K=3, STRUCTURE analysis matches with the data of Nei's genetic distance dimension results, indicating the presence of a common consistent pattern. CHU and EVK reindeer breeds are characterized by a closer genetic relationship in comparison with the EVN breed, which formed a separate cluster.

Reindeer (Rangifer tarandus L. 1758) are an essential element of the Russian Far North, providing a significant source of nutrition for the representatives of 18 ethnicities. The species has wild and domestic forms, which are in constant interaction. The aim of our study was to characterize the genetic structure of domestic and wild reindeer populations, using a genome-wide bovine genotyping array (BovineHD BeadChip). The wild reindeer samples were obtained from the western Taymyr Peninsula population and from the taiga and tundra populations in the Sakha Republic (Yakutia). The domestic populations included the Evenk, Even, and Chukotka-Khargin breeds of Yakutia and the Nenets breed from the Nenets Autonomous district and Murmansk region. The level of genetic diversity was higher for the wild population. Analyzing Neighbor-Net tree, multidimensional scaling, and Structure results, we observed strong genetic population structure and clear differentiation between domestic and wild populations. All regional populations of domestic reindeer were clearly separated, while wild reindeer showed similar genetic backgrounds. Nevertheless, we found contrasting patterns in the genetic structure of the tundra and taiga reindeer, in accordance with their morphological and ecological differences. Thus, our study revealed a clear genetic differentiation between domestic and wild reindeer populations. It provides novel insights into the genetic diversity and structure of reindeer populations, to support resource utilization and aid in the development of genetic improvement strategies and conservation programs for this species.

General senescence of the adult organism is closely connected with reproductive one. Meanwhile, the age-related reduction in the female fertility is primarily associated with a decline in the gamete quality. Molecular and cellular changes in oocytes of old mammalian females are very similar to those occurring during aging of matured ova of their young counterparts, suggesting similarities in underlying mechanisms. The aim of the present work was to study actions of two related pituitary hormones, prolactin (PRL) and growth hormone (GH), on age-associated modifications of metaphase-II (M-II) chromosomes in bovine oocytes using a model of the prolonged culture. We analyzed: (1) effects of PRL and GH on abnormal changes in the chromosome morphology in aging matured oocytes and the role of cumulus cells in these effects and (2) signaling pathways involved in the hormone actions. During the prolonged culture of oocytes, a gradual rise in the frequency of destructive modifications of M-II chromosomes was revealed. In the case of cumulus-enclosed oocytes (CEOs), PRL and GH exerted dose-dependent biphasic effects on the frequency of these modifications. Both PRL (50 ng/ml) and GH (10 ng/ml) decelerated the abnormal chromosome changes in CEOs, but did not affect the chromosome configuration in denuded oocytes. Concurrently, the presence of PRL and GH receptors in cumulus cells surrounding matured oocytes was demonstrated. Attenuating effects of both hormones on the chromosome modifications in aging CEOs were abolished by PP2 (an inhibitor of Src-family tyrosine kinases), triciribine (an inhibitor of Akt kinase), and calphostin C (a protein kinase C inhibitor). Our findings indicate that PRL and GH can exert the similar decelerating action on age-associated alterations in the M-II chromosome morphology in bovine ova, which is mediated by cumulus cells and may be related to activation of Src-family tyrosine kinases as well as Akt- and protein kinase C-dependent signal pathways.