Characterization of domain formation in complex membranes.

A widely known property of lipid membranes is their tendency to undergo a separation into disordered (Ld) and ordered (Lo) domains. This impacts the local structure of the membrane relevant for the physical (e.g., enhanced electroporation) and biological (e.g., protein sorting) significance of these regions. The increase in computing power, advancements in simulation software, and more detailed information about the composition of biological membranes shifts the study of these domains into the focus of classical molecular dynamics simulations. In this chapter, we present a versatile yet robust analysis pipeline that can be easily implemented and adapted for a wide range of lipid compositions. It employs Gaussian-based Hidden Markov Models to predict the hidden order states of individual lipids by describing their structure through the area per lipid and the average SCC order parameters per acyl chain. Regions of the membrane with a high correlation between ordered lipids are identified by employing the Getis-Ord local spatial autocorrelation statistic on a Voronoi tessellation of the lipids. As an example, the approach is applied to two distinct systems at a coarse-grained resolution, demonstrating either a strong tendency towards phase separation (1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), 1,2-dilinoleoyl-sn-glycero-3-phosphocholine (DIPC), cholesterol) or a weak tendency toward phase separation (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), 1-palmitoyl-2-docosahexaenoyl-sn-glycero-3-phosphocholine (PUPC), cholesterol). Explanations of the steps are complemented by coding examples written in Python, providing both a comprehensive understanding and practical guidance for a seamless integration of the workflow into individual projects.

Fuzzy random sensitivity analysis for the overall structure reliability of reinforced concrete freezing wellbores in deep alluvium based on hidden Markov model.

To address the shortcomings of traditional reliability theory in characterizing the stability of deep underground structures, the advanced first order second moment of reliability was improved to obtain fuzzy random reliability, which is more consistent with the working conditions. The traditional sensitivity analysis model was optimized using fuzzy random optimization, and an analytical calculation model of the mean and standard deviation of the fuzzy random reliability sensitivity was established. A big data hidden Markov model and expectation-maximization algorithm were used to improve the digital characteristics of fuzzy random variables. The fuzzy random sensitivity optimization model was used to confirm the effect of concrete compressive strength, thick-diameter ratio, reinforcement ratio, uncertainty coefficient of calculation model, and soil depth on the overall structural reliability of a reinforced concrete double-layer wellbore in deep alluvial soil. Through numerical calculations, these characteristics were observed to be the main influencing factors. Furthermore, while the soil depth was negatively correlated, the other influencing factors were all positively correlated with the overall reliability. This study provides an effective reference for the safe construction of deep underground structures in the future.

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.

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.

SignalP: The Evolution of a Web Server.

