Alternative cleavage and polyadenylation generates downstream uncapped RNA isoforms with translation potential.

The use of alternative promoters, splicing, and cleavage and polyadenylation (APA) generates mRNA isoforms that expand the diversity and complexity of the transcriptome. Here, we uncovered thousands of previously undescribed 5' uncapped and polyadenylated transcripts (5' UPTs). We show that these transcripts resist exonucleases due to a highly structured RNA and N6-methyladenosine modification at their 5' termini. 5' UPTs appear downstream of APA sites within their host genes and are induced upon APA activation. Strong enrichment in polysomal RNA fractions indicates 5' UPT translational potential. Indeed, APA promotes downstream translation initiation, non-canonical protein output, and consistent changes to peptide presentation at the cell surface. Lastly, we demonstrate the biological importance of 5' UPTs using Bcl2, a prominent anti-apoptotic gene whose entire coding sequence is a 5' UPT generated from 5' UTR-embedded APA sites. Thus, APA is not only accountable for terminating transcripts, but also for generating downstream uncapped RNAs with translation potential and biological impact.

A powerful approach to identify replicable variants in genome-wide association studies.

Replicability is the cornerstone of modern scientific research. Reliable identifications of genotype-phenotype associations that are significant in multiple genome-wide association studies (GWASs) provide stronger evidence for the findings. Current replicability analysis relies on the independence assumption among single-nucleotide polymorphisms (SNPs) and ignores the linkage disequilibrium (LD) structure. We show that such a strategy may produce either overly liberal or overly conservative results in practice. We develop an efficient method, ReAD, to detect replicable SNPs associated with the phenotype from two GWASs accounting for the LD structure. The local dependence structure of SNPs across two heterogeneous studies is captured by a four-state hidden Markov model (HMM) built on two sequences of p values. By incorporating information from adjacent locations via the HMM, our approach provides more accurate SNP significance rankings. ReAD is scalable, platform independent, and more powerful than existing replicability analysis methods with effective false discovery rate control. Through analysis of datasets from two asthma GWASs and two ulcerative colitis GWASs, we show that ReAD can identify replicable genetic loci that existing methods might otherwise miss.

On the core segmentation algorithms of copy number variation detection tools.

Shotgun sequencing is a high-throughput method used to detect copy number variants (CNVs). Although there are numerous CNV detection tools based on shotgun sequencing, their quality varies significantly, leading to performance discrepancies. Therefore, we conducted a comprehensive analysis of next-generation sequencing-based CNV detection tools over the past decade. Our findings revealed that the majority of mainstream tools employ similar detection rationale: calculates the so-called read depth signal from aligned sequencing reads and then segments the signal by utilizing either circular binary segmentation (CBS) or hidden Markov model (HMM). Hence, we compared the performance of those two core segmentation algorithms in CNV detection, considering varying sequencing depths, segment lengths and complex types of CNVs. To ensure a fair comparison, we designed a parametrical model using mainstream statistical distributions, which allows for pre-excluding bias correction such as guanine-cytosine (GC) content during the preprocessing step. The results indicate the following key points: (1) Under ideal conditions, CBS demonstrates high precision, while HMM exhibits a high recall rate. (2) For practical conditions, HMM is advantageous at lower sequencing depths, while CBS is more competitive in detecting small variant segments compared to HMM. (3) In case involving complex CNVs resembling real sequencing, HMM demonstrates more robustness compared with CBS. (4) When facing large-scale sequencing data, HMM costs less time compared with the CBS, while their memory usage is approximately equal. This can provide an important guidance and reference for researchers to develop new tools for CNV detection.

Inference of Locus-Specific Population Mixtures from Linked Genome-Wide Allele Frequencies.

