DomeTree: a canonical toolkit for mitochondrial DNA analyses in domesticated animals

Quantifying the effects of the mitochondrial genome on milk production traits in dairy cows: empirical results and modelling challenges

Significant advances in livestock traits have been achieved primarily through selection strategies targeting variation in the nuclear genome, with little attention given to mitogenome variation. We analyzed the influence of the mitogenome on milk production traits of Holstein cattle in Croatia based on strategically generated next-generation sequencing data for 109 cows pedigree-linked to 7115 milk production records (milk, fat and protein yield) from 3006 cows (first 5 lactations). Since little is known about the biology of the relationship between mitogenome variation and production traits, our quantitative genetic modeling was complex. Thus, the proportion of total variance explained by mitogenome inheritance was estimated using 5 different models: (1) cytoplasmic model with maternal lineages (CYTO), (2) haplotypic model with mitogenome sequences (HAPLO), (3) amino acid model with unique amino acid sequences (AMINO), (4) evolutionary model based on a phylogenetic analysis using Bayesian Evolutionary Analysis Sampling Trees phylogenetic analysis (EVOL), and (5) mitogenome SNP model (SNPmt). The polygenic autosomal and X chromosome additive genetic effects based on pedigree were modeled, together with the effects of herd-year-season interaction, permanent environment, location, and age at first calving. The estimated proportions of phenotypic variance explained by mitogenome in 4 different models (CYTO, HAPLO, AMINO, and SNPmt) were found to be substantial given the size of mitogenome, ranging from 5% to 7% for all 3 milk traits. At the same time, a negligible proportion of the phenotypic variance was explained by mitogenome with the EVOL model. Similarly, in all models, no proportion of phenotypic variance was explained by the X chromosome. Although our results should be confirmed in other dairy cattle populations, including a large number of sequenced mitogenomes and nuclear genomes, the potential of utilizing mitogenome information in animal breeding is promising, especially as the acquisition of complete genome sequences becomes cost-effective.

Genomic Reconstruction of the Successful Establishment of a Feralized Bovine Population on the Subantarctic Island of Amsterdam

The feral cattle of the subantarctic island of Amsterdam provide an outstanding case study of a large mammalian population that was established by a handful of founders and thrived within a few generations in a seemingly inhospitable environment. Here, we investigated the genetic history and composition of this population using genotyping and sequencing data. Our inference showed an intense but brief founding bottleneck around the late 19th century and revealed contributions from European taurine and Indian Ocean Zebu in the founder ancestry. Comparative analysis of whole-genome sequences further revealed a moderate reduction in genetic diversity despite high levels of inbreeding. The brief and intense bottleneck was associated with high levels of drift, a flattening of the site frequency spectrum and a slight relaxation of purifying selection on mildly deleterious variants. Unlike some populations that have experienced prolonged reductions in effective population size, we did not observe any significant purging of highly deleterious variants. Interestingly, the population's success in the harsh environment can be attributed to preadaptation from their European taurine ancestry, suggesting no strong bioclimatic challenge, and also contradicting evidence for insular dwarfism. Genome scan for footprints of selection uncovered a majority of candidate genes related to nervous system function, likely reflecting rapid feralization driven by behavioral changes and complex social restructuring. The Amsterdam Island cattle offers valuable insights into rapid population establishment, feralization, and genetic adaptation in challenging environments. It also sheds light on the unique genetic legacies of feral populations, raising ethical questions according to conservation efforts.

Ancient Sheep Genomes Reveal Four Millennia of North European Short-Tailed Sheep in the Baltic Sea Region

Sheep are among the earliest domesticated livestock species, with a wide variety of breeds present today. However, it remains unclear how far back this diversity goes, with formal documentation only dating back a few centuries. North European short-tailed (NEST) breeds are often assumed to be among the oldest domestic sheep populations, even thought to represent relicts of the earliest sheep expansions during the Neolithic period reaching Scandinavia <6,000 years ago. This study sequenced the genomes (up to 11.6X) of five sheep remains from the Baltic islands of Gotland and Åland, dating from the Late Neolithic (∼4,100 cal BP) to historical times (∼1,600 CE). Our findings indicate that these ancient sheep largely possessed the genetic characteristics of modern NEST breeds, suggesting a substantial degree of long-term continuity of this sheep type in the Baltic Sea region. Despite the wide temporal spread, population genetic analyses show high levels of affinity between the ancient genomes and they also exhibit relatively high genetic diversity when compared to modern NEST breeds, implying a loss of diversity in most breeds during the last centuries associated with breed formation and recent bottlenecks. Our results shed light on the development of breeds in Northern Europe specifically as well as the development of genetic diversity in sheep breeds, and their expansion from the domestication center in general.

Ancient Mitogenomes Reveal the Maternal Genetic History of East Asian Dogs

