Genome data uncover four synergistic key regulators for extremely small body size in horses

Genomic analyses of withers height and linear conformation traits in German Warmblood horses using imputed sequence-level genotypes

BACKGROUND

Body conformation, including withers height, is a major selection criterion in horse breeding and is associated with other important traits, such as health and performance. However, little is known about the genomic background of equine conformation. Therefore, the aim of this study was to use imputed sequence-level genotypes from up to 4891 German Warmblood horses to identify genomic regions associated with withers height and linear conformation traits. Furthermore, the traits were genetically characterised and putative causal variants for withers height were detected.

RESULTS

A genome-wide association study (GWAS) for withers height confirmed the presence of a previously known quantitative trait locus (QTL) on Equus caballus (ECA) chromosome 3 close to the LCORL/NCAPG locus, which explained 16% of the phenotypic variance for withers height. An additional significant association signal was detected on ECA1. Further investigations of the region on ECA3 identified a few promising candidate causal variants for withers height, including a nonsense mutation in the coding sequence of the LCORL gene. The estimated heritability for withers height was 0.53 and ranged from 0 to 0.34 for the conformation traits. GWAS identified significantly associated variants for more than half of the investigated conformation traits, among which 13 showed a peak on ECA3 in the same region as withers height. Genetic parameter estimation revealed high genetic correlations between these traits and withers height for the QTL on ECA3.

CONCLUSIONS

The use of imputed sequence-level genotypes from a large study cohort led to the discovery of novel QTL associated with conformation traits in German Warmblood horses. The results indicate the high relevance of the QTL on ECA3 for various conformation traits, including withers height, and contribute to deciphering causal mutations for body size in horses.

Selection Signatures Reveal Candidate Genes for the Cornish Rex Breed-Specific Phenotype

Many coat color, behavioral and morphological traits are specific and fixed across cat breeds, with several variants influencing these traits being common among different breeds. In the domestic cat, rexoid mutations have been documented in several breeds. In the Cornish Rex, four bp deletion in the LPAR6 gene has been found to cause a frame shift and a premature stop codon. In addition to the rexoid coat, Cornish Rex cats also have a characteristic head, ear shape and body type. Analysis of the selection signatures in the Cornish Rex genome revealed several regions that are under selective pressure. One of these is located in CFA B4, in the region where the ALX1 gene is located. The ALX1 gene in Burmese cats disrupts the cranial morphogenesis and causes brachycephaly in the heterozygous state. In our study, we confirmed the presence of a deletion in LPAR6 in 20 Cornish Rex and in four F1 hybrids between Cornish Rex and domestic cat. However, we did not confirm the presence of the deletion in ALX1 in Cornish Rex cats. Genome-wide selection signature analysis was performed using ROH islands and integrated haplotype score (iHS) statistics based on publicly available SNP array data of 11 Cornish Rex cats. The selection signatures were detected on chromosomes A1, A3, C2, B1, B4 and D1.

Genome-Wide Detection of Copy Number Variations Associated with Miniature Features in Horses

Copy number variations (CNVs) are crucial structural genomic variants affecting complex traits in humans and livestock animals. The current study was designed to conduct a comprehensive comparative copy number variation analysis among three breeds, Debao (DB), Baise (BS), and Warmblood (WB), with a specific focus on identifying genomic regions associated with miniature features in horses. Using whole-genome next-generation resequencing data, we identified 18,974 CNVs across 31 autosomes. Among the breeds, we found 4279 breed-specific CNV regions (CNVRs). Baise, Debao, and Warmblood displayed 2978, 986, and 895 distinct CNVRs, respectively, with 202 CNVRs shared across all three breeds. After removing duplicates, we obtained 1545 CNVRs from 26 horse genomes. Functional annotation reveals enrichment in biological functions, including antigen processing, cell metabolism, olfactory conduction, and nervous system development. Debao horses have 970 genes overlapping with CNVRs, possibly causing their small size and mountainous adaptations. We also found that the genes GHR, SOX9, and SOX11 may be responsible for the miniature features of the Debao horse by analyzing their overlapping CNVRs. Overall, this study offers valuable insights into the widespread presence of CNVs in the horse genome. The findings contribute to mapping horse CNVs and advance research on unique miniature traits observed in the Debao horse.

Single-cut osteotomy for correction of a complex multiplanar deformity of the radius in a Shetland pony foal

OBJECTIVE

To describe the surgical correction of a multiplanar deformity of the radius in a pony using a single-cut osteotomy.

STUDY DESIGN

Case report.

ANIMALS

A 9-week-old male Shetland pony foal with a bodyweight of 47 kg.

