Inherited genetic risk factors play a role in multiple myeloma (MM), yet considerable missing heritability exists. Rare risk variants at genome-wide association study (GWAS) loci are a new avenue to explore. Pleiotropy between lymphoid neoplasms (LNs) has been suggested in family history and genetic studies, but no studies have interrogated sequencing for pleiotropic genes or rare risk variants. Sequencing genetically enriched cases can help discover rarer variants. We analyzed exome sequencing in familial or early-onset MM cases to identify rare, functionally relevant variants near GWAS loci for a range of LNs. A total of 149 distinct and significant LN GWAS loci have been published. We identified six recurrent, rare, potentially deleterious variants within 5 kb of significant GWAS single nucleotide polymorphisms in 75 MM cases. Mutations were observed in BTNL2, EOMES, TNFRSF13B, IRF8, ACOXL and TSPAN32. All six genes replicated in an independent set of 255 early-onset MM or familial MM or precursor cases. Expansion of our analyses to the full length of these six genes resulted in a list of 39 rare and deleterious variants, seven of which segregated in MM families. Three genes also had significant rare variant burden in 733 sporadic MM cases compared with 935 control individuals: IRF8 (P = 1.0 × 10-6), EOMES (P = 6.0 × 10-6) and BTNL2 (P = 2.1 × 10-3). Together, our results implicate six genes in MM risk, provide support for genetic pleiotropy between LN subtypes and demonstrate the utility of sequencing genetically enriched cases to identify functionally relevant variants near GWAS loci.

Blastic plasmacytoid dendritic cell neoplasm (BPDCN) is a rare haematopoietic malignancy characterized by dismal prognosis and overall poor therapeutic response. Since the biology of BPDCN is barely understood, our study aims to shed light on the genetic make-up of these highly malignant tumors. Using targeted high-coverage massive parallel sequencing, we investigated 50 common cancer genes in 33 BPDCN samples. We detected point mutations in NRAS (27.3% of cases), ATM (21.2%), MET, KRAS, IDH2, KIT (9.1% each), APC and RB1 (6.1%), as well as in VHL, BRAF, MLH1, TP53 and RET1 (3% each). Moreover, NRAS-, KRAS- and ATM-mutations were found to be mutually exclusive and we observed recurrent mutations in NRAS, IDH2, APC and ATM. CDKN2A deletions were detected in 27.3% of the cases followed by deletions of RB1 (9.1%), PTEN and TP53 (3% each). The mutual exclusive distribution of some mutations may point to different subgroups of BPDCN whose biological significance remains to be explored.

JAK2V617F(+) myeloproliferative neoplasms (MPNs) frequently progress into leukemias, but the factors driving this process are not understood. Here, we find excess Hedgehog (HH) ligand secretion and loss of PTCH2 in myeloproliferative disease, which drives canonical and noncanonical HH-signaling. Interestingly, Ptch2(-/-) mice mimic dual pathway activation and develop a MPN-phenotype with leukocytosis (neutrophils and monocytes), strong progenitor and LKS mobilization, splenomegaly, anemia, and loss of lymphoid lineages. HSCs exhibit increased cell cycling with improved stress hematopoiesis after 5-FU treatment, and this results in HSC exhaustion over time. Cytopenias, LKS loss, and mobilization are all caused by loss of Ptch2 in the niche, whereas hematopoietic loss of Ptch2 drives leukocytosis and promotes LKS maintenance and replating capacity in vitro. Ptch2(-/-) niche cells show hyperactive noncanonical HH signaling, resulting in reduced production of essential HSC regulators (Scf, Cxcl12, and Jag1) and depletion of osteoblasts. Interestingly, Ptch2 loss in either the niche or in hematopoietic cells dramatically accelerated human JAK2V617F-driven pathogenesis, causing transformation of nonlethal chronic MPNs into aggressive lethal leukemias with >30% blasts in the peripheral blood. Our findings suggest HH ligand inhibitors as possible drug candidates that act on hematopoiesis and the niche to prevent transformation of MPNs into leukemias.

