A comparison of the TempO-Seq and Affymetrix microarray platform using RTqPCR validation.

Next-generation risk assessment relies on mechanistic data from new approach methods, including transcriptome data. Various technologies, such as high-throughput targeted sequencing methods and microarray technologies based on hybridization with complementary probes, are used to determine differentially expressed genes (DEGs). The integration of data from different technologies requires a good understanding of the differences arising from the use of various technologies.To better understand the differences between the TempO-Seq platform and Affymetrix chip technology, whole-genome data for the volatile compound dimethylamine were compared. Selected DEGs were also confirmed using RTqPCR validation. Although the overlap of DEGs between TempO-Seq and Affymetrix was no higher than 37%, a comparison of the gene regulation in terms of log2fold changes revealed a very high concordance. RTqPCR confirmed the majority of DEGs from either platform in the examined dataset. Only a few conflicts were found (11%), while 22% were not confirmed, and 3% were not detected.Despite the observed differences between the two platforms, both can be validated using RTqPCR. Here we highlight some of the differences between the two platforms and discuss their applications in toxicology.

Improving prediction of blood cancer using leukemia microarray gene data and Chi2 features with weighted convolutional neural network.

Blood cancer has emerged as a growing concern over the past decade, necessitating early diagnosis for timely and effective treatment. The present diagnostic method, which involves a battery of tests and medical experts, is costly and time-consuming. For this reason, it is crucial to establish an automated diagnostic system for accurate predictions. A particular field of focus in medical research is the use of machine learning and leukemia microarray gene data for blood cancer diagnosis. Even with a great deal of research, more improvements are needed to reach the appropriate levels of accuracy and efficacy. This work presents a supervised machine-learning algorithm for blood cancer prediction. This work makes use of the 22,283-gene leukemia microarray gene data. Chi-squared (Chi2) feature selection methods and the synthetic minority oversampling technique (SMOTE)-Tomek resampling is used to overcome issues with imbalanced and high-dimensional datasets. To balance the dataset for each target class, SMOTE-Tomek creates synthetic data, and Chi2 chooses the most important features to train the learning models from 22,283 genes. A novel weighted convolutional neural network (CNN) model is proposed for classification, utilizing the support of three separate CNN models. To determine the importance of the proposed approach, extensive experiments are carried out on the datasets, including a performance comparison with the most advanced techniques. Weighted CNN demonstrates superior performance over other models when coupled with SMOTE-Tomek and Chi2 techniques, achieving a remarkable 99.9% accuracy. Results from k-fold cross-validation further affirm the supremacy of the proposed model.

mLiftOver: harmonizing data across Infinium DNA methylation platforms.

MOTIVATION: Infinium DNA methylation BeadChips are widely used for genome-wide DNA methylation profiling at the population scale. Recent updates to probe content and naming conventions in the EPIC version 2 (EPICv2) arrays have complicated integrating new data with previous Infinium array platforms, such as the MethylationEPIC (EPIC) and the HumanMethylation450 (HM450) BeadChip. RESULTS: We present mLiftOver, a user-friendly tool that harmonizes probe ID, methylation level, and signal intensity data across different Infinium platforms. It manages probe replicates, missing data imputation, and platform-specific bias for accurate data conversion. We validated the tool by applying HM450-based cancer classifiers to EPICv2 cancer data, achieving high accuracy. Additionally, we successfully integrated EPICv2 healthy tissue data with legacy HM450 data for tissue identity analysis and produced consistent copy number profiles in cancer cells. AVAILABILITY AND IMPLEMENTATION: mLiftOver is implemented R and available in the Bioconductor package SeSAMe (version 1.21.13+): https://bioconductor.org/packages/release/bioc/html/sesame.html. Analysis of EPIC and EPICv2 platform-specific bias and high-confidence mapping is available at https://github.com/zhou-lab/InfiniumAnnotationV1/raw/main/Anno/EPICv2/EPICv2ToEPIC_conversion.tsv.gz. The source code is available at https://github.com/zwdzwd/sesame/blob/devel/R/mLiftOver.R under the MIT license.

