Genome sequence, comparative analysis, and population genetics of the domestic horse

Origin and spread of Thoroughbred racehorses inferred from complete mitochondrial genome sequences: Phylogenomic and Bayesian coalescent perspectives

The Thoroughbred horse breed was developed primarily for racing, and has a significant contribution to the qualitative improvement of many other horse breeds. Despite the importance of Thoroughbred racehorses in historical, cultural, and economical viewpoints, there was no temporal and spatial dynamics of them using the mitogenome sequences. To explore this topic, the complete mitochondrial genome sequences of 14 Thoroughbreds and two Przewalski’s horses were determined. These sequences were analyzed together along with 151 previously published horse mitochondrial genomes from a range of breeds across the globe using a Bayesian coalescent approach as well as Bayesian inference and maximum likelihood methods. The racing horses were revealed to have multiple maternal origins and to be closely related to horses from one Asian, two Middle Eastern, and five European breeds. Thoroughbred horse breed was not directly related to the Przewalski’s horse which has been regarded as the closest taxon to the all domestic horses and the only true wild horse species left in the world. Our phylogenomic analyses also supported that there was no apparent correlation between geographic origin or breed and the evolution of global horses. The most recent common ancestor of the Thoroughbreds lived approximately 8,100–111,500 years ago, which was significantly younger than the most recent common ancestor of modern horses (0.7286 My). Bayesian skyline plot revealed that the population expansion of modern horses, including Thoroughbreds, occurred approximately 5,500–11,000 years ago, which coincide with the start of domestication. This is the first phylogenomic study on the Thoroughbred racehorse in association with its spatio-temporal dynamics. The database and genetic history information of Thoroughbred mitogenomes obtained from the present study provide useful information for future horse improvement projects, as well as for the study of horse genomics, conservation, and in association with its geographical distribution.

Transcriptome analysis of immune genes in peripheral blood mononuclear cells of young foals and adult horses

During the neonatal period, the ability to generate immune effector and memory responses to vaccines or pathogens is often questioned. This study was undertaken to obtain a global view of the natural differences in the expression of immune genes early in life. Our hypothesis was that transcriptome analyses of peripheral blood mononuclear cells (PBMCs) of foals (on day 1 and day 42 after birth) and adult horses would show differential gene expression profiles that characterize natural immune processes. Gene ontology enrichment analysis provided assessment of biological processes affected by age, and a list of 897 genes with ≥2 fold higher (p<0.01) expression in day 42 when compared to day 1 foal samples. Up-regulated genes included B cell and T cell receptor diversity genes; DNA replication enzymes; natural killer cell receptors; granzyme B and perforin; complement receptors; immunomodulatory receptors; cell adhesion molecules; and cytokines/chemokines and their receptors. The list of 1,383 genes that had higher (p<0.01) expression on day 1 when compared to day 42 foal samples was populated by genes with roles in innate immunity such as antimicrobial proteins; pathogen recognition receptors; cytokines/chemokines and their receptors; cell adhesion molecules; co-stimulatory molecules; and T cell receptor delta chain. Within the 742 genes with increased expression between day 42 foal and adult samples, B cell immunity was the main biological process (p = 2.4E-04). Novel data on markedly low (p<0.0001) TLR3 gene expression, and high (p≤0.01) expression of IL27, IL13RA1, IREM-1, SIRL-1, and SIRPα on day 1 compared to day 42 foal samples point out potential mechanisms of increased susceptibility to pathogens in early life. The results portray a progression from innate immune gene expression predominance early in life to adaptive immune gene expression increasing with age with a putative overlay of immune suppressing genes in the neonatal phase. These results provide insight to the unique attributes of the equine neonatal and young immune system, and offer many avenues of future investigation.

Genome-wide analyses of the Jeju, Thoroughbred, and Jeju crossbred horse populations using the high density SNP array

The Jeju horse is an indigenous Korean horse breed that is currently registered with the Food and Agriculture Organization of the United Nations. However, there is severe lack of genomic studies on Jeju horse. This study was conducted to investigate genetic characteristics of horses including Jeju horse, Thoroughbred and Jeju crossbred (Jeju × Thoroughbred) populations. We compared the genomes of three horse populations using the Equine SNP70 Beadchip array. Short-range Linkage disequilibrium was the highest in Thoroughbred, whereas r² values were lowest in Jeju horse. Expected heterozygosity was the highest in Jeju crossbred (0.351), followed by the Thoroughbred (0.337) and Jeju horse (0.311). The level of inbreeding was slightly higher in Thoroughbred (- 0.009) than in Jeju crossbred (- 0.035) and Jeju horse (- 0.038). F_ST value was the highest between Jeju horse and Thoroughbred (0.113), whereas Jeju crossbred and Thoroughbred showed the lowest value (0.031). The genetic relationship was further assessed by principal component analysis, suggesting that Jeju crossbred is more genetically similar to Thoroughbred than Jeju horse population. Additionally, we detected potential selection signatures, for example, in loci located on LCORL/NCAPG and PROP1 genes that are known to influence body. Genome-wide analyses of the three horse populations showed that all the breeds had somewhat a low level of inbreeding within each population. In the population structure analysis, we found that Jeju crossbred was genetically closer to Thoroughbred than Jeju horse. Furthermore, we identified several signatures of selection which might be associated with traits of interest. To our current knowledge, this study is the first genomic research, analyzing genetic relationships of Jeju horse, Thoroughbred and Jeju crossbred.

Comparison of Oryza sativa and Oryza brachyantha Genomes Reveals Selection-Driven Gene Escape from the Centromeric Regions

Centromeres are dynamic chromosomal regions, and the genetic and epigenetic environment of the centromere is often regarded as oppressive to protein-coding genes. Here, we used comparative genomic and phylogenomic approaches to study the evolution of centromeres and centromere-linked genes in the genus Oryza We report a 12.4-Mb high-quality BAC-based pericentromeric assembly for Oryza brachyantha, which diverged from cultivated rice (Oryza sativa) ∼15 million years ago. The synteny analyses reveal seven medium (>50 kb) pericentric inversions in O. sativa and 10 in O. brachyantha Of these inversions, three resulted in centromere movement (Chr1, Chr7, and Chr9). Additionally, we identified a potential centromere-repositioning event, in which the ancestral centromere on chromosome 12 in O. brachyantha jumped ∼400 kb away, possibly mediated by a duplicated transposition event (>28 kb). More strikingly, we observed an excess of syntenic gene loss at and near the centromeric regions (P < 2.2 × 10^-16). Most (33/47) of the missing genes moved to other genomic regions; therefore such excess could be explained by the selective loss of the copy in or near centromeric regions after gene duplication. The pattern of gene loss immediately adjacent to centromeric regions suggests centromere chromatin dynamics (e.g., spreading or microrepositioning) may drive such gene loss.

Horse Y chromosome assembly displays unique evolutionary features and putative stallion fertility genes