Interferons (IFNs) are cytokines with potent antineoplastic and antiviral properties. IFNα has significant clinical activity in the treatment of myeloproliferative neoplasms (MPN), but the precise mechanisms by which it acts are not well understood. Here, we demonstrate that chromatin assembly factor 1 subunit B (CHAF1B), an Unc-51-like kinase 1 (ULK1)-interactive protein in the nuclear compartment of malignant cells, is overexpressed in patients with MPN. Remarkably, targeted silencing of CHAF1B enhances transcription of IFNα-stimulated genes and promotes IFNα-dependent antineoplastic responses in primary MPN progenitor cells. Taken together, our findings indicate that CHAF1B is a promising newly identified therapeutic target in MPN and that CHAF1B inhibition in combination with IFNα therapy might offer a novel strategy for treating patients with MPN.

Myeloproliferative neoplasms (MPNs) are blood cancers that are characterized by the excessive production of mature myeloid cells and arise from the acquisition of somatic driver mutations in haematopoietic stem cells (HSCs). Epidemiological studies indicate a substantial heritable component of MPNs that is among the highest known for cancers1. However, only a limited number of genetic risk loci have been identified, and the underlying biological mechanisms that lead to the acquisition of MPNs remain unclear. Here, by conducting a large-scale genome-wide association study (3,797 cases and 1,152,977 controls), we identify 17 MPN risk loci (P < 5.0 × 10-8), 7 of which have not been previously reported. We find that there is a shared genetic architecture between MPN risk and several haematopoietic traits from distinct lineages; that there is an enrichment for MPN risk variants within accessible chromatin of HSCs; and that increased MPN risk is associated with longer telomere length in leukocytes and other clonal haematopoietic states-collectively suggesting that MPN risk is associated with the function and self-renewal of HSCs. We use gene mapping to identify modulators of HSC biology linked to MPN risk, and show through targeted variant-to-function assays that CHEK2 and GFI1B have roles in altering the function of HSCs to confer disease risk. Overall, our results reveal a previously unappreciated mechanism for inherited MPN risk through the modulation of HSC function.

Pancreatic cystic neoplasms (PCNs) are recognized as precursor lesions of pancreatic cancer, with a marked increase in prevalence. Early detection of malignant PCNs is crucial for improving prognosis; however, current diagnostic methods are insufficient for accurately identifying malignant PCNs. Here, we utilized mass spectrometry (MS)-based glycosite- and glycoform-specific glycoproteomics, combined with proteomics, to explore potential cyst fluid diagnostic biomarkers for PCN. The glycoproteomic and proteomic landscape of pancreatic cyst fluid samples from PCN patients was comprehensively investigated, and its characteristics during the malignant transformation of PCN were analyzed. Under the criteria of screening specific cyst fluid biomarkers for the diagnosis of PCN, a group of cyst fluid glycoprotein biomarkers was identified. Through parallel reaction monitoring (PRM)-based targeted glycoproteomic analysis, we validated these chosen glycoprotein biomarkers in a second cohort, ultimately confirming N-glycosylated PHKB (Asn-935, H5N2F0S0; Asn-935, H4N4F0S0; Asn-935, H5N4F0S0), CEACAM5 (Asn-197, H5N4F0S0) and ATP6V0A4 (Asn-367, H6N4F0S0) as promising diagnostic biomarkers for distinguishing malignant PCNs. These glycoprotein biomarkers exhibited robust performance, with an area under the curve ranging from 0.771 to 0.948. In conclusion, we successfully established and conducted MS-based glycoproteomic analysis to identify novel cyst fluid glycoprotein biomarkers for PCN. These findings hold significant clinical implications, providing valuable insights for PCN decision-making, and potentially offering therapeutic targets for PCN treatment.

Primary cardiac tumours are rare and mostly benign lesions. Recent publications report that cardiac papillary fibroelastomas are the most common benign primary heart tumour, outnumbering myxomas. However, there is no consensus about their aetiology. We investigated the molecular profile of these tumours using next generation sequencing in a cohort of 16 cases. Eleven of 14 (79%) analysable tumours showed mutations of the KRAS oncogene. Our results provide unambiguous evidence that a significant proportion of these lesions are genuine neoplastic tumours caused by an oncogenic driver mutation.

