Journal List > J Pathol Transl Med > v.57(1) > 1516080847

Kim, Han, Kim, Jung, and Ryu: The proteomic landscape shows oncologic relevance in cystitis glandularis

Abstract

Background

The relationship between cystitis glandularis (CG) and bladder malignancy remains unclear.

Methods

We identified the oncologic significance of CG at the molecular level using liquid chromatography-tandem mass spectrometry-based proteomic analysis of 10 CG, 12 urothelial carcinoma (UC), and nine normal urothelium (NU) specimens. Differentially expressed proteins (DEPs) were identified based on an analysis of variance false discovery rate < 0.05, and their functional enrichment was analyzed using a network model, Gene Set Enrichment Analysis, and Gene Ontology annotation.

Results

We identified 9,890 proteins across all samples and 1,139 DEPs among the three entities. A substantial number of DEPs overlapped in CG/NU, distinct from UC. Interestingly, we found that a subset of DEP clusters (n = 53, 5%) was differentially expressed in NU but similarly between CG and UC. This “UC-like signature” was enriched for reactive oxygen species (ROS) and energy metabolism, growth and DNA repair, transport, motility, epithelial-mesenchymal transition, and cell survival. Using the top 10 shortlisted DEPs, including SOD2, PRKCD, CYCS, and HCLS1, we identified functional elements related to ROS metabolism, development, and transport using network analysis. The abundance of these four molecules in UC/CG than in NU was consistent with the oncologic functions in CG.

Conclusions

Using a proteomic approach, we identified a predominantly non-neoplastic landscape of CG, which was closer to NU than to UC. We also confirmed a small subset of common DEPs in UC and CG, suggesting that altered ROS metabolism might imply potential cancerous risks in CG.

Cystitis glandularis (CG), a common glandular lesion of the urinary bladder, is a cystic structure with a glandular lumen surrounded by outer urothelial cells that usually evolve from von Brunn’s nests after long-standing inflammation or irritation [1]. CG with intestinal metaplasia (IM) or the presence of mucin-containing goblet cells can produce abundant extracellular mucin that requires differentiation from primary adenocarcinoma, urachal adenocarcinoma, or metastatic adenocarcinoma from other organs [1,2]. CG can also develop a mass-like florid lesion that mimics bladder cancer on radiologic and cystoscopic examinations [3].
In addition to the clinical and histopathological disguise characteristics of CG, its relationship with malignancy has been discussed in studies with contrasting results [2]. For example, the co-occurrence of CG and adenocarcinoma has been frequently reported in the urinary bladder, which may indicate a connection between these two entities [1]. Furthermore, prior studies, using single gene or chromosome assays, underlaid the molecular basis supporting the premalignant nature of CG with or without IM by elucidating telomere shortening and chromosomal instability [4], p53 loss of heterozygosity and overexpression [5], and nuclear β-catenin expression [6]. In contrast, other studies failed to find a clear indication that CG, with or without IM, increased the future risk of developing bladder cancer during years of retrospective observation [7,8]. Comprehensive characterization would greatly help identify the pathobiology and clinical implications of CG.
In this study, we compared the proteomic landscape of CG with that of urothelial carcinoma (UC) and normal urothelium (NU). Differentially expressed proteins (DEPs) with functional enrichment were analyzed in CG, UC, and NU to identify the oncologic significance of CG at the molecular level.

MATERIALS AND METHODS

Patient selection

From the Seoul National University Hospital, 10 CG, 12 UC, and 9 NU specimens were collected after excluding patients who previously had any bladder tumor or intravesical treatment. Experienced pathologists (M.J. and H.S.R.) reviewed the diagnoses using hematoxylin-eosin slides based on the 2022 World Health Organization (WHO) classification [9]. All samples were obtained by transurethral resection of the urinary bladder, except for one NU, which was obtained from the ureter.

Liquid chromatography-tandem mass spectrometry