METHODS

The foal presented with a complex multiplanar deformity of the right radius. A 3-dimensional model of the bone was created based on computed tomography (CT) imaging. To correct the deformity, the cutting plane for a single-cut osteotomy was calculated following the mathematical approach described by Sangeorzan et al. After osteotomy, the bone was realigned and stabilized with two 4.5 locking compression plates (LCPs).

RESULTS

Recovery from surgery was uneventful, and the foal remained comfortable. A CT exam 15 weeks after surgery revealed that diaphyseal deformities improved substantially in procurvatum (from 8° to 1°), varus (from 27° to 0°), and rotation (30° to 5°). The operated radius was 2.1 cm shorter than the left. Eighteen-month follow up confirmed a functionally and cosmetically acceptable outcome.

CONCLUSION

The single-cut osteotomy resulted in the successful correction of a multiplanar equine long-bone deformity with a favorable outcome in a Shetland pony.

CLINICAL SIGNIFICANCE

Single-cut osteotomy is an alternative surgical technique for the correction of complex diaphyseal long-bone equine deformities. Computed tomography data and the possibility of printing 3D models provides a significant advantage for rehearsing the procedure and for evaluating the correction that was achieved.

DIAPH3 deficiency links microtubules to mitotic errors, defective neurogenesis, and brain dysfunction

Diaphanous (DIAPH) three (DIAPH3) is a member of the formin proteins that have the capacity to nucleate and elongate actin filaments and, therefore, to remodel the cytoskeleton. DIAPH3 is essential for cytokinesis as its dysfunction impairs the contractile ring and produces multinucleated cells. Here, we report that DIAPH3 localizes at the centrosome during mitosis and regulates the assembly and bipolarity of the mitotic spindle. DIAPH3-deficient cells display disorganized cytoskeleton and multipolar spindles. DIAPH3 deficiency disrupts the expression and/or stability of several proteins including the kinetochore-associated protein SPAG5. DIAPH3 and SPAG5 have similar expression patterns in the developing brain and overlapping subcellular localization during mitosis. Knockdown of SPAG5 phenocopies DIAPH3 deficiency, whereas its overexpression rescues the DIAHP3 knockdown phenotype. Conditional inactivation of Diaph3 in mouse cerebral cortex profoundly disrupts neurogenesis, depleting cortical progenitors and neurons, leading to cortical malformation and autistic-like behavior. Our data uncover the uncharacterized functions of DIAPH3 and provide evidence that this protein belongs to a molecular toolbox that links microtubule dynamics during mitosis to aneuploidy, cell death, fate determination defects, and cortical malformation.

A QTL for conformation of back and croup influences lateral gait quality in Icelandic horses

BACKGROUND

The back plays a vital role in horse locomotion, where the spine functions as a spring during the stride cycle. A complex interaction between the spine and the muscles of the back contribute to locomotion soundness, gait ability, and performance of riding and racehorses. Conformation is commonly used to select horses for breeding and performance in multiple horse breeds, where the back and croup conformation plays a significant role. The conformation of back and croup plays an important role on riding ability in Icelandic horses. However, the genes behind this trait are still unknown. Therefore, the aim of this study was to identify genomic regions associated with conformation of back and croup in Icelandic horses and to investigate their effects on riding ability. One hundred seventy-seven assessed Icelandic horses were included in the study. A genome-wide association analysis was performed using the 670 K+ Axiom Equine Genotyping Array, and the effects of different haplotypes in the top associated region were estimated for riding ability and additional conformation traits assessed during breeding field tests.

RESULTS

A suggestive quantitative trait loci (QTL) for the score of back and croup was detected on Equus caballus (ECA) 22 (p-value = 2.67 × 10- 7). Haplotype analysis revealed two opposite haplotypes, which resulted in higher and lower scores of the back and croup, respectively (p-value < 0.001). Horses with the favorable haplotype were more inclined to have a well-balanced backline with an uphill conformation and had, on average, higher scores for the lateral gaits tölt (p-value = 0.02) and pace (p-value = 0.004). This genomic region harbors three genes: C20orf85, ANKRD60 and LOC100056167. ANKRD60 is associated with body height in humans. C20orf85 and ANKRD60 are potentially linked to adolescent idiopathic scoliosis in humans.

CONCLUSIONS

Our results show that the detected QTL for conformation of back and croup is of importance for quality of lateral gaits in Icelandic horses. These findings could result in a genetic test to aid in the selection of breeding horses, thus they are of major interest for horse breeders. The results may also offer a gateway to comparative functional genomics by potentially linking both motor laterality and back inclination in horses with scoliosis in humans.