Cancer Sample Analysis Utilizing Single-Nucleotide Polymorphism Array and Array Comparative Genomic Hybridization.

Chromosomal microarray, including single-nucleotide polymorphism (SNP) array and array comparative genomic hybridization (aCGH), enables the detection of DNA copy-number loss and/or gain associated with unbalanced chromosomal aberrations. In addition, SNP array and aCGH with SNP component also detect copy-neutral loss of heterozygosity (CN-LOH). Here we describe the chromosomal microarray procedure from the sample preparation using extracted DNA to the scanning of the array chip.

Differential methylation region detection via an array-adaptive normalized kernel-weighted model.

A differentially methylated region (DMR) is a genomic region that has significantly different methylation patterns between biological conditions. Identifying DMRs between different biological conditions is critical for developing disease biomarkers. Although methods for detecting DMRs in microarray data have been introduced, developing methods with high precision, recall, and accuracy in determining the true length of DMRs remains a challenge. In this study, we propose a normalized kernel-weighted model to account for similar methylation profiles using the relative probe distance from "nearby" CpG sites. We also extend this model by proposing an array-adaptive version in attempt to account for the differences in probe spacing between Illumina's Infinium 450K and EPIC bead array respectively. We also study the asymptotic results of our proposed statistic. We compare our approach with a popular DMR detection method via simulation studies under large and small treatment effect settings. We also discuss the susceptibility of our method in detecting the true length of the DMRs under these two settings. Lastly, we demonstrate the biological usefulness of our method when combined with pathway analysis methods on oral cancer data. We have created an R package called idDMR, downloadable from GitHub repository with link: https://github.com/DanielAlhassan/idDMR, that allows for the convenient implementation of our array-adaptive DMR method.

Comparing preprocessing strategies for 3D-Gene microarray data of extracellular vesicle-derived miRNAs.

BACKGROUND: Extracellular vesicle-derived (EV)-miRNAs have potential to serve as biomarkers for the diagnosis of various diseases. miRNA microarrays are widely used to quantify circulating EV-miRNA levels, and the preprocessing of miRNA microarray data is critical for analytical accuracy and reliability. Thus, although microarray data have been used in various studies, the effects of preprocessing have not been studied for Toray's 3D-Gene chip, a widely used measurement method. We aimed to evaluate batch effect, missing value imputation accuracy, and the influence of preprocessing on measured values in 18 different preprocessing pipelines for EV-miRNA microarray data from two cohorts with amyotrophic lateral sclerosis using 3D-Gene technology. RESULTS: Eighteen different pipelines with different types and orders of missing value completion and normalization were used to preprocess the 3D-Gene microarray EV-miRNA data. Notable results were suppressed in the batch effects in all pipelines using the batch effect correction method ComBat. Furthermore, pipelines utilizing missForest for missing value imputation showed high agreement with measured values. In contrast, imputation using constant values for missing data exhibited low agreement. CONCLUSIONS: This study highlights the importance of selecting the appropriate preprocessing strategy for EV-miRNA microarray data when using 3D-Gene technology. These findings emphasize the importance of validating preprocessing approaches, particularly in the context of batch effect correction and missing value imputation, for reliably analyzing data in biomarker discovery and disease research.

Development and verification of a 10K liquid chip for Hainan black goat based on genotyping by pinpoint sequencing of liquid captured targets.

