Advertisement
Review
Open Access | 10.1172/JCI186702
1Computational Cancer Genomics Team, Genomic Epidemiology Branch, International Agency for Research on Cancer (IARC-WHO), Lyon, France.
2Department of Pulmonary Medicine, Erasmus MC Cancer institute, University Medical Center, Rotterdam, Netherlands.
3GROW School for Oncology and Developmental Biology, Maastricht University Medical Centre, Maastricht, Netherlands.
4University of Bristol Medical School, Bristol, United Kingdom.
5European Molecular Biology Laboratory Barcelona, Tissue Biology and Disease Modeling, Barcelona, Spain.
6Institute of Tissue Medicine and Pathology, University of Bern, Bern, Switzerland.
7Service d’Oncologie Médicale, Groupement Hospitalier Centre, Institut de Cancérologie des Hospices Civils de Lyon, Lyon, France.
Address correspondence to: Lynnette Fernandez-Cuesta, International Agency for Research on Cancer (IARC-WHO), 25 Av. Tony Garnier, 69007 Lyon, France. Email: fernandezcuestal@iarc.who.int.
Find articles by Fernandez-Cuesta, L. in: JCI | PubMed | Google Scholar
1Computational Cancer Genomics Team, Genomic Epidemiology Branch, International Agency for Research on Cancer (IARC-WHO), Lyon, France.
2Department of Pulmonary Medicine, Erasmus MC Cancer institute, University Medical Center, Rotterdam, Netherlands.
3GROW School for Oncology and Developmental Biology, Maastricht University Medical Centre, Maastricht, Netherlands.
4University of Bristol Medical School, Bristol, United Kingdom.
5European Molecular Biology Laboratory Barcelona, Tissue Biology and Disease Modeling, Barcelona, Spain.
6Institute of Tissue Medicine and Pathology, University of Bern, Bern, Switzerland.
7Service d’Oncologie Médicale, Groupement Hospitalier Centre, Institut de Cancérologie des Hospices Civils de Lyon, Lyon, France.
Address correspondence to: Lynnette Fernandez-Cuesta, International Agency for Research on Cancer (IARC-WHO), 25 Av. Tony Garnier, 69007 Lyon, France. Email: fernandezcuestal@iarc.who.int.
Find articles by Alcala, N. in: JCI | PubMed | Google Scholar
1Computational Cancer Genomics Team, Genomic Epidemiology Branch, International Agency for Research on Cancer (IARC-WHO), Lyon, France.
2Department of Pulmonary Medicine, Erasmus MC Cancer institute, University Medical Center, Rotterdam, Netherlands.
3GROW School for Oncology and Developmental Biology, Maastricht University Medical Centre, Maastricht, Netherlands.
4University of Bristol Medical School, Bristol, United Kingdom.
5European Molecular Biology Laboratory Barcelona, Tissue Biology and Disease Modeling, Barcelona, Spain.
6Institute of Tissue Medicine and Pathology, University of Bern, Bern, Switzerland.
7Service d’Oncologie Médicale, Groupement Hospitalier Centre, Institut de Cancérologie des Hospices Civils de Lyon, Lyon, France.
Address correspondence to: Lynnette Fernandez-Cuesta, International Agency for Research on Cancer (IARC-WHO), 25 Av. Tony Garnier, 69007 Lyon, France. Email: fernandezcuestal@iarc.who.int.
Find articles by Mathian, E. in: JCI | PubMed | Google Scholar
1Computational Cancer Genomics Team, Genomic Epidemiology Branch, International Agency for Research on Cancer (IARC-WHO), Lyon, France.
2Department of Pulmonary Medicine, Erasmus MC Cancer institute, University Medical Center, Rotterdam, Netherlands.
3GROW School for Oncology and Developmental Biology, Maastricht University Medical Centre, Maastricht, Netherlands.
4University of Bristol Medical School, Bristol, United Kingdom.
5European Molecular Biology Laboratory Barcelona, Tissue Biology and Disease Modeling, Barcelona, Spain.
6Institute of Tissue Medicine and Pathology, University of Bern, Bern, Switzerland.
7Service d’Oncologie Médicale, Groupement Hospitalier Centre, Institut de Cancérologie des Hospices Civils de Lyon, Lyon, France.
Address correspondence to: Lynnette Fernandez-Cuesta, International Agency for Research on Cancer (IARC-WHO), 25 Av. Tony Garnier, 69007 Lyon, France. Email: fernandezcuestal@iarc.who.int.
Find articles by Derks, J. in: JCI | PubMed | Google Scholar
1Computational Cancer Genomics Team, Genomic Epidemiology Branch, International Agency for Research on Cancer (IARC-WHO), Lyon, France.
2Department of Pulmonary Medicine, Erasmus MC Cancer institute, University Medical Center, Rotterdam, Netherlands.
3GROW School for Oncology and Developmental Biology, Maastricht University Medical Centre, Maastricht, Netherlands.
4University of Bristol Medical School, Bristol, United Kingdom.
5European Molecular Biology Laboratory Barcelona, Tissue Biology and Disease Modeling, Barcelona, Spain.
6Institute of Tissue Medicine and Pathology, University of Bern, Bern, Switzerland.
7Service d’Oncologie Médicale, Groupement Hospitalier Centre, Institut de Cancérologie des Hospices Civils de Lyon, Lyon, France.
Address correspondence to: Lynnette Fernandez-Cuesta, International Agency for Research on Cancer (IARC-WHO), 25 Av. Tony Garnier, 69007 Lyon, France. Email: fernandezcuestal@iarc.who.int.
Find articles by Thirlwell, C. in: JCI | PubMed | Google Scholar
1Computational Cancer Genomics Team, Genomic Epidemiology Branch, International Agency for Research on Cancer (IARC-WHO), Lyon, France.
2Department of Pulmonary Medicine, Erasmus MC Cancer institute, University Medical Center, Rotterdam, Netherlands.
3GROW School for Oncology and Developmental Biology, Maastricht University Medical Centre, Maastricht, Netherlands.
4University of Bristol Medical School, Bristol, United Kingdom.
5European Molecular Biology Laboratory Barcelona, Tissue Biology and Disease Modeling, Barcelona, Spain.
6Institute of Tissue Medicine and Pathology, University of Bern, Bern, Switzerland.
7Service d’Oncologie Médicale, Groupement Hospitalier Centre, Institut de Cancérologie des Hospices Civils de Lyon, Lyon, France.
Address correspondence to: Lynnette Fernandez-Cuesta, International Agency for Research on Cancer (IARC-WHO), 25 Av. Tony Garnier, 69007 Lyon, France. Email: fernandezcuestal@iarc.who.int.
Find articles by Dayton, T. in: JCI | PubMed | Google Scholar
1Computational Cancer Genomics Team, Genomic Epidemiology Branch, International Agency for Research on Cancer (IARC-WHO), Lyon, France.
2Department of Pulmonary Medicine, Erasmus MC Cancer institute, University Medical Center, Rotterdam, Netherlands.
3GROW School for Oncology and Developmental Biology, Maastricht University Medical Centre, Maastricht, Netherlands.
4University of Bristol Medical School, Bristol, United Kingdom.
5European Molecular Biology Laboratory Barcelona, Tissue Biology and Disease Modeling, Barcelona, Spain.
6Institute of Tissue Medicine and Pathology, University of Bern, Bern, Switzerland.
7Service d’Oncologie Médicale, Groupement Hospitalier Centre, Institut de Cancérologie des Hospices Civils de Lyon, Lyon, France.
Address correspondence to: Lynnette Fernandez-Cuesta, International Agency for Research on Cancer (IARC-WHO), 25 Av. Tony Garnier, 69007 Lyon, France. Email: fernandezcuestal@iarc.who.int.
Find articles by Marinoni, I. in: JCI | PubMed | Google Scholar
1Computational Cancer Genomics Team, Genomic Epidemiology Branch, International Agency for Research on Cancer (IARC-WHO), Lyon, France.
2Department of Pulmonary Medicine, Erasmus MC Cancer institute, University Medical Center, Rotterdam, Netherlands.
3GROW School for Oncology and Developmental Biology, Maastricht University Medical Centre, Maastricht, Netherlands.
4University of Bristol Medical School, Bristol, United Kingdom.
5European Molecular Biology Laboratory Barcelona, Tissue Biology and Disease Modeling, Barcelona, Spain.
6Institute of Tissue Medicine and Pathology, University of Bern, Bern, Switzerland.
7Service d’Oncologie Médicale, Groupement Hospitalier Centre, Institut de Cancérologie des Hospices Civils de Lyon, Lyon, France.
Address correspondence to: Lynnette Fernandez-Cuesta, International Agency for Research on Cancer (IARC-WHO), 25 Av. Tony Garnier, 69007 Lyon, France. Email: fernandezcuestal@iarc.who.int.
Find articles by
Perren, A.
in:
JCI
|
PubMed
|
Google Scholar
|
1Computational Cancer Genomics Team, Genomic Epidemiology Branch, International Agency for Research on Cancer (IARC-WHO), Lyon, France.
2Department of Pulmonary Medicine, Erasmus MC Cancer institute, University Medical Center, Rotterdam, Netherlands.
3GROW School for Oncology and Developmental Biology, Maastricht University Medical Centre, Maastricht, Netherlands.
4University of Bristol Medical School, Bristol, United Kingdom.
5European Molecular Biology Laboratory Barcelona, Tissue Biology and Disease Modeling, Barcelona, Spain.
6Institute of Tissue Medicine and Pathology, University of Bern, Bern, Switzerland.
7Service d’Oncologie Médicale, Groupement Hospitalier Centre, Institut de Cancérologie des Hospices Civils de Lyon, Lyon, France.
Address correspondence to: Lynnette Fernandez-Cuesta, International Agency for Research on Cancer (IARC-WHO), 25 Av. Tony Garnier, 69007 Lyon, France. Email: fernandezcuestal@iarc.who.int.
Find articles by Walter, T. in: JCI | PubMed | Google Scholar
1Computational Cancer Genomics Team, Genomic Epidemiology Branch, International Agency for Research on Cancer (IARC-WHO), Lyon, France.
2Department of Pulmonary Medicine, Erasmus MC Cancer institute, University Medical Center, Rotterdam, Netherlands.
3GROW School for Oncology and Developmental Biology, Maastricht University Medical Centre, Maastricht, Netherlands.
4University of Bristol Medical School, Bristol, United Kingdom.
5European Molecular Biology Laboratory Barcelona, Tissue Biology and Disease Modeling, Barcelona, Spain.
6Institute of Tissue Medicine and Pathology, University of Bern, Bern, Switzerland.
7Service d’Oncologie Médicale, Groupement Hospitalier Centre, Institut de Cancérologie des Hospices Civils de Lyon, Lyon, France.
Address correspondence to: Lynnette Fernandez-Cuesta, International Agency for Research on Cancer (IARC-WHO), 25 Av. Tony Garnier, 69007 Lyon, France. Email: fernandezcuestal@iarc.who.int.
Find articles by Foll, M. in: JCI | PubMed | Google Scholar
Published March 3, 2025 - More info
Neuroendocrine tumors (NETs) are a diverse group of malignancies that can occur in various organs, with a notable prevalence in the lungs and gastrointestinal tract, which are the focus of this Review. Although NETs are rare in individual organs, their incidence has increased over recent decades, highlighting the urgent need for current classification systems to evolve by incorporating recent advances in the understanding of NET biology. Several omics studies have revealed molecular subtypes, which, when integrated into existing classification frameworks, may provide more clinically relevant insights for patients with NETs. This Review examines recent progress in elucidating the biology of NETs, with a particular emphasis on the tumor microenvironment and cells of origin. The existence of different cells of origin, which may contribute to distinct molecular groups, along with profiles of immune infiltration — despite being generally low — could explain the emergence of more aggressive cases and the potential for metastatic progression. Given the molecular heterogeneity of NETs and the diversity of their microenvironments and different cells of origin, there is an urgent need to develop morphomolecular classification systems. Such systems would make it possible to better characterize tumor progression, identify new therapeutic targets, and, ultimately, guide the development of personalized therapies.
Neuroendocrine tumors (NETs) are malignant, well-differentiated, epithelial neuroendocrine neoplasms (NENs) (1). NETs can be found in the upper and lower airways, thymus, digestive system, urinary tract, reproductive organs of both sexes, breast, and skin (2–6) (Figure 1A). The digestive tract and lungs are the most common sites (7–9), and this Review will focus on these. The incidence rate of NETs is around seven new cases per 100,000 people annually, and it has been rising over the past two decades (Figure 1B) (7–10). The prevalence is nearly equal between sexes, with a slight male predominance in cases affecting the digestive tract (9, 10). While most patients are in their sixties to eighties, younger individuals can also be affected, particularly in association with hereditary syndromes, appendiceal tumors, and ovarian NETs linked to dermoid cysts. NETs are rarely seen in infants and children (11, 12). The cause of NETs remains unclear. Nonetheless, hereditary genetic mutations, such as MEN1, VHL, and NF1, are linked to NETs in the thorax and upper digestive tract. Risk factors include a family history of cancer, older age, high body mass index, and specific risk factors shared with nonneuroendocrine cancers, such as smoking and alcohol use (13–15).
Epidemiology of NETs. (A) Anatomic sites where NETs originate. (B) Temporal trend of NET incidence in the United States from the Surveillance, Epidemiology, and End Results (SEER; https://seer.cancer.gov/) database (data from ref. 7, Supplemental Table 1). (C) Five-year survival rate of patients with NETs with distant metastases as a function of anatomic site (data from ref. 7, Supplemental Table 3).
Low-grade (G1–G2) NETs usually follow an indolent course, with up to 90% of patients surviving for five years, depending on the site and stage. High-grade (G3) NETs behave more aggressively and are a recognized entity in the pancreas (16). While NETs generally have low rates of distant metastasis, those that do metastasize significantly impact the prognosis of the patients (Figure 1C). The potential to metastasize or invade nearby tissues varies based on their location, grade, and stage at diagnosis. Approximately, 20% of G1 and 50%–60% of G2 ileal and pancreatic NETs metastasize, although these numbers are a bit lower for lung NETs (9% and 25%, respectively) (17). Lung NETs have a five-year survival rate of over 90% for typical carcinoids (G1) and 60% for atypical carcinoids (G2). Pancreatic low-grade NETs have 10-year survival rates of 40%–50% (18–20). Gastroenteropancreatic NETs (GEP-NETs) have the highest rates of metastasis, with 20% showing liver metastases, followed by 5% with bone metastases and 2% with lung metastases. Lung NETs also tend to metastasize to the liver, albeit less frequently than GEP-NETs (4%–5%), and to the bones (3%). They also have a higher incidence of brain metastases compared with GEP-NETs (2% vs. 0%) (17, 21, 22).
Tumor histology remains the gold standard for the diagnosis and clinical decision-making for NETs, with the WHO classification being the most powerful prognostic tool currently available (23, 24). NETs are diagnosed based on histomorphological characteristics, including organoid growth patterns and “salt-and-pepper” chromatin appearance, combined with evidence of a neuroendocrine (NE) phenotype through the expression of NE markers. However, while the NET system is used in GEP organs, the carcinoid system (typical and atypical) is still used in thoracic organs (Table 1). The grade is determined by both the proliferation rate (measured by mitotic count or Ki-67 index) and the presence of necrosis. While grading is effective for prognosis, it falls short in aiding therapy selection and identifying relapse risk. It has been shown that recurrence risk increases with higher tumor grade and TNM stage (25, 26), but establishing precise cut-offs for current markers has been challenging for GEP-NETs (27–30) and lung NETs. A recent study including over 300 samples from the lungNENomics project (31) assessed the current and emerging criteria for the classification of lung NETs using traditional pathology techniques as well as innovative deep learning approaches on whole-slide images. The authors concluded that, while mitotic criteria can be complemented by manual or automated assessment of Ki-67 or PHH3 proteins expression, these markers do not significantly improve the prognostic value of the current classification and remain highly unspecific for aggressive cases (32). In addition, there are also discrepancies in nomenclature use between different organs that impair communication among clinicians (33). For example, in the pancreas, G3-NETs are a recognized entity, while in the lung, these tumors are still termed as large-cell NE carcinomas (LCNECs) with NET-like morphology (34). It is clear that we have exhausted the potential of current morphological features in assisting the clinical management of NETs and need to investigate using further approaches. Consequently, most patients with NETs undergo surgery followed by prolonged follow-up periods (5–10 years) with conventional or functional imaging tests for monitoring (35), which are burdensome and costly. In addition, there are no defined strategies for adjuvant systemic therapies, not even for more aggressive NETs.
Although many recent molecular findings still need thorough validation, a growing body of evidence suggests that their incorporation into the current classification system via a morphomolecular classification might be more clinically relevant for patients with NETs. A summary of the molecular features of different organs is shown in Table 2.
Lung. Lung NETs, unlike small-cell lung cancer (SCLC) and LCNEC, exhibit a low mutation rate, with frequent alterations in chromatin-remodeling genes and rare mutations in TP53 and RB1 (36–38). Significant mutations include somatic inactivation of MEN1, EIF1AX, and ARID1A, with around 5% of patients having MEN1 syndrome, a rare disorder that predisposes patients to developing tumors in endocrine glands and organs (39, 40). Despite being traditionally viewed as a single entity, a comprehensive multi-omic study by Alcala et al. identified three distinct molecular subgroups within lung NETs (41). Additional studies have identified similar subtypes (42, 43), or subtypes that were highly concordant with these observations (44, 45), validating the findings. Carcinoid A1 (LC1) tumors are characterized by high levels of ASCL1 and DLL3 expression and frequently harbor EIF1AX mutations. Carcinoid A2 (LC3) tumors have low levels of SLIT1 and ROBO1 expression. Both A1 and A2 subgroups have favorable prognoses, with more than 80% of patients surviving beyond 10 years. Carcinoid B (LC2) tumors exhibit high expression levels of UDP-glucuronosyltransferase (UGT), cytochrome P450 (CYP) family gene, ANGPTL3, and ERBB4, along with low levels of OTP and TTF1. These tumors often have MEN1 alterations and a poorer prognosis, with only 60% of patients surviving beyond 10 years. These molecular subtypes do not align with morphological classifications but show a trend toward typical (A1 and A2) and atypical (B) carcinoids. Each subgroup has unique clinical features and potential therapeutic targets, underscoring the need for further research into their biology, risk factors, and responses to treatment. One study also reported the existence of a group of uncommon carcinoids — named supracarcinoids — that exhibit genuine carcinoid-like morphology but molecular and clinical features of highly aggressive and poorly differentiated LCNECs (41). This finding has been replicated in an independent study by Simbolo et al. (46). Further analyses are underway to better characterized this new biological entity. Emerging molecular markers, such as protein expression of CD44, ASCL1, and OTP and TERT gene expression have shown prognostic value in lung NETs, alongside traditional markers like Ki-67 and somatostatin receptors (47–49). Studies have proposed combining these markers to stratify patients into more clinically relevant categories than the typical/atypical classification (50). Therefore, the use of comprehensive IHC panels can further characterize these molecular groups and potentially lead to personalized treatment strategies for patients with lung NETs. Advances in nuclear medicine have also provided new methods for tumor characterization, such as somatostatin receptor targeting and metabolic imaging with 18F-FDG or 18F-DOPA (51). These approaches have led to alternative classifications and new grading schemes (52).
Pancreas. A large whole-genome sequencing study from Scarpa et al. provided comprehensive insights into the genetic landscape of pancreatic NETs (53). The study identified driver genetic alterations that converged into four main pathways: chromatin remodeling, DNA damage repair, activation of mTOR signaling, and telomere maintenance. The authors also found that pancreatic NETs exhibit recurrent genetic inactivation of MEN1, ATRX, and DAXX and the activation of the PI3K/mTOR pathway (53). Approximately 40% of sporadic nonfunctioning pancreatic NETs (tumors that do not produce hormones) harbor mutations both in MEN1 and DAXX or ATRX (53–55). These mutations suggest that there is profound epigenetic dysregulation in pancreatic NET development. DAXX and ATRX mutations correlate with loss of nuclear expression and higher relapse risk (56). Indeed, a subgroup with recurrent loss of heterozygosity (LOH) on 10 specific chromosomes, and enriched for MEN1, DAXX, and ATRX mutations, showed worse prognosis and higher metastatic risk (57). Consequently, DAXX/ATRX immunohistochemistry has been suggested as a prognostic biomarker in the pancreas. In addition, the mutational status of DAXX and ATRX is strongly correlated with the alternative lengthening of telomeres (ALT), with ALT emerging as a reliable indicator of increased risk of metastasis for primary pancreatic NETs, supporting its introduction in clinical practice (58). Interestingly, once metastasized, the loss of DAXX/ATRX and the presence of ALT are associated with longer survival, though the reason for this is unknown (54). TP53 and KRAS mutations were also shown to be more frequent in pancreatic NET metastases than in primary tumors, but TP53 mutations were on the contrary less frequent in liver metastasis from lung NETs than in their primary counterparts (59). Insulinomas, NETs that produce insulin and are predominantly found in the pancreas, frequently exhibit mutations in the transcription factor YY1, which is absent in nonfunctioning pancreatic NETs (60). A notable proportion of clinically sporadic pancreatic NETs were found to have germline mutations. These include previously unreported mutations in DNA repair genes such as MUTYH, CHEK2, and BRCA2 as well as mutations in MEN1 and VHL, occurring in approximately 17% of patients (53).
For G3 pancreatic NETs, the limited data indicate that these tumors often harbor mutations in DAXX/ATRX and TP53, along with the loss of RB1 (61). Their epigenetic profile appears to be similar to that of G1 and G2 tumors (62). Notably, progression from G2 to G3 is not uncommon. While this progression can occasionally occur in naive tumors, it is more frequently observed under therapy (63–65).
Several transcriptomic analyses have proposed various molecular subclassifications, yet their clinical implications remain unclear. Initial analyses identified three main subtypes of pancreatic NETs, including less aggressive insulinomas and two nonfunctional clusters associated with metastatic disease (66). The metastasis-like primary 1 (MLP-1) and MLP-2 subtypes share gene signatures related to fibroblasts, stem cells, and hypoxia, with MLP-1 showing an immune-suppressive profile (67). These profiles could pave the way for new immunotherapy approaches for these tumors. Scarpa et al. also identified a subgroup of tumors associated with hypoxia and HIF signaling, suggesting that stemness, hypoxia, and metabolic changes, along with immune profile alterations, are major phenotypes associated with aggressive pancreatic NETs (53). However, transcriptome profiles do not clearly distinguish between small indolent tumors and those with a high relapse risk. Despite the low mutational burden (68), the frequent mutations in chromatin regulation genes highlight the significant role of epigenetic dysregulation in pancreatic NETs. Chromatin immunoprecipitation assays identified two major subtypes of pancreatic NETs based on H3K27ac enhancer profiles, one with an α signature expressing ARX and one with a β signature expressing PDX1 (69). DNA methylation profiles further stratify tumors by cell of origin, genetic background, and prognosis, identifying α-like tumors (MEN1 mutated and indolent) and intermediate tumors (MEN1 and either DAXX or ATRX mutated and with high risk of relapse), suggesting a progression from α to intermediate tumors (70).
MEN1 in pancreatic and lung NETs. A common occurrence among pancreatic NETs and lung NETs is the presence of both sporadic and germline MEN1 mutations. MEN1 is associated with histone lysine methyltransferase activity and regulation of the cell-cycle pathway through CDKN1B (p27) and CDKN2C (p18). In lung NETs, the identification of a somatic MEN1 mutation correlates with a poorer prognosis (41, 71), though it may be indolent in α-like and sporadic pancreatic NETs with isolated MEN1 deficiency (65). When combined with a DAXX or ATRX mutation, MEN1-inactivated pancreatic NETs have a poor prognosis. Interestingly, patients with germline MEN1 mutations and lung NETs have a favorable prognosis, possibly owing to early screening that identifies 86% of tumors, allowing timely curative treatment (72, 73). The heritability of MEN1 alterations among siblings of patients with MEN1-related lung NET and pancreatic NET is lower compared with other NET types (pituitary, adrenal, thymic) (74), though the reason for this is unclear. Thus, while MEN1 is clearly linked to the development of pancreatic and lung NETs, its precise role in tumor aggressiveness across different anatomical sites and its molecular mechanisms require further investigation.
Small intestine. Small-intestinal NETs also have a very low mutational rate, with the only recurrent mutations occurring in the CDKN1B gene in 8% of cases and loss of chromosome 18 observed in 50%–80% cases (75, 76). Three molecular subtypes have been described based on DNA methylation and genetic alterations. One group harbors chromosome 18 LOH (18LOH group, 55% of tumors), another group shows no large copy number variations (19%), and the third group exhibits multiple copy number variations, including gains of chromosomes 4, 5, and 20 (26%). Notably, tumors with CDKN1B mutations are found within the 18LOH group. A difference in progression-free survival was identified among these three subgroups, with the 18LOH group having longer progression-free survival compared with the other two groups (77). Additionally, small intestine NETs also exhibit DNA methylation clusters associated with different prognoses, CNV variations, and genetic backgrounds (77).
Cells of origin. The heterogeneity of NE cells throughout the body, including pancreatic islet cells and enteroendocrine cells in the intestine, is well established. Each of these NE cell populations is subdivided into specific subtypes primarily defined by the hormones they express. There is emerging evidence of heterogeneity among NE cells in the lung (78–81). Beyond hormonal or neuropeptide expression, heterogeneity may also arise concerning innervation (82). These cells might have different functions compared with the noninnervated ones. Furthermore, a third level of heterogeneity in NE cells can be attributed to their spatial distribution within tissues. For example, NE cells in the trachea likely differ in function from those in the bronchi or bronchioles in the lung, with the latter usually consisting of multiple NE cells clustering as a NE body. In the lung, tumorlets (small carcinoid lesions <5 mm) are considered preneoplastic lesions. These tumorlets are abundant in a condition called diffuse idiopathic NE cell hyperplasia, which is generally diagnosed in older female individuals with peripheral nodules and seems to be associated with the A1 molecular carcinoid subtype and EGFR expression (49, 83). NE cell hyperplasia may also develop more centrally in the airways and could contribute to the formation of centrally located carcinoids, which seems to be more often observed in younger patients and linked to the A2 molecular subtype. Therefore, the spatially distinct locations of NE cells may correlate with different molecularly defined lung NET subtypes, warranting further investigation.
Most NETs originate from NE cells of epithelial origin. These cells arise from local pluripotent stem cells in the endoderm, which undergo NE differentiation under the influence of transcription factors like PDX1 and NGN3 in the GEP system and ASCL1 in the lung. Examples include NETs from enterochromaffin cells in the gut and insulinomas from pancreatic islet cells. Insulinomas are hypothesized to originate from β cells in the islets of Langerhans, as suggested by epigenetic data (84). Genetic mutations in nonfunctioning NETs differ from those in insulinomas, indicating varying susceptibility among different cell types (85).
Owing to the scarcity of NE cells in most organs, our understanding of their subtypes and their association with specific tumor types remains limited. For example, a recent single-cell cancer atlas of the lung identified just 500 NE cells of 2.4 million lung cells, accounting for only 0.02% of the total (86). Several lines of evidence highlight the pivotal role of the cell of origin in determining tumor type following specific oncogenic mutations, particularly for SCLC (87, 88) and GEP-NEC (89, 90). But cells of origin have not been as extensively explored in NETs, partly due to the limited number of available models.
Most NETs exhibit low levels of immune infiltration but may present heterogeneous immune subtypes (41, 91), alongside significant abnormal vascularization (92, 93). Supporting these observations, a recent small single-cell study of three lung NETs revealed a distinct microenvironment characterized by the presence of noninflammatory monocyte-derived myeloid cells, vascular smooth muscle cells, pericytes, and a small proportion of potentially prognostic cancer-associated myofibroblasts (94). Transcriptomic analyses have also found lung supracarcinoids to have high infiltration levels (41). These tumors exhibited increased expression of immune checkpoint genes such as PDL1, a phenomenon observed occasionally in a small subset of NETs across various sites (95, 96). The microenvironment of pancreatic NETs shows a lower density of immune cells (97) than that observed in pancreatic NECs. An increased presence of CD8+ T cells combined with a reduced number of macrophages is associated with better outcomes (98, 99).
Recent pan-cancer analyses have classified immune microenvironments into 6–12 dominant archetypes, which may share evolutionary traits and vulnerabilities, particularly in relation to immunotherapy (100, 101). Among these archetypes, several immune-desert and myeloid-centric environments might be pertinent to NETs, although the very small sample size for NETs (n = 7) allows only for the formulation of hypotheses that require further validation in large cohorts, and the link with known NET molecular subtypes has not been explored. Specifically, pancreatic NETs (4 G1, 1 G2) and a colorectal (G1) NET were categorized as either immune-desert or myeloid-centric archetypes (101) (Figure 2A). The immune-desert monocyte archetype found in pancreatic NETs is associated with higher rates of cancer-associated fibroblasts, aligning with preliminary single-cell data on lung NETs; the only colorectal NET included in the study presented another immune-desert archetype, biased for CD4+ T cells and macrophages. Myeloid-centric environments can be mistaken for immune deserts due to the lack of T cells, explaining why NETs are often considered uniformly “cold” tumors. The myeloid-centric type 1 dendritic cell (DC1) archetype is characterized by elevated levels of type 1 conventional dendritic cells, which might also occur in subsets of lung NETs, in particular in the A1 molecular group, based on indirect evidence from transcriptomic deconvolution (41). This archetype may also feature higher levels of neutrophils compared with other immune archetypes. The existence of lung NETs with high infiltration, in particular in supracarcinoids (41), suggests that some tumors may correspond to the immune-rich or immune-rich stromal archetypes, although the precise archetype remains unclear.
The microenvironment of NETs. (A) Immune archetypes observed in NETs. Data shown are from Combes and colleagues (101). The hot-intermediate-cold classification of archetypes corresponds to that of Galon et al. (163), as reported Combes and colleagues (100) (representation of archetypes adapted with permission of Springer Nature Limited, which retains rights to the reference image). (B) Competing hypotheses about the crosstalk between the cell of origin of NETs and their microenvironment. (Top left) Under the spontaneous parallel evolution hypothesis, the cell of origin is solely responsible for the observed tumor groups (histopathological or molecular types and subtypes), and crosstalk between NE cells of origin and their microenvironment does not influence carcinogenesis. (Top right) Under the spontaneous divergent evolution hypothesis, somatic alterations (genetic or epigenetic) are responsible for the observed tumor groups, putting the NE cells of origin on different evolutionary trajectories. (Bottom left) Under the immune-mediated parallel evolution hypothesis, the cell of origin determines the observed tumor groups, but crosstalk between NE cells of origin and their microenvironment is crucial to initiate carcinogenesis. (Bottom right) Under the immune-mediated divergent evolution hypothesis, different microenvironments are responsible for the observed tumor groups by induction of different selective pressures that put the NE cells of origin on different evolutionary trajectories.
Interaction between cell of origin and microenvironment. NE cells from different tissues show evidence of dynamic crosstalk with their microenvironment, such as the response of NE cells to environmental and cellular derived stimuli and the direct influence they can exert both locally and systemically through the bioactive compounds they secrete (91). However, the interactions between the tumor microenvironment, NE cells, and tumor cells during carcinogenesis and NET progression are not well understood, resulting in unproven hypotheses (Figure 2B). The limited level of immune infiltration is generally interpreted as evidence of the limited impact of immune infiltration on NET evolution. Nevertheless, there is evidence that both immune-desert and myeloid-rich archetypes observed in pancreatic NETs and a colorectal NET strongly influence tumor formation and progression. The myeloid-rich DC1 archetype is speculated to evolve from a fibroblast-macrophage-monocyte axis. It exhibits features similar to those of visceral adipose tissue, where Tregs interact with macrophages to regulate adiposity. In this environment, conventional dendritic cells acquire a tolerogenic phenotype with decreased antigen-presenting functions. This phenotype is sustained by PPARγ, IL-10, and steroid signaling (100). The immune-desert monocyte archetype is speculated to resemble the immune system’s response during midstage wound healing. This stage is characterized by the presence of few T cells but many neutrophils and immature myeloid cells. This pattern aligns with skin wound healing, where neutrophil and monocyte infiltration is initially driven by IL-1 and TGF-β. Subsequently, monocytes differentiate into macrophages, which clear cellular debris and interact with fibroblasts to promote tissue healing through extracellular matrix remodeling (100). Interestingly, immune-desert and myeloid-centric tumors had the highest expression of Ki-67 and of cell-cycle–associated genes in general, suggesting that their lack of T cells might be associated with cell-cycle checkpoint avoidance.
Little is known about the progression from primary to metastatic in lung and pancreatic NETs. Most studies focus on changes in proliferation rates, with over 35% of lung NET and 50% of pancreatic NET metastases showing increased rates compared with primary tumors (22, 102). Few studies assess molecular changes in consecutive specimens, but low genomic heterogeneity between primary tumors and synchronous liver metastases has been observed (103). However, other studies report increased genomic imbalances in metastases (104). Mechanisms associated with metastasis formation include hypoxia, metabolic changes, and a stem-cell-like phenotype (66, 67). A recent study identified immune escape, stem cell signaling, and cell reprogramming as key pathways in progression (105). Additionally, an increased number of tumor-infiltrating T cells in metastatic pancreatic NETs suggests immune pathway activation (106). Comparing primary and metastatic transcriptomes has revealed potential treatment targets, though underlying progression mechanisms remain unclear (107). In the case of small intestine NETs, DNA methylation appears to change progressively from normal tissue to primary tumor and then to metastasis, suggesting that tumor evolution may be driven by epigenetic mechanisms (108).
In addition to well-differentiated NETs, NENs also include poorly differentiated and more aggressive NECs. NE cells are believed to be the origin of most NENs. Their plastic and dynamic nature is highlighted by their ability to respond to external stimuli and microenvironmental signals in a context-specific manner. NE cells of the stomach, intestine, pancreas, and lung have all been shown to be altered in number under certain pathological conditions. However, although they share characteristics, NETs and NECs are usually considered separate entities with distinct evolutionary histories and cells of origin (84, 109, 110). Indeed, normal epithelial cells can directly acquire a small-cell NEC phenotype by acquiring drivers that can reprogram cells (111). Balanis et al. also identified a convergence to a small-cell NE state across various epithelial cancers, often linked to poor prognosis (112). Similarly, treatment is known to influence the transition toward NEC (113), for example, adenocarcinoma of the lung and prostate under prolonged exposure to EGFR blockade or treatment-induced suppression of the androgen receptor in prostate cancer (114–117). The specific genomic profile of NECs across organs, almost always TP53 and often RB1 mutations, especially in SCLC (37, 84), also suggests fundamental differences between NETs and NECs and the necessary mechanisms for a convergent NEC phenotype. Nevertheless, although concurrent inactivation of TP53 and RB1 in prostate cancer was shown to lead to a NE phenotype in preclinical models (118), similar inactivation did not lead to the expression of NE markers in the colon (89), and it seems the situation might be the same in the lung. The existence of TP53 and RB1 in lung adenocarcinoma further supports the fact that these two genes might not be sufficient for NEC formation in all organs (119). It should also be noted that there is speculation of transdifferentiation from NETs to adenocarcinoma in the pancreas, in particular related to the acquisition of KRAS mutations (59).
However, there is some evidence suggesting that the separation between NETs and NECs might be subtler and that the progression of NETs toward a more aggressive NEC-like molecular profile exists. The discovery of supracarcinoids (41), the existence of which was confirmed in other cohorts (43, 46), further supports this link. Although rare, these highly aggressive entities resemble NECs, making it crucial to understand them for identifying patients with poor prognosis and specific treatment needs as well as to explore the link between NETs and NECs. The absence of typical NEC alterations in these samples suggests phenotypic convergence toward NECs through alternative drivers. Their unique microenvironment, composed of myeloid cells and macrophages, like that of LCNEC rather than NETs, may contribute to the ecological niche necessary for this convergence (120). One mechanism through which lung NETs may acquire NEC features is chromothripsis, typically affecting chromosomes 3, 11, and 12 (41, 53, 121). This phenomenon has been identified in up to 3% of SCLCs, primarily in never-smokers, and is associated with a unique molecular profile characterized by intact RB1 and TP53 genes. These tumors were enriched for mutations in cell-cycle gene ATM and chromatin remodeling genes ARID1A, MEN1, and EIF1AX, which are most observed in lung NETs (36). Some patients with this profile were indeed later identified as having lung NETs based on a different tumor sample. The increased oncogenicity was partially attributed to cell-cycle pathway disruption through amplification of CCND1 or CCND2/CKD4/MDM2. Notably, similar amplification of chromosomes 3 and 11, along with enrichment for ATM mutations, has been previously observed in lung NETs/LCNECs with higher proliferation rates (122).
Similar evidence is observed in pancreatic NETs. During progression, α-like tumors can evolve into intermediate tumors upon mutations in DAXX and ATRX, leading to genomic instability and activation of ALT (70). Rarely, metastatic insulinomas show ARX positivity with concurrent loss of DAXX/ATRX and ALT activation, suggesting a distinct tumorigenic mechanism in malignant insulinomas similar to nonfunctional pancreatic NETs through transdifferentiation from α cell tumors (123). In a Tp53- and Rb-mutated mouse model, metastatic and primary tumors arise from low-grade insulinomas via dedifferentiation along the β cell developmental pathway, resulting in downregulation of mature β cell markers and expression of pancreatic progenitor markers (124).
Among the systemic therapies for metastatic NETs, two somatostatin analogues (SSAs) (125, 126) and one peptide receptor radionuclide therapy (177Lu-edotreotide) (127) are accepted (Table 3). Regarding molecularly targeted therapies, everolimus (a mammalian target of rapamycin mTOR inhibitor, refs. 128–130), three tyrosine kinase inhibitors with antiangiogenic activity (sunitinib for pancreatic NETs, refs. 131; surufatinib in China, refs. 132, 133 for all NETs; and probably soon cabozantinib, ref. 134), and belzutifan (for patients with von Hippel–Lindau (VHL) disease, targeting HIF-2α, ref. 135) have been approved. Systemic cytotoxic chemotherapies include alkylating agents such as streptozotocin, temozolomide (combined with capecitabine), dacarbazine, and oxaliplatin.
Despite the several prognostic factors available for NETs management, almost no predictive factors of response to a specific treatment are prospectively validated (136, 137). Except for 177Lu-edotreotide, which requires the expression of somatostatin receptor 2 (SSTR2) on somatostatin receptor imaging, and belzutifan, which requires a germline mutation in the VHL gene, all other treatments are prescribed without driver biomarkers. Promoter methylation or low expression of methylguanine-methyltransferase (MGMT) is the best-proven predictive factor of response to alkylating agents including temozolomide; evaluating its status could help in choosing chemotherapy for NETs (138, 139). However, the assessment of MGMT methylation and expression by validated platforms must become more widely available before larger clinically use.
Through recent advances in our understanding of the molecular alterations occurring in NETs, the options for molecularly targeted therapy of driver mutations remain limited. This is because some NET subtypes harbor very few molecular alterations (small intestine NETs) and the more frequently occurring alterations are not actionable (MEN1, DAXX/ATRX, and ARID1A in pancreatic and lung NETs). Subsequently, achieving effective personalized medicine in NETs remains a challenge compared with other cancers. In a recent study from Boilève and colleagues, 19 patients (four with a NEC, 15 with a NET) were treated with molecularly targeted therapy as follows: immunotherapy (n = 3), tipifarnib (n = 1), NOTCH inhibitor (n = 1), EGFR inhibitor (n = 2), HER2 inhibitor (n = 1), and everolimus (n = 11); clinical benefit was seen in 67% of cases (140). In addition, positive outcomes have been seen in rare subtypes of NETs treated with an agnostic approach targeting rare gene rearrangements, including ALK, ROS, RET, NRG1, and NTRK1 (141, 142). The prevalence of these molecular alterations may be associated with the organ of origin. For example, BRAF-V600E mutation is more common in colon NECs than in colon adenocarcinoma (143), and it has also been identified in a supracarcinoid (143). In lung NETs, the presence of a targetable driver (e.g., ALK or RET) may be associated with a mucin-enriched or combined adenocarcinoma tumor, warranting further investigation (145, 146). The discovery of gene fusions and oncogenic driver mutations in a subset of NETs underscores the importance of genomic evaluation in these tumors. Finally, using “liquid biopsy” to identify the predominant NET clones to determine actionable alterations may help identify patients for personalized medicine (147, 148).
In lung NET, new insights may indicate drug susceptibility according to epigenetic/transcriptional molecular subtypes. The recent observation that some lung NETs express EGFR and that lung NET patient-derived tumor organoids require EGF for their growth suggests the need for clinical studies to determine whether EGFR could be a predictive biomarker for the response of a subset of lung NETs to EGFR-targeted therapies (144). In addition to EGFR, the expression of delta-like ligand 3 (DLL3), an inhibitory Notch pathway ligand, has been identified in lung NETs (149). DLL3 is significantly higher expressed in A1 but not in A2/B lung NET subtypes (41, 49). Proof of concept was recently highlighted in a case report of a patient with atypical lung NET showing clinical efficacy on treatment with a DLL3 bispecific T cell engager (150). Finally, a subset of lung NETs also shows high expression of hepatocyte nuclear factor 1α and 4α (HNF1a/HNF4a). In vitro such HNF+ lung NETs have been correlated with a response to FGFR3 and FGFR4 inhibitors (45).
There is evidence suggesting that most of these approved drugs also target the tumor microenvironment in addition to NET cells. SSAs, a cornerstone in NET management, bind to SSTRs overexpressed on many NET cells, inhibiting the release of various hormones and growth factors that alter the microenvironment. Therefore, SSAs are effective not only in controlling secretory syndrome induced by NETs, but also in stabilizing the disease with an antiproliferative effect (151). NETs often exhibit high levels of VEGF, which promotes angiogenesis by binding to VEGFR on endothelial cells. Key signaling pathways, such as PI3K/AKT/mTOR and Notch, regulate angiogenesis in NETs. Therefore, VEGF inhibitors, like sunitinib (131), cabozantinib (134), and bevacizumab, and mTOR inhibitors, like everolimus, which target angiogenesis, are potential treatments (130, 152). Understanding the molecular mechanisms of angiogenesis in NETs is essential for developing effective new drugs. Belzutifan is a novel agent targeting the HIF pathway, which plays a crucial role in cellular response to hypoxia and is dysregulated in some NETs (135). Though its current indication is in patients with VHL disease, belzutifan’s mechanism should also be identified outside the rare situation of germinal VHL disease.
Regarding other treatments targeting the microenvironment, immunotherapy using checkpoint inhibitors (CPIs) alone has been disappointing in NETs; this is especially due to its lack of predictive factors of response (153). Nevertheless, there is a subpopulation of NETs that show a signal of susceptible for response to CPI, enriched for atypical lung carcinoids (153–156). Tumor mutational burden and microsatellite instability–high status have been associated with better responses to CPI, as they may increase neoantigen load, enhancing immune recognition. However, they remain rare (<5%) in NETs (157). PD-L1 expression seems less predictive for response to CPI in NETs than in other cancers (153). Whether temozolomide-induced high tumor-mutation burden is a suitable molecular selection marker for immunotherapy in NETs is yet to be determined, as conflicting data have emerged from glioblastoma and colon cancers (158, 159). Combining immunotherapy with targeted therapies or peptide receptor radionuclide therapy may enhance its effectiveness and overcome resistance mechanisms, but the first results reporting this association were conflicting (160, 161).
The nomenclature gap across organs seems to be slowly closing, as in 2022, the WHO proposed a uniform NET nomenclature (1, 162). Still there is an urgent need for a clear, clinically relevant definition of high-grade NETs that incorporates current concepts of NEN evolution. Similarly, moving toward a morphomolecular classification is likely to be more relevant for patients with NETs. In parallel, gaining a better understanding of the cell of origin, evolutionary history, and the role of the microenvironment is essential for effectively tackling these diseases. The wide range of clinical behaviors displayed by NETs — from slow progression to metastatic — likely reflects an underappreciated heterogeneity in the NE cells of origin and their microenvironments. For example, pulmonary NETs can be found in different anatomic locations within the lung — peripheral or central (in the bronchi) — which are associated with distinct molecular features. It is very likely that these tumors arise from distinct NE cells depending on their location and exposure to different environmental stimuli. To accelerate the discovery of novel therapeutic targets, further comprehensive molecular studies are needed. Only such studies will guide researchers in testing biologically grounded therapeutic options that will improve patient outcomes and drive personalized therapy. The current trend toward agnostic oncology for rare cancers, which combines them based on single actionable molecular alterations or their common NE nature, is likely to fail without robust molecular basis.
We extend our gratitude to all rare cancers genomics initiative collaborators (https://rarecancersgenomics.com/) for their valuable contributions to the discussions, which greatly enriched the Review manuscript. Where authors are identified as personnel of the International Agency for Research on Cancer/WHO, the authors alone are responsible for the views expressed in this article and they do not necessarily represent the decisions, policy, or views of the International Agency for Research on Cancer/WHO.
Address correspondence to: Lynnette Fernandez-Cuesta, International Agency for Research on Cancer (IARC-WHO), 25 Av. Tony Garnier, 69007 Lyon, France. Email: fernandezcuestal@iarc.who.int.
Conflict of interest: CT is the secretary of the UK and Ireland Neuroendocrine Tumour Society and cochair of the Neuroendocrine Tumor Research Foundation Board of Scientific Advisors (USA). TW leads the scientific committee of the Groupe d’Etude des Tumeurs Endocrines and the European Neuroendocrine Tumor Society lung NET task force.
Copyright: © 2025, Fernandez-Cuesta et al. This is an open access article published under the terms of the Creative Commons Attribution 4.0 International License.
Reference information: J Clin Invest. 2025;135(5):e186702. https://doi.org/10.1172/JCI186702.