Admixture between populations and species is common in nature. Since the influx of new genetic material might be either facilitated or hindered by selection, variation in mixture proportions along the genome is expected in organisms undergoing recombination. Various graph-based models have been developed to better understand these evolutionary dynamics of population splits and mixtures. However, current models assume a single mixture rate for the entire genome and do not explicitly account for linkage. Here, we introduce TreeSwirl, a novel method for inferring branch lengths and locus-specific mixture proportions by using genome-wide allele frequency data, assuming that the admixture graph is known or has been inferred. TreeSwirl builds upon TreeMix that uses Gaussian processes to estimate the presence of gene flow between diverged populations. However, in contrast to TreeMix, our model infers locus-specific mixture proportions employing a hidden Markov model that accounts for linkage. Through simulated data, we demonstrate that TreeSwirl can accurately estimate locus-specific mixture proportions and handle complex demographic scenarios. It also outperforms related D- and f-statistics in terms of accuracy and sensitivity to detect introgressed loci.

Field model for multistate lateral diffusion of various transmembrane proteins observed in living Dictyostelium cells.

The lateral diffusion of transmembrane proteins on plasma membranes is a fundamental process for various cellular functions. Diffusion properties specific for individual protein species have been extensively studied, but the common features among protein species are poorly understood. Here, we systematically studied the lateral diffusion of various transmembrane proteins in the lower eukaryote Dictyostelium discoideum cells using a hidden Markov model for single-molecule trajectories obtained experimentally. As common features, all membrane proteins that had from one to ten transmembrane regions adopted three free diffusion states with similar diffusion coefficients regardless of their structural variability. All protein species reduced their mobility similarly upon the inhibition of microtubule or actin cytoskeleton dynamics, or myosin II. The relationship between protein size and the diffusion coefficient was consistent with the Saffman-Delbrück model, meaning that membrane viscosity is a major determinant of lateral diffusion, but protein size is not. These protein species-independent properties of multistate free diffusion were explained simply and quantitatively by free diffusion on the three membrane regions with different viscosities, which is in sharp contrast to the complex diffusion behavior of transmembrane proteins in higher eukaryotes.

rcCAE: a convolutional autoencoder method for detecting intra-tumor heterogeneity and single-cell copy number alterations.

Intra-tumor heterogeneity (ITH) is one of the major confounding factors that result in cancer relapse, and deciphering ITH is essential for personalized therapy. Single-cell DNA sequencing (scDNA-seq) now enables profiling of single-cell copy number alterations (CNAs) and thus aids in high-resolution inference of ITH. Here, we introduce an integrated framework called rcCAE to accurately infer cell subpopulations and single-cell CNAs from scDNA-seq data. A convolutional autoencoder (CAE) is employed in rcCAE to learn latent representation of the cells as well as distill copy number information from noisy read counts data. This unsupervised representation learning via the CAE model makes it convenient to accurately cluster cells over the low-dimensional latent space, and detect single-cell CNAs from enhanced read counts data. Extensive performance evaluations on simulated datasets show that rcCAE outperforms the existing CNA calling methods, and is highly effective in inferring clonal architecture. Furthermore, evaluations of rcCAE on two real datasets demonstrate that it is able to provide a more refined clonal structure, of which some details are lost in clonal inference based on integer copy numbers.

An FPGA-based hardware accelerator supporting sensitive sequence homology filtering with profile hidden Markov models.

BACKGROUND: Sequence alignment lies at the heart of genome sequence annotation. While the BLAST suite of alignment tools has long held an important role in alignment-based sequence database search, greater sensitivity is achieved through the use of profile hidden Markov models (pHMMs). Here, we describe an FPGA hardware accelerator, called HAVAC, that targets a key bottleneck step (SSV) in the analysis pipeline of the popular pHMM alignment tool, HMMER. RESULTS: The HAVAC kernel calculates the SSV matrix at 1739 GCUPS on a ∼  $3000 Xilinx Alveo U50 FPGA accelerator card, ∼  227× faster than the optimized SSV implementation in nhmmer. Accounting for PCI-e data transfer data processing, HAVAC is 65× faster than nhmmer's SSV with one thread and 35× faster than nhmmer with four threads, and uses ∼  31% the energy of a traditional high end Intel CPU. CONCLUSIONS: HAVAC demonstrates the potential offered by FPGA hardware accelerators to produce dramatic speed gains in sequence annotation and related bioinformatics applications. Because these computations are performed on a co-processor, the host CPU remains free to simultaneously compute other aspects of the analysis pipeline.

kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R.

BACKGROUND: Approximating the recent phylogeny of N phased haplotypes at a set of variants along the genome is a core problem in modern population genomics and central to performing genome-wide screens for association, selection, introgression, and other signals. The Li & Stephens (LS) model provides a simple yet powerful hidden Markov model for inferring the recent ancestry at a given variant, represented as an N × N distance matrix based on posterior decodings. RESULTS: We provide a high-performance engine to make these posterior decodings readily accessible with minimal pre-processing via an easy to use package kalis, in the statistical programming language R. kalis enables investigators to rapidly resolve the ancestry at loci of interest and developers to build a range of variant-specific ancestral inference pipelines on top. kalis exploits both multi-core parallelism and modern CPU vector instruction sets to enable scaling to hundreds of thousands of genomes. CONCLUSIONS: The resulting distance matrices accessible via kalis enable local ancestry, selection, and association studies in modern large scale genomic datasets.

Inference of genomic landscapes using ordered Hidden Markov Models with emission densities (oHMMed).

BACKGROUND: Genomes are inherently inhomogeneous, with features such as base composition, recombination, gene density, and gene expression varying along chromosomes. Evolutionary, biological, and biomedical analyses aim to quantify this variation, account for it during inference procedures, and ultimately determine the causal processes behind it. Since sequential observations along chromosomes are not independent, it is unsurprising that autocorrelation patterns have been observed e.g., in human base composition. In this article, we develop a class of Hidden Markov Models (HMMs) called oHMMed (ordered HMM with emission densities, the corresponding R package of the same name is available on CRAN): They identify the number of comparably homogeneous regions within autocorrelated observed sequences. These are modelled as discrete hidden states; the observed data points are realisations of continuous probability distributions with state-specific means that enable ordering of these distributions. The observed sequence is labelled according to the hidden states, permitting only neighbouring states that are also neighbours within the ordering of their associated distributions. The parameters that characterise these state-specific distributions are inferred. RESULTS: We apply our oHMMed algorithms to the proportion of G and C bases (modelled as a mixture of normal distributions) and the number of genes (modelled as a mixture of poisson-gamma distributions) in windows along the human, mouse, and fruit fly genomes. This results in a partitioning of the genomes into regions by statistically distinguishable averages of these features, and in a characterisation of their continuous patterns of variation. In regard to the genomic G and C proportion, this latter result distinguishes oHMMed from segmentation algorithms based in isochore or compositional domain theory. We further use oHMMed to conduct a detailed analysis of variation of chromatin accessibility (ATAC-seq) and epigenetic markers H3K27ac and H3K27me3 (modelled as a mixture of poisson-gamma distributions) along the human chromosome 1 and their correlations. CONCLUSIONS: Our algorithms provide a biologically assumption free approach to characterising genomic landscapes shaped by continuous, autocorrelated patterns of variation. Despite this, the resulting genome segmentation enables extraction of compositionally distinct regions for further downstream analyses.

Estimating Temporally Variable Selection Intensity from Ancient DNA Data.

Novel technologies for recovering DNA information from archaeological and historical specimens have made available an ever-increasing amount of temporally spaced genetic samples from natural populations. These genetic time series permit the direct assessment of patterns of temporal changes in allele frequencies and hold the promise of improving power for the inference of selection. Increased time resolution can further facilitate testing hypotheses regarding the drivers of past selection events such as the incidence of plant and animal domestication. However, studying past selection processes through ancient DNA (aDNA) still involves considerable obstacles such as postmortem damage, high fragmentation, low coverage, and small samples. To circumvent these challenges, we introduce a novel Bayesian framework for the inference of temporally variable selection based on genotype likelihoods instead of allele frequencies, thereby enabling us to model sample uncertainties resulting from the damage and fragmentation of aDNA molecules. Also, our approach permits the reconstruction of the underlying allele frequency trajectories of the population through time, which allows for a better understanding of the drivers of selection. We evaluate its performance through extensive simulations and demonstrate its utility with an application to the ancient horse samples genotyped at the loci for coat coloration. Our results reveal that incorporating sample uncertainties can further improve the inference of selection.

