Implication of FKBP6 for Male Fertility in Horses

Transcriptomic analysis of testis and epididymis tissues from Banna mini-pig inbred line boars with single-molecule long-read sequencing†

In mammals, testis and epididymis are critical components of the male reproductive system for androgen production, spermatogenesis, sperm transportation, as well as sperm maturation. Here, we report single-molecule real-time sequencing data from the testis and epididymis of the Banna mini-pig inbred line (BMI), a promising laboratory animal for medical research. We obtained high-quality full-length transcriptomes and identified 9879 isoforms and 8761 isoforms in the BMI testis and epididymis, respectively. Most of the isoforms we identified have novel exon structures that will greatly improve the annotation of testis- and epididymis-expressed genes in pigs. We also found that 3055 genes (over 50%) were shared between BMI testis and epididymis, indicating widespread expression profiles of genes related to reproduction. We characterized extensive alternative splicing events in BMI testis and epididymis and showed that 96 testis-expressed genes and 79 epididymis-expressed genes have more than six isoforms, revealing the complexity of alternative splicing. We accurately defined the transcribed isoforms in BMI testis and epididymis by combining Pacific Biotechnology Isoform-sequencing (PacBio Iso-Seq) and Illumina RNA Sequencing (RNA-seq) techniques. The refined annotation of some key genes governing male reproduction will facilitate further understanding of the molecular mechanisms underlying BMI male sterility. In addition, the high-confident identification of 548 and 669 long noncoding RNAs (lncRNAs) in these two tissues has established a candidate gene set for future functional investigations. Overall, our study provides new insights into the role of the testis and epididymis during BMI reproduction, paving the path for further studies on BMI male infertility.

Genome-wide association study for frozen-thawed sperm motility in stallions across various horse breeds

OBJECTIVE

The semen quality of stallions including sperm motility is an important target of selection as it has a high level of individual variability. However, effects of the molecular architecture of the genome on the mechanisms of sperm formation and their preservation after thawing have been poorly investigated. Here, we conducted a genome-wide association study (GWAS) for the sperm motility of cryopreserved semen in stallions of various breeds.

METHODS

Semen samples were collected from the stallions of 23 horse breeds. The following semen characteristics were examined: progressive motility (PM), progressive motility after freezing (FPM), and the difference between PM and FPM. The respective DNA samples from these stallions were genotyped using Axiom Equine Genotyping Array.

RESULTS

We performed a GWAS search for single nucleotide polymorphism (SNP) markers and potential genes related to motility properties of frozen-thawed semen in the stallions of various breeds. As a result of the GWAS analysis, two SNP markers, rs1141327473 and rs1149048772, were identified that were associated with preservation of the frozen-thawed stallion sperm motility, the relevant putative candidate genes being NME/NM23 family member 8 (NME8), olfactory receptor family 2 subfamily AP member 1 (OR2AP1), and olfactory receptor family 6 subfamily C member 4 (OR6C4). Potential implications of effects of these genes on sperm motility are herein discussed.

CONCLUSION

The GWAS results enabled us to localize novel SNPs and candidate genes for sperm motility in stallions. Implications of the study for horse breeding and genetics are a better understanding of genomic regions and candidate genes underlying stallion sperm quality, and improvement in horse reproduction and breeding techniques. The identified markers and genes for sperm cryotolerance and the respective genomic regions are promising candidates for further studying the biological processes in the formation and function of the stallion reproductive system.

Thoroughbred stallion fertility is significantly associated with FKBP6 genotype but not with inbreeding or the contribution of a leading sire

This is a follow-up study to validate the previously detected association of the FKBP6 gene with stallion subfertility. Using a select cohort of 150 Thoroughbred stallions with detailed breeding records, we confirm significant association (P < 0.0001) between low per-cycle pregnancy rates (≤50%) and a combined A/A-A/A genotype of SNPs chr13:11 353 372G>A and chr13:11 353 436A>C in FKBP6 exon 5. We also show that stallion subfertility and the combined genotype A/A-A/A are not associated with the level of genetic diversity based on 12 autosomal microsatellite markers, or with pedigree-based inbreeding rate, or the extent of contribution of a leading Thoroughbred sire, Northern Dancer, in a stallion's pedigree. We develop a TaqMan allelic discrimination assay for the two SNPs to facilitate accurate and high-throughput genotyping. We determine allele, genotype and combined genotype frequencies of FKBP6 exon 5 SNPs in a global cohort of 518 Thoroughbreds (76% stallions or geldings and 24% mares) and show that the frequency of the A/A-A/A genotype is 4%. Because there is no similar association between the FKBP6 exon 5 genotype and stallion subfertility in Hanoverians, we suggest that the two SNPs are not causative but rather tagging a breed-specific haplotype with genetic variants unique to Thoroughbreds. Further WGS-based research is needed to identify the molecular causes underlying the observed genotype-phenotype association in Thoroughbred stallions.

Impaired Reproductive Function in Equines: From Genetics to Genomics

