Mitochondrial Diversity and Phylogenetic Relationship of Eight Native Bulgarian Sheep Breeds.

The geographical, geomorphological, and climatic characteristics of Bulgaria are particularly favorable for animal breeding and, above all, for pastoral farming and sheep breeding. These conditions created prerequisites for the creation of about 30 unique local breeds of sheep. In this study we investigated the genetic diversity of eight of the most popular Bulgarian native breeds, based on the sequence analysis of a part of the mitochondrial D-loop region. An almost entire mitochondrial DNA (mtDNA) D-loop region (1180 bp) was amplified and sequenced. The obtained results showed the presence of a large number of haplotypes-225, belonging to two main haplogroups. The majority of samples showed a high prevalence of the European haplogroup B (95.2%) while the remaining individuals were assigned to haplogroup A (4.8%). None of the other reported mitochondrial haplogroups were observed. The number of polymorphic sites, nucleotide and haplotype diversity was high (240, 0.01237, and 0.9968, respectively), which is evidence for multiple maternal origins in all populations. The Tajima D-test value in all the study populations was -1.905 (p < 0.05), indicating that the abundance of rare alleles was most likely due to population expansion after a recent bottleneck. The Median joining network showed that almost all haplotypes belonging to haplogroup B formed a star-like network, which revealed a weak genetic differentiation and a large gene flow between the Bulgarian native breeds.

Composition and diversity of bacterial communities associated with honey bee foragers from two contrasting environments.

The honey bee is associated with a diverse community of microbes (viruses, bacteria, fungi, and protists), commonly known as the microbiome. Here, we present data on honey bee microbiota from two localities having different surrounding landscapes - mountain (the Rhodope Mountains) and lowland (the Danube plain). The bacterial communities of abdomen of adult bees were studied using amplicon sequencing of the 16S rRNA gene. The composition and dominance structure and their variability within and between localities, alpha and beta diversity, and core and differential taxa were compared at different hierarchical levels (operational taxonomic units to phylum). Seven genera (Lactobacillus, Gilliamella, Bifidobacterium, Commensalibacter, Bartonella, Snodgrassella, and Frischella), known to include core gut-associated phylotypes or species clusters, dominated (92-100%) the bacterial assemblages. Significant variations were found in taxa distribution across both geographical regions and within each apiary. Lactobacillus (Firmicutes) prevailed significantly in the mountain locality followed by Gilliamella and Bartonella (Proteobacteria). Bacteria of four genera, core (Bartonella and Lactobacillus) and non-core (Pseudomonas and Morganella), dominated the bee-associated assemblages of the Danube plain locality. Several ubiquitous bacterial genera (e.g., Klebsiella, Serratia, and Providencia), some species known also as potential and opportunistic bee pathogens, had been found in the lowland locality. Beta diversity analyses confirmed the observed differences in the bacterial communities from both localities. The occurrence of non-core taxa contributes substantially to higher microbial richness and diversity in bees from the Danube plain locality. We assume that the observed differences in the microbiota of honey bees from both apiaries are due to a combination of factors specific for each region. The surrounding landscape features of both localities and related vegetation, anthropogenic impact and land use intensity, the beekeeping management practices, and bee health status might all contribute to observed differences in bee microbiota traits.

Molecular Detection and Phylogenetic Analysis of Deformed Wing Virus and Sacbrood Virus Isolated from Pollen.

Among many pathogens and pests, honey bee viruses are known as one of the most common cause of diseases in honey bee colonies. In this study, we demonstrate that pollen grains and bee bread are potential sources of viral DNA. We extracted DNA from 3 types of pollen samples: directly provided by beekeepers (n = 12), purchased from trade markets (n = 5), and obtained from honeycombs (bee bread, n = 10). The extracted DNA was used for molecular detection (RT-PCR analysis) of six of the most widely distributed honey bee viruses: deformed wing virus, sacbrood virus, acute bee paralysis virus, black queen cell virus, Kashmir bee virus, Israeli acute paralysis virus, and chronic bee paralysis virus. We successfully managed to establish only the deformed wing virus (DWV) and the sacbrood virus (SBV), with different distribution frequencies depending on the territory of the country. The phylogenetic analyses of Bulgarian isolates were performed with the most similar sequences available in molecular databases from other countries. Phylogenies of Bulgarian viral strains demonstrated genetically heterogeneous populations of DWV and relatively homogenous populations of SBV. In conclusion, the results obtained from the current study have shown that pollen is a valuable source for molecular detection of honey bee pathogens. This allows epidemiological monitoring of honey bee diseases at a regional and a national level.

