Next Article in Journal
Unveiling Niaprazine’s Potential: Behavioral Insights into a Re-Emerging Anxiolytic Agent
Previous Article in Journal
Polyunsaturated Fatty Acids as Potential Treatments for COVID-19-Induced Anosmia
Previous Article in Special Issue
Prognosis Prediction in Head and Neck Squamous Cell Carcinoma by Radiomics and Clinical Information
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evolutionary Mechanism Based Conserved Gene Expression Biclustering Module Analysis for Breast Cancer Genomics

School of Biomedical Engineering, Guangzhou Medical University, Guangzhou 511436, China
*
Authors to whom correspondence should be addressed.
Submission received: 23 June 2024 / Revised: 23 August 2024 / Accepted: 2 September 2024 / Published: 12 September 2024
(This article belongs to the Special Issue Bioinformatics: From Methods to Applications)

Abstract

:
The identification of significant gene biclusters with particular expression patterns and the elucidation of functionally related genes within gene expression data has become a critical concern due to the vast amount of gene expression data generated by RNA sequencing technology. In this paper, a Conserved Gene Expression Module based on Genetic Algorithm (CGEMGA) is proposed. Breast cancer data from the TCGA database is used as the subject of this study. The p-values from Fisher’s exact test are used as evaluation metrics to demonstrate the significance of different algorithms, including the Cheng and Church algorithm, CGEM algorithm, etc. In addition, the F-test is used to investigate the difference between our method and the CGEM algorithm. The computational cost of the different algorithms is further investigated by calculating the running time of each algorithm. Finally, the established driver genes and cancer-related pathways are used to validate the process. The results of 10 independent runs demonstrate that CGEMGA has a superior average p-value of 1.54 × 10−4 ± 3.06 × 10−5 compared to all other algorithms. Furthermore, our approach exhibits consistent performance across all methods. The F-test yields a p-value of 0.039, indicating a significant difference between our approach and the CGEM. Computational cost statistics also demonstrate that our approach has a significantly shorter average runtime of 5.22 × 100 ± 1.65 × 10−1 s compared to the other algorithms. Enrichment analysis indicates that the genes in our approach are significantly enriched for driver genes. Our algorithm is fast and robust, efficiently extracting co-expressed genes and associated co-expression condition biclusters from RNA-seq data.

1. Introduction

The quest to understand cancer at the molecular level reveals a complex landscape where cancer cells demonstrate unique gene expression patterns diverging significantly from their healthy counterparts. The advent of high-throughput sequencing technologies, such as RNA sequencing (RNA-seq), has opened new vistas in cancer research by furnishing an unparalleled richness of transcriptomic data [1,2]. These technologies have catapulted us into an era where the voluminous and sophisticated nature of genomic data presents a formidable challenge, necessitating the use of computational techniques to sift through this extensive dataset to unearth genes that exhibit coordinated expression under specific conditions [3,4]. Concurrently, it has become increasingly clear that cancer is not a consequence of anomalies in single genes but emerges from complex interactions among multiple co-expressed RNAs.
Amidst this complexity, it is crucial to identify the specific conditions under which cancer genes are co-expressed. The specific conditions of co-expression, which leads to the research question of biclustering. Biclustering allows us to identify clusters of genes that exhibit coordinated expression under certain conditions, offering insights that conventional clustering methods, which assume co-expression across all conditions, might miss [5,6,7,8,9,10,11,12]. The quest for deciphering these patterns has spurred the development of various heuristic algorithms since the pioneering CC algorithm by Cheng and Church in 2000 [13]. These include the Modified Cheng and Church’s algorithm (MCC) [14], the Large Average Submatrix algorithm (LAS) [15], the Relative Density based Biclustering Method (RelDenClu) [16], and the Connectedness-based subspace clustering (CBSC) [17], all designed to tackle the biclustering challenge as an optimization problem, striving for efficient extraction of meaningful genomic insights from RNA-seq data [18,19]. In the comparative analysis of biclustering algorithms, each method presents distinct characteristics that allow it to be applied to specific research needs in a tailored manner. The spectral CC method is particularly effective in identifying co-expressed gene and condition submatrices, which is a valuable capability in the analysis of gene expression data. However, its high time complexity and sensitivity to noise can act as limitations, particularly when dealing with large datasets. MCC offers an ensemble learning approach that enhances stability and accuracy across various data types, with robustness against noise, although this is at the cost of higher computational requirements. LAS proposes a statistical method that is effective in detecting numerically significant submatrices within high-dimensional data, demonstrating strong resistance to noise and scalability to large datasets, despite its computational intensity. Conserved Gene Expression Motif (CGEM) employs a graph-based approach to identify conserved gene expression motifs, offering insights into gene regulatory networks [20]. However, its utility is more suited to medium to small-scale datasets, and its effectiveness depends on the stability of the graph structure. RelDenClu and CBSC are density-based subspace clustering algorithms that aim to partition data points into multiple subspaces, wherein data points within the same subspace exhibit higher connectivity compared to those in different subspaces. RelDenClu identifies subsets of observations exhibiting dependence between features by comparing joint and marginal densities, and subsequently groups data points based on these features. In contrast, CBSC leverages connectivity scores to identify subspaces. Both algorithms exhibit robustness, although CBSC is more computationally intensive. RelDenClu is capable of uncovering feature relationships based on nonlinear dependencies, whereas CBSC is more suited for linear relationships. In conclusion, the selection of an appropriate biclustering algorithm should be based on the characteristics of the data set, the specific research objectives, and the available computational resources. Each algorithm possesses distinctive advantages in identifying the intrinsic patterns within biological data. Table 1 provides a comparative analysis of the aforementioned six algorithms.
Although these algorithms represent a significant advance in the field, they are not without limitations. In particular, their reliance on exhaustive traversals to identify genetic traits associated with cancer renders them somewhat inefficient. This process is both time-consuming and a significant consumer of computing resources [21,22,23]. Furthermore, the complex interrelationships between diverse gene combinations and cancer progression present a significant challenge to optimization, particularly when using traditional methods [24]. However, recent advances have demonstrated the potential of evolutionary algorithms (EAs) in overcoming these challenges. They exhibit robust global optimization capabilities, parallel processing properties, adaptability and resilience, maintenance of population diversity, adaptive tuning of parameters, scalability of algorithms, and a search process that does not depend on predefined thresholds. These characteristics assist the genetic algorithm in circumventing local optimal solutions, thereby enhancing the algorithm’s efficacy in processing extensive gene expression data and offering a new paradigm for efficient bicluster identification.
Leveraging this insight, we introduce the Conserved Gene Expression Module algorithm enhanced by a Genetic Algorithm (CGEMGA). This novel approach not only capitalizes on the strengths of the CGEM algorithm but also harnesses the power of EAs to navigate the optimization landscape more effectively. By incorporating the Mean Squared Residual Score (MSR) criterion, CGEMGA stands out in its ability to efficiently identify the most optimal gene combinations pertinent to breast cancer, as validated by rigorous testing against breast cancer datasets, and Fisher’s exact test comparisons with other algorithms. Remarkably, this approach not only enhances reliability but also drastically reduces computational overhead, setting a new benchmark in biclustering algorithm efficiency.

2. Materials and Methods

2.1. Data Acquisition of Breast Cancer Samples