Fertility is one of the key factors in the economic and productive success of the equine industry. Despite this, studies on the genetic causes affecting reproductive performance are scarce, especially in mares, where the genetic architecture of the reproductive traits is extremely complex. Today, with the increasing availability of new genomic methodologies for this species, we are presented with an interesting opportunity to understand the genetic basis of equine reproductive disorders. These include, among others, novel techniques for detecting chromosomal abnormalities, whose association with infertility in horses was established over 50 years ago; new sequencing technologies permitting an accurate detection of point mutations influencing fertility, as well as the study of inbreeding and molecular homozygosity, which has been widely suggested as one of the main causes of low reproductive performance in horses. Finally, over the last few years, reproductive performance has also been associated with copy number variants and candidate genes detected by genome-wide association studies on fertility traits. However, such studies are still scarce, probably because they depend on the existence of large and accurate phenotypic datasets of reproductive and/or fertility traits, which are still difficult to obtain in equines.

Genetics of Equine Reproductive Diseases

Reproductive disorders are genetically heterogeneous and complex; available genetic tests are limited to chromosome analysis and 1 susceptibility gene. Cytogenetic analysis should be the first test to confirm or rule out chromosomal aberrations. No causative genes/mutations are known. The only available genetic test for stallion subfertility is based on a susceptibility gene FKBP6. The ongoing progress in equine genomics will improve the status of genetic testing. However, because subfertile phenotypes do not facilitate collection of large numbers of samples or pedigrees, and clinical causes of many cases remain unknown, further progress requires constructive cross-talk between geneticists, clinicians, breeders, and owners.

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.

Screening of whole genome sequences identified high-impact variants for stallion fertility

Background

Stallion fertility is an economically important trait due to the increase of artificial insemination in horses. The availability of whole genome sequence data facilitates identification of rare high-impact variants contributing to stallion fertility. The aim of our study was to genotype rare high-impact variants retrieved from next-generation sequencing (NGS)-data of 11 horses in order to unravel harmful genetic variants in large samples of stallions.

Methods

Gene ontology (GO) terms and search results from public databases were used to obtain a comprehensive list of human und mice genes predicted to participate in the regulation of male reproduction. The corresponding equine orthologous genes were searched in whole genome sequence data of seven stallions and four mares and filtered for high-impact genetic variants using SnpEFF, SIFT and Polyphen 2 software. All genetic variants with the missing homozygous mutant genotype were genotyped on 337 fertile stallions of 19 breeds using KASP genotyping assays or PCR-RFLP. Mixed linear model analysis was employed for an association analysis with de-regressed estimated breeding values of the paternal component of the pregnancy rate per estrus (EBV-PAT).

Results

We screened next generation sequenced data of whole genomes from 11 horses for equine genetic variants in 1194 human and mice genes involved in male fertility and linked through common gene ontology (GO) with male reproductive processes. Variants were filtered for high-impact on protein structure and validated through SIFT and Polyphen 2. Only those genetic variants were followed up when the homozygote mutant genotype was missing in the detection sample comprising 11 horses. After this filtering process, 17 single nucleotide polymorphism (SNPs) were left. These SNPs were genotyped in 337 fertile stallions of 19 breeds using KASP genotyping assays or PCR-RFLP. An association analysis in 216 Hanoverian stallions revealed a significant association of the splice-site disruption variant g.37455302G>A in NOTCH1 with the de-regressed estimated breeding values of the paternal component of the pregnancy rate per estrus (EBV-PAT). For 9 high-impact variants within the genes CFTR, OVGP1, FBXO43, TSSK6, PKD1, FOXP1, TCP11, SPATA31E1 and NOTCH1 (g.37453246G>C) absence of the homozygous mutant genotype in the validation sample of all 337 fertile stallions was obvious. Therefore, these variants were considered as potentially deleterious factors for stallion fertility.

Conclusions

In conclusion, this study revealed 17 genetic variants with a predicted high damaging effect on protein structure and missing homozygous mutant genotype. The g.37455302G>A NOTCH1 variant was identified as a significant stallion fertility locus in Hanoverian stallions and further 9 candidate fertility loci with missing homozygous mutant genotypes were validated in a panel including 19 horse breeds. To our knowledge this is the first study in horses using next generation sequencing data to uncover strong candidate factors for stallion fertility.

Electronic supplementary material

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

Fertility Assessment in Sorraia Stallions by Sperm-Fish and Fkbp6 Genotyping

The Sorraia, a critically endangered indigenous Iberian horse breed, is characterized by low genetic variability, high rate of inbreeding, bad sperm quality and subfertility. Here, we studied 11 phenotypically normal but subfertile Sorraia stallions by karyotyping, sex chromosome sperm-FISH and molecular analysis of FKBP6 - a susceptibility locus for impaired acrosome reaction (IAR). The stallions had normal sperm concentration (>300 million cells/ml), but the numbers of progressively motile sperm (21%) and morphologically normal sperm (28%) were invariably low. All stallions had a normal 64,XY karyotype. The majority of sperm (89%) had normal haploid sex chromosome content, although 11% of sperm carried various sex chromosome aneuploidies. No correlation was found between the percentage of sperm sex chromosome abnormalities and inbreeding, sperm morphology or stallion age. Direct sequencing of FKBP6 exon 4 for SNPs g.11040315G>A and g.11040379C>A revealed that none of the stallions had the susceptibility genotype (A/A-A/A) for IAR. Instead, all animals had a G/G-A/A genotype - a testimony of low genetic variability. The findings ruled out chromosomal abnormalities and genetic predisposition for IAR as contributing factors for subfertility. However, low fertility of the Sorraia stallions could be partly attributed to relatively higher rate of sex chromosome aneuploidies in the sperm.