Recent studies have suggested that dogs were domesticated during the Last Glacial Maximum (LGM) in Siberia, which contrasts with previous proposed domestication centers (e.g. Europe, the Middle East, and East Asia). Ancient DNA provides a powerful resource for the study of mammalian evolution and has been widely used to understand the genetic history of domestic animals. To understand the maternal genetic history of East Asian dogs, we have made a complete mitogenome dataset of 120 East Asian canids from 38 archaeological sites, including 102 newly sequenced from 12.9 to 1 ka BP (1,000 years before present). The majority (112/119, 94.12%) belonged to haplogroup A, and half of these (55/112, 49.11%) belonged to sub-haplogroup A1b. Most existing mitochondrial haplogroups were present in ancient East Asian dogs. However, mitochondrial lineages in ancient northern dogs (northeastern Eurasia and northern East Asia) were deeper and older than those in southern East Asian dogs. Results suggests that East Asian dogs originated from northeastern Eurasian populations after the LGM, dispersing in two possible directions after domestication. Western Eurasian (Europe and the Middle East) dog maternal ancestries genetically influenced East Asian dogs from approximately 4 ka BP, dramatically increasing after 3 ka BP, and afterwards largely replaced most primary maternal lineages in northern East Asia. Additionally, at least three major mitogenome sub-haplogroups of haplogroup A (A1a, A1b, and A3) reveal at least two major dispersal waves onto the Qinghai-Tibet Plateau in ancient times, indicating eastern (A1b and A3) and western (A1a) Eurasian origins.

Phylogenetic Analysis of Russian Native Sheep Breeds Based on mtDNA Sequences

Eurasia is represented by all climatic zones and various environments. A unique breed variety of farm animals has been developed in Russia, whose territory covers a large area of the continent. A total of 69 local breeds and types of dairy, wool, and meat sheep (Ovis aries) are maintained here. However, the genetic diversity and maternal origin of these local breeds have not been comprehensively investigated. In this study, we describe the diversity and phylogeny of Russian sheep breeds inhabiting different geographical regions based on the analysis of complete sequences of mitochondrial genomes (mtDNA). Complete mtDNA sequences of the studied sheep were obtained using next-generation sequencing technology (NGS). All investigated geographical groups of sheep were characterized by high haplotype (Hd = 0.9992) and nucleotide diversity (π = 0.00378). Analysis of the AMOVA results showed that genetic diversity was majorly determined by within-population differences (77.87%). We identified 128 haplotypes in all studied sheep. Haplotypes belonged to the following haplogroups: B (64.8%), A (28.9%), C (5.5%), and D (0.8%). Haplogroup B was predominant in the western part of Russia. A high level of mtDNA polymorphism in the studied groups of local sheep indicates the presence of a significant reserve of unique genotypes in Russia, which is to be explored.

Mitochondrial DNA diversity of D-loop region in three native Turkish cattle breeds

This study aimed to reveal the genetic variability of mitochondrial DNA (mtDNA) displacement-loop (D-loop) region in 62 animals belonging to three native Turkish cattle breeds, namely Anatolian Black (AB), East Anatolian Red (EAR) and Zavot (ZAV), and to conduct phylogenetic relationship analyses to obtain deeper information on their genetic origin and breeding history by comparison of 6 taurine and 11 indicine breeds, together with 66 polymorphic sites, a total of 31 haplotypes, of which 15, 10 and 6 were detected in AB, EAR and ZAV, respectively. Mean nucleotide and haplotype diversity were 0.01 and 0.891, respectively, whereas the genetic differentiation derived from Wright's F ST index was 0.174 across the breeds. A significant level of total variation (17.42 %) was observed among breeds in molecular variance analysis. Six main haplogroups (T, T1, T2, T3, Q and I2) were detected in Anatolian cattle populations, where T3 was the most frequent among breeds (43.55 %), whereas I2, an indicine specific haplogroup, was observed only in ZAV. At the breed level, phylogenetic analyses supported by 198 sequences of 17 cattle breeds and 3 outgroup species retrieved from the GenBank clustered native Turkish cattle breeds with the taurine group rather than the indicine one, as expected. However, indicine admixture at low frequency (8.89 %) was detected in the ZAV breed for the first time due to more likely gene flow from indicine cattle breeds raised in neighbour countries, particularly Iran. This finding should be further investigated in all native Turkish and indicine cattle breeds from nearby countries to clarify gene flow and indicine admixture in Anatolian cattle.

Unveiling new perspective of phylogeography, genetic diversity, and population dynamics of Southeast Asian and Pacific chickens

