Prehistoric genomes reveal the genetic foundation and cost of horse domestication

An endothelial regulatory module links blood pressure regulation with elite athletic performance

The control of transcription is crucial for homeostasis in mammals. A previous selective sweep analysis of horse racing performance revealed a 19.6 kb candidate regulatory region 50 kb downstream of the Endothelin3 (EDN3) gene. Here, the region was narrowed to a 5.5 kb span of 14 SNVs, with elite and sub-elite haplotypes analyzed for association to racing performance, blood pressure and plasma levels of EDN3 in Coldblooded trotters and Standardbreds. Comparative analysis of human HiCap data identified the span as an enhancer cluster active in endothelial cells, interacting with genes relevant to blood pressure regulation. Coldblooded trotters with the sub-elite haplotype had significantly higher blood pressure compared to horses with the elite performing haplotype during exercise. Alleles within the elite haplotype were part of the standing variation in pre-domestication horses, and have risen in frequency during the era of breed development and selection. These results advance our understanding of the molecular genetics of athletic performance and vascular traits in both horses and humans.

Genome-wide analysis of population structure, effective population size and inbreeding in Iranian and exotic horses

Population structure and genetic diversity are the key parameters to study the breeding history of animals. This research aimed to provide a characterization of the population structure and to compare the effective population size (Ne), LD decay, genetic diversity, and genomic inbreeding in Iranian native Caspian (n = 38), Turkmen (n = 24) and Kurdish (n = 29) breeds and some other exotic horses consisting of Arabian (n = 24), Fell pony (n = 21) and Akhal-Teke (n = 20). A variety of statistical population analysis techniques, such as principal component analysis (PCA), discriminant analysis of principal component (DAPC) and model-based method (STRUCTURE) were employed. The results of the population analysis clearly demonstrated a distinct separation of native and exotic horse breeds and clarified the relationships between studied breeds. The effective population size (Ne) for the last six generations was estimated 54, 49, 37, 35, 27 and 26 for the Caspian, Kurdish, Arabian, Turkmen, Akhal-Teke and Fell pony breeds, respectively. The Caspian breed showed the lowest LD with an average r2 value of 0.079, while the highest was observed in Fell pony (0.148). The highest and lowest average observed heterozygosity were found in the Kurdish breeds (0.346) and Fell pony (0.290) breeds, respectively. The lowest genomic inbreeding coefficient based on run of homozygosity (FROH) and excess of homozygosity (FHOM) was in the Caspian and Kurdish breeds, respectively, while based on genomic relationship matrix) FGRM) and correlation between uniting gametes) FUNI) the lowest genomic inbreeding coefficient was found in the Kurdish breed. The estimation of genomic inbreeding rates in the six breeds revealed that FROH yielded lower estimates compared to the other three methods. Additionally, the Iranian breeds displayed lower levels of inbreeding compared to the exotic breeds. Overall, the findings of this study provide valuable insights for the development of effective breeding management strategies aimed at preserving these horse breeds.

Genomic evidence for domestication selection in three hatchery populations of Chinook salmon, Oncorhynchus tshawytscha

Fish hatcheries are widely used to enhance fisheries and supplement declining wild populations. However, substantial evidence suggests that hatchery fish are subject to differential selection pressures compared to their wild counterparts. Domestication selection, or adaptation to the hatchery environment, poses a risk to wild populations if traits specific to success in the hatchery environment have a genetic component and there is subsequent introgression between hatchery and wild fish. Few studies have investigated domestication selection in hatcheries on a genomic level, and even fewer have done so in parallel across multiple hatchery-wild population pairs. In this study, we used low-coverage whole-genome sequencing to investigate signals of domestication selection in three separate hatchery populations of Chinook salmon, Oncorhynchus tshawytscha, after approximately seven generations of divergence from their corresponding wild progenitor populations. We sequenced 192 individuals from populations across Southeast Alaska and estimated genotype likelihoods at over six million loci. We discovered a total of 14 outlier peaks displaying high genetic differentiation (F ST) between hatchery-wild pairs, although no peaks were shared across the three comparisons. Peaks were small (53 kb on average) and often displayed elevated absolute genetic divergence (D xy) and linkage disequilibrium, suggesting some level of domestication selection has occurred. Our study provides evidence that domestication selection can lead to genetic differences between hatchery and wild populations in only a few generations. Additionally, our data suggest that population-specific adaptation to hatchery environments likely occurs through different genetic pathways, even for populations with similar standing genetic variation. These results highlight the need to collect paired genotype-phenotype data to understand how domestication may be affecting fitness and to identify potential management practices that may mitigate genetic risks despite multiple pathways of domestication.

Koban culture genome-wide and archeological data open the bridge between Bronze and Iron Ages in the North Caucasus

The North Caucasus played a key role during the ancient colonization of Eurasia and the formation of its cultural and genetic ancestry. Previous archeogenetic studies described a relative genetic and cultural continuity of ancient Caucasus societies, since the Eneolithic period. The Koban culture, which formed in the Late Bronze Age on the North Caucasian highlands, is considered as a cultural "bridge" between the ancient and modern autochthonous peoples of the Caucasus. Here, we discuss the place of this archeological culture and its representatives in the genetic orbit of Caucasian cultures using genome-wide SNP data from five individuals of the Koban culture and one individual of the early Alanic culture as well as previously published genomic data of ancient and modern North Caucasus individuals. Ancient DNA analysis shows that an ancient individual from Klin-Yar III, who was previously described as male, was in fact a female. Additional studies on well-preserved ancient human specimens are necessary to determine the level of local mobility and kinship between individuals in ancient societies of North Caucasus. Further studies with a larger sample size will allow us gain a deeper understanding of this topic.

The genomic landscape of mammal domestication might be orchestrated by selected transcription factors regulating brain and craniofacial development

Domestication transforms once wild animals into tamed animals that can be then exploited by humans. The process entails modifications in the body, cognition, and behavior that are essentially driven by differences in gene expression patterns. Although genetic and epigenetic mechanisms were shown to underlie such differences, less is known about the role exerted by trans-regulatory molecules, notably transcription factors (TFs) in domestication. In this paper, we conducted extensive in silico analyses aimed to clarify the TF landscape of mammal domestication. We first searched the literature, so as to establish a large list of genes selected with domestication in mammals. From this list, we selected genes experimentally demonstrated to exhibit TF functions. We also considered TFs displaying a statistically significant number of targets among the entire list of (domestication) selected genes. This workflow allowed us to identify 5 candidate TFs (SOX2, KLF4, MITF, NR3C1, NR3C2) that were further assessed in terms of biochemical and functional properties. We found that such TFs-of-interest related to mammal domestication are all significantly involved in the development of the brain and the craniofacial region, as well as the immune response and lipid metabolism. A ranking strategy, essentially based on a survey of protein-protein interactions datasets, allowed us to identify SOX2 as the main candidate TF involved in domestication-associated evolutionary changes. These findings should help to clarify the molecular mechanics of domestication and are of interest for future studies aimed to understand the behavioral and cognitive changes associated to domestication.

Single nucleotide polymorphisms (SNPs) and indels identified from whole-genome re-sequencing of four Chinese donkey breeds