In our quest to unravel the complex genetic architecture of breast cancer, we have procured gene expression profiles from The Cancer Genome Atlas (TCGA). This global repository offers an extensive compendium of human cancer genomes, serving as a beacon for researchers worldwide [25,26,27,28]. Our analysis is anchored in the rich dataset of 421 breast cancer samples, encompassing a diverse array of 12,129 genes, carefully selected for their relevance to our study.
To further reinforce the empirical foundation of our study, we refer to the Catalogue of Somatic Mutations in Cancer (COSMIC) as a reference point. The Cancer Gene Census (CGC) on the Sanger Institute’s website provides access to the comprehensive and meticulously curated catalogue of cancer-driving genes, COSMIC. This database, which is instrumental in both cancer genetics research and drug development initiatives, has recently been expanded to encompass 739 genes that are crucial to the progression of cancer [29,30]. The judicious choice of these datasets shows how seriously we’re taking the molecular intricacies of breast cancer.

2.2. Overview of Biclustering

The concept of biclustering, a term coined by Cheng and Church in the groundbreaking introduction of their CC algorithm in 2000, has evolved significantly over the years. This pioneering contribution heralded a new era in the exploration of gene expression dynamics across diverse conditions through the application of biclustering techniques. These sophisticated algorithms, by framing the biclustering challenge as an optimization problem (COP), set out to uncover patterns within the genetic matrix that elude traditional analysis methods. At the heart of this endeavor is the quest for submatrices within the gene expression matrix that exhibit unique patterns of interest, encapsulated in the specially defined score functions and heuristic solutions of the COP [31,32].
To illustrate, let’s consider the data matrix X = { E , F } , where E = { e 1 , e N } symbolizes the set of N genes and F = { f 1 , , f M } represents the set of M conditions. Within this framework, a bicluster B emerges as a subset, defined by B = { E B , F B ; E B E ,   F B F } . This subset transcends to the status of a bicluster only if it adheres to specific patterns, typically quantified by the Mean Squared Residual (MSR) score.
The MSR is defined as:
H I , J = 1 I J i I , j J a i j a i J a I j + a I J 2 ,
where a i j denotes the matrix element, with the row and column means and the overall mean of the submatrix B calculated as follows:
a i J = 1 J j J a i j , a I j = 1 I i I a i j , and   a I J = 1 I [ J ] i I , j J a i j = 1 I i I a i J = 1 J j J a I j .
This mathematical rigor provides a robust framework to unravel the complex patterns of gene expression specific to breast cancer, laying the groundwork for a deeper understanding of the disease’s genetic architecture.

2.3. Conserved Biclustering Algorithm

At the heart of our exploration into the genetic underpinnings of breast cancer lies the Conserved Gene Expression Motif (CGEM), an innovative biclustering algorithm. This method is grounded in a simple yet profound principle: if a gene maintains consistent expression across a subset of samples, its expression level is considered conserved within that specific subset [33]. This approach diverges markedly from traditional algorithms such as CC, which rely heavily on scoring schemes to identify significant patterns. Instead, CGEM seeks out conserved gene modules across the entirety of the dataset, employing constraints to uncover these vital connections [34].
The CGEM algorithm distinguishes itself by its ability to unearth the largest xmotif, representing a conserved gene expression pattern of paramount significance. This endeavor involves a meticulous, iterative process composed of three fundamental steps. Initially, the algorithm identifies the most extensive xmotif within the dataset, guided by a designated seed s. This identification process prioritizes the discovery of motifs with a substantial number of conserved genes. Following this, it strategically removes the samples aligning with this motif from consideration, thereby refining the dataset. The subsequent step involves a renewed search within this pared-down dataset, aiming to locate the next largest motif. This iterative loop continues until the algorithm selects the submatrix with the largest row-column size among all potential biclusters.
This iterative, constraint-based methodology sets CGEM apart, offering a novel lens through which to view the complex landscape of gene expression in breast cancer. By focusing on the conservation of gene expression across samples, CGEM provides invaluable insights into the genetic consistency that may underpin this disease, offering promising avenues for further research and potential therapeutic targets.

2.4. Evolutionary Mechanism-Based Conserved Gene Expression Biclustering Module

While the CGEM algorithm marks a significant leap forward in identifying conserved gene modules, it harbors certain limitations that hinder its potential for global optimization. These challenges include a predefined iteration count lacking mathematical justification, an absence of a fitness function for optimal module evaluation, and the arbitrary selection of the initial seed s , which compromises the pursuit of the global optimum [35]. Additionally, despite the high accuracy of many biclustering algorithms, their extensive exploratory capabilities contribute to increased computational demands [21,23].
To surmount these obstacles, we introduce the CGEMGA (Conserved Gene Expression Module Genetic Algorithm), an innovative approach that incorporates an evolutionary algorithm to navigate toward globally optimal biclustering. The genetic algorithm (GA), inspired by Darwin’s principle of natural selection, serves as the cornerstone of our method, utilizing population-based strategies and the survival of the fittest principle to refine the search for optimal gene modules [36]. Unlike the CGEM algorithm, which relies on randomly selected seed s , our approach optimizes the initial seed selection through GA mechanisms, including population genetics, crossover, and mutation. This strategy enables the identification of an optimal seed s by evaluating the MSR of submatrices generated from each seed and selecting the one with the minimal MSR value, thereby ensuring the selection of the most stable module.
Outlined in Table 2, our algorithm’s process begins with the input of genes, samples, and their expression values, alongside intervals representing gene expression states. Assuming distinct intervals for each gene’s states, the algorithm uniformly selects n s sample groups from the entire sample set. In contrast to the CGEM algorithm’s random seed generation, our method employs GA to ascertain optimal seeds through a meticulously defined procedure:
1.
Initiate with a population of n chromosomes C i ( i = 1,2 , , n ) as potential seeds s .
2.
For each chromosome C i , identify a subset of samples D i of size s d .
3.
Include gene-sample pairs ( g , s ) in set G i j if gene g exists in state s across all samples in D i , and also incorporate samples matching c across all gene states in G i j .
4.
Compute the MSR fitness value for each chromosome.
5.
Apply GA’s selection, mutation, and crossover operations to optimize based on the MSR fitness value, thereby deriving the optimal solution.
6.
Exclude any C i j representing less than a fraction α of the samples.
7.
Select the module with the lowest MSR from all C i as the final choice.
Employing GA to refine the initial seed s and evaluating modules via MSR during each iteration not only circumvents the fixed iteration constraint of the CGEM algorithm but also significantly reduces computational time, showcasing the efficacy of our evolutionary approach. Figure 1 illustrates the flowchart of our algorithm.

2.5. Evaluation Metrics