Post-task responses following working memory and movement are driven by transient spectral bursts with similar characteristics.

The post-movement beta rebound has been studied extensively using magnetoencephalography (MEG) and is reliably modulated by various task parameters as well as illness. Our recent study showed that rebounds, which we generalise as "post-task responses" (PTRs), are a ubiquitous phenomenon in the brain, occurring across the cortex in theta, alpha, and beta bands. Currently, it is unknown whether PTRs following working memory are driven by transient bursts, which are moments of short-lived high amplitude activity, similar to those that drive the post-movement beta rebound. Here, we use three-state univariate hidden Markov models (HMMs), which can identify bursts without a priori knowledge of frequency content or response timings, to compare bursts that drive PTRs in working memory and visuomotor MEG datasets. Our results show that PTRs across working memory and visuomotor tasks are driven by pan-spectral transient bursts. These bursts have very similar spectral content variation over the cortex, correlating strongly between the two tasks in the alpha (R2 = .89) and beta (R2 = .53) bands. Bursts also have similar variation in duration over the cortex (e.g., long duration bursts occur in the motor cortex for both tasks), strongly correlating over cortical regions between tasks (R2 = .56), with a mean over all regions of around 300 ms in both datasets. Finally, we demonstrate the ability of HMMs to isolate signals of interest in MEG data, such that the HMM probability timecourse correlates more strongly with reaction times than frequency filtered power envelopes from the same brain regions. Overall, we show that induced PTRs across different tasks are driven by bursts with similar characteristics, which can be identified using HMMs. Given the similarity between bursts across tasks, we suggest that PTRs across the cortex may be driven by a common underlying neural phenomenon.

Brain state dynamics differ between eyes open and eyes closed rest.

The human brain exhibits spatio-temporally complex activity even in the absence of external stimuli, cycling through recurring patterns of activity known as brain states. Thus far, brain state analysis has primarily been restricted to unimodal neuroimaging data sets, resulting in a limited definition of state and a poor understanding of the spatial and temporal relationships between states identified from different modalities. Here, we applied hidden Markov model (HMM) to concurrent electroencephalography-functional magnetic resonance imaging (EEG-fMRI) eyes open (EO) and eyes closed (EC) resting-state data, training models on the EEG and fMRI data separately, and evaluated the models' ability to distinguish dynamics between the two rest conditions. Additionally, we employed a general linear model approach to identify the BOLD correlates of the EEG-defined states to investigate whether the fMRI data could be used to improve the spatial definition of the EEG states. Finally, we performed a sliding window-based analysis on the state time courses to identify slower changes in the temporal dynamics, and then correlated these time courses across modalities. We found that both models could identify expected changes during EC rest compared to EO rest, with the fMRI model identifying changes in the activity and functional connectivity of visual and attention resting-state networks, while the EEG model correctly identified the canonical increase in alpha upon eye closure. In addition, by using the fMRI data, it was possible to infer the spatial properties of the EEG states, resulting in BOLD correlation maps resembling canonical alpha-BOLD correlations. Finally, the sliding window analysis revealed unique fractional occupancy dynamics for states from both models, with a selection of states showing strong temporal correlations across modalities. Overall, this study highlights the efficacy of using HMMs for brain state analysis, confirms that multimodal data can be used to provide more in-depth definitions of state and demonstrates that states defined across different modalities show similar temporal dynamics.

A novel starch-active lytic polysaccharide monooxygenase discovered with bioinformatics screening and its application in textile desizing.