The complex geographic and temporal origins of chicken domestication have attracted wide interest in molecular phylogeny and phylogeographic studies as they continue to be debated up to this day. In particular, the population dynamics and lineage-specific divergence time estimates of chickens in Southeast Asia (SEA) and the Pacific region are not well studied. Here, we analyzed 519 complete mitochondrial DNA control region sequences and identified 133 haplotypes with 70 variable sites. We documented 82.7% geographically unique haplotypes distributed across major haplogroups except for haplogroup C, suggesting high polymorphism among studied individuals. Mainland SEA (MSEA) chickens have higher overall genetic diversity than island SEA (ISEA) chickens. Phylogenetic trees and median-joining network revealed evidence of a new divergent matrilineage (i.e., haplogroup V) as a sister-clade of haplogroup C. The maximum clade credibility tree estimated the earlier coalescence age of ancestral D-lineage (i.e., sub-haplogroup D2) of continental chickens (3.7 kya; 95% HPD 1985-4835 years) while island populations diverged later at 2.1 kya (95% HPD 1467-2815 years). This evidence of earlier coalescence age of haplogroup D ancestral matriline exemplified dispersal patterns to the ISEA, and thereafter the island clade diversified as a distinct group.

Recovery of mitogenomes from whole genome sequences to infer maternal diversity in 1883 modern taurine and indicine cattle

Maternal diversity based on a sub-region of mitochondrial genome or variants were commonly used to understand past demographic events in livestock. Additionally, there is growing evidence of direct association of mitochondrial genetic variants with a range of phenotypes. Therefore, this study used complete bovine mitogenomes from a large sequence database to explore the full spectrum of maternal diversity. Mitogenome diversity was evaluated among 1883 animals representing 156 globally important cattle breeds. Overall, the mitogenomes were diverse: presenting 11 major haplogroups, expanding to 1309 unique haplotypes, with nucleotide diversity 0.011 and haplotype diversity 0.999. A small proportion of African taurine (3.5%) and indicine (1.3%) haplogroups were found among the European taurine breeds and composites. The haplogrouping was largely consistent with the population structure derived from alternate clustering methods (e.g. PCA and hierarchical clustering). Further, we present evidence confirming a new indicine subgroup (I1a, 64 animals) mainly consisting of breeds originating from China and characterised by two private mutations within the I1 haplogroup. The total genetic variation was attributed mainly to within-breed variance (96.9%). The accuracy of the imputation of missing genotypes was high (99.8%) except for the relatively rare heteroplasmic genotypes, suggesting the potential for trait association studies within a breed.

Genetic diversity and population genetic structure of Cambodian indigenous chickens

Objective

Cambodia is located within the distribution range of the red junglefowl, the common ancestor of domestic chickens. Although a variety of indigenous chickens have been reared in Cambodia since ancient times, their genetic characteristics have yet to be sufficiently defined. Here, we conducted a large-scale population genetic study to investigate the genetic diversity and population genetic structure of Cambodian indigenous chickens and their phylogenetic relationships with other chicken breeds and native chickens worldwide.

Methods

A Bayesian phylogenetic tree was constructed based on 625 mitochondrial DNA D-loop sequences, and Bayesian clustering analysis was performed for 666 individuals with 23 microsatellite markers, using samples collected from 28 indigenous chicken populations in 24 provinces and three commercial chicken breeds.

Results

A total of 92 haplotypes of mitochondrial D-loop sequences belonging to haplogroups A to F and J were detected in Cambodian chickens; in the indigenous chickens, haplogroup D (44.4%) was the most common, and haplogroups A (21.0%) and B (13.2%) were also dominant. However, haplogroup J, which is rare in domestic chickens but abundant in Thai red junglefowl, was found at a high frequency (14.5%), whereas the frequency of haplogroup E was considerably lower (4.6%). Population genetic structure analysis based on microsatellite markers revealed the presence of three major genetic clusters in Cambodian indigenous chickens. Their genetic diversity was relatively high, which was similar to findings reported for indigenous chickens from other Southeast Asian countries.

Conclusion

Cambodian indigenous chickens are characterized by mitochondrial D-loop haplotypes that are common to indigenous chickens throughout Southeast Asia, and may retain many of the haplotypes that originated from wild ancestral populations. These chickens exhibit high population genetic diversity, and the geographical distribution of three major clusters may be attributed to inter-regional trade and poultry transportation routes within Cambodia or international movement between Cambodia and other countries.

Mitochondrial DNA diversity and demographic history of Black-boned chickens in China