Dynamic evolutionary processes and complex structure make the Y chromosome among the most diverse and least understood regions in mammalian genomes. Here, we present an annotated assembly of the male specific region of the horse Y chromosome (eMSY), representing the first comprehensive Y assembly in odd-toed ungulates. The eMSY comprises single-copy, equine specific multi-copy, PAR transposed, and novel ampliconic sequence classes. The eMSY gene density approaches that of autosomes with the highest number of retained X-Y gametologs recorded in eutherians, in addition to novel Y-born and transposed genes. Horse, donkey and mule testis RNAseq reveals several candidate genes for stallion fertility. A novel testis-expressed XY ampliconic sequence class, ETSTY7, is shared with the parasite Parascaris genome, providing evidence for eukaryotic horizontal transfer and inter-chromosomal mobility. Our study highlights the dynamic nature of the Y and provides a reference sequence for improved understanding of equine male development and fertility.

Understanding eukaryotic chromosome segregation from a comparative biology perspective

A long-appreciated variation in fundamental cell biological processes between different species is becoming increasingly tractable due to recent breakthroughs in whole-genome analyses and genome editing techniques. However, the bulk of our mechanistic understanding in cell biology continues to come from just a few well-established models. In this Review, I use the highly diverse strategies of chromosome segregation in eukaryotes as an instrument for a more general discussion on phenotypic variation, possible rules underlying its emergence and its utility in understanding conserved functional relationships underlying this process. Such a comparative approach, supported by modern molecular biology tools, might provide a wider, holistic view of biology that is difficult to achieve when concentrating on a single experimental system.

Identification of genetic variation in equine collagenous lectins using targeted resequencing

Collagenous lectins are a family of soluble pattern recognition receptors that play an important role in innate immune resistance to infectious disease. Through recognition of carbohydrate motifs on the surface of pathogens, some collagenous lectins can activate the lectin pathway of complement, providing an effective means of host defense. Genetic polymorphisms in collagenous lectins have been shown in several species to predispose animals to a variety of infectious diseases. Infectious diseases are an important cause of morbidity in horses, however little is known regarding the role of equine collagenous lectins. Using a high-throughput, targeted re-sequencing approach, the relationship between genetic variation in equine collagenous lectin genes and susceptibility to disease was investigated. DNA was isolated from tissues obtained from horses submitted for post-mortem examination. Animals were divided into two populations, those with infectious or autoinflammatory diseases (n = 37) and those without (n = 52), and then subdivided by dominant pathological process for a total of 21 pools, each containing 4-5 horses. DNA was extracted from each horse and pooled in equimolar amounts, and the exons, introns, upstream (approximately 50 kb) and downstream (approximately 3 kb) regulatory regions for the 11 equine collagenous lectin genes and related MASP genes were targeted for re-sequencing. A custom target capture kit was used to prepare a sequencing library, which was sequenced on an Illumina MiSeq. After implementing quality control filters, 4559 variants were identified. Of these, 92 were present in the coding regions (43 missense, 1 nonsense, and 48 synonymous), 1414 in introns, 3029 in the upstream region, and 240 in the downstream region. In silico analysis of the missense short nucleotide variants identified 12 mutations with potential to disrupt collagenous lectin protein structure or function, 280 mutations located within predicted transcription factor binding sites, and 95 mutations located within predicted microRNA binding elements. Analysis of allelic association identified 113 mutations that segregated between the infectious/autoinflammatory and non-infectious populations. The variants discovered in this experiment represent potential genetic contributors to disease susceptibility of horses, and will serve as candidates for further population-level genotyping. This study contributes to the growing body of evidence that pooled, high-throughput sequencing is a viable strategy for cost-effective variant discovery.

Detecting the Population Structure and Scanning for Signatures of Selection in Horses (Equus caballus) From Whole-Genome Sequencing Data

Animal domestication gives rise to gradual changes at the genomic level through selection in populations. Selective sweeps have been traced in the genomes of many animal species, including humans, cattle, and dogs. However, little is known regarding positional candidate genes and genomic regions that exhibit signatures of selection in domestic horses. In addition, an understanding of the genetic processes underlying horse domestication, especially the origin of Chinese native populations, is still lacking. In our study, we generated whole genome sequences from 4 Chinese native horses and combined them with 48 publicly available full genome sequences, from which 15 341 213 high-quality unique single-nucleotide polymorphism variants were identified. Kazakh and Lichuan horses are 2 typical Asian native breeds that were formed in Kazakh or Northwest China and South China, respectively. We detected 1390 loss-of-function (LoF) variants in protein-coding genes, and gene ontology (GO) enrichment analysis revealed that some LoF-affected genes were overrepresented in GO terms related to the immune response. Bayesian clustering, distance analysis, and principal component analysis demonstrated that the population structure of these breeds largely reflected weak geographic patterns. Kazakh and Lichuan horses were assigned to the same lineage with other Asian native breeds, in agreement with previous studies on the genetic origin of Chinese domestic horses. We applied the composite likelihood ratio method to scan for genomic regions showing signals of recent selection in the horse genome. A total of 1052 genomic windows of 10 kB, corresponding to 933 distinct core regions, significantly exceeded neutral simulations. The GO enrichment analysis revealed that the genes under selective sweeps were overrepresented with GO terms, including “negative regulation of canonical Wnt signaling pathway,” “muscle contraction,” and “axon guidance.” Frequent exercise training in domestic horses may have resulted in changes in the expression of genes related to metabolism, muscle structure, and the nervous system.

A potential regulatory region near the EDN3 gene may control both harness racing performance and coat color variation in horses

The Swedish-Norwegian Coldblooded trotter and the heavier North-Swedish draught horse both descend from the North-Swedish horse, but the Coldblooded trotters have been selected for racing performance while the North-Swedish draught horse is mainly used for agricultural and forestry work. By comparing the genomes of Coldblooded trotters, North-Swedish draught horses and Standardbreds for a large number of single-nucleotide polymorphisms (SNPs), the aim of the study was to identify genetic regions that may be under selection for racing performance. We hypothesized that the selection for racing performance, in combination with unauthorized crossbreeding of Coldblooded trotters and Standardbreds, has created regions in the genome where the Coldblooded trotters and Standardbreds are similar, but differ from the North-Swedish draught horse. A fixation index (Fst) analysis was performed and sliding window Delta Fst values were calculated across the three breeds. Five windows, where the average Fst between Coldblooded trotters and Standardbreds was low and the average Fst between Coldblooded trotters and North-Swedish draught horses was high, were selected for further investigation. Associations between the most highly ranked SNPs and harness racing performance were analyzed in 400 raced Coldblooded trotters with race records. One SNP showed a significant association with racing performance, with the CC genotype appearing to be negatively associated. The SNP identified was genotyped in 1915 horses of 18 different breeds. The frequency of the TT genotype was high in breeds typically used for racing and show jumping while the frequency of the CC genotype was high in most pony breeds and draught horses. The closest gene in this region was the Endothelin3 gene (EDN3), a gene mainly involved in melanocyte and enteric neuron development. Both functional genetic and physiological studies are needed to fully understand the possible impacts of the gene on racing performance.

