Explore the vast range of titles available.

The domestication of horses was very important in the history of humankind. However, the ancestry of modern horses and the location and timing of their emergence remain unclear. Gaunitz et al. generated 42 ancient-horse genomes. Their source samples included the Botai archaeological site in Central Asia, considered to include the earliest domesticated horses. Unexpectedly, Botai horses were the ancestors not of modern domestic horses, but rather of modern Przewalski's horses. Thus, in contrast to current thinking on horse domestication, modern horses may have been domesticated in other, more Western, centers of origin. Science , this issue p. 111 The earliest herded horses were ancestors of feral Przewalski’s horses but not of modern domesticated horses. The Eneolithic Botai culture of the Central Asian steppes provides the earliest archaeological evidence for horse husbandry, ~5500 years ago, but the exact nature of early horse domestication remains controversial. We generated 42 ancient-horse genomes, including 20 from Botai. Compared to 46 published ancient- and modern-horse genomes, our data indicate that Przewalski’s horses are the feral descendants of horses herded at Botai and not truly wild horses. All domestic horses dated from ~4000 years ago to present only show ~2.7% of Botai-related ancestry. This indicates that a massive genomic turnover underpins the expansion of the horse stock that gave rise to modern domesticates, which coincides with large-scale human population expansions during the Early Bronze Age.

Background To date, genome-scale analyses in the domestic horse have been limited by suboptimal single nucleotide polymorphism (SNP) density and uneven genomic coverage of the current SNP genotyping arrays. The recent availability of whole genome sequences has created the opportunity to develop a next generation, high-density equine SNP array. Results Using whole genome sequence from 153 individuals representing 24 distinct breeds collated by the equine genomics community, we cataloged over 23 million de novo discovered genetic variants. Leveraging genotype data from individuals with both whole genome sequence, and genotypes from lower-density, legacy SNP arrays, a subset of ~5 million high-quality, high-density array candidate SNPs were selected based on breed representation and uniform spacing across the genome. Considering probe design recommendations from a commercial vendor (Affymetrix, now Thermo Fisher Scientific) a set of ~2 million SNPs were selected for a next-generation high-density SNP chip (MNEc2M). Genotype data were generated using the MNEc2M array from a cohort of 332 horses from 20 breeds and a lower-density array, consisting of ~670 thousand SNPs (MNEc670k), was designed for genotype imputation. Conclusions Here, we document the steps taken to design both the MNEc2M and MNEc670k arrays, report genomic and technical properties of these genotyping platforms, and demonstrate the imputation capabilities of these tools for the domestic horse.

Background The sample ascertainment bias due to complex population structures remains a major challenge in genome-wide investigations of complex traits. In this study we derived the high-resolution population structure and levels of autozygosity of 377 Lipizzan horses originating from five different European stud farms utilizing the SNP genotype information of the high density 700 k Affymetrix Axiom™ Equine genotyping array. Scanning the genome for overlapping runs of homozygosity (ROH) shared by more than 50% of horses, we identified homozygous regions (ROH islands) in order to investigate the gene content of those candidate regions by gene ontology and enrichment analyses. Results The high-resolution population network approach revealed well-defined substructures according to the origin of the horses (Austria, Slovakia, Croatia and Hungary). The highest mean genome coverage of ROH (S ROH ) was identified in the Austrian (S ROH  = 342.9), followed by Croatian (S ROH  = 214.7), Slovakian (S ROH  = 205.1) and Hungarian (S ROH  = 171.5) subpopulations. ROH island analysis revealed five common islands on ECA11 and ECA14, hereby confirming a closer genetic relationship between the Hungarian and Croatian as well as between the Austrian and Slovakian samples. Private islands were detected for the Hungarian and the Austrian Lipizzan subpopulations. All subpopulations shared a homozygous region on ECA11, nearly identical in position and length containing among other genes the homeobox-B cluster, which was also significantly ( p  < 0.001) highlighted by enrichment analysis. Gene ontology terms were mostly related to biological processes involved in embryonic morphogenesis and anterior/posterior specification. Around the STX17 gene (causative for greying), we identified a ROH island harbouring the genes NR4A3 , STX17 , ERP44 and INVS . Within further islands on ECA14, ECA16 and ECA20 we detected the genes SPRY4 , NDFIP1 , IMPDH2 , HSP90AB1 , whereas SPRY4 and HSP90AB1 are involved in melanoma metastasis and survival rate of melanoma patients in humans. Conclusions We demonstrated that the assessment of high-resolution population structures within one single breed supports the downstream genetic analyses (e.g. the identification of ROH islands). By means of ROH island analyses, we identified the genes SPRY4 , NDFIP1 , IMPDH2 , HSP90AB1 , which might play an important role for further studies on equine melanoma. Furthermore, our results highlighted the impact of the homeobox-A and B cluster involved in morphogenesis of Lipizzan horses.

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 ( F ST ) 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.

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.

