{"paper_id":"35d84ca7-e9be-4733-8e4a-d73b59b57d45","body_text":"Mammalian disorders of sex development (DSDs) encompass a broad range of complex congenital conditions that affect sex determination, sexual differentiation, and the development and function of gonads [ 1 , 2 , 3 , 4 , 5 , 6 ]. The majority of DSDs are characterized by an atypical combination of chromosomal, gonadal and phenotypic sex and are broadly classified as XX and XY DSDs [ 5 , 7 ]. Though, in some forms of XY DSDs, such as cryptorchidism and hypospadias, the genetic and gonadal/phenotypic sex are in concordance [ 8 ]. The prevalence of DSDs in humans is approximately 1.7% of all live births [ 7 , 9 ], whereas the frequency of non-syndromic cryptorchidism alone can be as high as 8% among full-term male births depending on the geographical area [ 8 , 10 ].\nSex development is a complex multi-step process controlled by many genes acting in networks [ 4 , 5 , 11 ], which complicates the search for the molecular causes of DSDs. Only a handful of human and mouse DSDs can be attributed to point mutations or copy number variants (CNVs) in a single gene, others are thought to be oligogenic with the involvement of additional modifier genes, while molecular causes of many conditions are not known. Furthermore, variations in the same gene can result in different DSD phenotypes and similar DSD phenotypes may have different genetic causes [ 4 , 11 ].\nMultiple forms of XX and XY DSDs have also been described in horses [ 12 , 13 , 14 , 15 ] and all forms affect equine health and reproduction, as well as the economy of breeders and owners. Among equine XY DSDs, the most common is cryptorchidism with a prevalence of 2–8% among male births [ 16 ], followed by a spectrum of XY and XX DSDs, all characterized by various degrees of discrepancy between genetic, gonadal and phenotypic sex and ambiguous external sex phenotypes [ 12 , 13 , 14 , 17 ].\nKnowledge about the molecular causes of equine DSDs is sparse. For example, no genes or mutations are known for any of the XX DSDs with intersex, hermaphrodite, pseudohermaphrodite or ambiguous sex phenotypes [ 12 , 13 , 14 ]. Among XY DSDs, the  SRY -negative forms are associated with various Y chromosome deletions [ 18 ], while genetic causes of  SRY -positive XY DSDs are known only for a few cases and are heterogeneous. In three families of different breeds,  SRY -positive XY DSDs have been associated with different mutations in the androgen receptor gene resulting in androgen insensitivity syndrome [ 19 , 20 , 21 ]. Earlier studies identified a large deletion in horse ( Equus caballus , ECA) chromosome 29 (ECA29) that was found in two related Standardbred horses with  SRY -positive XY DSDs and male pseudohermaphrodite phenotypes [ 22 ]. The deletion involves steroid hormone biosynthesis genes belonging to the aldo-keto reductase family 1 C ( AKR1C ), of which mutations have been associated with XY DSDs, such as cryptorchidism, undervirilized external genitalia, dihydrotestosterone (DHT) deficiency and sex reversal conditions in humans [ 1 , 3 , 4 ].\nIn this study, we will further investigate the molecular landscape of the deletion in horse chromosome 29 by PCR, fluorescence in situ hybridization (FISH) and whole genome sequencing (WGS). We characterize the deletion breakpoints, develop a simple PCR test for the detection of deletion homozygotes, and determine their frequency in general equine populations, as well as among horses with various DSDs and reproductive disorders.\n\nProcurement of blood and hair samples for DNA isolation followed the United States Government Principles for the Utilization and Care of Vertebrate Animals Used in Testing, Research and Training. These protocols were approved as AUPs #2012-0250 and #2018-0342 CA at Texas A&M University. Owners of the Swedish samples gave permission for their animals to be used in the study and the study was approved by the Ethics Committee for Animal Experiments in Uppsala, Sweden (#C 121/14 and #5.8.18–15453/2017).\nSamples of 940 horses of over 45 diverse breeds and breed mixes were used for experiments in this study ( Supplementary Table S1 ). These included 364 horses with various DSDs of which 159 were cryptorchids; 258 phenotypically normal but reproductively or hormonally abnormal male and female horses, and 318 control horses with normal male or female phenotypes, but with limited information about fertility. Control horses included Thoroughbred mare  Twilight , the DNA donor for horse reference genome assemblies EquCab2 [ 23 ] and EquCab3 [ 24 ], and her half-sibling Thoroughbred stallion  Bravo , the DNA donor of bacterial artificial chromosome (BAC) library CHORI 241 ( https://bacpacresources.org/ ). Both horses are normal with proven fertility. The majority of samples available from the depository of Texas A&M Molecular Cytogenetics laboratory also had karyotype data.\nIn addition, whole genome sequence (WGS) and/or high density 670K SNP genotyping data from publications [ 25 , 26 , 27 ] were used for population analysis. Horses from these cohorts were not specifically phenotyped for sex development or fertility.\nGenomic DNA was isolated from hair or peripheral blood using Gentra Puregene Blood Kit (Qiagen) following the manufacturer’s protocol. The DNA was checked for quality and quantity with Nanodrop 2000 spectrophotometer (Thermo Scientific).\nThe IDs and sequence coordinates of 4 partially overlapping genomic clones spanning the deletion in ECA29 and a control clone outside the deletion ( Table 1 ) were retrieved from NCBI ( https://www.ncbi.nlm.nih.gov/genome/?term=horse ). The clones were picked from the CHORI 241 BAC library ( https://bacpacresources.org/ ). High molecular weight BAC DNA was isolated with High Pure Plasmid Isolation Kit (Roche Diagnostics) per manufacturer instructions and checked for molecular weight, quality and quantity with Agilent 2100 TapeStation system.\nThe ECA29 deletion was analyzed by qualitative and quantitative PCR. Primers for sequence tagged sites (STSs) around putative deletion breakpoints and for genes inside the deletion ( Supplementary Table S2 ) were designed using horse reference genome EquCab3 [ 24 ] (NCBI:  https://www.ncbi.nlm.nih.gov/genome/?term=horse ) and Primer3 software [ 28 ]. Each primer was validated by in silico PCR in UCSC ( http://genome.ucsc.edu/ ) to produce a single product unique for ECA29.\nQualitative PCR was performed with JumpStart Taq DNA polymerase (Sigma) and the results were analyzed by agarose gel electrophoresis. Qualitative PCR with DELin1 primers to detect horses with homozygous deletion were carried out in duplicate reactions and duplicate experiments.\nQuantitative PCR (qPCR) experiments were performed with LightCycler® 480 (Roche Diagnostics) in triplicate assays and triplicate 20 µL reactions, each containing 50 ng of template DNA, 10 μM primers, and the Hot FirePol® EvaGreen® qPCR Mix Plus (Solis BioDyne). Relative copy numbers were determined in comparison to the reference sample (Thoroughbred  Bravo ) and normalized to an autosomal reference gene  GAPDH . Copy numbers were quantified by the comparative C T  method (2 ΔΔ Ct) [ 29 , 30 ] with  p  < 0.05 as a cut-off threshold for statistical significance.\nThe DNA (1 µg) of BACs 23N13 spanning the deletion and a control BAC 76H13 outside the deletion ( Table 1 ) were labeled with biotin-16-dUTP or digoxigenin-11-dUTP using Biotin- or DIG-Nick Translation Mix (Roche Diagnostics). Labeled probes were hybridized to metaphase chromosomes of horses where analysis by PCR suggested homozygous or heterozygous deletion, and to control horses without deletion. FISH was done following standard protocols [ 31 ]. Images for a minimum of 20 metaphase cells were captured for each experiment and analyzed with a Zeiss Axioplan2 fluorescent microscope equipped with Isis v5.2 (MetaSystems GmbH) software.\nWhole genomes of two horses with homozygous deletion in ECA29 and XY DSD phenotypes (H369, a male pseudohermaphrodite and H354, a bilateral cryptorchid;  Supplementary Table S1 ) were sequenced using 1x100 bp libraries prepared with TruSeq DNA PCR-Free Sample Preparation Kit (Illumina), and an Illumina HiSeq2500 platform producing approximately 6-10X average coverage per animal. All data were aligned to EquCab3 using bwa [ 32 ], sorted with samtools 1.9 [ 33 ], and duplicates marked with sambamba [ 34 ]. Publicly available control data from Jagannathan et al. 2019 (QH070 and HN001 for visualization) were mapped in the same fashion.\nThe 4 BAC clones ( Table 1 ) spanning the deletion were sequenced by Pacific Biosciences (PacBio) single molecule real time (SMRT) technology. Sequencing libraries with a molecule size range of 20 kb were generated using the RS DNA Template Preparation Kit (Pacific Biosciences), loaded on a SMRT cell and sequenced for long-reads with C2 chemistry on the PacBio RS [ 35 ], with an estimated yield of approximately 5 Gb data per SMRT cell. Resulting filtered subreads were assembled using canu 1.8 [ 36 ] for each individual BAC, as well as a separate assembly for all data. For reference based alignment, PacBio subreads were aligned to EquCab3 using minimap2.1 [ 37 ]. Furthermore, mapping of individual contigs generated by canu were aligned using minimap for visualization.\nThe content and functions of genes affected by the deletion were retrieved from NCBI ( http://genome.ucsc.edu/ ), UCSC ( http://genome.ucsc.edu/ ) and Ensembl ( http://www.ensembl.org/index.html ) Genome Browsers, and GeneCards ( http://www.genecards.org/ ).\nProtein sequences for horse, human and mouse  AKR1C  genes were retrieved from NCBI ( http://genome.ucsc.edu/ ) and aligned with the Clustal Omega multiple sequence alignment tool ( https://www.ebi.ac.uk/Tools/msa/clustalo/ ).\n\nThe initial discovery of the deletion in ECA29 was by array comparative genomic hybridization (aCGH), which identified a ~200 kb homozygous deletion in two related horses (H348 and H369;  Supplementary Table S1 ) with XY DSDs [ 22 ]. To determine, whether or not the deletion is an isolated event we used the published [ 22 ] deletion-specific primers ( DelIn1 ,  Supplementary Table S2 ) to screen a small selection of horses with XY DSDs by regular PCR. As a result, additional animals that were homozygous for the deleted sequence were identified ( Figure 1 ), indicating that the deletion is not limited to the two discovery horses and requires more detailed analysis.\nIn order to demarcate deletion breakpoints, a set of primers were designed for PCR walking around the putative deletion start (primers S1-S11) and end sites (primers E1–E9) ( Supplementary Table S2 ,  Figure 2 ), which were determined based on aCGH data [ 22 ]. Walking by PCR using DNA from horses with known homozygous deletion (D/D) and controls without deletion suggested that the deletion starts in a 1019 bp region at chr29: 29,740,911-29,741,930 between primers S7 and S8, and ends in a 95 bp region at chr29:29,944,184-29,944,279 between primers E8 and E9 ( Supplementary Table S2 ,  Figure 2 ). However, demarcation of the end of the deletion by PCR remained dubious because sequences corresponding to the amplicons of primers E10–E12 showed multiple hits in EquCab3. Furthermore, primers E4 did not amplify from D/D or controls ( Supplementary Table S2 ), suggesting that the region is complex and the reference assembly has gaps and/or errors in this region.\nTo determine the deletion boundaries more accurately, we generated WGS for two D/D horses: H354 and H369 ( Figure 1 ,  Supplementary Table S1 ). Sequence alignment with EquCab3 and published WGSs of 2 normal horses [ 27 ] showed that the deletion in ECA29 starts at position 29,741,143 bp in H354 and at 29,741,147 bp in H369, near a documented repetitive region in EquCab3, and consistent with the region determined by PCR ( Figure 3 ). In horse H354, the putative end was around 29,944,505 bp and in horse H369, at 29,952,177 bp, thus over 7.6 kb more distal. Both putative deletion end sites determined by WGS ( Figure 3 ) extended further compared to the one determined by PCR ( Figure 2 ). Based on WGS data, we estimated the approximate size of the deletion to be 200–210 kb. However, due to multi-mapping short read alignment from 29,944,163 to 29,952,267 roughly 20–30X average depth, the precise end breakpoint was not precisely determined. This ~8,100 bp region indicates a putative tandem duplication or complex rearrangement at the end of the deletion, potentially inducing the structural change. No duplication is evident in 88 control horses from the WGS data obtained from Jagannathan et al. 2019, thus genotypically normal horses do not possess the putative additional structural variant(s) associated with the deletion. Resolution is unattainable with short-read data, and should be explored using newer PacBio or Oxford Nanopore long-read data.\nIn an attempt to improve sequence assembly in this region, we de novo assembled long-read PacBio sequence data of 4 partially overlapping BACs spanning the deletion ( Table 1 ,  Figure 3 ). Previously, we have shown that Thoroughbreds  Bravo  and  Twilight , the DNA donors of the BAC library and genome assembly, respectively, differ by just 3 CNVs across the genome with no CNVs in ECA29 [ 22 ]. Thus, the BAC sequences should closely resemble the reference genome and can be used to improve the assembly in this region. Independent assembly of all four BACs produced one long contig with high read coverage (>3000) each. Assembly of all data together produced one 338,300 bp contig spanning the deletion. This confirmed that there is no evidence for duplication in the region, other than in D/D individuals ( Figure 3 ).\nTaken together, both PCR analysis and WGS confirmed the presence of a large, approximately 200 kb deletion in some horses with XY DSDs. While precise demarcation of deletion boundaries remained tentative due to sequence complexity in the region, the most proximal and most distal points of the deletion by WGS were chr29:29,741,143–29,952,177 spanning 211,034 bp.\nThe deleted region chr29:29,741,143–29,952,177 contains 2 predicted protein-coding genes and one RNA gene by Ensembl ( Figure 3 ). However, the Ensembl gene prediction ENSECAG00000021292 corresponds to 5 non-overlapping fragments of the same gene in NCBI ( Table 2 ), suggesting that the region contains members of a gene family which have not yet been accurately identified and annotated. Despite this, all predicted coding genes in the deleted region belong to aldo-keto reductase family 1, known to play essential roles in the metabolism of all steroid hormones [ 38 ], including the production of extra-testicular androgens [ 39 ]. Deletion of all genes in this region in select XY DSD horses with homozygous deletion was confirmed by PCR with gene-specific primers E661, E668, E458, E332 and ED2 ( Figure 3 ;  Supplementary Table S2 ).\nOnce approximate deletion start and end sites had been determined ( Figure 2  and  Figure 3 ), conventional qualitative PCR with primers designed inside the deletion ( Supplementary Table S2 ) was a fast and efficient approach for the detection of horses homozygous for the deletion ( Figure 1  and  Figure 2 ). However, because of the large size (~200 kb) of the deleted region and ambiguity in defining deletion boundaries, we were not able to detect heterozygous carriers by qualitative PCR. Therefore, we designed gene-specific primers G_2879 ( Figure 3 ,  Supplementary Table S2 ) for quantitative PCR (qPCR) and analyzed selected 12 horses with DSDs (including 8 cryptorchids) and 11 normal control horses for possible heterozygotes ( Figure 4 ). As expected, previously known D/D horses had zero copies. Other horses, both with DSDs and normal controls, showed considerable CNV in a range from 20 to 120 relative copies, which is consistent with previous CNV studies in the region [ 22 ].\nThe broad range of variation in qPCR results made calling deletion heterozygotes ambiguous and we further validated the results by FISH. We selected two horses with relative copy numbers below 50 ( Bravo —normal fertile male and H657—an XX intersex), one horse with > 50 copies (H395), and one with homozygous deletion and zero copies (H369) as a negative control. FISH was done with BAC 23N13, which spans the deletion ( Figure 3 ) and a control BAC 76H13, which maps to the proximal part of ECA29 ( Table 1 ). The results showed that H395 (>50 copies) and  Bravo  (<50 copies) did not have the deletion, horse H657 (<50 copies) with similar copy numbers by qPCR as  Bravo  was heterozygous, and horse H369, as expected, homozygous for the deletion ( Figure 5 ). These results indicate that FISH analysis can clearly discriminate between horses with no deletion, deletion heterozygotes and horses with homozygous deletion. However, FISH is a method of too low throughput for routine screening of large number of horses. On the other hand, relative copy number estimation by qPCR to discriminate between horses with no deletion and those with heterozygous deletion remained ambiguous. Furthermore, a broad range in copy numbers as detected by qPCR analysis ( Figure 4 ) suggests that this CNV may be multi-allelic like several common complex CNVs described in human diseases [ 40 ]. This also implies that qPCR-based analysis can confidently detect zero copies in homozygous deletion, maybe also 1 copy in simple heterozygotes, but will provide no information about the distribution of >2 copies on homologous chromosomes.\nBecause the first identified homozygotes for ECA29 deletion were horses with DSDs [ 22 ] ( Figure 1 ), and because the deletion involved genes necessary for steroid hormone biosynthesis ( Table 2 ) and sex development [ 38 ], we further investigated by PCR the frequency of deletion homozygotes among 622 developmentally or reproductively abnormal horses. These included horses with various XX and XY DSDs and other congenital disorders, horses with normal male or female sex development but with unexplained infertility or subfertility, and horses producing offspring with DSDs. We also tested 318 phenotypically normal male and female horses, though with limited information about their fertility (63 with proven fertility; 255 with no information on fertility). Details about all 940 individual horses are presented in  Supplementary Table S1 .\nAltogether, we identified 67 horses with homozygous deletion of which 53 horses (79%) were diagnosed with DSDs or fertility problems. The frequency of homozygous deletion was almost two times higher in the developmentally and reproductively abnormal cohort (8.1%) compared to controls (4.7%), with the highest incidence found in cryptorchid horses (9.4%) ( Table 3 ). We did not observe any bias of the deletion between bi- or unilateral cryptorchids. There was also no difference between the occurrence of the deletion across different breeds and D/D individuals were present in 17 different breeds and breed mixes out of 45 breeds involved.\nIn order to estimate the frequency of the homozygous deletion in the general equine population, we analyzed publicly available WGS and 670K SNP array data for the region chr29: 29,741,143—29,952,177 ( Table 4 ). The frequency of homozygous deletion carriers varies considerably between the studies being 0–1% among global warmblood and pony breeds [ 26 , 27 ], but reaching over 7% in a recent study, which is strongly biased towards draught horses [ 25 ]. Collectively, the frequency of homozygous deletion in general equine population is lower (5.6%) than that among horses with DSDs and fertility problems (8.1%). Because the sequence assembly of the region has ambiguities ( Figure 3 ), and the region likely represents a complex multi-allelic CNV, no efforts were made to estimate the number of putative deletion heterozygotes.\n\nHere we characterized a ~200 kb CNV region (CNVR) in horse chromosome 29 at 29.7–29.9 Mb, which contains genes essential for steroid hormone biosynthesis and was deleted from both homologs in 8–9% of horses with reproductive problems and DSDs. On the other hand, 79% (53/67) of the horses carrying the homozygous deletion were developmentally or reproductively abnormal.\nThe deleted region in ECA29 is embedded in a larger common CNVR that has been found in the general equine population of diverse origins by several previous studies [ 22 , 25 , 41 , 42 ]. The most extensive recent genome-wide analysis of 939 CNVRs in 1755 horses of 8 European breeds, detected the larger 714 kb common deletion-duplication CNVR (CNVR915) in 1204 (68.6%) horses [ 25 ]. As a notable contrast, the homozygous deletion of the smaller 200 kb region is rare and was found only in about 4.7% of phenotypically normal horses used in this study, and in 0% to 7.6% of general equine population from other studies ( Table 3  and  Table 4 ). Based on published data, the frequency of the homozygous deletion is extremely low (0% to 1%) among 26 global warmblood and pony breeds [ 26 , 27 ]. For example, the larger CNVR is present in 84% (216/256) of Exmoor ponies studied for CNVs [ 25 ], but none of these carry the homozygous deletion of the smaller region [ 26 ] ( Table 4 ). At the same time, the frequency of D/D individuals seems to be considerably higher (7.6%) among European draught breeds [ 25 ] and it would be of future interest to compare the genotyping data with reproductive phenotypes in these breeds.\nGiven that 79% of the carriers of the homozygous deletion in ECA29, as determined in this study, were developmentally or reproductively abnormal, the content of genes in the region is of particular interest. According to gene predictions, both deleted protein-coding genes belong to the aldo-keto reductase gene family 1C ( AKR1C ) ( Figure 3 ,  Table 2 ) members of which carry out various functions in steroid hormone metabolism. More specifically, AKR1C enzymes (AKR1C1–AKR1C4) function in vivo as 3-keto-, 17-keto- and 20-ketosteroid reductases and metabolize a broad spectrum of natural and synthetic therapeutic steroids [ 6 ]. They are involved both in androgen and estrogen metabolism by controlling concentrations of active androgens and estrogens, and regulating occupancy and trans-activation of androgen and estrogen receptors [ 38 ]. All pathways to the potent androgens testosterone and dihydrotestosterone (DHT) proceed through family members  AKR1C2-AKR1C4  [ 38 ], while  AKR1C2  and  AKR1C4  also participate in ‘the backdoor pathway’ that leads to DHT synthesis without testosterone intermediate [ 43 , 44 , 45 ].\nAltered expression of individual  AKR1C  genes is related to the development of prostate and breast cancer, and endometriosis in humans due to changes in progesterone, estrogen and androgen metabolism [ 38 , 46 ]. Mutations in  AKR1C  genes are associated with human XY DSDs, including cryptorchidism, sex reversal and defective androgen synthesis phenotypes [ 3 , 4 , 38 , 44 , 47 ].\nMouse models are less informative because there are no one-to-one mouse orthologs for human or horse  AKR1C  genes (Ensembl:  http://www.ensembl.org/index.html ). While amino acid sequence identity between human  AKR1C1-AKR1C4  paralogs is 82–99% ( Supplementary Table S3 ) [ 38 ], their one-to-one orthologs have been found only in primates, while orthologs in other mammals, including mouse, are many-to-many type (Ensembl:  http://www.ensembl.org/index.html ). On the one hand, this suggests rapid evolutionary divergence of  AKR1C  orthologs across species, and on the other hand, difficulties to annotate paralogs within a species. The closest murine homolog for human  AKR1C1-AKR1C4  genes is  Akr1c18  (MGI:  http://www.informatics.jax.org/marker ) sharing only 65–68% amino acid sequence identity with human and 51–65% with equine orthologs ( Supplementary Table S3 ). Nevertheless, it is noteworthy that  Akr1c18  (-/-) live knockouts are exclusively female with reduced fetal viability, low adult pregnancy rates and increased embryonic loss, whereas male knockouts are likely not viable [ 48 , 49 ].\nThe human and mouse data strongly suggest that  AKR1C  genes are necessary for normal sex development and reproduction and their mutations have consequences. Furthermore, in the CNV map of the human genome, the  AKR1C  genes are not among the ~100 “expendable” CNV genes which can be homozygously deleted with no deleterious effect [ 50 ]. This is consistent with the findings of the present study where the majority (79%) of horses with the homozygous deletion were abnormal regarding sex development and/or reproduction. Though, there is no single clear explanation for the remaining 21% of horses also carrying the homozygous deletion, in part because the available fertility data for the “normal” population were sparse ( Supplementary Table S1 ). Furthermore, the assembly of the region in ECA29 is incomplete even in the recently improved reference genome EquCab3 [ 24 ], and annotation of individual  AKR1C  family members is ambiguous ( Table 2 ). The region is structurally complex, containing CNVs, segmental duplications and  AKR1C  paralogs that share over 86% sequence identity with each other [ 6 , 51 ]. These features complicated accurate demarcation of deletion boundaries. For example, by combining PCR and WGS data, we confined the deletion between 29.7–29.9 Mb ( Figure 2  and  Figure 3 ). In EquCab3 in UCSC Genome Browser ( https://genome.ucsc.edu/ ), the cluster of  AKR1C  paralogs is located at chr29: 29.7–30.3 Mb, thus some members of the gene family may escape the deletion and compensate the loss of others. It is also possible that the size of the deletion differs between individual horses, but was not captured by PCR analysis of a short target region. However, perhaps the most plausible explanation for why not all animals with the deletion showed abnormal phenotypes, as well as why those with the deletion had a broad spectrum of abnormal phenotypes, stays in the very nature of the genetic regulation of steroid metabolism. It is well documented that among all CNV genes in humans and animals, the most abundant are G-protein coupled receptors, and immunity and steroid metabolism genes [ 22 , 50 , 52 , 53 ]. This implies that the gene networks and biochemical pathways governed by these genes are so critically important for the development, viability and adaptive plasticity that they are duplicated and backed up by alternatives in the genome [ 54 ]. We therefore hypothesize that phenotypic consequences of the deletion in ECA29 varies from case to case depending on the overall genetic background of each individual.\nFurther, in the large cohort of 622 developmentally and reproductively abnormal horses, only 8%-9% carried the homozygous deletion and phenotypes of the carriers were diverse. Phenotypically, the most uniform group were cryptorchids of which 9.4% were homozygous for the deletion ( Table 3 ), though cryptorchid phenotype too varied in a broad range. Therefore, it was not possible to statistically associate the deletion with any specific form of equine cryptorchidism or other DSD. Our observations are consistent with those in humans where studies on DSDs and disorders of reproduction are more advanced and abundant [ 1 , 2 , 3 , 4 ], but face similar challenges for uncovering the genetic causes. One is the very nature of these conditions which prevents reproduction and transmission, making almost every case a singleton. Another is the complex gene network underlying sex development and reproduction. In addition, most of these conditions have a spectrum of phenotypes. Together these factors exclude association studies based on large, well-phenotyped case-control cohorts, and require genetic studies on a case-by-case basis.\nNevertheless, based on the size and gene content of the deletion, the abnormal phenotype of the majority (79%) of deletion homozygotes, and the low frequency (~5%) of the homozygous deletion in general equine population, we consider this deletion as a possible causative, but definitely a risk factor, for some forms of equine DSDs, cryptorchidism, and reproductive disorders.\n\nWe characterized a homozygous ~200 kb deletion in ECA29 as a risk factor for equine DSDs and reproductive disorders. These conclusions were based on the following notions: i) the deletion eliminates functionally relevant  AKR1C  genes, which have been associated with DSDs in humans; ii) the frequency of the deletion was higher in abnormal horses (8–9%) compared to normal population (1–5%), and iii) the majority (79%) of horses with the homozygous deletion were phenotypically abnormal. Even though denoted as a  risk factor  because of difficulties to establish statistically sound causative genotype–phenotype relationship, the findings add to the currently very short list of genomic regions underlying equine disorders of sex development and reproduction (reviewed by [ 55 ]). A simple PCR analysis to identify ECA29 deletion homozygotes can be easily included in clinical genetic testing of problem horses, similarly to the routinely conducted test for the Y-linked  SRY  gene [ 18 ]. At the same time, for any further progress in the study of this deletion, it is important to improve phenotype descriptions for all new cases. For example, since the  AKR1C  genes are involved in steroid hormone biosynthesis, it would be important to include data on serum hormone levels. Further, the fact that only 8–9% of the 622 reproductively and developmentally abnormal horses in this study carried the ECA29 deletion also implies the need of continuing efforts to decipher the molecular causes of equine DSDs and reproductive disorders. High expectations are put on personalized genomics by generating WGS data for large numbers of individual horses, as well as on the ongoing equine initiative of the Functional Annotation of Animal Genomes (FAANG) project (reviewed by [ 55 ]) to identify important functional elements in the horse genome, particularly those underlying complex traits such as development and reproduction.","source_license":"CC-BY-4.0","license_restricted":false}