Myeloproliferative neoplasms (MPN) are clonal stem cell associated disorders inclusive of chronic myeloid leukemia (CML), Polycythaemia vera (PV), myelofibrosis (MF), and essential thrombocythemia (ET). They are characterized by increased production of myeloid cells with minimal effects on terminal differentiation but can undergo transformation to acute leukemias. PV is the most common chronic myeloproliferative neoplasm and in the majority of cases is characterized by a V617F point mutation in JAK2. This JAK2 activating mutation is also found in about half the patients with MF and ET. Such aberrant proteins offer great potential for the treatment of these diseases however inhibitors to JAK2 have had limited success in the clinic in terms of curing the disease. We have previously used advanced proteomic techniques to identify drug targets and thus develop novel treatment strategies to distinguish the leukemic clone in both CML and PV. Here, we build on our proteomic data sets to characterize a new target, the receptor tyrosine kinase AXL. AXL is overexpressed in acute myeloid leukemia and importantly small molecule inhibitors have been developed which are currently in clinical trial hence offer the opportunity to repurpose this drug for the treatment of MPNs. We demonstrate that AXL is upregulated and activated in JAK2 associated MPNs. Further we show that inhibition of AXL preferentially kills early hemopoietic stem cells from PV patients and as such represents a promising therapeutic approach for JAK2 driven MPNs.

RNA splicing and epigenetic gene mutations are the most frequent genetic lesions found in patients with myelodysplastic neoplasm (MDS). About 25% of patients present concomitant mutations in such pathways, suggesting a cooperative role in MDS pathogenesis. Importantly, mutations in the splicing factor ZRSR2 frequently associate with alterations in the epigenetic regulator TET2. However, the impact of these concurrent mutations in hematopoiesis and MDS remains unclear. Using CRISPR/Cas9 genetically engineered mice, we demonstrate that Zrsr2m/mTet2-/- promote MDS with reduced penetrance. Animals presented peripheral blood cytopenia, splenomegaly, extramedullary hematopoiesis, and multi-lineage dysplasia, signs consistent with MDS. We identified a myelo-erythroid differentiation block accompanied by an expansion of LT-HSC and MPP2 progenitors. Transplanted animals presented a similar phenotype, thus indicating that alterations were cell-autonomous. Whole-transcriptome analysis in HSPC revealed key alterations in ribosome, inflammation, and migration/motility processes. Moreover, we found the MAPK pathway as the most affected target by mRNA aberrant splicing. Collectively, this study shows that concomitant Zrsr2 mutation and Tet2 loss are sufficient to initiate MDS in mice. Understanding this mechanistic interplay will be crucial for the identification of novel therapeutic targets in the spliceosome/epigenetic MDS subgroup.

Gastroenteropancreatic neuroendocrine neoplasms (GEP-NEN), considered a heterogeneous neoplasia, exhibit ill-defined pathobiology and protean symptomatology and are ubiquitous in location. They are difficult to diagnose, challenging to manage, and outcome depends on cell type, secretory product, histopathologic grading, and organ of origin. A morphologic and molecular genomic review of these lesions highlights tumor characteristics that can be used clinically, such as somatostatin-receptor expression, and confirms features that set them outside the standard neoplasia paradigm. Their unique pathobiology is useful for developing diagnostics using somatostatin-receptor targeted imaging or uptake of radiolabeled amino acids specific to secretory products or metabolism. Therapy has evolved via targeting of protein kinase B signaling or somatostatin receptors with drugs or isotopes (peptide-receptor radiotherapy). With DNA sequencing, rarely identified activating mutations confirm that tumor suppressor genes are relevant. Genomic approaches focusing on cancer-associated genes and signaling pathways likely will remain uninformative. Their uniquely dissimilar molecular profiles mean individual tumors are unlikely to be easily or uniformly targeted by therapeutics currently linked to standard cancer genetic paradigms. The prevalence of menin mutations in pancreatic NEN and P27KIP1 mutations in small intestinal NEN represents initial steps to identifying a regulatory commonality in GEP-NEN. Transcriptional profiling and network-based analyses may define the cellular toolkit. Multianalyte diagnostic tools facilitate more accurate molecular pathologic delineations of NEN for assessing prognosis and identifying strategies for individualized patient treatment. GEP-NEN remain unique, poorly understood entities, and insight into their pathobiology and molecular mechanisms of growth and metastasis will help identify the diagnostic and therapeutic weaknesses of this neoplasia.

