Comparative Genomics Approach to Identify Candidate Genetic Loci for Male Fertility

Towards a Multi-Omics of Male Infertility

Infertility is a common problem affecting one in six couples and in 30% of infertile couples, the male factor is a major cause. A large number of genes are involved in spermatogenesis and a significant proportion of male infertility phenotypes are of genetic origin. Studies on infertility have so far primarily focused on chromosomal abnormalities and sequence variants in protein-coding genes and have identified a large number of disease-associated genes. However, it has been shown that a multitude of factors across various omics levels also contribute to infertility phenotypes. The complexity of male infertility has led to the understanding that an integrated, multi-omics analysis may be optimal for unravelling this disease. While there is a vast array of different factors across omics levels associated with infertility, the present review focuses on known factors from the genomics, epigenomics, transcriptomics, proteomics, metabolomics, glycomics, lipidomics, miRNomics, and integrated omics levels. These include: repeat expansions in AR, POLG, ATXN1, DMPK, and SHBG, multiple SNPs, copy number variants in the AZF region, disregulated miRNAs, altered H3K9 methylation, differential MTHFR, MEG3, PEG1, and LIT1 methylation, altered protamine ratios and protein hypo/hyperphosphorylation. This integrative review presents a step towards a multi-omics approach to understanding the complex etiology of male infertility. Currently only a few genetic factors, namely chromosomal abnormalities and Y chromosome microdeletions, are routinely tested in infertile men undergoing intracytoplasmic sperm injection. A multi-omics approach to understanding infertility phenotypes may yield a more holistic view of the disease and contribute to the development of improved screening methods and treatment options. Therefore, beside discovering as of yet unknown genetic causes of infertility, integrating multiple fields of study could yield valuable contributions to the understanding of disease development. Future multi-omics studies will enable to synthesise fragmented information and facilitate biomarker discovery and treatments in male infertility.

Top Trends in Multiomics Research: Evaluation of 52 Published Studies and New Ways of Thinking Terminology and Visual Displays

Multiomics study designs have significantly increased understanding of complex biological systems. The multiomics literature is rapidly expanding and so is their heterogeneity. However, the intricacy and fragmentation of omics data are impeding further research. To examine current trends in multiomics field, we reviewed 52 articles from PubMed and Web of Science, which used an integrated omics approach, published between March 2006 and January 2021. From studies, data regarding investigated loci, species, omics type, and phenotype were extracted, curated, and streamlined according to standardized terminology, and summarized in a previously developed graphical summary. Evaluated studies included 21 omics types or applications of omics technology such as genomics, transcriptomics, metabolomics, epigenomics, environmental omics, and pharmacogenomics, species of various phyla including human, mouse, Arabidopsis thaliana, Saccharomyces cerevisiae, and various phenotypes, including cancer and COVID-19. In the analyzed studies, diverse methods, protocols, results, and terminology were used and accordingly, assessment of the studies was challenging. Adoption of standardized multiomics data presentation in the future will further buttress standardization of terminology and reporting of results in systems science. This shall catalyze, we suggest, innovation in both science communication and laboratory medicine by making available scientific knowledge that is easier to grasp, share, and harness toward medical breakthroughs.

Two mouse lines selected for large litter size display different lifetime fecundities

We recently described two outbred mouse lines that were selected for large litter size at first delivery. However, lifetime fecundity appears to be economically more important for the husbandry of many polytocous species for which mouse lines might serve as bona fide animal models (e.g. for pigs). In the present study, we compared the lifetime fecundities of two highly fertile mouse lines (FL1 and FL2: >20 offspring/litter at first delivery) with those of an unselected control line (ctrl) and two lines that were selected for high body weight (DU6) and high protein mass (DU6P) without selection pressure on fertility. We tested the hypothesis that selection for large litter size at first parturition would also increase lifetime fecundity in mice, and we observed very large differences between lines. Whereas FL1 and ctrl delivered up to nine and ten litters, none of the DU6 and DU6P females gave birth to more than five litters. In line with this observation, FL1 delivered the most pups per lifetime (85.7/female). FL2 females produced the largest average litter sizes (20.4 pups/litter) in the first four litters; however, they displayed a reduced number of litters. With the exception of ctrl, litter sizes declined from litter to litter. Repeated delivery of litters with high offspring numbers did not affect the general health of FL females. The presented data demonstrate that two biodiverse, highly fertile mouse lines selected for large litter size at first delivery show different lifetime reproductive fitness levels. Thus, these mouse lines might serve as valuable mouse models for investigating lifetime productivity and longevity in farm animals.

Association of α/β-Hydrolase D16B with Bovine Conception Rate and Sperm Plasma Membrane Lipid Composition

