Rapid evolution and molecular convergence in cryptorchidism-related genes associated with inherently undescended testes in mammals
BMC Ecology and Evolution volume 21, Article number: 22 (2021)
The mammalian testis is an important male exocrine gland and spermatozoa-producing organ that usually lies in extra-abdominal scrotums to provide a cooler environment for spermatogenesis and sperm storage. Testicles sometimes fail to descend, leading to cryptorchidism. However, certain groups of mammals possess inherently ascrotal testes (i.e. testes that do not descend completely or at all) that have the same physiological functions as completely descended scrotal testes. Although several anatomical and hormonal factors involved in testicular descent have been studied, there is still a paucity of comprehensive research on the genetic mechanisms underlying the evolution of testicular descent in mammals and how mammals with ascrotal testes maintain their reproductive health.
We performed integrative phenotypic and comparative genomic analyses of 380 cryptorchidism-related genes and found that the mammalian ascrotal testes trait is derived from an ancestral scrotal state. Rapidly evolving genes in ascrotal mammals were enriched in the Hedgehog pathway—which regulates Leydig cell differentiation and testosterone secretion—and muscle development. Moreover, some cryptorchidism-related genes in ascrotal mammals had undergone positive selection and contained specific mutations and indels. Genes harboring convergent/parallel amino acid substitutions between ascrotal mammals were enriched in GTPase functions.
Our results suggest that the scrotal testis is an ancestral state in mammals, and the ascrotal phenotype was derived multiple times in independent lineages. In addition, the adaptive evolution of genes involved in testicular descent and the development of the gubernaculum contributed to the evolution of ascrotal testes. Accurate DNA replication, the proper segregation of genetic material, and appropriate autophagy are the potential mechanisms for maintaining physiological normality during spermatogenesis in ascrotal mammals. Furthermore, the molecular convergence of GTPases is probably a mechanism in the ascrotal testes of different mammals. This study provides novel insights into the evolution of the testis and scrotum in mammals and contributes to a better understanding of the pathogenesis of cryptorchidism in humans.
The testis is a consequential male exocrine gland that produces spermatozoa and an endocrine gland that secretes sex hormones. For most mammals, both the testis and epididymis are located in the scrotum, which is outside of the body and protects the testis. The multi-version ‘cooling hypothesis’ [1, 2] suggests that the scrotum provides an environment 2–4 °C cooler than the normal body temperature; germinal epithelium and spermatozoa are acutely sensitive to heat , so this environment is optimal for spermatogenesis and sperm storage [4,5,6].
Some other hypotheses have been raised to explain the evolutionary origin of testicular descent and the function of the scrotum—e.g. the ‘training hypothesis’ argues that the scrotum exposes the sperm to a hostile environment to “train” it for further fertilization , Portmann  contends that the scrotum serves as a sexual signal in some mammals, and the ‘galloping hypothesis’ states that the scrotum originated in mammals that gallop or jump and protects spermatogenesis and sperm storage from consequent fluctuations in intra-abdominal pressure [9, 10].
Cryptorchidism (from Greek, meaning “hidden testicle”) is a failure of the testis to descend into the scrotal sac. For mammals with completely descended testes (CDT), cryptorchidism is a developmental defect that causes severe dysfunctions such as germ cell maldevelopment, asthenospermia, and imbalances in hormones such as testosterone (T) and anti-Müllerian hormone (AMH) . Furthermore, cryptorchidism is associated with an increased risk of testicular malignancy and other diseases [5, 11]. Overall, cryptorchidism, a congenital malformation in most mammals, can affect normal male physiological functions.
Anatomical and physiological adaptations of the scrotum—e.g. thin skin, no subcutaneous fat, absence of hair/fur, and cremaster muscle—keep the epididymis and testis cool to maintain male fertility and health . However, some mammals possess natural and properly functioning undescended testes (testes that do not descend completely or at all). For example, the platypus (Ornithorhynchus anatinus), a monotreme, has high intra-abdominal undescended testes (UDT) in the same initial position as the ovary in females and does not develop a scrotum . Most species in Afrotheria—e.g. elephant (Loxodonta africana), cape golden mole (Chrysochloris asiatica), and manatee (Trichechus manatus latirostris)—have UDT [14,15,16], whereas the aardvark (Orycteropus afer afer) has incompletely descended testes (IDT) and lacks a scrotum . In addition, armadillo (Dasypus novemcinctus; Xenarthra) and several lineages of Boreoeutheria (e.g. cetaceans, flying foxes, some pinnipeds, eulipotyphlans, and certain rodents) possess IDT .
In general, the CDT forms through a distinct and sequential two-phase descent , each phase of which involves multiple mechanical and hormonal factors. During the first phase, or the transabdominal phase, the testis is anchored to the inguinal region from a high abdominal position with the help of a swelling gubernaculum, the cranial suspensory ligament (CSL), T, and AMH [12, 19]. Early in the embryo development process, the gubernaculum primarily consists of a mesenchymal core and muscular outer layer, and has the potential to develop into a striated muscle bundle . Later on, the gubernaculum can be found as a striated muscle bundle that eventually connects to the abdominal wall and scrotum [21, 22]. The gubernacular swelling reaction is mainly controlled by the insulin-like 3 protein (INSL3) and its receptor leucine-rich repeat-containing G protein coupled receptor 8 (RXFP2) .
During the second—or the inguino-scrotal—phase, the testis migrates into the scrotum. This process is regulated by a combination of gubernaculum, testosterone, calcitonin gene-related peptide (CGRP) released by the genitofemoral nerve (GFN), and intra-abdominal pressure [12, 24, 25]. A recent study associated the absence of testicular descent to the inactivation of INSL3 and RXFP2 in several UDT mammalian lineages like tenrec (Echinops telfairi), cape golden mole, cape elephant shrew (Elephantulus edwardii), and manatee . However, the same study also showed that many other ascrotal mammals have full INSL3 and RXFP2 gene sequences, such as the elephant. This implies that the development and evolution of testicular descent in mammals likely involve polygene. Overall, both the etiology of cryptorchidism and the genetic mechanisms related to the evolution of testicular descent remain largely unknown.
This study analyzes cryptorchidism-related genes in mammals to investigate the mechanisms underlying developmental testicular descent and its adaptive importance in ascrotal mammals. This is important not only to understand the evolutionary mechanisms of testicular descent in mammals, but also to reveal the genetic mechanisms underlying cryptorchidism; thus, this study may be used to improve reproductive health in humans.
Scrotal CDT is the ancestral state and was lost independently in many lineages
We used a recent species tree  to map the phenotypic states of testicular position and the presence of a scrotum. The mapping shows that the testis and scrotum phenotypes are diverse across mammals (Fig. 1, Additional file 1: Table S1). A comparison of the evolutionary models found that ARD (All Rates Different), which treats all rates as different during trait evolution, was the best fitting model for testis and scrotum evolution in mammals (Additional file 1: Table S2). Reconstruction of the ancestral state suggested that scrotal CDT had a higher probability of appearing in deep ancestral nodes than did the undescended testis, e.g. the ancestor of mammals (node 1: 0.963), therians (node 2: 0.991), and placental mammals (node 5: 0.999) (Fig. 1). Moreover, the testis and scrotum phenotypes recovered several shifts from scrotal CDT to ascrotal UDT in internal nodes within Afrotheria, and to the ascrotal IDT in certain ancestral nodes of rodents, e.g. in the respective ancestors of cetaceans, megabats, and pinnipeds, plus the ancestral nodes within Eulipotyphla (Fig. 1).
Molecular evolution of cryptorchidism-related genes involved in testicular descent
Branch model analysis in PAML detected 36 genes that evolved significantly (adjusted p < 0.01) increased molecular substitution rates in ascrotal IDT and UDT branches compared to CDT species (Fig. 2 and Additional file 1: Table S3). Among them, substitution rates in ascrotal mammals were up to 10.7 times greater than they were in scrotal ones. In addition, 33 out of these 36 genes were still rapidly evolved in ascrotal species in an expanded 62-mammal data set (Additional file 2: Fig S1), confirming the validity of the result (Additional file 1: Table S4). A further functional and pathway enrichment analysis revealed that 36 rapidly evolving genes were significantly enriched in the Hedgehog signaling pathway of Kyoto Encyclopedia of Genes and Genomes (KEGG) (Fig. 2 and Additional file 1: Table S5) and terms related to reproductive development (reproductive structure development, reproductive system development, male gonad development, development of primary male sexual characteristics, and male sex differentiation), muscle (muscle tissue development, striated muscle tissue development, and muscle organ development), and sex hormone receptors (protein-hormone receptor activity and hormone binding) of Gene Ontology (GO) (adjusted p < 0.05) (Fig. 2 and Additional file 1: Table S6).
Next, we detected genes with increased substitution rates in UDT mammals, which evolved as the most unique and extreme case of testis-scrotum phenotype. Forty-six genes had rapidly evolved in UDT species (Additional file 2: Fig S2 and Additional file 1: Table S7).
Phylogenetic generalized least squares (PGLS) regression was applied to detect any potential relationship between the evolutionary rates of cryptorchidism-related genes and testicular descent. Thirteen genes were found to be significantly associated: MKX, TMEM74, BTBD1, SPPL2C, ATRX, NANOS1, AP3B2, IRF6, PRRG4, RAB3GAP2, DSCC1, EIF3A, and STS (Table 1). Among them, chromatin-remodeling protein ATRX is involved in transcriptional regulation and telomere replication . The gene NANOS1 (Nanos C2HC-Type Zinc Finger 1), a member of the nanos family, is associated with spermatogenic impairment and translational regulation . Further functional enrichment using GO annotations revealed terms related to sister chromatid cohesion (adjusted p < 0.05) (Additional file 1: Table S8 and Additional file 2: Fig S3).
We used a branch-site model to test positively selected genes and amino acids in ascrotal IDT and UDT mammals. Evidence of positive selection was found in 13 ascrotal mammal genes (FLNA, DOCK1, CSMD3, MCMBP, FANCE, FBXL18, ARMC4, DEPTOR, ACTA1, JAG1, AIMP2, RAF1, and CCDC73) (Table 2). Enrichment analysis indicated that the functions of these 13 genes were significantly related to muscle (actin filament, sarcomere, contractile fiber part, myofibril, and contractile fiber), DNA replication (MCM complex), and GTPase (guanyl-nucleotide exchange factor complex) after multiple testing (adjusted p < 0.05) (Additional file 2: Fig S4 and Additional file 1: Table S9).
UDT species-specific amino acid mutations and indels
We failed to detect any specific mutations in ascrotal mammals (UDT or IDT species), but we did successfully identify 17 specific substitutions in 14 proteins of UDT mammals (Fig. 3). The protein AMER1 had three specific substitutions, FAT1 had two, and the other 12 proteins had one each. In addition, two indels were detected mainly in UDT mammals (Fig. 4): (1) a 108-bp deletion in AXIN1 caused a 36-amino-acid fragment to be missing in all the UDT mammals and a small number of the IDT and CDT species and (2) a six-bp-long deletion in all UDT and some CDT species.
Molecular convergence in ascrotal mammals
We initially observed convergent/parallel amino acid substitutions in 164 proteins in 176 IDT/UDT ascrotal branches (Additional file 1: Table S10 and Additional file 2: Fig S5). The Poisson test detected three parallel substitutions in the protein CPEB1 among different UDT branches, and 18 parallel changes were found in 13 proteins—AHSA2, FHL3, GNRHR, PROSER2, RALBP1, SERTAD2, WDPCP, C15orf40, CDC42EP4, CDT1, CPEB1, PGLS, and WHAMM—among the IDT branches (Additional file 1: Table S11). A functional enrichment analysis found that these 13 genes were related to small GTPase functions, such as GTP-Rho binding, Rho GTPase binding, and GTPase activator activity (Additional file 1: Table S12 and Fig. 5).
Most mammals have completely descended testes (CDT). Although the failure of testicular descent leads to cryptorchidism, certain groups of mammals possess incompletely descended (IDT) or undescended testes (UDT) that function normally. The evolution of the scrotum and testicular descent in mammals has long been the subject of scientific interest in a wide range of fields, including medicine, developmental biology, and evolutionary biology. However, little is known about this topic from an evolutionary perspective. Here, comparative genomics and evolutionary analyses of cryptorchidism-related genes provide some novel evidence for the evolutionary trajectory of testicular descent and potential mechanisms driving normal physiological functions for IDT and UDT in mammals.
The ascrotal testis evolved multiple times independently in mammals
It has been suggested that descended testes and the scrotum provide an optimal cooling environment for spermatogenesis and sperm storage, since they are cooler in temperature than the core body [1, 2, 29]. However, propitious testicular descent into the scrotum could be a sophisticated process that is physiologically, developmentally, and evolutionarily costly . Mammals evolved divergent scrotal CDT as well as ascrotal IDT and UDT. However, unlike in hard tissues, it is difficult to trace the evolutionary history of these organs because they do not fossilize. Nevertheless, inferring ancestral character states by mapping various phenotypes of living taxa onto a phylogeny could suggest the evolutionary trajectory of a certain trait .
Our ancestral state reconstruction showed that several deep ancestors of mammals possessed scrotal CDT (Fig. 1). The derived characters of ascrotal IDT and UDT were inferred to occur in multiple lineages—e.g. Monotremata, Afrotheria, Cingulata, Rodentia, Chiroptera, Cetartiodactyla, Pinnipedia, and Eulipotyphla. Although Kleisner et al.  and Lovegrove  suggested that the scrotal testis is a derived state in mammals, Werdelin and Nilsonne  and Sharma et al.  argue that the scrotal CDT is the ancestral state and was subsequently lost in separate lineages. Overall, our reconstruction of the evolutionary history of testicular descent and the presence of the scrotum relied on a high order-level species coverage and well-accepted phylogeny, providing more credible information about the plesiomorphic scrotal testis in mammals. Moreover, the finding that the ascrotal testis evolved independently in mammals is a fascinating instance of convergent evolution.
Genetic mechanisms driving the evolution of testicular descent in mammals
The two phases of testicular descent are regulated by a combination of mechanical and hormonal factors . Dysplasia of the gubernacula, a pair of structures needed to steer the descent, yielded abnormal testicular descent. Previous studies have also found that mesenchymal, fibroblast, and muscle cells are involved in gubernaculum development [21, 22].
In the present study, we identified a series of rapidly evolving genes in ascrotal IDT and UDT mammals that were enriched in functions related to muscle and striated muscle development (Fig. 2, Additional file 1: Table S3, and Table S6). Positively selected genes in the group of ascrotal species were also overrepresented in GO terms related to muscular components (Additional file 1: Table S9 and Additional file 2: Fig S4). These findings imply that evolutionary changes in genes related to the muscular gubernaculum were involved in testicular descent. This is similar to the finding from Barthold et al.  that the differentially expressed genes of the gubernaculum in wild-type and cryptorchid rat fetuses were enriched in categories related to muscle development. Combined, this evidence suggests that muscle-related genes played some important role in the evolution of testicular descent in mammals, specifically by contributing to mechanical traction.
The process of testicular descent also involves a series of hormones. T, which is produced by Leydig cells, could induce CSL regression and cause the testis to descend into the first phase. More importantly, during the second phase, androgens (including T) are the predominant hormonal controllers . Moreover, T plays a large role in a paracrine manner to stabilize the Wolffian ducts and further systemically masculinize the external genitalia .
Desert hedgehog (Dhh) (Hedgehog protein family) has been reported to play an important role in regulating T secretion and testis development. Dhh is produced by Sertoli cells and can regulate the proliferation and differentiation of Leydig cells and functions of T secretion [34, 35]. A missense mutation in the rat Hedgehog pathway resulted in androgen deficiency and a decrease in the number of Leydig cells with impaired functions . Our finding that the rapidly evolving genes in ascrotal mammals were overrepresented in the Hedgehog signaling pathway (Fig. 2, Additional file 1: Table S3, and Table S5) aligns with the critical role of the Hedgehog signaling pathway in T secretion and testis development.
Possible mechanisms for maintaining health in IDT and UDT mammals
The testis is a crucial producer of spermatozoa and sex hormones in males. However, cryptorchidism caused by testicular maldescent is a heterogeneous disorder associated with not only the macroscopic abnormal location of testes, but also further postnatal abnormalities. For example, cryptorchidism is one of the most common causes of infertility in men; it is characterized by depression of spermatogenesis  and reduction in the number of germ cells . Abnormal autophagy  and DNA damage  were suggested to be involved in the impairment of spermatozoa in the cryptorchid testis. The present study uncovered several lines of evidence for strategies by which ascrotal IDT and UDT mammals keep normal reproductive capacity under “cryptorchid conditions.”
First, genes involved in spermatogenesis were associated with healthy UDT. G1021V in the protein GLI1 and E/S90N in KITLG were identified as UDT mammal-specific amino acid mutations. The zinc finger transcription factor GLI1 is a member of the GLI family, whose overexpression or mutation causes disease in humans and mice [41, 42]. There is evidence that GLI1 mediates Desert hedgehog (Dhh) signaling in male mouse testes during spermatogenesis . The KIT/KITLG signaling system is essential for the proliferation, meiosis, migration, survival, and maturation of germ cells in testes [44, 45]. It has been reported that polymorphisms in KITLG are likely associated with germ cell tumors . Although residue 1021 of GLI1 is not located in a putative functional region and nonsynonymous mutation of residue 90 is an unreported polymorphism case, both might be targeted in response to spermatogenesis via unknown cis-regulations.
We also found that the genes ATRX and NANOS1, which function in spermatogenesis in the testis, evolved an association with testicular descent. ATRX (on the X-chromosome, responsible for alpha-thalassaemia and mental retardation) is believed to play a role in testicular development, since the majority of mutations in ATRX result in genital abnormalities . Tang et al.  suggested that ATRX might also contribute to adult spermatogenesis in human. Mutations in the gene NANOS1 are associated with spermatogenic failure and oligoasthenoteratozoospermia .
Of the genes with UDT mammal-specific indels, T-complex 11 like 1 (TCP11L1) is involved in cryptorchidism. Seabra et al.  described a cryptorchidism patient with azoospermia presenting a microdeletion at 11p13 in TCP11L1. Liu et al.  found that TCP11L1 shares its functional domain with and has similar subcellular localization to the TCP11 protein, suggesting that the two proteins have a similar role in spermatogenesis. Several studies have demonstrated that, similar to mutations and insertions, deletions in a protein region might enhance that protein’s molecular function if the region effects it; this also occurs with cis- or trans-acting regulation [52, 53]. It is likely that this unique deletion of TCP11L1 in UDT species enhances spermatogenic capability, but further functional assays are needed to confirm this.
Second, genes contributing to DNA repair and genome stability maintenance were found to be involved in potential mechanisms maintaining health in UDT mammals. Germ cells pass genetic information onto descendants, and this requires that replication is accurate and their genome is perfectly stable . Obstacles to replication might cause genomic instability and cancer formation . However, the physiological environment provided by the undescended testis challenges the enzymes and cellular mechanisms that appear to be well adapted to the lower temperature in the CDT . We found that PMS2 possesses a UDT species-specific mutation at position 258. A previous study found that PMS2 plays a crucial role not only in the post-replicative DNA mismatch repair system, but also in a process that induces cell cycle arrest and could lead to apoptosis in the case of major DNA damages—this is known as DNA damage signaling . Along with DNA replication repair, sister chromatid cohesion in the meiotic process is another crucial mechanism for maintaining genome stability . Our results showed that genes that evolved in association with testicular descent were significantly enriched in functions that maintain sister chromatid cohesion (Additional file 1: Table S8 and Additional file 2: Fig S3), a biological process in which sister chromatids of a replicated chromosome become tethered to each other. Cohesion in eukaryotic cells appears to lie at the heart of the meiotic process because it is compulsory for the repair of recombinogenic lesions and for chromosome segregation in dividing cells during meiotic anaphase [59, 60].
Third, genes involved in autophagy might help maintain health in UDT mammals. Well-adjusted germ cell proliferation and death are highly ordered in testicular spermatogenesis [61, 62]. However, testicular heating suppresses spermatogenesis and leads to increasing germ cell degeneration and death [63, 64]. Autophagy involves a biological process of self-cannibalization via lysosomal degradation, namely nonapoptotic cellular demise . Zheng et al.  found that spermatogenesis was impaired in a surgery-induced cryptorchid mouse model; they also found cryptorchidism-induced autophagy and apoptosis synchronously promoting germ cell death. Notably, there have hardly been any reports of over-autophagy in naturally ascrotal mammals. Therefore, we believe that, to maintain homeostasis in testes and germ cells, ascrotal species might adjust their autophagy levels to adapt to the high ambient temperature of the undescended testes. This is supported by the significant correlation of autophagy associated genes TMEM74 and RAB3GAP2 in the PGLS analysis and a rapidly evolving gene RAB7A in UDT mammals (Table 1 and Fig. 2). Transmembrane protein 74 (TMEM74) might be an important element promoting autophagy under cell stress conditions, because the knockdown of this gene hampers the cell’s autophagy function when starvation is imposed . RAB3GAP2 and RAB7A are encoded by family members of RAB, which is a group of GTP-binding proteins that regulate vesicular transport. Both genes are involved in the autophagy pathway [68, 69].
Molecular convergence of the ascrotal IDT and UDT in GTPase
For decades, researches have demonstrated that phenotypic convergence can result from convergent molecular mechanisms . In our study, phylogenetic reconstruction revealed that the ascrotal testis evolved in different mammalian lineages independently, and that this could be regarded as convergent evolution. Genes harboring convergent/parallel substitutions in the ascrotal mammals were found to be significantly enriched in small GTPase-related terms (Additional file 1: Table S12 and Fig. 5). Small GTPases are a large family of molecular switches that play a pivotal role in various cellular processes . Specifically, recent advances have suggested that Ras GTPases are involved in testicular descent . Syndromes such as Noonan, Cardiofaciocutaneous, LEOPARD, and Costello, which are characterized as cryptorchidism, are RASopathies . Moreover, GTPase-mediated signal transduction was identified in the functional analysis of the differentially expressed genes of wild-type and cryptorchid rats . A previous study argues that Rho GTPases are critical for cytoskeletal reorganization and myogenesis . This is consistent with the roles of INSL3/RXFP2 and other candidate genes in regulating myogenesis and muscle development in the development of gubernaculum during testicular descent . Hence, it could be hypothesized that small GTPases drive the convergent evolution of the ascrotal testis in different mammals, playing an important role in the development of gubernaculum via myogenesis.
Nevertheless, the development of testicular descent is a complicated process involving many genes and pathways. Future studies using genome-wide scans are needed to elucidate the molecules and mechanisms dictating the evolution of testicular descent in mammals and test the candidate genes by functional analyses. Additionally, primates, which have a relatively low prevalence in the evolutionary history of the Mammalia , were overrepresented in this analysis. A more precise and comprehensive sampling that is proportional to the distribution of species in the Mammalia class is needed to remedy any potential ascertainment bias in the present study.
Our study combined phenotypic evolution and comparative genomics investigations of 380 cryptorchidism-related genes in mammals and found that the scrotal testis is the ancestral state in mammals and the ascrotal testis evolved multiple times independently. More importantly, the rapidly evolving and positively selected genes that we found suggest that the derived status of the ascrotal IDT and UDT phenotypes in mammals can be attributed to the adaptive evolution of genes involved in testicular descent and muscle development. Moreover, we demonstrated that accurate DNA replication, high genome stability, and appropriate autophagy are likely the mechanisms by which ascrotal mammals maintain normal spermatogenesis and physiological health. Our results suggest that small GTPases are associated with molecular mechanisms that contribute to the convergent phenotype of ascrotal testes in different mammals. Our study provides some novel insights into the evolution of testicular descent in mammals, contributes to a better understanding of the pathogenesis of cryptorchidism in humans, and offers further experimental validation for these candidate genes and residues.
Data on the presence and testicular position of the scrotums of 62 representative mammals in the UCSC 100-way multiple alignments (http://hgdownload.soe.ucsc.edu/downloads.html) and the monotreme echidna (Tachyglossus aculeatus) were taken from publications (Additional file 1: Table S1). Three categories of testicular position were classified: completely descended testis (CDT), incompletely descended testis (IDT), and undescended testis (UDT) (Fig. 1 and Additional file 1: Table S1). 380 cryptorchidism-related protein coding genes screened from The Cryptorchidism Gene Database version 3  were collected. Homologous exon alignments of these cryptorchidism-related genes were gathered from UCSC human 100-way multiple alignments at http://genome.ucsc.edu/. The entire transcripts were concatenated by exon-employing custom perl scripts. The longest transcript was retained for genes with multiple splice variants. To improve the quality of sequence alignments and subsequent evolutionary analyses, the nucleotide alignments were realigned by PRANK v.170427 in codon mode . Sequences were then treated with Gblocks . Incomplete codons and premature stop codons were prohibited. Based on the criteria that (1) at least one representative species for each main mammalian order was chosen, (2) all the focused “cryptorchidism” (IDT + UDT) species were included, and (3) species with higher quality genomes were kept, several CDT species which were mainly from the overrepresented primates were removed. Sequences from 49 mammals (28 CDT species, 15 IDT, and 6 UDT species) covering 16 orders (Primates, Rodentia, Lagomorpha, Cetartiodactyla, Perissodactyla, Carnivora, Chiroptera, Eulipotyphla, Proboscidea, Sirenia, Afrosoricida, Macroscelidea, Tubulidentata, Cingulata, Marsupialia, and Monotremata) (Additional file 1: Table S1, Additional file 2: Fig S5) were used in the evolutionary analyses described below.
Ancestral state reconstruction
The discrete morphologies of the scrotum and testis used in this study for extant taxa were assigned according to previous studies (Additional file 1: Table S1). Ancestral character reconstructions of mammalian scrotum presence and testis position were performed via the phytools and ape packages in R statistical software [78, 79]. We first tested the fitness of the following models: Equal Rate model (ER): discrete character evolution in which a single parameter governs all transition rates; All Rates Different model (ARD): all possible transitions can occur at different rates; and Symmetrical Rates model (SYM): forward and reverse transitions share the same parameter. We then used the model with the best AIC score to reconstruct the ancestral state of mammalian scrotum presence and testis position. We used a well-supported mammal phylogeny from TimeTree (http://www.timetree.org/)  as the input tree.
Selective pressure test
Rapidly evolved and positively selected genes generally contribute to adaptive evolution in response to new and changed environments . To test whether heterogeneity in the evolutionary rates of cryptorchidism-related genes among mammals are correlated with different testis positions, we undertook a suite of analyses to detect differences in selection intensity and positive selection using codeml in the PAML package .
A false discovery rate (FDR) correction was applied to conservatively account for multiple testing. First, the branch model was recruited to examine ω (the ratio of nonsynonymous to synonymous substitution rates) values between each of two groups of branches. A nested model comparison between ascrotal and scrotal testis mammals—including a two-ratio model (in which scrotal and ascrotal mammals have different ω values) and a one-ratio model (in which all mammals have one ω value)—was used to evaluate rapidly evolving genes in ascrotal IDT and UDT mammals. Another branch-specific comparison was made between IDT and UDT species within ascrotal mammals.
Second, the branch-site model was used to detect positively selected genes and amino acids across the combined IDT and UDT branches. The ω value was used to describe selection pressure: ω < 1, ω = 1, and ω > 1 indicated purified selection, neutral evolution, and positive selection, respectively.
Association analysis between sequence evolution and phenotype
To investigate the potential relationship between gene evolutionary rates and testicular and scrotal evolution in mammals, we employed a PGLS regression in the caper package of R . We used the root-to-tip ω value to access the evolutionary rate of each coding sequence because it contains more evolutionary history than does the terminal-branch ω value . The free ratio model of the codeml program in PAML was used to estimate the ω values of internal and terminal branches. We used a binary state combination of the scrotal CDT and ascrotal IDT and UDT, with relevant characters shown in Additional file 1: Table S1. The lambda (λ) value, used as a quantitative measure of phylogenetic signals, was estimated by the maximum likelihood method .
Shared specific amino acids and convergent/parallel substitutions in UDT species
Species-specific amino acid mutations and peptide fragment indel (insertion and deletion) events have been inferred to be associated with functional, physiological, and phenotypic changes [85, 86]. For each amino acid site within a protein, the shared UDT species-specific mutations refer to the amino acids along all UDT branches that are different from other species. The shared specific substitutions were identified by strict identity. To identify shared specific amino acid mutations in the UDT species, in-house perl scripts were employed on each column of the trimmed amino acid alignments.
It has been suggested that convergent phenotypes can be mediated by parallel and convergent amino acid substitutions [87,88,89]. Convergent substitution occurs when different amino acids at a specific amino acid site from two distantly-related branches converge, whereas parallel substitution occurs when the same amino acid in two branches are derived from another same amino acid in independent ancestors . Convergent/parallel amino acid substitutions in ascrotal branches were detected via the method described in Zou and Zhang .
For each gene, the ancestral node sequence was first reconstructed using the codeml program in the PAML package . Second, a tree with branch lengths was extracted from the output file from the abovementioned calculation. Then, the relative evolutionary rates of all amino acid sites within one gene were calculated using the aaml model in codeml. The frequency of each amino acid at each site was counted for each gene. Using this method , all observed and expected cases of convergent and parallel amino acid substitutions were calculated under the recommended JTT-fgene matrix. We focused on and filtered the convergent/parallel substitutions between each of two ascrotal branches. Finally, a Poisson test was employed to identify the significance between the observed and expected numbers of substitutions in ascrotal branches.
We used the over-representation tests enrichGO and enrichKEGG incorporated into clusterProfiler version 3.6.0  to analyze Gene Ontology (GO) and KEGG enrichment, respectively. Benjamini and Hochberg (BH) multiple test correction  was performed to adjust p values and decrease the likelihood of false positives.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its supplementary information files.
Completely descended testis
Incompletely descended testis
Calcitonin gene-related peptide
Equal rate model
All rates different model
Symmetrical rates model
University of California, Santa Cruz
False discovery rate
Phylogenetic generalized least squares
Kyoto Encyclopedia of Genes and Genomes
Benjamini and Hochberg
Cranial suspensory ligament
Moore CR. The biology of the mammalian testis and scrotum. Q Rev Biol. 1926;1:4–50.
Bedford JM. Anatomical evidence for the epididymis as the prime mover in the evolution of the scrotum. Am J Anat. 1978;152:483–507.
Zorgniotti AW. Temperature and environmental effects on the testis. New York: Plenum Press; 1991.
Klonisch T, Fowler PA, Hombach-Klonisch S. Molecular and genetic regulation of testis descent and external genitalia development. Dev Biol. 2004;270:1–18.
Durairajanayagam D, Agarwal A, Ong C. Causes, effects and molecular mechanisms of testicular heat stress. Reprod Biomed Online. 2015;30:14–27.
Simoni M, Huhtaniemi IT. Endocrinology of the testis and male reproduction. Cham: Springer; 2017.
Freeman S. The evolution of the scrotum: a new hypothesis. J Theor Biol. 1990;145:429–45.
Portmann A. Animal forms and patterns. 1967.
Frey VR. Zur Ursache des Hodenabstiegs (Descensus testiculorum) bei Säugetieren 1. J Zool Syst Evol Res. 1991;29:40–65.
Chance MRA. Reason for externalization of the testis of mammals. J Zool. 1996;239:691–5.
Toppari J, Rodprasert W, Virtanen HE. Cryptorchidism-disease or symptom? Ann Endocrinol. 2014;75:72–6.
Hutson JM, Thorup JM, Beasley SW. Descent of the testis. London: Edward Arnold; 1992.
Grant T, Fanning D. The platypus. Kensington: New South Wales University Press; 1984.
Werdelin L, Nilsonne Å. The evolution of the scrotum and testicular descent in mammals: a phylogenetic view. J Theor Biol. 1999;196:61–72.
Kleisner K, Ivell R, Flegr J. The evolutionary history of testicular externalization and the origin of the scrotum. J Biosci. 2010;35:27–37.
Lovegrove BG. Cool sperm: why some placental mammals have a scrotum. J Evol Biol. 2014;27:801–14.
Sharma V, Lehmann T, Stuckas H, Funke L, Hiller M. Loss of RXFP2 and INSL3 genes in Afrotheria shows that testicular descent is the ancestral condition in placental mammals. PLoS Biol. 2018;16:e2005293.
Hutson JM. A biphasic model for the hormonal control of testicular descent. Lancet. 1985;326:419–21.
Hutson JM, Hasthorpe S, Heyns CF. Anatomical and functional aspects of testicular descent and cryptorchidism. Endocr Rev. 1997;18:259–80.
Emmen JM, McLuskey A, Adham IM, Engel W, Grootegoed JA, Brinkmann AO. Hormonal control of gubernaculum development during testis descent: gubernaculum outgrowth in vitro requires both insulin-like factor and androgen. Endocrinology. 2000;141:4720–7.
Heyns CF. The gubernaculum during testicular descent in the human fetus. J Anat. 1987;153:93.
Costa WS, Sampaio FJ, Favorito LA, Cardoso LE. Testicular migration: remodeling of connective tissue and muscle cells in human gubernaculum testis. J Urology. 2002;167:2171–6.
Niedzielski JK, Oszukowska E, Słowikowska-Hilczer J. Undescended testis–current trends and guidelines: a review of the literature. Arch Med Sci. 2016;12:667.
Lewis LG. Cryptorchism. J Urology. 1948;60:345–56.
Larkins SL, Hutson JM, Williams MP. Localisation of calcitonin gene-related peptide immunoreactivity within the spinal nucleus of the genitofemoral nerve. Pediatr Surg Int. 1991;6:176–9.
Kumar S, Stecher G, Suleski M, Hedges SB. TimeTree: a resource for timelines, timetrees, and divergence times. Mol Biol Evol. 2017;34:1812–9.
Bagheri-Fam S, Argentaro A, Svingen T, Combes AN, Sinclair AH, Koopman P, Harley V. Defective survival of proliferating Sertoli cells and androgen receptor function in a mouse model of the ATR-X syndrome. Hum Mol Genet. 2011;20:2213–24.
Lu C, Xu M, Wang Y, Qin Y, Du G, Wu W, Han X, Ji C, Yang Y, Gu A, Xia Y, Song L, Wang S, Wang X. Genetic variants in meiotic program initiation pathway genes are associated with spermatogenic impairment in a Han Chinese population. PLoS ONE. 2013;8:e53443.
Short RV. The testis: the witness of the mating system, the site of mutation and the engine of desire. Acta Paediatr. 1997;86:3–7.
Moore KL, Persaud TVN, Torchia MG. The developing human-e-book: clinically oriented embryology. Amsterdam: Elsevier Health Sciences; 2018.
Bokma F. Detection of “punctuated equilibrium” by Bayesian estimation of speciation and extinction rates, ancestral character states, and rates of anagenetic and cladogenetic evolution on a molecular phylogeny. Evolution. 2008;62:2718–26.
Barthold JS, McCahan SM, Singh AV, Knudsen TB, Si X, Campion L, Akins RE. Altered expression of muscle-and cytoskeleton-related genes in a rat strain with inherited cryptorchidism. J Androl. 2008;29:352–66.
Hughes IA, Acerini CL. Factors controlling testis descent. Eur J Endocrinol. 2008;159(suppl 1):S75–82.
Szczepny A, Hime GR, Loveland KL. Expression of hedgehog signaling components in adult mouse testis. Dev Dyn. 2006;235:3063–70.
Varjosalo M, Taipale J. Hedgehog: functions and mechanisms. Gene Dev. 2008;22:2454–72.
Kawai Y, Noguchi J, Akiyama K, Takeno Y, Fujiwara Y, Kajita S, Tsuji T, Kikuchi K, Kaneko H, Kunieda T. A missense mutation of the Dhh gene is associated with male pseudohermaphroditic rats showing impaired Leydig cell development. Reproduction. 2011;141:217–25.
Hadziselimovic F. Early successful orchidopexy does not prevent from developing azoospermia. Int Braz J Urol. 2006;32:570–3.
Cortes D, Clasen-Linde E, Hutson JM, Li R, Thorup J. The Sertoli cell hormones inhibin-B and anti Müllerian hormone have different patterns of secretion in prepubertal cryptorchid boys. J Pediatr Surg. 2016;51:475–80.
Yefimova MG, Buschiazzo A, Burel A, Lavault MT, Pimentel C, Jouve G, Jaillard S, Jegou B, Bourmeyster N, Ravel C. Autophagy is increased in cryptorchid testis resulting in abnormal spermatozoa. Asian J Androl. 2019;21:570.
Song SH, Lim JJ, Bang JK, Cha SK, Lee DR, Kim YS, Ahn TY, Yoon TK. Evaluation of sperm deoxyribonucleic acid (dna) damage and effects on embryo development using a mouse cryptorchidism model. Urology. 2013;82:743.e17-743.e23.
Altaba AR. Gli proteins and Hedgehog signaling: development and cancer. Trends Genet. 1999;15:418–25.
Villavicencio EH, Walterhouse DO, Iannaccone PM. The sonic hedgehog-patched-gli pathway in human development and disease. Am J Hum Genet. 2000;67:1047–54.
Kroft TL, Patterson J, Won Yoon J, Doglio L, Walterhouse DO, Iannaccone PM, Goldberg E. GLI1 localization in the germinal epithelial cells alternates between cytoplasm and nucleus: upregulation in transgenic mice blocks spermatogenesis in pachytene. Biol Reprod. 2001;65:1663–71.
Runyan C, Schaible K, Molyneaux K, Wang Z, Levin L, Wylie C. Steel factor controls midline cell death of primordial germ cells and is essential for their normal proliferation and migration. Development. 2006;133:4861–9.
Cheng P, Chen H, Liu S, Pu X. SNPs in KIT and KITLG genes may be associated with oligospermia in Chinese population. Biomarkers. 2013;18:650–4.
Qin J, Cai Z, Yang Y, Xing J, Zhuang X. Association between KITLG gene polymorphisms and testicular germ cell tumors: a systematic review and meta-analysis. Curr Med Sci. 2019;39:702–6.
Gibbons RJ, Higgs DR. Molecular-clinical spectrum of the ATR-X syndrome. Am J Med Genet. 2000;97:204–12.
Tang P, Argentaro A, Pask AJ, O’Donnell L, Marshall-Graves J, Familari M, Harley VR. Localization of the chromatin remodelling protein, ATRX in the adult testis. J Reprod Dev. 2011;57:317–21.
Kusz-Zamelczyk K, Sajek M, Spik A, Glazar R, Jędrzejczak P, Latos-Bieleńska A, Kotecki M, Pawelczyk L, Jaruzelska J. Mutations of NANOS1, a human homologue of the Drosophila morphogen, are associated with a lack of germ cells in testes or severe oligo-astheno-teratozoospermia. J Med Genet. 2013;50:187–93.
Seabra CM, Quental S, Neto AP, Carvalho F, Gonçalves J, Oliveira JP, Fernandes S, Sousa M, Barros A, Amorim A, Lopes AM. A novel Alu-mediated microdeletion at 11p13 removes WT1 in a patient with cryptorchidism and azoospermia. Reprod Biomed Online. 2014;29:388–91.
Liu Y, Wang X, Sun H, Lu Y, Zeng M, Tao D, et al. Expression and structure analysis of a novel mouse gene Tcp11l1. Chin J Appl Environ Biol. 2008;14:783–6.
Ozaki T, Naka M, Takada N, Tada M, Sakiyama S, Nakagawara A. Deletion of the COOH-terminal region of p73α enhances both its transactivation function and DNA-binding activity but inhibits induction of apoptosis in mammalian cells. Cancer Res. 1999;59:5902–7.
Niu XL, Peters KG, Kontos CD. Deletion of the carboxyl terminus of Tie2 enhances kinase activity, signaling, and function evidence for an autoinhibitory mechanism. J Biol Chem. 2002;277:31768–73.
Helle F. Germ cell DNA-repair systems—possible tools in cancer research? Cancer Gene Ther. 2012;19:299.
Hsieh HJ, Peng G. Cellular responses to replication stress: implications in cancer biology and therapy. DNA Repair. 2016. https://doi.org/10.1016/j.dnarep.2016.11.002.
Steinberger A. Effects of temperature on the biochemistry of the testis. In: Temperature and environmental effects on the testis. Boston: Springer; 1991. p. 33–47.
Borràs E, Pineda M, Cadiñanos J, del Valle J, Brieger A, Hinrichsen I, et al. Refining the role of PMS2 in Lynch syndrome: germline mutational analysis improved by comprehensive assessment of variants. J Med Genet. 2013;50:552–63.
Morrison C, Vagnarelli P, Sonoda E, Takeda S, Earnshaw WC. Sister chromatid cohesion and genome stability in vertebrate cells. Biochem Soc T. 2003;31:263–5.
Klein F, Mahr P, Galova M, Buonomo SB, Michaelis C, Nairz K, Nasmyth K. A central role for cohesins in sister chromatid cohesion, formation of axial elements, and recombination during yeast meiosis. Cell. 1999;98:91–103.
Peters JM, Tedeschi A, Schmitz J. The cohesin complex and its roles in chromosome biology. Gene Dev. 2008;22:3089–114.
Print CG, Loveland KL. Germ cell suicide: new insights into apoptosis during spermatogenesis. BioEssays. 2000;22:423–30.
Zhang M, Jiang M, Bi Y, Zhu H, Zhou Z, Sha J. Autophagy and apoptosis act as partners to induce germ cell death after heat stress in mice. PLoS ONE. 2012;7:e41412.
Rockett JC, Mapp FL, Garges JB, Luft JC, Mori C, Dix DJ. Effects of hyperthermia on spermatogenesis, apoptosis, gene expression, and fertility in adult male mice. Biol Reprod. 2001;65:229–39.
Lue Y, Wang C, Liu Y, Hikim APS, Zhang X, Ng C, Hu Z, Li Y, Leung A, Swerdloff RS. Transient testicular warming enhances the suppressive effect of testosterone on spermatogenesis in adult cynomolgus monkeys (Macaca fascicularis). J Clin Endocr Metab. 2006;91:539–45.
Levine B, Yuan J. Autophagy in cell death: an innocent convict? J Clin Invest. 2005;115:2679–88.
Zheng Y, Zhang P, Zhang C, Zeng W. Surgery-induced cryptorchidism induces apoptosis and autophagy of spermatogenic cells in mice. Zygote. 2019;27:101–10.
Yu C, Wang L, Lv B, Lu Y, Chen Y, Ma D, Shi T, Wang L. TMEM74, a lysosome and autophagosome protein, regulates autophagy. Biochem Bioph Res Co. 2008;369:622–9.
Zhang M, Chen L, Wang S, Wang T. Rab7: roles in membrane trafficking and disease. Biosci Rep. 2009;29:193–209.
Spang N, Feldmann A, Huesmann H, Bekbulat F, Schmitt V, Hiebel C, et al. RAB3GAP1 and RAB3GAP2 modulate basal and rapamycin-induced autophagy. Autophagy. 2014;10:2297–309.
Foote AD, Liu Y, Thomas GW, Vinař T, Alföldi J, Deng J, et al. Convergent evolution of the genomes of marine mammals. Nat Genet. 2015;47:272.
Bourne HR, Sanders DA, McCormick F. The GTPase superfamily: conserved structure and molecular mechanism. Nature. 1991;349:117–27.
Cannistraci CV, Ogorevc J, Zorc M, Ravasi T, Dovc P, Kunej T. Pivotal role of the muscle-contraction pathway in cryptorchidism and evidence for genomic connections with cardiomyopathy pathways in RASopathies. BMC Med Genomics. 2013;6:5.
Bryan BA, Li D, Wu X, Liu M. The Rho family of small GTPases: crucial regulators of skeletal myogenesis. Cell Mol Life Sci. 2005;62:1547–55.
Burgin CJ, Colella JP, Kahn PL, Upham NS. How many species of mammals are there? J Mammal. 2018;99:1–14.
Urh K, Kunej T. Genome-wide screening for smallest regions of overlaps in cryptorchidism. Reprod Biomed Online. 2018;37:85–99.
Löytynoja A, Goldman N. Phylogeny-aware gap placement prevents errors in sequence alignment and evolutionary analysis. Science. 2008;320:1632–5.
Castresana J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000;17:540–52.
Paradis E, Claude J, Strimmer K. APE: analyses of phylogenetics and evolution in R language. Bioinformatics. 2004;20:289–90.
Revell LJ. Phytools: an R package for phylogenetic comparative biology (and other things). Methods Ecol Evol. 2012;3:217–23.
Swanson WJ, Vacquier VD. The rapid evolution of reproductive proteins. Nat Rev Genet. 2002;3:137.
Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007;24:1586–91.
Orme D, Freckleton R, Thomas G, Petzoldt T. The caper package: comparative analysis of phylogenetics and evolution in R. R package version. 2013;5:1–36.
Montgomery SH, Capellini I, Venditti C, Barton RA, Mundy NI. Adaptive evolution of four microcephaly genes and the evolution of brain size in anthropoid primates. Mol Biol Evol. 2010;28:625–38.
Pagel M. Inferring the historical patterns of biological evolution. Nature. 1999;401:877.
Frazer KA, Chen X, Hinds DA, Pant PK, Patil N, Cox DR. Genomic DNA insertions and deletions occur frequently between humans and nonhuman primates. Genome Res. 2003;13:341–6.
Wierer M, Schrey AK, Kühne R, Ulbrich SE, Meyer HH. A single glycine-alanine exchange directs ligand specificity of the elephant progestin receptor. PLoS ONE. 2012;7:e50350.
Feldman CR, Brodie ED, Pfrender ME. Constraint shapes convergence in tetrodotoxin-resistant sodium channels of snakes. Proc Natl Acad Sci USA. 2012;109:4556–61.
Storz JF. Causes of molecular convergence and parallelism in protein evolution. Nat Rev Genet. 2016;17:239.
Hu Y, Wu Q, Ma S, Ma T, Shan L, Wang X, Nie Y, Ning Z, Yan L, Xiu Y, Wei F. Comparative genomics reveals convergent evolution between the bamboo-eating giant and red pandas. Proc Natl Acad Sci USA. 2017;114:1081–6.
Zou Z, Zhang J. Are convergent and parallel amino acid substitutions in protein evolution more prevalent than neutral expectations? Mol Biol Evol. 2015;32:2085–96.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16:284–7.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc B. 1995;57:289–300.
We thank all of our lab members and the members of Jiangsu Key Laboratory for Biodiversity and Biotechnology for their suggestions and support during this project. We also thank Noah Last of Third Draft Editing for his English language editing.
This research was supported by the State Key Program of National Science of China to Guang Yang (Grant Number 31630071), the National Natural Science Foundation of China to Wenhua Ren (Grant Number 31872219), the National Key Research and Development Project to Guang Yang (Grant Number 2016YFC0503204), and the National Natural Science Foundation of China to Shixia Xu (Grant Number 31772448). The funding bodies for this study had no role in the study design, collection of data, data analysis or interpretation, or writing the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The presence of scrotum and the testicular position of 63 representative mammals. Table S2. Model comparison of the evolution of testicular descent in the ancestral state reconstruction. Table S3. Rapidly evolved genes in ascrotal testis mammals. Table S4. Rapidly evolving genes in a 62-species mammalian data set. Table S5. KEGG enrichment of rapidly evolving genes in ascrotal mammals. Table S6. GO enrichment of rapidly evolving genes in ascrotal mammals. Table S7. Rapidly evolving genes in UDT mammals. Table S8. GO enrichment of genes evolved with significant regression with testicular descent. Table S9. GO enrichment of positively selected genes in ascrotal IDT and UDT mammals. Table S10. The observed convergent/parallel amino acid substitutions in the ascrotal IDT and UDT branches. Table S11. Parallel and convergent amino acid substitutions in ascrotal UDT and IDT mammals. Table S12. GO enrichment of genes exhibit parallel/convergent substitutions in IDT mammals.
The expanded 62-mammal data set. Fig S2. Rapidly evolving genes in UDT mammals. Fig S3. Functional enrichment of genes that evolved significantly correlated to character of testis position (p value < 0.05). Fig S4. Top 10 GO function enrichment of positive selected genes in ascrotal mammals. Fig S5. The working tree with numbered internal nodes for identifying convergent/parallel amino acid substitutions.
About this article
Cite this article
Chai, S., Tian, R., Bi, J. et al. Rapid evolution and molecular convergence in cryptorchidism-related genes associated with inherently undescended testes in mammals. BMC Ecol Evo 21, 22 (2021). https://doi.org/10.1186/s12862-021-01753-5
- Testicular descent
- Cryptorchidism-related genes
- Molecular convergence
- Rapid evolution