Genetic Diversity and Structure of the Main Danubian Horse Paternal Genealogical Lineages Based on Microsatellite Genotyping.

The Danubian horse, together with the Pleven and the Eastern Bulgarian horse breeds, is one of the modern breeds in Bulgaria. The objective of this study was to compare the genetic structure and genetic diversity of six paternal genealogical lineages of the Danubian horse breed (Zdravko, NONIUS XVII-30, Torpedo, Lider, Kalifa, and Hrabar). In total, 166 individuals from the six genealogical lines were investigated, based on 15 STR markers (short tandem repeats, also known as microsatellites). In total, 184 alleles were found in the six populations, using 15 microsatellite loci. The mean number of alleles, the effective number of alleles, and the polymorphic information content (PIC) values per locus were 12.28, 9.48, and 0.73, respectively. In a comparison of the allelic diversity among sire lineages, the highest genetic diversity (Na) was observed in Lider and Kalifa (14.60 ± 0.21), while the lowest value of this parameter was observed in the Zdravko lineage 4.20 ± 0.35. The largest genetic diversity was found in loci HMS3 and HMS7, with 13 alleles, and the smallest polymorphism was noted for the locus ASB17, with 10 alleles. The level of observed heterozygosity was in the range of 0.65 ± 0.069 for the Zdravko lineage to 0.93 ± 0.01 for the Torpedo lineage. The expected heterozygosity level range was from 0.57 ± 0.048 to 0.91 ± 0.01 for all horse lineages. Structure analysis revealed three main gene pools in the study population. The first pool included the Zdravko lineage; the second had the NONIUS XVII-30, Torpedo, Lider, and Kalifa lineages; and the third defined the Hrabar lineage, which was significantly differentiated from the other genealogical lineages.

The Sequence Analysis of Mitochondrial DNA Revealed Some Major Centers of Horse Domestications: The Archaeologist's Cut.

The question about the time and the place of horse domestication, a process which had a profound impact on the progress of mankind, is disputable. According to the most widely accepted hypothesis, the earliest domestication of the horse happened in the western parts of the Eurasian steppes, between the Northern Black Sea region and present-day Kazakhstan and Turkmenistan. It seems that it occurred not earlier than the first half and most probably during the middle (even the last third) of the fourth millennium BC (from ∼ 5.5 kya). The next steps of large-scale horse breeding occurred almost simultaneously in Eurasia and North Africa due to the development of the social structure of human communities. On the other hand, the morphological differences between wild and domestic animals are rather vague and the genetic introgression between them is speculative. In this review, we have tried to gather all available scientific data on the existing possible hypotheses for the earliest domestication of the horse, as well as to highlight some data on the most plausible ones. This is due to the frequency of some significant data on the frequency of strictly defined mitotypes in different historical periods of human civilizations existing in the same periods.

Sequence analysis of the mitochondrial D-loop region throws a new light on the origin of Hungarian Nonius, Danubian Horse and Serbian Nonius.