BACKGROUND: Lytic polysaccharide monooxygenases (LPMOs) catalyzing the oxidative cleavage of different types of polysaccharides have potential to be used in various industries. However, AA13 family LPMOs which specifically catalyze starch substrates have relatively less members than AA9 and AA10 families to limit their application range. Amylase has been used in enzymatic desizing treatment of cotton fabric for semicentury which urgently need for new assistant enzymes to improve reaction efficiency and reduce cost so as to promote their application in the textile industry. RESULTS: A total of 380 unannotated new genes which probably encode AA13 family LPMOs were discovered by the Hidden Markov model scanning in this study. Ten of them have been successfully heterologous overexpressed. AlLPMO13 with the highest activity has been purified and determined its optimum pH and temperature as pH 5.0 and 50 °C. It also showed various oxidative activities on different substrates (modified corn starch > amylose > amylopectin > corn starch). The results of enzymatic textile desizing application showed that the best combination of amylase (5 g/L), AlLPMO13 (5 mg/L), and H2O2 (3 g/L) made the desizing level and the capillary effects increased by 3 grades and more than 20%, respectively, compared with the results treated by only amylase. CONCLUSION: The Hidden Markov model constructed basing on 34 AA13 family LPMOs was proved to be a valid bioinformatics tool for discovering novel starch-active LPMOs. The novel enzyme AlLPMO13 has strong development potential in the enzymatic textile industry both concerning on economy and on application effect.

PHR-search: a search framework for protein remote homology detection based on the predicted protein hierarchical relationships.

Protein remote homology detection is one of the most fundamental research tool for protein structure and function prediction. Most search methods for protein remote homology detection are evaluated based on the Structural Classification of Proteins-extended (SCOPe) benchmark, but the diverse hierarchical structure relationships between the query protein and candidate proteins are ignored by these methods. In order to further improve the predictive performance for protein remote homology detection, a search framework based on the predicted protein hierarchical relationships (PHR-search) is proposed. In the PHR-search framework, the superfamily level prediction information is obtained by extracting the local and global features of the Hidden Markov Model (HMM) profile through a convolution neural network and it is converted to the fold level and class level prediction information according to the hierarchical relationships of SCOPe. Based on these predicted protein hierarchical relationships, filtering strategy and re-ranking strategy are used to construct the two-level search of PHR-search. Experimental results show that the PHR-search framework achieves the state-of-the-art performance by employing five basic search methods, including HHblits, JackHMMER, PSI-BLAST, DELTA-BLAST and PSI-BLASTexB. Furthermore, the web server of PHR-search is established, which can be accessed at http://bliulab.net/PHR-search.

Construct a variable-length fragment library for de novo protein structure prediction.

Although remarkable achievements, such as AlphaFold2, have been made in end-to-end structure prediction, fragment libraries remain essential for de novo protein structure prediction, which can help explore and understand the protein-folding mechanism. In this work, we developed a variable-length fragment library (VFlib). In VFlib, a master structure database was first constructed from the Protein Data Bank through sequence clustering. The hidden Markov model (HMM) profile of each protein in the master structure database was generated by HHsuite, and the secondary structure of each protein was calculated by DSSP. For the query sequence, the HMM-profile was first constructed. Then, variable-length fragments were retrieved from the master structure database through dynamically variable-length profile-profile comparison. A complete method for chopping the query HMM-profile during this process was proposed to obtain fragments with increased diversity. Finally, secondary structure information was used to further screen the retrieved fragments to generate the final fragment library of specific query sequence. The experimental results obtained with a set of 120 nonredundant proteins show that the global precision and coverage of the fragment library generated by VFlib were 55.04% and 94.95% at the RMSD cutoff of 1.5 Å, respectively. Compared with the benchmark method of NNMake, the global precision of our fragment library had increased by 62.89% with equivalent coverage. Furthermore, the fragments generated by VFlib and NNMake were used to predict structure models through fragment assembly. Controlled experimental results demonstrate that the average TM-score of VFlib was 16.00% higher than that of NNMake.

Comparison of actigraphy-measured and parent-reported sleep in association with weight status among preschool children.