Mutations in the BRCA1 tumor suppressor gene, such as 5382insC (BRCA1insC), give carriers an increased risk for breast, ovarian, prostate, and pancreatic cancers. We have previously reported that, in mice, Brca1 deficiency in the hematopoietic system leads to pancytopenia and, as a result, early lethality. We explored the cellular consequences of Brca1-null and BRCA1insC alleles in combination with Trp53 deficiency in the murine hematopoietic system. We found that Brca1 and Trp53 codeficiency led to a highly penetrant erythroproliferative disorder that is characterized by hepatosplenomegaly and by expanded megakaryocyte erythroid progenitor (MEP) and immature erythroid blast populations. The expanded erythroid progenitor populations in both BM and spleen had the capacity to transmit the disease into secondary mouse recipients, suggesting that Brca1 and Trp53 codeficiency provides a murine model of hematopoietic neoplasia. This Brca1/Trp53 model replicated Poly (ADP-ribose) polymerase (PARP) inhibitor olaparib sensitivity seen in existing Brca1/Trp53 breast cancer models and had the benefits of monitoring disease progression and drug responses via peripheral blood analyses without sacrificing experimental animals. In addition, this erythroid neoplasia developed much faster than murine breast cancer, allowing for increased efficiency of future preclinical studies.

Myeloproliferative neoplasms (MPNs) such as chronic myelogenous (CML) and chronic myelomonocytic leukemias (CMML) are frequently induced by tyrosine kinase oncogenes. Although these MPNs are sensitive to tyrosine kinase inhibitors such as imatinib, patients often relapse upon withdrawal of therapy. We used a model of MPN, which is induced by co-expression of the oncoproteins HIP1/PDGFβR (H/P) and AML1/ETO from their endogenous loci, to examine the mechanisms of disease development and recurrence following imatinib withdrawal. Although the MPN displayed a full hematologic response to imatinib, 100% of the diseased mice relapsed upon drug withdrawal. MPN persistence was not due to imatinib resistance mutations in the H/P oncogene or massive gene expression changes. Within 1 week of imatinib treatment, more than 98% of gene expression changes induced by the oncogenes in isolated hematopoietic stem and progenitor cells (lineage(-)Sca-1(+)c-Kit(+) immunophenotype) normalized. Supplementation of imatinib with granulocyte colony-stimulating factor or arsenic trioxide reduced MPN-initiating cell frequencies and the combination of imatinib with arsenic trioxide cured a large fraction of mice with MPNs. In contrast, no mice in the imatinib-treated control cohorts were cured. These data suggest that treatment with a combination of arsenic trioxide and imatinib can eliminate refractory MPN-initiating cells and reduce disease relapse.

Chronic myeloproliferative neoplasms (MPN) frequently evolve to a blast phase (BP) that is almost uniformly resistant to induction chemotherapy or hypomethylating agents. We explored the functional properties, genomic architecture, and cell of origin of MPN-BP initiating cells (IC) using a serial NSG mouse xenograft transplantation model. Transplantation of peripheral blood mononuclear cells (MNC) from 7 of 18 patients resulted in a high degree of leukemic cell chimerism and recreated clinical characteristics of human MPN-BP. The function of MPN-BP ICs was not dependent on the presence of JAK2V617F, a driver mutation associated with the initial underlying MPN. By contrast, multiple MPN-BP IC subclones coexisted within MPN-BP MNCs characterized by different myeloid malignancy gene mutations and cytogenetic abnormalities. MPN-BP ICs in 4 patients exhibited extensive proliferative and self-renewal capacity, as demonstrated by their ability to recapitulate human MPN-BP in serial recipients. These MPN-BP IC subclones underwent extensive continuous clonal competition within individual xenografts and across multiple generations, and their subclonal dynamics were consistent with functional evolution of MPN-BP IC. Finally, we show that MPN-BP ICs originate from not only phenotypically identified hematopoietic stem cells, but also lymphoid-myeloid progenitor cells, which were each characterized by differences in MPN-BP initiating activity and self-renewal capacity.