In the complex process of validating our evolutionary mechanism-based conserved gene expression biclustering module, precision in detecting breast cancer genes stands paramount. To achieve this, we juxtapose our method against both the CGEM algorithm and other prevailing algorithms through a series of ten meticulously conducted independent runs. Post each execution, the generated results undergo a rigorous comparison with the breast cancer gene data cataloged in COSMIC. This comparative analysis serves not only as a validation of our method’s effectiveness but also as a critical link that bridges theoretical advancement with empirical confirmation.
The cornerstone of our evaluation lies in the application of Fisher’s exact test, a statistical method renowned for its precision in assessing the association between the outcomes of our algorithm and the presence of breast cancer genes. The Fisher exact test embarks on this task by meticulously calculating p values, predicated on the analysis of all conceivable configurations of 2 × 2 contingency tables that manifest marginal totals equivalent to or surpassing those observed. The test’s foundation is built upon a sample drawn from a population of size N , with m objects exhibiting trait A ( a within the sample and c outside the sample) and n objects not exhibiting trait A ( b within the sample and d outside the sample). Here, the aggregates of a and b form r , while those of c and d constitute s . The calculation of the p-value unfurls through the following formula:
p = m ! n ! r ! s ! a ! b ! c ! d ! N ! .
Delving deeper into the specifics, our application of Fisher’s exact test involves the enumeration of genes identified by the biclustering algorithm and present in the CGC dataset as a , juxtaposed with genes identified by the algorithm but absent in the CGC dataset as b . Concurrently, we account for genes present in the CGC dataset yet overlooked by the biclustering algorithm as c , and genes neither detected by the algorithm nor listed in the CGC dataset as d . A threshold of p < 0.05 delineates the boundary for statistical significance, guiding us in discerning meaningful associations.
The essence of our study is encapsulated within the realm of unsupervised learning models, which inherently operate without the crutch of predefined labels, thus obviating the hurdles of data partitioning. This attribute, while liberating, mandates the necessity for external validation through third-party corroborations. In our quest for validation, we anchor our trust in established cancer driver genes and related pathways, leveraging these benchmarks not just as a means of validation but as a beacon guiding our exploratory voyage through the genomic landscape of breast cancer.

3. Results

3.1. Experiment Setup

To examine the effectiveness of biclustering data searching of our CGEMGA algorithm, we conducted several experiments to analyze our approach from multiple perspectives: (1) investigating the superiority between our approach and CGEM algorithms by comparing the p-values of Fisher’s test with the MSR values [37,38,39]; (2) comparing our approach with widely-used existing biclustering approaches; (3) testing the computational cost of the entire algorithm’s runtime; (4) studying the quantitative differences between our approach and the CGEM algorithm by the F-test [40]; (5) unscrambling gene functions by enrichment analysis.
The experimental apparatus was manufactured by Lenovo in Shenzhen, China, and was powered by an Intel Xeon(R) CPU E3-1225 v6 @ 3.30 GHz × 2 processors, backed by 32 GB of RAM. The datasets comprised 421 breast cancer samples, intricately mapped out across 12,129 genes. Leveraging unsupervised methods, we plunged directly into the data, eschewing any division between training, testing, or validation datasets, thus ensuring an unadulterated analysis.

3.2. Ablation Study

3.2.1. Evolutionary Effect

To showcase our method’s superiority, we pitted it against the CGEM algorithm in ten independent trials, comparing their statistical significance and MSR values. Table 3 reveals a striking contrast: while the CGEM algorithm’s p-values fluctuate significantly, ours remain remarkably stable, underscoring the robustness of our approach. Specifically, CGEM’s p-values span from 9.62 × 10−6 to 9.71 × 10−2 averaging at 1.18 × 10−2 ± 3.01 × 10−2. The MSR variations mirror this volatility. In stark contrast, our method consistently delivers p values with minimal variance, showcasing not only our algorithm’s precision but also its reliability, a testament to the evolutionary mechanics at its core.

3.2.2. Stability Analysis

Venturing further, we employed the F-test to quantitatively assess the differences between our algorithm and CGEM, seeking to cement the stability and reliability of our findings. The F-test, executed with SPSS, yields an F-value of 4.940, decisively surpassing the threshold of 3.18. This, coupled with a p-value of 0.039, significantly below the conventional benchmark of 0.05, unequivocally confirms the superior performance and stability of our approach over CGEM.

3.2.3. Computational Cost

In the realm of computational efficiency, our method shines bright, dramatically outpacing CGEM and other evaluated algorithms. An analysis of runtimes across ten trials showcased CGEM’s considerably longer durations, with an average runtime dwarfing ours. Our method clocked in at an astonishingly low average of 5.22 ± 1.65 × 10−1 s (Figure 2 and Figure S1 and Table 4), a mere fraction of CGEM’s, thus not only illustrating our approach’s swiftness but its unparalleled efficiency and stability in the face of complex genomic data.

3.3. Comparison Study

In order to obtain a comprehensive assessment of the identification performance, an analysis of the most widely used biclustering approaches is conducted, with a comparison to our approach made using Fisher’s exact test p-values. Through meticulous analysis, juxtaposed against the rigorous benchmarks set by Fisher’s exact test p-values, we sought to illuminate the distinctive prowess of our approach. The tableau of results, as detailed in Table 5, unveils a panorama of statistical variance across the algorithms. Notably, the MCC algorithm’s p-values oscillated broadly, marking a contrast against the more consistent yet equally varied performance of the CC and LAS algorithms. The RelDenClu and CBSC algorithms yield p-values that are relatively consistent and demonstrate superior performance compared to the MCC, CC, LAS, and CGEM algorithms. Amidst these statistical results, our CGEMGA algorithm emerged as a beacon of stability and precision, boasting the lowest average p-value with unparalleled consistency.
To visually articulate these findings, Figure 3 unfolds in a duo of plots, each casting our algorithm in a comparative light against both the CGEM algorithm and the ensemble of established biclustering methods. Further enriched by plots C and D, which transform these p-values into a more visually impactful negative log10 scale, our narrative of superiority is vividly underscored. Figure 4, echoing this theme, presents a compelling graphical representation of our method’s statistical dominance, showcasing the smallest p-values amidst a backdrop of high stability and consistency, significantly outshining the comparative algorithms.
Moreover, in order to examine the convergence of the GA proposed by our approach, we undertake an analysis of the number of iterations required for convergence. The algorithm is designed to generate, on average, two modules per iteration. As illustrated in Figure 5, both modules converge to the global optimal solution value of 9.37 × 10−2 in less than five iterations, indicating a rapid convergence rate.

3.4. Functional Enrichment Analysis

While our algorithm’s prowess in navigating the complex genomic landscape of breast cancer is undeniably impressive, the true essence of this journey lies in unraveling the biological narratives of the genes thus identified. Functional enrichment analysis emerges as a vital tool in this quest, bridging the gap between gene clusters and their biological functions. By leveraging the DAVID platform, a beacon in the bioinformatics realm, we not only validate the biological significance of the identified genes but also illuminate the pathways they orchestrate. Our analysis, capturing forty-five genes within the TCGA breast cancer dataset (Table 6), reveals a rich tapestry of functional associations, underscoring the enriched biological relevance of these gene modules (Table 7).
A closer examination reveals further fascinating insights. For instance, the role of the CDH1 gene in lobular breast cancer is well documented [41,42], while the FEN1 gene is associated with poor prognoses [43]. Additionally, the complex interactions between POLQ and TP53 have been extensively studied [44,45]. These accounts, substantiated by meticulous research, not only corroborate the accuracy of our approach but also illuminate the complex interplay between genetics and cancer pathology.

4. Discussion and Conclusion