Birth, evolution, and transmission of satellite-free mammalian centromeric domains

Mammalian centromeres are associated with highly repetitive DNA (satellite DNA), which has so far hindered molecular analysis of this chromatin domain. Centromeres are epigenetically specified, and binding of the CENPA protein is their main determinant. In previous work, we described the first example of a natural satellite-free centromere on Equus caballus Chromosome 11. Here, we investigated the satellite-free centromeres of Equus asinus by using ChIP-seq with anti-CENPA antibodies. We identified an extraordinarily high number of centromeres lacking satellite DNA (16 of 31). All of them lay in LINE- and AT-rich regions. A subset of these centromeres is associated with DNA amplification. The location of CENPA binding domains can vary in different individuals, giving rise to epialleles. The analysis of epiallele transmission in hybrids (three mules and one hinny) showed that centromeric domains are inherited as Mendelian traits, but their position can slide in one generation. Conversely, centromere location is stable during mitotic propagation of cultured cells. Our results demonstrate that the presence of more than half of centromeres void of satellite DNA is compatible with genome stability and species survival. The presence of amplified DNA at some centromeres suggests that these arrays may represent an intermediate stage toward satellite DNA formation during evolution. The fact that CENPA binding domains can move within relatively restricted regions (a few hundred kilobases) suggests that the centromeric function is physically limited by epigenetic boundaries.

Downregulation of MicroRNA eca-mir-128 in Seminal Exosomes and Enhanced Expression of CXCL16 in the Stallion Reproductive Tract Are Associated with Long-Term Persistence of Equine Arteritis Virus

Equine arteritis virus (EAV) can establish long-term persistent infection in the reproductive tract of stallions and is shed in the semen. Previous studies showed that long-term persistence is associated with a specific allele of the CXCL16 gene (CXCL16S) and that persistent infection is maintained despite the presence of a local inflammatory and humoral and mucosal antibody responses. In this study, we demonstrated that equine seminal exosomes (SEs) are enriched in a small subset of microRNAs (miRNAs). Most importantly, we demonstrated that long-term EAV persistence is associated with the downregulation of an SE-associated miRNA (eca-mir-128) and with an enhanced expression of CXCL16 in the reproductive tract, a putative target of eca-mir-128. The findings presented here suggest that SE eca-mir-128 is implicated in the regulation of the CXCL16/CXCR6 axis in the reproductive tract of persistently infected stallions, a chemokine axis strongly implicated in EAV persistence. This is a novel finding and warrants further investigation to identify its specific mechanism in modulating the CXCL16/CXCR6 axis in the reproductive tract of the EAV long-term carrier stallion.IMPORTANCE Equine arteritis virus (EAV) has the ability to establish long-term persistent infection in the stallion reproductive tract and to be shed in semen, which jeopardizes its worldwide control. Currently, the molecular mechanisms of viral persistence are being unraveled, and these are essential for the development of effective therapeutics to eliminate persistent infection. Recently, it has been determined that long-term persistence is associated with a specific allele of the CXCL16 gene (CXCL16S) and is maintained despite induction of local inflammatory, humoral, and mucosal antibody responses. This study demonstrated that long-term persistence is associated with the downregulation of seminal exosome miRNA eca-mir-128 and enhanced expression of its putative target, CXCL16, in the reproductive tract. For the first time, this study suggests complex interactions between eca-mir-128 and cellular elements at the site of EAV persistence and implicates this miRNA in the regulation of the CXCL16/CXCR6 axis in the reproductive tract during long-term persistence.

Transformation of animal genomics by next-generation sequencing technologies: a decade of challenges and their impact on genetic architecture

For more than a quarter of a century, sequencing technologies from Sanger's method to next-generation high-throughput techniques have provided fascinating opportunities in the life sciences. The continuing upward trajectory of sequencing technologies will improve livestock research and expedite the development of various new genomic and technological studies with farm animals. The use of high-throughput technologies in livestock research has increased interest in metagenomics, epigenetics, genome-wide association studies, and identification of single nucleotide polymorphisms and copy number variations. Such studies are beginning to provide revolutionary insights into biological and evolutionary processes. Farm animals, such as cattle, swine, and horses, have played a dual role as economically and agriculturally important animals as well as biomedical research models. The first part of this study explores the current state of sequencing methods, many of which are already used in animal genomic studies, and the second part summarizes the state of cattle, swine, horse, and chicken genome sequencing and illustrates its achievements during the last few years. Finally, we describe several high-throughput sequencing approaches for the improved detection of known, unknown, and emerging infectious agents, leading to better diagnosis of infectious diseases. The insights from viral metagenomics and the advancement of next-generation sequencing will strongly support specific and efficient vaccine development and provide strategies for controlling infectious disease transmission among animal populations and/or between animals and humans. However, prospective sequencing technologies will require further research and in-field testing before reaching the marketplace.

Function of Junk: Pericentromeric Satellite DNA in Chromosome Maintenance

Satellite DNAs are simple tandem repeats that exist at centromeric and pericentromeric regions on eukaryotic chromosomes. Unlike the centromeric satellite DNA that comprises the vast majority of natural centromeres, function(s) for the much more abundant pericentromeric satellite repeats are poorly understood. In fact, the lack of coding potential allied with rapid divergence of repeat sequences across eukaryotes has led to their dismissal as "junk DNA" or "selfish parasites." Although implicated in various biological processes, a conserved function for pericentromeric satellite DNA remains unidentified. We have addressed the role of satellite DNA through studying chromocenters, a cytological aggregation of pericentromeric satellite DNA from multiple chromosomes into DNA-dense nuclear foci. We have shown that multivalent satellite DNA-binding proteins cross-link pericentromeric satellite DNA on chromosomes into chromocenters. Disruption of chromocenters results in the formation of micronuclei, which arise by budding off the nucleus during interphase. We propose a model that satellite DNAs are critical chromosome elements that are recognized by satellite DNA-binding proteins and incorporated into chromocenters. We suggest that chromocenters function to preserve the entire chromosomal complement in a single nucleus, a fundamental and unquestioned feature of eukaryotic genomes. We speculate that the rapid divergence of satellite DNA sequences between closely related species results in discordant chromocenter function and may underlie speciation and hybrid incompatibility.

Improved de novo genomic assembly for the domestic donkey

Donkeys and horses share a common ancestor dating back to about 4 million years ago. Although a high-quality genome assembly at the chromosomal level is available for the horse, current assemblies available for the donkey are limited to moderately sized scaffolds. The absence of a better-quality assembly for the donkey has hampered studies involving the characterization of patterns of genetic variation at the genome-wide scale. These range from the application of genomic tools to selective breeding and conservation to the more fundamental characterization of the genomic loci underlying speciation and domestication. We present a new high-quality donkey genome assembly obtained using the Chicago HiRise assembly technology, providing scaffolds of subchromosomal size. We make use of this new assembly to obtain more accurate measures of heterozygosity for equine species other than the horse, both genome-wide and locally, and to detect runs of homozygosity potentially pertaining to positive selection in domestic donkeys. Finally, this new assembly allowed us to identify fine-scale chromosomal rearrangements between the horse and the donkey that likely played an active role in their divergence and, ultimately, speciation.

