Trio-binning of a hinny refines the comparative organization of the horse and donkey X chromosomes and reveals novel species-specific features.

We generated single haplotype assemblies from a hinny hybrid which significantly improved the gapless contiguity for horse and donkey autosomal genomes and the X chromosomes. We added over 15 Mb of missing sequence to both X chromosomes, 60 Mb to donkey autosomes and corrected numerous errors in donkey and some in horse reference genomes. We resolved functionally important X-linked repeats: the DXZ4 macrosatellite and ampliconic Equine Testis Specific Transcript Y7 (ETSTY7). We pinpointed the location of the pseudoautosomal boundaries (PAB) and determined the size of the horse (1.8 Mb) and donkey (1.88 Mb) pseudoautosomal regions (PARs). We discovered distinct differences in horse and donkey PABs: a testis-expressed gene, XKR3Y, spans horse PAB with exons1-2 located in Y and exon3 in the X-Y PAR, whereas the donkey XKR3Y is Y-specific. DXZ4 had a similar ~ 8 kb monomer in both species with 10 copies in horse and 20 in donkey. We assigned hundreds of copies of ETSTY7, a sequence horizontally transferred from Parascaris and massively amplified in equids, to horse and donkey X chromosomes and three autosomes. The findings and products contribute to molecular studies of equid biology and advance research on X-linked conditions, sex chromosome regulation and evolution in equids.

Refining the evolutionary tree of the horse Y chromosome.

The Y chromosome carries information about the demography of paternal lineages, and thus, can prove invaluable for retracing both the evolutionary trajectory of wild animals and the breeding history of domesticates. In horses, the Y chromosome shows a limited, but highly informative, sequence diversity, supporting the increasing breeding influence of Oriental lineages during the last 1500 years. Here, we augment the primary horse Y-phylogeny, which is currently mainly based on modern horse breeds of economic interest, with haplotypes (HT) segregating in remote horse populations around the world. We analyze target enriched sequencing data of 5 Mb of the Y chromosome from 76 domestic males, together with 89 whole genome sequenced domestic males and five Przewalski's horses from previous studies. The resulting phylogeny comprises 153 HTs defined by 2966 variants and offers unprecedented resolution into the history of horse paternal lineages. It reveals the presence of a remarkable number of previously unknown haplogroups in Mongolian horses and insular populations. Phylogenetic placement of HTs retrieved from 163 archaeological specimens further indicates that most of the present-day Y-chromosomal variation evolved after the domestication process that started around 4200 years ago in the Western Eurasian steppes. Our comprehensive phylogeny significantly reduces ascertainment bias and constitutes a robust evolutionary framework for analyzing horse population dynamics and diversity.

Copy number variation of horse Y chromosome genes in normal equine populations and in horses with abnormal sex development and subfertility: relationship of copy number variations with Y haplogroups.

Structural rearrangements like copy number variations in the male-specific Y chromosome have been associated with male fertility phenotypes in human and mouse but have been sparsely studied in other mammalian species. Here, we designed digital droplet PCR assays for 7 horse male-specific Y chromosome multicopy genes and SRY and evaluated their absolute copy numbers in 209 normal male horses of 22 breeds, 73 XY horses with disorders of sex development and/or infertility, 5 Przewalski's horses and 2 kulans. This established baseline copy number for these genes in horses. The TSPY gene showed the highest copy number and was the most copy number variable between individuals and breeds. SRY was a single-copy gene in most horses but had 2-3 copies in some indigenous breeds. Since SRY is flanked by 2 copies of RBMY, their copy number variations were interrelated and may lead to SRY-negative XY disorders of sex development. The Przewalski's horse and kulan had 1 copy of SRY and RBMY. TSPY and ETSTY2 showed significant copy number variations between cryptorchid and normal males (P < 0.05). No significant copy number variations were observed in subfertile/infertile males. Notably, copy number of TSPY and ETSTY5 differed between successive male generations and between cloned horses, indicating germline and somatic mechanisms for copy number variations. We observed no correlation between male-specific Y chromosome gene copy number variations and male-specific Y chromosome haplotypes. We conclude that the ampliconic male-specific Y chromosome reference assembly has deficiencies and further studies with an improved male-specific Y chromosome assembly are needed to determine selective constraints over horse male-specific Y chromosome gene copy number and their relation to stallion reproduction and male biology.