In the complex field of breast cancer genomics, we propose a state-of-the-art biclustering algorithm that draws inspiration from the adaptive power of genetic algorithms (GA). At the core of our exploration, utilizing the robust analytical frameworks of Fisher’s exact test and the F-test, we meticulously scrutinized the performance of various algorithms against a backdrop of TCGA breast cancer data. This rigorous examination not only confirmed the superior significance of our method but also highlighted its computational agility and the depth of biological insights it unveils through functional enrichment analysis. With the F-test result (F = 4.940, p = 0.039) underscoring significant distinctions from the CGEM algorithm and a swift convergence within a mere five iterations, our method’s computational efficiency shines brightly, clocking an average processing time of just 5.22 s.
Elevating the foundational CGEM algorithm, our CGEMGA incarnation introduces a strategic optimization of the initial seed s through GA, coupled with the Mean Squared Residual (MSR) serving as a continual benchmark for each iteration. This innovation ensures not only a globally optimal output but also a refined selection of candidate gene modules, adeptly navigating the vast seas of RNA-seq data to pinpoint co-expressed genomes and their co-expression conditions with unprecedented efficiency and stability. The development of CGEMGA is shown to be superior to CGEM in terms of speed and its ability to link driver genes with breast cancer.
However, our voyage is not without its navigational buoys and potential horizons for expansion. The reliance on TCGA as the sole data harbor introduces a need for broader validation across diverse genomic databases. Moreover, the exclusive use of MSR as the guiding criterion beckons the exploration of additional metrics to enrich our biclustering search compass. Historical beacons, such as the scaling MSR (SMSR) introduced by Mukhopadhyay et al. in 2009 [46], and the biclustering with iterative sorting of weighted coefficients (BISWC) approach employed by Teng and Chan in 2008 [47], which meticulously prioritize and filter features based on their significance, hint at the vast potential of integrating multiple criteria to further refine our algorithm’s accuracy and relevance.
Furthermore, in their foundational work, Gunnar Carlsson and colleagues introduced the concept of topological data analysis (TDA), which offers insights into the underlying data structures and key learning processes, thereby facilitating improvements in deep-learning performance and generalization [48]. Tianyu Zhang and colleagues have proposed a multimodal deep-learning model that fuses mammography and ultrasound images with the objective of enhancing the precision of breast cancer molecular subtype prediction [49]. This approach focuses on the most pertinent features for the prediction task through an attention mechanism. Paul Gamble has developed a deep-learning system for the direct prediction of biomarker status in breast cancer tissues [50]. Moreover, he presents a methodology for enhancing the efficacy and precision of biomarker detection by elucidating the correlation between the morphological characteristics discerned by the model and the biomarker status through interpretable analysis. Xinmin Zhang provides a comprehensive overview of the current status and future direction of molecular classification in breast cancer, emphasizing the crucial role of molecular classification in individualized therapy [51]. Moreover, he presents the development of more accurate, reliable, and straightforward molecular classification methods. Lehmann examines two variants of the rs12976445 polymorphism of the miR-125a gene in breast cancer patients, investigating their correlation with breast cancer [52]. Furthermore, he elucidates the mechanism by which the U variant may diminish the expression level of miR-125a through online simulation. The findings of this research will provide new insights and ideas for the diagnosis and treatment of breast cancer. In their review, Yasmin Cura and colleagues considered the clinical relevance of genetic polymorphisms affecting the efficacy and safety of breast cancer treatments [53]. Moreover, they emphasized the importance of pharmacogenetic guidelines based on these polymorphisms and explored the development of more precise predictive models for individuals. The aforementioned methods illustrate that breast cancer prediction research is exploring innovative avenues that integrate bioinformatics concepts, including multimodal information fusion, deep-learning models, interpretability analysis, and molecular marker discovery, with bioinformatics techniques such as data mining, machine learning, network analysis, and computational biology. This represents the central focus of our future research endeavors.
In essence, the introduced CGEMGA algorithm serves not only as a model of efficiency and robustness in the pursuit of deciphering the genomic landscape of breast cancer but also as a validation of the transformative power of evolutionary algorithms in deciphering the intricate harmonies of gene expression patterns. In the future, our research will focus on incorporating a range of fitness functions as criteria for bicluster identification. This will allow us to refine the accuracy and enhance the pragmatic utility of our computational methodology. As we pursue this course of research, the practical implications of our findings encourage further investigation, suggesting a multitude of potential applications in the ongoing fight against breast cancer.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biomedicines12092086/s1. Figure S1: Runtime of CGEM and CGEMGA in seconds.

Author Contributions