Dangerous damage to mitochondrial DNA (mtDNA) can be ameliorated during mammalian development through a highly debated mechanism called the mtDNA bottleneck. Uncertainty surrounding this process limits our ability to address inherited mtDNA diseases. We produce a new, physically motivated, generalisable theoretical model for mtDNA populations during development, allowing the first statistical comparison of proposed bottleneck mechanisms. Using approximate Bayesian computation and mouse data, we find most statistical support for a combination of binomial partitioning of mtDNAs at cell divisions and random mtDNA turnover, meaning that the debated exact magnitude of mtDNA copy number depletion is flexible. New experimental measurements from a wild-derived mtDNA pairing in mice confirm the theoretical predictions of this model. We analytically solve a mathematical description of this mechanism, computing probabilities of mtDNA disease onset, efficacy of clinical sampling strategies, and effects of potential dynamic interventions, thus developing a quantitative and experimentally-supported stochastic theory of the bottleneck. Mitochondria are structures that provide vital sources of energy in our cells. DNA contained within mitochondria encodes important mitochondrial machinery, and most human cells contain hundreds or thousands of mitochondrial DNA molecules in addition to the DNA that is stored in the nucleus. Mitochondrial DNA is inherited from mothers via the egg, and the details of this inheritance are poorly understood. This question is important because inherited mistakes in mitochondrial DNA can have detrimental consequences on health, with links to fatal diseases and many other conditions. An unfertilised egg cell contains many copies of mitochondrial DNA molecules; some may have mutations and some may not. After fertilisation, the egg divides, the number of cells in the developing embryo increases, and the number of mitochondrial DNA molecules per cell changes. If the original egg cell contained defective mitochondrial DNA, some of these new cells end up containing more defective copies than others, leading to cell-to-cell differences in the developing embryo. This potentially allows cells with the greatest number of defective mitochondria to be eliminated. The increase in this cell-to-cell variability is called ‘bottlenecking’, and its mechanism remains highly debated. Johnston et al. have now used tools from maths, statistics and new experiments to address this debate, in the light of several studies that measured the mitochondrial DNA content in developing mice. This approach allowed a new theoretical model of mitochondrial DNA during the growth of an organism to be produced, which encompasses a wide range of existing theories and allows them to be compared. This model starts from the viewpoint that the hundreds or thousands of mitochondrial DNA molecules in a cell can be thought of as a population undergoing random ‘birth’ and ‘death’, and it allows the first statistical comparison of the many proposed bottleneck mechanisms. Johnston et al. find support for two ways that cells segregate mitochondria as they multiply, and show that the decrease in the number of mitochondrial DNA molecules during bottlenecking is flexible. This reconciles a debate amongst previous studies. These findings are confirmed using new experimental data from mice, which are genetically distinct from existing studies, illustrating the generality of the model's findings. Furthermore, an analytic mathematical description that describes in detail how bottlenecking might work is produced. Finally, Johnston et al. provide examples using this new theoretical model to suggest therapeutic strategies for diseases caused by mitochondrial DNA mutations. Future work will need to test these suggestions, and link mathematical understanding of mitochondria with healthcare data.