Formalin-fixed paraffin-embedded slides were scrapped and lysed in a sodium dodecyl sulfate-extraction buffer. After protein isolation, samples were sonicated and precipitated using acetone. Proteins were digested according to the filter-aided sample preparation procedure [10]. Tandem mass tag 6-plex labeling was employed for the peptide samples according to the manufacturer’s instructions, with modifications. After pooling the labeled peptides, the sample was separated into 12 fractions using Agilent 1290 bioinert high-pH reverse-phase liquid chromatography (Agilent, Santa Clara, CA, USA). For each peptide fraction, liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis was conducted using a Q-ExactivePlus mass spectrometer (Thermo Fisher, Waltham, MA) with an Ultimate 3000 RSLC system (Dionex, Sunnyvale, CA, USA), as previously described, with modifications [11]. Peptides were identified using a false discovery rate (FDR) of 1% as cutoff. Raw data are available in the ProteomeXchange Consortium (PXD027602) [12,13].

Proteomic data analysis

The MS raw data were processed using MaxQuant ver. 1.5.3.1 (Max Planck Institute of Biochemistry, Munich, Germany) [14] with the Andromeda engine [15]. For label-free quantification, the iBAQ algorithm was used on MaxQuant platform [16]. Using normalized abundance, we identified DEPs among CG, UC, and NU based on an analysis of variance (ANOVA) [14] test as a cutoff for permutation-based FDR < 0.05. Subsequently, the DEPs were hierarchically clustered based on Euclidean distance. Following a review of their abundance, we enlisted “UC-like” DEP clusters as “UC-like signature”, for which CG and UC overlapped but NU did not. All analyses were conducted using Perseus ver. 1.6.14.0 (Max Planck Institute of Biochemistry).

Functional enrichment analysis

To determine the biological processes represented by the lists of DEPs, we investigated protein-protein interactions (PPIs) using String database [17], Gene Ontology-biologic process (GOBP) annotation using Toppgene Suite [18], and the Molecular Signatures Database (MSigDB) Hallmark gene sets using pre-ranked Gene Set Enrichment Analysis (GSEA) [19]. To examine the associations among GOBPs, we reconstructed a network model using REVIGO [20]. Network models of PPI and GOBP were visualized using Cytoscape ver. 3.7.2 [21]. The z-scores of the selected proteins were compared between the group with and without IM using a two-tailed t-test.

RESULTS

We investigated 31 urothelial specimens, including 10 CGs, 12 UCs, and nine NUs, using a proteomic approach. A schematic outline of this study is shown in Fig. 1. The median ages were 40, 70, and 63.5 years, respectively, and the male-to-female ratios were 1.0, 0.4, and 0.9, respectively. Half of the specimens that contained CG were accompanied by IM. All UCs were confined to the mucosa (stage Ta), and eight (67%) were of WHO high grade.
From LC-MS/MS analysis, we identified 9,890 proteins across all samples and 1,139 DEPs among the three entities (ANOVA FDR < 0.05). With unsupervised hierarchical clustering, a substantial number of DEPs were discovered to overlap with CG/NU in a way distinct from UC (Fig. 2). CG was not differentiated by the presence of IM (Fig. 2, asterisks). Interestingly, we found a subset of DEP clusters (n = 53, 5%) that were differentially expressed in NU but similarly expressed between CG and UC (Fig. 2, arrows); thus, these proteins were named as “UC-like signature” (Fig. 3A). Since the signature were similar within each group, the patterns of these signatures were thought to represent the general profile of CG.
We hypothesized that this signature might represent the UC-like characteristics of the CG. Using PPI analysis, we explored cellular processes represented by the “UC-like signature”. This resulted in sets of DEP networks within the “UC-like signature”, which were enriched for reactive oxygen species (ROS) and energy metabolism [2224], growth and DNA repair [25,26], transport [27,28], motility and epithelial-mesenchymal transition (EMT) [29,30], and cell survival [31,32] (Fig. 3B). The unit associated with ROS and energy metabolism, which was connected to the one with growth and DNA repair functions, included the top-ranked DEPs (e.g., SOD2, PRKCD, CYCS, and MRPL23) (Fig. 3B). Using the official pre-ranked GSEA of the “UC-like signature”, we also identified significant enrichment of ultraviolet (UV) response (normalized enrichment score = 1.772, FDR = 0.03) for SOD2, PRKCD, and MRPL23 (Fig. 3C). UV radiation stimulates carcinogenic sequences by promoting ROS generation and DNA damage [33]. These data suggest that CG is predominantly a benign lesion, but oxidative stress might be a gateway theme representing oncogenic potential in CG.
To further characterize the ”UC-like signature”, we shortlisted the top 10 most significant DEPs (Fig. 3A, asterisk), followed by GOBP analysis (Fig. 4A). Subsequently, we constructed a network of connections between the GOBP terms. Consistent with the aforementioned single-protein networks, functional elements related to ROS metabolism, development, and transport were highlighted by GOBPs (Fig. 4B), which suggests that these functions could collectively substantiate the cancerous risks in CG.
Of note, GOBPs were represented almost exclusively by four molecules, i.e., SOD2, PRKCD, CYCS, and HCLS1 (Fig. 4A), and their abundance in UC/CG compared to NU accordantly pointed toward oncologic functions. For example, previous studies have shown an association of decreased level of SOD2 and increased level of PRKCD, as observed in CG/UC versus NU, with ROS production [22,24]. Overt oxidative stress by activating oncogenic signaling pathways, DNA mutations, EMT, and stromal remodeling promotes bladder cancer development and progression [22]. In line with this, altered expression of CYCS and HCLS1 might deregulate oxidative respiration and the RAS signaling pathway, thereby enhancing the malignant behavior of bladder cancer [23,34]. We further examined whether these four statistically and functionally significant molecules were differentially expressed with IM and found no significant differences related to IM (Fig. 4C).

