T-S2Inet: Transformer-based sequence-to-image network for accurate nanopore sequence recognition.

MOTIVATION: Nanopore sequencing is a new macromolecular recognition and perception technology that enables high-throughput sequencing of DNA, RNA, even protein molecules. The sequences generated by nanopore sequencing span a large time frame, and the labor and time costs incurred by traditional analysis methods are substantial. Recently, research on nanopore data analysis using machine learning algorithms has gained unceasing momentum, but there is often a significant gap between traditional and deep learning methods in terms of classification results. To analyze nanopore data using deep learning technologies, measures such as sequence completion and sequence transformation can be employed. However, these technologies do not preserve the local features of the sequences. To address this issue, we propose a sequence-to-image (S2I) module that transforms sequences of unequal length into images. Additionally, we propose the Transformer-based T-S2Inet model to capture the important information and improve the classification accuracy. RESULTS: Quantitative and qualitative analysis shows that the experimental results have an improvement of around 2% in accuracy compared to previous methods. The proposed method is adaptable to other nanopore platforms, such as the Oxford nanopore. It is worth noting that the proposed method not only aims to achieve the most advanced performance, but also provides a general idea for the analysis of nanopore sequences of unequal length. AVAILABILITY AND IMPLEMENTATION: The main program is available at https://github.com/guanxiaoyu11/S2Inet.

Automated individual cortical parcellation via consensus graph representation learning.

Cortical parcellation plays a pivotal role in elucidating the brain organization. Despite the growing efforts to develop parcellation algorithms using functional magnetic resonance imaging, achieving a balance between intra-individual specificity and inter-individual consistency proves challenging, making the generation of high-quality, subject-consistent cortical parcellations particularly elusive. To solve this problem, our paper proposes a fully automated individual cortical parcellation method based on consensus graph representation learning. The method integrates spectral embedding with low-rank tensor learning into a unified optimization model, which uses group-common connectivity patterns captured by low-rank tensor learning to optimize subjects' functional networks. This not only ensures consistency in brain representations across different subjects but also enhances the quality of each subject's representation matrix by eliminating spurious connections. More importantly, it achieves an adaptive balance between intra-individual specificity and inter-individual consistency during this process. Experiments conducted on a test-retest dataset from the Human Connectome Project (HCP) demonstrate that our method outperforms existing methods in terms of reproducibility, functional homogeneity, and alignment with task activation. Extensive network-based comparisons on the HCP S900 dataset reveal that the functional network derived from our cortical parcellation method exhibits greater capabilities in gender identification and behavior prediction than other approaches.

Dynamic cerebellar and sensorimotor network compensation in tremor-dominated Parkinson's disease.

AIM: Parkinson's disease (PD) tremor is associated with dysfunction in the basal ganglia (BG), cerebellum (CB), and sensorimotor networks (SMN). We investigated tremor-related static functional network connectivity (SFNC) and dynamic functional network connectivity (DFNC) in PD patients. METHODS: We analyzed the resting-state functional MRI data of 21 tremor-dominant Parkinson's disease (TDPD) patients and 29 healthy controls. We compared DFNC and SFNC between the three networks and assessed their associations with tremor severity. RESULTS: TDPD patients exhibited increased SFNC between the SMN and BG networks. In addition, they spent more mean dwell time (MDT) in state 2, characterized by sparse connections, and less MDT in state 4, indicating stronger connections. Furthermore, enhanced DFNC between the CB and SMN was observed in state 2. Notably, the MDT of state 2 was positively associated with tremor scores. CONCLUSION: The enhanced dynamic connectivity between the CB and SMN in TDPD patients suggests a potential compensatory mechanism. However, the tendency to remain in a state of sparse connectivity may contribute to the severity of tremor symptoms.

Constructing fine-grained spatiotemporal neonatal functional atlases with spectral functional network learning.

The early stages of human development are increasingly acknowledged as pivotal in laying the groundwork for subsequent behavioral and cognitive development. Spatiotemporal (4D) brain functional atlases are important in elucidating the development of human brain functions. However, the scarcity of such atlases for early life stages stems from two primary challenges: (1) the significant noise in functional magnetic resonance imaging (fMRI) that complicates the generation of high-quality atlases for each age group, and (2) the rapid and complex changes in the early human brain that hinder the maintenance of temporal consistency in 4D atlases. This study tackles these challenges by integrating low-rank tensor learning with spectral embedding, thereby proposing a novel, data-driven 4D functional atlas generation framework based on spectral functional network learning (SFNL). This method utilizes low-rank tensor learning to capture common functional connectivity (FC) patterns across different ages, thus optimizing FCs for each age group to improve the temporal consistency of functional networks. Incorporating spectral embedding aids in mitigating potential noise in FC networks derived from fMRI data by reconstructing networks in the spectral space. Utilizing SFNL-generated functional networks enables the creation of consistent and highly qualified spatiotemporal functional atlases. The framework was applied to the developing Human Connectome Project (dHCP) dataset, generating the first neonatal 4D functional atlases with fine-grained temporal and spatial resolutions. Experimental evaluations focusing on functional homogeneity, reliability, and temporal consistency demonstrated the superiority of our framework compared to existing methods for constructing 4D atlases. Additionally, network analysis experiments, including individual identification, functional systems development, and local efficiency assessments, further corroborate the efficacy and robustness of the generated atlases. The 4D atlases and related codes will be made publicly accessible (https://github.com/zhaoyunxi/neonate-atlases).

Cross-cohort replicable resting-state functional connectivity in predicting symptoms and cognition of schizophrenia.

Schizophrenia (SZ) is a debilitating mental illness characterized by adolescence or early adulthood onset of psychosis, positive and negative symptoms, as well as cognitive impairments. Despite a plethora of studies leveraging functional connectivity (FC) from functional magnetic resonance imaging (fMRI) to predict symptoms and cognitive impairments of SZ, the findings have exhibited great heterogeneity. We aimed to identify congruous and replicable connectivity patterns capable of predicting positive and negative symptoms as well as cognitive impairments in SZ. Predictable functional connections (FCs) were identified by employing an individualized prediction model, whose replicability was further evaluated across three independent cohorts (BSNIP, SZ = 174; COBRE, SZ = 100; FBIRN, SZ = 161). Across cohorts, we observed that altered FCs in frontal-temporal-cingulate-thalamic network were replicable in prediction of positive symptoms, while sensorimotor network was predictive of negative symptoms. Temporal-parahippocampal network was consistently identified to be associated with reduced cognitive function. These replicable 23 FCs effectively distinguished SZ from healthy controls (HC) across three cohorts (82.7%, 90.2%, and 86.1%). Furthermore, models built using these replicable FCs showed comparable accuracies to those built using the whole-brain features in predicting symptoms/cognition of SZ across the three cohorts (r = .17-.33, p < .05). Overall, our findings provide new insights into the neural underpinnings of SZ symptoms/cognition and offer potential targets for further research and possible clinical interventions.

S2Snet: deep learning for low molecular weight RNA identification with nanopore.

Ribonucleic acid (RNA) is a pivotal nucleic acid that plays a crucial role in regulating many biological activities. Recently, one study utilized a machine learning algorithm to automatically classify RNA structural events generated by a Mycobacterium smegmatis porin A nanopore trap. Although it can achieve desirable classification results, compared with deep learning (DL) methods, this classic machine learning requires domain knowledge to manually extract features, which is sophisticated, labor-intensive and time-consuming. Meanwhile, the generated original RNA structural events are not strictly equal in length, which is incompatible with the input requirements of DL models. To alleviate this issue, we propose a sequence-to-sequence (S2S) module that transforms the unequal length sequence (UELS) to the equal length sequence. Furthermore, to automatically extract features from the RNA structural events, we propose a sequence-to-sequence neural network based on DL. In addition, we add an attention mechanism to capture vital information for classification, such as dwell time and blockage amplitude. Through quantitative and qualitative analysis, the experimental results have achieved about a 2% performance increase (accuracy) compared to the previous method. The proposed method can also be applied to other nanopore platforms, such as the famous Oxford nanopore. It is worth noting that the proposed method is not only aimed at pursuing state-of-the-art performance but also provides an overall idea to process nanopore data with UELS.

Active learning for efficient analysis of high-throughput nanopore data.

MOTIVATION: As the third-generation sequencing technology, nanopore sequencing has been used for high-throughput sequencing of DNA, RNA, and even proteins. Recently, many studies have begun to use machine learning technology to analyze the enormous data generated by nanopores. Unfortunately, the success of this technology is due to the extensive labeled data, which often suffer from enormous labor costs. Therefore, there is an urgent need for a novel technology that can not only rapidly analyze nanopore data with high-throughput, but also significantly reduce the cost of labeling. To achieve the above goals, we introduce active learning to alleviate the enormous labor costs by selecting the samples that need to be labeled. This work applies several advanced active learning technologies to the nanopore data, including the RNA classification dataset (RNA-CD) and the Oxford Nanopore Technologies barcode dataset (ONT-BD). Due to the complexity of the nanopore data (with noise sequence), the bias constraint is introduced to improve the sample selection strategy in active learning. Results: The experimental results show that for the same performance metric, 50% labeling amount can achieve the best baseline performance for ONT-BD, while only 15% labeling amount can achieve the best baseline performance for RNA-CD. Crucially, the experiments show that active learning technology can assist experts in labeling samples, and significantly reduce the labeling cost. Active learning can greatly reduce the dilemma of difficult labeling of high-capacity nanopore data. We hope active learning can be applied to other problems in nanopore sequence analysis. AVAILABILITY AND IMPLEMENTATION: The main program is available at https://github.com/guanxiaoyu11/AL-for-nanopore. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Identify connectome between genotypes and brain network phenotypes via deep self-reconstruction sparse canonical correlation analysis.

MOTIVATION: As a rising research topic, brain imaging genetics aims to investigate the potential genetic architecture of both brain structure and function. It should be noted that in the brain, not all variations are deservedly caused by genetic effect, and it is generally unknown which imaging phenotypes are promising for genetic analysis. RESULTS: In this work, genetic variants (i.e. the single nucleotide polymorphism, SNP) can be correlated with brain networks (i.e. quantitative trait, QT), so that the connectome (including the brain regions and connectivity features) of functional brain networks from the functional magnetic resonance imaging data is identified. Specifically, a connection matrix is firstly constructed, whose upper triangle elements are selected to be connectivity features. Then, the PageRank algorithm is exploited for estimating the importance of different brain regions as the brain region features. Finally, a deep self-reconstruction sparse canonical correlation analysis (DS-SCCA) method is developed for the identification of genetic associations with functional connectivity phenotypic markers. This approach is a regularized, deep extension, scalable multi-SNP-multi-QT method, which is well-suited for applying imaging genetic association analysis to the Alzheimer's Disease Neuroimaging Initiative datasets. It is further optimized by adopting a parametric approach, augmented Lagrange and stochastic gradient descent. Extensive experiments are provided to validate that the DS-SCCA approach realizes strong associations and discovers functional connectivity and brain region phenotypic biomarkers to guide disease interpretation. AVAILABILITY AND IMPLEMENTATION: The Matlab code is available at https://github.com/meimeiling/DS-SCCA/tree/main. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Three-way parallel group independent component analysis: Fusion of spatial and spatiotemporal magnetic resonance imaging data.

Advances in imaging acquisition techniques allow multiple imaging modalities to be collected from the same subject. Each individual modality offers limited yet unique views of the functional, structural, or dynamic temporal features of the brain. Multimodal fusion provides effective ways to leverage these complementary perspectives from multiple modalities. However, the majority of current multimodal fusion approaches involving functional magnetic resonance imaging (fMRI) are limited to 3D feature summaries that do not incorporate its rich temporal information. Thus, we propose a novel three-way parallel group independent component analysis (pGICA) fusion method that incorporates the first-level 4D fMRI data (temporal information included) by parallelizing group ICA into parallel ICA via a unified optimization framework. A new variability matrix was defined to capture subject-wise functional variability and then link it to the mixing matrices of the other two modalities. Simulation results show that the three-way pGICA provides highly accurate cross-modality linkage estimation under both weakly and strongly correlated conditions, as well as comparable source estimation under different noise levels. Results using real brain imaging data identified one linked functional-structural-diffusion component associated to differences between schizophrenia and controls. This was replicated in an independent cohort, and the identified components were also correlated with major cognitive domains. Functional network connectivity revealed visual-subcortical and default mode-cerebellum pairs that discriminate between schizophrenia and controls. Overall, both simulation and real data results support the use of three-way pGICA to identify multimodal spatiotemporal links and to pursue the study of brain disorders under a single unifying multimodal framework.

Links between electroconvulsive therapy responsive and cognitive impairment multimodal brain networks in late-life major depressive disorder.

BACKGROUND: Although electroconvulsive therapy (ECT) is an effective treatment for depression, ECT cognitive impairment remains a major concern. The neurobiological underpinnings and mechanisms underlying ECT antidepressant and cognitive impairment effects remain unknown. This investigation aims to identify ECT antidepressant-response and cognitive-impairment multimodal brain networks and assesses whether they are associated with the ECT-induced electric field (E-field) with an optimal pulse amplitude estimation. METHODS: A single site clinical trial focused on amplitude (600, 700, and 800 mA) included longitudinal multimodal imaging and clinical and cognitive assessments completed before and immediately after the ECT series (n = 54) for late-life depression. Another two independent validation cohorts (n = 84, n = 260) were included. Symptom and cognition were used as references to supervise fMRI and sMRI fusion to identify ECT antidepressant-response and cognitive-impairment multimodal brain networks. Correlations between ECT-induced E-field within these two networks and clinical and cognitive outcomes were calculated. An optimal pulse amplitude was estimated based on E-field within antidepressant-response and cognitive-impairment networks. RESULTS: Decreased function in the superior orbitofrontal cortex and caudate accompanied with increased volume in medial temporal cortex showed covarying functional and structural alterations in both antidepressant-response and cognitive-impairment networks. Volume increases in the hippocampal complex and thalamus were antidepressant-response specific, and functional decreases in the amygdala and hippocampal complex were cognitive-impairment specific, which were validated in two independent datasets. The E-field within these two networks showed an inverse relationship with HDRS reduction and cognitive impairment. The optimal E-filed range as [92.7-113.9] V/m was estimated to maximize antidepressant outcomes without compromising cognitive safety. CONCLUSIONS: The large degree of overlap between antidepressant-response and cognitive-impairment networks challenges parameter development focused on precise E-field dosing with new electrode placements. The determination of the optimal individualized ECT amplitude within the antidepressant and cognitive networks may improve the treatment benefit-risk ratio. TRIAL REGISTRATION: ClinicalTrials.gov Identifier: NCT02999269.

[Unsupervised deep learning for identifying the O 6-carboxymethyl guanine by nanopore sequencing].

O 6-carboxymethyl guanine(O 6-CMG) is a highly mutagenic alkylation product of DNA that causes gastrointestinal cancer in organisms. Existing studies used mutant Mycobacterium smegmatis porin A (MspA) nanopore assisted by Phi29 DNA polymerase to localize it. Recently, machine learning technology has been widely used in the analysis of nanopore sequencing data. But the machine learning always need a large number of data labels that have brought extra work burden to researchers, which greatly affects its practicability. Accordingly, this paper proposes a nano-Unsupervised-Deep-Learning method (nano-UDL) based on an unsupervised clustering algorithm to identify methylation events in nanopore data automatically. Specially, nano-UDL first uses the deep AutoEncoder to extract features from the nanopore dataset and then applies the MeanShift clustering algorithm to classify data. Besides, nano-UDL can extract the optimal features for clustering by joint optimizing the clustering loss and reconstruction loss. Experimental results demonstrate that nano-UDL has relatively accurate recognition accuracy on the O 6-CMG dataset and can accurately identify all sequence segments containing O 6-CMG. In order to further verify the robustness of nano-UDL, hyperparameter sensitivity verification and ablation experiments were carried out in this paper. Using machine learning to analyze nanopore data can effectively reduce the additional cost of manual data analysis, which is significant for many biological studies, including genome sequencing.

[Cross subject personality assessment based on electroencephalogram functional connectivity and domain adaptation].

The research shows that personality assessment can be achieved by regression model based on electroencephalogram (EEG). Most of existing researches use event-related potential or power spectral density for personality assessment, which can only represent the brain information of a single region. But some research shows that human cognition is more dependent on the interaction of brain regions. In addition, due to the distribution difference of EEG features among subjects, the trained regression model can not get accurate results of cross subject personality assessment. In order to solve the problem, this research proposes a personality assessment method based on EEG functional connectivity and domain adaption. This research collected EEG data from 45 normal people under different emotional pictures (positive, negative and neutral). Firstly, the coherence of 59 channels in 5 frequency bands was taken as the original feature set. Then the feature-based domain adaptation was used to map the feature to a new feature space. It can reduce the distribution difference between training and test set in the new feature space, so as to reduce the distribution difference between subjects. Finally, the support vector regression model was trained and tested based on the transformed feature set by leave-one-out cross-validation. What's more, this paper compared the methods used in previous researches. The results showed that the method proposed in this paper improved the performance of regression model and obtained better personality assessment results. This research provides a new method for personality assessment.

GACDN: generative adversarial feature completion and diagnosis network for COVID-19.

BACKGROUND: The outbreak of coronavirus disease 2019 (COVID-19) causes tens of million infection world-wide. Many machine learning methods have been proposed for the computer-aided diagnosis between COVID-19 and community-acquired pneumonia (CAP) from chest computed tomography (CT) images. Most of these methods utilized the location-specific handcrafted features based on the segmentation results to improve the diagnose performance. However, the prerequisite segmentation step is time-consuming and needs the intervention by lots of expert radiologists, which cannot be achieved in the areas with limited medical resources. METHODS: We propose a generative adversarial feature completion and diagnosis network (GACDN) that simultaneously generates handcrafted features by radiomic counterparts and makes accurate diagnoses based on both original and generated features. Specifically, we first calculate the radiomic features from the CT images. Then, in order to fast obtain the location-specific handcrafted features, we use the proposed GACDN to generate them by its corresponding radiomic features. Finally, we use both radiomic features and location-specific handcrafted features for COVID-19 diagnosis. RESULTS: For the performance of our generated location-specific handcrafted features, the results of four basic classifiers show that it has an average of 3.21% increase in diagnoses accuracy. Besides, the experimental results on COVID-19 dataset show that our proposed method achieved superior performance in COVID-19 vs. community acquired pneumonia (CAP) classification compared with the state-of-the-art methods. CONCLUSIONS: The proposed method significantly improves the diagnoses accuracy of COVID-19 vs. CAP in the condition of incomplete location-specific handcrafted features. Besides, it is also applicable in some regions lacking of expert radiologists and high-performance computing resources.

Discovering network phenotype between genetic risk factors and disease status via diagnosis-aligned multi-modality regression method in Alzheimer's disease.

MOTIVATION: Neuroimaging genetics is an emerging field to identify the associations between genetic variants [e.g. single-nucleotide polymorphisms (SNPs)] and quantitative traits (QTs) such as brain imaging phenotypes. However, most of the current studies focus only on the associations between brain structure imaging and genetic variants, while neglecting the connectivity information between brain regions. In addition, the brain itself is a complex network, and the higher-order interaction may contain useful information for the mechanistic understanding of diseases [i.e. Alzheimer's disease (AD)]. RESULTS: A general framework is proposed to exploit network voxel information and network connectivity information as intermediate traits that bridge genetic risk factors and disease status. Specifically, we first use the sparse representation (SR) model to build hyper-network to express the connectivity features of the brain. The network voxel node features and network connectivity edge features are extracted from the structural magnetic resonance imaging (sMRI) and resting-state functional magnetic resonance imaging (fMRI), respectively. Second, a diagnosis-aligned multi-modality regression method is adopted to fully explore the relationships among modalities of different subjects, which can help further mine the relation between the risk genetics and brain network features. In experiments, all methods are tested on the Alzheimer's Disease Neuroimaging Initiative (ADNI) database. The experimental results not only verify the effectiveness of our proposed framework but also discover some brain regions and connectivity features that are highly related to diseases. AVAILABILITY AND IMPLEMENTATION: The Matlab code is available at http://ibrain.nuaa.edu.cn/2018/list.htm.

Topological correction of infant white matter surfaces using anatomically constrained convolutional neural network.

Reconstruction of accurate cortical surfaces without topological errors (i.e., handles and holes) from infant brain MR images is very important in early brain development studies. However, infant brain MR images typically suffer extremely low tissue contrast and dynamic imaging appearance patterns. Thus, it is inevitable to have large amounts of topological errors in the segmented infant brain tissue images, which lead to inaccurately reconstructed cortical surfaces with topological errors. To address this issue, inspired by recent advances in deep learning, we propose an anatomically constrained network for topological correction on infant cortical surfaces. Specifically, in our method, we first locate regions of potential topological defects by leveraging a topology-preserving level set method. Then, we propose an anatomically constrained network to correct those candidate voxels in the located regions. Since infant cortical surfaces often contain large and complex handles or holes, it is difficult to completely correct all errors using one-shot correction. Therefore, we further enroll these two steps into an iterative framework to gradually correct large topological errors. To the best of our knowledge, this is the first work to introduce deep learning approach for topological correction of infant cortical surfaces. We compare our method with the state-of-the-art methods on both simulated topological errors and real topological errors in human infant brain MR images. Moreover, we also validate our method on the infant brain MR images of macaques. All experimental results show the superior performance of the proposed method.

Iron homeostasis in pregnancy and spontaneous abortion.

During pregnancy, iron requirements are increased to support maternal erythropoietic expansion and fetal growth and development. To meet these requirements, dietary iron absorption increases, and available iron stores are mobilized. These adjustments are thought to be in large part mediated by the iron-regulatory hormone hepcidin, which controls the concentrations of ferroportin, the sole exporter of iron into the extracellular fluid and blood plasma. Hepcidin regulation of iron availability during healthy and abnormal pregnancies is not well understood. In our cross-sectional study, we compared hepcidin, iron and hematological parameters between nonpregnant control women, healthy pregnant women in the first and second trimester, and women with spontaneous abortion in the first trimester. We found that in healthy pregnancy, hepcidin increased in the first trimester compared with nonpregnant women, but then decreased during the second trimester. The second trimester hepcidin levels decreased despite stable serum iron concentrations, suggesting active suppression of hepcidin, presumably to enhance iron availability as iron demand increases. In women with spontaneous abortion during the first trimester, hepcidin, serum iron, and ferritin concentrations were all increased compared with the first trimester healthy pregnancy. Although the specific mechanisms remain to be determined, our findings demonstrate that maternal hepcidin is regulated by signals related to the progression of pregnancy, and that pregnancy loss is associated with profound changes in maternal iron metabolism. These observations highlight the existence of fetoplacental signals that modulate maternal iron homeostasis.

Multi-Atlas Based Methods in Brain MR Image Segmentation.

Brain region-of-interesting (ROI) segmentation is an important prerequisite step for many computer-aid brain disease analyses. However, the human brain has the complicated anatomical structure. Meanwhile, the brain MR images often suffer from the low intensity contrast around the boundary of ROIs, large inter-subject variance and large inner-subject variance. To address these issues, many multi-atlas based segmentation methods are proposed for brain ROI segmentation in the last decade. In this paper, multi-atlas based methods for brain MR image segmentation were reviewed regarding several registration toolboxes which are widely used in the multi-atlas methods, conventional methods for label fusion, datasets that have been used for evaluating the multi-atlas methods, as well as the applications of multi-atlas based segmentation in clinical researches. We propose that incorporating the anatomical prior into the end-to-end deep learning architectures for brain ROI segmentation is an important direction in the future.

MiR-30a-5p inhibits osteosarcoma cell proliferation and migration by targeting FOXD1.

Despite a number of studies have emphasized the extensive role of microRNA (miRNA) in the development of multiple cancers, the role of miR-30a-5p in the progression of osteosarcoma (OS) and the underlying mechanism are still limited. We detected the expression level of MiR-30a-5p and forkhead box D1 (FOXD1) in Clinical OS specimens and found that miR-30a-5p was significantly decreased while FOXD1 was markedly increased. Dual luciferase assay confirmed that FOXD1 was directly regulated by miR-30a-5p. In vitro assay showed that inhibitior of FOXD1 suppressed cell proliferation, migration and invasion in MG63 and U2OS cells, while overexpression of FOXD1 promoted OS cell proliferation and migration. In vivo assay further showed the inhibition of tumor growth after knockdown of FOXD1. These results suggested that FOXD1 might play key roles in OS development and progression, and was negatively regulated by miR-30a-5p in OS.

Identification of associations between genotypes and longitudinal phenotypes via temporally-constrained group sparse canonical correlation analysis.

MOTIVATION: Neuroimaging genetics identifies the relationships between genetic variants (i.e., the single nucleotide polymorphisms) and brain imaging data to reveal the associations from genotypes to phenotypes. So far, most existing machine-learning approaches are widely used to detect the effective associations between genetic variants and brain imaging data at one time-point. However, those associations are based on static phenotypes and ignore the temporal dynamics of the phenotypical changes. The phenotypes across multiple time-points may exhibit temporal patterns that can be used to facilitate the understanding of the degenerative process. In this article, we propose a novel temporally constrained group sparse canonical correlation analysis (TGSCCA) framework to identify genetic associations with longitudinal phenotypic markers. RESULTS: The proposed TGSCCA method is able to capture the temporal changes in brain from longitudinal phenotypes by incorporating the fused penalty, which requires that the differences between two consecutive canonical weight vectors from adjacent time-points should be small. A new efficient optimization algorithm is designed to solve the objective function. Furthermore, we demonstrate the effectiveness of our algorithm on both synthetic and real data (i.e., the Alzheimer's Disease Neuroimaging Initiative cohort, including progressive mild cognitive impairment, stable MCI and Normal Control participants). In comparison with conventional SCCA, our proposed method can achieve strong associations and discover phenotypic biomarkers across multiple time-points to guide disease-progressive interpretation. AVAILABILITY AND IMPLEMENTATION: The Matlab code is available at https://sourceforge.net/projects/ibrain-cn/files/ . CONTACT: dqzhang@nuaa.edu.cn or shenli@iu.edu.