Black-boned chickens (Gallus domesticus, herein abbreviated BBCs) are well known for their unique appearance and medicinal properties and have a long breeding history in China. However, the genetic diversity and demographic history of BBCs remain unclear. In this study, we analyzed 844 mitochondrial DNA D-loop sequences, including 346 de novo sequences and 498 previously published sequences from 20 BBC breeds. We detected a generally high level of genetic diversity among the BBCs, with average haplotype and nucleotide diversities of 0.917 ± 0.0049 and 0.01422, respectively. Nucleotide diversity was highest in populations from Southwest China (0.01549 ± 0.00026), particularly in Yunnan Province (0.01624 ± 0.00025). Significant genetic divergence was detected between most breeds, particularly between Yunnan chickens and those from all other provinces. Haplogroups F and G had the highest levels of genetic diversity and were restricted to Southwest China, particularly Yunnan Province. Based on neutrality tests and mismatch distribution analyses, we did not obtain evidence for rapid population expansions and observed similar demographic histories in BBCs and local non-BBCs. Our results suggest that Chinese BBCs have complex breeding histories and may be selected in situ from local domestic chickens. These results improve our understanding of the genetic heritage and breeding histories of these desirable chickens.

Ancient DNA Evidence from China Reveals the Expansion of Pacific Dogs

The ancestral homeland of Australian dingoes and Pacific dogs is proposed to be in South China. However, the location and timing of their dispersal and relationship to dog domestication is unclear. Here, we sequenced 7,000- to 2,000-year-old complete mitochondrial DNA (mtDNA) genomes of 27 ancient canids (one gray wolf and 26 domestic dogs) from the Yellow River and Yangtze River basins (YYRB). These are the first complete ancient mtDNA of Chinese dogs from the cradle of early Chinese civilization. We found that most ancient dogs (18/26) belong to the haplogroup A1b lineage that is found in high frequency in present-day Australian dingoes and precolonial Pacific Island dogs but low frequency in present-day China. Particularly, a 7,000-year-old dog from the Tianluoshan site in Zhejiang province possesses a haplotype basal to the entire haplogroup A1b lineage. We propose that A1b lineage dogs were once widely distributed in the YYRB area. Following their dispersal to South China, and then into Southeast Asia, New Guinea and remote Oceania, they were largely replaced by dogs belonging to other lineages in the last 2,000 years in present-day China, especially North China.

Genetic variation of Nigerian cattle inferred from maternal and paternal genetic markers

The African cattle provide unique genetic resources shaped up by both diverse tropical environmental conditions and human activities, the assessment of their genetic diversity will shade light on the mechanism of their remarkable adaptive capacities. We therefore analyzed the genetic diversity of cattle samples from Nigeria using both maternal and paternal DNA markers. Nigerian cattle can be assigned to 80 haplotypes based on the mitochondrial DNA (mtDNA) D-loop sequences and haplotype diversity was 0.985 + 0.005. The network showed two major matrilineal clustering: the dominant cluster constituting the Nigerian cattle together with other African cattle while the other clustered Eurasian cattle. Paternal analysis indicates only zebu haplogroup in Nigerian cattle with high genetic diversity 1.000 ± 0.016 compared to other cattle. There was no signal of maternal genetic structure in Nigerian cattle population, which may suggest an extensive genetic intermixing within the country. The absence of Bos indicus maternal signal in Nigerian cattle is attributable to vulnerability bottleneck of mtDNA lineages and concordance with the view of male zebu genetic introgression in African cattle. Our study shades light on the current genetic diversity in Nigerian cattle and population history in West Africa.

Population Dynamics in Italian Canids between the Late Pleistocene and Bronze Age

Dog domestication is still largely unresolved due to time-gaps in the sampling of regions. Ancient Italian canids are particularly understudied, currently represented by only a few specimens. In the present study, we sampled 27 canid remains from Northern Italy dated between the Late Pleistocene and Bronze Age to assess their genetic variability, and thus add context to dog domestication dynamics. They were targeted at four DNA fragments of the hypervariable region 1 of mitochondrial DNA. A total of 11 samples had good DNA preservation and were used for phylogenetic analyses. The dog samples were assigned to dog haplogroups A, C and D, and a Late Pleistocene wolf was set into wolf haplogroup 2. We present our data in the landscape of ancient and modern dog genetic variability, with a particular focus on the ancient Italian samples published thus far. Our results suggest there is high genetic variability within ancient Italian canids, where close relationships were evident between both a ~24,700 years old Italian canid, and Iberian and Bulgarian ancient dogs. These findings emphasize that disentangling dog domestication dynamics benefits from the analysis of specimens from Southern European regions.

Geographic Origin and Genetic Characteristics of Japanese Indigenous Chickens Inferred from Mitochondrial D-Loop Region and Microsatellite DNA Markers