This paper represents the fundamental report of the survey of genome-wide changes of four Chinese indigenous donkey breeds, Dezhou (DZ), Guangling (GL), North China (NC), and Shandong Little donkey (SDL), and the findings will prove usefully for identification of biomarkers that perhaps predict or characterize the growth and coat color patterns. Three genomic regions in CYP3A12, TUBGCP5, and GSTA1 genes, were identified as putative selective sweeps in all researched donkey populations. The loci of candidate genes that may have contributed to the phenotypes in body size (ACSL4, MSI2, ADRA1B, and CDKL5) and coat color patterns (KITLG and TBX3) in donkey populations would be found in underlying strong selection signatures when compared between large and small donkey types, and between different coat colors. The results of the phylogenetic analysis, FST, and principal component analysis (PCA) supported that each population cannot clearly deviate from each other, showing no obvious population structure. We can conclude from the population history that the formation processes between DZS and NC, GL, and SDL are completely different. The genetic variants discovered here provide a rich resource to help identify potential genomic markers and their associated molecular mechanisms that impact economically important traits for Chinese donkey breeding programs.

Allele biased transcription factor binding across human brain regions gives mechanistic insight into eQTLs

Transcription Factors (TFs) influence gene expression by facilitating or disrupting the formation of transcription initiation machinery at particular genomic loci. Because genomic localization of TFs is in part driven by TF recognition of DNA sequence, variation in TF binding sites can disrupt TF-DNA associations and affect gene regulation. To identify variants that impact TF binding in human brain tissues, we quantified allele bias for 93 TFs analyzed with ChIP-seq experiments of multiple structural brain regions from two donors. Using graph genomes constructed from phased genomic sequence data, we compared ChIP-seq signal between alleles at heterozygous variants within each tissue sample from each donor. Comparison of results from different brain regions within donors and the same regions between donors provided measures of allele bias reproducibility. We identified thousands of DNA variants that show reproducible bias in ChIP-seq for at least one TF. We found that alleles that are rarer in the general population were more likely than common alleles to exhibit large biases, and more frequently led to reduced TF binding. Combining ChIP-seq with RNA-seq, we identified TF-allele interaction biases with RNA bias in a phased allele linked to 6,709 eQTL variants identified in GTEx data, 3,309 of which were found in neural contexts. Our results provide insights into the effects of both common and rare variation on gene regulation in the brain. These findings can facilitate mechanistic understanding of cis-regulatory variation associated with biological traits, including disease.

Genetic variation and domestication of horses revealed by 10 chromosome-level genomes and whole-genome resequencing

Understanding the genetic variations of the horse (Equus caballus) genome will improve breeding conservation and welfare. However, genetic variations in long segments, such as structural variants (SVs), remain understudied. We de novo assembled 10 chromosome-level three-dimensional horse genomes, each representing a distinct breed, and analysed horse SVs using a multi-assembly approach. Our findings suggest that SVs with the accumulation of mammalian-wide interspersed repeats related to long interspersed nuclear elements might be a horse-specific mechanism to modulate genome-wide gene regulatory networks. We found that olfactory receptors were commonly loss and accumulated deleterious mutations, but no purge of deleterious mutations occurred during horse domestication. We examined the potential effects of SVs on the spatial structure of chromatin via topologically associating domains (TADs). Breed-specific TADs were significantly enriched by breed-specific SVs. We identified 4199 unique breakpoint-resolved novel insertions across all chromosomes that account for 2.84 Mb sequences missing from the reference genome. Several novel insertions might have potential functional consequences, as 519 appeared to reside within 449 gene bodies. These genes are primarily involved in pathogen recognition, innate immune responses and drug metabolism. Moreover, 37 diverse horses were resequenced. Combining this with public data, we analysed 97 horses through a comparative population genomics approach to identify the genetic basis underlying breed characteristics using Thoroughbreds as a case study. We provide new scientific evidence for horse domestication, an understanding of the genetic mechanism underlying the phenotypic evolution of horses, and a comprehensive genetic variation resource for further genetic studies of horses.

Characterization of the Caudal Ventral Tubercle in the Sixth Cervical Vertebra in Modern Equus ferus caballus

This study examined the anomalous variations of the ventral process of C6 in modern E. ferus caballus. The aim was to provide an incremental grading protocol measuring the absence of the caudal ventral tubercle (CVT) in this ventral process. The findings revealed the most prevalent absent CVT (aCVT) was left unilateral (n = 35), with bilateral (n = 29) and right unilateral (n = 12). Grading was determined in equal increments of absence 1/4, 2/4, 3/4, with 4/4 representing a complete aCVT in 56/76, with a significance of p = 0.0013. This also applied to bilateral specimens. In those C6 osseous specimens displaying a 4/4 grade aCVT, 41/56 had a partial absence of the caudal aspect of the cranial ventral tubercle (CrVT). Here, grading absent CrVTs (aCrVT) followed similarly to aCVTs, though 4/4 was not observed. The significance between 4/4 grade aCVTs and the presentation of an aCrVT was left p = 0.00001 and right p = 0.00018. In bilateral specimens, C6 morphologically resembled C5, implying a homeotic transformation that limited the attachment sites for the cranial and thoracal longus colli muscle. This potentially diminishes function and caudal cervical stability. Therefore, it is recommended that further studies examine the morphological extent of this equine complex vertebral malformation (ECVM) as well as its interrelationships and genetic code/blueprint.

Imputed genomes of historical horses provide insights into modern breeding

Historical genomes can provide important insights into recent genomic changes in horses, especially the development of modern breeds. In this study, we characterized 8.7 million genomic variants from a panel of 430 horses from 73 breeds, including newly sequenced genomes from 20 Clydesdales and 10 Shire horses. We used this modern genomic variation to impute the genomes of four historically important horses, consisting of publicly available genomes from 2 Przewalski's horses, 1 Thoroughbred, and a newly sequenced Clydesdale. Using these historical genomes, we identified modern horses with higher genetic similarity to those in the past and unveiled increased inbreeding in recent times. We genotyped variants associated with appearance and behavior to uncover previously unknown characteristics of these important historical horses. Overall, we provide insights into the history of Thoroughbred and Clydesdale breeds and highlight genomic changes in the endangered Przewalski's horse following a century of captive breeding.

High-Resolution Genotyping of Expressed Equine MHC Reveals a Highly Complex MHC Structure

The Major Histocompatibility Complex (MHC) genes play a key role in a number of biological processes, most notably in immunological responses. The MHCI and MHCII genes incorporate a complex set of highly polymorphic and polygenic series of genes, which, due to the technical limitations of previously available technologies, have only been partially characterized in non-model but economically important species such as the horse. The advent of high-throughput sequencing platforms has provided new opportunities to develop methods to generate high-resolution sequencing data on a large scale and apply them to the analysis of complex gene sets such as the MHC. In this study, we developed and applied a MiSeq-based approach for the combined analysis of the expressed MHCI and MHCII repertoires in cohorts of Thoroughbred, Icelandic, and Norwegian Fjord Horses. The approach enabled us to generate comprehensive MHCI/II data for all of the individuals (n = 168) included in the study, identifying 152 and 117 novel MHCI and MHCII sequences, respectively. There was limited overlap in MHCI and MHCII haplotypes between the Thoroughbred and the Icelandic/Norwegian Fjord horses, showcasing the variation in MHC repertoire between genetically divergent breeds, and it can be inferred that there is much more MHC diversity in the global horse population. This study provided novel insights into the structure of the expressed equine MHC repertoire and highlighted unique features of the MHC in horses.

