<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//TaxonX//DTD Taxonomic Treatment Publishing DTD v0 20100105//EN" "https://foliamedica.bg/nlm/tax-treatment-NS0.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:tp="http://www.plazi.org/taxpub" article-type="research-article" dtd-version="3.0" xml:lang="en">
  <front>
    <journal-meta>
      <journal-id journal-id-type="publisher-id">87</journal-id>
      <journal-id journal-id-type="index">urn:lsid:arphahub.com:pub:A116C711-4C18-5A38-8F1E-5E97753A8A64</journal-id>
      <journal-title-group>
        <journal-title xml:lang="en">Folia Medica</journal-title>
        <abbrev-journal-title xml:lang="en">FM</abbrev-journal-title>
      </journal-title-group>
      <issn pub-type="ppub">0204-8043</issn>
      <issn pub-type="epub">1314-2143</issn>
      <publisher>
        <publisher-name>Plovdiv Medical University</publisher-name>
      </publisher>
    </journal-meta>
    <article-meta>
      <article-id pub-id-type="doi">10.3897/folmed.67.e168906</article-id>
      <article-id pub-id-type="publisher-id">168906</article-id>
      <article-categories>
        <subj-group subj-group-type="heading">
          <subject>Research Article</subject>
        </subj-group>
        <subj-group subj-group-type="scientific_subject">
          <subject>Diagnostic medicine</subject>
          <subject>General Pathology</subject>
          <subject>Molecular biology</subject>
          <subject>Pediatrics &amp; Genetic diseases</subject>
        </subj-group>
      </article-categories>
      <title-group>
        <article-title>﻿Identifying potential novel biomarkers for varicocele: A bioinformatics approach to genomics analysis</article-title>
      </title-group>
      <contrib-group content-type="authors">
        <contrib contrib-type="author" corresp="yes">
          <name name-style="western">
            <surname>Farid</surname>
            <given-names>Muhammad</given-names>
          </name>
          <email xlink:type="simple">muhammad2100034023@webmail.uad.ac.id</email>
          <uri content-type="orcid">https://orcid.org/0009-0002-5301-0854</uri>
          <xref ref-type="aff" rid="A1">1</xref>
        </contrib>
        <contrib contrib-type="author" corresp="no">
          <name name-style="western">
            <surname>Aqilanisa</surname>
            <given-names>Hanin Fitri</given-names>
          </name>
          <xref ref-type="aff" rid="A1">1</xref>
        </contrib>
      </contrib-group>
      <aff id="A1">
        <label>1</label>
        <addr-line content-type="verbatim">Faculty of Medicine, Ahmad Dahlan University, Yogyakarta, Indonesia</addr-line>
        <institution>Faculty of Medicine, Ahmad Dahlan University</institution>
        <addr-line content-type="city">Yogyakarta</addr-line>
        <country>Indonesia</country>
      </aff>
      <author-notes>
        <fn fn-type="corresp">
          <p>Corresponding author: Muhammad Farid, Faculty of Medicine, Ahmad Dahlan University, Yogyakarta, Indonesia; Email: <email xlink:type="simple">muhammad2100034023@webmail.uad.ac.id</email></p>
        </fn>
      </author-notes>
      <pub-date pub-type="collection">
        <year>2025</year>
      </pub-date>
      <pub-date pub-type="epub">
        <day>18</day>
        <month>12</month>
        <year>2025</year>
      </pub-date>
      <volume>67</volume>
      <issue>6</issue>
      <elocation-id>e168906</elocation-id>
      <uri content-type="arpha" xlink:href="http://openbiodiv.net/92DF6751-6692-5E1E-959F-933595EE1CAD">92DF6751-6692-5E1E-959F-933595EE1CAD</uri>
      <history>
        <date date-type="received">
          <day>15</day>
          <month>08</month>
          <year>2025</year>
        </date>
        <date date-type="accepted">
          <day>02</day>
          <month>09</month>
          <year>2025</year>
        </date>
      </history>
      <permissions>
        <copyright-statement>Muhammad Farid, Hanin Fitri Aqilanisa</copyright-statement>
        <license license-type="creative-commons-attribution" xlink:href="http://creativecommons.org/licenses/by/4.0/" xlink:type="simple">
          <license-p>This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</license-p>
        </license>
      </permissions>
      <abstract>
        <p>﻿<bold>Abstract</bold></p>
        <p><bold>Introduction</bold>: Varicocele, characterized by the enlargement of scrotal veins, is a common contributor to male infertility, but its genetic underpinnings remain largely unknown.</p>
        <p><bold>Aim</bold>: The goal of this study is to identify potential biomarkers associated with varicocele in order to better understand its molecular mechanisms.</p>
        <p><bold>Materials and methods</bold>: Using the three primary databases, NCBI, DisGeNET, and OpenTarget, we analyzed gene variants and found 79 pertinent genes associated with varicocele. Protein-protein interaction analysis was performed using STRING and visualized with Cytoscape. Molecular Complex Detection (<abbrev xlink:title="Molecular Complex Detection" id="ABBRID0EPD">MCODE</abbrev>) and CytoHubba tools helped identify significant protein clusters.</p>
        <p><bold>Results</bold>: The gene ontology analysis shows that there are 79 proteins involved in the inflammatory process, the regulation of gene expression, and cellular components that play a role in oxidative stress and angiogenesis. Our results revealed three key biomarkers: Interleukin-1 beta (<abbrev xlink:title="Interleukin-1 beta" id="ABBRID0EXD">IL1B</abbrev>), B-cell lymphoma 2 (<abbrev xlink:title="B-cell lymphoma 2" id="ABBRID0E2D">BCL2</abbrev>), and matrix metalloproteinase-9 (<abbrev xlink:title="matrix metalloproteinase-9" id="ABBRID0E6D">MMP-9</abbrev>). These proteins are involved in critical processes, such as inflammation, oxidative stress, angiogenesis, and vascular damage, that are central to the pathophysiology of varicocele.</p>
        <p><bold>Conclusion</bold>: The identification of <abbrev xlink:title="Interleukin-1 beta" id="ABBRID0EHE">IL1B</abbrev>, <abbrev xlink:title="B-cell lymphoma 2" id="ABBRID0ELE">BCL2</abbrev>, and <abbrev xlink:title="matrix metalloproteinase-9" id="ABBRID0EPE">MMP-9</abbrev> offers new insights into varicocele’s molecular mechanisms and suggests potential targets for diagnostic and therapeutic strategies, advancing personalized treatment approaches for fertility restoration.</p>
      </abstract>
      <kwd-group>
        <label>Keywords</label>
        <kwd>BCL2</kwd>
        <kwd>bioinformatics</kwd>
        <kwd>IL1B</kwd>
        <kwd>MMP-9</kwd>
        <kwd>varicocele</kwd>
      </kwd-group>
      <funding-group>
        <funding-statement>No</funding-statement>
      </funding-group>
    </article-meta>
    <notes>
      <sec sec-type="Citation" id="SECID0E2E">
        <title>Citation</title>
        <p>Farid M, Aqilanisa HF. Identifying potential novel biomarkers for varicocele: A bioinformatics approach to genomics analysis. Folia Med (Plovdiv) 2025;67(6):е168906. doi: <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.3897/folmed.67.e168906">10.3897/folmed.67.e168906</ext-link>.</p>
      </sec>
    </notes>
  </front>
  <body>
    <sec sec-type="﻿Introduction" id="SECID0EHF">
      <title>﻿Introduction</title>
      <p>Varicocele is a medical condition in which men’s scrotum veins widen, potentially interfering with testicular hemodynamics.<sup>[<xref ref-type="bibr" rid="B1">1</xref>]</sup> It is estimated to affects 15% to 20% of the world’s male population, and is thought to be responsible for more than 40% of male infertility.<sup>[<xref ref-type="bibr" rid="B2">2</xref>, <xref ref-type="bibr" rid="B3">3</xref>]</sup> Varicocele has a significant impact on the quality of sperm and semen. This condition is caused by instability of the hemodynamics.<sup>[<xref ref-type="bibr" rid="B4">4</xref>, <xref ref-type="bibr" rid="B5">5</xref>]</sup> DNA fragmentation is also known to play a big role in the development of varicocele, which is caused by the process of apoptosis, packaging abnormal chromatin, as well as oxidative stress.<sup>[<xref ref-type="bibr" rid="B6">6</xref>]</sup> In addition, varicocele is known to damage sperm DNA function through improved oxidative stress and hypoxic effects on the testes.<sup>[<xref ref-type="bibr" rid="B7">7</xref>]</sup></p>
      <p>Various risk factors have been identified in the development of varicocele, such as damage to the venous valves, physical activity, weight, hormonal imbalance, smoking habits, and genetic factors.<sup>[<xref ref-type="bibr" rid="B8">8</xref>, <xref ref-type="bibr" rid="B9">9</xref>]</sup> Patients with varicocele are often underdiagnosed, either when presenting with pain or during fertility evaluations.<sup>[<xref ref-type="bibr" rid="B10">10</xref>]</sup> Although most cases of varicocele are asymptomatic, some patients present with a feeling of heaviness in the scrotum or even testicular atrophy.<sup>[<xref ref-type="bibr" rid="B1">1</xref>]</sup></p>
      <p>Mutation genetics, such as single nucleotide polymorphism (<abbrev xlink:title="single nucleotide polymorphism" id="ABBRID0ETH">SNP</abbrev>), are also linked with varicocele and infertility.<sup>[<xref ref-type="bibr" rid="B11">11</xref>]</sup> Previous research has shown that varicocele is the only component of strong genetics, with evidence that the condition can be naturally derived, indicating that genetics play an important role in the development of this disease.‌<sup>[<xref ref-type="bibr" rid="B12">12</xref>, <xref ref-type="bibr" rid="B13">13</xref>]</sup> One of the interesting findings is the improvement in the expression of the enzyme nitric oxide synthase (<abbrev xlink:title="nitric oxide synthase" id="ABBRID0EKAAC">NOS</abbrev>) in patients’ varicocele, which causes relaxation in cells, muscles, and blood vessels.<sup>[<xref ref-type="bibr" rid="B14">14</xref>]</sup> Furthermore, the development of varicoceles is influenced by endothelial vascular overload, which is mediated by prokineticin 2 (<italic><abbrev xlink:title="prokineticin 2" id="ABBRID0EWAAC">PK2</abbrev></italic>) and a number of other genes that contribute to varicocele pathophysiology.<sup>[<xref ref-type="bibr" rid="B15">15</xref>, <xref ref-type="bibr" rid="B16">16</xref>]</sup></p>
      <p>Despite these findings, the genetic basis of varicocele remains poorly understood. Current studies have mainly focused on clinical manifestations and general pathophysiological mechanisms, while the role of specific genetic biomarkers, particularly genes involved in oxidative stress and angiogenesis, is still underexplored and not clearly validated.<sup>[<xref ref-type="bibr" rid="B17">17</xref>]</sup> This lack of precise molecular understanding contributes to delayed diagnosis and limited options for targeted therapy. Therefore, further investigation into genetics associated with varicocele is essential to establish reliable biomarkers that may improve early detection, risk stratification, and personalized treatment strategies. Various genome database platforms are used to investigate single nucleotide polymorphisms (<abbrev xlink:title="single nucleotide polymorphisms" id="ABBRID0ENBAC">SNPs</abbrev>) that have the potential to regulate angiogenesis and oxidative stress, making them biomarkers and treatment targets.<sup>[<xref ref-type="bibr" rid="B18 B19 B20">18–20</xref>]</sup></p>
    </sec>
    <sec sec-type="﻿Aim" id="SECID0EXBAC">
      <title>﻿Aim</title>
      <p>While the results from this study are still limited, the findings provide valuable insights into the potential role of genetic biomarkers in the diagnosis and treatment of varicocele. The aim of this study was to look into genes that may influence the development of varicoceles in the hopes of laying the groundwork for future research. This could eventually lead to the discovery of new biomarkers for diagnosing varicocele that can be used as therapeutic targets in the future.</p>
    </sec>
    <sec sec-type="materials|methods" id="SECID0E3BAC">
      <title>﻿Materials and methods</title>
      <sec sec-type="﻿Identification of variant genes associa-ted with varicocele" id="SECID0EACAC">
        <title>﻿Identification of variant genes associa-ted with varicocele</title>
        <p>This study used three different database platforms, all of which were accessed on February 22, 2025: OpenTarget (<ext-link xlink:href="https://platform.opentargets.org/" ext-link-type="uri" xlink:type="simple">https://platform.opentargets.org/</ext-link>), DisGeNET (<ext-link xlink:href="https://www.disgenet.com/" ext-link-type="uri" xlink:type="simple">https://www.disgenet.com/</ext-link>), and NCBI. The gene variants associated with varicocele were identified by searching the term “varicocele” on each platform to gather relevant genes. The screening process was based on specific criteria: genes with an OpenTarget score greater than 0.03, a DisGeNET score greater than 0.3, and genes expressed in Homo sapiens for NCBI.<sup>[<xref ref-type="bibr" rid="B21">21</xref>]</sup></p>
      </sec>
      <sec sec-type="﻿Discovery of varicocele gene biomarkers" id="SECID0EWCAC">
        <title>﻿Discovery of varicocele gene biomarkers</title>
        <p>Following the identification of key genes, protein-protein interaction (<abbrev xlink:title="protein-protein interaction" id="ABBRID0E3CAC">PPI</abbrev>) analysis was conducted to further explore the biological relationships among the encoded proteins. This analysis utilized the STRING database (<ext-link xlink:href="https://string-db.org/" ext-link-type="uri" xlink:type="simple">https://string-db.org/</ext-link>), accessed on February 22, 2025, which provides comprehensive information on known and predicted functional protein interactions. The resulting <abbrev xlink:title="protein-protein interaction" id="ABBRID0EFDAC">PPI</abbrev> networks were visualized using Cytoscape software, enabling a clear representation of the biological interactions. To identify significant protein clusters, the Molecular Complex Detection (<abbrev xlink:title="Molecular Complex Detection" id="ABBRID0EJDAC">MCODE</abbrev>) plugin in Cytoscape was applied with the following parameters: degree cutoff = 2, node score cutoff = 0.2, k-core = 2, and max depth = 100. Subsequently, to determine the most central and potentially influential proteins within the network, the CytoHubba plugin was employed using the Maximal Clique Centrality (<abbrev xlink:title="Maximal Clique Centrality" id="ABBRID0ENDAC">MCC</abbrev>) algorithm, which is recognized for its accuracy in identifying hub proteins. The top ten ranked nodes were selected, and the three highest-scoring proteins were proposed as potential biomarkers associated with varicocele due to their central role in the interaction network.</p>
      </sec>
      <sec sec-type="﻿Biological pathway analysis and functional roles" id="SECID0ERDAC">
        <title>﻿Biological pathway analysis and functional roles</title>
        <p>To investigate the biological mechanisms involved, the STRING database was accessed on February 22, 2025, to analyze protein-protein interactions. Biological pathway analysis was conducted using WebGestalt, which integrates the KEGG database to link genomic data with high-level functional insights, with significant enrichment determined at a q-value (FDR) cutoff of 0.05. The enrichment results were visualized using a bubble node graph to facilitate interpretation of the connections between various biological components. Additionally, gene function was assessed through Gene Ontology (<abbrev xlink:title="Gene Ontology" id="ABBRID0EXDAC">GO</abbrev>) analysis, which classifies genes into three main categories: Molecular Function (<abbrev xlink:title="Molecular Function" id="ABBRID0E2DAC">MF</abbrev>), Biological Process (<abbrev xlink:title="Biological Process" id="ABBRID0E6DAC">BP</abbrev>), and Cellular Component (<abbrev xlink:title="Cellular Component" id="ABBRID0EDEAC">CC</abbrev>). The outcomes of the <abbrev xlink:title="Gene Ontology" id="ABBRID0EHEAC">GO</abbrev> analysis were presented in bar charts to illustrate the functional distribution of the identified target genes in this study.</p>
      </sec>
    </sec>
    <sec sec-type="﻿Results" id="SECID0ELEAC">
      <title>﻿Results</title>
      <sec sec-type="﻿Identification of variant genes associa­ted with varicocele" id="SECID0EPEAC">
        <title>﻿Identification of variant genes associa­ted with varicocele</title>
        <p>Gene screening was carried out using specific inclusion criteria: genes with a relevance score greater than 0.03 in OpenTarget, greater than 0.3 in DisGeNET, and genes restricted to Homo sapiens, according to NCBI. The identification of gene variants associated with varicocele was conducted through three major databases: OpenTarget, DisGeNET, and NCBI, accessed on February 22, 2025. A search using the keyword “varicocele” yielded 415 genes from OpenTarget, 204 from DisGeNET, and 93 from NCBI. After applying the filtering criteria, 25 genes from OpenTarget, 4 from DisGeNET, and 58 from NCBI met the specified thresholds and were specific to Homo sapiens. Removing duplicates, a total of 79 unique genes were identified <bold>(Table <xref ref-type="table" rid="T1">1</xref>)</bold>. These findings provide an initial set of candidate genes potentially associated with varicocele, serving as a basis for further functional and network analyses.</p>
        <table-wrap id="T1" position="float" orientation="portrait">
          <label>Table 1.</label>
          <caption>
            <p>Varicocele-associated genes</p>
          </caption>
          <table id="TID0EBDAE" rules="all">
            <tbody>
              <tr>
                <td rowspan="1" colspan="1"><italic>DNMT</italic>1</td>
                <td rowspan="1" colspan="1">
                  <italic>ARE</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>
                    <abbrev xlink:title="matrix metalloproteinase-9" id="ABBRID0E1FAC">MMP-9</abbrev>
                  </italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>CST3</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>CACNA1H</italic>
                </td>
              </tr>
              <tr>
                <td rowspan="1" colspan="1">
                  <italic>NLRP3</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>NOS2</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>PTGS2</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>CHI3L1</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>PEPD</italic>
                </td>
              </tr>
              <tr>
                <td rowspan="1" colspan="1">
                  <italic>BRD2</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>SLC7A11</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>GSTM1</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>FASLG</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>GNLY</italic>
                </td>
              </tr>
              <tr>
                <td rowspan="1" colspan="1">
                  <italic>TRPV1</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>INSL3</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>
                    <abbrev xlink:title="B-cell lymphoma 2" id="ABBRID0E1IAC">BCL2</abbrev>
                  </italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>SHH</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>HSPA2</italic>
                </td>
              </tr>
              <tr>
                <td rowspan="1" colspan="1">
                  <italic>TLR4</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>IL-1A</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>MTOR</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>AMH</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>ARG2</italic>
                </td>
              </tr>
              <tr>
                <td rowspan="1" colspan="1">
                  <italic>PROK2</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>CATSPER1</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>GSTT1</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>MIR210</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>SERPINA5</italic>
                </td>
              </tr>
              <tr>
                <td rowspan="1" colspan="1">
                  <italic>PPARA</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>DNMT3B</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>IGF1</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>NOS1</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>SRD5A1</italic>
                </td>
              </tr>
              <tr>
                <td rowspan="1" colspan="1">
                  <italic>
                    <abbrev xlink:title="Interleukin-1 beta" id="ABBRID0ENMAC">IL1B</abbrev>
                  </italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>VEGFA</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>GSTP1</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>RARA</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>PLCZ1</italic>
                </td>
              </tr>
              <tr>
                <td rowspan="1" colspan="1">
                  <italic>PLOD1</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>FASN</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>HMOX1</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>NT5E</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>CACNA1G</italic>
                </td>
              </tr>
              <tr>
                <td rowspan="1" colspan="1">
                  <italic>EGF</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>CASP9</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>PARP1</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>BAK1</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>FGL2</italic>
                </td>
              </tr>
              <tr>
                <td rowspan="1" colspan="1">
                  <italic>FN1</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>TP53</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>ESR2</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>PTGS1</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>NAIP</italic>
                </td>
              </tr>
              <tr>
                <td rowspan="1" colspan="1">
                  <italic>PHGDH</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>TGFB1</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>FAS</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>EDNRB</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>PRM1</italic>
                </td>
              </tr>
              <tr>
                <td rowspan="1" colspan="1">
                  <italic>DNMT3A</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>MTHFR</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>EDN1</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>AGTR2</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>PRM2</italic>
                </td>
              </tr>
              <tr>
                <td rowspan="1" colspan="1">
                  <italic>LEP</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>NOS3</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>TNFSF10</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>EDNRA</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>MIR210HG</italic>
                </td>
              </tr>
              <tr>
                <td rowspan="1" colspan="1">
                  <italic>HIF1A</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>ESR1</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>BAX</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>TXNRD1</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>ROPN1</italic>
                </td>
              </tr>
              <tr>
                <td rowspan="1" colspan="1">
                  <italic>CASP1</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>ACE</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>AHR</italic>
                </td>
                <td rowspan="1" colspan="1">
                  <italic>INHBB</italic>
                </td>
                <td rowspan="1" colspan="1"/>
              </tr>
            </tbody>
          </table>
        </table-wrap>
      </sec>
      <sec sec-type="﻿Discovery of varicocele gene biomarkers" id="SECID0E6EAE">
        <title>﻿Discovery of varicocele gene biomarkers</title>
        <p>Following the identification of genes associated with varicocele, their corresponding proteins were analyzed using the STRING database to investigate functional associations. The resulting Protein-Protein Interaction (<abbrev xlink:title="Protein-Protein Interaction" id="ABBRID0EFFAE">PPI</abbrev>) network consisted of 76 nodes and 594 edges, with an average node degree of 15.6 and an average local clustering coefficient of 0.68. The network showed statistically significant enrichment (<italic>p</italic>&lt;1.0e-16), suggesting that the observed interactions are unlikely to occur by chance and are functionally meaningful. To identify highly interconnected sub-networks, the <abbrev xlink:title="Molecular Complex Detection" id="ABBRID0ELFAE">MCODE</abbrev> plugin in Cytoscape was utilized, resulting in five distinct protein clusters, each exhibiting unique topological characteristics in <bold>Fig. <xref ref-type="fig" rid="F1">1</xref></bold>. These clusters coded visually as red diamonds (cluster 1), orange ellipses (cluster 2), green polygons (cluster 3), gray octagons (cluster 4), and purple triangles (cluster 5) are potential biomarker modules in varicocele-related molecular networks.</p>
        <fig id="F1" position="float" orientation="portrait">
          <object-id content-type="arpha">2D98E120-56D8-5784-BCA0-8E465A133395</object-id>
          <label>Figure 1.</label>
          <caption>
            <p>Construction of a Protein-Protein interaction (<abbrev xlink:title="protein-protein interaction" id="ABBRID0E4FAE">PPI</abbrev>) network for genes associated with varicocele.</p>
          </caption>
          <graphic xlink:href="foliamedica-67-6-e168906-g001.jpg" position="float" orientation="portrait" xlink:type="simple" id="oo_1494656.jpg">
            <uri content-type="original_file">https://binary.pensoft.net/fig/1494656</uri>
          </graphic>
        </fig>
        <p>Further prioritization of key regulatory proteins was carried out using the CytoHubba plugin with the Maximal Clique Centrality (<abbrev xlink:title="Maximal Clique Centrality" id="ABBRID0EIGAE">MCC</abbrev>) algorithm, which identified the top ten hub proteins <bold>(Fig. <xref ref-type="fig" rid="F2">2</xref>)</bold>. These were visualized using color gradients from dark red to light red, representing decreasing centrality values. Among these, <abbrev xlink:title="Interleukin-1 beta" id="ABBRID0ETGAE">IL1B</abbrev>, <abbrev xlink:title="B-cell lymphoma 2" id="ABBRID0EXGAE">BCL2</abbrev>, and <abbrev xlink:title="matrix metalloproteinase-9" id="ABBRID0E2GAE">MMP-9</abbrev> emerged as the most central proteins and are proposed as potential biomarkers for varicocele diagnosis or targeted therapeutic development.</p>
        <fig id="F2" position="float" orientation="portrait">
          <object-id content-type="arpha">E2900532-CAEA-50DA-8FF6-79817DFA0A9D</object-id>
          <label>Figure 2.</label>
          <caption>
            <p>Higher degree values are represented by red to yellow colors. The rank of genes was identified from CytoHubba. The higher the score, the higher the rank of the biomarker.</p>
          </caption>
          <graphic xlink:href="foliamedica-67-6-e168906-g002.jpg" position="float" orientation="portrait" xlink:type="simple" id="oo_1494657.jpg">
            <uri content-type="original_file">https://binary.pensoft.net/fig/1494657</uri>
          </graphic>
        </fig>
      </sec>
      <sec sec-type="﻿Biological pathway analysis and functional roles" id="SECID0ELHAE">
        <title>﻿Biological pathway analysis and functional roles</title>
        <p>To further explore the biological significance of these proteins, functional and pathway enrichment analyses were conducted. STRING-based KEGG pathway analysis revealed 115 enriched biological pathways, with the top three being pathways in cancer, platinum drug resistance, and the HIF-1 signaling pathway. These findings suggest a link between varicocele and oxidative stress responses, as well as cellular mechanisms commonly implicated in tumorigenesis. Gene Ontology analysis provided additional insights into the functional roles of the proteins. A total of 1,205 biological processes (<abbrev xlink:title="Biological Process" id="ABBRID0ERHAE">BP</abbrev>) were enriched, with the most significant being blood pressure regulation, response to hypoxia, and response to oxygen levels highlighting the relevance of these pathways in varicocele-related vascular and metabolic dysregulation. In the Molecular Function (<abbrev xlink:title="Molecular Function" id="ABBRID0EVHAE">MF</abbrev>) category, 55 functions were identified, notably including transcription coactivator binding, nuclear receptor activity, and DNA methyltransferase activity, which indicate involvement in gene expression control and epigenetic regulation. Meanwhile, the Cellular Component (<abbrev xlink:title="Cellular Component" id="ABBRID0EZHAE">CC</abbrev>) analysis identified 29 significant components, such as the inflammasome complex, voltage-gated calcium channel complex, and platelet alpha granules, pointing to the role of inflammation and vascular pressure modulation in varicocele pathophysiology <bold>(Fig. <xref ref-type="fig" rid="F3">3</xref>)</bold>.</p>
        <fig id="F3" position="float" orientation="portrait">
          <object-id content-type="arpha">7E3255C3-7F45-593C-B735-A815BB35E140</object-id>
          <label>Figure 3.</label>
          <caption>
            <p>Visualization of the top 10 biological pathways and functional roles: (<bold>A</bold>) KEGG pathway, (<bold>B</bold>) biological processes, (<bold>C</bold>) molecular function, (<bold>D</bold>) cellular components.</p>
          </caption>
          <graphic xlink:href="foliamedica-67-6-e168906-g003.jpg" position="float" orientation="portrait" xlink:type="simple" id="oo_1494658.jpg">
            <uri content-type="original_file">https://binary.pensoft.net/fig/1494658</uri>
          </graphic>
        </fig>
      </sec>
    </sec>
    <sec sec-type="﻿Discussion" id="SECID0EZIAE">
      <title>﻿Discussion</title>
      <p>Varicocele is still a medical condition with unknown causes, though genetic factors are becoming more widely recognized as influencing its development.<sup>[<xref ref-type="bibr" rid="B20">20</xref>]</sup> This study aimed to identify key genetic biomarkers associated with varicocele using integrated bioinformatics approaches, including STRING, <abbrev xlink:title="Molecular Complex Detection" id="ABBRID0EGJAE">MCODE</abbrev>, and CytoHubba, which have been effectively applied in related conditions such as male infertility and liver fibrosis.<sup>[<xref ref-type="bibr" rid="B21">21</xref>, <xref ref-type="bibr" rid="B22">22</xref>]</sup> Our analysis identified three significant candidate biomarkers: IL-1B, <abbrev xlink:title="B-cell lymphoma 2" id="ABBRID0EVJAE">BCL2</abbrev>, and <abbrev xlink:title="matrix metalloproteinase-9" id="ABBRID0EZJAE">MMP-9</abbrev>.</p>
      <p>Interleukin-1 beta (IL-1B) is a proinflammatory cytokine critically involved in vascular pathology, including varicocele, characterized by abnormal venous dilation in the scrotum. It mediates inflammatory signaling through the NF-κB pathway, contributing to endothelial dysfunction and vascular inflammation observed in related diseases like atherosclerosis.<sup>[<xref ref-type="bibr" rid="B23">23</xref>]</sup> Such mechanisms plausibly underpin varicocele progression via chronic vascular inflammation and increased endothelial permeability, which may exacerbate venous damage.<sup>[<xref ref-type="bibr" rid="B24">24</xref>]</sup> Additionally, IL-1B promotes oxidative stress by stimulating reactive oxygen species (ROS) production, impairing endothelial function and exacerbating venous hypertension typical of varicocele.<sup>[<xref ref-type="bibr" rid="B25">25</xref>]</sup> Its role in angiogenesis, through upregulation of VEGF in macrophages, further implicates IL-1β in abnormal vessel formation.<sup>[<xref ref-type="bibr" rid="B26">26</xref>]</sup> These multifaceted functions position IL-1B as a promising biomarker and potential therapeutic target, especially for inflammation-driven varicocele.</p>
      <p>B-cell lymphoma 2 (<abbrev xlink:title="B-cell lymphoma 2" id="ABBRID0E4KAE">Bcl-2</abbrev>), an anti-apoptotic protein, regulates cell survival and is linked to angiogenesis, oxidative stress, and autophagy. Increased <abbrev xlink:title="B-cell lymphoma 2" id="ABBRID0EBLAE">BCL2</abbrev> expression has been observed in spermatic veins of varicocele patients, suggesting enhanced endothelial cell survival that may contribute to pathological angiogenesis and venous abnormalities.<sup>[<xref ref-type="bibr" rid="B27">27</xref>]</sup> The imbalance between <abbrev xlink:title="B-cell lymphoma 2" id="ABBRID0EMLAE">BCL2</abbrev> and pro-apoptotic BAX exacerbates oxidative damage to vascular and testicular tissues, potentially impairing sperm quality and fertility.‌<sup>[<xref ref-type="bibr" rid="B28">28</xref>]</sup> Furthermore, <abbrev xlink:title="B-cell lymphoma 2" id="ABBRID0EXLAE">BCL2</abbrev>’s interaction with hypoxia-inducible factors such as HIF-1α<sup>[<xref ref-type="bibr" rid="B29">29</xref>]</sup> reflects a cellular adaptive response to hypoxic stress caused by disrupted blood flow, stimulating angiogenesis and cell proliferation.<sup>[<xref ref-type="bibr" rid="B30">30</xref>]</sup> These insights highlight <abbrev xlink:title="B-cell lymphoma 2" id="ABBRID0EJMAE">BCL2</abbrev>’s critical role in maintaining venous homeostasis under varicocele-associated stress.</p>
      <p><abbrev xlink:title="matrix metalloproteinase-9" id="ABBRID0EPMAE">MMP-9</abbrev> is a matrix metalloproteinase involved in extracellular matrix remodeling critical for vascular integrity.‌<sup>[<xref ref-type="bibr" rid="B31">31</xref>]</sup> Genetic polymorphisms in <abbrev xlink:title="matrix metalloproteinase-9" id="ABBRID0E1MAE">MMP-9</abbrev> have been linked to varicose veins, a pathophysiologically related condition, and may similarly predispose individuals to varicocele through altered extracellular matrix degradation leading to vessel dilation and valve incompetence.<sup>[<xref ref-type="bibr" rid="B32">32</xref>]</sup><abbrev xlink:title="matrix metalloproteinase-9" id="ABBRID0EFNAE">MMP-9</abbrev> activation by ROS generated under heat stress and venous congestion damages collagen and elastin fibers, reducing vascular elasticity and exacerbating varicocele progression.‌<sup>[<xref ref-type="bibr" rid="B32">32</xref>]</sup> Moreover, dysregulated <abbrev xlink:title="matrix metalloproteinase-9" id="ABBRID0EQNAE">MMP-9</abbrev> activity disrupts angiogenesis, leading to inefficient collateral vessel formation and persistent venous congestion. The imbalance between <abbrev xlink:title="matrix metalloproteinase-9" id="ABBRID0EUNAE">MMP-9</abbrev> and its inhibitor TIMP2 further aggravates this pathological process, underscoring <abbrev xlink:title="matrix metalloproteinase-9" id="ABBRID0EYNAE">MMP-9</abbrev>’s central role in varicocele pathogenesis.</p>
      <p>Previous studies have investigated the genetic underpinnings of male infertility, either through targeted analysis of specific polymorphisms, such as the NOS3 variants linked to oxidative stress in varicocele patients, or by identifying broader infertility-associated genes using genome-wide bioinformatics approaches.<sup>[<xref ref-type="bibr" rid="B33">33</xref>]</sup> However, these efforts did not specifically isolate or prioritize genetic biomarkers uniquely relevant to varicocele as a distinct pathological entity. The bioinformatics study, for instance, focused on general male infertility without distinguishing varicocele-associated molecular mechanisms, identifying genes like TEX11 and SYCP3 primarily involved in meiosis and spermatogenesis. While previous studies have noted the involvement of these proteins in vascular disorders, our integrative bioinformatics approach uniquely prioritizes <abbrev xlink:title="Interleukin-1 beta" id="ABBRID0EFOAE">IL1B</abbrev>, <abbrev xlink:title="B-cell lymphoma 2" id="ABBRID0EJOAE">BCL2</abbrev>, and <abbrev xlink:title="matrix metalloproteinase-9" id="ABBRID0ENOAE">MMP-9</abbrev> as varicocele-specific biomarkers by combining multi-database gene screening with network and functional analyses. This provides novel insights into their potential interplay in varicocele development, bridging gaps between genetic predisposition and molecular mechanisms. However, this study’s reliance on in silico data limits direct clinical applicability without experimental validation.</p>
      <p>This study’s limitations include its reliance solely on computational bioinformatics analysis without experimental validation, which limits the ability to confirm the actual expression or functional impact of the identified biomarkers <abbrev xlink:title="Interleukin-1 beta" id="ABBRID0ETOAE">IL1B</abbrev>, <abbrev xlink:title="B-cell lymphoma 2" id="ABBRID0EXOAE">BCL2</abbrev>, and <abbrev xlink:title="matrix metalloproteinase-9" id="ABBRID0E2OAE">MMP-9</abbrev> in varicoceles. The absence of patient-derived clinical samples and functional assays means that the biological roles of these genes in varicocele pathogenesis remain hypothetical. Additionally, the cross-sectional design limits interpretation of causality and dynamic gene expression changes over time or in response to treatment. To address these limitations, future studies should include molecular validation using patient tissues, assess biomarker specificity through case-control comparisons, and apply longitudinal designs to evaluate changes post-varicocelectomy. Incorporating predictive modeling tools such as ROC analysis and regression frameworks would also help determine the diagnostic and prognostic utility of these candidate genes in clinical practice.</p>
    </sec>
    <sec sec-type="﻿Conclusion" id="SECID0E6OAE">
      <title>﻿Conclusion</title>
      <p>This study successfully identified 79 varicocele-associated genes by integrating data from OpenTarget, DisGeNET, and NCBI databases. Subsequent protein-protein interaction network analysis and module detection pinpointed three key biomarkers, <abbrev xlink:title="Interleukin-1 beta" id="ABBRID0EFPAE">IL1B</abbrev>, <abbrev xlink:title="B-cell lymphoma 2" id="ABBRID0EJPAE">BCL2</abbrev>, and <abbrev xlink:title="matrix metalloproteinase-9" id="ABBRID0ENPAE">MMP-9</abbrev>, which are critically involved in varicocele pathogenesis through mechanisms such as inflammation, oxidative stress, angiogenesis, and vascular damage. These biomarkers hold promise for improving varicocele diagnosis and developing targeted therapeutic strategies. Overall, our findings enhance the molecular understanding of varicocele and lay the groundwork for future research aimed at more precise and effective clinical interventions.</p>
    </sec>
    <sec sec-type="﻿Funding" id="SECID0ERPAE">
      <title>﻿Funding</title>
      <p>This research was conducted without any external financial support.</p>
    </sec>
    <sec sec-type="﻿Author contributions" id="SECID0EWPAE">
      <title>﻿Author contributions</title>
      <p>M.F. conceptualized and designed the study, conducted the majority of data analysis, interpreted the results, and wrote the original draft of the manuscript. H.F.A. assisted in data collection, contributed to the bioinformatics analysis, and participated in reviewing and editing the manuscript. Both authors read and approved the final manuscript and agree to be accountable for all aspects of the work.</p>
    </sec>
    <sec sec-type="﻿Competing interests" id="SECID0E2PAE">
      <title>﻿Competing interests</title>
      <p>The authors have declared that no competing interests exist.</p>
    </sec>
  </body>
  <back>
    <ack>
      <title>﻿Acknowledgements</title>
      <p>Special thanks to those who continuously inspired and motivated the authors throughout the research process. The data supporting the findings of this study can be obtained from the authors upon reasonable request.</p>
    </ack>
    <ref-list>
      <title>﻿References</title>
      <ref id="B1">
        <mixed-citation xlink:type="simple">1. Leslie SW, Sajjad H, Siref LE. Varicocele. [Updated 2023 Nov 13]. In: StatPearls [Internet]. Treasure Island (FL): StatPearls Publishing; 2025 Jan-. Available from: <ext-link xlink:href="https://www.ncbi.nlm.nih.gov/books/NBK448113/" ext-link-type="uri" xlink:type="simple">https://www.ncbi.nlm.nih.gov/books/NBK448113/</ext-link>.</mixed-citation>
      </ref>
      <ref id="B2">
        <mixed-citation xlink:type="simple">2. Alsaikhan B, Alrabeeah K, Delouya G, et al. Epidemiology of varicocele. Asian J Androl 2016; 18(2):179–81.</mixed-citation>
      </ref>
      <ref id="B3">
        <mixed-citation xlink:type="simple">3. Lomboy JR, Coward RM. The varicocele: clinical presentation, evaluation, and surgical management. Semin Intervent Radiol 2016; 33(3):163–9.</mixed-citation>
      </ref>
      <ref id="B4">
        <mixed-citation xlink:type="simple">4. Jensen CFS, Østergren P, Dupree JM, et al. Varicocele and male infertility. Nat Rev Urol 2017; 14(9):523–33.</mixed-citation>
      </ref>
      <ref id="B5">
        <mixed-citation xlink:type="simple">5. Cannarella R, Shah R, Hamoda TA, et al. Does varicocele repair improve conventional semen parameters? A meta-analytic study of before-after data. World J Mens Health 2023; 42(1):92.</mixed-citation>
      </ref>
      <ref id="B6">
        <mixed-citation xlink:type="simple">6. Esfahani MHN, Tavalaee M. Origin and role of DNA damage in varicocele. Int J Fertil Steril 2012; 6(3):141–6. PMCID: <ext-link xlink:type="simple" ext-link-type="pmcid" xlink:href="PMC3850301">PMC3850301</ext-link>.</mixed-citation>
      </ref>
      <ref id="B7">
        <mixed-citation xlink:type="simple">7. Brannigan RE. Introduction: varicoceles: a contemporary perspective. Fertil Steril 2017; 108(3):361–3.</mixed-citation>
      </ref>
      <ref id="B8">
        <mixed-citation xlink:type="simple">8. Biyani CS, Cartledge J, Janetschek G. Varicocele. Clin Evid [Internet] 2008; Available from: <ext-link xlink:type="simple" ext-link-type="uri" xlink:href="http://www.clinicalevidence.com">www.clinicalevidence.com</ext-link></mixed-citation>
      </ref>
      <ref id="B9">
        <mixed-citation xlink:type="simple">9. Tanrikut C, Goldstein M, Rosoff JS, et al. Varicocele as a risk factor for androgen deficiency and effect of repair. BJU Int 2011; 108(9):1480–4.</mixed-citation>
      </ref>
      <ref id="B10">
        <mixed-citation xlink:type="simple">10. Shridharani A, Owen RC, Elkelany OO, et al. The significance of clinical practice guidelines on adult varicocele detection and management. Asian J Androl 2016; 18(2):269–75.</mixed-citation>
      </ref>
      <ref id="B11">
        <mixed-citation xlink:type="simple">11. Heidari MM, Khatami M, Talebi AR, et al. Mutation analysis of TNP1 gene in infertile men with varicocele. Iran J Reprod Med 2014; 12(4):236.</mixed-citation>
      </ref>
      <ref id="B12">
        <mixed-citation xlink:type="simple">12. Raman JD, Walmsley K, Goldstein M. Inheritance of varicoceles. Urology 2005; 65(6):1186–9.</mixed-citation>
      </ref>
      <ref id="B13">
        <mixed-citation xlink:type="simple">13. Mokhtari G, Pourreza F, Falahatkar S, et al. Comparison of prevalence of varicocele in first-degree relatives of patients with varicocele and male kidney donors. Urology 2008; 71(4):666–8.</mixed-citation>
      </ref>
      <ref id="B14">
        <mixed-citation xlink:type="simple">14. Kahraman CY, Tasdemir S, Sahin I, et al. The relationship between endothelial nitric oxide synthase gene (NOS3) polymorphisms, NOS3 expression, and varicocele. Genet Test Mol Biomarkers 2016; 20(4):191–6.</mixed-citation>
      </ref>
      <ref id="B15">
        <mixed-citation xlink:type="simple">15. Lecouter J, Lin R, Tejada M, et al. The endocrine-gland-derived VEGF homologue Bv8 promotes angiogenesis in the testis: Localization of Bv8 receptors to endothelial cells. Proc Natl Acad Sci USA [Internet] 2003; Available from: <ext-link xlink:type="simple" ext-link-type="uri" xlink:href="http://www.pnas.org">www.pnas.org</ext-link></mixed-citation>
      </ref>
      <ref id="B16">
        <mixed-citation xlink:type="simple">16. Santana VP, Miranda-Furtado CL, de Oliveira-Gennaro FG, et al. Genetics and epigenetics of varicocele pathophysiology: an overview. J Assist Reprod Genet 2017; 34(7):839–47.</mixed-citation>
      </ref>
      <ref id="B17">
        <mixed-citation xlink:type="simple">17. Chiba K, Ramasamy R, Lamb DJ, et al. The varicocele: Diagnostic dilemmas, therapeutic challenges and future perspectives. Asian J Androl 2016; 18(2):276–81.</mixed-citation>
      </ref>
      <ref id="B18">
        <mixed-citation xlink:type="simple">18. Selvam MKP, Samanta L, Agarwal A. Functional analysis of differentially expressed acetylated spermatozoal proteins in infertile men with unilateral and bilateral varicocele. Int J Mol Sci 2020; 21(9).</mixed-citation>
      </ref>
      <ref id="B19">
        <mixed-citation xlink:type="simple">19. Zhu J, Ou N, Song Y, et al. Identification and verification of key genes in varicocele rats through high-throughput sequencing and bioinformatics analysis. Andrologia 2020; 52(9).</mixed-citation>
      </ref>
      <ref id="B20">
        <mixed-citation xlink:type="simple">20. Nasr-Esfahani A, Valipour Motlagh A, Adib M, et al. Integrative bioinformatics analysis of the cell division cycle and ribosomal pathways in the rat varicocele: implications for drug discovery. Cell J 2023; 25(10):727–37.</mixed-citation>
      </ref>
      <ref id="B21">
        <mixed-citation xlink:type="simple">21. Junahar D, Dwiputri R, Adikusuma W, et al. Potential biomarker signatures in male infertility: integrative genomic analysis. Egypt J Med Hum Genet 2024; 25(1):39. doi: <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1186/s43042-024-00512-7">10.1186/s43042-024-00512-7</ext-link></mixed-citation>
      </ref>
      <ref id="B22">
        <mixed-citation xlink:type="simple">22. Paulina AJ, Mutiara YV, Irham LM, et al. Genetic-driven biomarkers for liver fibrosis through bioinformatic approach. Egypt J Med Hum Genet 2024;25(1):58. doi: <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1186/s43042-024-00528-z">10.1186/s43042-024-00528-z</ext-link></mixed-citation>
      </ref>
      <ref id="B23">
        <mixed-citation xlink:type="simple">23. Shen J, Zhao M, Zhang C, et al. IL-1β in atherosclerotic vascular calcification: From bench to bedside. Int J Biol Sci 2021; 17(15):4353–64.</mixed-citation>
      </ref>
      <ref id="B24">
        <mixed-citation xlink:type="simple">24. Fahey E, Doyle SL. IL-1 family cytokine regulation of vascular permeability and angiogenesis. Front Immunol 2019; 10:1426. doi: <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.3389/fimmu.2019.01426/full">10.3389/fimmu.2019.01426/full</ext-link></mixed-citation>
      </ref>
      <ref id="B25">
        <mixed-citation xlink:type="simple">25. Mantsounga CS, Lee C, Neverson J, et al. Macrophage IL-1β promotes arteriogenesis by autocrine STAT3- and NF-κB-mediated transcription of pro-angiogenic VEGF-A. Cell Rep 2022;38(5). doi: <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1016/j.celrep.2022.110309">10.1016/j.celrep.2022.110309</ext-link></mixed-citation>
      </ref>
      <ref id="B26">
        <mixed-citation xlink:type="simple">26. Dinarello CA. A clinical perspective of IL-1β as the gatekeeper of inflammation. Eur J Immunol 2011; 41(5):1203–17.</mixed-citation>
      </ref>
      <ref id="B27">
        <mixed-citation xlink:type="simple">27. Lee JD, Jeng SY, Lee TH. Bcl-2 overexpression in the internal spermatic vein of patients with varicocele. J Formos Med Assoc 2007; 106(4):308–12.</mixed-citation>
      </ref>
      <ref id="B28">
        <mixed-citation xlink:type="simple">28. Mostafa T, Rashed L, Nabil N, et al. Seminal BAX and BCL2 gene and protein expressions in infertile men with varicocele. Urology 2014; 84(3):590–5.</mixed-citation>
      </ref>
      <ref id="B29">
        <mixed-citation xlink:type="simple">29. Lee JD, Yang WK, Lee TH. Increased expression of hypoxia-inducible factor-1 alpha and Bcl-2 in varicocele and varicose veins. Ann Vasc Surg 2012; 26(8):1100–5.</mixed-citation>
      </ref>
      <ref id="B30">
        <mixed-citation xlink:type="simple">30. Zhu SM, Rao T, Yang X, et al. Autophagy may play an important role in varicocele. Mol Med Rep 2017; 16(4):5471–9.</mixed-citation>
      </ref>
      <ref id="B31">
        <mixed-citation xlink:type="simple">31. Mazuchová J, Halašová E, Mazuch J, et al. Investigation of association between genetic polymorphisms of MMP2, MMP8, MMP9 and TIMP2 and development of varicose veins in the Slovak population – pilot study. Physiol Res 2021; 69(Supplement 3):S443–54.</mixed-citation>
      </ref>
      <ref id="B32">
        <mixed-citation xlink:type="simple">32. Serra R, Buffone G, Costanzo G, et al. Altered metalloproteinase-9 expression as least common denominator between varicocele, inguinal hernia, and chronic venous disorders. Ann Vasc Surg 2014; 28(3):705–9.</mixed-citation>
      </ref>
      <ref id="B33">
        <mixed-citation xlink:type="simple">33. Omar SS, Mahfouz W, Dawood W, et al. Relation of nitric oxide synthase gene (NOS3) polymorphisms to varicocele risk and post-varicocelectomy seminal oxidative stress reduction. Andrologia 2020; 52(3).</mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>