We have encountered pancreatic tumors with unique histologic features, which do not conform to any of the known tumors of the pancreas or other anatomical sites. We aimed to define their clinicopathologic features and whether they are characterized by recurrent molecular signatures. Eight cases were identified; studied histologically and by immunohistochemistry. Selected cases were also subjected to whole-exome sequencing (WES; n = 4), RNA-sequencing (n = 6), Archer FusionPlex assay (n = 5), methylation profiling using the Illumina MethylationEPIC (850k) array platform (n = 6), and TERT promoter sequencing (n = 5). Six neoplasms occurred in females. The mean age was 43 years (range: 26-75). Five occurred in the head/neck of the pancreas. All patients were treated surgically; none received neoadjuvant/adjuvant therapy. All patients are free of disease after 53 months of median follow-up (range: 8-94). The tumors were well-circumscribed, and the median size was 1.8 cm (range: 1.3-5.8). Microscopically, the unencapsulated tumors had a geographic pattern of epithelioid cell nests alternating with spindle cell fascicles. Some areas showed dense fibrosis, in which enmeshed tumor cells imparted a slit-like pattern. The predominant epithelioid cells had scant cytoplasm and round-oval nuclei with open chromatin. The spindle cells displayed irregular, hyperchromatic nuclei. Mitoses were rare. No lymph node metastases were identified. All tumors were positive for vimentin, CD99 and cytokeratin (patchy), while negative for markers of solid pseudopapillary neoplasm, neuroendocrine, acinar, myogenic/rhabdoid, vascular, melanocytic, or lymphoid differentiation, gastrointestinal stromal tumor as well as MUC4. Whole-exome sequencing revealed no recurrent somatic mutations or amplifications/homozygous deletions in any known oncogenes or tumor suppressor genes. RNA-sequencing and the Archer FusionPlex assay did not detect any recurrent likely pathogenic gene fusions. Single sample gene set enrichment analysis revealed that these tumors display a likely mesenchymal transcriptomic program. Unsupervised analysis (t-SNE) of their methylation profiles against a set of different mesenchymal neoplasms demonstrated a distinct methylation pattern. Here, we describe pancreatic neoplasms with unique morphologic/immunophenotypic features and a distinct methylation pattern, along with a lack of abnormalities in any of key genetic drivers, supporting that these neoplasms represent a novel entity with an indolent clinical course. Given their mesenchymal transcriptomic features, we propose the designation of "sclerosing epithelioid mesenchymal neoplasm" of the pancreas.

PTEN acts as a phosphatase for PIP3 and negatively regulates the PI3K/AKT pathway, and p27(KIP1) is a cyclin-dependent kinase inhibitor that regulates the G1 to S-phase transition by binding to and regulating the activity of cyclin-dependent kinases. Genetic alterations of PTEN or CDKN1B (p27(KIP1)) are common in hematological malignancies. To better understand how mutations in these two genes might cooperate in leukemogenesis, we inactivated both genes in the hematological compartment in mice. Here, we show that the combined inactivation of Pten and Cdkn1b results in a more severe myeloproliferative neoplasm phenotype associated with lower hemoglobin, enlarged spleen and liver, and shorter lifespan compared to inactivation of Pten alone. More severe anemia and increased myeloid infiltration and destruction of the spleen contributed to the earlier death of these mice, and elevated p-AKT, cyclin D1, and cyclin D3 might contribute to the development of this phenotype. In conclusion, PTEN and p27(KIP1) cooperate in tumor suppression in the hematological compartment.

Interferon α (IFNα) is used to treat JAK2V617F-driven myeloproliferative neoplasms (MPNs) but rarely clears the disease. We investigated the IFNα mechanism of action focusing on PML, an interferon target and key senescence gene whose targeting by arsenic trioxide (ATO) drives eradication of acute promyelocytic leukemia. ATO sharply potentiated IFNα-induced growth suppression of JAK2V617F patient or mouse hematopoietic progenitors, which required PML and was associated with features of senescence. In a mouse MPN model, combining ATO with IFNα enhanced and accelerated responses, eradicating MPN in most mice by targeting disease-initiating cells. These results predict potent clinical efficacy of the IFNα+ATO combination in patients and identify PML as a major effector of therapy, even in malignancies with an intact PML gene.

Myeloproliferative neoplasms (MPNs) are characterized by the activated JAK2/STAT pathway. Pleckstrin-2 (Plek2) is a downstream target of the JAK2/STAT5 pathway and is overexpressed in patients with MPNs. We previously revealed that Plek2 plays critical roles in the pathogenesis of JAK2-mutated MPNs. The nonessential roles of Plek2 under physiologic conditions make it an ideal target for MPN therapy. Here, we identified first-in-class Plek2 inhibitors through an in silico high-throughput screening approach and cell-based assays, followed by the synthesis of analogs. Plek2-specific small-molecule inhibitors showed potent inhibitory effects on cell proliferation. Mechanistically, Plek2 interacts with and enhances the activity of Akt through the recruitment of downstream effector proteins. The Plek2-signaling complex also includes Hsp72, which protects Akt from degradation. These functions were blocked by Plek2 inhibitors via their direct binding to the Plek2 dishevelled, Egl-10 and pleckstrin (DEP) domain. The role of Plek2 in activating Akt signaling was further confirmed in vivo using a hematopoietic-specific Pten-knockout mouse model. We next tested Plek2 inhibitors alone or in combination with an Akt inhibitor in various MPN mouse models, which showed significant therapeutic efficacies similar to that seen with the genetic depletion of Plek2. The Plek2 inhibitor was also effective in reducing proliferation of CD34-positive cells from MPN patients. Our studies reveal a Plek2/Akt complex that drives cell proliferation and can be targeted by a class of antiproliferative compounds for MPN therapy.