Pathways involved in pony body size development

BACKGROUND

The mechanism of body growth in mammals is poorly understood. Here, we investigated the regulatory networks involved in body growth through transcriptomic analysis of pituitary and epiphyseal tissues of smaller sized Debao ponies and Mongolian horses at the juvenile and adult stages.

RESULTS

We found that growth hormone receptor (GHR) was expressed at low levels in long bones, although growth hormone (GH) was highly expressed in Debao ponies compared with Mongolian horses. Moreover, significant downregulated of the GHR pathway components m-RAS and ATF3 was found in juvenile ponies, which slowed the proliferation of bone osteocytes. However, WNT2 and PLCβ2 were obviously upregulated in juvenile Debao ponies, which led to premature mineralization of the bone extracellular matrix. Furthermore, we found that the WNT/Ca²⁺ pathway may be responsible for regulating body growth. GHR was demonstrated by q-PCR and Western blot analyses to be expressed at low levels in long bones of Debao ponies. Treatment with WNT antagonistI decreased the expression of WNT pathway components (P < 0.05) in vitro. Transduction of ATDC5 cells with a GHR-RNAi lentiviral vector decreased the expression of the GHR pathway components (P < 0.05). Additionally, the expression of the IGF-1 gene in the liver was lower in Debao ponies than in Mongolian horses at the juvenile and adult stages. Detection of plasma hormone concentrations showed that Debao ponies expressed higher levels of IGF-1 as juveniles and higher levels of GH as adults than Mongolian horses, indicating that the hormone regulation in Debao ponies differs from that in Mongolian horses.

CONCLUSION

Our work provides insights into the genetic regulation of short stature growth in mammals and can provide useful information for the development of therapeutic strategies for small size.

Evaluation of a new variant in the aggrecan gene potentially associated with chondrodysplastic dwarfism in Miniature horses

Chondrodysplastic dwarfism in Miniature horses is an autosomal recessive disorder previously associated with four mutations (D1, D2, D3*, and D4) in the aggrecan (ACAN) gene. The aim of this study was to identify additional variants in the candidate ACAN gene associated with chondrodysplastic dwarfism in Miniature horses. Fifteen dwarf Miniature horses were found to possess only one of the dwarfism-causing variants, and two possessed none of the variants. The ACAN exons (EquCab3.0) of seven dwarf Miniature horses were sequenced. A missense SNP in coding exon 11 (g.95271115A > T, c.6465A > T-RefSeq XM_005602799.2), which resulted in the amino acid substitution p.Leu2155Phe (RefSeq XP_005602856.2), was initially associated with the dwarf phenotype. The variant was tested and found present in 14 dwarf foals as well as one parent of each, and both parents of a dwarf possessing two copies. Genetic testing of 347 phenotypically normal Miniature horses demonstrated that none had more than one of the dwarf alleles or c.6465A > T. However, a study of large breeds revealed the presence of c.6465A > T, which was present in homozygosis in two Mangalarga Marchador horses. We suggest that c.6465A > T as a marker of disequilibrium or complex interactions in the Miniature horse genome could contribute to the associated dwarfism.

Signatures of selection analysis using whole-genome sequence data reveals novel candidate genes for pony and light horse types

Natural selection and domestication have shaped modern horse populations, resulting in a vast range of phenotypically diverse breeds. Horse breeds are classified into three types (pony, light, and draft) generally based on their body type. Understanding the genetic basis of horse type variation and selective pressures related to the evolutionary trend can be particularly important for current selection strategies. Whole-genome sequences were generated for 14 pony and 32 light horses to investigate the genetic signatures of selection of the horse type in pony and light horses. In the overlapping extremes of the fixation index and nucleotide diversity results, we found novel genomic signatures of selective sweeps near key genes previously implicated in body measurements including C4ORF33, CRB1, CPN1, FAM13A, and FGF12 that may influence variation in pony and light horse types. This study contributes to a better understanding of the genetic background of differences between pony and light horse types.

Heterozygotes for ACAN dwarfism alleles in horses have reduced stature

Homozygous and compound heterozygous Miniature horses for ACAN alleles D1, D2, D3* and D4 exhibit chondrodysplastic dwarfism (OMIA 001271-9796). In a previous study, the carrier rate for these four alleles, combined, was 26.2%. The purpose of this study was to investigate whether carriers of these dwarfism-causing alleles had a shorter withers height than non-carriers. A total of 245 Miniature horses were tested for these four ACAN alleles and also were measured for withers height. Of these horses, 98 were carriers and 147 were non-carriers. A statistically significant difference of 1.43 inches was observed with the carriers being shorter (P = 1.72E - 11). The range of heights for the two groups overlapped, indicating that other factors, including genes, have an impact on withers height. However, the high carrier rate of these dwarfism-causing variants may be due to selection for decreased height.

Identification of Novel lncRNAs Differentially Expressed in Placentas of Chinese Ningqiang Pony and Yili Horse Breeds

As a nutrient sensor, the placenta plays a key role in regulating fetus growth and development. Long non-coding RNAs (lncRNAs) have been shown to regulate growth-related traits. However, the biological function of lncRNAs in horse placentas remains unclear. To compare the expression patterns of lncRNAs in the placentas of the Chinese Ningqiang (NQ) and Yili (YL) breeds, we performed a transcriptome analysis using RNA sequencing (RNA-seq) technology. NQ is a pony breed with an average adult height at the withers of less than 106 cm, whereas that of YL is around 148 cm. Based on 813 million high-quality reads and stringent quality control procedures, 3011 transcripts coding for 1464 placental lncRNAs were identified and mapped to the horse reference genome. We found 107 differentially expressed lncRNAs (DELs) between NQ and YL, including 68 up-regulated and 39 down-regulated DELs in YL. Six (TBX3, CACNA1F, EDN3, KAT5, ZNF281, TMED2, and TGFB1) out of the 233 genes targeted by DELs were identified as being involved in limb development, skeletal myoblast differentiation, and embryo development. Two DELs were predicted to target the TBX3 gene, which was found to be under strong selection and associated with small body size in the Chinese Debao pony breed. This finding suggests the potential functional significance of placental lncRNAs in regulating horse body size.

Clinical, cytogenetic and molecular genetic characterization of a tandem fusion translocation in a male Holstein cattle with congenital hypospadias and a ventricular septal defect

Hypospadias, disorder of sex development (DSD), is a sporadic congenital abnormality of the genital region in male ruminants, which is characterized by a non-fused urethra during fetal development. Detailed clinical examination classified the hypospadias phenotype of a male Holstein calf studied here as the perineal type. In combined use of cytogenetic analysis and whole genome sequencing, a non-mosaic, pseudo-monosomy 59, XY + tan(18;27) was detected. This chromosomal aberration had its origin in a tandem fusion translocation of the bovine autosomes (BTA) 18 and 27 with an accompanying loss of genomic sequences mainly in the distal end of BTA 18 and the proximal end of BTA 27. The resulting phenotype included hypospadias, growth retardation and ventricular septal defect.

Whole-Genome Signatures of Selection in Sport Horses Revealed Selection Footprints Related to Musculoskeletal System Development Processes

Selective breeding has led to gradual changes at the genome level of horses. Deciphering selective pressure patterns is progressive to understand how breeding strategies have shaped the sport horse genome; although, little is known about the genomic regions under selective pressures in sport horse breeds. The major goal of this study was to shed light on genomic regions and biological pathways under selective pressures in sport horses. In this study, whole-genome sequences of 16 modern sport and 35 non-sport horses were used to investigate the genomic selective signals of sport performance, by employing fixation index, nucleotide diversity, and Tajima's D approaches. A total number of 49 shared genes were identified using these approaches. The functional enrichment analysis for candidate genes revealed novel significant biological processes related to musculoskeletal system development, such as limb development and morphogenesis, having been targeted by selection in sport breeds.

Ten years of the horse reference genome: insights into equine biology, domestication and population dynamics in the post-genome era

The horse reference genome from the Thoroughbred mare Twilight has been available for a decade and, together with advances in genomics technologies, has led to unparalleled developments in equine genomics. At the core of this progress is the continuing improvement of the quality, contiguity and completeness of the reference genome, and its functional annotation. Recent achievements include the release of the next version of the reference genome (EquCab3.0) and generation of a reference sequence for the Y chromosome. Horse satellite‐free centromeres provide unique models for mammalian centromere research. Despite extremely low genetic diversity of the Y chromosome, it has been possible to trace patrilines of breeds and pedigrees and show that Y variation was lost in the past approximately 2300 years owing to selective breeding. The high‐quality reference genome has led to the development of three different SNP arrays and WGSs of almost 2000 modern individual horses. The collection of WGS of hundreds of ancient horses is unique and not available for any other domestic species. These tools and resources have led to global population studies dissecting the natural history of the species and genetic makeup and ancestry of modern breeds. Most importantly, the available tools and resources, together with the discovery of functional elements, are dissecting molecular causes of a growing number of Mendelian and complex traits. The improved understanding of molecular underpinnings of various traits continues to benefit the health and performance of the horse whereas also serving as a model for complex disease across species.