Pre-Diagnostic Circulating RNAs Networks Identify Testicular Germ Cell Tumour Susceptibility Genes
==================================================================================================

* Joshua Burton
* Trine B. Rounge
* Trine B. Haugen
* Marcin W. Wojewodzic

## Abstract

Testicular germ cell tumour (TGCT) is a malignancy with known inherited risk factors, affecting young men. We have previously identified several hundred circulating RNAs that were differentially expressed in pre-diagnostic serum samples from TGCT cases when compared to healthy controls. In this study we performed network preservation analyses of pre-diagnostic serum mRNA and miRNA. Hub genes, enriched functional pathways, and regulatory feature prediction were identified for all TGCT, seminoma, and non-seminoma cases separately, compared to controls. We identified *UBCA1, RCC1, FMR1, OSA3,* and *UBE2W* as hub genes associated with TGCT. The genes *OSA3* and *UBE2W* have previously been associated with testicular dysgenesis syndrome (TDS) disorders. Previously described TGCT susceptibility genes *TEX14*, *NARS2,* and *G3BP2* were identified as hub genes in both seminoma and non-seminoma networks. Furthermore, network module analysis showed prediction of transcription factors for oestrogen-related receptors. The overlap between network hub genes and TGCT susceptibility genes indicates a role in the progression from germ cell neoplasia in situ (GCNIS) to TGCT that should be further investigated.

Key words
*   Network analysis
*   pre-diagnostic
*   TGCT
*   testicular cancer
*   RNA

## Introduction

Testicular germ cell tumour (TGCT) is the most common form of testicular cancer affecting primarily younger males. Throughout the 20th and 21st century there have been rising rates of TGCT in developed countries. However, in recent years the incidence rates have begun to plateau in Northern and Western European countries with high incidence, in Norway at around 11 cases per 100,000 person-years (1–3). It is generally agreed that TGCT develop from a pre-malignancy, from germ cell neoplasia in situ (GCNIS), obtained through a loss of differentiation of gonocytes during foetal life which progress to a malignancy during postnatal puberty (4). There are around 80 risk loci for TGCT as well as environmental risk factors, and ethnicity also affects the likelihood of developing TGCT (5–8). Further disorders commonly associated with the testicular dysgenesis syndrome (TDS) include cryptorchidism, poor semen quality, and hypospadias (9). These disorders are known to be associated with higher risk of TGCT to varying degrees (10). However, the aetiology behind TGCT is still mainly unknown.

Advances in TGCT treatment have led to an increased 5-year survival rate, up to 95% (11), however, survivors of TGCT are still at risk of secondary cancers, cardiovascular issues, epigenetic alterations, and reduced longevity associated with cisplatin treatment (12–15). It is therefore important to diagnose TGCT early to reduce the rounds and dosage of cisplatin-based treatment needed and therefore potentially reduce the long-term side effects of cisplatin. Novel methods of analysing existing data are therefore needed.

Analyses of pre-diagnostic samples may increase understanding of the progression of TGCT. The pre-diagnostic serum samples from the Janus Serum Bank are well suited for such studies (16, 17). Our previous studies of pre-diagnostic samples from this biobank using traditional linear models with negative binomial distribution employed on count data, identified genes associated with the development of lung and TGCT (18, 19). A total of 818 circulating RNAs were found to be differentially expressed in pre-diagnostic TGCT cases when compared to controls. For instance, differences in expression were observed for the male fertility-related gene, *TEX101*, and the X-chromosome inactivation gene, *TSIX* (18).

Weighted gene correlation network analysis (WGCNA) is a systems biology approach that allows for increased comprehension of large and multidimensional transcriptomic datasets (20). WGCNA identifies gene co-expression patterns between multiple samples from various states. Gene co-expression can be used to identify a variety of elements such as future therapeutic targets or biomarkers as well as screening datasets for networks of genes relating to certain traits (21). WGCNA also allows for the comparison of such networks to show nodes and clusters that are preserved between differing states, such as histological subtypes. This is one of the major differences between standard RNA-Seq approaches, where linear models are deployed (limma, edgeR, or DESeq2) (22–24). WGCNA does not only consider the differentially expressed genes, instead opting to utilise the comprehensive gene list to search for inter-connectivity between clusters of genes (20). Using these methods, previous research has found related genes to specific cancers through mRNA and miRNA analyses, including lung adenocarcinoma (25), bladder cancer (26), and breast cancer (27).

The WGCNA can be used to obtain new information about the mechanisms underlying TGCT development. Alterations in the networks may potentially indicate the presence of a disease or precede the development of one (28, 29). Observed changes within gene networks also have the potential to be used for biomarker discovery (30). Extensive networks relying on multiple different biological components, such as transcription factors, mRNAs, miRNAs, miRNA targets, and have been demonstrated in multiple diseases, including TGCT development (31). The data in Mallik’s work originated from post-diagnostic samples of seminoma and non-seminomas and compared only TGCT cases without controls. Subtype-specific networks were constructed, and mRNAs, miRNA and TFs were identified. In both the seminoma and non-seminoma modules, the top miRNAs were onco-related miRNAs for several types of cancers. Common hub genes in Mallik’s work, *SP1* and *MYC,* were also identified between the networks constructed for seminomas and non-seminomas.

Post-diagnostic TGCT networks contains Phosphatase and Tensin Homolog (*PTEN*), targeted by three miRNAs, hsa-miR-141, hsa-miR-222, and hsa-miR-21, with hsa-miR-222 also regulating KIT and Tumor Protein 53 (*TP53*) (32). The complexity of the networks can often be overlooked when only differential expression is studied, therefore a more in-depth method of comparing disease states should be pursued in order to obtain a more complete understanding of what drives TGCT development.

We aimed to identify mRNA and miRNA networks associated with TGCT development in pre-diagnostic serum samples by using WGCNA. We compared networks from different TGCT histological types. We hypothesised that by using network analysis, we could propose genes and mechanisms related to TGCT development.

## Materials and Methods

### Data and sample description

We proceeded with the data in our previous study of pre-diagnostic RNA profile in serum from TGCT patients (18). In brief, we obtained samples from the Janus Serum Bank that were linked to registry data from the Cancer Registry of Norway to identify matched cases and controls. We included patients who developed TGCT up to 10 years after sample donation and control samples who were cancer-free within 10 years after sample donation. The serum samples were collected between 1972 – 2004 from participants in health surveys and red cross blood donors.

In total, 79 TGCT cases were identified that had blood donated up to 10 years prior to TGCT analyses. Of these, 52 were identified as seminoma, and 27 as non-seminoma. In addition, 111 matched controls that were cancer-free 10 years after blood donation were also retrieved. This sample set is part of a larger data collection named JanusRNA, described in detail by Langseth., et al (33). RNA profiles were generated by extracting RNAs from 400 µl serum, and small RNA sequencing libraries were produced with NEBNext kit (Cat. No E7300, New England Biolabs Inc.) and sequenced on the HiSeq 2500 (Illumina) as previously described (34).

### Bioinformatics analysis

Pre-processing of raw transcriptomics data was run using a high-performance computing cluster at the Cancer Registry of Norway. Detailed description can be found in a previously published article (34) ([https://github.com/sinanugur/sncRNA-workflow](https://github.com/sinanugur/sncRNA-workflow)). In brief, the reads generated an average depth of 18.4 million raw reads per sample. Adapters were removed with AdapterRemoval (v.2.1.7) with collapsed reads then being mapped to the human genome (hg38) with Bowtie2 (v2.2.9). Annotation set for RNAs was GENCODE (v.26) and for miRNAs was miRbase (v.22). We included mRNAs with at least five reads in more than 20% of the samples, and in the miRNA analysis we removed miRNAs that had no reads in more than 10% of the samples.

Controls and cases for the network construction were matched using age at sampling, time of sampling, and blood donation group, including separate matching for samples in the ‘0-7 years before diagnosis’ and ‘0-4 years before diagnosis’ groups as part of a sensitivity analysis (Figure 1 I). The same matching methods were performed between the subtypes, seminomas and non-seminomas during the network analysis based on TGCT histology. The *optmatch* R package (v0.9–11) (github.com/markmfredrickson/optmatch) was used to find these matching sets. Full overview of these procedures was given in a previous article (Burton et al., 2020).

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/01/26/2022.12.16.22283563/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2023/01/26/2022.12.16.22283563/F1)

Figure 1. 
I) Design of the study. Samples were identified from the Janus Serum Bank (JSB) in periods prior to diagnosis. II) Workflow for analysis. a) RNAs were initially sequenced from JSB serum sample, and sncRNAs were quantified. b) Samples were split based on cancer diagnosis and histology and matched with controls. c) mRNA and miRNA datasets were separated for different analyses. d, e & f) Quality control was performed on the datasets, and data was prepared for network analyses. g) WGCNA pipelines were used to perform network analyses. h) The preservation analyses were then performed on comparable networks. i) Enrichment analysis was then performed on relevant modules.

Statistical analyses were performed on the processed datasets. Firstly, counts were normalised using counts per million from the EdgeR package (22). The R package ComBat (35) was then used to determine if there was any batch effect in our data. A simplified WGCNA analysis on mRNA cases data which had undergone batch correction was run. Power graph analyses were compared with the full mRNA cases analysis which had not undergone batch correction.