We have identified a Holstein sire named Tarantino who had been approved for artificial insemination that is based on normal semen characteristics (i.e., morphology, thermoresistance, motility, sperm concentration), but had no progeny after 412 first inseminations, resulting in a non-return rate (NRdev) of -29. Using whole genome association analysis and next generation sequencing, an associated nonsense variant in the α/β-hydrolase domain-containing 16B gene (ABHD16B) on bovine chromosome 13 was identified. The frequency of the mutant allele in the German Holstein population was determined to be 0.0018 in 222,645 investigated cattle specimens. The mutant allele was traced back to Whirlhill Kingpin (bornFeb. 13th, 1959) as potential founder. The expression of ABHD16B was detected by Western blotting and immunohistochemistry in testis and epididymis of control bulls. A lipidome comparison of the plasma membrane of fresh semen from carriers and controls showed significant differences in the concentration of phosphatidylcholine (PC), diacylglycerol (DAG), ceramide (Cer), sphingomyelin (SM), and phosphatidylcholine (-ether) (PC O-), indicating that ABHD16B plays a role in lipid biosynthesis. The altered lipid contents may explain the reduced fertilization ability of mutated sperms.

Reproductive performance primarily depends on the female genotype in a two-factorial breeding experiment using high-fertility mouse lines

Mouse models showing an improved fertility phenotype are barely described in the literature. In the present study, we further characterized two outbred mouse models that have been selected for the phenotype 'high fertility' for more than 177 generations (fertility lines (FL) 1 and 2). In order to delineate the impact of males and females on fertility parameters, we performed a two-factorial breeding experiment by mating males and females of the three different genotypes (FL1, FL2, unselected control (Ctrl)) in all 9 possible combinations. Reproductive performance, such as number of offspring per litter or total birth weight of the entire pup, mainly depends on the female genotype. Although the reproductive performance of FL1 and FL2 is very similar, their phenotypes differ. FL2 animals of both genders are larger compared to FL1 and control animals. Females of the control line delivered offspring earlier compared to FL1 and FL2 dams. Males of FL1 are the lightest and the only ones who gained weight during the two weeks mating period. To address whether this effect is correlated with differing serum androgen levels, we measured the concentrations of testosterone, dehydroepiandrosterone, 4-androstenedione, androstanediol and dihydrotestosterone in males of all three lines by GC-MS. We measured serum testosterone between 5.0 and 6.4 ng/mL, whereas the concentrations of the other androgens were at least one order of magnitude lower, with no significant differences between the lines. Our data indicate that reproductive outcome largely depends on the genotype of the female in a two-factorial breeding experiment and supports previous findings that the phenotype 'high fertility' is warranted by using different physiological strategies.

Initiative for standardization of reporting genetics of male infertility

The number of publications on research of male infertility is increasing. Technologies used in research of male infertility generate complex results and various types of data that need to be appropriately managed, arranged, and made available to other researchers for further use. In our previous study, we collected over 800 candidate loci for male fertility in seven mammalian species. However, the continuation of the work towards a comprehensive database of candidate genes associated with different types of idiopathic human male infertility is challenging due to fragmented information, obtained from a variety of technologies and various omics approaches. Results are published in different forms and usually need to be excavated from the text, which hinders the gathering of information. Standardized reporting of genetic anomalies as well as causative and risk factors of male infertility therefore presents an important issue. The aim of the study was to collect examples of diverse genomic loci published in association with human male infertility and to propose a standardized format for reporting genetic causes of male infertility. From the currently available data we have selected 75 studies reporting 186 representative genomic loci which have been proposed as genetic risk factors for male infertility. Based on collected and formatted data, we suggested a first step towards unification of reporting the genetics of male infertility in original and review studies. The proposed initiative consists of five relevant data types: 1) genetic locus, 2) race/ethnicity, number of participants (infertile/controls), 3) methodology, 4) phenotype (clinical data, disease ontology, and disease comorbidity), and 5) reference. The proposed form for standardized reporting presents a baseline for further optimization with additional genetic and clinical information. This data standardization initiative will enable faster multi-omics data integration, database development and sharing, establishing more targeted hypotheses, and facilitating biomarker discovery.

Sperm superoxide dismutase is associated with bull fertility