The objective of our study was to investigate the genetic structure of yet uninvestigated populations of three closely related horse breeds - the Danubian Horse, the Hungarian Nonius and the Serbian Nonius - in order to clarify their origin and genetic diversity. A 640-bp-long fragment of the mtDNA D-loop region was amplified and sequenced. The results showed that the investigated breeds have different genetic profiles although they share some common characteristics. We identified nine of the 17 haplogroups described in modern horses. Most of the obtained sequences fall into the M, L, G, and O'P lineages, which is indicative of the genetic profile of the ancestral mares that had probably been used at the initial stages of the formation of the breeds. The population of the Danubian Horse is characterised by a high prevalence of the Anatolian specific haplogroup G (45%), followed by the Western Eurasian specific haplogroups L and M (both about 21%). In the Hungarian Nonius breed we found the highest frequency of the Western Eurasian haplogroup M (44%), followed by the Middle Eastern O'P (26%) and the Central Asian specific E (13%) and G (13%). The Serbian Nonius showed a distinct genetic profile, characterised by a high prevalence of the rare European haplogroup D (67%), followed by the Central Asian specific haplogroup G (17%). The high percentage of haplogroups shared especially between the Danubian and the Hungarian Nonius indicates the possibility of a common origin of the two breeds. In contrast, the Serbian Nonius showed a specific genetic profile, which can be explained by a different and independent origin.

New data on the evolutionary history of the European bison (Bison bonasus) based on subfossil remains from Southeastern Europe.

The origin and evolutionary history of the European bison Bison bonasus (wisent) have become clearer after several morphological, genomic, and paleogenomic studies in the last few years, but these paleogenomic studies have raised new questions about the evolution of the species. Here, we present additional information about the population diversity of the species based on the analysis of new subfossil Holocene remains from the Balkan Peninsula. Seven ancient samples excavated from caves in Western Stara Planina in Bulgaria were investigated by mitochondrial D-loop (HVR1) sequence analysis. The samples were dated to 3,800 years BP by radiocarbon analysis. Additionally, a phylogenetic analysis was performed to investigate the genetic relationship among the investigated samples and all mitochondrial DNA sequences from the genus Bison available in GenBank. The results clustered with the sequences from the extinct Holocene South-Eastern (Balkan) wisent to the fossil Alpine population from France, Austria, and Switzerland, but not with those from the recent Central European (North Sea) one and the now extinct Caucasian population. In conclusion, these data indicate that the Balkan wisent that existed in historical time represented a relict and probably an isolated population of the Late Pleistocene-Holocene South-Western mountainous population of the wisent.

Molecular Detection of Nosema spp. in Honey in Bulgaria.

Environmental DNA (eDNA) analysis is related to screening genetic material of various organisms in environmental samples. Honey represents a natural source of exogenous DNA, which allows for the detection of different honey bee pathogens and parasites. In the present study, we extracted DNA from 20 honey samples from different regions in Bulgaria and tested for the presence of DNA of the ectoparasitic mite Varroa destructor, as well as Nosema apis and Nosema ceranae. Only Nosema ceranae was detected, showing up in 30% of all samples, which confirms the widespread prevalence of this pathogen. All positive samples were found in plain regions of the country, while this pathogen was not detected in mountainous parts. None of the samples gave positive amplifications for the Nosema apis and Varroa mite. The obtained results from this study confirm previous observations that eDNA contained in honey is a potent source for effective biomonitoring of actual diseases in the honey bee.

Mitochondrial Profiles of the East Bulgarian and the Pleven Horse Breeds.

It is well known that horse breeding in Bulgaria is a cultural heritage in Bulgaria, dating from prehistoric and historic times. Until now, molecular data on Bulgarian horses from the plain regions of the country were not available. Therefore, for the first time, we have collected genetic information about some modern horse breeds from the plain regions in Bulgaria. A total of 50 horses originating from different families from two different breeds were investigated: the first one was the Pleven horse (n = 11, breeding in the Danubian Plain), and the second one was the East Bulgarian horse breed (n = 39, Shumen district). These breeds were genotyped according to the mitochondrial D-loop region. The results showed that the Pleven horse particularly carries the European haplogroup L (45.45%), followed by the Middle East haplogroup C (27.27%). In contrast to the Pleven horse, the East Bulgarian horse breed revealed almost equal frequencies of the European specific haplogroup L (33.33%) and the Central Asiatic haplogroup Q (35.90%). Analyses of these two horse breeds revealed a specific genetic profile, but it is obvious that the East Bulgarian horse showed an unusual, mixed profile-a massive admixture with the Asiatic-type haplogroup Q and a high haplogroup diversity. In conclusion, the differences in genetic structure of the two plain horse breeds may be explained with the various horse breeds involved in their formation.