### WGCNA pipeline steps to module preservation analysis

Weighted Gene Correlation Network Analysis (WGCNA) package (1.70-3) was used on the pre-processed data. First, we ensured scale free topology that is necessary assumption for network analysis by choosing an appropriate power value for all analyses. WGCNA was then run through the following steps described in the flowchart (Figure 1 II). Firstly, feature (mRNA or miRNA expression) correlations were used to construct the individual networks, followed by clusters identification of interconnected nodes hereafter referred to as modules. Then the nodes within the modules are summarised through a highly connected hub node (‘eigengenes’ or ‘eigenmirs’ for mRNA and miRNA network, respectively), alleviating multiple testing issues by reducing the dimensions of the datasets. The distinct individual modules, named by colours, can then be identified, and all nodes within the network can be annotated depending on how close they are to the modules. In addition, TGCT samples were divided by seminoma and non-seminoma samples, and network construction was performed as above, with only mRNA counts. This was due to the lack of sufficient miRNAs needed for network analysis in each histology.

### Preservation network analysis

Network preservation analysis was used for the comparison of two networks and the modules preserved across networks. We identified modules that are preserved across the two constructed networks and the degree to which they are preserved. For the network preservation analysis, we made several comparisons between groups, initially, contrasts between mRNA serum levels in cases and controls, followed by miRNA serum levels between cases and controls and mRNA serum levels between seminomas and non-seminomas. For each of these subsets we constructed separate networks, whilst also assuring the same scale-free topology. Module preservation statistics were computed using the modulePreservation function (200 permutations) implemented in WGCNA (36) and further implemented from Frisch’s work, (37). Network module preservation statistics quantify how density and connectivity patterns of modules defined in a reference dataset (control samples) are preserved in a test dataset (TGCT samples). Network adjacency comparisons are superior estimates of module preservation to standard cross-tabulation techniques (36). We used network adjacency comparison to assess preservation of gene co-expression patterns in control and TGCT modules, and between histology types. The overall significance of the preservation was assessed using Zsummary and median rank statistics (36). Based on the thresholds proposed in the documentation, Zsummary <2 indicates no preservation, 2 < Zsummary < 10 weak to moderate evidence of preservation, and Zsummary >10 strong evidence of module preservation across networks. Median rank statistics further consider module size for preservation assessment (36).

### Pathway analysis

Functional annotation was performed using the enrichment analysis tool Enrichr (38–40), by inputting gene lists from modules generated by WGCNA. Of particular interest were modules that were not preserved between cases and controls and those which were preserved between seminomas and non-seminomas. Enrichr databases were then reviewed for significantly enriched pathways and pathways associated with TGCT or TDS conditions.

### Motive finder, transcription factor mRNA