Decreasing mammalian fertility and sperm quality have created an urgent need to find effective methods to distinguish non-viable from viable fertilising spermatozoa. The aims of the present study were to evaluate expression levels of ?-tubulin 2C (TUBB2C), heat shock protein 10 (HSP10), hexokinase 1 (HXK1) and superoxide dismutase 1 (SOD1) in spermatozoa from Holstein bulls with varying fertility using western blotting and to analyse the biological networks of these key sperm proteins using a bioinformatics software (Metacore; Thomson-Reuters, Philadelphia, PA, USA). The rationales behind this study were that the sperm proteins play crucial roles in fertilisation and early embryonic development in mammals and ascertaining the biological networks of the proteins helps us better understand sperm physiology and early mammalian development. The results showed that expression of SOD1 was higher in spermatozoa from high fertility bulls (PPin vivo bull fertility. The findings are important because they illuminate molecular and cellular determinants of sperm viability and the identified protein markers can be used to determine bull fertility.

High-fertility phenotypes: two outbred mouse models exhibit substantially different molecular and physiological strategies warranting improved fertility

Animal models are valuable tools in fertility research. Worldwide, there are more than 400 transgenic or knockout mouse models available showing a reproductive phenotype; almost all of them exhibit an infertile or at least subfertile phenotype. By contrast, animal models revealing an improved fertility phenotype are barely described. This article summarizes data on two outbred mouse models exhibiting a 'high-fertility' phenotype. These mouse lines were generated via selection over a time period of more than 40 years and 161 generations. During this selection period, the number of offspring per litter and the total birth weight of the entire litter nearly doubled. Concomitantly with the increased fertility phenotype, several endocrine parameters (e.g. serum testosterone concentrations in male animals), physiological parameters (e.g. body weight, accelerated puberty, and life expectancy), and behavioral parameters (e.g. behavior in an open field and endurance fitness on a treadmill) were altered. We demonstrate that the two independently bred high-fertility mouse lines warranted their improved fertility phenotype using different molecular and physiological strategies. The fertility lines display female- as well as male-specific characteristics. These genetically heterogeneous mouse models provide new insights into molecular and cellular mechanisms that enhance fertility. In view of decreasing fertility in men, these models will therefore be a precious information source for human reproductive medicine. Translated abstract A German translation of abstract is freely available at http://www.reproduction-online.org/content/147/4/427/suppl/DC1.

Obesity gene atlas in mammals

Obesity in humans has increased at an alarming rate over the past two decades and has become one of the leading public health problems worldwide. Studies have revealed a large number of genes/markers that are associated with obesity and/or obesity-related phenotypes, indicating an urgent need to develop a central database for helping the community understand the genetic complexity of obesity. In the present study, we collected a total of 1,736 obesity associated loci and created a freely available obesity database, including 1,515 protein-coding genes and 221 microRNAs (miRNAs) collected from four mammalian species: human, cattle, rat, and mouse. These loci were integrated as orthologs on comparative genomic views in human, cattle, and mouse. The database and genomic views are freely available online at: http://www.integratomics-time.com/fat_deposition. Bioinformatics analyses of the collected data revealed some potential novel obesity related molecular markers which represent focal points for testing more targeted hypotheses and designing experiments for further studies. We believe that this centralized database on obesity and adipogenesis will facilitate development of comparative systems biology approaches to address this important health issue in human and their potential applications in animals.

Pivotal role of the muscle-contraction pathway in cryptorchidism and evidence for genomic connections with cardiomyopathy pathways in RASopathies

BACKGROUND

Cryptorchidism is the most frequent congenital disorder in male children; however the genetic causes of cryptorchidism remain poorly investigated. Comparative integratomics combined with systems biology approach was employed to elucidate genetic factors and molecular pathways underlying testis descent.

METHODS

Literature mining was performed to collect genomic loci associated with cryptorchidism in seven mammalian species. Information regarding the collected candidate genes was stored in MySQL relational database. Genomic view of the loci was presented using Flash GViewer web tool (http://gmod.org/wiki/Flashgviewer/). DAVID Bioinformatics Resources 6.7 was used for pathway enrichment analysis. Cytoscape plug-in PiNGO 1.11 was employed for protein-network-based prediction of novel candidate genes. Relevant protein-protein interactions were confirmed and visualized using the STRING database (version 9.0).

RESULTS

The developed cryptorchidism gene atlas includes 217 candidate loci (genes, regions involved in chromosomal mutations, and copy number variations) identified at the genomic, transcriptomic, and proteomic level. Human orthologs of the collected candidate loci were presented using a genomic map viewer. The cryptorchidism gene atlas is freely available online: http://www.integratomics-time.com/cryptorchidism/. Pathway analysis suggested the presence of twelve enriched pathways associated with the list of 179 literature-derived candidate genes. Additionally, a list of 43 network-predicted novel candidate genes was significantly associated with four enriched pathways. Joint pathway analysis of the collected and predicted candidate genes revealed the pivotal importance of the muscle-contraction pathway in cryptorchidism and evidence for genomic associations with cardiomyopathy pathways in RASopathies.

CONCLUSIONS

The developed gene atlas represents an important resource for the scientific community researching genetics of cryptorchidism. The collected data will further facilitate development of novel genetic markers and could be of interest for functional studies in animals and human. The proposed network-based systems biology approach elucidates molecular mechanisms underlying co-presence of cryptorchidism and cardiomyopathy in RASopathies. Such approach could also aid in molecular explanation of co-presence of diverse and apparently unrelated clinical manifestations in other syndromes.

Genome-wide association study implicates testis-sperm specific FKBP6 as a susceptibility locus for impaired acrosome reaction in stallions

Impaired acrosomal reaction (IAR) of sperm causes male subfertility in humans and animals. Despite compelling evidence about the genetic control over acrosome biogenesis and function, the genomics of IAR is as yet poorly understood, providing no molecular tools for diagnostics. Here we conducted Equine SNP50 Beadchip genotyping and GWAS using 7 IAR–affected and 37 control Thoroughbred stallions. A significant (P<6.75E-08) genotype–phenotype association was found in horse chromosome 13 in FK506 binding protein 6 (FKBP6). The gene belongs to the immunophilins FKBP family known to be involved in meiosis, calcium homeostasis, clathrin-coated vesicles, and membrane fusions. Direct sequencing of FKBP6 exons in cases and controls identified SNPs g.11040315G>A and g.11040379C>A (p.166H>N) in exon 4 that were significantly associated with the IAR phenotype both in the GWAS cohort (n = 44) and in a large multi-breed cohort of 265 horses. All IAR stallions were homozygous for the A-alleles, while this genotype was found only in 2% of controls. The equine FKBP6 was exclusively expressed in testis and sperm and had 5 different transcripts, of which 4 were novel. The expression of this gene in AC/AG heterozygous controls was monoallelic, and we observed a tendency for FKBP6 up-regulation in IAR stallions compared to controls. Because exon 4 SNPs had no effect on the protein structure, it is likely that FKBP6 relates to the IAR phenotype via regulatory or modifying functions. In conclusion, FKBP6 was considered a susceptibility gene of incomplete penetrance for IAR in stallions and a candidate gene for male subfertility in mammals. FKBP6 genotyping is recommended for the detection of IAR–susceptible individuals among potential breeding stallions. Successful use of sperm as a source of DNA and RNA propagates non-invasive sample procurement for fertility genomics in animals and humans.

Impaired acrosomal reaction (IAR) of sperm causes male subfertility in humans and animals, and currently the molecular causes of the condition are not known. Here we report the mapping, identification, and functional analysis of a susceptibility locus for IAR in stallions. The candidate region was mapped to horse chromosome 13 by SNP genotyping and GWAS of 7 IAR affected and 44 control Thoroughbred stallions. Re-sequencing and case-control analysis of functionally relevant candidate genes in the region identified FKBP6 gene as a significantly associated locus. The association was confirmed by genotyping 265 male horses of multiple breeds. FKBP6 belongs to the immunophilins FKBP family known to be involved in meiosis, calcium homeostasis, clathrin-coated vesicles, and membrane fusions. We showed that the equine FKBP6 is exclusively and monoallelically expressed in testis and sperm and has 5 different transcripts, of which 4 were novel. Overall, FKBP6 was considered a susceptibility gene of incomplete penetrance for IAR in stallions and a candidate gene for male subfertility in other mammals. Successful use of sperm as a source of DNA and RNA propagates non-invasive sample procurement for fertility genomics in animals and humans.

Chemocentric informatics approach to drug discovery: identification and experimental validation of selective estrogen receptor modulators as ligands of 5-hydroxytryptamine-6 receptors and as potential cognition enhancers

We have devised a chemocentric informatics methodology for drug discovery integrating independent approaches to mining biomolecular databases. As a proof of concept, we have searched for novel putative cognition enhancers. First, we generated Quantitative Structure-Activity Relationship (QSAR) models of compounds binding to 5-hydroxytryptamine-6 receptor (5-HT(6)R), a known target for cognition enhancers, and employed these models for virtual screening to identify putative 5-HT(6)R actives. Second, we queried chemogenomics data from the Connectivity Map ( http://www.broad.mit.edu/cmap/ ) with the gene expression profile signatures of Alzheimer's disease patients to identify compounds putatively linked to the disease. Thirteen common hits were tested in 5-HT(6)R radioligand binding assays and ten were confirmed as actives. Four of them were known selective estrogen receptor modulators that were never reported as 5-HT(6)R ligands. Furthermore, nine of the confirmed actives were reported elsewhere to have memory-enhancing effects. The approaches discussed herein can be used broadly to identify novel drug-target-disease associations.