Evidence for Early European Neolithic Dog Dispersal: New Data on Southeastern European Subfossil Dogs from the Prehistoric and Antiquity Ages.

The history of dog domestication is still under debate, but it is doubtless the process of an ancient partnership between dogs (Canis familiaris) and humans. Although data on ancient DNA for dog diversity are still incomplete, it is clear that several regional dog populations had formed in Eurasia up to the Holocene. During the Neolithic Revolution and the transition from hunter-gatherer to farmer societies, followed by civilization changes in the Antiquity period, the dog population structure also changed. This process was due to replacement with newly formed dog populations. In this study, we present for the first time mitochondrial data of ancient dog remains from the Early Neolithic (8000 years before present (BP)) to Late Antiquity (up to 3th century AD) from southeastern Europe (the Balkans). A total of 16 samples were analyzed, using the mitochondrial D-loop region (HVR1). The results show the presence of A (70%) and B (25%) clades throughout the Early and Late Neolithic Period. In order to clarify the position of our results within the ancient dog population in Eneolithic Eurasia, we performed phylogenetic analysis with the available genetic data sets. This data showed a similarity of the ancient Bulgarian dogs to Italian (A, B, and C clades) and Iberian (clades A and C) dogs' populations. A clear border can be seen between southern European genetic dog structure, on the one hand, and on the other hand, central-western (clade C), eastern (clade D) and northern Europe (clades A and C). This corresponds to genetic data for European humans during the same period, without admixture between dog populations. Also, our data have shown the presence of clade B in ancient Eurasia. This is not unexpected, as the B haplogroup is widely distributed in extant Balkan dogs and wolves. The presence of this clade both in dogs and in wolves on the Balkans may be explained with hybridization events before the Neolithic period. The spreading of this clade across Europe, together with the A clade, is related to the possible dissemination of newly formed dog breeds from Ancient Greece, Thrace, and the Roman Empire.

Molecular and morphological characterisation of Longidoruspolyae sp. n. and L.pisi Edward, Misra & Singh, 1964 (Dorylaimida, Longidoridae) from Bulgaria.

Longidoruspolyae sp. n., a bisexual nematode species found in the rhizosphere of pear tree (Pyruscommunis L.), is described and characterised using an integrative approach. The new species has a female body length of 6.8-9.1 mm; a comparatively long odontostyle (114.0-127.5 µm); a narrow lip region (14.0-15.5 µm), anteriorly flattened and almost continuous with the body profile; pocket-like amphidial pouches long, deeply bilobed, and slightly asymmetrical, a guide ring at 37-42 µm from the anterior end; normal arrangement of pharyngeal glands; and a short bluntly rounded to hemispherical tail. Four juvenile stages identified: the first stage with a digitate tail, and the second and subsequent stages with a bluntly rounded tail. Males have one adcloacal pair and a row of 10 or 11 single ventromedian supplements; spicules 71.0-74.5 µm long. Based on morphometric data, the new species belongs to a group of species spread over Europe (L.arthensis, L.silvae, L.uroshis,), Iran (L.kheirii), and Syria (L.pauli), which share common characters such as amphidial fovea, lip region and tail shapes, similar odontostyle and body length, and similar first-stage juvenile tail shape. Codes for identifying the new species are A5, B2, C34, D3, E3, F45, G12, H1, I2, J1, K7. The phylogenetic analysis based on D2-D3 expansion domains of the rRNA gene revealed that the new species has the closest relationships with L.athesinus from Italy and three unidentified Longidorus spp. from USA (Longidorus sp. 1, Longidorus sp. 2, and Longidorus sp. 6). New morphometric and molecular data (18S rRNA gene, ITS1-5.8S-ITS2 regions and D2-D3 28S rRNA gene sequences) for three populations of L.pisi from Bulgaria were obtained and variations between populations are discussed.