The regulatory patterns of the mRNAs included in the modules were investigated using a promoter motif analysis tool i-cisTarget ([https://gbiomed.kuleuven.be/apps/lcb/i-cisTarget/](https://gbiomed.kuleuven.be/apps/lcb/i-cisTarget/) Last accessed 6th June 2022). A full analysis with standard parameters was run with gene symbols as the input and RefSeq r45 and v.6 of the i-cisTarget database, with only TF binding site databases included in the analysis. This analysis could reveal if the module is regulated by specific transcription factors for genes related to male fertility or cancer development.

### Determining miRNA targets on mRNA

New modules predicting miRNA as preserved between cases and controls underwent a miRNA prediction analysis to determine which mRNAs the preserved miRNAs targeted. By translating miRNA to mRNA targets it also allows for easier comparison between the miRNA and mRNA networks. MiRNA targets are determined through sequence matching. Targets were extracted from miRDB (v5.0) with a cut off score of > 70 allowing only targets of high confidence to be included in the analysis.

### Enrichment analyses

Enrichment analyses were performed on the network and module data, firstly, KEGG analysis was performed, using Enrichr once more, on the top ten hub genes in each of the mRNA cases, seminoma and non-seminoma modules. A parallel enrichment analysis was also performed on all the genes in the module to give an overall prediction for enriched pathways. For miRNA cases, mRNA targets were used for enrichment analysis, with the cut off score used previously (>70).

### Cytoscape visualisation

Cytoscape (v. 3.9.0) was used to visualise the network modules. Cytoscape files were generated using the R package WGCNA’s export to Cytoscape function, generating both the node and the edge files necessary to visualise. Additional layout packages from yFiles Layout Algorithms (v. 1.1.2) were used for visualisation of the networks.

### Sensitivity analysis

In order to determine if the time to diagnosis of each sample affected network construction, a sensitivity analysis was performed by excluding the samples with a longer time period between sampling and diagnosis. As the data was split into four separate time frames from a previous study ([0-2], [2-5], [5-8], [8-10]), we first excluded the 8-10 years between sampling and diagnosis samples, which also gives us smaller datasets to work with. As a further sensitivity step, the 5-8 years between sampling and diagnosis samples were removed for a second round.

## Results

### Power calculations

To ensure scale-free topology, we performed power analysis and selected a power value equal to 10 for mRNA networks and miRNA networks (Sup F1).

#### Main characteristics of the networks

Network statistics showed that the mRNA cases network consisted of 19 modules, containing a total of 10596 mRNAs, with an average module density of 0.026 and an average heterogeneity of 0.99. The control network for mRNA consisted of 11 modules containing 8462 mRNAs, with an average module density of 0.028 and an average heterogeneity of 1.06.

The miRNA cases network consisted of five modules containing 403 miRNAs, with an average module density of 0.09 and average heterogeneity of 0.52, whilst the miRNA control network consisted of four modules containing a total of 403 miRNAs, with an average module density of 0.1 and an average heterogeneity of 0.49.

The seminoma network consisted of 22 modules containing 11032 mRNAs, with an average module density of 0.05 and average heterogeneity of 0.72, whilst the non-seminoma network consisted of 26 modules containing 11026 mRNAs, with an average module density of 0.06 and an average heterogeneity of 0.72.

### Modules of interest in mRNA cases vs. control

Modules with TGCT-related genes in mRNA cases include the black(cases), magenta(cases), purple(cases), salmon(cases), and the midnightblue(cases) modules (Table 5). The midnightblue(cases) module consisted of 169 genes with a density of 0.03, and the top three eigengenes for this module were *RTCB*, *JAKMIP1,* and *MAP3K12* (Table 1). The black(cases) module consisted of 580 genes with a density of 0.019, and the top three eigengenes in this module were *UBAC1*, *UBE2W* and *OAS3*. The magenta(cases) module consisted of 343 genes with a density of 0.025, and the top three eigengenes in this module were *BEST3, SLC16A2,* and *RCC1*. Purple(cases) module consisted of 303 genes with a density of 0.028, and the top three eigengenes in this module were *CPNE5*, *FMR1,* and *RBFOX2*. Salmon(cases) module consisted of only 194 genes with a density of 0.038, and the top three eigengenes in this module are *KIAA2026*, *IGF2BP3,* and *SBF2*.

View this table:
[Table 1.](http://medrxiv.org/content/early/2023/01/26/2022.12.16.22283563/T1)

Table 1. 
The top three eigengenes shared between the new midnightblue(cases) module and the preserved red(NS) module, with their gene name and gene ontology annotation. These modules were chosen as they shared the same eigengenes and were specific to TGCT cases.

### mRNA cases vs. controls preservation analysis

Preservation analysis was performed on the mRNA case/control networks to see which modules were preserved and which were not. Modules that were present in cases but absent from controls, are called new modules, and in total there were 11 new modules in this analysis. New modules are of particular interest as they could show important genes involved in TGCT development. Modules where a significant number of mRNAs are shared between paired modules, are called preserved modules, seven in total. Finally, modules where several mRNAs are shared between paired modules but not to a significant degree, are known here as non-preserved modules, and there is only one.

Blue(cases)/blue(controls) showed significant preservation, with 603 preserved genes, and a p-value < 0.001. The case module midnightblue(cases) was identified as being a new module, as well as black(cases), magenta(cases), purple(cases) and salmon(cases) modules.

### Modules of interest in Seminoma (SE) and Non-seminoma (NS)

The modules of interest for TGCT development in the SE network were the blue(SE), lightyellow(SE), and midnightblue(SE). Blue(SE) consists of 1251 genes, with a density of 0.044. The top three eigengenes for this module are *UBE2W*, *OAS3,* and *TNFAIP2*. Lightyellow(SE) consists of 168 genes, with a density of 0.06. The top three eigengenes for this module are *UVSSA*, *KLHL21,* and *FIGNL1*. The module mightnightblue(SE) consists of 231 genes, with a density of 0.097. The top three eigengenes for this module include *ZC4H2*, *ATP1B1,* and *GRID2*.

The modules of interest in the NS network were yellow(NS), which consists of 947 genes, with a density of 0.032. The top three eigengenes for this module are *CROT*, *HELB,* and *PRELID3A*. The second module of interest was red(NS), which consisted of 636 genes, with a density of 0.044. The top three hub genes for this module were *JAKMIP1*, *MAP3K12,* and *RTCB* (Table 1).

### SE vs. NS network’s preservation analysis

Preservation analysis of histology-related networks identified several modules of interest that were preserved between seminomas and non-seminomas networks. Preserved modules show commonalities between seminomas and non-seminomas. Of the 22 modules in the seminoma network and the 26 in the non-seminoma network, 20 were preserved between the two. One such module that was significantly preserved was the blue(SE)/red(NS) module, which shared 114 genes between both histology states. There were also two new modules present for seminomas and six new modules present for non-seminomas, showing possible histology-specific genes of interest. Of note is the recurring presence of this module between all case networks.

### Modules of interest in miRNA cases

The module of interest within miRNA cases was the brown(miCases) module, which consisted of 89 genes and had a density of 0.065. The top three eigenmiRs for this module were hsa-miR-30e-5p, hsa-miR-191-5p, and hsa-miR-199a-3p (Table 2).

View this table:
[Table 2.](http://medrxiv.org/content/early/2023/01/26/2022.12.16.22283563/T2)

Table 2. 
The top three eigenmiRs for the miRNA module brown(miCases) with top ten highest confidence mRNA targets IDs.

### miRNA cases vs. controls preservation analysis

Preservation analysis of miRNA networks identified four modules that were preserved between cases and controls and one new module, specific to cases. mRNA target prediction for the top three eigenmiRs within this module was performed using > 70 as the cut off value for prediction score (Table 2). A total of 1310 unique mRNA targets were identified. EigenmiRs targeting prediction showed that all three eigengenes were predicted to target *MAP3K* genes as well as related pathway genes. Hsa-miR-30e-5p was predicted to target *KRAS* and other members of the RAS oncogene family, as well the gene *MAP3K12*, with a confidence score of 78.

### Enrichment Analyses

KEGG analysis of the mRNA TGCT cases network showed enrichment of the Ubiquitin mediated proteolysis pathway in lightcyan(cases). Purple(cases) showed significant enrichment of the Notch signalling pathway. Magenta(cases), a new module, showed significant enrichment of the AMPK signalling pathway. Finally, the preserved green(cases) module also showed significant enrichment of pathways in cancer (Table 3). Furthermore, the red(SE) module showed significant enrichment of the T cell receptor signalling pathway (Table 4). Salmon(NS), a preserved module between seminomas and non-seminomas, shows significant enrichment of the non-small cell lung cancer associated pathway (Table 5).

View this table:
[Table 3.](http://medrxiv.org/content/early/2023/01/26/2022.12.16.22283563/T3)

Table 3. 
mRNA case network modules specific to mRNA cases when compared to mRNA controls. Top three eigengenes in each module listed, as well as significant enrichment analysis results for top 10 eigengenes and for all genes included in the module. Highlighted rows contain genes of interest within the top three eigengenes.

View this table:
[Table 4.](http://medrxiv.org/content/early/2023/01/26/2022.12.16.22283563/T4)

Table 4. 
mRNA network modules for seminomas, with preservation state shows mRNA seminoma vs non-seminoma analysis results. Included are the top three eigengenes in each module as well as enrichment analysis results for top 10 eigengenes and for all genes included in the module. Modules that were preserved show number of genes preserved between seminomas and non-seminomas. Highlighted rows contain genes of interest within the eigengenes.

View this table:
[Table 5.](http://medrxiv.org/content/early/2023/01/26/2022.12.16.22283563/T5)

Table 5. 
mRNA network modules for non-seminomas, with preservation state shows mRNA non-seminoma vs seminoma analysis results. Included are the top three hub genes in each module as well as enrichment analysis results for all 10 hub genes and for all genes included in the module. Modules that were preserved show number of genes preserved between seminomas and non-seminomas. Highlighted rows contain genes of interest within the eigengenes.

View this table:
[Table 6.](http://medrxiv.org/content/early/2023/01/26/2022.12.16.22283563/T6)

Table 6. 
miRNA case network modules, including the top three hub miRNAs in each module as well as enrichment analysis results for top mRNA targets for the top three hub miRNAs in the module. Preservation state shows miRNA cases vs controls preservation analysis results. Modules that were preserved show number of genes preserved between miRNA cases and miRNA controls. Highlighted rows contain miRNAs of interest within the top three eigenmiRs.

Pathway enrichment analysis was then performed using the eigenmiR’s mRNA targets. The following KEGG pathways showed significant enrichment in the brown(miCases) module: long term potentiation, axon guidance, regulation of the actin cytoskeleton, ErbB signalling pathway, T cell receptor signalling pathway, neurotrophin signalling pathway, MAPK signalling pathway, choline metabolism in cancer, FoxO signalling pathway, glioma, cGMP-PKG signalling pathway, renal cell carcinoma, cellular senescence, and B cell receptor signalling pathway. A further enrichment analysis was performed on the top ten eigenmiRs in brown(miCases), within results from the database Jensen COMPARTMENTS, showing significant enrichment of the following entries: extracellular exosome complex, lysosomal multienzyme complex, microvesicle, EmrE multidrug transporter complex, PTEN phosphatase complex, Phosphatidylinositol phosphate phosphatase complex, and apoptotic body.

### Network figures

Modelling of the network with Cytoscape elucidated interesting characteristics within the networks. The mightnightblue(SE) visualised network (Figure 2) contained 12 genes, including the top 10 hub genes from the network, excluding *ELF1*. In addition to the hub genes, the genes *TEX14, ABCB1, and DARS2* appear in the network. An enrichment analysis of these 12 genes showed results for “Carcinoma testes” in the DisGeNET database (adjusted p value < 0.1, OR = 114.02), as well as significant (p value < 0.05) enrichment in several male fertility-related disorders in MGI Mammalian Phenotype Level 4 2021, including male infertility (MP:0001925), decreased testis weight (MP:0004852), and abnormal spermatocyte morphology (MP:0006379).

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/01/26/2022.12.16.22283563/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2023/01/26/2022.12.16.22283563/F2)

Figure 2. 
Network model of modules within a,b) mRNA cases, c) seminomas, d) non-seminomas. Connections between genes show a weighted correlation in expression. The distance measurement threshold for this network was set to 0.3. Genes with coloured backgrounds are eigengenes for each module. Sequence motifs show cisbp regulatory element with highest normalised enrichment score for all genes in the module. Module was selected for visualisation due to specificity to TGCT, as well the nodes passing the adjacency threshold set in the analysis.

### Modules with shared eigengenes

Analysis of the eigengenes in modules revealed a pattern emerging in most networks related to TGCT cases. The repeated appearance of the gene MAP3K12 within the eigengenes, and the presence of MAPK signalling pathways within enrichment analysis was further investigated by isolating the modules in which this appeared and determining mRNA commonality between the three. The brown(miCases) module was included through its predicted targets with the same confidence scores used previously. The 21 genes common to all modules were the following: *ATRX, BEND7, DCAF12, FYCO1, GMEB2, GRM3, ICK, IGLON5, KLHL14, MAP3K12, MAPK4, MASTL, NRIP3, PER2, PPP1R37, PRKCE, SCN3A, SLC25A36, STYX, TTLL2, UCHL5* (Figure 3).

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/01/26/2022.12.16.22283563/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2023/01/26/2022.12.16.22283563/F3)

Figure 3. 
Common genes between each of the three modules containing, either the *MAP3K12* mRNA, or miRNAs that are known to target the *MAP3K12* gene with a confidence score of 70 in miRbase. These modules were midnightblue(cases), brown(miCases) and red(NS). Common genes between these modules were analysed due to the presence of similar eigengenes.

### Regulatory feature prediction

Analysis of the several modules’ genes using the prediction tool i-cis Target showed regulatory elements with possible links to TGCT development. For both the modules greenyellow(cases) and midnightblue(SE), the highest scoring regulatory elements from the cisbp database was oestrogen-related receptor alpha gene transcription factor (*ESRRA*), with *ESRRG* also amongst the highest scoring for midnightblue(SE) (Figure 3). The red(NS) module’s highest scoring regulatory element was associated with *CTCF*, and for the magenta(cases) module, the highest score regulatory feature was associated with *TCF3*.

## Discussion

In our previous analysis of RNA profiles from the Janus Serum Bank samples using linear models we identified RNAs which were differentially expressed in patients prior to diagnosis of TGCT. These included genes such as *TEX101,* a male fertility-related gene, and *TSIX,* X-chromosome inactivation genes (18, 41). As gene products and regulatory elements often operate in close relationships, we used network analysis framework and the tool WGCNA to identify interplay between genes and ncRNAs that could be involved in TGCT development in the years leading up to diagnosis. We combined miRNA and mRNA analyses by using miRNA targets, to further validate both types of RNA. When comparing the networks constructed from mRNA case and mRNA control data, we were able to identify several modules specific to cases that contained genes that were previously known to be associated with TGCT.

### mRNA Cases Network

Within the mRNA cases network there were several modules with eigengenes associated with both TGCT and male subfertility. The new black(cases) module’s top eigengene, *UBCA1*, has been previously seen to be downregulated in TGCT tumour samples compared to controls (42). This suggest that *UBCA1* is involved in early TGCT development and could be used as a potential biomarker. In addition, links to early defects in the male reproductive system of mice have been described with *UBE2W*, the second hub gene in this module (43). In the study, knock-out mice were used to investigate the role UBE2W plays in testis development. It was observed that the testis is vulnerable to the loss of UBE2W as seen through the high rate of male infertility after loss. Alongside this, UBE2W knock-out causes an increase in testicular vacuolation defects, pointing to its role in the maintenance of the testis structure. Defects in the male reproductive system of humans were also linked to the third hub gene, *OSA3*. Missense variants of the *OSA3* gene were found in members of a 22-family group with recurrent disorders of sexual development and hypospadias, a component of the testicular dysgenesis syndrome, often associated with TGCT (44). *UBE2W* and *OSA3* were also present in the blue(SE) module’s eigengenes, and *UBCA1*, *UBE2W,* and *OSA3* appeared in the red(NS) module’s eigengenes. This indicates that these genes may play an important role in the formation of both seminomas and non-seminomas due to their repeated appearance as eigengenes in both the overall cases, and in both histological types.

Within the new magenta(cases) module, the eigengene *RCC1* has known associations with TGCT. *RCC1* can form a fusion gene with *ABHD12B,* and this fusion, *RCC1-ABHD12B,* has been shown to be present in 9% of TGCT tissues but was not present in any control testis tissues (45). However, another study showed expression of *RCC1-ABHD12B* in 60% of GCNIS samples, in 80% of seminoma samples, and in all embryonal carcinoma cell lines (46). In this study, another fusion, *RCC1-HENMT1,* was also observed in a significant number of TGCT tissues of undifferentiated histological subtypes of TGCT, including all GCNIS samples, and all other TGCT samples, excluding on yolk sac tumour.

The MAPK signalling pathway-related gene *MAP3K12*, was observed as an eigengene for the new midnightblue(cases) module. Previously *MAP3K12* has been associated with prostate cancer, and regulation of *MAP3K12* by miRNAs could suppress prostate cancer progression (47). *MAP3K12* was also observed as one of the genes associated with activation of oncogene-induced angiogenesis in hepatocellular carcinomas (48). MAPK signalling pathways have also been labelled as one of the dominant functional pathways in spermatogenesis, alongside the AMPK pathways (49). We observed both MAPK-related genes and significant association with AMPK-related pathways in the new modules of the mRNA cases network. The involvement of these pathways in TGCT development should be investigated further as they may play an important role.

In addition, the new purple(cases) module eigengene *FMR1* has been investigated in TGCT patients and could be used to determine overall prognosis (50). It was observed that expression levels of *FMR1* were positively correlated with the clinical outcome of TGCT and therefore could act as tumour suppressors alongside *AR* and *GPC3* genes (50). The eigengene, *RBFOX2*, has been previously studied in cryptorchidism susceptibility. Paralogs of this gene are expressed in the gubernaculum; failure of this fibrous cord to elongate during development causes the birth defect commonly associated with TGCT (51, 52).

The new Salmon(cases) module contained the eigengene *IGF2BP3*, also known as *IMP3*. The IMP3 protein is an oncofoetal protein which has known roles in both embryogenesis and carcinogenesis (53). Within teratomas, IMP3 expression has been seen in 100% of mature teratoma components and in 96% of all metastatic testicular teratomas. In this study, it was also observed to be expressed in 99% of all other TGCT components (53). Therefore, it should also be investigated further for possible role in development of GCNIS and its progression to TGCT, due to its shared roles in embryogenesis and carcinogenesis.

### Histology-Related Networks

Within the seminoma network, one of the modules, the new lightyellow(SE) module, contained the eigengene *G3BP2*. Studies into TGCT have identified the cytogenetic band 4q21.1 as being a possible location for a susceptibility locus, the genes positioned within this band include *G3BP2,* whilst not fully labelled as a TGCT susceptibility gene, its location within this band should be noted (8). Our observation of this gene within seminomas could explain its identification within genome-wide associations studies (GWAS), and its role within seminoma histology should be investigated further.

The preserved midnightblue module(SE) also contained an eigengene that has previously been identified in TGCT studies (Figure 2). The *GRID2* gene is an extremely large gene (1.47Mbs within 4q22.3) that could be involved in hepatocellular carcinomas (54). Focal deletions in *GRID2* have been seen in TGCT, though it has primarily been seen within non-seminomas (55). Finally, within the non-seminoma networks, the module yellow(NS) contained the eigengene *NARS2,* which has previously been identified as a TGCT susceptibility gene through GWAS within the cytogenetic band 11q14.1 (8), and is also involved in ovarian serous cystadenocarcinomas as a possible progression associated gene (56).

Eigengenes related to TGCT, subfertility, and disorders of sexual development in the case networks are potential liquid biopsy biomarkers. With samples being taken up to 10 years before diagnosis and still exhibiting male reproductive health-related genes, it could be beneficial to look more closely into the pathways associated with these RNAs and the reason for their presence in serum. Although there is not one common pathway or mechanism between these eigengenes, we know that TGCT and cancer in general are multifaceted and complex. Even within histological subtypes of TGCT, tumour growth stems from multiple cell types with multiple mechanisms, and therefore we should expect many shared pathways for all the eigengenes related to TGCT.

TGCT analysis using network methods has been performed previously on TGCT tissue, which has allowed the identification of genes that are differentially expressed between histological subtypes (31). In particular, the transcription factor-target gene pairs, *MYC* and *SP1,* were found to be common eigengenes in both seminomas and non-seminomas, as well as the miRNA, miR-182-5p (31). *MYC* have an oncogenic association with TGCT (57), *SP1* is thought to be involved in the process of aerobic glycolysis in TGCT cells (58), and miR-182-5p and its targets are associated with cancer cell proliferation and migration (59, 60). However, within our pre-diagnostic data, we did not see *MYC* or *SP1* as eigengenes in any of the seminoma or non-seminoma modules. This could mean that these genes are not associated with progression of GCNIS to TGCT but could be involved in the maintenance of the malignancies via the previously described processes.

### Inter-network features

The mRNAs, *RTCB*, *JAKMIP1,* and *MAP3K12*, as well as the miRNA hsa-miR-30e-5p were present as eigengenes in nearly all case related modules, and hsa-miR-30e-5p is thought to potentially target both *MAP3K12* and *RTCB*. When investigating the common mRNAs shared between all three, the presence of fertility- and cancer-related genes was observed. Early spermatid production relies on ion channels for both overall fertility and normal sperm physiology, and the gene *SCNA3* is involved in this process and defects within these channels are associated with male infertility (61). The precise role of *ATRX* in human gonadal development is still not entirely clear, however, links between the gene *ATRX* and *DMRT1* have been previously noted. It is also thought that ATRX acts downstream of the genes *SPY* and *SOX9*, which are crucial for Sertoli cell development and upregulation of anti-Müllerian hormone during testicular development *in utero* (62, 63). Within mice models, *Styx* has been seen to be strongly involved in spermatid development, with ablation of *Styx* in mice causing a 1,000-fold reduction in spermatozoa production (64). These shared genes between three case-specific modules show common male fertility-related associations, indicating again the presence of possible markers for TGCT progression.

Another shared component between networks was the midnightblue(SE) module. When comparing the seminoma and non-seminoma networks together, it was noted that this midnightblue module was significantly preserved between the two histological states. Modules that were shared between seminomas and non-seminomas could offer insights into common pathways and genes of interest common to most types of TGCT. Initially the gene GRID2 had been identified as a gene of interest due to its associations with non-seminomas. However, when visualising the preserved midnightblue(SE) module in Cytoscape (Figure 2), further genes of interest were noted, such as the gene *TEX14,* which was associated with the hub genes for the module, though not an eigengene itself. *TEX14* has been identified as a susceptibility locus for TGCT (8). As shown previously, other TGCT susceptibility genes also appeared amongst the hub genes of modules in both the seminoma and non-seminoma modules, including *NARS2* in the yellow(NS) module and *G3BP2* in the lightyellow(SE) module. Both these modules are significantly preserved between both histological states. Our findings are in line with the GWAS identifying these loci regardless of histology, however, our findings propose that changes at these loci may be observable earlier than previously thought.

In the midnightblue(SE) and the greenyellow(cases) modules, we also saw the enrichment of EGRR genes during the regulatory elements analysis. *ERR* are a group of genes that all code for oestrogen-related receptor proteins, and previously they have been found to be crucial in controlling energy metabolism for both normal and for cancerous cells, particularly ERRA and ERRG (65). In some cancers, such as endometrial, colorectal, and prostate cancer, ERR expression is low, which suggests that ERRs have an overall negative effect on tumour progression (66, 67). In the formation of Leydig cells, ERRs are also crucial, as lack of ERRs has been shown to affect Leydig cell morphology and physiology in cultures (68). ERRA has also shown to mediate the adaption of metabolic systems to the tumour microenvironment and could be involved in the direct activation of Wnt11/β-catenin pathway, and this in turn leads to an increase in the capacity of cells to migrate (68). Overall, these genes are strongly involved in the molecular characteristics of Leydig cell tumours. Although this is a less common form of TGCT, the presence of ERR transcription factor enrichment in our seminoma samples could indicate a joint aetiology between the two tumour types.

Further analysis of regulatory elements identified transcription factors for the gene *CTCF* as being associated with the genes within the red(NS) module. CTCF and its paralogues have been previously identified in testicular tumour cells. It has been suggested that the co-expression of CTCF and BORIS could be responsible for epigenetic deregulation leading to cancer. Furthermore, the mechanisms surrounding CTCF/BORIS could be an essential factor leading to the immortalisation of testicular cancer cells (69).

### Strengths and Limitations

The strength of WGCNA is the reduction in the complexity of the data, allowing for a more accurate interpretation. With this approach we have detected TGCT-related genes in addition to those previously identified with linear models. The size and scope of the dataset is also advantageous and allows us to investigate pre-diagnostic TGCT patients, whilst also considering confounding factors from health surveys.

However, sample size is still a limitation with this study. Within this dataset, there is a distribution amongst different RNA types, whilst we were able to perform a network analysis for different histological types on the mRNA from our data, the scarcity of other RNA types limited sub-analysis on seminoma and non-seminoma for other types, including miRNA. The mRNA levels were derived from small RNA sequencing of long-term stored serum samples. We assume the short mRNA fragments, up to 47 nucleotides in length, represent the full-length mRNA profiles. An increase in statistical power would benefit the WGCNA analyses across histological types as well as alternate RNA types.

### Conclusion

The primary insights gained from this work include multiple pathways and genes with possible associations to TGCT initiation or growth, providing better understanding of the development of this disease. We see TGCT-related genes with high degrees of connectivity in TGCT exclusive networks, such as *TEX14*, *NARS2*, and *G3BP2*. In addition, we observe the enrichment of cancer-related mismatch repair KEGG pathway and have identified predicted transcriptional factors such as oestrogen-related receptors (ERR) in multiple TGCT modules. Enrichment analysis also showed “Carcinoma testes” and several male fertility-related disorders as significantly enriched within a seminoma module.

Furthermore, there is the potential for RNAs discussed here to be utilised as biomarkers for earlier detection. Possible future studies could involve the use of knock-out cell lines, including 3D models (70), for eigengenes discovered in the modules of our networks. Systematic investigation of these genes and their roles in cell lines would allow for validation of the altered phenotypes for each of these genes and propose a possible order of molecular events in tumorigenesis.

## Data Availability

The datasets generated for this article are not readily available because of the principles and conditions set out in articles 6 (1) (e) and 9 (2) (j) of the General Data Protection Regulation (GDPR). National legal basis as per the Regulations on population-based health surveys and ethical approval from the Norwegian Regional Committee for Medical and Health Research Ethics (REC) are also required. Requests to access the datasets should be directed to the corresponding authors.

## Data Availability

The datasets generated for this article are not readily available because of the principles and conditions set out in articles 6 (1) (e) and 9 (2) (j) of the General Data Protection Regulation (GDPR). National legal basis as per the Regulations on population-based health surveys and ethical approval from the Norwegian Regional Committee for Medical and Health Research Ethics (REC) are also required. Requests to access the datasets should be directed to the corresponding authors.

## Author Contributions

MW conceptualized the network studies. All authors designed the study. JB and MW performed the analyses of the data. All authors drafted the manuscript. All authors interpreted and discussed the results, contributed to the writing, and approved the final manuscript.

## Funding

This work was supported by internal funding of OsloMet–Oslo Metropolitan University and Cancer Registry of Norway.

## Ethics Statement

The studies involving human participants were reviewed and approved by the Norwegian Regional Committee for Medical and Health Research Ethics (REC no: 2016/1290). The patients/participants provide their written informed consent to participate in this study.

## Data Availability Statement

The datasets generated for this article are not readily available because of the principles and conditions set out in articles 6 (1) (e) and 9 (2) (j) of the General Data Protection Regulation (GDPR). National legal basis as per the Regulations on population-based health surveys and ethical approval from the Norwegian Regional Committee for Medical and Health Research Ethics (REC) are also required. Requests to access the datasets should be directed to the corresponding authors.

## Online Appendix

### Appendix 1

#### Sensitivity Analysis

Sensitivity analysis of the mRNA networks was performed through the exclusion of samples diagnosed from 8-10 years after sampling. Weighted correlation gene networks were constructed using the remaining samples diagnosed from 0-7 years after sampling with matched controls. For the sensitivity analysis of 0-7 years, the cases network consisted of 26 modules with an average module size of 412. The largest module, turquoise(cases 0-7) consisted of 1948 mRNAs and had a density of 0.012 and a heterogeneity of 0.529. Control network with 8–10-year samples excluded consisted of a total of 14 modules with an average module size of 765. The largest module here was turquoise(controls 0-7) which consisted of 2200 mRNAs with a density of 0.033 and a heterogeneity of 1.324.

Enrichment analysis of 0-7 years showed several new modules in both cases and control networks. The turquoise(cases 0-7) module finds both cancer-related and germ cell development-related pathways during enrichment analysis.

Further sensitivity analysis was performed through the exclusion of all samples diagnosed 6-10 years after sampling. Weighted gene correlation networks showed that using 0–5-year samples would allow for the construction of a network with 28 modules in total. The average module size was 382.

![Supplementary Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/01/26/2022.12.16.22283563/F4.medium.gif)

[Supplementary Figure 1.](http://medrxiv.org/content/early/2023/01/26/2022.12.16.22283563/F4)

Supplementary Figure 1. 
Power analysis for mRNA (A) and miRNA (B), showing 10 as the appropriate soft threshold for weighted gene correlation network analysis.

## Acknowledgements

We would like to acknowledge Sinan U. Umu for his contribution with the pre-processed analysis and management of the Janus Serum Bank sncRNA data.

## Footnotes

*   Small edits

*   Received December 16, 2022.
*   Revision received January 26, 2023.
*   Accepted January 26, 2023.


*   © 2023, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial 4.0 International), CC BY-NC 4.0, as described at [http://creativecommons.org/licenses/by-nc/4.0/](http://creativecommons.org/licenses/by-nc/4.0/)

## References

1.  1).Trabert, B., Chen, J., Devesa, S. S., Bray, F., & McGlynn, K. A. (2015). International patterns and trends in testicular cancer incidence, overall and by histologic subtype, 1973-2007. Andrology, 3(1), 4–12. [https://doi.org/10.1111/andr.293](https://doi.org/10.1111/andr.293)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/andr.293&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25331326&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F01%2F26%2F2022.12.16.22283563.atom) 

2.  2).Znaor, A., Skakkebaek, N. E., Rajpert-De Meyts, E., Kuliš, T., Laversanne, M., Gurney, J., Sarfati, D., McGlynn, K. A., & Bray, F. (2022). Global patterns in testicular cancer incidence and mortality in 2020. International Journal of Cancer, 151(5), 692–698. [https://doi.org/10.1002/ijc.33999](https://doi.org/10.1002/ijc.33999)
    
    

3.  3).Cancer Registry of Norway. Cancer in Norway 2021 - Cancer incidence, mortality, survival and prevalence in Norway. Oslo: Cancer Registry of Norway, 2022.
    
    

4.  4).Skakkebaek, N.E., Berthelsen, J.G., Giwercman, A., Muller, J. (1987). Carcinoma-in-situ of the testis: possible origin from gonocytes and precursor of all types of germ cell tumours except spermatocytoma. International Journal of Andrology, 10, 19–28. [https://doi.org/10.1111/j.1365-2605.1987.tb00161.x](https://doi.org/10.1111/j.1365-2605.1987.tb00161.x)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1365-2605.1987.tb00161.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=3034791&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F01%2F26%2F2022.12.16.22283563.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1987G683400003&link_type=ISI) 

5.  5).Litchfield, K., Mitchell, J. S., Shipley, J., Huddart, R., Rajpert-De Meyts, E., Skakkebæk, N. E., Houlston, R. S., & Turnbull, C. (2015). Polygenic susceptibility to testicular cancer: implications for personalised health care. British Journal of Cancer, 113(10), 1512–1518. [https://doi.org/10.1038/bjc.2015.334](https://doi.org/10.1038/bjc.2015.334)
    
    

6.  6).Wang, Z., McGlynn, K. A., Rajpert-De Meyts, E., Bishop, D. T., Chung, C. C., Dalgaard, M. D., Greene, M. H., Gupta, R., Grotmol, T., Haugen, T. B., Karlsson, R., Litchfield, K., Mitra, N., Nielsen, K., Pyle, L. C., Schwartz, S. M., Thorsson, V., Vardhanabhuti, S., Wiklund, F., … Nathanson, K. L. (2017). Meta-analysis of five genome-wide association studies identifies multiple new loci associated with testicular germ cell tumor. Nature Genetics, 49(7), 1141–1147. [https://doi.org/10.1038/ng.3879](https://doi.org/10.1038/ng.3879)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3879&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28604732&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F01%2F26%2F2022.12.16.22283563.atom) 

7.  7).Pyle, L. C., & Nathanson, K. L. (2016). Genetic changes associated with testicular cancer susceptibility. Seminars in Oncology, 43(5), 575–581. [https://doi.org/10.1053/j.seminoncol.2016.08.004](https://doi.org/10.1053/j.seminoncol.2016.08.004)
    
    

8.  8).Pluta, J., Pyle, L. C., Nead, K. T., Wilf, R., Li, M., Mitra, N., … Consortium, T. T. C. (2021). Identification of 22 susceptibility loci associated with testicular germ cell tumors. Nature Communications, 12(1), 4487. [http://doi.org/10.1038/s41467-021-24334-y](http://doi.org/10.1038/s41467-021-24334-y)
    
    

9.  9).Skakkebæk, N. E., Rajpert-De Meyts, E., & Main, K. M. (2001). Testicular dysgenesis syndrome: an increasingly common developmental disorder with environmental aspects: Opinion. Human Reproduction, 16(5), 972–978. [https://doi.org/10.1093/humrep/16.5.972](https://doi.org/10.1093/humrep/16.5.972)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/humrep/16.5.972&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11331648&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F01%2F26%2F2022.12.16.22283563.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000168646500027&link_type=ISI) 

10. 10).Dieckmann, K.-P., & Pichlmeier, U. (2004). Clinical epidemiology of testicular germ cell tumors. World Journal of Urology, 22(1), 2–14. [https://doi.org/10.1007/s00345-004-0398-8](https://doi.org/10.1007/s00345-004-0398-8)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00345-004-0398-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15034740&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F01%2F26%2F2022.12.16.22283563.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000222089800002&link_type=ISI) 

11. 11).Kvammen, Ø., Myklebust, T. Å., Solberg, A., Møller, B., Klepp, O. H., Fosså, S. D., & Tandstad, T. (2016). Long-term Relative Survival after Diagnosis of Testicular Germ Cell Tumor. Cancer Epidemiology, Biomarkers & Prevention, 25(5), 773–779. [https://doi.org/10.1158/1055-9965.EPI-15-1153](https://doi.org/10.1158/1055-9965.EPI-15-1153)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiY2VicCI7czo1OiJyZXNpZCI7czo4OiIyNS81Lzc3MyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzAxLzI2LzIwMjIuMTIuMTYuMjIyODM1NjMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

12. 12).Bucher-Johannessen, C., Page, C. M., Haugen, T. B., Wojewodzic, M. W., Fosså, S. D., Grotmol, T., Haugnes, H. S., & Rounge, T. B. (2019). Cisplatin treatment of testicular cancer patients introduces long-term changes in the epigenome. Clinical Epigenetics, 11(1), 179. [https://doi.org/10.1186/s13148-019-0764-4](https://doi.org/10.1186/s13148-019-0764-4)
    
    

13. 13).Haugnes, H. S., Negaard, H. F., Jensvoll, H., Wilsgaard, T., Tandstad, T., & Solberg, A. (2021). Thromboembolic Events During Treatment with Cisplatin-based Chemotherapy in Metastatic Testicular Germ-cell Cancer 2000-2014: A Population-based Cohort Study. European Urology Open Science, 32, 19–27. [https://doi.org/10.1016/j.euros.2021.07.007](https://doi.org/10.1016/j.euros.2021.07.007)
    
    

14. 14).Raphael, M. J., Wei, X., Karim, S., Robinson, A. G., Bedard, P. L., & Booth, C. M. (2019). Neurotoxicity Among Survivors of Testicular Cancer: A Population-based Study. Clinical Oncology (Royal College of Radiologists (Great Britain)), 31(9), 653– 658. [https://doi.org/10.1016/j.clon.2019.04.008](https://doi.org/10.1016/j.clon.2019.04.008)
    
    

15. 15).Hellesnes, R., Kvammen, Ø., Myklebust, T. Å., Bremnes, R. M., Karlsdottir, Á., Negaard, H. F. S., Tandstad, T., Wilsgaard, T., Fosså, S. D., & Haugnes, H. S. (2020). Continuing increased risk of second cancer in long-term testicular cancer survivors after treatment in the cisplatin era. International Journal of Cancer, 147(1), 21–32. [https://doi.org/10.1002/ijc.32704](https://doi.org/10.1002/ijc.32704)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ijc.32704&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31597192&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F01%2F26%2F2022.12.16.22283563.atom) 

16. 16).Hjerkind, K. V., Gislefoss, R. E., Tretli, S., Nystad, W., Bjørge, T., Engeland, A., … Langseth, H. (2017). Cohort Profile Update: The Janus Serum Bank Cohort in Norway. International Journal of Epidemiology, 46(4), 1101–1102f. [http://doi.org/10.1093/ije/dyw302](http://doi.org/10.1093/ije/dyw302)
    
    

17. 17).Langseth, H., Gislefoss, R. E., Martinsen, J. I., Dillner, J., & Ursin, G. (2017). Cohort Profile: The Janus Serum Bank Cohort in Norway. International Journal of Epidemiology, 46(2), 403–404g. [http://doi.org/10.1093/ije/dyw027](http://doi.org/10.1093/ije/dyw027)
    
    

18. 18).Burton, J., Umu, S. U., Langseth, H., Grotmol, T., Grimsrud, T. K., Haugen, T. B., & Rounge, T. B. (2020). Serum RNA Profiling in the 10-Years Period Prior to Diagnosis of Testicular Germ Cell Tumor. Frontiers in Oncology, 10, 2207. [http://doi.org/10.3389/fonc.2020.574977](http://doi.org/10.3389/fonc.2020.574977)
    
    

19. 19).Umu, S. U., Langseth, H., Zuber, V., Helland, Å., Lyle, R., & Rounge, T. B. (2022). Serum RNAs can predict lung cancer up to 10 years prior to diagnosis. eLife, 11, e71035. [http://doi.org/10.7554/eLife.71035](http://doi.org/10.7554/eLife.71035)
    
    

20. 20).Langfelder, P., & Horvath, S. (2008). WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics, 1, 559. [https://bmcbioinformatics.biomedcentral.com/articles/10.1186/1471-2105-9-559](https://bmcbioinformatics.biomedcentral.com/articles/10.1186/1471-2105-9-559)
    
    

21. 21).Yuan, L., Qian, G., Chen, L., Wu, C.-L., Dan, H. C., Xiao, Y., & Wang, X. (2018). Co-expression Network Analysis of Biomarkers for Adrenocortical Carcinoma. Frontiers in Genetics, 9, 328. [https://doi.org/10.3389/fgene.2018.00328](https://doi.org/10.3389/fgene.2018.00328)
    
    

22. 22).Robinson, M. D., McCarthy, D. J., & Smyth, G. K. (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics (Oxford, England), 26(1), 139–140. [https://doi.org/10.1093/bioinformatics/btp616](https://doi.org/10.1093/bioinformatics/btp616)
    
    

23. 23).Love, M. I., Huber, W., & Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biology, 15(12), 550. [https://doi.org/10.1186/s13059-014-0550-8](https://doi.org/10.1186/s13059-014-0550-8)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13059-014-0550-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25516281&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F01%2F26%2F2022.12.16.22283563.atom) 

24. 24).Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W., & Smyth, G. K. (2015). limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Research, 43(7), e47. [https://doi.org/10.1093/nar/gkv007](https://doi.org/10.1093/nar/gkv007)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkv007&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25605792&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F01%2F26%2F2022.12.16.22283563.atom) 

25. 25).Liao, Y., Wang, Y., Cheng, M., Huang, C., & Fan, X. (2020). Weighted Gene Coexpression Network Analysis of Features That Control Cancer Stem Cells Reveals Prognostic Biomarkers in Lung Adenocarcinoma. Frontiers in Genetics, 11. [https://www.frontiersin.org/articles/10.3389/fgene.2020.00311](https://www.frontiersin.org/articles/10.3389/fgene.2020.00311)
    
    

26. 26).Di, Y., Chen, D., Yu, W., & Yan, L. (2019). Bladder cancer stage-associated hub genes revealed by WGCNA co-expression network analysis. Hereditas, 156, 7. [https://doi.org/10.1186/s41065-019-0083-y](https://doi.org/10.1186/s41065-019-0083-y)
    
    

27. 27).Shi, G., Shen, Z., Liu, Y., & Yin, W. (2020). Identifying Biomarkers to Predict the Progression and Prognosis of Breast Cancer by Weighted Gene Co-expression Network Analysis. Frontiers in Genetics, 11. [https://www.frontiersin.org/articles/10.3389/fgene.2020.597888](https://www.frontiersin.org/articles/10.3389/fgene.2020.597888)
    
    

28. 28).Pu, M., Chen, J., Tao, Z., Miao, L., Qi, X., Wang, Y., & Ren, J. (2019). Regulatory network of miRNA on its target: coordination between transcriptional and post-transcriptional regulation of gene expression. Cellular and Molecular Life Sciences : CMLS, 76(3), 441–451. [https://doi.org/10.1007/s00018-018-2940-7](https://doi.org/10.1007/s00018-018-2940-7)
    
    

29. 29).Boguslawska, J., Kryst, P., Poletajew, S., & Piekielko-Witkowska, A. (2019). TGF-β and microRNA Interplay in Genitourinary Cancers. Cells, 8(12). [https://doi.org/10.3390/cells8121619](https://doi.org/10.3390/cells8121619)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/cells8121584&link_type=DOI) 

30. 30).Chen, Y., Liao, R., Yao, Y., Wang, Q., & Fu, L. (2022). Machine learning to identify immune-related biomarkers of rheumatoid arthritis based on WGCNA network. Clinical Rheumatology, 41(4), 1057–1068. [https://doi.org/10.1007/s10067-021-05960-9](https://doi.org/10.1007/s10067-021-05960-9)
    
    

31. 31).Mallik, S., Qin, G., Jia, P., & Zhao, Z. (2020). Molecular signatures identified by integrating gene expression and methylation in non-seminoma and seminoma of testicular germ cell tumours. Epigenetics, 16(2), 162–176. [https://doi.org/10.1080/15592294.2020.1790108](https://doi.org/10.1080/15592294.2020.1790108)
    
    

32. 32).Zhao, Y., Xu, Z., Wang, N., & Wang, S. (2016). Regulatory network of microRNAs and genes in testicular cancer. Oncology Letters, 12(5), 3640–3646. [http://doi.org/10.3892/ol.2016.5043](http://doi.org/10.3892/ol.2016.5043)
    
    

33. 33).Langseth, H., Umu, S. U., Bucher-Johannessen, C., Babigumira, R., Leithaug, M., Lauritzen, M., … Rounge, T. B. (2021). Data Resource Profile: thousands of circulating RNA profiles of pre-clinical samples from the Janus Serum Bank Cohort. (Pre-print) medRxiv, 2021.01.22.20243154. [http://doi.org/10.1101/2021.01.22.20243154](http://doi.org/10.1101/2021.01.22.20243154)
    
    

34. 34).Umu, S. U., Langseth, H., Bucher-Johannessen, C., Fromm, B., Keller, A., Meese, E., … Rounge, T. B. (2018). A comprehensive profile of circulating RNAs in human serum. RNA Biology, 15(2), 242–250. [http://doi.org/10.1080/15476286.2017.1403003](http://doi.org/10.1080/15476286.2017.1403003)
    
    

35. 35).Johnson, W. E., Li, C., & Rabinovic, A. (2007). Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics (Oxford, England), 8(1), 118–127. [http://doi.org/10.1093/biostatistics/kxj037](http://doi.org/10.1093/biostatistics/kxj037)
    
    

36. 36).Langfelder, P., Luo, R., Oldham, MC., Horvath S. (2011). Is My Network Module Preserved and Reproducible? PLoS Computational Biology 7(1): e1001057. [https://doi.org/10.1371/journal.pcbi.1001057](https://doi.org/10.1371/journal.pcbi.1001057)
    
    

37. 37).Frisch, D., Becker, D., & Wojewodzic, M. W. (2020). Dissecting the Transcriptomic Basis of Phenotypic Evolution in an Aquatic Keystone Grazer. Molecular Biology and Evolution, 37(2), 475–487. [http://doi.org/10.1093/molbev/msz234](http://doi.org/10.1093/molbev/msz234)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/molbev/msz234&link_type=DOI) 

38. 38).Chen, E. Y., Tan, C. M., Kou, Y., Duan, Q., Wang, Z., Meirelles, G. V., Clark, N. R., & Ma’ayan, A. (2013). Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinformatics, 14, 128. [https://doi.org/10.1186/1471-2105-14-128](https://doi.org/10.1186/1471-2105-14-128)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2105-14-128&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23586463&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F01%2F26%2F2022.12.16.22283563.atom) 

39. 39).Kuleshov, M. V, Jones, M. R., Rouillard, A. D., Fernandez, N. F., Duan, Q., Wang, Z., Koplev, S., Jenkins, S. L., Jagodnik, K. M., Lachmann, A., McDermott, M. G., Monteiro, C. D., Gundersen, G. W., & Ma’ayan, A. (2016). Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Research, 44(W1), W90–7. [https://doi.org/10.1093/nar/gkw377](https://doi.org/10.1093/nar/gkw377)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkw377&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27141961&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F01%2F26%2F2022.12.16.22283563.atom) 

40. 40).Xie, Z., Bailey, A., Kuleshov, M. V, Clarke, D. J. B., Evangelista, J. E., Jenkins, S. L., Lachmann, A., Wojciechowicz, M. L., Kropiwnicki, E., Jagodnik, K. M., Jeon, M., & Ma’ayan, A. (2021). Gene Set Knowledge Discovery with Enrichr. Current Protocols, 1(3), e90. https://doi.org/[https://doi.org/10.1002/cpz1.90](https://doi.org/10.1002/cpz1.90)
    
    

41. 41).Burton, J., Wojewodzic, M. W., Rounge, T. B., & Haugen, T. B. (2022). A Role of the TEX101 Interactome in the Common Aetiology Behind Male Subfertility and Testicular Germ Cell Tumor. Frontiers in Oncology. Retrieved from [https://www.frontiersin.org/article/10.3389/fonc.2022.892043](https://www.frontiersin.org/article/10.3389/fonc.2022.892043)
    
    

42. 42).Zhang, C., Zhang, W., Cui, H., Zhang, B., Miao, P., Yang, Q., Bai, M., Jiao, H., & Chang, D. (2022). Role of Hub Genes in the Occurrence and Development of Testicular Cancer Based on Bioinformatics. International Journal of General Medicine, 15, 645–660. [https://doi.org/10.2147/IJGM.S342611](https://doi.org/10.2147/IJGM.S342611)
    
    

43. 43).Wang, B., Merillat, S. A., Vincent, M., Huber, A. K., Basrur, V., Mangelberger, D., Zeng, L., Elenitoba-Johnson, K., Miller, R. A., Irani, D. N., Dlugosz, A. A., Schnell, S., Scaglione, K. M., & Paulson, H. L. (2016). Loss of the Ubiquitin-conjugating Enzyme UBE2W Results in Susceptibility to Early Postnatal Lethality and Defects in Skin, Immune, and Male Reproductive Systems. The Journal of Biological Chemistry, 291(6), 3030–3042. [https://doi.org/10.1074/jbc.M115.676601](https://doi.org/10.1074/jbc.M115.676601)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEwOiIyOTEvNi8zMDMwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDEvMjYvMjAyMi4xMi4xNi4yMjI4MzU2My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

44. 44).Tsai, C.-L., Tsai, C.-N., Lee, Y.-S., Wang, H.-S., Lee, L.-Y., Lin, C.-Y., Yang, S. Y., & Chao, A. (2020). Genetic analysis of a Taiwanese family identifies a DMRT3-OAS3 interaction that is involved in human sexual differentiation through the regulation of ESR1 expression. Fertility and Sterility, 114(1), 133–143. https://doi.org/[https://doi.org/10.1016/j.fertnstert.2020.03.008](https://doi.org/10.1016/j.fertnstert.2020.03.008)
    
    

45. 45).Zhao, S., Hoff, A. M., & Skotheim, R. I. (2020). ScaR-a tool for sensitive detection of known fusion transcripts: establishing prevalence of fusions in testicular germ cell tumors. NAR Genomics and Bioinformatics, 2(1), lqz025. [https://doi.org/10.1093/nargab/lqz025](https://doi.org/10.1093/nargab/lqz025)
    
    

46. 46).Hoff, A. M., Alagaratnam, S., Zhao, S., Bruun, J., Andrews, P. W., Lothe, R. A., & Skotheim, R. I. (2016). Identification of Novel Fusion Genes in Testicular Germ Cell Tumors. Cancer Research, 76(1), 108–116. [https://doi.org/10.1158/0008-5472.CAN-15-1790](https://doi.org/10.1158/0008-5472.CAN-15-1790)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiY2FucmVzIjtzOjU6InJlc2lkIjtzOjg6Ijc2LzEvMTA4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDEvMjYvMjAyMi4xMi4xNi4yMjI4MzU2My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

47. 47).Yu, J., Feng, Y., Wang, Y., & An, R. (2018). Aryl hydrocarbon receptor enhances the expression of miR-150-5p to suppress in prostate cancer progression by regulating MAP3K12. Archives of Biochemistry and Biophysics, 654, 47–54. https://doi.org/[https://doi.org/10.1016/j.abb.2018.07.010](https://doi.org/10.1016/j.abb.2018.07.010)
    
    

48. 48).Sun, L., Xi, S., Zhou, Z., Zhang, F., Hu, P., Cui, Y., Wu, S., Wang, Y., Wu, S., Wang, Y., Du, Y., Zheng, J., Yang, H., Chen, M., Yan, Q., Yu, D., Shi, C., Zhang, Y., Xie, D., … Li, Y. (2022). Elevated expression of RIT1 hyperactivates RAS/MAPK signal and sensitizes hepatocellular carcinoma to combined treatment with sorafenib and AKT inhibitor. Oncogene, 41(5), 732–744. [https://doi.org/10.1038/s41388-021-02130-8](https://doi.org/10.1038/s41388-021-02130-8)
    
    

49. 49).Ni, F.-D., Hao, S.-L., & Yang, W.-X. (2019). Multiple signaling pathways in Sertoli cells: recent findings in spermatogenesis. Cell Death & Disease, 10(8), 541. [https://doi.org/10.1038/s41419-019-1782-z](https://doi.org/10.1038/s41419-019-1782-z)
    
    

50. 50).Cong, R., Ji, C., Zhang, J., Zhang, Q., Zhou, X., Yao, L., Luan, J., Meng, X., & Song, N. (2021). m6A RNA methylation regulators play an important role in the prognosis of patients with testicular germ cell tumor. Translational Andrology and Urology, 10(2), 662–679. [https://doi.org/10.21037/tau-20-963](https://doi.org/10.21037/tau-20-963)
    
    

51. 51).McGlynn, K. A. (2001). Environmental and host factors in testicular germ cell tumors. Cancer Investigation, 19(8), 842–853. [https://doi.org/10.1081/cnv-100107746](https://doi.org/10.1081/cnv-100107746)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1081/CNV-100107746&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11768038&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F01%2F26%2F2022.12.16.22283563.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000172523000011&link_type=ISI) 

52. 52).Wang, Y., Gray, D. R., Robbins, A. K., Crowgey, E. L., Chanock, S. J., Greene, M. H., McGlynn, K. A., Nathanson, K., Turnbull, C., Wang, Z., Devoto, M., & Barthold, J. S. (2018). Subphenotype meta-analysis of testicular cancer genome-wide association study data suggests a role for RBFOX family genes in cryptorchidism susceptibility. Human Reproduction (Oxford, England), 33(5), 967–977. [https://doi.org/10.1093/humrep/dey066](https://doi.org/10.1093/humrep/dey066)
    
    

53. 53).Goodman, S., Zhang, L., Cheng, L., & Jiang, Z. (2014). Differential expression of IMP3 between male and female mature teratomas—immunohistochemical evidence of malignant nature. Histopathology, 65(4), 483–489. https://doi.org/[https://doi.org/10.1111/his.12409](https://doi.org/10.1111/his.12409)
    
    