Contrasted patterns of evolution of the LINE-1 retrotransposon in perissodactyls: the history of a LINE-1 extinction

Background

LINE-1 (L1) is the dominant autonomously replicating non-LTR retrotransposon in mammals. Although our knowledge of L1 evolution across the tree of life has considerably improved in recent years, what we know of L1 evolution in mammals is biased and comes mostly from studies in primates (mostly human) and rodents (mostly mouse). It is unclear if patterns of evolution that are shared between those two groups apply to other mammalian orders. Here we performed a detailed study on the evolution of L1 in perissodactyls by making use of the complete genome of the domestic horse and of the white rhinoceros. This mammalian order offers an excellent model to study the extinction of L1 since the rhinoceros is one of the few mammalian species to have lost active L1.

Results

We found that multiple L1 lineages, carrying different 5’UTRs, have been simultaneously active during the evolution of perissodactyls. We also found that L1 has continuously amplified and diversified in horse. In rhinoceros, L1 was very prolific early on. Two successful families were simultaneously active until ~20my ago but became extinct suddenly at exactly the same time.

Conclusions

The general pattern of L1 evolution in perissodactyls is very similar to what was previously described in mouse and human, suggesting some commonalities in the way mammalian genomes interact with L1. We confirmed the extinction of L1 in rhinoceros and we discuss several possible mechanisms.

Electronic supplementary material

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

Identification and expression analysis of cDNA encoding insulin-like growth factor 2 in horses

Insulin-like growth factor 2 (IGF2) is responsible for a broad range of physiological processes during fetal development and adulthood, but genomic analyses of IGF2 containing the 5ʹ- and 3ʹ-untranslated regions (UTRs) in equines have been limited. In this study, we characterized the IGF2 mRNA containing the UTRs, and determined its expression pattern in the fetal tissues of horses. The complete equine IGF2 mRNA sequence harboring another exon approximately 2.8 kb upstream from the canonical transcription start site was identified as a new transcript variant. As this upstream exon did not contain the start codon, the amino acid sequence was identical to the canonical variant. Analysis of the deduced amino acid sequence revealed that the protein possessed two major domains, IlGF and IGF2_C, and analysis of IGF2 sequence polymorphism in fetal tissues of Hokkaido native horse and Thoroughbreds revealed a single nucleotide polymorphism (T to C transition) at position 398 in Thoroughbreds, which caused an amino acid substitution at position 133 in the IGF2 sequence. Furthermore, the expression pattern of the IGF2 mRNA in the fetal tissues of horses was determined for the first time, and was found to be consistent with those of other species. Taken together, these results suggested that the transcriptional and translational products of the IGF2 gene have conserved functions in the fetal development of mammals, including horses.

Ancient genomes revisit the ancestry of domestic and Przewalski’s 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.

Construction of two whole genome radiation hybrid panels for dromedary (Camelus dromedarius): 5000RAD and 15000RAD

The availability of genomic resources including linkage information for camelids has been very limited. Here, we describe the construction of a set of two radiation hybrid (RH) panels (5000RAD and 15000RAD) for the dromedary (Camelus dromedarius) as a permanent genetic resource for camel genome researchers worldwide. For the 5000RAD panel, a total of 245 female camel-hamster radiation hybrid clones were collected, of which 186 were screened with 44 custom designed marker loci distributed throughout camel genome. The overall mean retention frequency (RF) of the final set of 93 hybrids was 47.7%. For the 15000RAD panel, 238 male dromedary-hamster radiation hybrid clones were collected, of which 93 were tested using 44 PCR markers. The final set of 90 clones had a mean RF of 39.9%. This 15000RAD panel is an important high-resolution complement to the main 5000RAD panel and an indispensable tool for resolving complex genomic regions. This valuable genetic resource of dromedary RH panels is expected to be instrumental for constructing a high resolution camel genome map. Construction of the set of RH panels is essential step toward chromosome level reference quality genome assembly that is critical for advancing camelid genomics and the development of custom genomic tools.

The Hidden Genomic and Transcriptomic Plasticity of Giant Marker Chromosomes in Cancer

Genome amplification in the form of rings or giant rod-shaped marker chromosomes (RGMs) is a common genetic alteration in soft tissue tumors. The mitotic stability of these structures is often rescued by perfectly functioning analphoid neocentromeres, which therefore significantly contribute to cancer progression. Here, we disentangled the genomic architecture of many neocentromeres stabilizing marker chromosomes in well-differentiated liposarcoma and lung sarcomatoid carcinoma samples. In cells carrying heavily rearranged RGMs, these structures were assembled as patchworks of multiple short amplified sequences, disclosing an extremely high level of complexity and definitely ruling out the existence of regions prone to neocentromere seeding. Moreover, by studying two well-differentiated liposarcoma samples derived from the onset and the recurrence of the same tumor, we documented an expansion of the neocentromeric domain that occurred during tumor progression, which reflects a strong selective pressure acting toward the improvement of the neocentromeric functionality in cancer. In lung sarcomatoid carcinoma cells we documented, extensive "centromere sliding" phenomena giving rise to multiple, closely mapping neocentromeric epialleles on separate coexisting markers occur, likely due to the instability of neocentromeres arising in cancer cells. Finally, by investigating the transcriptional activity of neocentromeres, we came across a burst of chimeric transcripts, both by extremely complex genomic rearrangements, and cis/trans-splicing events. Post-transcriptional editing events have been reported to expand and variegate the genetic repertoire of higher eukaryotes, so they might have a determining role in cancer. The increased incidence of fusion transcripts, might act as a driving force for the genomic amplification process, together with the increased transcription of oncogenes.

A genome-wide association study for body weight in Japanese Thoroughbred racehorses clarifies candidate regions on chromosomes 3, 9, 15, and 18

Body weight is an important trait to confirm growth and development in humans and animals. In Thoroughbred racehorses, it is measured in the postnatal, training, and racing periods to evaluate growth and training degrees. The body weight of mature Thoroughbred racehorses generally ranges from 400 to 600 kg, and this broad range is likely influenced by environmental and genetic factors. Therefore, a genome-wide association study (GWAS) using the Equine SNP70 BeadChip was performed to identify the genomic regions associated with body weight in Japanese Thoroughbred racehorses using 851 individuals. The average body weight of these horses was 473.9 kg (standard deviation: 28.0) at the age of 3, and GWAS identified statistically significant SNPs on chromosomes 3 (BIEC2_808466, P=2.32E-14), 9 (BIEC2_1105503, P=1.03E-7), 15 (BIEC2_322669, P=9.50E-6), and 18 (BIEC2_417274, P=1.44E-14), which were associated with body weight as a quantitative trait. The genomic regions on chromosomes 3, 9, 15, and 18 included ligand-dependent nuclear receptor compressor-like protein (LCORL), zinc finger and AT hook domain containing (ZFAT), tribbles pseudokinase 2 (TRIB2), and myostatin (MSTN), respectively, as candidate genes. LCORL and ZFAT are associated with withers height in horses, whereas MSTN affects muscle mass. Thus, the genomic regions identified in this study seem to affect the body weight of Thoroughbred racehorses. Although this information is useful for breeding and growth management of the horses, the production of genetically modified animals and gene doping (abuse/misuse of gene therapy) should be prohibited to maintain horse racing integrity.