Phylogeny of hymenolepidid cestodes (Cestoda: Cyclophyllidea) from mammalian hosts based on partial 28S rDNA, with focus on parasites from shrews.

The aims of the study are to enrich the partial 28S rDNA dataset for hymenolepidids by adding new sequences for species parasitic in the genera Sorex, Neomys and Crocidura (Soricidae) and to propose a new hypothesis for the relationships among mammalian hymenolepidids. New sequences were obtained for Coronacanthus integrus, C. magnihamatus, C. omissus, C. vassilevi, Ditestolepis diaphana, Lineolepis scutigera, Spasskylepis ovaluteri, Staphylocystis tiara, S. furcata, S. uncinata, Vaucherilepis trichophorus and Neoskrjabinolepis sp. The phylogenetic analysis (based on 56 taxa) confirmed the major clades identified by Haukisalmi et al. (Zool Scr 39:631-641, 2010) based on analysis of 31 species: Ditestolepis clade, Hymenolepis clade, Rodentolepis clade and Arostrilepis clade; however, the support was weak for the early divergent lineages of the tree and for the Arostrilepis clade. Novelties revealed include the molecular evidence for the monophyly of Coronacanthus, the non-monophyletic status of Staphylocystis and the polyphyly of Staphylocystoides. The analysis has confirmed the monophyly of Hymenolepis, the monophyly of hymenolepidids from glirids, the position of Pararodentolepis and Nomadolepis as sister taxa, the polyphyly of Rodentolepis, the position of Neoskrjabinolepis and Lineolepis as sister taxa, and the close relationship among the genera with the entire reduction of rostellar apparatus. Resolved monophyletic groups are supported by the structure of the rostellar apparatus. The diversification of the Ditestolepis clade is associated with soricids. The composition of the other major clades suggests multiple evolutionary events of host switching, including between different host orders. The life cycles of Coronacanthus and Vaucherilepis are recognised as secondarily aquatic as these taxa are nested in terrestrial groups.

Mitochondrial diversity of Bulgarian native dogs suggests dual phylogenetic origin.

The dog has been the first domesticated animal to have a central role in human society from ancient times to present day. Although there have been numerous investigations of dog phylogeny and origin, genetic data of dogs in the region of the Balkan Peninsula (South-Eastern Europe) are still scarce. Therefore, the aim of the present study was to perform phylogenetic analysis of three native Bulgarian dog breeds. A total of 130 samples were analyzed at HVR1 (hypervariable region, D-loop region). The samples were taken from two hunting dog breeds (Bulgarian Hound Dog: Barak, n = 34; Bulgarian Scenthound Dog: Gonche, n = 45) as well as from a Bulgarian Shepherd Dog (n = 51). The first two breeds are reared in a flat region of the country (the Northern part of Bulgaria, the Danubian Plain), while the last breed is a typical representative of the mountainous part of the country. The results have shown the presence of almost all main clades-A, B, C and D-in the three dog breeds taken together, except clades E and F, as expected. With regard to haplogroups distribution, there are clear differences among investigated breeds. While hunting breeds exhibit a prevalence of clade C, the mountainous Shepherd dog shows presence of the D2 haplogroup but absence of the C clade. In conclusion, the present study has been the first to investigate the mitochondrial DNA diversity of native dog breeds in Bulgaria. The results have revealed a clear difference of haplogroups dissemination in native hunting and shepherd dogs, which suggests a dual independent phylogenetic origin, without hybridization events between these dogs.

Molecular detection and phylogenetic assessment of six honeybee viruses in Apis mellifera L. colonies in Bulgaria.

