Introduction
Male infertility is a complex, often multifactorial pathological condition affecting ~7% of the global male population (1). While the etiology of infertility is wide-ranging, up to 15% of males with infertility have an underlying causative genetic defect (1–3). Mouse knockout models have identified >400 genes, including ATM, SYCP1-3, SYCE1 and CADM, leading to monogenic infertility in men (4). Studies on the molecular mechanisms of cell-cell interactions have provided a better understanding of the causes of fertility issues due to fertilization defects (5,6). While there are established guidelines and practices for diagnostic karyotyping, azoospermia factor (AZF) deletion screening and cystic fibrosis transmembrane conductance regulator testing, the overall diagnostic yield is as low as 4% (7). However, rapid advancements in genomic medicine achieved by next-generation sequencing (NGS) technologies lack of results standardization, creating a knowledge gap in the clinical area of male infertility (8).
A literature review on monogenic forms of male infertility identified 78 genes associated with male infertility in 1,337 publications in 2019 (9), with this number rising by 33% to 104 genes by 2022 (10). In the era of genomics and personalized medicine, this number is low compared with other clinical fields, such as intellectual disability (8). To improve biological understanding as well as diagnostic yield and clinical relevance of genetic testing, further genes leading to male infertility must be identified.
Our previous study outlined a possible cause for asthenoteratozoospermia in a patient harboring a deletion of ~8-Mb in the 5q22.2q23.1 locus, including the testis-specific serine kinase 1B (TSSK1B) gene (11). Genes belonging to kinase and phosphatase families are responsible for activation and deactivation of intra- and extracellular transduction via phosphorylation and dephosphorylation (12). These processes are key for correct regulation and metabolic processes of the cell that are mediated by different receptors and enzymes (13). TSSKs are part of the AMP-activated protein kinase family (14). They comprise six genes, which are almost exclusively present in the testes (>1,000-fold concentration compared with other organs). These genes serve a role in spermatogenesis and are responsible for correct morphogenesis and differentiation following meiosis when spermatid elongation occurs. This has been demonstrated using recombinant mouse models, where sterile phenotypes have been observed in double Tssk1 and Tssk2 knockout (KO) and Tssk6 KO mice (15). A sub-fertile phenotype with reduced TSSK1 and TSSK2 levels has also been observed in a Tssk4 KO model (16,17). These results highlight the importance of TSSKs in research, as they can be targeted by both drugs and inhibitors, suggesting that TSSKs may not only cause infertility, but can also provide a route for the development of male non-hormonal contraceptive drugs. Originally, TSSK genes were found in mouse tissue using degenerate oligonucleotide primers while searching for novel kinases, which led to the discovery of TSSK1. Subsequently, TSSK2 was identified via low stringency screening due to its close proximity and genetic linkage to TSSK1 (18,19). Following the discovery of these genes, yeast two-hybrid and co-immunoprecipitation experiments were performed to detect proteins that may interact with these kinases. A novel 65 kDa protein was identified, namely testis-specific kinase substrate (TSKS), that interacts with both TSSK1 and TSSK2 (19).
According to studies conducted on a diverse array of species, it was proposed that the TSSK family originated with Tssk5 in amphibians ~380 million years ago (MYA). Following the Paleocene-Eocene radiation period, a novel gene appeared in primates and humans, known as Tssk1b (Fig. 1). Moreover, in most species Tssk1/2 are linked on one chromosome and their activity is hypothesized to be parallel, as their combined absence has been shown to cause infertility in KO mice (17). However, with the appearance of Tssk1b (duplicated on another chromosome and not linked to Tssk2), Tssk1 was inactivated through negative selection and mutation to become a pseudogene known as Tssk1 (Fig. 1). Even though the genes diverge from one another, they exhibit sequence conservation and similarity. The N-terminal domain is represented by the 1–272 amino acid sequence, which bears the kinase domain, whereas the C-terminal is present in the 273-end sequence. In humans, N-terminal similarity between Tssk1b and Tssk2 is 82.0%, which is similar to other mammalian species, while C-terminal similarity between the two genes is notably lower at 14.7%. This may be due to the regulatory function of the C-terminal domain, which is specific for each kinase. This could also mean that each kinase serves specific functions, in addition to exhibiting overlapping effects. With regards to evolutionary conservation, the C-terminal of Tssk1 shows 65.7% mean sequence conservation compared with the C-terminal of Tssk2, which is 87.6%. When considering retrogenes and duplications, newer genes seem to have the ability to evolve faster, meaning that their conservation rate will be lower. Based on these observations on conservation rates, it is hypothesized that Tssk2 was the first gene to appear on the chromosome following the retroposition of Tssk1 (20). Compared with non-mammalian species, the C-terminal conservation of human Tssk1b is lower (61.4%) than that of Tssk2 (86.8%), whereas N-terminal conservation is similar (92.5 and 94.9%, respectively) (20). This suggests that the two domains are affected by different selective evolutionary pressures.
Spermatogenesis is divided into three main phases: Mitotic division, which generates a pool of spermatocytes; meiosis to generate haploid spermatids and spermiogenesis, in which spermatids differentiate into spermatozoa (21). Nayak et al (22) tested mouse tissue obtained from maturing seminiferous tubules at different periods after birth. The procedure included immuno-fluorescent staining to visualize difference in expression of both Tssk1 and Tssk2. Tssk1 expression was first observed at low levels in week 3, concurrent with the appearance of round spermatids, predominantly in cells in meiotic metaphase. By contrast, Tssk2 was detected in cells undergoing spermiogenesis at weeks 4–5 and was not detected in metaphase dividing cells (22). This highlights the importance of TSSK1 (the first gene to be expressed) in both early sperm development and mature sperm formation (22).
In spermatogenesis, following meiosis, a cloud-like organelle (nuage) appears in the round spermatid known as chromatoid body and is localized around the haploid nucleus (23). The chromatoid body is the RNA-controlling center responsible for organization and regulation of mRNA and RNA pathways associated with the haploid genome of the spermatid (Fig. 2A) (24). Its function is rather short-term and during the transition from round to elongated spermatid, the chromatoid body loses certain functions and enzymes (such as murine P-element induced wimpy testis family member proteins) to split into two well-defined structures, namely satellite and ring (Fig. 2B) (25). It has been observed that Tssk1 and Tssk2, as well as TSKS, accumulate in both the ring around the flagellum and in the satellite in the cytoplasm. Experiments in both wild-type and double Tssk1,2 КО mice revealed that КО results in disturbed mitochondrial sheath formation, rendering the sperm cells non-functional (17).
The sheath-forming complex includes not only TSSK1B, but also testis-specific phosphatase, Ppp1cc2, which serves a key role in completion of spermatogenesis in mice (26). Pull-down assays using GST-Ppp1cc2 expressed and purified from bacteria as bait against protein lysates from mouse testis tissue have been performed to determine the interaction between Ppp1cc2 and other proteins (27). An indirect interaction of Ppp1cc2 with Tssk1 was shown to be mediated through TSKS via RVxF motif (at amino acid position 51–55), which interacts with both proteins to form a complex. To test the activity of Ppp1cc, deletion via mutagenesis was performed in male mice; this resulted in germ cell reduction primarily occurring during spermatid elongation (post-meiosis) phase at the time of chromatoid body dissociation, causing the spermatogenic cycle to halt (28). Thus, the testicular kinase/phosphatase complex is key for formation of the mitochondrial sheath during spermatogenesis.
Considering the literature and our previous experience regarding the role of the TSSK1B gene in a patient with asthenoteratozoospermia, the aim of the present study was to perform TSSK1B genetic screening on a selected group of mutation-negative male patients with infertility. The study aimed to assess overall genomic variability within TSSK1B compared with previously reported variants in population databases, investigate the association between missense mutations in TSSK1B and clinical phenotype with regards to semen quality and to assess findings using in silico predictors and novel protein folding algorithms.
Discussion
Our previous report presented the first known case of a TSSK1B deletion associated with a clinical phenotype of asthenoteratozoospermia and male infertility (11). Since then, the TSSK family has emerged as a potential cause for male infertility and may be a target for the development of novel male contraceptive solutions through targeted inhibition of TSSK1B (15). By performing comprehensive targeted patient screening to investigate the role of TSSK1B in male infertility, the present study discovered both known and novel missense variants that should be evaluated in vivo to determine the extent of their clinical manifestation. Mutations discovered toward the C-terminus of TSSK1B, namely p.3D>N and p.52F>L, did not exhibit the same phenotype as in our initial report (asthenoteratozoospermia), yet yielded pathogenic scores via in silico prediction methods. Of note, p.66M>V and p.237R>C mutations were associated with the clinical phenotype discovered in our initial study. Finally, the p.293G>E mutation manifested as two clinical phenotypes, azoospermia and teratozoospermia. While classification via ACMG guidelines shows either that mutations are of unknown significance or possibly benign, this may be due to >0.01 population frequency in reference databases (34). Infertility-causing mutations might not be functionally investigated because of their high representation in a population. Such sequence alterations may not directly affect the quality of life of the individual; hence they are treated as benign.
While in silico predictors such as PolyPhen and SIFT (31) have been used to determine the potential outcome of a missense mutation, the next step in the assessment of novel pathological variants is to build and compare protein models. Although AlphaFold2 does not provide a solution to the protein folding problem (34), it has an accuracy level comparable to that of X-ray crystallography (35). This novel artificial intelligence program achieved a milestone level of accuracy in the biannual Critical Assessment of Protein Structure Prediction experiment in 2020 and is used to predict protein structures for databases such as the European Bioinformatics Institute. While it is still not validated for clinical use, the results produced are notable, with a margin of error of only 1.6 angstroms when simulating the folding of the protein (36). Using this method, altered amino acid sequence caused by the missense mutations identified in our patient cohort were simulated; AlphaFold2 predicted a small helical motif in four of the altered sequences. This fold occurred seemingly with no regard to the location of the change yet it was predicted to affect the kinase domain of TSSK1B. As previous studies show, damaging or removing one of the copies of this gene can cause infertility (14,17), which indicates that small structure alterations lead to a pathological changes.
To support the present findings, the evolutionary background of TSSKs and integration of the Ppp1cc2-Tssk1 complex were investigated, which suggested that this gene is key for sperm maturity. TSSK1B was the first of the family to evolve in humans and the complex with Ppp1cc2 is essential for proceeding with the spermatid elongation stage and completing spermatogenesis. Furthermore, TSSK1B is the earliest kinase of this family to be expressed, at 3 weeks opposed to other TSSKs expressed at 3.5 weeks (22,37), which further outlines its importance in spermatogenesis.
In the present cohort, 11 of 100 patients carried a missense mutation in TSSK1B; although this frequency is biased by the cohort size and patient selection, an 11% carrier rate warrants further study. Not all mutations are equal: In silico predictors showed that 3 of the 5 protein-altering variants have scores that are interpreted as damaging. As the N-terminus of TSSK1 is highly conserved, this suggests a potential damaging outcome caused by the mutations reported in the present study. Further functional studies are required to verify the pathological consequences of these alterations.
Multiple studies have highlighted the need for novel male contraceptives (38,39). Current methods are limited to condoms and vasectomy, compared with the variety of available female contraceptives; the limitations of these methods are high failure rate of condoms and irreversibility of surgical vasectomy (39). While trials of hormonal methods are successful, they come with side effects, such as weight gain, acne, mood changes and changes in libido (39,40). Previous studies demonstrate promising results in targeting TSSKs using kinase inhibitors and the potential application of TSSK allosteric inhibitors (15,41), thereby emphasizing not only the important role of this family of kinases in spermatogenesis, but also how they can be targeted pharmaceutically.
The only putative pathological manifestation of TSSK1B has been studied in mice, where haploinsufficiency results in decreased sperm maturation and offspring carrying exclusively the wild-type TSSK1B allele (14). The present patient cohort was referred for a number of clinical phenotypes of infertility and some individuals carried a heterozygous variant of TSSK1B. Functional study investigating TSSK1B function in human spermatozoids is warranted. Second, through the advent of NGS, carrier screening can be performed in cases with previously unidentified genetic causes. Lastly, research into this gene may lead to novel male contraceptive solutions.