DISCUSSION

The malignant risk of CG and its association with bladder cancer have been disputed [1,2]. Considering the high prevalence of CG, there is a dire need to clarify its oncologic implications. Using LC-MS/MS, we revealed the proteomic landscape of CG in relation to malignant (UC) and normal (NU) urothelial tissues. CG generally presented an NU-like profile, in contrast to UC, indicating the overall benign nature of CG. By retrospectively observing a handful of patients with CG, consistent with this result, previous studies failed to confirm clear association between UC and CG [7,8]. Conversely, we confirmed that some of the clustered DEPs showed an opposite pattern across the diagnosis; these proteins were shared by CG and UC but not by NU. Using functional and network analyses, we found that these DEPs coded for ROS and energy metabolism, growth and DNA repair, transport, motility and EMT, and cell survival, further supporting the relation of this “UC-like signature” with the malignant risk in CG. To our knowledge, this is the first study to determine the oncogenic characteristics of CG using a comprehensive proteomic analysis. Further large-scale studies are required to determine the risks of UC development in patients with CG.
Using an interconnected functional network analysis of the shortlisted proteins, including SOD2, PRKCD, CYCS, and HCLS1, ROS metabolism, structural development, and transport functions were concordantly found to be significantly enriched in the “UC-like signature.” ROS are tightly regulated in normal cellular environments, and increased ROS act as versatile transducers for the generation and progression of bladder cancer [22]. For example, both UC and CG showed downregulation of SOD2, a scavenging molecule, and upregulation of PRKCD, an oxidative enzyme associated with ROS accumulation. ROS imbalance induces oxidative DNA damage, mutations, propagation of oncogenic signals including RAS, mitogen-activated protein kinases, phosphoinositide 3-kinase, and nuclear factor κB pathways, EMT, and stromal modification in bladder cancer [22,24]. Consistent with this, a previous study suggested that alterations in ROS metabolism participate in inflammation-associated cancer sequences [35,36]. In addition, upregulation of PRKCD can further support a malignant phenotype in bladder cancer by promoting migration and invasion [37]. CYCS is instrumental in regulating oxidative phosphorylation and energy metabolism, and high maintenance of CYCS, as identified in UC/CG, could induce metabolic modification in bladder cancer, or Warburg effect [23]. In addition, UC and CG showed lower HCLS1 levels than NU. HCLS1, an actin-binding molecule supporting cellular transport and trafficking, has been shown to result in adverse outcomes in bladder cancer [34,38]. We believe that the “UC-like signature” might reflect the oncogenic pathway related to altered ROS and energy metabolism and structural modification in CG, and the expression of these signature proteins deserves further study to determine the risk of aggression in CG.
This study lacked specimens with adenocarcinoma in the analysis. The malignant association of CG with IM has been conjectured, typically in terms of bladder adenocarcinoma, but definite evidence has not been confirmed [1]. Instead, altered ROS and energy metabolism have been implicated in the transformation of dysplasia to adenocarcinoma in the esophagus, consistent with the functions enriched in the “UC-like signature” [35]. In addition, the overall proteomic profile and biomarkers relevant to such functions showed no difference regardless of the presence of IM in CG. However, previous studies have suggested that CG with IM may be more advanced than CG without IM regarding tumorigenic potential, as exemplified by more robust telomere shortening or β-catenin activation [4,6]. Therefore, we reasonably speculate that the “UC-like signature” might reflect the malignancy risk encoded in CG irrespective of IM at a global level.
In conclusion, using comprehensive proteomic profiling, we identified a predominantly non-neoplastic landscape of CG that is closer to NU than to UC. Furthermore, we confirmed a small subset of common DEPs in UC and CG, suggesting altered functions of ROS metabolism that might imply potential cancerous risks in CG.