Conceptualization, Z.B. and W.Y.; methodology, J.X. (Jianing Xi); software, J.X. (Jinnan Xie); validation, Y.L. and J.C.; formal analysis, W.Y.; investigation, Z.H.; resources, Y.C.; data curation, W.Y.; writing—original draft preparation, J.X. (Jianing Xi); writing—review and editing, W.Y., J.X. (Jinnan Xie) and Z.B.; visualization, Y.L.; supervision, J.X. (Jianing Xi); project administration, J.X. (Jianing Xi); funding acquisition, J.X. (Jianing Xi) and Z.B. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Guangzhou Education Science Planning Project (Grant No. 202214340), partially by the Special Foundation in Department of Higher Education of Guangdong (Grant No. 2022ZDX2053), partially by the Discipline Construction Project of Guangzhou Medical University (Grant No. 02-445-2301244XM), partially by the Guangzhou Basic and Applied Basic Research Foundation (Grant No. 2023A04J0386), partially by the Guangdong Basic and Applied Basic Research Foundation (Grant No. 2024A1515010851), and partially by the Industry-university Cooperative Education Project (Grant No. 201902120032). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets analyzed for this study can be found in the TCGA and cBioPortal, https://www.cancer.gov/ccg/research/genome-sequencing/tcga (accessed on 15 November 2022) and http://www.cbioportal.org/ (accessed on 15 November 2022).

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Zhang, Z.; Jing, Y.; Chen, B.; Zhang, H.; Liu, T.; Dong, S.; Zhang, L.; Yan, X.; Yang, S.; Chen, L.; et al. The application of targeted RNA sequencing for the analysis of fusion genes, gene mutations, IKZF1 intragenic deletion, and CRLF2 overexpression in acute lymphoblastic leukemia. Int. J. Lab. Hematol. 2024, 46, 670–677. [Google Scholar] [CrossRef] [PubMed]
  2. Sun, P.; Luan, Y.; Cai, X.; Liu, Q.; Ren, P.; Xin, P.; Yu, Y.; Song, B.; Wang, Y.; Chang, H.; et al. Predicting mechanism of immune response in microsatellite instability colorectal cancer. Heliyon 2024, 10, e28120. [Google Scholar] [CrossRef] [PubMed]
  3. Ye, G.; Zhang, C.; Zhuang, Y.; Liu, H.; Song, E.; Li, K.; Liao, Y. An advanced nomogram model using deep learning radiomics and clinical data for predicting occult lymph node metastasis in lung adenocarcinoma. Transl. Oncol. 2024, 44, 101922. [Google Scholar] [CrossRef]
  4. Wang, L.; Hong, C.; Song, J.; Yao, J. CTEC: A cross-tabulation ensemble clustering approach for single-cell RNA sequencing data analysis. Bioinformatics 2024, 40, btae130. [Google Scholar] [CrossRef] [PubMed]
  5. Watts, J.; Allen, E.; Mitoubsi, A.; Khojandi, A.; Eales, J.; Jalali-Najafabadi, F.; Papamarkou, T. Adapting Random Forests to Predict Obesity-Associated Gene Expression. In Proceedings of the 2022 44th Annual International Conference of the IEEE Engineering in Medicine & Biology Society (EMBC), Glasgow, UK, 11–15 July 2022; pp. 4407–4410. [Google Scholar]
  6. Resmini, R.; Silva, L.; Araujo, A.S.; Medeiros, P.; Muchaluat-Saade, D.; Conci, A. Combining Genetic Algorithms and SVM for Breast Cancer Diagnosis Using Infrared Thermography. Sensors 2021, 21, 4802. [Google Scholar] [CrossRef]
  7. Seifert, S.; Gundlach, S.; Junge, O.; Szymczak, S. Integrating biological knowledge and gene expression data using pathway-guided random forests: A benchmarking study. Bioinformatics 2020, 36, 4301–4308. [Google Scholar] [CrossRef]
  8. Kim, W.-J.; Choi, B.R.; Noh, J.J.; Lee, Y.-Y.; Kim, T.-J.; Lee, J.-W.; Kim, B.-G.; Choi, C.H. Comparison of RNA-Seq and microarray in the prediction of protein expression and survival prediction. Front. Genet. 2024, 15, 1342021. [Google Scholar] [CrossRef]
  9. Huang, Y.; Zhao, Y.; Capstick, A.; Palermo, F.; Haddadi, H.; Barnaghi, P. Analyzing entropy features in time-series data for pattern recognition in neurological conditions. Artif. Intell. Med. 2024, 150, 102821. [Google Scholar] [CrossRef]
  10. Ha, C.T.; Tageldein, M.M.; Harding, S.M. The entanglement of DNA damage and pattern recognition receptor signaling. DNA Repair 2024, 133, 103595. [Google Scholar] [CrossRef]
  11. Hauschild, A.-C.; Lemanczyk, M.; Matschinske, J.; Frisch, T.; Zolotareva, O.; Holzinger, A.; Baumbach, J.; Heider, D. Federated Random Forests can improve local performance of predictive models for various healthcare applications. Bioinformatics 2022, 38, 2278–2286. [Google Scholar] [CrossRef]
  12. Chu, H.-M.; Liu, J.-X.; Zhang, K.; Zheng, C.-H.; Wang, J.; Kong, X.-Z. A binary biclustering algorithm based on the adjacency difference matrix for gene expression data analysis. BMC Bioinform. 2022, 23, 381. [Google Scholar] [CrossRef]
  13. Cheng, Y.; Church, G.M. Biclustering of expression data. In Proceedings of the 8th International Conference on Intelligent Systems for Molecular Biology (ISMB 2000), San Diego, CA, USA, 19–23 August 2000; Volume 8, pp. 93–103. [Google Scholar]
  14. Hanczar, B.; Nadif, M. Ensemble methods for biclustering tasks. Pattern Recognit. 2012, 45, 3938–3949. [Google Scholar] [CrossRef]
  15. Andrey, A.S.; Weigman, V.J.; Perou, C.M.; Nobel, A.B. Finding large average submatrices in high dimensional data. Ann. Appl. Stat. 2009, 3, 985–1012. [Google Scholar]
  16. Jain, N.; Ghosh, S.; Murthy, C.A. RelDenClu: A Relative Density based Biclustering Method for identifying non-linear feature relations. arXiv 2021, arXiv:1811.04661. [Google Scholar]
  17. Jain, N.; Murthy, C.A. Connectedness-based subspace clustering. Knowl. Inf. Syst. 2019, 58, 9–34. [Google Scholar] [CrossRef]
  18. Xi, J.; Wang, M.; Li, A. DGPathinter: A novel model for identifying driver genes via knowledge-driven matrix factorization with prior knowledge from interactome and pathways. PeerJ Comput. Sci. 2017, 3, e133. [Google Scholar] [CrossRef]
  19. Xi, J.; Li, A.; Wang, M. HetRCNA: A Novel Method to Identify Recurrent Copy Number Alternations from Heterogeneous Tumor Samples Based on Matrix Decomposition Framework. IEEE/ACM Trans. Comput. Biol. Bioinform. 2020, 17, 422–434. [Google Scholar] [CrossRef]
  20. Murali, T.M.; Kasif, S. Extracting conserved gene expression motifs from gene expression data. In Proceedings of the Pacific Symposium on Biocomputing (PSB), Kauai, HI, USA, 3–7 January 2003; pp. 77–88. [Google Scholar]
  21. Williams, A.; Halappanavar, S. Application of bi-clustering of gene expression data and gene set enrichment analysis methods to identify potentially disease causing nanomaterials. Data Brief 2017, 15, 933–940. [Google Scholar] [CrossRef]
  22. Madeira, S.C.; Oliveira, A.L. Biclustering algorithms for biological data analysis: A survey. IEEE/ACM Trans. Comput. Biol. Bioinform. 2004, 1, 24–45. [Google Scholar] [CrossRef]
  23. Tang, C.Y.; Hung, C.L.; Zheng, H.; Lin, C.Y.; Jiang, H. Novel Computational Technologies for Next-Generation Sequencing Data Analysis and Their Applications. Int. J. Genom. 2015, 2015, 254685. [Google Scholar] [CrossRef]
  24. Xi, J.; Yuan, X.; Wang, M.; Li, A.; Li, X.; Huang, Q. Inferring subgroup-specific driver genes from heterogeneous cancer samples via subspace learning with subgroup indication. Bioinformatics 2020, 36, 1855–1863. [Google Scholar] [CrossRef] [PubMed]
  25. Craven, K.E.; Gökmen-Polar, Y.; Badve, S.S. CIBERSORT analysis of TCGA and METABRIC identifies subgroups with better outcomes in triple negative breast cancer. Sci. Rep. 2021, 11, 4691. [Google Scholar] [CrossRef]
  26. Thennavan, A.; Beca, F.; Xia, Y.; Garcia-Recio, S.; Allison, K.; Collins, L.C.; Tse, G.M.; Chen, Y.-Y.; Schnitt, S.J.; Hoadley, K.A.; et al. Molecular analysis of TCGA breast cancer histologic types. Cell Genom. 2021, 1, 100067. [Google Scholar] [CrossRef] [PubMed]
  27. Linehan, W.M.; Ricketts, C.J. The Cancer Genome Atlas of renal cell carcinoma: Findings and clinical implications. Nat. Rev. Urol. 2019, 16, 539–552. [Google Scholar] [CrossRef] [PubMed]
  28. Rau, A.; Flister, M.; Rui, H.; Auer, P.L. Exploring drivers of gene expression in the Cancer Genome Atlas. Bioinformatics 2019, 35, 62–68. [Google Scholar] [CrossRef]
  29. Malhotra, S.; Alsulami, A.F.; Heiyun, Y.; Ochoa, B.M.; Jubb, H.; Forbes, S.; Blundell, T.L. Understanding the impacts of missense mutations on structures and functions of human cancer-related genes: A preliminary computational analysis of the COSMIC Cancer Gene Census. PLoS ONE 2019, 14, e0219935. [Google Scholar] [CrossRef]
  30. Alsulami, A.F.; Torres, P.H.; Moghul, I.; Arif, S.M.; Chaplin, A.; Vedithi, S.; Blundell, T. COSMIC Cancer Gene Census 3D database: Understanding the impacts of mutations on cancer targets. Brief. Bioinform. 2021, 22, bbab220. [Google Scholar] [CrossRef]
  31. De Jong, K. Learning with genetic algorithms: An overview. Mach. Learn. 1988, 3, 121–138. [Google Scholar] [CrossRef]
  32. Liu, X.; Yu, T.; Zhao, X.; Long, C.; Han, R.; Su, Z.; Li, G. ARBic: An all-round biclustering algorithm for analyzing gene expression data. NAR Genom. Bioinform. 2023, 5, lqad009. [Google Scholar] [CrossRef]
  33. Hastie, T.; Tibshirani, R.; Eisen, M.B.; Alizadeh, A.; Levy, R.; Staudt, L.; Chan, W.C.; Botstein, D.; Brown, P. ‘Gene shaving’ as a method for identifying distinct sets of genes with similar expression patterns. Genome Biol. 2000, 1, research0003.1. [Google Scholar] [CrossRef]
  34. Li, J.; Wong, L. Identifying good diagnostic gene groups from gene expression profiles using the concept of emerging patterns. Bioinformatics 2002, 18, 725–734. [Google Scholar] [CrossRef]
  35. Castanho, E.N.; Aidos, H.; Madeira, S.C. Biclustering fMRI time series: A comparative study. BMC Bioinform. 2022, 23, 192. [Google Scholar] [CrossRef] [PubMed]
  36. Colin, P.J.; Eleveld, D.J.; Thomson, A.H. Genetic Algorithms as a Tool for Dosing Guideline Optimization: Application to Intermittent Infusion Dosing for Vancomycin in Adults. CPT Pharmacomet. Syst. Pharmacol. 2020, 9, 294–302. [Google Scholar] [CrossRef] [PubMed]
  37. Connelly, L.M. Fisher’s Exact Test. Medsurg. Nurs. 2016, 25, 58–61. [Google Scholar]
  38. Lopez-Gimenez, M.R.; Garcia Gomez, J.J. The Fisher’s test. Med. Clin. 1993, 101, 156–157. [Google Scholar]
  39. Blevins, L.; McDonald, C.J. Fisher’s Exact Test: An easy-to-use statistical test for comparing outcomes. MD Comput. 1985, 2, 15–19, 68. [Google Scholar]
  40. Wang, S.; Cui, H. Generalized F test for high dimensional linear regression coefficients. J. Multivar. Anal. 2013, 117, 134–149. [Google Scholar] [CrossRef]
  41. Sammons, S.; Elliott, A.; Barroso-Sousa, R.; Chumsri, S.; Tan, A.R.; Sledge, G.W., Jr.; Tolaney, S.M.; Torres, E.T.R. Concurrent predictors of an immune responsive tumor microenvironment within tumor mutational burden-high breast cancer. Front. Oncol. 2023, 13, 1235902. [Google Scholar] [CrossRef] [PubMed]
  42. Perez-Duran, J.; Luna, A.; Portilla, A.; Martínez, P.; Ceballos, G.; Ortíz-Flores, M.Á.; Solis-Paredes, J.M.; Nájera, N. (-)-Epicatechin Inhibits Metastatic-Associated Proliferation, Migration, and Invasion of Murine Breast Cancer Cells In Vitro. Molecules 2023, 28, 6229. [Google Scholar] [CrossRef]
  43. Xu, L.; Shen, J.M.; Qu, J.L.; Song, N.; Che, X.F.; Hou, K.Z.; Shi, J.; Zhao, L.; Shi, S.; Liu, Y.P.; et al. FEN1 is a prognostic biomarker for ER+ breast cancer and associated with tamoxifen resistance through the ERα/cyclin D1/Rb axis. Ann. Transl. Med. 2021, 9, 258. [Google Scholar] [CrossRef]
  44. Kim, J.; Jeong, K.; Jun, H.; Kim, K.; Bae, J.M.; Song, M.G.; Yi, H.; Park, S.; Woo, G.U.; Lee, D.W.; et al. Mutations of TP53 and genes related to homologous recombination repair in breast cancer with germline BRCA1/2 mutations. Hum. Genom. 2023, 17, 2. [Google Scholar] [CrossRef]
  45. Grote, I.; Poppe, A.; Lehmann, U.; Christgen, M.; Kreipe, H.; Bartels, S. Frequency of genetic alterations differs in advanced breast cancer between metastatic sites. Genes Chromosomes Cancer 2023, 63, e23199. [Google Scholar] [CrossRef]
  46. Mukhopadhyay, A.; Maulik, U.; Bandyopadhyay, S. A novel coherence measure for discovering scaling biclusters from gene expression data. J. Bioinform. Comput. Biol. 2009, 7, 853–868. [Google Scholar] [CrossRef]
  47. Teng, L.; Chan, L. Discovering Biclusters by Iteratively Sorting with Weighted Correlation Coefficient in Gene Expression Data. J. Signal Process. Syst. 2008, 50, 267–280. [Google Scholar] [CrossRef]
  48. Carlsson, G.E.; Gabrielsson, R.B. Topological Approaches to Deep Learning. arXiv 2018, arXiv:1811.01122. [Google Scholar]
  49. Zhang, T.; Tan, T.; Han, L.; Appelman, L.; Veltman, J.; Wessels, R.; Duvivier, K.M.; Loo, C.; Gao, Y.; Wang, X.; et al. Predicting breast cancer types on and beyond molecular level in a multi-modal fashion. npj Breast Cancer 2023, 9, 16. [Google Scholar] [CrossRef]
  50. Gamble, P.; Jaroensri, R.; Wang, H.; Tan, F.; Moran, M.; Brown, T.; Flament-Auvigne, I.; Rakha, E.A.; Toss, M.; Dabbs, D.J.; et al. Determining breast cancer biomarker status and associated morphological features using deep learning. Commun. Med. 2021, 1, 14. [Google Scholar] [CrossRef] [PubMed]
  51. Zhang, X. Molecular Classification of Breast Cancer: Relevance and Challenges. Arch. Pathol. Lab. Med. 2022, 147, 46–51. [Google Scholar] [CrossRef] [PubMed]
  52. Lehmann, T.P.; Miskiewicz, J.; Szostak, N.; Szachniuk, M.; Grodecka-Gazdecka, S.; Jagodziński, P.P. In Vitro and in Silico Analysis of miR-125a with rs12976445 Polymorphism in Breast Cancer Patients. Appl. Sci. 2020, 10, 7275. [Google Scholar] [CrossRef]
  53. Cura, Y.; Ramírez, C.P.; Martín, A.S.; Martínez, F.M.; Hernández, M.C.; Tortosa, M.d.C.R.; Morales, A.J. Genetic polymorphisms on the effectiveness or safety of breast cancer treatment: Clinical relevance and future perspectives. Mutat. Res./Rev. Mutat. Res. 2021, 788, 108391. [Google Scholar] [CrossRef]