BACKGROUND: China has thousands years of goat breeding and abundant goat genetic resources. Additionally, the Hainan black goat is one of the high-quality local goat breeds in China. In order to conserve the germplasm resources of the Hainan black goat, facilitate its genetic improvement and further protect the genetic diversity of goats, it is urgent to develop a single nucleotide polymorphism (SNP) chip for Hainan black goat. RESULTS: In this study, we aimed to design a 10K liquid chip for Hainan black goat based on genotyping by pinpoint sequencing of liquid captured targets (cGPS). A total of 45,588 candidate SNP sites were obtained, 10,677 of which representative SNP sites were selected to design probes, which finally covered 9,993 intervals and formed a 10K cGPS liquid chip for Hainan black goat. To verify the 10K cGPS liquid chip, some southern Chinese goat breeds and a sheep breed with similar phenotype to the Hainan black goat were selected. A total of 104 samples were used to verify the clustering ability of the 10K cGPS liquid chip for Hainan black goat. The results showed that the detection rate of sites was 97.34% -99.93%. 84.5% of SNP sites were polymorphic. The heterozygosity rate was 3.08%-36.80%. The depth of more than 99.4% sites was above 10X. The repetition rate was 99.66%-99.82%. The average consistency between cGPS liquid chip results and resequencing results was 85.58%. In addition, the phylogenetic tree clustering analysis verified that the SNP sites on the chip had better clustering ability. CONCLUSION: These results indicate that we have successfully realized the development and verification of the 10K cGPS liquid chip for Hainan black goat, which provides a useful tool for the genome analysis of Hainan black goat. Moreover, the 10K cGPS liquid chip is conducive to the research and protection of Hainan black goat germplasm resources and lays a solid foundation for its subsequent breeding work.

Development and validation of the Axiom-MaruPri SNP chip for genetic analyses of domesticated old world camelids.

Camels play a crucial socio-economic role in sustaining the livelihoods of millions in arid and semi-arid regions. They possess remarkable physiological attributes which enable them to thrive in extreme environments, and provide a source of meat, milk and transportation. With their unique traits, camels embody an irreplaceable source of untapped genomic knowledge. This study introduces Axiom-MaruPri, a medium-density SNP chip meticulously designed and validated for both Camelus bactrianus and Camelus dromedarius. Comprising of 182,122 SNP markers, derived from the re-sequenced data of nine Indian dromedary breeds and the double-humped Bactrian camel, this SNP chip offers 34,894 markers that display polymorphism in both species. It achieves an estimated inter-marker distance of 14 Kb, significantly enhancing the coverage of the camel genome. The medium-density chip has been successfully genotyped using 480 camel samples, achieving an impressive 99 % call rate, with 96 % of the 182,122 SNPs being highly reliable for genotyping. Phylogenetic analysis and Discriminant Analysis of Principal Components yield clear distinctions between Bactrian camels and dromedaries. Moreover, the discriminant functions substantially enhance the classification of dromedary camels into different breeds. The clustering of various camel breeds reveals an apparent correlation between geographical and genetic distances. The results affirm the efficacy of this SNP array, demonstrating high genotyping precision and clear differentiation between Bactrian and dromedary camels. With an enhanced genome coverage, accuracy and economic efficiency the Axiom_MaruPri SNP chip is poised to advance genomic breeding research in camels. It holds the potential to serve as an invaluable genetic resource for investigating population structure, genome-wide association studies and implementing genomic selection in domesticated camelid species.

DNA Microarray and Bioinformatic Analysis Reveals the Potential of Whale Oil in Enhancing Hair Growth in a C57BL/6 Mice Dorsal Skin Model.

Much research has been conducted to determine how hair regeneration is regulated, as this could provide therapeutic, cosmetic, and even psychological interventions for hair loss. The current study focused on the hair growth effect and effective utilization of fatty oil obtained from Bryde's whales through a high-throughput DNA microarray approach in conjunction with immunohistochemical observations. The research also examined the mechanisms and factors involved in hair growth. In an experiment using female C57BL/6J mice, the vehicle control group (VC: propylene glycol: ethanol: water), the positive control group (MXD: 3% minoxidil), and the experimental group (WO: 20% whale oil) were topically applied to the dorsal skin of the mouse. The results showed that 3% MXD and 20% WO were more effective than VC in promoting hair growth, especially 20% WO. Furthermore, in hematoxylin and eosin-stained dorsal skin tissue, an increase in the number of hair follicles and subcutaneous tissue thickness was observed with 20% WO. Whole-genome transcriptome analysis also confirmed increases for 20% WO in filaggrin (Flg), a gene related to skin barrier function; fibroblast growth factor 21 (Fgf21), which is involved in hair follicle development; and cysteine-rich secretory protein 1 (Crisp1), a candidate gene for alopecia areata. Furthermore, the results of KEGG pathway analysis indicated that 20% WO may have lower stress and inflammatory responses than 3% MXD. Therefore, WO is expected to be a safe hair growth agent.