Notes

Ethics Statement

All procedures performed in the current study were approved by the Institutional Review Board (IRB) (H-2009-163-1160; 05 October 2020) per the 1964 Helsinki Declaration and its later amendments. Formal written informed consent was not required, with a waiver by the appropriate IRB.

Availability of Data and Material

The datasets generated or analyzed during the current study are available in the ProteomeXchange Consortium (http://www.proteomexchange.org/) via the PRIDE database (accession: PXD027602).

Code Availability

Not applicable.

Author Contributions

Conceptualization: HSR. Data curation: HSR. Formal analysis: JYK, MJ, HSR. Funding acquisition: MJ. Investigation: JYK, MJ, HSR. Methodology: HSR. Project administration: HSR. Resources: HSR. Software: DH, HK, HSR. Supervision: HSR. Validation: DH, HK, MJ, HSR. Visualization: JYK, MJ, HSR. Writing—original draft: JYK, MJ. Writing—review and editing: MJ, HSR.

Conflicts of Interest

The authors declare that they have no potential conflicts of interest.

Funding Statement

This research was supported by the Korean Society of Pathologists Grant (No. KSPG2020-04).

References

1. Williamson SR, Lopez-Beltran A, Montironi R, Cheng L. Glandular lesions of the urinary bladder: clinical significance and differential diagnosis. Histopathology. 2011; 58:811–34.
2. Samaratunga H, Martignoni G, Egevad L, Delahunt B. Premalignant lesions of the urinary bladder. Pathology. 2013; 45:243–50.
crossref
3. Song YS, Jang KS, Jang SH, et al. The intestinal type of florid cystitis glandularis mimics bladder tumor: a case report. Korean J Pathol. 2007; 41:116–8.
4. Morton MJ, Zhang S, Lopez-Beltran A, et al. Telomere shortening and chromosomal abnormalities in intestinal metaplasia of the urinary bladder. Clin Cancer Res. 2007; 13:6232–6.
crossref
5. Srivastava S, Yan B, Chin SY, Muliana T, Salto-Tellez M, Teh M. Nuclear p53 expression is associated with allelic imbalance (TP53) in glandular dysplasia and typical cystitis glandularis: a LCM-based molecular analysis. Clin Genitourin Cancer. 2012; 10:57–9.
crossref
6. Bryan RT, Nicholls JH, Harrison RF, Jankowski JA, Wallace DM. The role of beta-catenin signaling in the malignant potential of cystitis glandularis. J Urol. 2003; 170:1892–6.
7. Corica FA, Husmann DA, Churchill BM, et al. Intestinal metaplasia is not a strong risk factor for bladder cancer: study of 53 cases with long-term follow-up. Urology. 1997; 50:427–31.
crossref
8. Smith AK, Hansel DE, Jones JS. Role of cystitis cystica et glandularis and intestinal metaplasia in development of bladder carcinoma. Urology. 2008; 71:915–8.
crossref
9. WHO Classification of Tumours Editorial Board. Urinary and male genital tumours. Lyon: International Agency for Research on Cancer;2022. p. 131–92.
10. Wisniewski JR, Zougman A, Nagaraj N, Mann M. Universal sample preparation method for proteome analysis. Nat Methods. 2009; 6:359–62.
crossref
11. Han D, Moon S, Kim Y, Kim J, Jin J, Kim Y. In-depth proteomic analysis of mouse microglia using a combination of FASP and StageTip-based, high pH, reversed-phase fractionation. Proteomics. 2013; 13:2984–8.
crossref
12. Jung M, Lee C, Han D, et al. Proteomic-based machine learning analysis reveals PYGB as a novel immunohistochemical biomarker to distinguish inverted urothelial papilloma from low-grade papillary urothelial carcinoma with inverted growth. Front Oncol. 2022; 12:841398.
crossref
13. Kim B, Jung M, Moon KC, et al. Quantitative proteomics identifies TUBB6 as a biomarker of muscle-invasion and poor prognosis in bladder cancer. Int J Cancer. 2022. Sep. 10. [Epub]. https://doi.org/10.1002/ijc.34265 .
crossref
14. Tyanova S, Temu T, Cox J. The MaxQuant computational platform for mass spectrometry-based shotgun proteomics. Nat Protoc. 2016; 11:2301–19.
crossref
15. Cox J, Neuhauser N, Michalski A, Scheltema RA, Olsen JV, Mann M. Andromeda: a peptide search engine integrated into the MaxQuant environment. J Proteome Res. 2011; 10:1794–805.
crossref
16. Schwanhausser B, Busse D, Li N, et al. Global quantification of mammalian gene expression control. Nature. 2011; 473:337–42.
crossref
17. Szklarczyk D, Gable AL, Lyon D, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019; 47:D607–13.
crossref
18. Chen J, Bardes EE, Aronow BJ, Jegga AG. ToppGene Suite for gene list enrichment analysis and candidate gene prioritization. Nucleic Acids Res. 2009; 37:W305–311.
crossref
19. Liberzon A, Birger C, Thorvaldsdottir H, Ghandi M, Mesirov JP, Tamayo P. The Molecular Signatures Database (MSigDB) hallmark gene set collection. Cell Syst. 2015; 1:417–25.
20. Supek F, Bosnjak M, Skunca N, Smuc T. REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS One. 2011; 6:e21800.
crossref
21. Shannon P, Markiel A, Ozier O, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003; 13:2498–504.
crossref
22. Mendes F, Pereira E, Martins D, et al. Oxidative stress in bladder cancer: an ally or an enemy? Mol Biol Rep. 2021; 48:2791–802.
crossref
23. Herrmann PC, Herrmann EC. Oxygen metabolism and a potential role for cytochrome c oxidase in the Warburg effect. J Bioenerg Biomembr. 2007; 39:247–50.
crossref
24. Talior I, Tennenbaum T, Kuroki T, Eldar-Finkelman H. PKC-delta-dependent activation of oxidative stress in adipocytes of obese and insulin-resistant mice: role for NADPH oxidase. Am J Physiol Endocrinol Metab. 2005; 288:E405–11.
25. Bywater MJ, Pearson RB, McArthur GA, Hannan RD. Dysregulation of the basal RNA polymerase transcription apparatus in cancer. Nat Rev Cancer. 2013; 13:299–314.
crossref
26. Yuniati L, Lauriola A, Gerritsen M, et al. Ubiquitylation of the ER-shaping protein lunapark via the CRL3 (KLHL12) ubiquitin ligase complex. Cell Rep. 2020; 31:107664.
27. Chukkapalli S, Amessou M, Dekhil H, et al. Ehd3, a regulator of vesicular trafficking, is silenced in gliomas and functions as a tumor suppressor by controlling cell cycle arrest and apoptosis. Carcinogenesis. 2014; 35:877–85.
crossref
28. Nakamura N, Yamamoto A, Wada Y, Futai M. Syntaxin 7 mediates endocytic trafficking to late endosomes. J Biol Chem. 2000; 275:6523–9.
crossref
29. Wang G, Li Y, Li J, et al. microRNA-199a-5p suppresses glioma progression by inhibiting MAGT1. J Cell Biochem. 2019; 120:15248–54.
crossref
30. Moamer A, Hachim IY, Binothman N, Wang N, Lebrun JJ, Ali S. A role for kinesin-1 subunits KIF5B/KLC1 in regulating epithelial mesenchymal plasticity in breast tumorigenesis. EBioMedicine. 2019; 45:92–107.
crossref
31. Sharma A, Almasan A. USP14 regulates DNA damage response and is a target for radiosensitization in non-small cell lung cancer. Int J Mol Sci. 2020; 21:6383.
crossref
32. Wang L, Li Z, Sievert D, et al. Loss of NARS1 impairs progenitor proliferation in cortical brain organoids and leads to microcephaly. Nat Commun. 2020; 11:4038.
crossref
33. Bosch R, Philips N, Suarez-Perez JA, et al. Mechanisms of photoaging and cutaneous photocarcinogenesis, and photoprotective strategies with phytochemicals. Antioxidants (Basel). 2015; 4:248–68.
crossref
34. Tikoo A, Shakri R, Connolly L, et al. Treatment of ras-induced cancers by the F-actin-bundling drug MKT-077. Cancer J. 2000; 6:162–8.
35. O’Farrell NJ, Phelan JJ, Feighery R, et al. Differential expression profiles of oxidative stress levels, 8-oxo-dG and 4-HNE, in Barrett’s esophagus compared to esophageal adenocarcinoma. Int J Mol Sci. 2019; 20:4449.
crossref
36. Wigner P, Grebowski R, Bijak M, Saluk-Bijak J, Szemraj J. The interplay between oxidative stress, inflammation and angiogenesis in bladder cancer development. Int J Mol Sci. 2021; 22:4483.
crossref
37. Garg R, Benedetti LG, Abera MB, Wang H, Abba M, Kazanietz MG. Protein kinase C and cancer: what we know and what we do not. Oncogene. 2014; 33:5225–37.
crossref
38. Sanchez-Carbayo M, Socci ND, Lozano J, Saint F, Cordon-Cardo C. Defining molecular profiles of poor outcome in patients with invasive bladder cancer using oligonucleotide microarrays. J Clin Oncol. 2006; 24:778–89.
crossref

Fig. 1
Schematic outline of the study. (Top) Liquid chromatography-tandem mass spectrometry (LC-MS/MS)–based proteomic analysis of 12 urothelial carcinoma (UC),10 cystitis glandularis (CGs), and nine normal urothelial (NU) specimens identified 9,890 proteins. Formalin-fixed paraffin-embedded (FFPE) slides were scrapped and proteins were isolated. After sonication and precipitation of samples, proteins were digested according to the filter-aided sample preparation (FASP) procedure. (Bottom) Differentially expressed proteins (DEPs) were discovered based on analysis of variance (false discovery rate < 0.05), and DEPs with UC-like signatures were revealed in CG. Functional enrichment analysis was performed for the DEPs and top 10 DEPs, respectively.
jptm-2022-10-24f1.gif
Fig. 2
Unsupervised hierarchical clustering of 1,139 differentially expressed proteins identified by an analysis of variance test of urothelial carcinoma (UC), cystitis glandularis (CG), and normal urothelium (NU). Arrows indicate a “UC-like signature” that showed disparity between UC/CG and NU. Asterisks indicate CG samples showing intestinal metaplasia.
jptm-2022-10-24f2.gif
Fig. 3
Pathobiologic characteristics of the “urothelial carcinoma (UC)-like signature”. (A) The UC-like signature proteins. Asterisks denote the 10 top-listed proteins by false discovery rate (FDR). (B) Protein-protein interaction networks of the “UC-like signature” proteins with their common functions. (C) Gene Set Enrichment Analysis of the “UC-like signature” shows enrichment of response to ultraviolet in molecular hallmark function. CG, cystitis glandularis; EMT, epithelial-mesenchymal transition; FWER, family wise error rate; NES, normalized enrichment score; NU, normal urothelium; ROS, reactive oxygen species.
jptm-2022-10-24f3.gif
Fig. 4
Imperative functions enriched in the “urothelial carcinoma-like signature”. (A) Gene Ontology-biologic processes (GOBPs) represented by the top 10 shortlisted proteins of “urothelial carcinoma-like signature”. (B) Network analysis of the GOBPs identifies reactive oxygen species (ROS) metabolism, structure development, and transport as common functional themes. (C) SOD2, PRKCD, CYCS, and HCLS1 were comparable between cystitis glandularis with and without intestinal metaplasia (IM). The value is expressed with z-scores of the proteome abundance. FDR, false discovery rate.
jptm-2022-10-24f4.gif
TOOLS
Similar articles