SignalP ( https://services.healthtech.dtu.dk/services/SignalP-6.0/ ) is a very popular prediction method for signal peptides, the intrinsic signals that make proteins secretory. The SignalP web server has existed since 1995 and is now in its sixth major version. In this historical account, we (three authors who have taken part in the entire journey plus the first author of the latest version) describe the differences between the versions and discuss the various decisions taken along the way.

Abnormal dynamic functional connectivity and topological properties of cerebellar network in male obstructive sleep apnea.

PURPOSE: To investigate dynamic functional connectivity (dFC) within the cerebellar-whole brain network and dynamic topological properties of the cerebellar network in obstructive sleep apnea (OSA) patients. METHODS: Sixty male patients and 60 male healthy controls were included. The sliding window method examined the fluctuations in cerebellum-whole brain dFC and connection strength in OSA. Furthermore, graph theory metrics evaluated the dynamic topological properties of the cerebellar network. Additionally, hidden Markov modeling validated the robustness of the dFC. The correlations between the abovementioned measures and clinical assessments were assessed. RESULTS: Two dynamic network states were characterized. State 2 exhibited a heightened frequency, longer fractional occupancy, and greater mean dwell time in OSA. The cerebellar networks and cerebrocerebellar dFC alterations were mainly located in the default mode network, frontoparietal network, somatomotor network, right cerebellar CrusI/II, and other networks. Global properties indicated aberrant cerebellar topology in OSA. Dynamic properties were correlated with clinical indicators primarily on emotion, cognition, and sleep. CONCLUSION: Abnormal dFC in male OSA may indicate an imbalance between the integration and segregation of brain networks, concurrent with global topological alterations. Abnormal default mode network interactions with high-order and low-level cognitive networks, disrupting their coordination, may impair the regulation of cognitive, emotional, and sleep functions in OSA.

Place Recognition through Multiple LiDAR Scans Based on the Hidden Markov Model.

Autonomous driving systems for unmanned ground vehicles (UGV) operating in enclosed environments strongly rely on LiDAR localization with a prior map. Precise initial pose estimation is critical during system startup or when tracking is lost, ensuring safe UGV operation. Existing LiDAR-based place recognition methods often suffer from reduced accuracy due to only matching descriptors from individual LiDAR keyframes. This paper proposes a multi-frame descriptor-matching approach based on the hidden Markov model (HMM) to address this issue. This method enhances the place recognition accuracy and robustness by leveraging information from multiple frames. Experimental results from the KITTI dataset demonstrate that the proposed method significantly enhances the place recognition performance compared with the scan context-based single-frame descriptor-matching approach, with an average performance improvement of 5.8% and with a maximum improvement of 15.3%.

Predicting nodal metastasis progression of oral tongue cancer using a hidden Markov model in MRI.

Objectives: The presence of occult nodal metastases in patients with oral tongue squamous cell carcinomas (OTSCCs) has implications for treatment. More than 30% of patients will have occult nodal metastases, yet a considerable number of patients undergo unnecessary invasive neck dissection to confirm nodal status. In this work, we propose a probabilistic model for lymphatic metastatic spread that can quantify the risk of microscopic involvement at the lymph node level (LNL) given the location of macroscopic metastases and the tumor stage using the MRI method. Materials and methods: A total of 108 patients of OTSCCs were included in the study. A hidden Markov model (HMM) was used to compute the probabilities of transitions between states over time based on MRI. Learning of the transition probabilities was performed via Markov chain Monte Carlo sampling and was based on a dataset of OTSCC patients for whom involvement of individual LNLs was reported. Results: Our model found that the most common involvement was that of level I and level II, corresponding to a high probability of 𝑝b1 = 0.39 ± 0.05, 𝑝b2 = 0.53 ± 0.09; lymph node level I had metastasis, and the probability of metastasis in lymph node II was high (93.79%); lymph node level II had metastasis, and the probability of metastasis in lymph node III was small (7.88%). Lymph nodes progress faster in the early stage and slower in the late stage. Conclusion: An HMM can produce an algorithm that is able to predict nodal metastasis evolution in patients with OTSCCs by analyzing the macroscopic metastases observed in the upstream levels, and tumor category.

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.

Understanding Human Cognition Through Computational Modeling.

One important goal of cognitive science is to understand the mind in terms of its representational and computational capacities, where computational modeling plays an essential role in providing theoretical explanations and predictions of human behavior and mental phenomena. In my research, I have been using computational modeling, together with behavioral experiments and cognitive neuroscience methods, to investigate the information processing mechanisms underlying learning and visual cognition in terms of perceptual representation and attention strategy. In perceptual representation, I have used neural network models to understand how the split architecture in the human visual system influences visual cognition, and to examine perceptual representation development as the results of expertise. In attention strategy, I have developed the Eye Movement analysis with Hidden Markov Models method for quantifying eye movement pattern and consistency using both spatial and temporal information, which has led to novel findings across disciplines not discoverable using traditional methods. By integrating it with deep neural networks (DNN), I have developed DNN+HMM to account for eye movement strategy learning in human visual cognition. The understanding of the human mind through computational modeling also facilitates research on artificial intelligence's (AI) comparability with human cognition, which can in turn help explainable AI systems infer humans' belief on AI's operations and provide human-centered explanations to enhance human-AI interaction and mutual understanding. Together, these demonstrate the essential role of computational modeling methods in providing theoretical accounts of the human mind as well as its interaction with its environment and AI systems.

A Bayesian hierarchical hidden Markov model for clustering and gene selection: Application to kidney cancer gene expression data.

We introduce a Bayesian approach for biclustering that accounts for the prior functional dependence between genes using hidden Markov models (HMMs). We utilize biological knowledge gathered from gene ontologies and the hidden Markov structure to capture the potential coexpression of neighboring genes. Our interpretable model-based clustering characterized each cluster of samples by three groups of features: overexpressed, underexpressed, and irrelevant features. The proposed methods have been implemented in an R package and are used to analyze both the simulated data and The Cancer Genome Atlas kidney cancer data.

Spatial-Temporal Characteristics of Brain Activity in Autism Spectrum Disorder Based on Hidden Markov Model and Dynamic Graph Theory: A Resting-State fMRI Study.

Autism spectrum disorder (ASD) is a common neurodevelopmental disorder. Functional magnetic resonance imaging (fMRI) can be used to measure the temporal correlation of blood-oxygen-level-dependent (BOLD) signals in the brain to assess the brain's intrinsic connectivity and capture dynamic changes in the brain. In this study, the hidden Markov model (HMM) and dynamic graph (DG) theory are used to study the spatial-temporal characteristics and dynamics of brain networks based on dynamic functional connectivity (DFC). By using HMM, we identified three typical brain states for ASD and healthy control (HC). Furthermore, we explored the correlation between HMM time-varying properties and clinical autism scale scores. Differences in brain topological characteristics and dynamics between ASD and HC were compared by DG analysis. The experimental results indicate that ASD is more inclined to enter a strongly connected HMM brain state, leading to the isolation of brain networks and alterations in the topological characteristics of brain networks, such as default mode network (DMN), ventral attention network (VAN), and visual network (VN). This work suggests that using different data-driven methods based on DFC to study brain network dynamics would have better information complementarity, which can provide a new direction for the extraction of neuro-biomarkers in the early diagnosis of ASD.

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.

Auditory cues modulate the short timescale dynamics of STN activity during stepping in Parkinson's disease.

BACKGROUND: Gait impairment has a major impact on quality of life in patients with Parkinson's disease (PD). It is believed that basal ganglia oscillatory activity at ß frequencies (15-30 Hz) may contribute to gait impairment, but the precise dynamics of this oscillatory activity during gait remain unclear. Additionally, auditory cues are known to lead to improvements in gait kinematics in PD. If the neurophysiological mechanisms of this cueing effect were better understood they could be leveraged to treat gait impairments using adaptive Deep Brain Stimulation (aDBS) technologies. OBJECTIVE: We aimed to characterize the dynamics of subthalamic nucleus (STN) oscillatory activity during stepping movements in PD and to establish the neurophysiological mechanisms by which auditory cues modulate gait. METHODS: We studied STN local field potentials (LFPs) in eight PD patients while they performed stepping movements. Hidden Markov Models (HMMs) were used to discover transient states of spectral activity that occurred during stepping with and without auditory cues. RESULTS: The occurrence of low and high ß bursts was suppressed during and after auditory cues. This manifested as a decrease in their fractional occupancy and state lifetimes. Interestingly, α transients showed the opposite effect, with fractional occupancy and state lifetimes increasing during and after auditory cues. CONCLUSIONS: We show that STN oscillatory activity in the α and ß frequency bands are differentially modulated by gait-promoting oscillatory cues. These findings suggest that the enhancement of α rhythms may be an approach for ameliorating gait impairments in PD.

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.

Combining acoustic telemetry with archival tagging to investigate the spatial dynamic of the understudied pollack, Pollachius pollachius.

Combining fish tracking methods is a promising way of leveraging the strengths of each approach while mitigating their individual weaknesses. Acoustic telemetry provides presence information as the fish move within receiver range, eliminating the need for tag recovery. Archival tags, on the other hand, record environmental variables on tag retrieval, enabling continuous path reconstruction of a fish beyond coastal regions. This study capitalizes on the combination of both methods for geolocating pollack, Pollachius pollachius, an understudied species of the northeast Atlantic, where declining stocks are raising concern. Essential knowledge of population structure and connectivity between essential habitats is critically lacking and could help inform stock assessment and management. The aims of the study were (1) to evaluate the feasibility of double-tagging pollack, known for being challenging to tag, and (2) to track seasonal movements across the Channel to gain first insights into pollack spatial ecology. In 2022, an extensive network of acoustic receivers was been deployed in the Channel along the French, English, and Belgian coasts as part of the Fish Intel project. We tagged 83 pollack with acoustic transmitters, among which 48 were double-tagged with data storage tags. Post-tagging survival assessment, conducted on a subset of 35 individuals, revealed a successful procedure with a 97% short-term survival rate. By October 2023, the acoustic telemetry network detected 30 out of 83 pollack at least once, with no large-scale movements observed across the Channel. Presence in the network fluctuates seasonally, peaking in summer, particularly among immature fish. Integrating acoustic detections with temperature and depth time series in a geolocation model enabled trajectory reconstruction of 10 recaptured pollack, seven of which were detected by the network. This combined tracking approach revealed coastal movements along the coast of Brittany in France, highlighting the ecological significance of the Iroise Sea for pollack throughout the year, particularly in summer. The geolocation model also suggested movements towards the entrance of the western Channel. This study highlights the complementarity of acoustic telemetry and archival tagging in reconstructing fish movements in their natural environment. As data accumulate, these innovative tracking methods promise to continually unveil new insights into the spatial ecology of the understudied pollack, which is essential for the species' management.

Auxiliary two-filter particle smoothing for one generalized hidden Markov model.

This paper develops two-filter particle smoothing (TFPS) algorithms for the nonlinear fixed-interval smoothing problem of one generalized hidden Markov model (GHMM), where the current observation depends not only on the current state, but also on one-step previous state. Firstly, by Bayesian approach, the two-filter smoothing (TFS) formula for GHMM is established to calculate smoothing densities. In this TFS formula, the backward information prediction density is generally not a density of the state. This results in a difficulty that the normal sequential Monte Carlo (SMC) sampling technique cannot be directly applied to design corresponding TFPS algorithms based on the TFS formula. To solve this difficulty, a generalized TFS formula for GHMM is then proposed by introducing a sequence of artificial densities. By combining this generalized TFS formula, SMC, and the auxiliary variable sampling technique, a basic auxiliary TFPS (ATFPS) algorithm with quadratic computational complexity is proposed, and a simplified ATFPS algorithm with linear computational complexity is further developed. Finally, the effectiveness and superiority of the two proposed ATFPS algorithms for GHMM are verified via simulation examples and real experimental data.

Noise reduction of electron holography observations for a thin-foiled Nd-Fe-B specimen using the wavelet hidden Markov model.

In this study, we investigate the effectiveness of noise reduction in electron holography, based on the wavelet hidden Markov model (WHMM), which allows the reasonable separation of weak signals from noise. Electron holography observations from a Nd2Fe14B thin foil showed that the noise reduction method suppressed artificial phase discontinuities generated by phase retrieval. From the peak signal-to-noise ratio, it was seen that the impact of denoising was significant for observations with a narrow spacing of interference fringes, which is a key parameter for the spatial resolution of electron holography. These results provide essential information for improving the precision of electron holography studies.

Use of DNA forceps to measure receptor-ligand dissociation equilibrium constants in a single-molecule competition assay.

The ability of biophysicists to decipher the behavior of individual biomolecules has steadily improved over the past thirty years. However, it still remains unclear how an ensemble of data acquired at the single-molecule level compares with the data acquired on an ensemble of the same molecules. We here propose an assay to tackle this question in the context of dissociation equilibrium constant measurements. A sensor is built by engrafting a receptor and a ligand onto a flexible dsDNA scaffold and mounting this assembly on magnetic tweezers. This way, looking at the position of the magnetic bead enables one to determine in real-time if the two molecular partners are associated or not. Next, to quantify the affinity of the scrutinized single-receptor for a given competitor, various amounts of the latter molecule are introduced in solution and the equilibrium response of the sensor is monitored throughout the titration protocol. Proofs of concept are established for the binding of three rapamycin analogs to the FKBP12 cis-trans prolyl isomerase. For each of these drugs the mean affinity constant obtained on a ten of individual receptors agrees with the one previously determined in a bulk assay. Furthermore, experimental contingencies are sufficient to explain the dispersion observed over the single-molecule values.

A Hybrid Approach for Energy Consumption and Improvement in Sensor Network Lifespan in Wireless Sensor Networks.

In this paper, we propose an improved clustering algorithm for wireless sensor networks (WSNs) that aims to increase network lifetime and efficiency. We introduce an enhanced fuzzy spider monkey optimization technique and a hidden Markov model-based clustering algorithm for selecting cluster heads. Our approach considers factors such as network cluster head energy, cluster head density, and cluster head position. We also enhance the energy-efficient routing strategy for connecting cluster heads to the base station. Additionally, we introduce a polling control method to improve network performance while maintaining energy efficiency during steady transmission periods. Simulation results demonstrate a 1.2% improvement in network performance using our proposed model.