Myelodysplastic neoplasm (MDS) is a hematopoietic stem cell disorder that may evolve into acute myeloid leukemia. Fatal infection is among the most common cause of death in MDS patients, likely due to myeloid cell cytopenia and dysfunction in these patients. Mutations in genes that encode components of the spliceosome represent the most common class of somatically acquired mutations in MDS patients. To determine the molecular underpinnings of the host defense defects in MDS patients, we investigated the MDS-associated spliceosome mutation U2AF1-S34F using a transgenic mouse model that expresses this mutant gene. We found that U2AF1-S34F causes a profound host defense defect in these mice, likely by inducing a significant neutrophil chemotaxis defect. Studies in human neutrophils suggest that this effect of U2AF1-S34F likely extends to MDS patients as well. RNA-seq analysis suggests that the expression of multiple genes that mediate cell migration are affected by this spliceosome mutation and therefore are likely drivers of this neutrophil dysfunction.

Tumor heterogeneity presents a challenge for inferring clonal evolution and driver gene identification. Here, we describe a method for analyzing the cancer genome at a single-cell nucleotide level. To perform our analyses, we first devised and validated a high-throughput whole-genome single-cell sequencing method using two lymphoblastoid cell line single cells. We then carried out whole-exome single-cell sequencing of 90 cells from a JAK2-negative myeloproliferative neoplasm patient. The sequencing data from 58 cells passed our quality control criteria, and these data indicated that this neoplasm represented a monoclonal evolution. We further identified essential thrombocythemia (ET)-related candidate mutations such as SESN2 and NTRK1, which may be involved in neoplasm progression. This pilot study allowed the initial characterization of the disease-related genetic architecture at the single-cell nucleotide level. Further, we established a single-cell sequencing method that opens the way for detailed analyses of a variety of tumor types, including those with high genetic complex between patients.

Solid pseudopapillary neoplasm (SPN) of the pancreas is a low-grade malignant neoplasm that accounts for ~5% of cystic pancreatic tumors and ~0.9‑2.7% of exocrine pancreatic tumors. The transcription profiling data (GSE43795) of 14 SPN and 6 control samples were downloaded from the Gene Expression Omnibus (GEO) database. Using the Limma package, Student's t‑tests were performed to identify differentially expressed genes (DEGs) between SPN and control samples [with the following criterion: False discovery rate (FDR)<0.01 and log2 fold‑change (FC)≥3]. Pathway and functional enrichment analyses were performed to investigate the biological processes that the DEGs were involved in. Protein‑protein interaction (PPI) network and sub‑network analyses were conducted to comprehensively understand the interactions between DEGs. The screened DEGs were further annotated according to information relating to transcription factors and tumor associated genes (TAGs). A total of 710 upregulated and 710 downregulated DEGs were observed, including 74 transcriptional factors and 124 TAGs. Membrane metallo‑endopeptidase (MME), matrix metalloproteinase (MMP)-2 and MMP‑9 were also identified as key TAGs. Following PPI network analysis, hub nodes of epidermal growth factor receptor (EGFR), proto‑oncogene tyrosine protein kinase Fyn (FYN), c‑JUN (JUN), glucagon (GCG), c‑Myc (MYC) and CD44 were identified, the majority of which participate in the epidermal growth factor receptor (ErbB) and gonadotropin-releasing hormone (GnRH) signaling pathways. A sub‑network involving 70 gene nodes was also identified, with EGFR as the central gene. MME, MMP‑2 and MMP‑9 contribute to proliferative diabetic retinopathy and also involved in SPN. The genes EGFR, FYN, JUN, GCG, MYC and CD44 may therefore be key genes in SPN, and the ErbB and GnRH signaling pathways may be an important contributor to SPN progression.