Effect of genotyping errors on linkage map construction based on repeated chip analysis of two recombinant inbred line populations in wheat (Triticum aestivum L.).

Linkage maps are essential for genetic mapping of phenotypic traits, gene map-based cloning, and marker-assisted selection in breeding applications. Construction of a high-quality saturated map requires high-quality genotypic data on a large number of molecular markers. Errors in genotyping cannot be completely avoided, no matter what platform is used. When genotyping error reaches a threshold level, it will seriously affect the accuracy of the constructed map and the reliability of consequent genetic studies. In this study, repeated genotyping of two recombinant inbred line (RIL) populations derived from crosses Yangxiaomai × Zhongyou 9507 and Jingshuang 16 × Bainong 64 was used to investigate the effect of genotyping errors on linkage map construction. Inconsistent data points between the two replications were regarded as genotyping errors, which were classified into three types. Genotyping errors were treated as missing values, and therefore the non-erroneous data set was generated. Firstly, linkage maps were constructed using the two replicates as well as the non-erroneous data set. Secondly, error correction methods implemented in software packages QTL IciMapping (EC) and Genotype-Corrector (GC) were applied to the two replicates. Linkage maps were therefore constructed based on the corrected genotypes and then compared with those from the non-erroneous data set. Simulation study was performed by considering different levels of genotyping errors to investigate the impact of errors and the accuracy of error correction methods. Results indicated that map length and marker order differed among the two replicates and the non-erroneous data sets in both RIL populations. For both actual and simulated populations, map length was expanded as the increase in error rate, and the correlation coefficient between linkage and physical maps became lower. Map quality can be improved by repeated genotyping and error correction algorithm. When it is impossible to genotype the whole mapping population repeatedly, 30% would be recommended in repeated genotyping. The EC method had a much lower false positive rate than did the GC method under different error rates. This study systematically expounded the impact of genotyping errors on linkage analysis, providing potential guidelines for improving the accuracy of linkage maps in the presence of genotyping errors.

Critical evaluation of the reliability of DNA methylation probes on the Illumina MethylationEPIC v1.0 BeadChip microarrays.

DNA methylation (DNAm) plays a crucial role in a number of complex diseases. However, the reliability of DNAm levels measured using Illumina arrays varies across different probes. Previous research primarily assessed probe reliability by comparing duplicate samples between the 450k-450k or 450k-EPIC platforms, with limited investigations on Illumina EPIC v1.0 arrays. We conducted a comprehensive assessment of the EPIC v1.0 array probe reliability using 69 blood DNA samples, each measured twice, generated by the Alzheimer's Disease Neuroimaging Initiative study. We observed higher reliability in probes with average methylation beta values of 0.2 to 0.8, and lower reliability in type I probes or those within the promoter and CpG island regions. Importantly, we found that probe reliability has significant implications in the analyses of Epigenome-wide Association Studies (EWAS). Higher reliability is associated with more consistent effect sizes in different studies, the identification of differentially methylated regions (DMRs) and methylation quantitative trait locus (mQTLs), and significant correlations with downstream gene expression. Moreover, blood DNAm measurements obtained from probes with higher reliability are more likely to show concordance with brain DNAm measurements. Our findings, which provide crucial reliability information for probes on the EPIC v1.0 array, will serve as a valuable resource for future DNAm studies.

Robustness in population-structure and demographic-inference results derived from the Aedes aegypti genotyping chip and whole-genome sequencing data.