54. 54).Smith, D. I., Zhu, Y., McAvoy, S., & Kuhn, R. (2006). Common fragile sites, extremely large genes, neural development and cancer. Cancer Letters, 232(1), 48–57. [https://doi.org/10.1016/j.canlet.2005.06.049](https://doi.org/10.1016/j.canlet.2005.06.049)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.canlet.2005.06.049&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16221525&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F01%2F26%2F2022.12.16.22283563.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000235251600006&link_type=ISI) 

55. 55).Shen, H., Shih, J., Hollern, D. P., Wang, L., Bowlby, R., Tickoo, S. K., Thorsson, V., Mungall, A. J., Newton, Y., Hegde, A. M., Armenia, J., Sánchez-Vega, F., Pluta, J., Pyle, L. C., Mehra, R., Reuter, V. E., Godoy, G., Jones, J., Shelley, C. S., … Hoadley, K. A. (2018). Integrated Molecular Characterization of Testicular Germ Cell Tumors. Cell Reports, 23(11), 3392–3406. https://doi.org/[https://doi.org/10.1016/j.celrep.2018.05.039](https://doi.org/10.1016/j.celrep.2018.05.039)
    
    

56. 56).Zhang, X., & Wang, Y. (2019). Identification of hub genes and key pathways associated with the progression of gynecological cancer. Oncology Letters, 18(6), 6516–6524. [https://doi.org/10.3892/ol.2019.11004](https://doi.org/10.3892/ol.2019.11004)
    
    

57. 57).Ventelä, S., Mäkelä, J.-A., Sears, R. C., Toppari, J., & Westermarck, J. (2016). MYC is not detected in highly proliferating normal spermatogonia but is coupled with CIP2A in testicular cancers. Matters, 2016. [https://doi.org/10.19185/matters.201602000040](https://doi.org/10.19185/matters.201602000040)
    
    

58. 58).Zhou, S., Min, Z., Sun, K., Qu, S., Zhou, J., Duan, H., Liu, H., Liu, X., Gong, Z., & Li, D. (2018). miR-199a-3p/Sp1/LDHA axis controls aerobic glycolysis in testicular tumor cells. Int J Mol Med, 42(4), 2163–2174. [https://doi.org/10.3892/ijmm.2018.3771](https://doi.org/10.3892/ijmm.2018.3771)
    
    

59. 59).Wang, F., Wu, D., Xu, Z., Chen, J., Zhang, J., Li, X., Chen, S., He, F., Xu, J., Su, L., Luo, D., Zhang, S., & Wang, W. (2019). miR-182-5p affects human bladder cancer cell proliferation, migration and invasion through regulating Cofilin 1. Cancer Cell International, 19, 42. [https://doi.org/10.1186/s12935-019-0758-5](https://doi.org/10.1186/s12935-019-0758-5)
    
    

60. 60).Zhao, Y. S., Yang, W. C., Xin, H. W., Han, J. X., & Ma, S. G. (2019). MiR-182-5p Knockdown Targeting PTEN Inhibits Cell Proliferation and Invasion of Breast Cancer Cells. Yonsei Medical Journal, 60(2), 148–157. [https://doi.org/10.3349/ymj.2019.60.2.148](https://doi.org/10.3349/ymj.2019.60.2.148)
    
    

61. 61).Shukla, K. K., Mahdi, A. A., & Rajender, S. (2012). Ion channels in sperm physiology and male fertility and infertility. Journal of Andrology, 33(5), 777–788. [https://doi.org/10.2164/jandrol.111.015552](https://doi.org/10.2164/jandrol.111.015552)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2164/jandrol.111.015552&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22441763&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F01%2F26%2F2022.12.16.22283563.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000308782500002&link_type=ISI) 

62. 62).Pask, A., Renfree, M. B., & Marshall Graves, J. A. (2000). The human sex-reversing ATRX gene has a homologue on the marsupial Y chromosome, ATRY: implications for the evolution of mammalian sex determination. Proceedings of the National Academy of Sciences of the United States of America, 97(24), 13198–13202. [https://doi.org/10.1073/pnas.230424497](https://doi.org/10.1073/pnas.230424497)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiOTcvMjQvMTMxOTgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wMS8yNi8yMDIyLjEyLjE2LjIyMjgzNTYzLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

63. 63).Huyhn, K., Renfree, M. B., Graves, J. A., & Pask, A. J. (2011). ATRX has a critical and conserved role in mammalian sexual differentiation. BMC Developmental Biology, 11, 39. [https://doi.org/10.1186/1471-213X-11-39](https://doi.org/10.1186/1471-213X-11-39)
    
    

64. 64).Wishart, M. J., & Dixon, J. E. (2002). The archetype STYX/dead-phosphatase complexes with a spermatid mRNA-binding protein and is essential for normal sperm production. Proceedings of the National Academy of Sciences of the United States of America, 99(4), 2112–2117. [https://doi.org/10.1073/pnas.251686198](https://doi.org/10.1073/pnas.251686198)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czo5OiI5OS80LzIxMTIiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wMS8yNi8yMDIyLjEyLjE2LjIyMjgzNTYzLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

65. 65).Vernier, M., Dufour, C. R., McGuirk, S., Scholtes, C., Li, X., Bourmeau, G., … Giguère, V. (2020). Estrogen-related receptors are targetable ROS sensors. Genes & Development, 34(7–8), 544–559. [http://doi.org/10.1101/gad.330746.119](http://doi.org/10.1101/gad.330746.119)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZ2VuZXNkZXYiO3M6NToicmVzaWQiO3M6MTA6IjM0LzctOC81NDQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wMS8yNi8yMDIyLjEyLjE2LjIyMjgzNTYzLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

66. 66).Bianco, S., Sailland, J., & Vanacker, J.-M. (2012). ERRs and cancers: Effects on metabolism and on proliferation and migration capacities. The Journal of Steroid Biochemistry and Molecular Biology, 130(3), 180–185. https://doi.org/[https://doi.org/10.1016/j.jsbmb.2011.03.014](https://doi.org/10.1016/j.jsbmb.2011.03.014)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jsbmb.2011.03.014&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21414406&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F01%2F26%2F2022.12.16.22283563.atom) 