Analysis of the Y chromosome is the best-established way to reconstruct paternal family history in humans. Here, we applied fine-scaled Y-chromosomal haplotyping in horses with biallelic markers and demonstrate the potential of our approach to address the ancestry of sire lines. We de novo assembled a draft reference of the male-specific region of the Y chromosome from Illumina short reads and then screened 5.8 million basepairs for variants in 130 specimens from intensively selected and rural breeds and nine Przewalski’s horses. Among domestic horses we confirmed the predominance of a young’crown haplogroup’ in Central European and North American breeds. Within the crown, we distinguished 58 haplotypes based on 211 variants, forming three major haplogroups. In addition to two previously characterised haplogroups, one observed in Arabian/Coldblooded and the other in Turkoman/Thoroughbred horses, we uncovered a third haplogroup containing Iberian lines and a North African Barb Horse. In a genealogical showcase, we distinguished the patrilines of the three English Thoroughbred founder stallions and resolved a historic controversy over the parentage of the horse ‘Galopin’, born in 1872. We observed two nearly instantaneous radiations in the history of Central and Northern European Y-chromosomal lineages that both occurred after domestication 5,500 years ago.

Vital mitochondrial DNA (mtDNA) populations exist in cells and may consist of heteroplasmic mixtures of mtDNA types. The evolution of these heteroplasmic populations through development, ageing, and generations is central to genetic diseases, but is poorly understood in mammals. Here we dissect these population dynamics using a dataset of unprecedented size and temporal span, comprising 1947 single-cell oocyte and 899 somatic measurements of heteroplasmy change throughout lifetimes and generations in two genetically distinct mouse models. We provide a novel and detailed quantitative characterisation of the linear increase in heteroplasmy variance throughout mammalian life courses in oocytes and pups. We find that differences in mean heteroplasmy are induced between generations, and the heteroplasmy of germline and somatic precursors diverge early in development, with a haplotype-specific direction of segregation. We develop stochastic theory predicting the implications of these dynamics for ageing and disease manifestation and discuss its application to human mtDNA dynamics. Mitochondrial populations in cells may consist of heteroplasmic mixtures of mtDNA types, and their evolution through development, aging and generations is central to genetic diseases. Here the authors dissect these population dynamics using a large mouse-based data set to characterise the dynamics of heteroplasmy mean and variance throughout life and across generations.

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.

The paternally inherited Y chromosome displays the population genetic history of males. While modern domestic horses (Equus caballus) exhibit abundant diversity within maternally inherited mitochondrial DNA, no significant Y-chromosomal sequence diversity has been detected. We used high throughput sequencing technology to identify the first polymorphic Y-chromosomal markers useful for tracing paternal lines. The nucleotide variability of the modern horse Y chromosome is extremely low, resulting in six haplotypes (HT), all clearly distinct from the Przewalski horse (E. przewalskii). The most widespread HT1 is ancestral and the other five haplotypes apparently arose on the background of HT1 by mutation or gene conversion after domestication. Two haplotypes (HT2 and HT3) are widely distributed at high frequencies among modern European horse breeds. Using pedigree information, we trace the distribution of Y-haplotype diversity to particular founders. The mutation leading to HT3 occurred in the germline of the famous English Thoroughbred stallion \"Eclipse\" or his son or grandson and its prevalence demonstrates the influence of this popular paternal line on modern sport horse breeds. The pervasive introgression of Thoroughbred stallions during the last 200 years to refine autochthonous breeds has strongly affected the distribution of Y-chromosomal variation in modern horse breeds and has led to the replacement of autochthonous Y chromosomes. Only a few northern European breeds bear unique variants at high frequencies or fixed within but not shared among breeds. Our Y-chromosomal data complement the well established mtDNA lineages and document the male side of the genetic history of modern horse breeds and breeding practices.