Case Report: Disorder of Sexual Development in a Chinese Crested Dog With XX/XY Leukocyte Chimerism and Mixed Cell Testicular Tumors.

A 10-year-old intact female Chinese Crested dog was presented for evaluation and further diagnostics due to persistent symptoms of vulvar swelling, vaginal discharge, and an 8-year history of acyclicity. At presentation, generalized hyperpigmentation and truncal alopecia were identified, with no aberrations of the female phenotype. Vaginal cytology confirmed the influence of estrogen at multiple veterinary visits, and hormonal screening of progesterone and anti-Mullerian hormone indicated gonadal presence. Based on findings from abdominal laparotomy and gonadectomy, the tissue was submitted for histopathology. Histopathologic evaluation identified the gonads to be abnormal testes containing multiple Sertoli and interstitial (Leydig) cell tumors. The histopathologic diagnosis of testes and concurrent normal external female phenotype in the patient lead to a diagnosis of a disorder of sexual development (DSD). Karyotype evaluation by conventional and molecular analysis revealed a two cell line chimeric pattern of 78,XX (80%) and 78,XY (20%) among blood leukocytes, as well as a positive PCR test for the Y-linked SRY gene. Cytogenetic analysis of skin fibroblasts revealed the presence of 78,XX cells exclusively, and PCR tests for the Y-linked SRY gene were negative in the hair and skin samples. These results are consistent with an XX/XY blood chimerism. This is one of the few case reports of a canine with the diagnosis of leukocyte chimerism with normal female phenotypic external genitalia. This case illustrates a distinct presentation for hormonally active Sertoli cell tumorigenesis and demonstrates surgery as a curative treatment option for clinically affected patients.

The role of impaired acrosomal exocytosis (IAE) in stallion subfertility: A retrospective analysis of the clinical condition, and an update on its diagnosis by high throughput technologies.

Acrosomal dysfunction has been considered as a cause of subfertility in males of different species, including stallions. A subset of subfertile stallions with acrosomal dysfunction is unique because they have normal sperm quality (motility, morphology, viability, and DNA quality). The current work aims to describe the clinical characteristics of subfertile stallions that were diagnosed with Impaired Acrosomal Exocytosis (IAE) by using two high throughput methods: flow cytometry and molecular genetic analysis, and to identify the prevalence of subfertility due to IAE in stallions evaluated at Texas A&M University. Clinical data from 1,128 stallions evaluated during 17 years at a Veterinary Teaching Hospital was retrospectively analyzed. Only stallions with a history of subfertility not explained following a breeding soundness examination and/or conventional semen analysis, were included. For those stallions, the acrosomal exocytosis test (AE test), in which sperm is incubated at 37 °C for up to 2 h in the presence of the calcium ionophore A23187, was used to determine IAE. The difference in AE-Rate (AE-Diff) between each pair of fertile control stallion and subfertile stallion was categorized as: Normal: AE-Diff < 14%; Questionable: AE-Diff 15-29%; Abnormal: AE-Diff > 30%. In selected cases, blood or hair was procured for identification of the susceptibility genotype for IAE, A/A-A/A, in the FKBP6 gene, exon 5. Twenty-one (21) stallions (1.86% total population analyzed) had reduced fertility despite having acceptable sperm quality. Sperm from these stallions were subjected to the AE Test. Of these, 8 stallions had reduced sperm AE-rate, based on the AE Test (8/21; 38.1%). Subsequently, blood or hair samples from these 8 stallions which had either questionable (AE-Diff 15 - 29%; n = 5) or abnormal (AE-Diff > 30%; n = 3) responses to the AE Test were analyzed for the susceptibility genotype for IAE, A/A-A/A (FKBP6 gene, exon 5). Seven out of the eight (7/8) stallions carried this susceptibility genotype. All of these were Thoroughbreds. After 2 h of incubation, the viability in fertile stallion sperm was lower than in A/A-A/A stallions (4% vs. 25%, respectively; P < 0.05), while the AE-rate was higher for fertile than for A/A-A/A stallions (85% vs. 56%, respectively; P < 0.05). The use of two high throughput tests (i.e., flow cytometry and molecular genetic analysis) may complement each other in the diagnosis of IAE in breeding stallions. In this study, 5/7 subfertile stallions diagnosed with the IAE susceptibility genotype would have been diagnosed as normal with the AE Test. This study introduces a subset of stallions with the IAE genotype with fertility higher than has been previously reported (i.e., <15% per-cycle pregnancy rate), suggesting that IAE manifests as a broader range of subfertility.