The mosquito Aedes aegypti is the primary vector of many human arboviruses such as dengue, yellow fever, chikungunya, and Zika, which affect millions of people worldwide. Population genetic studies on this mosquito have been important in understanding its invasion pathways and success as a vector of human disease. The Axiom aegypti1 SNP chip was developed from a sample of geographically diverse A. aegypti populations to facilitate genomic studies on this species. We evaluate the utility of the Axiom aegypti1 SNP chip for population genetics and compare it with a low-depth shotgun sequencing approach using mosquitoes from the native (Africa) and invasive ranges (outside Africa). These analyses indicate that results from the SNP chip are highly reproducible and have a higher sensitivity to capture alternative alleles than a low-coverage whole-genome sequencing approach. Although the SNP chip suffers from ascertainment bias, results from population structure, ancestry, demographic, and phylogenetic analyses using the SNP chip were congruent with those derived from low-coverage whole-genome sequencing, and consistent with previous reports on Africa and outside Africa populations using microsatellites. More importantly, we identified a subset of SNPs that can be reliably used to generate merged databases, opening the door to combined analyses. We conclude that the Axiom aegypti1 SNP chip is a convenient, more accurate, low-cost alternative to low-depth whole-genome sequencing for population genetic studies of A. aegypti that do not rely on full allelic frequency spectra. Whole-genome sequencing and SNP chip data can be easily merged, extending the usefulness of both approaches.

Theoretical Studies of DNA Microarray Present Potential Molecular and Cellular Interconnectivity of Signaling Pathways in Immune System Dysregulation.

Autoimmunity is defined as the inability to regulate immunological activities in the body, especially in response to external triggers, leading to the attack of the tissues and organs of the host. Outcomes include the onset of autoimmune diseases whose effects are primarily due to dysregulated immune responses. In past years, there have been cases that show an increased susceptibility to other autoimmune disorders in patients who are already experiencing the same type of disease. Research in this field has started analyzing the potential molecular and cellular causes of this interconnectedness, bearing in mind the possibility of advancing drugs and therapies for the treatment of autoimmunity. With that, this study aimed to determine the correlation of four autoimmune diseases, which are type 1 diabetes (T1D), psoriasis (PSR), systemic sclerosis (SSc), and systemic lupus erythematosus (SLE), by identifying highly preserved co-expressed genes among datasets using WGCNA. Functional annotation was then employed to characterize these sets of genes based on their systemic relationship as a whole to elucidate the biological processes, cellular components, and molecular functions of the pathways they are involved in. Lastly, drug repurposing analysis was performed to screen candidate drugs for repositioning that could regulate the abnormal expression of genes among the diseases. A total of thirteen modules were obtained from the analysis, the majority of which were associated with transcriptional, post-transcriptional, and post-translational modification processes. Also, the evaluation based on KEGG suggested the possible role of TH17 differentiation in the simultaneous onset of the four diseases. Furthermore, clomiphene was the top drug candidate for regulating overexpressed hub genes; meanwhile, prilocaine was the top drug for regulating under-expressed hub genes. This study was geared towards utilizing transcriptomics approaches for the assessment of microarray data, which is different from the use of traditional genomic analyses. Such a research design for investigating correlations among autoimmune diseases may be the first of its kind.

A rapid and low-cost platform for detection of bacterial based on microchamber PCR microfluidic chip.

Polymerase chain reaction (PCR) has been considered as the gold standard for detecting nucleic acids. The simple PCR system is of great significance for medical applications in remote areas, especially for the developing countries. Herein, we proposed a low-cost self-assembled platform for microchamber PCR. The working principle is rotating the chamber PCR microfluidic chip between two heaters with fixed temperature to solve the problem of low temperature variation rate. The system consists of two temperature controllers, a screw slide rail, a chamber array microfluidic chip and a self-built software. Such a system can be constructed at a cost of about US$60. The micro chamber PCR can be finished by rotating the microfluidic chip between two heaters with fixed temperature. Results demonstrated that the sensitivity of the temperature controller is 0.1℃. The relative error of the duration for the microfluidic chip was 0.02 s. Finally, we successfully finished amplification of the target gene of Porphyromonas gingivalis in the chamber PCR microfluidic chip within 35 min and on-site detection of its PCR products by fluorescence. The chip consisted of 3200 cylindrical chambers. The volume of reagent in each volume is as low as 0.628 nL. This work provides an effective method to reduce the amplification time required for micro chamber PCR.