Honey bee colonies suffer from various pathogens, including honey bee viruses. About 24 viruses have been reported so far. However, six of them are considered to cause severe infection which inflicts heavy losses on beekeeping. The aim of this study was to investigate incidence of six honey bee viruses: deformed wing virus (DWV), acute bee paralysis virus (ABPV), chronic bee paralysis virus (CBPV), sacbrood virus (SBV), kashmir bee virus (KBV), and black queen cell virus (BQCV) by a reverse transcription polymerase chain reaction (RT-PCR). A total of 250 adult honey bee samples were obtained from 50 colonies from eight apiaries situated in three different parts of the country (South, North and West Bulgaria). The results showed the highest prevalence of DWV followed by SBV and ABPV, and one case of BQCV. A comparison with homology sequences available in GenBank was performed by phylogenetic analysis, and phylogenetic relationships were discussed in the context of newly described genotypes in the uninvestigated South Eastern region of Europe. In conclusion, the present study has been the first to provide sequencing data and phylogenetics analyses of some honey bee viruses in Bulgaria.

The first report of the prevalence of Nosema ceranae in Bulgaria.

Nosema apis and Nosema ceranae are the two main microsporidian parasites causing nosematosis in the honey bee Apis mellifera. The aim of the present study is to investigate the presence of Nosema apis and Nosema ceranae in the area of Bulgaria. The 16S (SSU) rDNA gene region was chosen for analysis. A duplex PCR assay was performed on 108 honey bee samples from three different parts of the country (South, North and West Bulgaria). The results showed that the samples from the northern part of the country were with the highest prevalence (77.2%) for Nosema ceranae while those from the mountainous parts (the Rodopa Mountains, South Bulgaria) were with the lowest rate (13.9%). Infection with Nosema apis alone and co-infection N. apis/N. ceranae were not detected in any samples. These findings suggest that Nosema ceranae is the dominant species in the Bulgarian honey bee. It is not known when the introduction of Nosema ceranae in Bulgaria has occurred, but as in the rest of the world, this species has become the dominant one in Bulgarian Apis mellifera. In conclusion, this is the first report for molecular detection of Nosema infection of honey bee in Bulgaria. The results showed that N. ceranae is the main Nosema species in Bulgaria.

Balkan brachicerous cattle - the first domesticated cattle in Europe.

The present study aimed to compare mitochondrial diversity among Balkan Neolithic/Chalcolithic cattle and present day Shorthorn Rhodopean cattle (Busha) to throw a new insight into European cattle domestication. The results showed that both ancient- and present-day samples belonged to the macrohaplogroup T. From the 28 sequences (8 ancient and 20 modern), the T1 and T2 haplogroup represent about 3.6% (1/28; 1/28). The T3 haplogroup was with the highest frequency - 57% (16/28). Based on the SNPs on 16057A and 16133C, the new T6 haplogroup was proposed. This haplogroup represents 75% from the ancient and 20% from the present day Bulgarian brachicerous cattle population. The survey in GenBank data base did not find a similar motif, except for the recent Serbian Busha cattle. Overall, these results showed that: (i) The newly named T6 haplogroup is Balkan specific; (ii) The T6 haplogroup survives in present day Bulgarian rhodopean cattle; (iii) The Balkan brachicerous cattle is the oldest European cattle breed.

Biogeography and phylogenetic position of Enchodeloides signyensis (Loof, 1975), gen. n., comb. n. from Maritime Antarctic (Nematoda, Nordiidae).

The taxonomic position of the endemic Antarctic species Enchodeloides signyensis (Loof, 1975), gen. n., comb. n. (= Enchodelus signyensis Loof, 1975) is discussed on the basis of morphological study, including SEM, morphometric data, postembryonic observations, and sequence data of 18S rDNA and the D2-D3 expansion fragments of the large subunit rDNA. A number of characters such as the cuticle and stoma structures, including the presence of moderately developed cuticularised ring around the oral aperture, peculiarities of pharynx expansion, size and position of the posterior pair of pharyngeal nuclei, a less complex uterus, and the position of a posterior ventromedian supplement show that this species differs substantially from the other members of the genus Enchodelus. Furthermore, both the 18S and 28S rDNA-based phylogenetic trees of the Enchodelus sequences available in the GenBank formed two distinct clusters with E. signyensis being a part of a well-supported group with species of the genus Pungentus; therefore, it is proposed that its taxonomic position should be reconsidered.