Genomic Analysis of Demographic History and Ecological Niche Modeling in the Endangered Sumatran Rhinoceros Dicerorhinus sumatrensis

The vertebrate extinction rate over the past century is approximately 22-100 times greater than background extinction rates [1], and large mammals are particularly at risk [2, 3]. Quaternary megafaunal extinctions have been attributed to climate change [4], overexploitation [5], or a combination of the two [6]. Rhinoceroses (Family: Rhinocerotidae) have a rich fossil history replete with iconic examples of climate-induced extinctions [7], but current pressures threaten to eliminate this group entirely. The Sumatran rhinoceros (Dicerorhinus sumatrensis) is among the most imperiled mammals on earth. The 2011 population was estimated at ≤216 wild individuals [8], and currently the species is extirpated, or nearly so, throughout the majority of its former range [8-12]. Understanding demographic history is important in placing current population status into a broader ecological and evolutionary context. Analysis of the Sumatran rhinoceros genome reveals extreme changes in effective population size throughout the Pleistocene. Population expansion during the early to middle Pleistocene was followed by decline. Ecological niche modeling indicated that changing climate most likely played a role in the decline of the Sumatran rhinoceros, as less suitable habitat on an emergent Sundaland corridor isolated Sumatran rhinoceros populations. By the end of the Pleistocene, the Sundaland corridor was submerged, and populations were fragmented and consequently reduced to low Holocene levels from which they would never recover. Past events denuded the Sumatran rhinoceros of genetic diversity through population decline, fragmentation, or some combination of the two and most likely made the species even more susceptible to later exploitation and habitat loss. VIDEO ABSTRACT.

Differential Expression of Serum MicroRNAs Supports CD4⁺ T Cell Differentiation into Th2/Th17 Cells in Severe Equine Asthma

MicroRNAs (miRNAs) regulate post-transcriptional gene expression and may be exported from cells via exosomes or in partnership with RNA-binding proteins. MiRNAs in body fluids can act in a hormone-like manner and play important roles in disease initiation and progression. Hence, miRNAs are promising candidates as biomarkers. To identify serum miRNA biomarkers in the equine model of asthma we investigated small RNA derived from the serum of 34 control and 37 asthmatic horses. These samples were used for next generation sequencing, novel miRNA identification and differential miRNA expression analysis. We identified 11 significantly differentially expressed miRNAs between case and control horses: eca-miR-128, eca-miR-744, eca-miR-197, eca-miR-103, eca-miR-107a, eca-miR-30d, eca-miR-140-3p, eca-miR-7, eca-miR-361-3p, eca-miR-148b-3p and eca-miR-215. Pathway enrichment using experimentally validated target genes of the human homologous miRNAs showed a significant enrichment in the regulation of epithelial-to-mesenchymal transition (key player in airway remodeling in asthma) and the phosphatidylinositol (3,4,5)-triphosphate (PIP3) signaling pathway (modulator of CD4⁺ T cell maturation and function). Downregulated miR-128 and miR-744 supports a Th2/Th17 type immune response in severe equine asthma.

MHC haplotype diversity in Persian Arabian horses determined using polymorphic microsatellites

Previous research on the equine major histocompatibility complex (MHC) demonstrated strong correlations between haplotypes defined by polymorphic intra-MHC microsatellites and haplotypes defined using classical serology. Here, we estimated MHC diversity in a sample of 124 Arabian horses from an endangered strain native to Iran (Persian Asil Arabians), using a validated 10-marker microsatellite panel. In a group of 66 horses related as parent-offspring pairs or half-sibling groups, we defined 51 MHC haplotypes, 49 of which were new. In 47 of the remaining 58 unrelated horses, we could assign one previously identified MHC haplotype, and by default, we gave provisional haplotype status to the remaining constellation of microsatellite alleles. In these horses, we found 21 haplotypes that we had previously defined and 31 provisional haplotypes, two of which had been identified in an earlier study. This gave a total of 78 new MHC haplotypes. The final 11 horses were MHC heterozygotes that we could not phase using information from any of the previously validated or provisional haplotypes. However, we could determine that these horses carried a total of 22 different undefined haplotypes. In the overall population sample, we detected three homozygous horses and one maternally inherited recombinant from 21 informative segregations. Virtually all of the horses tested were MHC heterozygotes, and most unrelated horses (98%) were heterozygous for rare microsatellite-defined haplotypes found less than three times in the sampled horses. This is evidence for a very high level of MHC haplotype variation in the Persian Asil Arabian horse.

A missense variant in the coil1A domain of the keratin 25 gene is associated with the dominant curly hair coat trait (Crd) in horse

Background

Curly horses present a variety of curl phenotypes that are associated with various degrees of curliness of coat, mane, tail and ear hairs. Their origin is still a matter of debate and several genetic hypotheses have been formulated to explain the diversity in phenotype, including the combination of autosomal dominant and recessive alleles. Our purpose was to map the autosomal dominant curly hair locus and identify the causal variant using genome-wide association study (GWAS) and whole-genome sequencing approaches.

Results

A GWAS was performed using a Bayesian sparse linear mixed model, based on 51 curly and 19 straight-haired French and North American horses from 13 paternal families genotyped on the Illumina EquineSNP50 BeadChip. A single strong signal was observed on equine chromosome 11, in a region that encompasses the type I keratin gene cluster. This region was refined by haplotype analysis to a segment including 36 genes, among which are 10 keratin genes (KRT-10, -12, -20, -23, -24, -25, -26, -27, -28, -222). To comprehensively identify candidate causal variants within all these genes, whole-genome sequences were obtained for one heterozygous curly stallion and its straight-haired son. Among the four non-synonymous candidate variants identified and validated in the curly region, only variant g.21891160G>A in the KRT25 gene (KRT25:p.R89H) was in perfect agreement with haplotype status in the whole pedigree. Genetic association was then confirmed by genotyping a larger population consisting of 353 horses. However, five discordant curly horses were observed, which carried neither the variant nor the main haplotype associated with curliness. Sequencing of KRT25 for two discordant horses did not identify any other deleterious variant, which suggests locus rather than allelic heterogeneity for the curly phenotype.

Conclusions