Shared reproductive disruption, not neural crest or tameness, explains the domestication syndrome

Altered neural crest cell (NCC) behaviour is an increasingly cited explanation for the domestication syndrome in animals. However, recent authors have questioned this explanation, while others cast doubt on whether domestication syndrome even exists. Here, we review published literature concerning this syndrome and the NCC hypothesis, together with recent critiques of both. We synthesize these contributions and propose a novel interpretation, arguing shared trait changes under ancient domestication resulted primarily from shared disruption of wild reproductive regimes. We detail four primary selective pathways for 'reproductive disruption' under domestication and contrast these succinct and demonstrable mechanisms with cryptic genetic associations posited by the NCC hypothesis. In support of our perspective, we illustrate numerous important ways in which NCCs contribute to vertebrate reproductive phenotypes, and argue it is not surprising that features derived from these cells would be coincidentally altered under major selective regime changes, as occur in domestication. We then illustrate several pertinent examples of Darwin's 'unconscious selection' in action, and compare applied selection and phenotypic responses in each case. Lastly, we explore the ramifications of reproductive disruption for wider evolutionary discourse, including links to wild 'self-domestication' and 'island effect', and discuss outstanding questions.

Deleterious Variation in Natural Populations and Implications for Conservation Genetics

Deleterious mutations decrease reproductive fitness and are ubiquitous in genomes. Given that many organisms face ongoing threats of extinction, there is interest in elucidating the impact of deleterious variation on extinction risk and optimizing management strategies accounting for such mutations. Quantifying deleterious variation and understanding the effects of population history on deleterious variation are complex endeavors because we do not know the strength of selection acting on each mutation. Further, the effect of demographic history on deleterious mutations depends on the strength of selection against the mutation and the degree of dominance. Here we clarify how deleterious variation can be quantified and studied in natural populations. We then discuss how different demographic factors, such as small population size, nonequilibrium population size changes, inbreeding, and gene flow, affect deleterious variation. Lastly, we provide guidance on studying deleterious variation in nonmodel populations of conservation concern.

MetaDamage tool: Examining post-mortem damage in sedaDNA on a metagenomic scale

The use of metagenomic datasets to support ancient sedimentary DNA (sedaDNA) for paleoecological reconstruction has been demonstrated to be a powerful tool to understand multi-organism responses to climatic shifts and events. Authentication remains integral to the ancient DNA discipline, and this extends to sedaDNA analysis. Furthermore, distinguishing authentic sedaDNA from contamination or modern material also allows for a better understanding of broader questions in sedaDNA research, such as formation processes, source and catchment, and post-depositional processes. Existing tools for the detection of damage signals are designed for single-taxon input, require a priori organism specification, and require a significant number of input sequences to establish a signal. It is therefore often difficult to identify an established cytosine deamination rate consistent with ancient DNA across a sediment sample. In this study, we present MetaDamage, a tool that examines cytosine deamination on a metagenomic (all organisms) scale for multiple previously undetermined taxa and can produce a damage profile based on a few hundred reads. We outline the development and testing of the MetaDamage tool using both authentic sedaDNA sequences and simulated data to demonstrate the resolution in which MetaDamage can identify deamination levels consistent with the presence of ancient DNA. The MetaDamage tool offers a method for the initial assessment of the presence of sedaDNA and a better understanding of key questions of preservation for paleoecological reconstruction.

Common protein-coding variants influence the racing phenotype in galloping racehorse breeds

Selection for system-wide morphological, physiological, and metabolic adaptations has led to extreme athletic phenotypes among geographically diverse horse breeds. Here, we identify genes contributing to exercise adaptation in racehorses by applying genomics approaches for racing performance, an end-point athletic phenotype. Using an integrative genomics strategy to first combine population genomics results with skeletal muscle exercise and training transcriptomic data, followed by whole-genome resequencing of Asian horses, we identify protein-coding variants in genes of interest in galloping racehorse breeds (Arabian, Mongolian and Thoroughbred). A core set of genes, G6PC2, HDAC9, KTN1, MYLK2, NTM, SLC16A1 and SYNDIG1, with central roles in muscle, metabolism, and neurobiology, are key drivers of the racing phenotype. Although racing potential is a multifactorial trait, the genomic architecture shaping the common athletic phenotype in horse populations bred for racing provides evidence for the influence of protein-coding variants in fundamental exercise-relevant genes. Variation in these genes may therefore be exploited for genetic improvement of horse populations towards specific types of racing.

Developmental instability in domesticated mammals

Measures of fluctuating asymmetry (FA) have been adopted widely as an estimate of developmental instability. Arising from various sources of stress, developmental instability is associated with an organism's capacity to maintain fitness. The process of domestication has been framed as an environmental stress with human-specified parameters, suggesting that FA may manifest to a larger degree among domesticates compared to their wild relatives. This study used three-dimensional geometric morphometric landmark data to (a) quantify the amount of FA in the cranium of six domestic mammal species and their wild relatives and, (b) provide novel assessment of the commonalities and differences across domestic/wild pairs concerning the extent to which random variation arising from the developmental system assimilates into within-population variation. The majority of domestic mammals showed greater disparity for asymmetric shape, however, only two forms (Pig, Dog) showed significantly higher disparity as well as a higher degree of asymmetry compared to their wild counterparts (Wild Boar, Wolf). Contra to predictions, most domestic and wild forms did not show a statistically significant correspondence between symmetric shape variation and FA, however, a moderate correlation value was recorded for most pairs (r-partial least squares >0.5). Within pairs, domestic and wild forms showed similar correlation magnitudes for the relationship between the asymmetric and symmetric components. In domesticates, new variation may therefore retain a general, conserved pattern in the gross structuring of the cranium, whilst also being a source for response to selection on specific features.

Analysis of ancient and modern horse genomes reveals the critical impact of lncRNA-mediated epigenetic regulation on horse domestication

Background: The domestication of horses has played critical roles in human civilizations. The excavation of ancient horse DNA provides crucial data for studying horse domestication. Studies of horse domestication can shed light on the general mechanisms of animal domestication.

Objective: We wish to explore the gene transcription regulation by long noncoding RNAs (lncRNAs) that influence horse domestication.

Methods: First, we assembled the ancient DNA sequences of multiple horses at different times and the genomes of horses, donkeys, and Przewalski horses. Second, we extracted sequences of lncRNA genes shared in ancient horses and sequences of lncRNA genes and the promoter regions of domestication-critical genes shared in modern horses, modern donkeys, and Przewalski horses to form two sample groups. Third, we used the LongTarget program to predict potential regulatory interactions between these lncRNAs and these domestication-critical genes and analyzed the differences between the regulation in ancient/modern horses and between horses/donkeys/Przewalski horses. Fourth, we performed functional enrichment analyses of genes that exhibit differences in epigenetic regulation.