Molecular Cytogenetic and Y Copy Number Analysis of a Reciprocal ECAY-ECA13 Translocation in a Stallion with Complete Meiotic Arrest.

We present a detailed molecular cytogenetic analysis of a reciprocal translocation between horse (ECA) chromosomes Y and 13 in a Friesian stallion with complete meiotic arrest and azoospermia. We use dual-color fluorescence in situ hybridization with select ECAY and ECA13 markers and show that the translocation breakpoint in ECAY is in the multicopy region and in ECA13, at the centromere. One resulting derivative chromosome, Y;13p, comprises of ECAY heterochromatin (ETSTY7 array), a small single copy and partial Y multicopy region, and ECA13p. Another derivative chromosome 13q;Y comprises of ECA13q and most of the single copy ECAY, the pseudoautosomal region and a small part of the Y multicopy region. A copy number (CN) analysis of select ECAY multicopy genes shows that the Friesian stallion has significantly (p < 0.05) reduced CNs of TSPY, ETSTY1, and ETSTY5, suggesting that the translocation may not be completely balanced, and genetic material is lost. We discuss likely meiotic behavior of abnormal chromosomes and theorize about the possible effect of the aberration on Y regulation and the progression of meiosis. The study adds a unique case to equine clinical cytogenetics and contributes to understanding the role of the Y chromosome in male meiosis.

An 8.22 Mb Assembly and Annotation of the Alpaca (Vicugna pacos) Y Chromosome.

The unique evolutionary dynamics and complex structure make the Y chromosome the most diverse and least understood region in the mammalian genome, despite its undisputable role in sex determination, development, and male fertility. Here we present the first contig-level annotated draft assembly for the alpaca (Vicugna pacos) Y chromosome based on hybrid assembly of short- and long-read sequence data of flow-sorted Y. The latter was also used for cDNA selection providing Y-enriched testis transcriptome for annotation. The final assembly of 8.22 Mb comprised 4.5 Mb of male specific Y (MSY) and 3.7 Mb of the pseudoautosomal region. In MSY, we annotated 15 X-degenerate genes and two novel transcripts, but no transposed sequences. Two MSY genes, HSFY and RBMY, are multicopy. The pseudoautosomal boundary is located between SHROOM2 and HSFY. Comparative analysis shows that the small and cytogenetically distinct alpaca Y shares most of MSY sequences with the larger dromedary and Bactrian camel Y chromosomes. Most of alpaca X-degenerate genes are also shared with other mammalian MSYs, though WWC3Y is Y-specific only in alpaca/camels and the horse. The partial alpaca Y assembly is a starting point for further expansion and will have applications in the study of camelid populations and male biology.

Two Novel Cases of Autosomal Translocations in the Horse: Warmblood Family Segregating t(4;30) and a Cloned Arabian with a de novo t(12;25).

We report 2 novel autosomal translocations in the horse. In Case 1, a breeding stallion with a balanced t(4p;30) had produced normal foals and those with congenital abnormalities. Of his 9 phenotypically normal offspring, 4 had normal karyotypes, 4 had balanced t(4p;30), and 1 carried an unbalanced translocation with tertiary trisomy of 4p. We argue that unbalanced forms of t(4p;30) are more tolerated and result in viable congenital abnormalities, without causing embryonic death like all other known equine autosomal translocations. In Case 2, two stallions produced by somatic cell nuclear transfer from the same donor were karyotyped because of fertility issues. A balanced translocation t(12q;25) was found in one, but not in the other clone. The findings underscore the importance of routine cytogenetic screening of breeding animals and animals produced by assisted reproductive technologies. These cases will contribute to molecular studies of translocation breakpoints and their genetic consequences in the horse.

Characterization of A Homozygous Deletion of Steroid Hormone Biosynthesis Genes in Horse Chromosome 29 as A Risk Factor for Disorders of Sex Development and Reproduction.

Disorders of sex development (DSD) and reproduction are not uncommon among horses, though knowledge about their molecular causes is sparse. Here we characterized a ~200 kb homozygous deletion in chromosome 29 at 29.7-29.9 Mb. The region contains AKR1C genes which function as ketosteroid reductases in steroid hormone biosynthesis, including androgens and estrogens. Mutations in AKR1C genes are associated with human DSDs. Deletion boundaries, sequence properties and gene content were studied by PCR and whole genome sequencing of select deletion homozygotes and control animals. Deletion analysis by PCR in 940 horses, including 622 with DSDs and reproductive problems and 318 phenotypically normal controls, detected 67 deletion homozygotes of which 79% were developmentally or reproductively abnormal. Altogether, 8-9% of all abnormal horses were homozygous for the deletion, with the highest incidence (9.4%) among cryptorchids. The deletion was found in ~4% of our phenotypically normal cohort, ~1% of global warmblood horses and ponies, and ~7% of draught breeds of general horse population as retrieved from published data. Based on the abnormal phenotype of the carriers, the functionally relevant gene content, and the low incidence in general population, we consider the deletion in chromosome 29 as a risk factor for equine DSDs and reproductive disorders.

Azoospermia and Y Chromosome-Autosome Translocation in a Friesian Stallion.

This case report describes spermatogenic arrest and azoospermia in a stallion with a unique Y chromosome-autosome translocation. Clinical diagnosis of azoospermia was based on history of infertility and evaluation of ejaculates collected for artificial insemination. Clinical and ultrasonographic evaluation of the external and internal genitalia did not reveal any abnormalities except for smaller than normal testicular size. Azoospermia of testicular origin was confirmed by determining alkaline phosphatase concentration in semen. Histological evaluation of testicular tissue after castration confirmed early spermatogenic arrest. Cytogenetic evaluation showed the presence of translocation between the Y chromosome and chromosome 13. To the authors' knowledge, this is the first case of azoospermia with a cytogenetically detected Y chromosome abnormality, suggesting that the horse Y chromosome may carry sequences critical for normal spermatogenesis.

Population Genetic Analysis of the Estonian Native Horse Suggests Diverse and Distinct Genetics, Ancient Origin and Contribution from Unique Patrilines.

The Estonian Native Horse (ENH) is a medium-size pony found mainly in the western islands of Estonia and is well-adapted to the harsh northern climate and poor pastures. The ancestry of the ENH is debated, including alleged claims about direct descendance from the extinct Tarpan. Here we conducted a detailed analysis of the genetic makeup and relationships of the ENH based on the genotypes of 15 autosomal short tandem repeats (STRs), 18 Y chromosomal single nucleotide polymorphisms (SNPs), mitochondrial D-loop sequence and lateral gait allele in DMRT3. The study encompassed 2890 horses of 61 breeds, including 33 ENHs. We show that the expected and observed genetic diversities of the ENH are among the highest within 52 global breeds, and the highest among 8 related Northern European ponies. The genetically closest breeds to the ENH are the Finn Horse, and the geographically more distant primitive Hucul and Konik. ENH matrilines are diverse and relate to draught and Pontic-Caspian breeds. ENH patrilines relate to draught breeds, and to a unique haplogroup not described before. None of the 33 ENHs carried the "gait" mutation, but the mutation was found in 2 Huculs. The study demonstrates that the ENH is a genetically distinct and diverse breed of ancient origin with no notable pressure of selective breeding.