We identified the KRT25:p.R89H variant as responsible for the dominant curly trait, but a second dominant locus may also be involved in the shape of hairs within North American Curly horses.

Electronic supplementary material

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

Clinical veterinary proteomics: Techniques and approaches to decipher the animal plasma proteome

Over the last two decades, technological advancements in the field of proteomics have advanced our understanding of the complex biological systems of living organisms. Techniques based on mass spectrometry (MS) have emerged as powerful tools to contextualise existing genomic information and to create quantitative protein profiles from plasma, tissues or cell lines of various species. Proteomic approaches have been used increasingly in veterinary science to investigate biological processes responsible for growth, reproduction and pathological events. However, the adoption of proteomic approaches by veterinary investigators lags behind that of researchers in the human medical field. Furthermore, in contrast to human proteomics studies, interpretation of veterinary proteomic data is difficult due to the limited protein databases available for many animal species. This review article examines the current use of advanced proteomics techniques for evaluation of animal health and welfare and covers the current status of clinical veterinary proteomics research, including successful protein identification and data interpretation studies. It includes a description of an emerging tool, sequential window acquisition of all theoretical fragment ion mass spectra (SWATH-MS), available on selected mass spectrometry instruments. This newly developed data acquisition technique combines advantages of discovery and targeted proteomics approaches, and thus has the potential to advance the veterinary proteomics field by enhancing identification and reproducibility of proteomics data.

To pace or not to pace: a pilot study of four- and five-gaited Icelandic horses homozygous for the DMRT3 'Gait Keeper' mutation

The Icelandic horse is a breed known mainly for its ability to perform the ambling four-beat gait 'tölt' and the lateral two-beat gait pace. The natural ability of the breed to perform these alternative gaits is highly desired by breeders. Therefore, the discovery that a nonsense mutation (C>A) in the DMRT3 gene was the main genetic factor for horses' ability to perform gaits in addition to walk, trot and canter was of great interest. Although several studies have demonstrated that homozygosity for the DMRT3 mutation is important for the ability to pace, only about 70% of the homozygous mutant (AA) Icelandic horses are reported to pace. The aim of the study was to genetically compare four- and five-gaited (i.e. horses with and without the ability to pace) AA Icelandic horses by performing a genome-wide association (GWA) analysis. All horses (n = 55) were genotyped on the 670K Axiom Equine Genotyping Array, and a GWA analysis was performed using the genabel package in r. No SNP demonstrated genome-wide significance, implying that the ability to pace goes beyond the presence of a single gene variant. Despite its limitations, the current study provides additional information regarding the genetic complexity of pacing ability in horses. However, to fully understand the genetic differences between four- and five-gaited AA horses, additional studies with larger sample materials and consistent phenotyping are needed.

Satellite DNA: An Evolving Topic

Satellite DNA represents one of the most fascinating parts of the repetitive fraction of the eukaryotic genome. Since the discovery of highly repetitive tandem DNA in the 1960s, a lot of literature has extensively covered various topics related to the structure, organization, function, and evolution of such sequences. Today, with the advent of genomic tools, the study of satellite DNA has regained a great interest. Thus, Next-Generation Sequencing (NGS), together with high-throughput in silico analysis of the information contained in NGS reads, has revolutionized the analysis of the repetitive fraction of the eukaryotic genomes. The whole of the historical and current approaches to the topic gives us a broad view of the function and evolution of satellite DNA and its role in chromosomal evolution. Currently, we have extensive information on the molecular, chromosomal, biological, and population factors that affect the evolutionary fate of satellite DNA, knowledge that gives rise to a series of hypotheses that get on well with each other about the origin, spreading, and evolution of satellite DNA. In this paper, I review these hypotheses from a methodological, conceptual, and historical perspective and frame them in the context of chromosomal organization and evolution.

Characterization of gene expression and genetic variation of horse ERBB receptor feedback inhibitor 1 in Thoroughbreds

Objective

This study aimed to test the expression patterns of ERBB receptor feedback inhibitor 1 (ERRFI1) before and after exercise and the association of non-synonymous single-nucleotide polymorphisms (nsSNPs) of horse ERRFI1 with racing traits in Thoroughbreds.

Methods

We performed bioinformatics and gene expression analyses for horse ERRFI1. Transcription factor (TF) binding sites in the 5′-regulatory region of this gene were identified through a tool for prediction of TF-binding site (PROMO). A general linear model was used to detect the association between the nsSNP (LOC42830758 A to G) and race performance.

Results

Quantitative polymerase chain reaction analysis showed that expression level of ERRFI1 after exercise was 1.6 times higher than that before exercise. Ten transcription factors were predicted from the ERRFI1 regulatory region. A novel nsSNP (LOC42830758 A to G) was found in ERRFI1, which was associated with three racing traits including average prize money, average racing index, and 3-year-old starts percentile ranking.

Conclusion

Our analysis will be helpful as a basis for studying genes and SNPs that affect race performance in racehorses.

Molecular cloning and functional expression of the K+ channel KV7.1 and the regulatory subunit KCNE1 from equine myocardium

BACKGROUND

The voltage-gated K⁺-channel K_V7.1 and the subunit KCNE1, encoded by the KCNQ1 and KCNE1 genes, respectively, are responsible for termination of the cardiac action potential. In humans, mutations in these genes can predispose patients to arrhythmias and sudden cardiac death (SCD).

AIM

To characterize equine K_V7.1/KCNE1 currents and compare them to human K_V7.1/KCNE1 currents to determine whether K_V7.1/KCNE1 plays a similar role in equine and human hearts.

METHODS

mRNA encoding K_V7.1 and KCNE1 was isolated from equine hearts, sequenced, and cloned into expression vectors. The channel subunits were heterologously expressed in Xenopus laevis oocytes or CHO-K1 cells and characterized using voltage-clamp techniques.

RESULTS

Equine K_V7.1/KCNE1 expressed in CHO-K1 cells exhibited electrophysiological properties that are overall similar to the human orthologs; however, a slower deactivation was found which could result in more open channels at fast rates.

CONCLUSION

The results suggest that the equine K_V7.1/KCNE1 channel may be important for cardiac repolarization and this could indicate that horses are susceptible to SCD caused by mutations in KCNQ1 and KCNE1.

Single nucleotide polymorphisms for DNA typing in the domestic horse