This study compared weekday and weekend actigraphy-measured and parent-reported sleep in relation to weight status among preschool-aged children. Participants were 3-6 years old preschoolers from the cross-sectional DAGIS-study with sleep data for ≥2 weekday and ≥2 weekend nights. Parents-reported sleep onset and wake-up times were gathered alongside 24 h hip-worn actigraphy. An unsupervised Hidden-Markov Model algorithm provided actigraphy-measured night time sleep without the guidance of reported sleep times. Waist-to-height ratio and age-and-sex-specific body mass index characterised weight status. Comparison of methods were assessed with consistency in quintile divisions and Spearman correlations. Associations between sleep and weight status were assessed with adjusted regression models. Participants included 638 children (49% girls) with a mean ± SD age of 4.76 ± 0.89. On weekdays, 98%-99% of actigraphy-measured and parent-reported sleep estimates were classified in the same or adjacent quintile and were strongly correlated (rs = 0.79-0.85, p < 0.001). On weekends, 84%-98% of actigraphy-measured and parent-reported sleep estimates were respectively classified and correlations were moderate to strong (rs = 0.62-0.86, p < 0.001). Compared with actigraphy-measured sleep, parent-reported sleep had consistently earlier onset, later wake-up, and greater duration. Earlier actigraphy-measured weekday sleep onset and midpoint were associated with a higher body mass index (respective ß-estimates: -0.63, p < 0.01 and -0.75, p < 0.01) and waist-to-height ratio (-0.004, p = 0.03 and -0.01, p = 0.02). Though the sleep estimation methods were consistent and correlated, actigraphy measures should be favoured as they are more objective and sensitive to identifying associations between sleep timing and weight status compared with parent reports.

Flexible modeling of headache frequency fluctuations in migraine with hidden Markov models.

OBJECTIVE: To explore hidden Markov models (HMMs) as an approach for defining clinically meaningful headache-frequency-based groups in migraine. BACKGROUND: Monthly headache frequency in patients with migraine is known to vary over time. This variation has not been completely characterized and is not well accounted for in the classification of individuals as having chronic or episodic migraine, a diagnosis with potentially significant impacts on the individual. This study investigated variation in reported headache frequency in a migraine population and proposed a model for classifying individuals by frequency while accounting for natural variation. METHODS: The American Registry for Migraine Research (ARMR) was a longitudinal multisite study of United States adults with migraine. Study participants completed quarterly questionnaires and daily headache diaries. A series of HMMs were fit to monthly headache frequency data calculated from the diary data of ARMR. RESULTS: Changes in monthly headache frequency tended to be small, with 47% of transitions resulting in a change of 0 or 1 day. A substantial portion (24%) of months reflected daily headache with individuals ever reporting daily headache likely to consistently report daily headache. An HMM with four states with mean monthly headache frequency emissions of 3.52 (95% Prediction Interval [PI] 0-8), 10.10 (95% PI 4-17), 20.29 (95% PI 12-28), and constant 28 days/month had the best fit of the models tested. Of sequential month-to-month headache frequency transitions, 12% were across the 15-headache days chronic migraine cutoff. Under the HMM, 38.7% of those transitions involved a change in the HMM state, and the remaining 61.3% of the time, a change in chronic migraine classification was not accompanied by a change in the HMM state. CONCLUSION: A divide between the second and third states of this model aligns most strongly with the current episodic/chronic distinction, although there is a meaningful overlap between the states that supports the need for flexibility. An HMM has appealing properties for classifying individuals according to their headache frequency while accounting for natural variation in frequency. This empirically derived model may provide an informative classification approach that is more stable than the use of a single cutoff value.

Subthalamic stimulation modulates motor network in Parkinson's disease: recover, relieve and remodel.