Japanese indigenous chickens have a long breeding history, possibly beginning 2000 years ago. Genetic characterization of Japanese indigenous chickens has been performed using mitochondrial D-loop region and microsatellite DNA markers. Their phylogenetic relationships with chickens worldwide and genetic variation within breeds have not yet been examined. In this study, the genetic characteristics of 38 Japanese indigenous chicken breeds were assessed by phylogenetic analyses of mitochondrial D-loop sequences compared with those of indigenous chicken breeds overseas. To evaluate the genetic relationships among Japanese indigenous chicken breeds, a STRUCTURE analysis was conducted using 27 microsatellite DNA markers. D-loop sequences of Japanese indigenous chickens were classified into five major haplogroups, A-E, among 15 haplogroups found in chickens worldwide. The haplogroup composition suggested that Japanese indigenous chickens originated mainly from China, with some originating from Southeast Asia. The STRUCTURE analyses revealed that Japanese indigenous chickens are genetically differentiated from chickens overseas; Japanese indigenous chicken breeds possess distinctive genetic characteristics, and Jidori breeds, which have been reared in various regions of Japan for a long time, are genetically close to each other. These results provide new insights into the history of chickens around Asia in addition to novel genetic data for the conservation of Japanese indigenous chickens.

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.

Genetic analysis identifies the missing parchment of New Zealand's founding document, the Treaty of Waitangi

Genetic analyses provide a powerful tool with which to identify the biological components of historical objects. Te Tiriti o Waitangi | The Treaty of Waitangi is New Zealand’s founding document, intended to be a partnership between the indigenous Māori and the British Crown. Here we focus on an archived piece of blank parchment that has been proposed to be the missing portion of the lower parchment of the Waitangi Sheet of the Treaty. However, its physical dimensions and characteristics are not consistent with this hypothesis. We perform genetic analyses on the parchment membranes of the Treaty, plus the blank piece of parchment. We find that all three parchments were made from ewes and that the blank parchment is highly likely to be a portion cut from the lower membrane of the Waitangi Sheet because they share identical whole mitochondrial genomes, including an unusual heteroplasmic site. We suggest that the differences in size and characteristics between the two pieces of parchment may have resulted from the Treaty’s exposure to water in the early 20^th century and the subsequent repair work, light exposure during exhibition or the later conservation treatments in the 1970s and 80s. The blank piece of parchment will be valuable for comparison tests to study the effects of earlier treatments and to monitor the effects of long-term display on the Treaty.

Phylogenetic Studies on Red Junglefowl (Gallus gallus) and Native Chicken (Gallus gallus domesticus) in Samar Island, Philippines using the Mitochondrial DNA D-Loop Region

A study was conducted to provide genetic information on the matrilineal phylogeny and genetic diversity of Red junglefowl (RJF) and native chickens in Samar Island, Philippines and to identify the genetic distance between Philippine junglefowls and other RJF species in Southeast Asia using complete mitochondrial DNA D-loop sequences. A total of 5 RJFs and 43 native chickens from Samar Island were included in this study. The results showed that Samar RJFs had a nucleotide diversity of 0.0050±0.0016, which was lower than those of three subspecies of Gallus gallus: G. g. gallus, G. g. spadiceus, and G. g. jabouillei. Meanwhile, Samar native chickens showed lower nucleotide diversity (0.0056±0.0004) than domestic fowls in some neighboring Southeast Asian countries, but higher than those in African and European countries. Phylogenetic analysis showed that 3 haplotypes of Samar RJFs clustered to haplogroup D1, and that 2 haplotypes clustered to haplogroup D2. Chickens native to Samar Island showed 100% resemblance to those in the haplogroup shared by domestic chickens and RJFs. Haplogroups A and B and sub-haplogroups D1 and E1 were the more widely distributed matrilineal lineages in Samar Island. Phylogenetic analysis of Samar RJFs showed that they were closely related to Myanmar RJFs (99.6%), Indonesia RJFs (99.5%), and Thailand RJFs (99.1%). This study is an initial investigation estimating the matrilineal phylogeny and genetic diversity of chicken populations in Samar Island, Philippines for developing strategies aimed at the future conservation and improvement of valuable genetic resources.

Genetic diversity of 21 experimental chicken lines with diverse origins and genetic backgrounds

The genetic characteristics and diversity of 21 experimental chicken lines registered with the National BioResource Project of Japan were examined using mitochondrial D-loop sequences and 54 microsatellite DNA markers. A total of 12 haplotypes were detected in the 500-bp mitochondrial DNA sequences of the hypervariable segment I for 349 individuals of 21 lines. The 12 haplotypes belonged to three (A, D, and E) haplogroups, out of the eight (A‒H) common haplogroups in domestic chickens and red junglefowls. The haplogroups A and D were widely represented in indigenous chickens in the Asian and Pacific regions, and the haplogroup E was the most prevalent in domestic chickens. Genetic clustering by discriminant analysis of principal components with microsatellite markers divided 681 individuals of 21 lines into three groups that consisted of Fayoumi-, European-, and Asian- derived lines. In each of the cladograms constructed with Nei's genetic distances based on allele frequencies and the membership coefficients provided by STRUCTURE and with the genetic distance based on the proportion of shared alleles, the genetic relationships coincided well with the breeding histories of the lines. Microsatellite markers showed remarkably lower genetic heterozygosities (less than 0.1 observed heterozygosity) for eight lines (GSP, GSN/1, YL, PNP, BM-C, WL-G, BL-E, and #413), which have been maintained as closed colonies for more than 40 years (except for #413), indicating their usefulness as experimental chicken lines in laboratory animal science research.