Genetic markers are important resources for individual identification and parentage assessment. Although short tandem repeats (STRs) have been the traditional DNA marker, technological advances have led to single nucleotide polymorphisms (SNPs) becoming an attractive alternative. SNPs can be highly multiplexed and automatically scored, which allows for easier standardization and sharing among laboratories. Equine parentage is currently assessed using STRs. We obtained a publicly available SNP dataset of 729 horses representing 32 diverse breeds. A proposed set of 101 SNPs was analyzed for DNA typing suitability. The overall minor allele frequency of the panel was 0.376 (range 0.304-0.419), with per breed probability of identities ranging from 5.6 × 10^-35 to 1.86 × 10^-42 . When one parent was available, exclusion probabilities ranged from 0.9998 to 0.999996, although when both parents were available, all breeds had exclusion probabilities greater than 0.9999999. A set of 388 horses from 35 breeds was genotyped to evaluate marker performance on known families. The set included 107 parent-offspring pairs and 101 full trios. No horses shared identical genotypes across all markers, indicating that the selected set was sufficient for individual identification. All pairwise comparisons were classified using ISAG rules, with one or two excluding markers considered an accepted parent-offspring pair, two or three excluding markers considered doubtful and four or more excluding markers rejecting parentage. The panel had an overall accuracy of 99.9% for identifying true parent-offspring pairs. Our developed marker set is both present on current generation SNP chips and can be highly multiplexed in standalone panels and thus is a promising resource for SNP-based DNA typing.

Impaired response of the bronchial epithelium to inflammation characterizes severe equine asthma

BACKGROUND

Severe equine asthma is a naturally occurring lung inflammatory disease of mature animals characterized by neutrophilic inflammation, bronchoconstriction, mucus hypersecretion and airway remodeling. Exacerbations are triggered by inhalation of dust and microbial components. Affected animals eventually are unable of aerobic performance. In this study transcriptomic differences between asthmatic and non-asthmatic animals in the response of the bronchial epithelium to an inhaled challenge were determined.

RESULTS

Paired endobronchial biopsies were obtained pre- and post-challenge from asthmatic and non-asthmatic animals. The transcriptome, determined by RNA-seq and analyzed with edgeR, contained 111 genes differentially expressed (DE) after challenge between horses with and without asthma, and 81 of these were upregulated. Genes involved in neutrophil migration and activation were in central location in interaction networks, and related gene ontology terms were significantly overrepresented. Relative abundance of specific gene products as determined by immunohistochemistry was correlated with differential gene expression. Gene sets involved in neutrophil chemotaxis, immune and inflammatory response, secretion, blood coagulation and apoptosis were overrepresented among up-regulated genes, while the rhythmic process gene set was overrepresented among down-regulated genes. MMP1, IL8, TLR4 and MMP9 appeared to be the most important proteins in connecting the STRING protein network of DE genes.

CONCLUSIONS

Several differentially expressed genes and networks in horses with asthma also contribute to human asthma, highlighting similarities between severe human adult and equine asthma. Neutrophil activation by the bronchial epithelium is suggested as the trigger of the inflammatory cascade in equine asthma, followed by epithelial injury and impaired repair and differentiation. Circadian rhythm dysregulation and the sonic Hedgehog pathway were identified as potential novel contributory factors in equine asthma.

Foundations of performance – factors that contribute to excellence in equine exercise

Horses are renowned for their incredible capacity for a range of athletic activities, and participation in athletic events arguably represents the most critical strut of the equine industry. Successful performance is typically a primary focus during participation in competitive athletic events, and relies upon a variety of innate physiological and structural factors of the athlete. However, a wide range of external factors also influence performance, and many of these can be readily manipulated. Therefore, thorough assessment of the individual’s inherent capacity for a specific athletic discipline must be combined with optimisation of external factors including nutrition and training to promote excellent performance. Recent progress in methods of athlete selection and monitoring of training responses are assisting continued improvements in equine performance.

A curated catalog of canine and equine keratin genes