Figure 1. Flowchart of the CGEMGA algorithm.
Figure 1. Flowchart of the CGEMGA algorithm.
Biomedicines 12 02086 g001
Figure 2. Runtime of CGEM and CGEMGA in seconds (the solid red line in the figure represents the median value, while the blue dashed line represents the mean value).
Figure 2. Runtime of CGEM and CGEMGA in seconds (the solid red line in the figure represents the median value, while the blue dashed line represents the mean value).
Biomedicines 12 02086 g002
Figure 3. The comparison of p-values between different methods. Where (A) stands for p values comparison between our approach and CGEM, and (B) stands for p values comparison between our approach and RelDenClu, CBSC, LAS, MCC, and CC. (C,D) are negative log10 operations on the p values in the (A,B).
Figure 3. The comparison of p-values between different methods. Where (A) stands for p values comparison between our approach and CGEM, and (B) stands for p values comparison between our approach and RelDenClu, CBSC, LAS, MCC, and CC. (C,D) are negative log10 operations on the p values in the (A,B).
Biomedicines 12 02086 g003
Figure 4. Comparison of mean p-value of five methods.
Figure 4. Comparison of mean p-value of five methods.
Biomedicines 12 02086 g004
Figure 5. Iterative convergence curve of the genetic algorithm with MSR as a fitness function in our approach.
Figure 5. Iterative convergence curve of the genetic algorithm with MSR as a fitness function in our approach.
Biomedicines 12 02086 g005
Table 1. A comparative analysis of the characteristics of the six algorithms, namely CC, MCC, LAS, CGEM, RelDenClu, and CBSC.
Table 1. A comparative analysis of the characteristics of the six algorithms, namely CC, MCC, LAS, CGEM, RelDenClu, and CBSC.
CharacteristicsAlgorithms
CCMCCLASCGEMRelDenCluCBSC
Core IdeaIterative spectral method for finding co-expressed gene and condition submatricesEnsemble learning method combining multiple base biclustering algorithms for improved stabilityStatistical method for finding large average submatrices in high-dimensional data, focusing on numerical featuresGraph-based method for extracting conserved gene expression motifsFind sets of observations with high local densityFind subspaces with high connectivity
Algorithm TypeSpectral clusteringEnsemble learningStatistical methodGraph-based methodDensity clusteringConnectivity clustering
Time ComplexityHighModerate to high (depending on the number and type of base algorithms)HighModerateHighHigh
Space ComplexityModerateHigh (requires storage of results from multiple base algorithms)High (requires storage of extensive submatrix information)ModerateHighHigh
Applicable Data TypeExpression dataGeneral (applicable to various types of data)General (applicable to various types of data)Expression dataMicroarray dataHigh-dimensional data
RobustnessModerate (sensitive to noise)High (ensemble methods reduce the impact of noise)High (statistical methods have some resistance to noise)Moderate (depends on the stability of the graph structure)HighHigh
ScalabilityModerate (suitable for medium to small-scale data)Good (can be scaled to large-scale data)Good (suitable for large-scale data)Moderate (suitable for medium to small-scale data)Moderate (suitable for medium to small-scale data)Medium
Pattern Type DiscoveredCo-expression patternsDiverse patterns (depending on the base algorithms)Numerically significant submatricesConserved expression motifsNonlinear relationships between featuresSubspace structures based on connectivity
Real-world ApplicationsGene expression analysis in bioinformaticsWidely applied in bioinformatics and machine learningBioinformatics, image processing, and other fieldsGene network analysis in bioinformaticsGene functional grouping in microarray dataGene functional grouping in bioinformatics
Table 2. Pseudo-code of CGEMGA algorithm.
Table 2. Pseudo-code of CGEMGA algorithm.
Algorithm 1 FINDMODULE(): algorithm for computing the largest module.
1. for i = 1 to ns do
2. GA begin
3. Create an initial population of n chromosomes Ci (i = 1, 2, …, n) as seeds
4. Set iteration counter t = 0
5. Choose a subset Di of the samples with size sd
6. For every gene g in Di, include the pair (g, s) in the set Gij if g is in the state s in c and all Di samples
7. Cij = set of samples that agree with c in all the gene-states in Gij
8. Calculate the MSR fitness value for each chromosome
9. while (t < MAX)
10. Select a pair of chromosomes form initial population based on MSR fitness
11. Apply crossover operation on selected pair with crossover probability
12. Apply mutation on the offspring with mutation probability
13. Replace old population with newly generated population
14. Increment the current iteration t by 1.
15. end while
16. Discard (Cij, Gij) if Cij contains less than αn samples.
17. returen the best solution, Ci with min MSR
18. GA end
19. return the module (C*, G*) that maximises |Gij|, 1 ≤ ins
Table 3. Comparison of Fisher’s test p-values and MSR values of CGEM with our approach on 10 independent runs.
Table 3. Comparison of Fisher’s test p-values and MSR values of CGEM with our approach on 10 independent runs.
No.Method
CGEMCGEMGA
p ValueMSR Valuep ValueMSR Value
19.62 × 10−64.71 × 10−21.13 × 10−49.37 × 10−2
22.47 × 10−51.27 × 10−11.19 × 10−49.37 × 10−2
31.73 × 10−41.29 × 1001.25 × 10−49.37 × 10−2
41.84 × 10−41.44 × 1001.28 × 10−49.37 × 10−2
56.54 × 10−45.44 × 1001.55 × 10−49.37 × 10−2
61.20 × 10−36.89 × 1001.63 × 10−49.37 × 10−2
71.90 × 10−31.15 × 101.75 × 10−49.37 × 10−2
86.00 × 10−32.95 × 101.81 × 10−49.37 × 10−2
91.08 × 10−25.78 × 101.91 × 10−49.37 × 10−2
109.71 × 10−27.12 × 1021.91 × 10−49.37 × 10−2
Mean ± SD1.18 × 10−2 ± 3.01 × 10−28.26 × 10 ± 2.21 × 1021.54 × 10−4 ± 3.06 × 10−59.37 × 10−2 ± 0
Table 4. Runtime of CGEM and our approach in seconds.
Table 4. Runtime of CGEM and our approach in seconds.
MethodTimes (n = 10)Mean ± SD
12345678910
CGEM1.89 × 1031.79 × 1031.77 × 1031.23 × 1031.21 × 1031.16 × 1031.13 × 1031.01 × 1039.59 × 1029.32 × 1021.31 × 103 ± 3.66 × 102
CGEMGA4.94 × 1005.02 × 1005.11 × 1005.13 × 1005.25 × 1005.29 × 1005.30 × 1005.33 × 1005.39 × 1005.45 × 1005.22 × 100 ± 1.65 × 10−1
Table 5. Fisher’s test p-values for 10 independent runs of the seven methods.
Table 5. Fisher’s test p-values for 10 independent runs of the seven methods.
MethodFisher’ Test p Values (n = 10)Mean ± SD
12345678910
MCC1.50 × 10−37.30 × 10−39.20 × 10−31.03 × 10−11.75 × 10−14.31 × 10−15.09 × 10−16.67 × 10−18.98 × 10−19.03 × 10−13.70 × 10−1 ± 3.62 × 10−1
CC1.50 × 10−37.30 × 10−34.99 × 10−21.74 × 10−11.78 × 10−11.95 × 10−14.12 × 10−15.59 × 10−16.73 × 10−18.54 × 10−13.10 × 10−1 ± 3.00 × 10−1
LAS6.68 × 10−56.44 × 10−46.44 × 10−41.70 × 10−36.70 × 10−33.04 × 10−26.65 × 10−27.56 × 10−25.44 × 10−16.26 × 10−11.35 × 10−1 ± 2.40 × 10−1
CGEM9.62 × 10−62.47 × 10−51.73 × 10−41.84 × 10−46.54 × 10−41.20 × 10−31.90 × 10−36.00 × 10−31.08 × 10−29.71 × 10−21.18 × 10−2 ± 3.01 × 10−2
RelDenClu1.40 × 10−681.43 × 10−671.34 × 10−341.35 × 10−346.33 × 10−343.05 × 10−137.49 × 10−131.44 × 10−81.45 × 10−22.01 × 10−23.46 × 10−3 ± 6.98 × 10−3
CBSC1.54 × 10−131.57 × 10−131.48 × 10−123.35 × 10−111.48 × 10−106.97 × 10−108.23 × 10−61.58 × 10−41.60 × 10−42.21 × 10−22.24 × 10−3 ± 7.41 × 10−3
CGEMGA1.13 × 10−41.19 × 10−41.25 × 10−41.28 × 10−41.55 × 10−41.63 × 10−41.75 × 10−41.81 × 10−41.91 × 10−41.91 × 10−41.54 × 10−4 ± 3.06 × 10−5
Table 6. The gene find by our approach form the TCGA breast cancer data.
Table 6. The gene find by our approach form the TCGA breast cancer data.
No.Gene
Symbol
NameCytogenetic BandNo.Gene
Symbol
NameCytogenetic Band
1MTORmechanistic target of rapamycin1p36.2224SDHAF2succinate dehydrogenase complex assembly factor 211q12.2
2SF3B1splicing factor 3b, subunit 1, 155 kDa2q33.125KDM5Alysine (K)-specific demethylase 5A, JARID1A12p13.33
3POLQDNA polymerase theta3q13.3326PRPF40Bpre-mRNA processing factor 40 homolog B12q13.12
4MECOMMDS1 and EVI1 complex locus3q26.227NCOR2nuclear receptor corepressor 212q24.31
5TET2tet oncogene family member 24q2428RAD51BRAD51 paralog B14q24.1
6FAT1FAT atypical cadherin 14q35.229TCL1AT-cell leukemia/lymphoma 1A14q32.13
7TLX3T-cell leukemia, homeobox 3 (HOX11L2)5q35.130DROSHAdrosha ribonuclease III15p13.3
8SRSF3serine/arginine-rich splicing factor 36p21.3131CHD2chromodomain helicase DNA binding protein 215q26.1
9DEKDEK oncogene (DNA binding)6p22.332PRKCBprotein kinase C beta16p12.2
10SGK1serum/glucocorticoid regulated kinase 16q23.233RMI2RecQ mediated genome instability 216p13.13
11EZRezrin6q25.334CDH1cadherin 1, type 1, E-cadherin (epithelial) (ECAD)16q22.1
12MACC1MET transcriptional regulator MACC17p21.135TP53tumor protein p5317p13.1
13SBDSShwachman-Bodian-Diamond syndrome protein7q11.2136KAT7lysine acetyltransferase 717q21.33
14CUX1cut-like homeobox 17q22.137SRSF2serine/arginine-rich splicing factor 217q25.2
15KAT6AK(lysine) acetyltransferase 6A8p11.2138KDSR3-ketodihydrosphingosine reductase18q21.33
16GNAQguanine nucleotide binding protein (Gprotein), q polypeptide9q21.239CEP89centrosomal protein 89 kDa19q13.11
17CNTRLcentriolin9q33.240ARHGAP35Rho GTPase activating protein 3519q13.32
18LARP4BLa ribonucleoprotein domain family member 4B10p15.341TOP1topoisomerase (DNA) I20q12
19A1CFAPOBEC1 complementation factor10q11.2342KDM5Clysine (K)-specific demethylase 5C (JARID1C)Xp11.22
20KAT6BK(lysine) acetyltransferase 6B10q22.243KDM6Alysine (K)-specific demethylase 6A, UTXXp11.3
21NUP98nucleoporin 98kDa11p15.444TMSB4XThymosin Beta 4 X-LinkedXp22.2
22CLP1cleavage and polyadenylation factor I subunit 111q12.145CRLF2cytokine receptor-like factor 2Xp22.33
23FEN1flap structure-specific endonuclease 111q12.2
Table 7. The functional enrichment analysis results of the discovered drivers of our approach.
Table 7. The functional enrichment analysis results of the discovered drivers of our approach.
TermPercentagep-ValueFDR
hsa05205:Proteoglycans in cancer11.13.8 × 10−35.98 × 10−1
hsa05214:Glioma6.72.4 × 10−27.75 × 10−1
hsa04971:Gastric acid secretion6.72.4 × 10−27.75 × 10−1
hsa05200:Pathways in cancer13.32.4 × 10−27.75 × 10−1
h_pkcPathway:Activation of PKC through G protein coupled receptor4.43.0 × 10−27.53 × 10−1
hsa03040:Spliceosome8.93.1 × 10−27.75 × 10−1
hsa05163:Human cytomegalovirus infection8.93.4 × 10−27.75 × 10−1
hsa04670:Leukocyte transendothelial migration6.75.1 × 10−27.75 × 10−1
hsa04935:Growth hormone synthesis, secretion and action6.75.6 × 10−27.75 × 10−1
hsa04071:Sphingolipid signaling pathway6.75.6 × 10−27.75 × 10−1
hsa04919:Thyroid hormone signaling pathway6.75.6 × 10−27.75 × 10−1
h_myosinPathway:PKC-catalyzed phosphorylation of inhibitory phosphoprotein of myosin phosphatase4.45.9 × 10−27.53 × 10−1
h_ccr5Pathway:Pertussis toxin-insensitive CCR5 Signaling in Macrophage4.47.1 × 10−27.53 × 10−1
hsa04371:Apelin signaling pathway6.77.2 × 10−27.75 × 10−1
hsa05206:MicroRNAs in cancer8.97.4 × 10−27.75 × 10−1
hsa05017:Spinocerebellar ataxia6.77.6 × 10−27.75 × 10−1
h_calcineurinPathway:Effects of calcineurin in Keratinocyte Differentiation4.47.9 × 10−27.53 × 10−1
hsa05226:Gastric cancer6.78.1 × 10−27.75 × 10−1
hsa04150:mTOR signaling pathway6.78.8 × 10−27.75 × 10−1
h_par1pathway:Thrombin signaling and protease-activated receptors4.49.1 × 10−27.53 × 10−1
h_chemicalPathway:Apoptotic Signaling in Response to DNA Damage4.49.1 × 10−27.53 × 10−1
h_ccr3Pathway:CCR3 signaling in Eosinophils4.49.5 × 10−27.53 × 10−1
h_eif4Pathway:Regulation of eIF4e and p70 S6 Kinase4.49.9 × 10−27.53 × 10−1
h_cxcr4Pathway:CXCR4 Signaling Pathway4.49.9 × 10−27.53 × 10−1
hsa05225:Hepatocellular carcinoma6.71.0 × 10−17.75 × 10−1
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yuan, W.; Li, Y.; Han, Z.; Chen, Y.; Xie, J.; Chen, J.; Bi, Z.; Xi, J. Evolutionary Mechanism Based Conserved Gene Expression Biclustering Module Analysis for Breast Cancer Genomics. Biomedicines 2024, 12, 2086. https://doi.org/10.3390/biomedicines12092086

AMA Style

Yuan W, Li Y, Han Z, Chen Y, Xie J, Chen J, Bi Z, Xi J. Evolutionary Mechanism Based Conserved Gene Expression Biclustering Module Analysis for Breast Cancer Genomics. Biomedicines. 2024; 12(9):2086. https://doi.org/10.3390/biomedicines12092086

Chicago/Turabian Style

Yuan, Wei, Yaming Li, Zhengpan Han, Yu Chen, Jinnan Xie, Jianguo Chen, Zhisheng Bi, and Jianing Xi. 2024. "Evolutionary Mechanism Based Conserved Gene Expression Biclustering Module Analysis for Breast Cancer Genomics" Biomedicines 12, no. 9: 2086. https://doi.org/10.3390/biomedicines12092086

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop