Prevalence and clinical characteristics of venous thromboembolism in patients with lung cancer: a systematic review and meta-analysis.

Background: The prevalence of venous thromboembolism (VTE) is high in patients with cancer and can often present as the first symptom of malignancy. Cancer-associated VTE is one of the most important risk factors contributing to cancer mortality, making its prevention and treatment critical for patients with lung cancer. Methods: We systematically searched for observational studies that estimated the prevalence of VTE in patients with lung cancer. A comprehensive search of electronic databases, including PubMed, EMBASE and Cochrane Library, was systematically conducted from database inception through January 21, 2022. The qualities of included studies were assessed in three domains, including patient selection, comparison, and results. Random effects meta-analyses of the prevalence of VTE in lung cancer were conducted using the metaprop procedure. Chi-square test and I 2 value were used to evaluate study heterogeneity. Results: Thirty-five studies involving 742,156 patients were considered eligible for this study. The pooled prevalence of VTE among patients with lung cancer was 5% (95% CI: 0.043-0.056, P = 0.000). The regional prevalence of VTE was 7% (95% CI: 0.06-0.08; I2 = 99.2%) in North America, 8% (95% CI: 0.06-0.10; I2 = 97.6%) in Asia, 6% (95% CI: 0.04-0.09; I2 = 95.9%) in Europe and 11% (95% CI: 0.07-0.15) in Australasia. Conclusions: The prevalence of lung cancer-related VTE is high and region-specific. These results of this review emphasize the importance of understanding the incidence of lung cancer-related VTE and provide argue for VTE screening of patients with lung cancer. Systematic Review Registration: https://www.crd.york.ac.uk/prospero/, identifier PROSPERO (CRD42022306400).

A streamlined approach for intelligent ship object detection using EL-YOLO algorithm.

Maritime objects frequently exhibit low-quality and insufficient feature information, particularly in complex maritime environments characterized by challenges such as small objects, waves, and reflections. This situation poses significant challenges to the development of reliable object detection including the strategies of loss function and the feature understanding capabilities in common YOLOv8 (You Only Look Once) detectors. Furthermore, the widespread adoption and unmanned operation of intelligent ships have generated increasing demands on the computational efficiency and cost of object detection hardware, necessitating the development of more lightweight network architectures. This study proposes the EL-YOLO (Efficient Lightweight You Only Look Once) algorithm based on YOLOv8, designed specifically for intelligent ship object detection. EL-YOLO incorporates novel features, including adequate wise IoU (AWIoU) for improved bounding box regression, shortcut multi-fuse neck (SMFN) for a comprehensive analysis of features, and greedy-driven filter pruning (GDFP) to achieve a streamlined and lightweight network design. The findings of this study demonstrate notable advancements in both detection accuracy and lightweight characteristics across diverse maritime scenarios. EL-YOLO exhibits superior performance in intelligent ship object detection using RGB cameras, showcasing a significant improvement compared to standard YOLOv8 models.

SCRN: Single-cell Gene Regulatory Network Identification in Alzheimer's Disease.

Alzheimer's disease (AD) is the most common neurodegenerative disease, and it consumes considerable medical resources with increasing number of patients every year. Mounting evidence show that the regulatory disruptions altering the intrinsic activity of genes in brain cells contribute to AD pathogenesis. To gain insights into the underlying gene regulation in AD, we proposed a graph learning method, Single-Cell based Regulatory Network (SCRN), to identify the regulatory mechanisms based on single-cell data. SCRN implements the γ-decaying heuristic link prediction based on graph neural networks and can identify reliable gene regulatory networks using locally closed subgraphs. In this work, we first performed UMAP dimension reduction analysis on single-cell RNA sequencing (scRNA-seq) data of AD and normal samples. Then we used SCRN to construct the gene regulatory network based on three well-recognized AD genes (APOE, CX3CR1, and P2RY12). Enrichment analysis of the regulatory network revealed significant pathways including NGF signaling, ERBB2 signaling, and hemostasis. These findings demonstrate the feasibility of using SCRN to uncover potential biomarkers and therapeutic targets related to AD.

Harmonic Wavelet Neural Network for Discovering Neuropathological Propagation Patterns in Alzheimer's Disease.

Emerging research indicates that the degenerative biomarkers associated with Alzheimer's disease (AD) exhibit a non-random distribution within the cerebral cortex, instead following the structural brain network. The alterations in brain networks occur much earlier than the onset of clinical symptoms, thereby affecting the progression of brain disease. In this context, the utilization of computational methods to ascertain the propagation patterns of neuropathological events would contribute to the comprehension of the pathophysiological mechanism involved in the evolution of AD. Despite the encouraging findings achieved by existing graph-based deep learning approaches in analyzing irregular graph data, their applications in identifying the spreading pathway of neuropathology are limited due to two disadvantages. They include (1) lack of a common brain network as an unbiased reference basis for group comparison, and (2) lack of an appropriate mechanism for the identification of propagation patterns. To this end, we propose a proof-of-concept harmonic wavelet neural network (HWNN) to predict the early stage of AD and localize disease-related significant wavelets, which can be used to characterize the spreading pathways of neuropathological events across the brain network. The extensive experiments constructed on both synthetic and real datasets demonstrate that our proposed method achieves superior performance in classification accuracy and statistical power of identifying propagation patterns, compared with other representative approaches.

Pathology steered stratification network for subtype identification in Alzheimer's disease.

BACKGROUND: Alzheimer's disease (AD) is a heterogeneous, multifactorial neurodegenerative disorder characterized by three neurobiological factors beta-amyloid, pathologic tau, and neurodegeneration. There are no effective treatments for AD at a late stage, urging for early detection and prevention. However, existing statistical inference approaches in neuroimaging studies of AD subtype identification do not take into account the pathological domain knowledge, which could lead to ill-posed results that are sometimes inconsistent with the essential neurological principles. PURPOSE: Integrating systems biology modeling with machine learning, the study aims to assist clinical AD prognosis by providing a subpopulation classification in accordance with essential biological principles, neurological patterns, and cognitive symptoms. METHODS: We propose a novel pathology steered stratification network (PSSN) that incorporates established domain knowledge in AD pathology through a reaction-diffusion model, where we consider non-linear interactions between major biomarkers and diffusion along the brain structural network. Trained on longitudinal multimodal neuroimaging data, the biological model predicts long-term evolution trajectories that capture individual characteristic progression pattern, filling in the gaps between sparse imaging data available. A deep predictive neural network is then built to exploit spatiotemporal dynamics, link neurological examinations with clinical profiles, and generate subtype assignment probability on an individual basis. We further identify an evolutionary disease graph to quantify subtype transition probabilities through extensive simulations. RESULTS: Our stratification achieves superior performance in both inter-cluster heterogeneity and intra-cluster homogeneity of various clinical scores. Applying our approach to enriched samples of aging populations, we identify six subtypes spanning AD spectrum, where each subtype exhibits a distinctive biomarker pattern that is consistent with its clinical outcome. CONCLUSIONS: The proposed PSSN (i) reduces neuroimage data to low-dimensional feature vectors, (ii) combines AT[N]-Net based on real pathological pathways, (iii) predicts long-term biomarker trajectories, and (iv) stratifies subjects into fine-grained subtypes with distinct neurological underpinnings. PSSN provides insights into pre-symptomatic diagnosis and practical guidance on clinical treatments, which may be further generalized to other neurodegenerative diseases.

Uncovering the System Vulnerability and Criticality of Human Brain Under Dynamical Neuropathological Events in Alzheimer's Disease.

BACKGROUND: Despite the striking efforts in investigating neurobiological factors behind the acquisition of amyloid-ß (A), protein tau (T), and neurodegeneration ([N]) biomarkers, the mechanistic pathways of how AT[N] biomarkers spreading throughout the brain remain elusive. OBJECTIVE: To disentangle the massive heterogeneities in Alzheimer's disease (AD) progressions and identify vulnerable/critical brain regions to AD pathology. METHODS: In this work, we characterized the interaction of AT[N] biomarkers and their propagation across brain networks using a novel bistable reaction-diffusion model, which allows us to establish a new systems biology underpinning of AD progression. We applied our model to large-scale longitudinal neuroimages from the ADNI database and studied the systematic vulnerability and criticality of brains. RESULTS: Our model yields long term prediction that is statistically significant linear correlated with temporal imaging data, produces clinically consistent risk prediction, and captures the Braak-like spreading pattern of AT[N] biomarkers in AD development. CONCLUSIONS: Our major findings include (i) tau is a stronger indicator of regional risk compared to amyloid, (ii) temporal lobe exhibits higher vulnerability to AD-related pathologies, (iii) proposed critical brain regions outperform hub nodes in transmitting disease factors across the brain, and (iv) comparing the spread of neuropathological burdens caused by amyloid-ß and tau diffusions, disruption of metabolic balance is the most determinant factor contributing to the initiation and progression of AD.

Learning pyramidal multi-scale harmonic wavelets for identifying the neuropathology propagation patterns of Alzheimer's disease.

Previous studies have established that neurodegenerative disease such as Alzheimer's disease (AD) is a disconnection syndrome, where the neuropathological burdens often propagate across the brain network to interfere with the structural and functional connections. In this context, identifying the propagation patterns of neuropathological burdens sheds new light on understanding the pathophysiological mechanism of AD progression. However, little attention has been paid to propagation pattern identification by fully considering the intrinsic properties of brain-network organization, which plays an important role in improving the interpretability of the identified propagation pathways. To this end, we propose a novel harmonic wavelet analysis approach to construct a set of region-specific pyramidal multi-scale harmonic wavelets, it allows us to characterize the propagation patterns of neuropathological burdens from multiple hierarchical modules across the brain network. Specifically, we first extract underlying hub nodes through a series of network centrality measurements on the common brain network reference generated from a population of minimum spanning tree (MST) brain networks. Then, we propose a manifold learning method to identify the region-specific pyramidal multi-scale harmonic wavelets corresponding to hub nodes by seamlessly integrating the hierarchically modular property of the brain network. We estimate the statistical power of our proposed harmonic wavelet analysis approach on synthetic data and large-scale neuroimaging data from ADNI. Compared with the other harmonic analysis techniques, our proposed method not only effectively predicts the early stage of AD but also provides a new window to capture the underlying hub nodes and the propagation pathways of neuropathological burdens in AD.

Estimating Outlier-Immunized Common Harmonic Waves for Brain Network Analyses on the Stiefel Manifold.

Since brain network organization is essentially governed by the harmonic waves derived from the Eigen-system of the underlying Laplacian matrix, discovering the harmonic-based alterations provides a new window to understand the pathogenic mechanism of Alzheimer's disease (AD) in a unified reference space. However, current reference (common harmonic waves) estimation studies over the individual harmonic waves are often sensitive to outliers, which are obtained by averaging the heterogenous individual brain networks. To address this challenge, we propose a novel manifold learning approach to identify a set of outlier-immunized common harmonic waves. The backbone of our framework is calculating the geometric median of all individual harmonic waves on the Stiefel manifold, instead of Fréchet mean, thus improving the robustness of learned common harmonic waves to the outliers. A manifold optimization scheme with theoretically guaranteed convergence is tailored to solve our method. The experimental results on synthetic data and real data demonstrate that the common harmonic waves learned by our approach are not only more robust to the outliers than the state-of-the-art methods, but also provide a putative imaging biomarker to predict the early stage of AD.

scENT for Revealing Gene Clusters From Single-Cell RNA-Seq Data.

Recently, the fast development of single-cell RNA-seq (scRNA-seq) techniques has enabled high-resolution transcriptomic statistical analysis of individual cells in heterogeneous tissues, which can help researchers to explore the relationship between genes and human diseases. The emerging scRNA-seq data results in new analysis methods aiming to identify cell-level clustering and annotations. However, there are few methods developed to gain insights into the gene-level clusters with biological significance. This study proposes a new deep learning-based framework, scENT (single cell gENe clusTer), to identify significant gene clusters from single-cell RNA-seq data. We started with clustering the scRNA-seq data into multiple optimal groups, followed by a gene set enrichment analysis to identify classes of over-represented genes. Considering high-dimensional data with extensive zeros and dropout issues, scENT integrates perturbation in the learning process of clustering scRNA-seq data to improve its robustness and performance. Experimental results show that scENT outperformed other benchmarking methods on simulation data. To validate the biological insights of scENT, we applied it to the public experimental scRNA-seq data profiled from patients with Alzheimer's disease and brain metastasis. scENT successfully identified novel functional gene clusters and associated functions, facilitating the discovery of prospective mechanisms and the understanding of related diseases.

Dynamic Cross-Task Representation Adaptation for Clinical Targets Co-Segmentation in CT Image-Guided Post-Prostatectomy Radiotherapy.

Adjuvant and salvage radiotherapy after radical prostatectomy requires precise delineations of prostate bed (PB), i.e., the clinical target volume, and surrounding organs at risk (OARs) to optimize radiotherapy planning. Segmenting PB is particularly challenging even for clinicians, e.g., from the planning computed tomography (CT) images, as it is an invisible/virtual target after the operative removal of the cancerous prostate gland. Very recently, a few deep learning-based methods have been proposed to automatically contour non-contrast PB by leveraging its spatial reliance on adjacent OARs (i.e., the bladder and rectum) with much more clear boundaries, mimicking the clinical workflow of experienced clinicians. Although achieving state-of-the-art results from both the clinical and technical aspects, these existing methods improperly ignore the gap between the hierarchical feature representations needed for segmenting those fundamentally different clinical targets (i.e., PB and OARs), which in turn limits their delineation accuracy. This paper proposes an asymmetric multi-task network integrating dynamic cross-task representation adaptation (i.e., DyAdapt) for accurate and efficient co-segmentation of PB and OARs in one-pass from CT images. In the learning-to-learn framework, the DyAdapt modules adaptively transfer the hierarchical feature representations from the source task of OARs segmentation to match up with the target (and more challenging) task of PB segmentation, conditioned on the dynamic inter-task associations learned from the learning states of the feed-forward path. On a real-patient dataset, our method led to state-of-the-art results of PB and OARs co-segmentation. Code is available at https://github.com/ladderlab-xjtu/DyAdapt.

Dual-Stream Subspace Clustering Network for revealing gene targets in Alzheimer's disease.

The rapid development of scRNA-seq technology in recent years has enabled us to capture high-throughput gene expression profiles at single-cell resolution, reveal the heterogeneity of complex cell populations, and greatly advance our understanding of the underlying mechanisms in human diseases. Traditional methods for gene co-expression clustering are limited to discovering effective gene groups in scRNA-seq data. In this paper, we propose a novel gene clustering method based on convolutional neural networks called Dual-Stream Subspace Clustering Network (DS-SCNet). DS-SCNet can accurately identify important gene clusters from large scales of single-cell RNA-seq data and provide useful information for downstream analysis. Based on the simulated datasets, DS-SCNet successfully clusters genes into different groups and outperforms mainstream gene clustering methods, such as DBSCAN and DESC, across different evaluation metrics. To explore the biological insights of our proposed method, we applied it to real scRNA-seq data of patients with Alzheimer's disease (AD). DS-SCNet analyzed the single-cell RNA-seq data with 10,850 genes, and accurately identified 8 optimal clusters from 6673 cells. Enrichment analysis of these gene clusters revealed functional signaling pathways including the ILS signaling, the Rho GTPase signaling, and hemostasis pathways. Further analysis of gene regulatory networks identified new hub genes such as ELF4 as important regulators of AD, which indicates that DS-SCNet contributes to the discovery and understanding of the pathogenesis in Alzheimer's disease.

A Graph Convolutional Multiple Instance Learning on a Hypersphere Manifold Approach for Diagnosing Chronic Obstructive Pulmonary Disease in CT Images.

Chronic obstructive pulmonary disease (COPD) is a prevalent chronic disease with high morbidity and mortality. The early diagnosis of COPD is vital for clinical treatment, which helps patients to have a better quality of life. Because COPD can be ascribed to chronic bronchitis and emphysema, lesions in a computed tomography (CT) image can present anywhere inside the lung with different types, shapes and sizes. Multiple instance learning (MIL) is an effective tool for solving COPD discrimination. In this study, a novel graph convolutional MIL with the adaptive additive margin loss (GCMIL-AAMS) approach is proposed to diagnose COPD by CT. Specifically, for those early stage patients, the selected instance-level features can be more discriminative if they were learned by our proposed graph convolution and pooling with self-attention mechanism. The AAMS loss can utilize the information of COPD severity on a hypersphere manifold by adaptively setting the angular margins to improve the performance, as the severity can be quantified as four grades by pulmonary function test. The results show that our proposed GCMIL-AAMS method provides superior discrimination and generalization abilities in COPD discrimination, with areas under a receiver operating characteristic curve (AUCs) of 0.960 ± 0.014 and 0.862 ± 0.010 in the test set and external testing set, respectively, in 5-fold stratified cross validation; moreover, it demonstrates that graph learning is applicable to MIL and suggests that MIL may be adaptable to graph learning.

Two-stage deep learning network-based few-view image reconstruction for parallel-beam projection tomography.

Background: Projection tomography (PT) is a very important and valuable method for fast volumetric imaging with isotropic spatial resolution. Sparse-view or limited-angle reconstruction-based PT can greatly reduce data acquisition time, lower radiation doses, and simplify sample fixation modes. However, few techniques can currently achieve image reconstruction based on few-view projection data, which is especially important for in vivo PT in living organisms. Methods: A 2-stage deep learning network (TSDLN)-based framework was proposed for parallel-beam PT reconstructions using few-view projections. The framework is composed of a reconstruction network (R-net) and a correction network (C-net). The R-net is a generative adversarial network (GAN) used to complete image information with direct back-projection (BP) of a sparse signal, bringing the reconstructed image close to reconstruction results obtained from fully projected data. The C-net is a U-net array that denoises the compensation result to obtain a high-quality reconstructed image. Results: The accuracy and feasibility of the proposed TSDLN-based framework in few-view projection-based reconstruction were first evaluated with simulations, using images from the DeepLesion public dataset. The framework exhibited better reconstruction performance than traditional analytic reconstruction algorithms and iterative algorithms, especially in cases using sparse-view projection images. For example, with as few as two projections, the TSDLN-based framework reconstructed high-quality images very close to the original image, with structural similarities greater than 0.8. By using previously acquired optical PT (OPT) data in the TSDLN-based framework trained on computed tomography (CT) data, we further exemplified the migration capabilities of the TSDLN-based framework. The results showed that when the number of projections was reduced to 5, the contours and distribution information of the samples in question could still be seen in the reconstructed images. Conclusions: The simulations and experimental results showed that the TSDLN-based framework has strong reconstruction abilities using few-view projection images, and has great potential in the application of in vivo PT.

Characterizing the propagation pathway of neuropathological events of Alzheimer's disease using harmonic wavelet analysis.

Empirical imaging biomarkers such as the level of the regional pathological burden are widely used to measure the risk of developing neurodegenerative diseases such as Alzheimer's disease (AD). However, ample evidence shows that the brain network (wirings of white matter fibers) plays a vital role in the progression of AD, where neuropathological burdens often propagate across the brain network in a prion-like manner. In this context, characterizing the spreading pathway of AD-related neuropathological events sheds new light on understanding the heterogeneity of pathophysiological mechanisms in AD. In this work, we propose a manifold-based harmonic network analysis approach to explore a novel imaging biomarker in the form of the AD propagation pattern, which eventually allows us to identify the AD-related spreading pathways of neuropathological events throughout the brain. The backbone of this new imaging biomarker is a set of region-adaptive harmonic wavelets that represent the common network topology across individuals. We conceptualize that the individual's brain network and its associated pathology pattern form a unique system, which vibrates as do all natural objects in the universe. Thus, we can computationally excite such a brain system using selected harmonic wavelets that match the system's resonance frequency, where the resulting oscillatory wave manifests the system-level propagation pattern of neuropathological events across the brain network. We evaluate the statistical power of our harmonic network analysis approach on large-scale neuroimaging data from ADNI. Compared with the other empirical biomarkers, our harmonic wavelets not only yield a new imaging biomarker to potentially predict the cognitive decline in the early stage but also offer a new window to capture the in-vivo spreading pathways of neuropathological burden with a rigorous mathematics insight.

Characterizing Network Selectiveness to the Dynamic Spreading of Neuropathological Events in Alzheimer's Disease.

BACKGROUND: Mounting evidence shows that the neuropathological burdens manifest preference in affecting brain regions during the dynamic progression of Alzheimer's disease (AD). Since the distinct brain regions are physically wired by white matter fibers, it is reasonable to hypothesize the differential spreading pattern of neuropathological burdens may underlie the wiring topology, which can be characterized using neuroimaging and network science technologies. OBJECTIVE: To study the dynamic spreading patterns of neuropathological events in AD. METHODS: We first examine whether hub nodes with high connectivity in the brain network (assemble of white matter wirings) are susceptible to a higher level of pathological burdens than other regions that are less involved in the process of information exchange in the network. Moreover, we propose a novel linear mixed-effect model to characterize the multi-factorial spreading process of neuropathological burdens from hub nodes to non-hub nodes, where age, sex, and APOE4 indicators are considered as confounders. We apply our statistical model to the longitudinal neuroimaging data of amyloid-PET and tau-PET, respectively. RESULTS: Our meta-data analysis results show that 1) AD differentially affects hub nodes with a significantly higher level of pathology, and 2) the longitudinal increase of neuropathological burdens on non-hub nodes is strongly correlated with the connectome distance to hub nodes rather than the spatial proximity. CONCLUSION: The spreading pathway of AD neuropathological burdens might start from hub regions and propagate through the white matter fibers in a prion-like manner.

Group-Wise Hub Identification by Learning Common Graph Embeddings on Grassmannian Manifold.

Human brain is a complex yet economically organized system, where a small portion of critical hub regions support the majority of brain functions. The identification of common hub nodes in a population of networks is often simplified as a voting procedure on the set of identified hub nodes across individual brain networks, which ignores the intrinsic data geometry and partially lacks the reproducible findings in neuroscience. Hence, we propose a first-ever group-wise hub identification method to identify hub nodes that are common across a population of individual brain networks. Specifically, the backbone of our method is to learn common graph embedding that can represent the majority of local topological profiles. By requiring orthogonality among the graph embedding vectors, each graph embedding as a data element is residing on the Grassmannian manifold. We present a novel Grassmannian manifold optimization scheme that allows us to find the common graph embeddings, which not only identify the most reliable hub nodes in each network but also yield a population-based common hub node map. Results of the accuracy and replicability on both synthetic and real network data show that the proposed manifold learning approach outperforms all hub identification methods employed in this evaluation.

Harnessing the Power of Hybrid Light Propagation Model for Three-Dimensional Optical Imaging in Cancer Detection.

Optical imaging is an emerging technology capable of qualitatively and quantitatively observing life processes at the cellular or molecular level and plays a significant role in cancer detection. In particular, to overcome the disadvantages of traditional optical imaging that only two-dimensionally and qualitatively detect biomedical information, the corresponding three-dimensional (3D) imaging technology is intensively explored to provide 3D quantitative information, such as localization and distribution and tumor cell volume. To retrieve these information, light propagation models that reflect the interaction between light and biological tissues are an important prerequisite and basis for 3D optical imaging. This review concentrates on the recent advances in hybrid light propagation models, with particular emphasis on their powerful use for 3D optical imaging in cancer detection. Finally, we prospect the wider application of the hybrid light propagation model and future potential of 3D optical imaging in cancer detection.

Joint hub identification for brain networks by multivariate graph inference.

Recent developments in neuroimaging allow us to investigate the structural and functional connectivity between brain regions in vivo. Mounting evidence suggests that hub nodes play a central role in brain communication and neural integration. Such high centrality, however, makes hub nodes particularly susceptible to pathological network alterations and the identification of hub nodes from brain networks has attracted much attention in neuroimaging. Current popular hub identification methods often work in a univariate manner, i.e., selecting the hub nodes one after another based on either heuristic of the connectivity profile at each node or predefined settings of network modules. Since the topological information of the entire network (such as network modules) is not fully utilized, current methods have limited power to identify hubs that link multiple modules (connector hubs) and are biased toward identifying hubs having many connections within the same module (provincial hubs). To address this challenge, we propose a novel multivariate hub identification method. Our method identifies connector hubs as those that partition the network into disconnected components when they are removed from the network. Furthermore, we extend our hub identification method to find the population-based hub nodes from a group of network data. We have compared our hub identification method with existing methods on both simulated and human brain network data. Our proposed method achieves more accurate and replicable discovery of hub nodes and exhibits enhanced statistical power in identifying network alterations related to neurological disorders such as Alzheimer's disease and obsessive-compulsive disorder.

Learning dynamic graph embeddings for accurate detection of cognitive state changes in functional brain networks.

Mounting evidence shows that brain functions and cognitive states are dynamically changing even in the resting state rather than remaining at a single constant state. Due to the relatively small changes in BOLD (blood-oxygen-level-dependent) signals across tasks, it is difficult to detect the change of cognitive status without requiring prior knowledge of the experimental design. To address this challenge, we present a dynamic graph learning approach to generate an ensemble of subject-specific dynamic graph embeddings, which allows us to use brain networks to disentangle cognitive events more accurately than using raw BOLD signals. The backbone of our method is essentially a representation learning process for projecting BOLD signals into a latent vertex-temporal domain with the greater biological underpinning of brain activities. Specifically, the learned representation domain is jointly formed by (1) a set of harmonic waves that govern the topology of whole-brain functional connectivities and (2) a set of Fourier bases that characterize the temporal dynamics of functional changes. In this regard, our dynamic graph embeddings provide a new methodology to investigate how these self-organized functional fluctuation patterns oscillate along with the evolving cognitive status. We have evaluated our proposed method on both simulated data and working memory task-based fMRI datasets, where our dynamic graph embeddings achieve higher accuracy in detecting multiple cognitive states than other state-of-the-art methods.