Aberrant dynamic switches between internal brain states are believed to underlie motor dysfunction in Parkinson's disease. Deep brain stimulation of the subthalamic nucleus is a well-established treatment for the motor symptoms of Parkinson's disease, yet it remains poorly understood how subthalamic stimulation modulates the whole-brain intrinsic motor network state dynamics. To investigate this, we acquired resting-state functional magnetic resonance imaging time-series data from 27 medication-free patients with Parkinson's disease (mean age: 64.8 years, standard deviation: 7.6) who had deep brain stimulation electrodes implanted in the subthalamic nucleus, in both on and off stimulation states. Sixteen matched healthy individuals were included as a control group. We adopted a powerful data-driven modelling approach, known as a hidden Markov model, to disclose the emergence of recurring activation patterns of interacting motor regions (whole-brain intrinsic motor network states) via the blood oxygen level-dependent signal detected in the resting-state functional magnetic resonance imaging time-series data from all participants. The estimated hidden Markov model disclosed the dynamics of distinct whole-brain motor network states, including frequency of occurrence, state duration, fractional coverage and their transition probabilities. Notably, the data-driven decoding of whole-brain intrinsic motor network states revealed that subthalamic stimulation reshaped functional network expression and stabilized state transitions. Moreover, subthalamic stimulation improved motor symptoms by modulating key trajectories of state transition within whole-brain intrinsic motor network states. This modulation mechanism of subthalamic stimulation was manifested in three significant effects: recovery, relieving and remodelling effects. Significantly, recovery effects correlated with improvements in tremor and posture symptoms induced by subthalamic stimulation (P < 0.05). Furthermore, subthalamic stimulation was found to restore a relatively low level of fluctuation of functional connectivity in all motor regions to a level closer to that of healthy participants. Also, changes in the fluctuation of functional connectivity between motor regions were associated with improvements in tremor and gait symptoms (P < 0.05). These findings fill a gap in our knowledge of the role of subthalamic stimulation at the level of neural activity, revealing the regulatory effects of subthalamic stimulation on whole-brain inherent motor network states in Parkinson's disease. Our results provide mechanistic insight and explanation for how subthalamic stimulation modulates motor symptoms in Parkinson's disease.

Brain-heart interaction disruption in major depressive disorder: disturbed rhythm modulation of the cardiac cycle on brain transient theta bursts.

Brain neurons support arousal and cognitive activity in the form of spectral transient bursts and cooperate with the peripheral nervous system to adapt to the surrounding environment. However, the temporal dynamics of brain-heart interactions have not been confirmed, and the mechanism of brain-heart interactions in major depressive disorder (MDD) remains unclear. This study aimed to provide direct evidence for brain-heart synchronization in temporal dynamics and clarify the mechanism of brain-heart interaction disruption in MDD. Eight-minute resting-state (closed eyes) electroencephalograph and electrocardiogram signals were acquired simultaneously. The Jaccard index (JI) was used to measure the temporal synchronization between cortical theta transient bursts and cardiac cycle activity (diastole and systole) in 90 MDD patients and 44 healthy controls (HCs) at rest. The deviation JI was used to reflect the equilibrium of brain activity between diastole and systole. The results showed that the diastole JI was higher than the systole JI in both the HC and MDD groups; compared to HCs, the deviation JI attenuated at F4, F6, FC2, and FC4 in the MDD patients. The eccentric deviation JI was negatively correlated with the despair factor scores of the HAMD, and after 4 weeks of antidepressant treatment, the eccentric deviation JI was positively correlated with the despair factor scores of the HAMD. It was concluded that brain-heart synchronization existed in the theta band in healthy individuals and that disturbed rhythm modulation of the cardiac cycle on brain transient theta bursts at right frontoparietal sites led to brain-heart interaction disruption in MDD.

The frontooccipital interaction mechanism of high-frequency acoustoelectric signal.

Based on acoustoelectric effect, acoustoelectric brain imaging has been proposed, which is a high spatiotemporal resolution neural imaging method. At the focal spot, brain electrical activity is encoded by focused ultrasound, and corresponding high-frequency acoustoelectric signal is generated. Previous studies have revealed that acoustoelectric signal can also be detected in other non-focal brain regions. However, the processing mechanism of acoustoelectric signal between different brain regions remains sparse. Here, with acoustoelectric signal generated in the left primary visual cortex, we investigated the spatial distribution characteristics and temporal propagation characteristics of acoustoelectric signal in the transmission. We observed a strongest transmission strength within the frontal lobe, and the global temporal statistics indicated that the frontal lobe features in acoustoelectric signal transmission. Then, cross-frequency phase-amplitude coupling was used to investigate the coordinated activity in the AE signal band range between frontal and occipital lobes. The results showed that intra-structural cross-frequency coupling and cross-structural coupling co-occurred between these two lobes, and, accordingly, high-frequency brain activity in the frontal lobe was effectively coordinated by distant occipital lobe. This study revealed the frontooccipital long-range interaction mechanism of acoustoelectric signal, which is the foundation of improving the performance of acoustoelectric brain imaging.