Keratins represent a large protein family with essential structural and functional roles in epithelial cells of skin, hair follicles, and other organs. During evolution the genes encoding keratins have undergone multiple rounds of duplication and humans have two clusters with a total of 55 functional keratin genes in their genomes. Due to the high similarity between different keratin paralogs and species-specific differences in gene content, the currently available keratin gene annotation in species with draft genome assemblies such as dog and horse is still imperfect. We compared the National Center for Biotechnology Information (NCBI) (dog annotation release 103, horse annotation release 101) and Ensembl (release 87) gene predictions for the canine and equine keratin gene clusters to RNA-seq data that were generated from adult skin of five dogs and two horses and from adult hair follicle tissue of one dog. Taking into consideration the knowledge on the conserved exon/intron structure of keratin genes, we annotated 61 putatively functional keratin genes in both the dog and horse, respectively. Subsequently, curators in the RefSeq group at NCBI reviewed their annotation of keratin genes in the dog and horse genomes (Annotation Release 104 and Annotation Release 102, respectively) and updated annotation and gene nomenclature of several keratin genes. The updates are now available in the NCBI Gene database (https://www.ncbi.nlm.nih.gov/gene).

RNA Interference and BMP-2 Stimulation Allows Equine Chondrocytes Redifferentiation in 3D-Hypoxia Cell Culture Model: Application for Matrix-Induced Autologous Chondrocyte Implantation

As in humans, osteoarthritis (OA) causes considerable economic loss to the equine industry. New hopes for cartilage repair have emerged with the matrix-associated autologous chondrocyte implantation (MACI). Nevertheless, its limitation is due to the dedifferentiation occurring during the chondrocyte amplification phase, leading to the loss of its capacity to produce a hyaline extracellular matrix (ECM). To enhance the MACI therapy efficiency, we have developed a strategy for chondrocyte redifferentiation, and demonstrated its feasibility in the equine model. Thus, to mimic the cartilage microenvironment, the equine dedifferentiated chondrocytes were cultured in type I/III collagen sponges for 7 days under hypoxia in the presence of BMP-2. In addition, chondrocytes were transfected by siRNA targeting Col1a1 and Htra1 mRNAs, which are overexpressed during dedifferentiation and OA. To investigate the quality of the neo-synthesized ECM, specific and atypical cartilage markers were evaluated by RT-qPCR and Western blot. Our results show that the combination of 3D hypoxia cell culture, BMP-2 (Bone morphogenetic protein-2), and RNA interference, increases the chondrocytes functional indexes (Col2a1/Col1a1, Acan/Col1a1), leading to an effective chondrocyte redifferentiation. These data represent a proof of concept for this process of application, in vitro, in the equine model, and will lead to the improvement of the MACI efficiency for cartilage tissue engineering therapy in preclinical/clinical trials, both in equine and human medicine.

Developing a 670k genotyping array to tag ~2M SNPs across 24 horse breeds

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.

Electronic supplementary material

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

Genomewide association study reveals a risk locus for equine metabolic syndrome in the Arabian horse

Equine obesity can cause life-threatening secondary chronic conditions, similar to those in humans and other animal species. Equine metabolic syndrome (EMS), primarily characterized by hyperinsulinemia, is often present in obese horses and ponies. Due to clinical similarities to conditions such as pituitary pars intermedia dysfunction (formerly equine Cushing's disease), conclusive diagnosis of EMS often proves challenging. Aside from changes in diet and exercise, few targeted treatments are available for EMS, emphasizing the need for early identification of at-risk individuals to enable implementation of preventative measures. A genomewide association study (GWAS) using Arabian horses with a history of severe laminitis secondary to EMS revealed significant genetic markers near a single candidate gene () that may play a role in cholesterol homeostasis. The best marker, BIEC2-263524 (chr14:69276814 T > C), was correlated with elevated insulin values and increased frequency of laminitis ( = 0.0024 and = 9.663 × 10, respectively). In a second population of Arabian horses, the BIEC2-263524 marker maintained its associations with higher modified insulin-to-glucose ratio (MIRG) values ( = 0.0056) and BCS ( = 0.0063). Screening of the predicted coding regions by sequencing identified a polymorphic guanine homopolymer and 5 haplotypes in the 3' untranslated region (UTR). An 11 guanine (11-G) allele at was correlated with elevated insulin values in the GWAS population ( = 0.0008) and, in the second population, elevated MIRG and increased BCS > 6.5 ( = 0.0055 and = 0.0162, respectively). The BIEC2-263524-C and the 3' UTR -11(G) polymorphisms were correlated at a 98% frequency, indicating strong linkage disequilibrium across this 150-kb haplotype. Assays for these markers could diagnose horses with a genetic predisposition to develop obesity. Additionally, discovery of FAM174A function may improve our understanding of the etiology of this troubling illness in the horse and warrants investigation of this locus for a role in metabolic- and obesity-related disorders of other species.

Genomic landscape of copy number variation and copy neutral loss of heterozygosity events in equine sarcoids reveals increased instability of the sarcoid genome

Although they are the most common neoplasms in equids, sarcoids are not fully characterized at the molecular level. Therefore, the objective of this study was to characterize the landscape of structural rearrangements, such as copy number variation (CNV) and copy neutral loss of heterozygosity (cnLOH), in the genomes of sarcoid tumor cells. This information will not only broaden our understanding of the characteristics of this genome but will also improve the general knowledge of this tumor and the mechanisms involved in its generation. To this end, Equine SNP64K Illumina microarrays were applied along with bioinformatics tools dedicated for signal intensity analysis. The analysis revealed increased instability of the genome of sarcoid cells compared with unaltered skin tissue samples, which was manifested by the prevalence of CNV and cnLOH events. Many of the identified CNVs overlapped with the other research results, but the simultaneously observed variability in the number and sizes of detected aberrations indicated a need for further studies and the development of more reliable bioinformatics algorithms. The functional analysis of genes co-localized with the identified aberrations revealed that these genes are engaged in vital cellular processes. In addition, a number of these genes directly contribute to neoplastic transformation. Furthermore, large numbers of cnLOH events identified in the sarcoids suggested that they may play no less significant roles than CNVs in the carcinogenesis of this tumor. Thus, our results indicate the importance of cnLOH and CNV in equine sarcoid oncogenesis and present a direction of future research.

The wolf reference genome sequence (Canis lupus lupus) and its implications for Canis spp. population genomics

Background

An increasing number of studies are addressing the evolutionary genomics of dog domestication, principally through resequencing dog, wolf and related canid genomes. There is, however, only one de novo assembled canid genome currently available against which to map such data - that of a boxer dog (Canis lupus familiaris). We generated the first de novo wolf genome (Canis lupus lupus) as an additional choice of reference, and explored what implications may arise when previously published dog and wolf resequencing data are remapped to this reference.

Results

Reassuringly, we find that regardless of the reference genome choice, most evolutionary genomic analyses yield qualitatively similar results, including those exploring the structure between the wolves and dogs using admixture and principal component analysis. However, we do observe differences in the genomic coverage of re-mapped samples, the number of variants discovered, and heterozygosity estimates of the samples.

Conclusion

In conclusion, the choice of reference is dictated by the aims of the study being undertaken; if the study focuses on the differences between the different dog breeds or the fine structure among dogs, then using the boxer reference genome is appropriate, but if the aim of the study is to look at the variation within wolves and their relationships to dogs, then there are clear benefits to using the de novo assembled wolf reference genome.

Electronic supplementary material

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

Y Chromosome Uncovers the Recent Oriental Origin of Modern Stallions

The Y chromosome directly reflects male genealogies, but the extremely low Y chromosome sequence diversity in horses has prevented the reconstruction of stallion genealogies [1, 2]. Here, we resolve the first Y chromosome genealogy of modern horses by screening 1.46 Mb of the male-specific region of the Y chromosome (MSY) in 52 horses from 21 breeds. Based on highly accurate pedigree data, we estimated the de novo mutation rate of the horse MSY and showed that various modern horse Y chromosome lineages split much later than the domestication of the species. Apart from few private northern European haplotypes, all modern horse breeds clustered together in a roughly 700-year-old haplogroup that was transmitted to Europe by the import of Oriental stallions. The Oriental horse group consisted of two major subclades: the Original Arabian lineage and the Turkoman horse lineage. We show that the English Thoroughbred MSY was derived from the Turkoman lineage and that English Thoroughbred sires are largely responsible for the predominance of this haplotype in modern horses.

Endurance Exercise Ability in the Horse: A Trait with Complex Polygenic Determinism

Endurance horses are able to run at more than 20 km/h for 160 km (in bouts of 30-40 km). This level of performance is based on intense aerobic metabolism, effective body heat dissipation and the ability to endure painful exercise. The known heritabilities of endurance performance and exercise-related physiological traits in Arabian horses suggest that adaptation to extreme endurance exercise is influenced by genetic factors. The objective of the present genome-wide association study (GWAS) was to identify single nucleotide polymorphisms (SNPs) related to endurance racing performance in 597 Arabian horses. The performance traits studied were the total race distance, average race speed and finishing status (qualified, eliminated or retired). We used three mixed models that included a fixed allele or genotype effect and a random, polygenic effect. Quantile-quantile plots were acceptable, and the regression coefficients for actual vs. expected log10p-values ranged from 0.865 to 1.055. The GWAS revealed five significant quantitative trait loci (QTL) corresponding to 6 SNPs on chromosomes 6, 1, 7, 16, and 29 (two SNPs) with corrected p-values from 1.7 × 10-6 to 1.8 × 10-5. Annotation of these 5 QTL revealed two genes: sortilin-related VPS10-domain-containing receptor 3 (SORCS3) on chromosome 1 is involved in protein trafficking, and solute carrier family 39 member 12 (SLC39A12) on chromosome 29 is active in zinc transport and cell homeostasis. These two coding genes could be involved in neuronal tissues (CNS). The other QTL on chromosomes 6, 7, and 16 may be involved in the regulation of the gene expression through non-coding RNAs, CpG islands and transcription factor binding sites. On chromosome 6, a new candidate equine long non-coding RNA (KCNQ1OT1 ortholog: opposite antisense transcript 1 of potassium voltage-gated channel subfamily Q member 1 gene) was predicted in silico and validated by RT-qPCR in primary cultures of equine myoblasts and fibroblasts. This lncRNA could be one element of the cardiac rhythm regulation. Our GWAS revealed that equine performance during endurance races is a complex polygenic trait, and is partially governed by at least 5 QTL: two coding genes involved in neuronal tissues and three other loci with many regulatory functions such as slowing down heart rate.