The association of mitochondrial DNA haplotypes and phenotypic traits in pigs

Background

The mitochondrial genome (mtDNA) is an emerging determiner of phenotypic traits and disease. mtDNA is inherited in a strict maternal fashion from the population of mitochondria present in the egg at fertilisation. Individuals are assigned to mtDNA haplotypes and those with sequences that cluster closely have common origins and their migration patterns can be mapped. Previously, we identified five mtDNA haplotypes in the commercial breeding lines of Australian pigs, which defined their common origins, and showed how these mtDNA haplotypes influenced litter size and reproductive function in terms of egg and embryo quality and fertilisation efficiency.

Results

We have determined whether mtDNA haplotypes influence other phenotypic traits. These include fat density; muscle depth; fat to leanness ratios; lifetime daily gain; teat quality; muscle score; front and rear leg assessments; percentage offspring weaned; weaning to oestrus intervals; gilt age at selection; and gestational length. In all, we assessed 5687 pigs of which 2762 were females and 2925 were males. We assessed all animals together and then by gender. We further assessed by gender based on whether a sire had joined with females from only one haplotype or from more than one haplotype. We determined that fat density, muscle depth, fat to leanness ratios, lifetime daily gain and teat quality were influenced by mtDNA haplotype and that there were gender specific effects on teat quality.

Conclusions

Our data illustrate that mtDNA haplotypes are associated with a number of important phenotypic traits indicative of economic breeding values in breeding pigs with gender-specific differences. Interestingly, there are ‘trade offs’ whereby some mtDNA haplotypes perform better for one selection criterion, such as muscle depth, but less so for another, for example teat quality, indicating that pig mtDNA haplotypes are afforded an advantage in one respect but a disadvantage in another.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0629-4) contains supplementary material, which is available to authorized users.

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.

Some maternal lineages of domestic horses may have origins in East Asia revealed with further evidence of mitochondrial genomes and HVR-1 sequences

Objectives

There are large populations of indigenous horse (Equus caballus) in China and some other parts of East Asia. However, their matrilineal genetic diversity and origin remained poorly understood. Using a combination of mitochondrial DNA (mtDNA) and hypervariable region (HVR-1) sequences, we aim to investigate the origin of matrilineal inheritance in these domestic horses.

Methods

To investigate patterns of matrilineal inheritance in domestic horses, we conducted a phylogenetic study using 31 de novo mtDNA genomes together with 317 others from the GenBank. In terms of the updated phylogeny, a total of 5,180 horse mitochondrial HVR-1 sequences were analyzed.

Results

Eightteen haplogroups (Aw-Rw) were uncovered from the analysis of the whole mitochondrial genomes. Most of which have a divergence time before the earliest domestication of wild horses (about 5,800 years ago) and during the Upper Paleolithic (35-10 KYA). The distribution of some haplogroups shows geographic patterns. The Lw haplogroup contained a significantly higher proportion of European horses than the horses from other regions, while haplogroups Jw, Rw, and some maternal lineages of Cw, have a higher frequency in the horses from East Asia. The 5,180 sequences of horse mitochondrial HVR-1 form nine major haplogroups (A-I). We revealed a corresponding relationship between the haplotypes of HVR-1 and those of whole mitochondrial DNA sequences. The data of the HVR-1 sequences also suggests that Jw, Rw, and some haplotypes of Cw may have originated in East Asia while Lw probably formed in Europe.

Conclusions

Our study supports the hypothesis of the multiple origins of the maternal lineage of domestic horses and some maternal lineages of domestic horses may have originated from East Asia.

Analysis of the genetic variation in mitochondrial DNA, Y-chromosome sequences, and MC1R sheds light on the ancestry of Nigerian indigenous pigs

Background

The history of pig populations in Africa remains controversial due to insufficient evidence from archaeological and genetic data. Previously, a Western ancestry for West African pigs was reported based on loci that are involved in the determination of coat color. We investigated the genetic diversity of Nigerian indigenous pigs (NIP) by simultaneously analyzing variation in mitochondrial DNA (mtDNA), Y-chromosome sequence and the melanocortin receptor 1 (MC1R) gene.

Results