Ultrasimple size encoded microfluidic chip for rapid simultaneous multiplex detection of DNA sequences.

Simultaneous multiplexed analysis can provide comprehensive information for disease diagnosis. However, the current multiplex methods rely on sophisticated barcode technology, which hinders its wider application. In this study, an ultrasimple size encoding method is proposed for multiplex detection using a wedge-shaped microfluidic chip. Driving by negative pressure, microparticles are naturally arranged in distinct stripes based on their sizes within the chip. This size encoding method demonstrates a high level of precision, allowing for accuracy in distinguishing 3-5 sizes of microparticles with a remarkable accuracy rate of up to 99%, even the microparticles with a size difference as small as 0.5 µm. The entire size encoding process is completed in less than 5 min, making it ultrasimple, reliable, and easy to operate. To evaluate the function of this size encoding microfluidic chip, three commonly co-infectious viruses' nucleic acid sequences (including complementary DNA sequences of HIV and HCV, and DNA sequence of HBV) are employed for multiplex detection. Results indicate that all three DNA sequences can be sensitively detected without any cross-interference. This size-encoding microfluidic chip-based multiplex detection method is simple, rapid, and high-resolution, its successful application in serum samples renders it highly promising for potential clinical promotion.

Biporous silica nanostructure-induced nanovortex in microfluidics for nucleic acid enrichment, isolation, and PCR-free detection.

Efficient pathogen enrichment and nucleic acid isolation are critical for accurate and sensitive diagnosis of infectious diseases, especially those with low pathogen levels. Our study introduces a biporous silica nanofilms-embedded sample preparation chip for pathogen and nucleic acid enrichment/isolation. This chip features unique biporous nanostructures comprising large and small pore layers. Computational simulations confirm that these nanostructures enhance the surface area and promote the formation of nanovortex, resulting in improved capture efficiency. Notably, the chip demonstrates a 100-fold lower limit of detection compared to conventional methods used for nucleic acid detection. Clinical validations using patient samples corroborate the superior sensitivity of the chip when combined with the luminescence resonance energy transfer assay. The enhanced sample preparation efficiency of the chip, along with the facile and straightforward synthesis of the biporous nanostructures, offers a promising solution for polymer chain reaction-free detection of nucleic acids.

Microarray-Based DNA Methylation Profiling: Validation Considerations for Clinical Testing.

Microarray-based methylation profiling has emerged as a valuable tool for refining diagnoses and revealing novel tumor subtypes, particularly in central nervous system tumors. Despite the increasing adoption of this technique in clinical genomic laboratories, no technical standards have been published in establishing minimum criteria for test validation. A working group with experience and expertise in DNA-based methylation profiling tests on central nervous system tumors collaborated to develop practical discussion points and focus on important considerations for validating this test in clinical laboratory settings. The experience in validating this methodology in a clinical setting is summarized. Specifically, the advantages and challenges associated with utilizing an in-house classifier compared with a third-party classifier are highlighted. Additionally, experiences in demonstrating the assay's sensitivity and specificity, establishing minimum sample criteria, and implementing quality control metrics are described. As methylation profiling for tumor classification expands to other tumor types and continues to evolve for various other applications, the critical considerations described here are expected to serve as a guidance for future efforts in establishing professional guidelines for this assay.

Low-input and single-cell methods for Infinium DNA methylation BeadChips.