Comparative FISH-Mapping of MC1R, ASIP, and TYRP1 in New and Old World Camelids and Association Analysis With Coat Color Phenotypes in the Dromedary (Camelus dromedarius).

Melanocortin 1 receptor (MC1R), the agouti signaling protein (ASIP), and tyrosinase related protein 1 (TYRP1) are among the major regulators of pigmentation in mammals. Recently, MC1R and ASIP sequence variants were associated with white and black/dark brown coat colors, respectively, in the dromedary. Here we confirmed this association by independent sequencing and mutation discovery of MC1R and ASIP coding regions and by TaqMan genotyping in 188 dromedaries from Saudi Arabia and United States, including 38 black, 53 white, and 97 beige/brown/red animals. We showed that heterozygosity for a missense mutation c.901C > T in MC1R is sufficient for the white coat color suggesting a possible dominant negative effect. Likewise, we confirmed that the majority of black dromedaries were homozygous for a frameshift mutation in ASIP exon 2, except for 4 animals, which were heterozygous. In search for additional mutations underlying the black color, we identified another frameshift mutation in ASIP exon 4 and 6 new variants in MC1R including a significantly associated SNP in 3'UTR. In pursuit of sequence variants that may modify dromedary wild-type color from dark-reddish brown to light beige, we identified 4 SNPs and one insertion in TYRP1 non-coding regions. However, none of these were associated with variations in wild-type colors. Finally, the three genes were cytogenetically mapped in New World (alpaca) and Old World (dromedary and Bactrian camel) camelids. The MC1R was assigned to chr21, ASIP to chr19 and TYRP1 to chr4 in all 3 species confirming extensive conservation of camelid karyotypes. Notably, while the locations of ASIP and TYRP1 were in agreement with human-camelid comparative map, mapping MC1R identified a new evolutionary conserved synteny segment between camelid chromosome 21 and HSA16. The findings contribute to coat color genomics and the development of molecular tests in camelids and toward the chromosome level reference assemblies of camelid genomes.

Genetic diversity and origin of the feral horses in Theodore Roosevelt National Park.

Feral horses in Theodore Roosevelt National Park (TRNP) represent an iconic era of the North Dakota Badlands. Their uncertain history raises management questions regarding origins, genetic diversity, and long-term genetic viability. Hair samples with follicles were collected from 196 horses in the Park and used to sequence the control region of mitochondrial DNA (mtDNA) and to profile 12 autosomal short tandem repeat (STR) markers. Three mtDNA haplotypes found in the TRNP horses belonged to haplogroups L and B. The control region variation was low with haplotype diversity of 0.5271, nucleotide diversity of 0.0077 and mean pairwise difference of 2.93. We sequenced one mitochondrial genome from each haplotype determined by the control region. Two complete mtDNA sequences of haplogroup L were closely related to the mtDNA of American Paint horse. The TRNP haplotype B did not have close matches in GenBank. The phylogenetic test placed this sequence in a group consisting of two horses from China, one from Yakutia, and one from Italy raising a possibility of historical transportation of horses from Siberia and East Asia to North America. Autosomal STR loci were polymorphic and indicated that the TRNP horses were distinctly different from 48 major horse breeds. Heterozygosity, mean number of alleles, and other measures of diversity indicated that TRNP herd diversity was below that observed for most other feral herds and domestic breeds. Both mtDNA and STRs demonstrated that the existing genetic data sets of horses are insufficient to determine the exact origins of the TRNP horses. However, measures of nuclear and mitochondrial diversity have elucidated management needs. It is recommended that new genetic stock be introduced and that adaptive management principles are employed to ensure that unique mitochondrial lineages are preserved and genetic diversity is increased and maintained over time.