Median-joining network analysis of mtDNA D-loop sequences from 201 NIP and previously characterized loci clustered NIP with populations from the West (Europe/North Africa) and East/Southeast Asia. Analysis of partial sequences of the Y-chromosome in 57 Nigerian boars clustered NIP into lineage HY1. Finally, analysis of MC1R in 90 NIP resulted in seven haplotypes, among which the European wild boar haplotype was carried by one individual and the European dominant black by most of the other individuals (93%). The five remaining unique haplotypes differed by a single synonymous substitution from European wild type, European dominant black and Asian dominant black haplotypes.

Conclusions

Our results demonstrate a European and East/Southeast Asian ancestry for NIP. Analyses of MC1R provide further evidence. Additional genetic analyses and archaeological studies may provide further insights into the history of African pig breeds. Our findings provide a valuable resource for future studies on whole-genome analyses of African pigs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-017-0326-1) contains supplementary material, which is available to authorized users.

Endogenous Retroviral Insertions Indicate a Secondary Introduction of Domestic Sheep Lineages to the Caucasus and Central Asia between the Bronze and Iron Age

Sheep were one of the first livestock species domesticated by humans. After initial domestication in the Middle East they were spread across Eurasia. The modern distribution of endogenous Jaagsiekte sheep retrovirus insertions in domestic sheep breeds suggests that over the course of millennia, successive introductions of improved lineages and selection for wool quality occurred in the Mediterranean region and most of Asia. Here we present a novel ancient DNA approach using data of endogenous retroviral insertions in Bronze and Iron Age domestic sheep from the Caucasus and Pamir mountain areas. Our findings support a secondary introduction of wool sheep from the Middle East between the Late Bronze Age and Iron Age into most areas of Eurasia.

Mitochondrial DNA sequence variation in Iranian native dogs

The dog mtDNA diversity picture from wide geographical sampling but from a small number of individuals per region or breed, displayed little geographical correlation and high degree of haplotype sharing between very distant breeds. For a clear picture, we extensively surveyed Iranian native dogs (n = 305) in comparison with published European (n = 443) and Southwest Asian (n = 195) dogs. Twelve haplotypes related to haplogroups A, B and C were shared by Iranian, European, Southwest Asian and East Asian dogs. In Iran, haplotype and nucleotide diversities were highest in east, southeast and northwest populations while western population had the least. Sarabi and Saluki dog populations can be assigned into haplogroups A, B, C and D; Qahderijani and Kurdi to haplogroups A, B and C, Torkaman to haplogroups A, B and D while Sangsari and Fendo into haplogroups A and B, respectively. Evaluation of population differentiation using pairwise F_ST generally revealed no clear population structure in most Iranian dog populations. The genetic signal of a recent demographic expansion was detected in East and Southeast populations. Further, in accordance with previous studies on dog-wolf hybridization for haplogroup d2 origin, the highest number of d2 haplotypes in Iranian dog as compared to other areas of Mediterranean basin suggests Iran as the probable center of its origin. Historical evidence showed that Silk Road linked Iran to countries in South East Asia and other parts of the world, which might have probably influenced effective gene flow within Iran and these regions. The medium nucleotide diversity observed in Iranian dog calls for utilization of appropriate management techniques in increasing effective population size.

Genetic evidence from mitochondrial DNA corroborates the origin of Tibetan chickens

Chicken is the most common poultry species and is important to human societies. Tibetan chicken (Gallus gallus domesticus) is a breed endemic to China that is distributed mainly on the Qinghai-Tibet Plateau. However, its origin has not been well characterized. In the present study, we sequenced partial mitochondrial DNA (mtDNA) control region of 239 and 283 samples from Tibetan and Sichuan indigenous chickens, respectively. Incorporating 1091 published sequences, we constructed the matrilineal genealogy of Tibetan chickens to further document their domestication history. We found that the genetic structure of the mtDNA haplotypes of Tibetan chickens are dominated by seven major haplogroups (A-G). In addition, phylogenetic and network analyses showed that Tibetan chickens are not distinguishable from the indigenous chickens in surrounding areas. Furthermore, some clades of Tibetan chickens may have originated from game fowls. In summary, our results collectively indicated that Tibetan chickens may have diverged from indigenous chickens in the adjacent regions and hybridized with various chickens.

On the origins and genetic diversity of South American chickens: one step closer

Local chicken populations are a major source of food in the rural areas of South America. However, very little is known about their genetic composition and diversity. Here, we analyzed five populations from South America to investigate their maternal genetic origin and diversity, hoping to mitigate the lack of information on local chicken populations from this region. We also included three populations of chicken from the Iberian Peninsula and one from Easter Island, which are potential sources of the first chickens introduced in South America. The obtained sequencing data from South American chickens indicate the presence of four haplogroups (A, B, E and D) that can be further subdivided into nine sub-haplogroups. Of these, four (B1, D1a, E1a(b), E1b) were absent from local Iberian Peninsula chickens and one (D1a) was present only on Easter Island. The presence of the sub-haplogroups A1a(b) and E1a(b) in South America, previously only observed in Eastern Asia, and the significant population differentiation between Iberian Peninsula and South American populations, suggest a second maternal source of the extant genetic pool in South American chickens.

A cryptic mitochondrial DNA link between North European and West African dogs

Domestic dogs have an ancient origin and a long history in Africa. Nevertheless, the timing and sources of their introduction into Africa remain enigmatic. Herein, we analyse variation in mitochondrial DNA (mtDNA) D-loop sequences from 345 Nigerian and 37 Kenyan village dogs plus 1530 published sequences of dogs from other parts of Africa, Europe and West Asia. All Kenyan dogs can be assigned to one of three haplogroups (matrilines; clades): A, B, and C, while Nigerian dogs can be assigned to one of four haplogroups A, B, C, and D. None of the African dogs exhibits a matrilineal contribution from the African wolf (Canis lupus lupaster). The genetic signal of a recent demographic expansion is detected in Nigerian dogs from West Africa. The analyses of mitochondrial genomes reveal a maternal genetic link between modern West African and North European dogs indicated by sub-haplogroup D1 (but not the entire haplogroup D) coalescing around 12,000 years ago. Incorporating molecular anthropological evidence, we propose that sub-haplogroup D1 in West African dogs could be traced back to the late-glacial dispersals, potentially associated with human hunter-gatherer migration from southwestern Europe.

A whole mitochondria analysis of the Tyrolean Iceman's leather provides insights into the animal sources of Copper Age clothing

The attire of the Tyrolean Iceman, a 5,300-year-old natural mummy from the Ötzal Italian Alps, provides a surviving example of ancient manufacturing technologies. Research into his garments has however, been limited by ambiguity surrounding their source species. Here we present a targeted enrichment and sequencing of full mitochondrial genomes sampled from his clothes and quiver, which elucidates the species of production for nine fragments. Results indicate that the majority of the samples originate from domestic ungulate species (cattle, sheep and goat), whose recovered haplogroups are now at high frequency in today’s domestic populations. Intriguingly, the hat and quiver samples were produced from wild species, brown bear and roe deer respectively. Combined, these results suggest that Copper Age populations made considered choices of clothing material from both the wild and domestic populations available to them. Moreover, these results show the potential for the recovery of complete mitochondrial genomes from degraded prehistoric artefacts.

Genomic and archaeological evidence suggest a dual origin of domestic dogs

The geographic and temporal origins of dogs remain controversial. We generated genetic sequences from 59 ancient dogs and a complete (28x) genome of a late Neolithic dog (dated to ~4800 calendar years before the present) from Ireland. Our analyses revealed a deep split separating modern East Asian and Western Eurasian dogs. Surprisingly, the date of this divergence (~14,000 to 6400 years ago) occurs commensurate with, or several millennia after, the first appearance of dogs in Europe and East Asia. Additional analyses of ancient and modern mitochondrial DNA revealed a sharp discontinuity in haplotype frequencies in Europe. Combined, these results suggest that dogs may have been domesticated independently in Eastern and Western Eurasia from distinct wolf populations. East Eurasian dogs were then possibly transported to Europe with people, where they partially replaced European Paleolithic dogs.

Genetic diversity of the Hungarian Gidran horse in two mitochondrial DNA markers

The Gidran is a native Hungarian horse breed that has approached extinction several times. Phylogenetic analysis of two mitochondrial markers (D-loop and cytochrome-b) was performed to determine the genetic characterization of the Gidran for the first time as well as to detect errors in the management of the Gidran stud book. Sequencing of 686 bp of CYTB and 202 bp of the D-loop in 260 mares revealed 24 and 32 haplotypes, respectively, among 31 mare families. BLAST analysis revealed six novel CYTB and four D-loop haplotypes that have not been previously reported. The Gidran mares showed high haplotype (CYTB: 0.8735 ± 0.011; D-loop: 0.9136 ± 0.008) and moderate nucleotide (CYTB: 0.00472 ± 0.00017; D-loop: 0.02091 ± 0.00068) diversity. Of the 31 Gidran mare families, only 15 CYTB (48.4%) and 17 D-loop (54.8%) distinct haplotypes were formed using the two markers separately. Merged markers created 24 (77.4%) mare families, which were in agreement with the mare families in the stud book. Our key finding was that the Gidran breed still possesses high genetic diversity despite its history. The obtained haplotypes are mostly consistent with known mare families, particularly when the two mtDNA markers were merged. Our results could facilitate conservation efforts for preserving the genetic diversity of the Gidran.