The Infinium BeadChip is the most widely used DNA methylome assay technology for population-scale epigenome profiling. However, the standard workflow requires over 200 ng of input DNA, hindering its application to small cell-number samples, such as primordial germ cells. We developed experimental and analysis workflows to extend this technology to suboptimal input DNA conditions, including ultra-low input down to single cells. DNA preamplification significantly enhanced detection rates to over 50% in five-cell samples and ∼25% in single cells. Enzymatic conversion also substantially improved data quality. Computationally, we developed a method to model the background signal's influence on the DNA methylation level readings. The modified detection P-value calculation achieved higher sensitivities for low-input datasets and was validated in over 100 000 public diverse methylome profiles. We employed the optimized workflow to query the demethylation dynamics in mouse primordial germ cells available at low cell numbers. Our data revealed nuanced chromatin states, sex disparities, and the role of DNA methylation in transposable element regulation during germ cell development. Collectively, we present comprehensive experimental and computational solutions to extend this widely used methylation assay technology to applications with limited DNA.

Population characteristic exploitation-based multi-orientation multi-objective gene selection for microarray data classification.

Gene selection is a process of selecting discriminative genes from microarray data that helps to diagnose and classify cancer samples effectively. Swarm intelligence evolution-based gene selection algorithms can never circumvent the problem that the population is prone to local optima in the process of gene selection. To tackle this challenge, previous research has focused primarily on two aspects: mitigating premature convergence to local optima and escaping from local optima. In contrast to these strategies, this paper introduces a novel perspective by adopting reverse thinking, where the issue of local optima is seen as an opportunity rather than an obstacle. Building on this foundation, we propose MOMOGS-PCE, a novel gene selection approach that effectively exploits the advantageous characteristics of populations trapped in local optima to uncover global optimal solutions. Specifically, MOMOGS-PCE employs a novel population initialization strategy, which involves the initialization of multiple populations that explore diverse orientations to foster distinct population characteristics. The subsequent step involved the utilization of an enhanced NSGA-II algorithm to amplify the advantageous characteristics exhibited by the population. Finally, a novel exchange strategy is proposed to facilitate the transfer of characteristics between populations that have reached near maturity in evolution, thereby promoting further population evolution and enhancing the search for more optimal gene subsets. The experimental results demonstrated that MOMOGS-PCE exhibited significant advantages in comprehensive indicators compared with six competitive multi-objective gene selection algorithms. It is confirmed that the "reverse-thinking" approach not only avoids local optima but also leverages it to uncover superior gene subsets for cancer diagnosis.

MESBC: A novel mutually exclusive spectral biclustering method for cancer subtyping.

Many soft biclustering algorithms have been developed and applied to various biological and biomedical data analyses. However, few mutually exclusive (hard) biclustering algorithms have been proposed, which could better identify disease or molecular subtypes with survival significance based on genomic or transcriptomic data. In this study, we developed a novel mutually exclusive spectral biclustering (MESBC) algorithm based on spectral method to detect mutually exclusive biclusters. MESBC simultaneously detects relevant features (genes) and corresponding conditions (patients) subgroups and, therefore, automatically uses the signature features for each subtype to perform the clustering. Extensive simulations revealed that MESBC provided superior accuracy in detecting pre-specified biclusters compared with the non-negative matrix factorization (NMF) and Dhillon's algorithm, particularly in very noisy data. Further analysis of the algorithm on real datasets obtained from the TCGA database showed that MESBC provided more accurate (i.e., smaller p-value) overall survival prediction in patients with lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC) cancers when compared to the existing, gold-standard subtypes for lung cancers (integrative clustering). Furthermore, MESBC detected several genes with significant prognostic value in both LUAD and LUSC patients. External validation on an independent, unseen GEO dataset of LUAD showed that MESBC-derived clusters based on TCGA data still exhibited clear biclustering patterns and consistent, outstanding prognostic predictability, demonstrating robust generalizability of MESBC. Therefore, MESBC could potentially be used as a risk stratification tool to optimize the treatment for the patient, improve the selection of patients for clinical trials, and contribute to the development of novel therapeutic agents.