Results: First, genes associated with neural crest development and domestication syndrome are important targets of lncRNAs. Second, compared with undomesticated Przewalski horses, more lncRNAs participate in the epigenetic regulation in modern horses and donkeys, suggesting that domestication is linked to more epigenetic regulatory changes. Third, lncRNAs’ potential target genes in modern horses are mainly involved in two functional areas: 1) the nervous system, behavior, and cognition, and 2) muscle, body size, cardiac function, and metabolism.

Conclusion: Domestication is linked to substantial epigenetic regulatory changes. Genes associated with neural crest development and domestication syndrome underwent noticeable lncRNA-mediated epigenetic regulation changes during horse domestication.

The genomic history and global expansion of domestic donkeys

Donkeys transformed human history as essential beasts of burden for long-distance movement, especially across semi-arid and upland environments. They remain insufficiently studied despite globally expanding and providing key support to low- to middle-income communities. To elucidate their domestication history, we constructed a comprehensive genome panel of 207 modern and 31 ancient donkeys, as well as 15 wild equids. We found a strong phylogeographic structure in modern donkeys that supports a single domestication in Africa ~5000 BCE, followed by further expansions in this continent and Eurasia and ultimately returning to Africa. We uncover a previously unknown genetic lineage in the Levant ~200 BCE, which contributed increasing ancestry toward Asia. Donkey management involved inbreeding and the production of giant bloodlines at a time when mules were essential to the Roman economy and military.

Detection of common copy number of variation underlying selection pressure in Middle Eastern horse breeds using whole-genome sequence data

Dareshouri, Arabian, and Akhal-Teke are three Middle Eastern horse breeds that have been selected for endurance and adaptation to harsh climates. Deciphering the genetic characteristics of these horses by tracing selection footprints and copy number of variations will be helpful in improving our understanding of equine breeds' development and adaptation. For this purpose, we sequenced the whole-genome of four Dareshouri horses using Illumina Hiseq panels and compared them with publicly available whole-genome sequences of Arabian (n=3) and Akhal-Teke (n=3) horses . Three tests of FLK, hapFLK, and pooled heterozygosity were applied using a sliding window (window size=100kb, step size=50kb) approach to detect putative selection signals. Copy number variation analysis was applied to investigate copy number of variants (CNVs), and the results were used to suggest selection signatures involving CNVs. Whole-genome sequencing demonstrated 8,837,950 single nucleotide polymorphisms (SNPs) in autosomal chromosomes. We suggested 58 genes and three quantitative trait loci (QTLs), including some related to horse gait, insect bite hypersensitivity, and withers height, based on selective signals detected by adjusted p-value of Mahalanobis distance based on the rank-based P-values (Md-rank-P) method. We proposed 12 genomic regions under selection pressure involving CNVs which were previously reported to be associated with metabolism energy (SLC5A8), champagne dilution in horses (SLC36A1), and synthesis of polyunsaturated fatty acids (FAT2). Only 10 Middle Eastern horses were tested in this study; therefore, the conclusions are speculative. Our findings are useful to better understanding the evolution and adaptation of Middle Eastern horse breeds.

Radiocarbon and genomic evidence for the survival of Equus Sussemionus until the late Holocene

The exceptionally rich fossil record available for the equid family has provided textbook examples of macroevolutionary changes. Horses, asses, and zebras represent three extant subgenera of Equus lineage, while the Sussemionus subgenus is another remarkable Equus lineage ranging from North America to Ethiopia in the Pleistocene. We sequenced 26 archaeological specimens from Northern China in the Holocene that could be assigned morphologically and genetically to Equus ovodovi, a species representative of Sussemionus. We present the first high-quality complete genome of the Sussemionus lineage, which was sequenced to 13.4× depth of coverage. Radiocarbon dating demonstrates that this lineage survived until ~3500 years ago, despite continued demographic collapse during the Last Glacial Maximum and the great human expansion in East Asia. We also confirmed the Equus phylogenetic tree and found that Sussemionus diverged from the ancestor of non-caballine equids ~2.3–2.7 million years ago and possibly remained affected by secondary gene flow post-divergence. We found that the small genetic diversity, rather than enhanced inbreeding, limited the species’ chances of survival. Our work adds to the growing literature illustrating how ancient DNA can inform on extinction dynamics and the long-term resilience of species surviving in cryptic population pockets.

Genetic load: genomic estimates and applications in non-model animals

Genetic variation, which is generated by mutation, recombination and gene flow, can reduce the mean fitness of a population, both now and in the future. This 'genetic load' has been estimated in a wide range of animal taxa using various approaches. Advances in genome sequencing and computational techniques now enable us to estimate the genetic load in populations and individuals without direct fitness estimates. Here, we review the classic and contemporary literature of genetic load. We describe approaches to quantify the genetic load in whole-genome sequence data based on evolutionary conservation and annotations. We show that splitting the load into its two components - the realized load (or expressed load) and the masked load (or inbreeding load) - can improve our understanding of the population genetics of deleterious mutations.

Signatures of selection underpinning rapid coral adaptation to the world's warmest reefs

Coral populations in the world’s warmest reefs, the Persian/Arabian Gulf (PAG), represent an ideal model system to understand the evolutionary response of coral populations to past and present environmental change and to identify genomic loci that contribute to elevated thermal tolerance. Here, we use population genomics of the brain coral Platygyra daedalea to show that corals in the PAG represent a distinct subpopulation that was established during the Holocene marine transgression, and identify selective sweeps in their genomes associated with thermal adaptation. We demonstrate the presence of positive and disruptive selection and provide evidence for selection of differentially methylated haplotypes. While demographic analyses suggest limited potential for genetic rescue of neighboring Indian Ocean reefs, the presence of putative targets of selection in corals outside of the PAG offers hope that loci associated with thermal tolerance may be present in the standing genetic variation.

Genetic diversity and selection in Puerto Rican horses

Since the first Spanish settlers brought horses to America centuries ago, several local varieties and breeds have been established in the New World. These were generally a consequence of the admixture of the different breeds arriving from Europe. In some instances, local horses have been selectively bred for specific traits, such as appearance, endurance, strength, and gait. We looked at the genetics of two breeds, the Puerto Rican Non-Purebred (PRNPB) (also known as the "Criollo") horses and the Puerto Rican Paso Fino (PRPF), from the Caribbean Island of Puerto Rico. While it is reasonable to assume that there was a historic connection between the two, the genetic link between them has never been established. In our study, we started by looking at the genetic ancestry and diversity of current Puerto Rican horse populations using a 668 bp fragment of the mitochondrial DNA D-loop (HVR1) in 200 horses from 27 locations on the island. We then genotyped all 200 horses in our sample for the "gait-keeper" DMRT3 mutant allele previously associated with the paso gait especially cherished in this island breed. We also genotyped a subset of 24 samples with the Illumina Neogen Equine Community genome-wide array (65,000 SNPs). This data was further combined with the publicly available PRPF genomes from other studies. Our analysis show an undeniable genetic connection between the two varieties in Puerto Rico, consistent with the hypothesis that PRNPB horses represent the descendants of the original genetic pool, a mix of horses imported from the Iberian Peninsula and elsewhere in Europe. Some of the original founders of PRNRB population must have carried the "gait-keeper" DMRT3 allele upon arrival to the island. From this admixture, the desired traits were selected by the local people over the span of centuries. We propose that the frequency of the mutant "gait-keeper" allele originally increased in the local horses due to the selection for the smooth ride and other characters, long before the PRPF breed was established. To support this hypothesis, we demonstrate that PRNPB horses, and not the purebred PRPF, carry a signature of selection in the genomic region containing the DMRT3 locus to this day. The lack of the detectable signature of selection associated with the DMRT3 in the PRPF would be expected if this native breed was originally derived from the genetic pool of PRNPB horses established earlier and most of the founders already had the mutant allele. Consequently, selection specific to PRPF later focused on allels in other genes (including CHRM5, CYP2E1, MYH7, SRSF1, PAM, PRN and others) that have not been previously associated with the prized paso gait phenotype in Puerto Rico or anywhere else.

Sequence Variant in the TRIM39-RPP21 Gene Readthrough is Shared Across a Cohort of Arabian Foals Diagnosed with Juvenile Idiopathic Epilepsy

Juvenile idiopathic epilepsy (JIE) is a self-limiting neurological disorder with a suspected genetic predisposition affecting young Arabian foals of the Egyptian lineage. The condition is characterized by tonic-clonic seizures with intermittent post-ictal blindness, in which most incidents are sporadic and unrecognized. This study aimed to identify genetic components shared across a local cohort of Arabian foals diagnosed with JIE via a combined whole genome and targeted resequencing approach: Initial whole genome comparisons between a small cohort of nine diagnosed foals (cases) and 27 controls from other horse breeds identified variants uniquely shared amongst the case cohort. Further validation via targeted resequencing of these variants, that pertain to non-intergenic regions, on additional eleven case individuals revealed a single 19bp deletion coupled with a triple-C insertion (Δ19InsCCC) within the TRIM39-RPP21 gene readthrough that was uniquely shared across all case individuals, and absent from three additional Arabian controls. Furthermore, we have confirmed recent findings refuting potential linkage between JIE and other inherited diseases in the Arabian lineage, and refuted the potential linkage between JIE and genes predisposing a similar disorder in human newborns. This is the first study to report a genetic variant to be shared in a sub-population cohort of Arabian foals diagnosed with JIE. Further evaluation of the sensitivity and specificity of the Δ19InsCCC allele within additional cohorts of the Arabian horse is warranted in order to validate its credibility as a marker for JIE, and to ascertain whether it has been introduced into other horse breeds by Arabian ancestry.

Museum Genomics

Natural history collections are invaluable repositories of biological information that provide an unrivaled record of Earth's biodiversity. Museum genomics-genomics research using traditional museum and cryogenic collections and the infrastructure supporting these investigations-has particularly enhanced research in ecology and evolutionary biology, the study of extinct organisms, and the impact of anthropogenic activity on biodiversity. However, leveraging genomics in biological collections has exposed challenges, such as digitizing, integrating, and sharing collections data; updating practices to ensure broadly optimal data extraction from existing and new collections; and modernizing collections practices, infrastructure, and policies to ensure fair, sustainable, and genomically manifold uses of museum collections by increasingly diverse stakeholders. Museum genomics collections are poised to address these challenges and, with increasingly sensitive genomics approaches, will catalyze a future era of reproducibility, innovation, and insight made possible through integrating museum and genome sciences.

Role of genes related to performance and reproduction of Thoroughbreds in training and breeding - A review

Thoroughbreds have been selected for speed and stamina since the 1700s. This selection resulted in structural and functional system-wide adaptations that enhanced physiological characteristics for outstanding speed of 61-71 kph (38-44 mph) between 1,000 and 3,200 m (5 furlongs - 2 miles). At present, horseracing is still an economically important industrial sector, therefore intensive research is underway to explore genes that allow the utilisation of genetic abilities and are significant in breeding and training. This study aims to provide an overview of genetic research and its applicability related to Thoroughbreds.

Deleterious protein-coding variants in diverse cattle breeds of the world

The domestication of wild animals has resulted in a reduction in effective population sizes, which can affect the deleterious mutation load of domesticated breeds. In addition, artificial selection contributes to the accumulation of deleterious mutations because of an increased rate of inbreeding among domesticated animals. Since founder population sizes and artificial selection differ between cattle breeds, their deleterious mutation load can vary. We investigated this question by using whole-genome data from 432 animals belonging to 54 worldwide cattle breeds. Our analysis revealed a negative correlation between genomic heterozygosity and nonsynonymous-to-silent diversity ratio, which suggests a higher proportion of single nucleotide variants (SNVs) affecting proteins in low-diversity breeds. Our results also showed that low-diversity breeds had a larger number of high-frequency (derived allele frequency (DAF) > 0.51) deleterious SNVs than high-diversity breeds. An opposite trend was observed for the low-frequency (DAF ≤ 0.51) deleterious SNVs. Overall, the number of high-frequency deleterious SNVs was larger in the genomes of taurine cattle breeds than of indicine breeds, whereas the number of low-frequency deleterious SNVs was larger in the genomes of indicine cattle than in those of taurine cattle. Furthermore, we observed significant variation in the counts of deleterious SNVs within taurine breeds. The variations in deleterious mutation load between taurine and indicine breeds could be attributed to the population sizes of the wild progenitors before domestication, whereas the variations observed within taurine breeds could be due to differences in inbreeding level, strength of artificial selection, and/or founding population size. Our findings imply that the incidence of genetic diseases can vary between cattle breeds.

Ancient Faunal History Revealed by Interdisciplinary Biomolecular Approaches

Starting four decades ago, studies have examined the ecology and evolutionary dynamics of populations and species using short mitochondrial DNA fragments and stable isotopes. Through technological and analytical advances, the methods and biomolecules at our disposal have increased significantly to now include lipids, whole genomes, proteomes, and even epigenomes. At an unprecedented resolution, the study of ancient biomolecules has made it possible for us to disentangle the complex processes that shaped the ancient faunal diversity across millennia, with the potential to aid in implicating probable causes of species extinction and how humans impacted the genetics and ecology of wild and domestic species. However, even now, few studies explore interdisciplinary biomolecular approaches to reveal ancient faunal diversity dynamics in relation to environmental and anthropogenic impact. This review will approach how biomolecules have been implemented in a broad variety of topics and species, from the extinct Pleistocene megafauna to ancient wild and domestic stocks, as well as how their future use has the potential to offer an enhanced understanding of drivers of past faunal diversity on Earth.

Suggested Absence of Horizontal Transfer of Retrotransposons between Humans and Domestic Mammal Species

Transposable element sequences are usually vertically inherited but have also spread across taxa via horizontal transfer. Previous investigations of ancient horizontal transfer of transposons have compared consensus sequences, but this method resists detection of recent single or low copy number transfer events. The relationship between humans and domesticated animals represents an opportunity for potential horizontal transfer due to the consistent shared proximity and exposure to parasitic insects, which have been identified as plausible transfer vectors. The relatively short period of extended human-animal contact (tens of thousands of years or less) makes horizontal transfer of transposons between them unlikely. However, the availability of high-quality reference genomes allows individual element comparisons to detect low copy number events. Using pairwise all-versus-all megablast searches of the complete suite of retrotransposons of thirteen domestic animals against human, we searched a total of 27,949,823 individual TEs. Based on manual comparisons of stringently filtered BLAST search results for evidence of vertical inheritance, no plausible instances of HTT were identified. These results indicate that significant recent HTT between humans and domesticated animals has not occurred despite the close proximity, either due to the short timescale, inhospitable recipient genomes, a failure of vector activity, or other factors.

Modularity patterns in mammalian domestication: Assessing developmental hypotheses for diversification

The neural crest hypothesis posits that selection for tameness resulted in mild alterations to neural crest cells during embryonic development, which directly or indirectly caused the appearance of traits associated with the "domestication syndrome" (DS). Although representing an appealing unitary explanation for the generation of domestic phenotypes, support for this hypothesis from morphological data and for the validity of the DS remains a topic of debate. This study used the frameworks of morphological integration and modularity to assess patterns that concern the embryonic origin of the skull and issues around the neural crest hypothesis. Geometric morphometric landmarks were used to quantify cranial trait interactions between six pairs of wild and domestic mammals, comprising representatives that express between five and 17 of the traits included in the DS, and examples from each of the pathways by which animals entered into relationships with humans. We predicted the presence of neural crest vs mesoderm modular structure to the cranium, and that elements in the neural crest module would show lower magnitudes of integration and higher disparity in domestic forms compared to wild forms. Our findings support modular structuring based on tissue origin (neural crest, mesoderm) modules, along with low module integration magnitudes for neural crest cell derived cranial elements, suggesting differential capacity for evolutionary response among those elements. Covariation between the neural crest and mesoderm modules accounted for major components of shape variation for most domestic/wild pairs. Contra to our predictions, however, we find domesticates share similar integration magnitudes to their wild progenitors, indicating that higher disparity in domesticates is not associated with magnitude changes to integration among either neural crest or mesoderm derived elements. Differences in integration magnitude among neural crest and mesoderm elements across species suggest that developmental evolution preserves a framework that promotes flexibility under the selection regimes of domestication.

Recent selection of candidate genes for mammal domestication in Europeans and language change in Europe: a hypothesis

BACKGROUND AND AIM

Human evolution resulted from changes in our biology, behaviour, and culture. One source of these changes has been hypothesised to be our self-domestication (that is, the development in humans of features commonly found in domesticated strains of mammals, seemingly as a result of selection for reduced aggression). Signals of domestication, notably brain size reduction, have increased in recent times.

METHODS

In this paper, we compare whole-genome data between the Late Neolithic/Bronze Age individuals and modern Europeans.

RESULTS

We show that genes associated with mammal domestication and with neural crest development and function are significantly differently enriched in nonsynonymous single nucleotide polymorphisms between these two groups.

CONCLUSION

We hypothesise that these changes might account for the increased features of self-domestication in modern humans and, ultimately, for subtle recent changes in human cognition and behaviour, including language.

From the Eurasian Steppes to the Roman Circuses: A Review of Early Development of Horse Breeding and Management

Simple Summary

Horses were domesticated later than any other major livestock species. Their role in shaping ancient civilizations cannot be overestimated. As a primary means of transportation, an essential asset in warfare, and later one of the key elements of circus entertainment, horses quickly became luxurious goods. Vast amounts of money were invested in the horse industry resulted resulting in the rapid development of horse breeding and husbandry. This review examines paleogenetic, archeological, and classical studies on managing horses in antiquity. Many ancient approaches and practices in horse management are still relevant today and some of them, now abandoned, are worth re-examination.

Abstract

The domestication of the horse began about 5500 years ago in the Eurasian steppes. In the following millennia horses spread across the ancient world, and their role in transportation and warfare affected every ancient culture. Ownership of horses became an indicator of wealth and social status. The importance of horses led to a growing interest in their breeding and management. Many phenotypic traits, such as height, behavior, and speed potential, have been proven to be a subject of selection; however, the details of ancient breeding practices remain mostly unknown. From the fourth millennium BP, through the Iron Age, many literature sources thoroughly describe horse training systems, as well as various aspects of husbandry, many of which are still in use today. The striking resemblance of ancient and modern equine practices leaves us wondering how much was accomplished through four thousand years of horse breeding.

Large-scale genomic analysis reveals the genetic cost of chicken domestication

BACKGROUND

Species domestication is generally characterized by the exploitation of high-impact mutations through processes that involve complex shifting demographics of domesticated species. These include not only inbreeding and artificial selection that may lead to the emergence of evolutionary bottlenecks, but also post-divergence gene flow and introgression. Although domestication potentially affects the occurrence of both desired and undesired mutations, the way wild relatives of domesticated species evolve and how expensive the genetic cost underlying domestication is remain poorly understood. Here, we investigated the demographic history and genetic load of chicken domestication.

RESULTS

We analyzed a dataset comprising over 800 whole genomes from both indigenous chickens and wild jungle fowls. We show that despite having a higher genetic diversity than their wild counterparts (average π, 0.00326 vs. 0.00316), the red jungle fowls, the present-day domestic chickens experienced a dramatic population size decline during their early domestication. Our analyses suggest that the concomitant bottleneck induced 2.95% more deleterious mutations across chicken genomes compared with red jungle fowls, supporting the "cost of domestication" hypothesis. Particularly, we find that 62.4% of deleterious SNPs in domestic chickens are maintained in heterozygous states and masked as recessive alleles, challenging the power of modern breeding programs to effectively eliminate these genetic loads. Finally, we suggest that positive selection decreases the incidence but increases the frequency of deleterious SNPs in domestic chicken genomes.

CONCLUSION

This study reveals a new landscape of demographic history and genomic changes associated with chicken domestication and provides insight into the evolutionary genomic profiles of domesticated animals managed under modern human selection.

Genomics of Adaptations in Ungulates

Ungulates are a group of hoofed animals that have long interacted with humans as essential sources of food, labor, clothing, and transportation. These consist of domesticated, feral, and wild species raised in a wide range of habitats and biomes. Given the diverse and extreme environments inhabited by ungulates, unique adaptive traits are fundamental for fitness. The documentation of genes that underlie their genomic signatures of selection is crucial in this regard. The increasing availability of advanced sequencing technologies has seen the rapid growth of ungulate genomic resources, which offers an exceptional opportunity to understand their adaptive evolution. Here, we summarize the current knowledge on evolutionary genetic signatures underlying the adaptations of ungulates to different habitats.

An overview of current population genomics methods for the analysis of whole-genome resequencing data in eukaryotes

Characterizing the population history of a species and identifying loci underlying local adaptation is crucial in functional ecology, evolutionary biology, conservation and agronomy. The constant improvement of high-throughput sequencing techniques has facilitated the production of whole genome data in a wide range of species. Population genomics now provides tools to better integrate selection into a historical framework, and take into account selection when reconstructing demographic history. However, this improvement has come with a profusion of analytical tools that can confuse and discourage users. Such confusion limits the amount of information effectively retrieved from complex genomic data sets, and impairs the diffusion of the most recent analytical tools into fields such as conservation biology. It may also lead to redundancy among methods. To address these isssues, we propose an overview of more than 100 state-of-the-art methods that can deal with whole genome data. We summarize the strategies they use to infer demographic history and selection, and discuss some of their limitations. A website listing these methods is available at www.methodspopgen.com.

Ancient horse genomes reveal the timing and extent of dispersals across the Bering Land Bridge

The Bering Land Bridge (BLB) last connected Eurasia and North America during the Late Pleistocene. Although the BLB would have enabled transfers of terrestrial biota in both directions, it also acted as an ecological filter whose permeability varied considerably over time. Here we explore the possible impacts of this ecological corridor on genetic diversity within, and connectivity among, populations of a once wide-ranging group, the caballine horses (Equus spp.). Using a panel of 187 mitochondrial and eight nuclear genomes recovered from present-day and extinct caballine horses sampled across the Holarctic, we found that Eurasian horse populations initially diverged from those in North America, their ancestral continent, around 1.0-0.8 million years ago. Subsequent to this split our mitochondrial DNA analysis identified two bidirectional long-range dispersals across the BLB ~875-625 and ~200-50 thousand years ago, during the Middle and Late Pleistocene. Whole genome analysis indicated low levels of gene flow between North American and Eurasian horse populations, which probably occurred as a result of these inferred dispersals. Nonetheless, mitochondrial and nuclear diversity of caballine horse populations retained strong phylogeographical structuring. Our results suggest that barriers to gene flow, currently unidentified but possibly related to habitat distribution across Beringia or ongoing evolutionary divergence, played an important role in shaping the early genetic history of caballine horses, including the ancestors of living horses within Equus ferus.

Steller's sea cow genome suggests this species began going extinct before the arrival of Paleolithic humans

Anthropogenic activity is the top factor directly related to the extinction of several animal species. The last Steller's sea cow (Hydrodamalis gigas) population on the Commander Islands (Russia) was wiped out in the second half of the 18^th century due to sailors and fur traders hunting it for the meat and fat. However, new data suggests that the extinction process of this species began much earlier. Here, we present a nuclear de novo assembled genome of H. gigas with a 25.4× depth coverage. Our results demonstrate that the heterozygosity of the last population of this animal is low and comparable to the last woolly mammoth population that inhabited Wrangel Island 4000 years ago. Besides, as a matter of consideration, our findings also demonstrate that the extinction of this marine mammal starts along the North Pacific coastal line much earlier than the first Paleolithic humans arrived in the Bering sea region.

Genomics in the Horse Industry: Discovering New Questions at Every Turn

The sheer diversity of heritable physiological traits, and the ingenuity of genome derived research technologies, extends the study of genetics to impact diverse scientific fields. Equine science is no exception, experiencing a number of genome-enabled discoveries that spur further research in areas like nutrition, reproduction, and exercise physiology. Yet unexpected findings, especially those that over-turn commonly held beliefs in the horse industry, can create challenges in outreach, education and communication with stakeholders. For example, studies of ancient DNA revealed that the oldest domesticated equids in the archeological record were in fact another species, the Przewalski's horse, leaving the origins of our modern horses a mystery yet to be solved. Genomic analysis of ancestry can illuminate relationships older than our prized pedigree records, and in some cases, identify unexpected inconsistencies in those pedigrees. Even our interpretation of what constitutes a genetic disease is changing, as we re-examine common disease alleles; how these alleles impact equine physiology, and how they are perceived by breeders and professionals in the industry. Effectively translating genetic tools for utilization in horse management and preparing our community for the debate surrounding ethical questions that may arise from genomic studies, may be the next great challenges we face as scientists and educators.

Investigating the population structure and genetic diversity of Arabian horses in Oman using SNP markers

Arabian horses were selected for metabolic efficiency, beauty, efficiency and endurance. Therefore, Bedouins have for centuries traced their prized horses' ancestries. With the establishment of the World Arabian Horse Organization (WAHO), registration of Arabian horses became centralized and countries worldwide registered them in its database. Most existing Arabian horses in Oman today were imported after the 1970s and are predominantly flat-racing Arabians. This work aimed at revealing the genetic background and diversity of Omani Arabian horses by comparing them with Arabian horses from a diverse genetic background. To that end, we genotyped 63 randomly sampled Arabian horses from Oman using the Illumina Equine SNP70. For comparison, SNP genotypes of 12 Saudi Arabian horses, 27 French, 77 Egyptian, 11 Polish and 36 US Arabians were included in the study. We additionally included 17 Thoroughbred horses and 21 horses representing large and small breeds as an outgroup. Our MDS analysis and phylogenetic analysis showed that the Arabian horses in Oman cluster primarily with French Arabian horses, with a few horses clustering within the Polish/US Arabians. The French Arabian horse cluster was the closest to the Thoroughbred horses. Amongst the Arabian horses, plink average genomic inbreeding levels were highest in the Egyptian Arabian (0.169) followed by the Saudi Arabian horses (0.137) and lowest in the Omani and French Arabian horses, -0.041 and -0.079 respectively. To our knowledge, this is the first report on the genetic background and diversity of Arabian horses in Oman. Our results demonstrated a definite subpopulation structure among Arabian horses and this information should advise future decision-making on Arabian horse breeding.

Ancient DNA reveals the lost domestication history of South American camelids in Northern Chile and across the Andes

The study of South American camelids and their domestication is a highly debated topic in zooarchaeology. Identifying the domestic species (alpaca and llama) in archaeological sites based solely on morphological data is challenging due to their similarity with respect to their wild ancestors. Using genetic methods also presents challenges due to the hybridization history of the domestic species, which are thought to have extensively hybridized following the Spanish conquest of South America that resulted in camelids slaughtered en masse. In this study, we generated mitochondrial genomes for 61 ancient South American camelids dated between 3,500 and 2,400 years before the present (Early Formative period) from two archaeological sites in Northern Chile (Tulán-54 and Tulán-85), as well as 66 modern camelid mitogenomes and 815 modern mitochondrial control region sequences from across South America. In addition, we performed osteometric analyses to differentiate big and small body size camelids. A comparative analysis of these data suggests that a substantial proportion of the ancient vicuña genetic variation has been lost since the Early Formative period, as it is not present in modern specimens. Moreover, we propose a domestication hypothesis that includes an ancient guanaco population that no longer exists. Finally, we find evidence that interbreeding practices were widespread during the domestication process by the early camelid herders in the Atacama during the Early Formative period and predating the Spanish conquest.

Traces of Late Bronze and Early Iron Age Mongolian Horse Mitochondrial Lineages in Modern Populations