The morphological and molecular identity of Longidorus piceicola Lisková, Robbins & Brown, 1997 from Romania (Nematoda, Dorylaimida).

Longidorus piceicola, a new geographical and host record from Romania, was described and illustrated on the basis of two populations originating from a coniferous and a deciduous forest. The main morphological characters of specimens from Romania correspond very well with the type material collected from the soil around Picea abies L. (Slovakia) except for the shorter body and tail. The D2-D3 fragment of 28S rDNA from both populations was amplified and sequenced, and the sequences were identical to L. piceicola sequence from Slovakia. The partial 18S-ITS1-5.8S-ITS2 rDNA regions from one of the populations were sequenced for the first time. The evolutionary relationships between L. piceicola and the closest species L. intermedius based on D2-D3 sequence divergence and single-nucleotide polymorphisms are discussed. Although having very low sequence dissimilarity (0.3-0.9 %) both species have distinct morphology and biology. Longidorus piceicola differs from L. intermedius in having a much longer odontostyle, body, distance anterior end - guide ring, a wider lip region, more ventromedian supplements (11 vs 5-7) in the male, and develops through four rather than three juvenile stages. Furthermore, L. piceicola occurs more frequently in association with conifers, while L. intermedius is found mainly in oak forests.

Mitochondrial diversity in mountain horse population from the South-Eastern Europe.

It is believed that population structure of mountain horse breeds is preserved from any genetic introgression, because of their geographical isolation and specific semi-wild life style of husbandry. Till date there are no molecular data for the Balkan horses. In this study we try to give information about some autochthonous mountain horse breeds from Bulgaria. A total of 121 horses from three different mountain massifs are presented: Stara Planina (the Balkan mountain), the Rhodopes and Rila-Pirin massif were genotyped according to mitochondrial D-loop region. The results showed huge diversity of all known haplogroups with exception of C, F and R. West Eurasian haplogroups B, D, M and L were with the highest frequencies. Haplogroups A, J, I, O'P and Q were also observed with the highest frequencies, but not equally distributed among the three populations. Analyses of the horse breeds reveal preserved genetic profile of the Balkan and the Rhodopes mountains populations. In contrast, a Rila-Pirin breed unexpectedly showed mixed profile - a massive genetic introgression with an Asiatic-type haplogroups. A similar mixed Euro-Asiatic haplotype profile possessed the Carpathian mountain pony, although both populations are separated geographically and historically. The genetic pool of three Bulgarian mountain horse populations is highly heterogenic and because of that these breeds should be preserved.

An independent event of Neolithic cattle domestication on the South-eastern Balkans: evidence from prehistoric aurochs and cattle populations.

Neolithic/Chalcolithic livestock domestication is an important issue for understanding the mode of life and economics of ancient human communities. The Balkans appears to be a crucial point for clarifying the socio-economical interrelations between the Oldest Middle Eastern/Anatolian and newly formed cultures in Europe. Two main hypotheses regarding the early history of cattle domestication, from their ancestor - the aurochs, have been discussed: multipoint domestication centers or single point origin and subsequent worldwide dissemination. In this study, we provide molecular data about the Balkan aurochs for the first time as well as additional information for the Neolithic/Chalcolithic cattle populations in this geographic location. A total of seventeen samples from different ancient settlements were analyzed according to D-loop control region. The results did not show different genetic profile of wild and domestic populations. All haplotypes were found to belong to the basic macro-haplogroup T. The majority of specimens (n = 14) were defined to form a new Balkan-specific T6 haplogroup. Only two of the ancient samples analyzed were assigned to the T3 haplotype predominating in Europe. We attempt to throw new light on the earliest cattle domestication events in Europe, thus, the results presented are discussed in two directions: (a) The possibility of local independent domestication processes in Neolithic South-Eastern Europe; (b) The single point domestication in the Middle East and subsequent cattle dissemination in Europe. Our data does not exclude the possibility for independent domestication events followed by a second wave of parallel dissemination of cattle herds via the Mediterranean route.