67. 67).Misawa, A., & Inoue, S. (2015). Estrogen-Related Receptors in Breast Cancer and Prostate Cancer. Frontiers in Endocrinology, 6, 83. [https://doi.org/10.3389/fendo.2015.00083](https://doi.org/10.3389/fendo.2015.00083)
    
    

68. 68).Kotula-Balak, M., Milon, A., Pawlicki, P., Opydo-Chanek, M., Pacwa, A., Lesniak, K., Sekula, M., Zarzycka, M., Bubka, M., Tworzydlo, W., Bilinska, B., & Hejmej, A. (2018). Insights into the role of estrogen-related receptors α, β and γ in tumor Leydig cells. Tissue and Cell, 52, 78–91. https://doi.org/[https://doi.org/10.1016/j.tice.2018.04.003](https://doi.org/10.1016/j.tice.2018.04.003)
    
    

69. 69).Renaud S, Loukinov D, Alberti L, Vostrov A, Kwon YW, Bosman FT, Lobanenkov V, Benhattar J. (2011). BORIS/CTCFL-mediated transcriptional regulation of the hTERT telomerase gene in testicular and ovarian tumor cells. Nucleic Acids Research. Feb;39(3):862–73. [https://doi.org/10.1093/nar/gkq827](https://doi.org/10.1093/nar/gkq827).
    
    

70. 70).Power, R. J., Hearn, J., Gillis, C. J., Harvey, D., French, C., & Organ, M. (2021). Development of a 3D-printed testicular cancer model for testicular examination education. Canadian Urological Association Journal, 15(4), E221–E226. [https://doi.org/10.5489/cuaj.6675](https://doi.org/10.5489/cuaj.6675)