The Mongolian horse is one of the most ancient and relatively unmanaged horse breeds. The population history of the Mongolian horse remains poorly understood due to a lack of information on ancient and modern DNA. Here, we report nearly complete mitochondrial genome data obtained from five ancient Mongolian horse samples of the Khereksur and Deer Stone culture (late 2nd to 1st third of the 1st millennium BC) and one ancient horse specimen from the Xiongnu culture (1st century BC to 1st century AD) using target enrichment and high-throughput sequencing methods. Phylogenetic analysis involving ancient, historical, and modern mitogenomes of horses from Mongolia and other regions showed the presence of three mitochondrial haplogroups in the ancient Mongolian horse populations studied here and similar haplotype composition of ancient and modern horse populations of Mongolia. Our results revealed genetic continuity between the Mongolian horse populations of the Khereksur and Deer Stone culture and those of the Xiongnu culture owing to the presence of related mitotypes. Besides, we report close phylogenetic relationships between haplotypes of the Khereksur and Deer Stone horses and the horses of indigenous breeds of the Middle East (Caspian and Iranian), China (Naqu, Yunnan, and Jinjiang), and Italy (Giara) as well as genetic similarity between the Xiongnu Mongolian horses and those of the most ancient breeds of the Middle East (Arabian) and Central Asia (Akhal-Teke). Despite all the migrations of the Mongolian peoples over the past 3000 years, mitochondrial haplogroup composition of Mongolian horse populations remains almost unchanged.

Genome diversity in Ukraine

Background

The main goal of this collaborative effort is to provide genome-wide data for the previously underrepresented population in Eastern Europe, and to provide cross-validation of the data from genome sequences and genotypes of the same individuals acquired by different technologies. We collected 97 genome-grade DNA samples from consented individuals representing major regions of Ukraine that were consented for public data release. BGISEQ-500 sequence data and genotypes by an Illumina GWAS chip were cross-validated on multiple samples and additionally referenced to 1 sample that has been resequenced by Illumina NovaSeq6000 S4 at high coverage.

Results

The genome data have been searched for genomic variation represented in this population, and a number of variants have been reported: large structural variants, indels, copy number variations, single-nucletide polymorphisms, and microsatellites. To our knowledge, this study provides the largest to-date survey of genetic variation in Ukraine, creating a public reference resource aiming to provide data for medical research in a large understudied population.

Conclusions

Our results indicate that the genetic diversity of the Ukrainian population is uniquely shaped by evolutionary and demographic forces and cannot be ignored in future genetic and biomedical studies. These data will contribute a wealth of new information bringing forth a wealth of novel, endemic and medically related alleles.

The origins and spread of domestic horses from the Western Eurasian steppes

Domestication of horses fundamentally transformed long-range mobility and warfare¹. However, modern domesticated breeds do not descend from the earliest domestic horse lineage associated with archaeological evidence of bridling, milking and corralling^2-4 at Botai, Central Asia around 3500 BC³. Other longstanding candidate regions for horse domestication, such as Iberia⁵ and Anatolia⁶, have also recently been challenged. Thus, the genetic, geographic and temporal origins of modern domestic horses have remained unknown. Here we pinpoint the Western Eurasian steppes, especially the lower Volga-Don region, as the homeland of modern domestic horses. Furthermore, we map the population changes accompanying domestication from 273 ancient horse genomes. This reveals that modern domestic horses ultimately replaced almost all other local populations as they expanded rapidly across Eurasia from about 2000 BC, synchronously with equestrian material culture, including Sintashta spoke-wheeled chariots. We find that equestrianism involved strong selection for critical locomotor and behavioural adaptations at the GSDMC and ZFPM1 genes. Our results reject the commonly held association⁷ between horseback riding and the massive expansion of Yamnaya steppe pastoralists into Europe around 3000 BC^8,9 driving the spread of Indo-European languages¹⁰. This contrasts with the scenario in Asia where Indo-Iranian languages, chariots and horses spread together, following the early second millennium BC Sintashta culture^11,12.

Genetic Mapping of Climbing and Mimicry: Two Behavioral Traits Degraded During Silkworm Domestication

Behavioral changes caused by domestication in animals are an important issue in evolutionary biology. The silkworm, Bombyx mori, is an ideal fully domesticated insect model for studying both convergent domestication and behavior evolution. We explored the genetic basis of climbing for foraging and mimicry, two degraded behaviors during silkworm domestication, in combination of bulked segregant analysis (BSA) and selection sweep screening. One candidate gene, ASNA1, located in the 3-5 Mb on chromosome 19, harboring a specific non-synonymous mutation in domestic silkworm, might be involved in climbing ability. This mutation was under positive selection in Lepidoptera, strongly suggesting its potential function in silkworm domestication. Nine candidate domesticated genes related to mimicry were identified on chromosomes 13, 21, and 27. Most of the candidate domesticated genes were generally expressed at higher levels in the brain of the wild silkworm. This study provides valuable information for deciphering the molecular basis of behavioral changes associated with silkworm domestication.

Variability of ACOX1 Gene Polymorphisms across Different Horse Breeds with Regard to Selection Pressure

Simple Summary

The genetic mechanisms occurring in organisms are shaped by selection pressure. Features that ought to be useful under given conditions leave their marks on the genome in the form of mutations, thereby creating different alleles. In this study, five different horse breeds were examined to find the connection between an individual’s lifestyle and the presence of the peroxisomal acyl-coenzyme A oxidase 1 (ACOX1) gene, which is necessary for some metabolic pathways. Results indicated that different ACOX1 gene alleles play various roles in primitive breeds and domesticated horses. This led to the conclusion that the DNA profile can be rated on the basis of adaptation to living conditions, opening the gate for further investigation.

Abstract

The ACOX1 gene encodes peroxisomal acyl-coenzyme A oxidase 1, the first enzyme in the fatty acid β-oxidation pathway, which could be significant for organisms exposed to long periods of starvation and harsh living conditions. We hypothesized that variations within ACOX1, revealed by RNA Sequencing (RNA-Seq), might be based on adaptation to living conditions and had resulted from selection pressure. There were five different horse breeds used in this study, representing various utility types: Arabian, Thoroughbred, Polish Konik, draft horses, and Hucul. The single-nucleotide polymorphism (SNP) located in the ACOX1 (rs782885985) was used as a marker and was identified using the PCR restriction fragment length polymorphism method (PCR-RFLP). Results indicated extremely different genotype and allele distributions of the ACOX1 gene across breeds. A predominance of the G allele was exhibited in horses that had adapted to difficult environmental conditions, namely, Polish Konik and Huculs, which are considered to be primitive breeds. The prevalence of the T allele in Thoroughbreds indicated that ACOX1 is significant in energy metabolism during flat racing.

Selection in the Finnhorse, a native all-around horse breed

Selection by breeders modifies the morphology, behaviour and performance of domesticated species. Here, we examined signs of selection in Finnhorse, the only native horse breed in Finland. We first searched divergent genomic regions between Finnhorses and other breeds, as well as between different breeding sections of the Finnhorse with data from Illumina Equine SNP70 BeadChip, and then studied several of the detected regions in more detail. We found altogether 35 common outlier SNPs between Finnhorses and other breeds using two different selection tests. Many of the SNPs were located close to genes affecting coat colour, performance, size, sugar metabolism, immune response and olfaction. We selected genes affecting coat colour (KIT, MITF, PMEL), performance (MSTN) and locomotion (DMRT3) for a more detailed examination. In addition, we looked for, and found, associations with height at withers and SNPs located close to gene LCORL. Among the four breeding sections of Finnhorses (harness trotters, riding horses, draught horses and pony‐sized horses), a single SNP located close to the DMRT3 gene was significantly differentiated and only between harness trotters and pony‐sized horses.