Movement decisions driving metapopulation connectivity respond to social resources in a long-lived ungulate, bighorn sheep (Ovis canadensis).

The spatial availability of social resources is speculated to structure animal movement decisions, but the effects of social resources on animal movements are difficult to identify because social resources are rarely measured. Here, we assessed whether varying availability of a key social resource-access to receptive mates-produces predictable changes in movement decisions among bighorn sheep in Nevada, the United States. We compared the probability that males made long-distance 'foray' movements, a critical driver of connectivity, across three ecoregions with varying temporal duration of a socially mediated factor, breeding season. We used a hidden Markov model to identify foray events and then quantified the effects of social covariates on the probability of foray using a discrete choice model. We found that males engaged in forays at higher rates when the breeding season was short, suggesting that males were most responsive to the social resource when its existence was short lived. During the breeding season, males altered their response to social covariates, relative to the non-breeding season, though patterns varied, and age was associated with increased foray probability. Our results suggest that animals respond to the temporal availability of social resources when making the long-distance movements that drive connectivity. This article is part of the theme issue 'The spatial-social interface: a theoretical and empirical integration'.

Evidence of longitudinal differences in spring migration strategies of an Arctic-nesting goose.

During spring, migratory birds are required to optimally balance energetic costs of migration across heterogeneous landscapes and weather conditions to survive and reproduce successfully. Therefore, an individual's migratory performance may influence reproductive outcomes. Given large-scale changes in land use, climate, and potential carry-over effects, understanding how individuals migrate in relation to breeding outcomes is critical to predicting how future scenarios may affect populations. We used GPS tracking devices on 56 Greater White-fronted Geese (Anser albifrons) during four spring migrations to examine whether migration characteristics influenced breeding propensity and breeding outcome. We found a strong longitudinal difference in arrival to the breeding areas (18 days earlier), pre-nesting duration (90.9% longer), and incubation initiation dates (9 days earlier) between western- and eastern-Arctic breeding regions, with contrasting effects on breeding outcomes, but no migration characteristic strongly influenced breeding outcome. We found that breeding region influenced whether an individual likely pursued a capital or income breeding strategy. Where individuals fell along the capital-income breeding continuum was influenced by longitude, revealing geographic effects of life-history strategy among conspecifics. Factors that govern breeding outcomes likely occur primarily upon arrival to breeding areas or are related to individual quality and previous breeding outcome, and may not be directly tied to migratory decision-making across broad scales.

Comparison between EEG and MEG of static and dynamic resting-state networks.

The characterisation of resting-state networks (RSNs) using neuroimaging techniques has significantly contributed to our understanding of the organisation of brain activity. Prior work has demonstrated the electrophysiological basis of RSNs and their dynamic nature, revealing transient activations of brain networks with millisecond timescales. While previous research has confirmed the comparability of RSNs identified by electroencephalography (EEG) to those identified by magnetoencephalography (MEG) and functional magnetic resonance imaging (fMRI), most studies have utilised static analysis techniques, ignoring the dynamic nature of brain activity. Often, these studies use high-density EEG systems, which limit their applicability in clinical settings. Addressing these gaps, our research studies RSNs using medium-density EEG systems (61 sensors), comparing both static and dynamic brain network features to those obtained from a high-density MEG system (306 sensors). We assess the qualitative and quantitative comparability of EEG-derived RSNs to those from MEG, including their ability to capture age-related effects, and explore the reproducibility of dynamic RSNs within and across the modalities. Our findings suggest that both MEG and EEG offer comparable static and dynamic network descriptions, albeit with MEG offering some increased sensitivity and reproducibility. Such RSNs and their comparability across the two modalities remained consistent qualitatively but not quantitatively when the data were reconstructed without subject-specific structural MRI images.

A topical VAEGAN-IHMM approach for automatic story segmentation.

Feature representations with rich topic information can greatly improve the performance of story segmentation tasks. VAEGAN offers distinct advantages in feature learning by combining variational autoencoder (VAE) and generative adversarial network (GAN), which not only captures intricate data representations through VAE's probabilistic encoding and decoding mechanism but also enhances feature diversity and quality via GAN's adversarial training. To better learn topical domain representation, we used a topical classifier to supervise the training process of VAEGAN. Based on the learned feature, a segmentor splits the document into shorter ones with different topics. Hidden Markov model (HMM) is a popular approach for story segmentation, in which stories are viewed as instances of topics (hidden states). The number of states has to be set manually but it is often unknown in real scenarios. To solve this problem, we proposed an infinite HMM (IHMM) approach which utilized an HDP prior on transition matrices over countably infinite state spaces to automatically infer the state's number from the data. Given a running text, a Blocked Gibbis sampler labeled the states with topic classes. The position where the topic changes was a story boundary. Experimental results on the TDT2 corpus demonstrated that the proposed topical VAEGAN-IHMM approach was significantly better than the traditional HMM method in story segmentation tasks and achieved state-of-the-art performance.

Dynamic functional connectivity and gene expression correlates in temporal lobe epilepsy: insights from hidden markov models.

BACKGROUD: Temporal lobe epilepsy (TLE) is associated with abnormal dynamic functional connectivity patterns, but the dynamic changes in brain activity at each time point remain unclear, as does the potential molecular mechanisms associated with the dynamic temporal characteristics of TLE. METHODS: Resting-state functional magnetic resonance imaging (rs-fMRI) was acquired for 84 TLE patients and 35 healthy controls (HCs). The data was then used to conduct HMM analysis on rs-fMRI data from TLE patients and an HC group in order to explore the intricate temporal dynamics of brain activity in TLE patients with cognitive impairment (TLE-CI). Additionally, we aim to examine the gene expression profiles associated with the dynamic modular characteristics in TLE patients using the Allen Human Brain Atlas (AHBA) database. RESULTS: Five HMM states were identified in this study. Compared with HCs, TLE and TLE-CI patients exhibited distinct changes in dynamics, including fractional occupancy, lifetimes, mean dwell time and switch rate. Furthermore, transition probability across HMM states were significantly different between TLE and TLE-CI patients (p < 0.05). The temporal reconfiguration of states in TLE and TLE-CI patients was associated with several brain networks (including the high-order default mode network (DMN), subcortical network (SCN), and cerebellum network (CN). Furthermore, a total of 1580 genes were revealed to be significantly associated with dynamic brain states of TLE, mainly enriched in neuronal signaling and synaptic function. CONCLUSIONS: This study provides new insights into characterizing dynamic neural activity in TLE. The brain network dynamics defined by HMM analysis may deepen our understanding of the neurobiological underpinnings of TLE and TLE-CI, indicating a linkage between neural configuration and gene expression in TLE.

Absolute direction in organelle movement.

In movement analysis, correlated random walk (CRW) models often use so-called turning angles, which are measured relative to the previous movement direction. To segregate between different movement modes, hidden Markov models (HMMs) describe movements as piecewise stationary CRWs in which the distributions of turning angles and step sizes depend on the underlying state. This typically allows for the segregation of movement modes that show different movement speeds. We show that in some cases, it may be interesting to investigate absolute angles, that is, biased random walks (BRWs) instead of turning angles. In particular, while discrimination between states in the turning angle setting can only rely on movement speed, models with absolute angles can be used to discriminate between sections of different movement directions. A preprocessing algorithm is provided that enables the analysis of absolute angles in the existing R package moveHMM. In a data set of movements of cell organelles, models using not the turning angle but the absolute angle could capture interesting additional properties. Goodness-of-fit was increased for HMMs with absolute angles, and HMMs with absolute angles tended to choose a higher number of states, suggesting the existence and relevance of prominent directional changes in the present data set. These results suggest that models with absolute angles can provide important information in the analysis of movement patterns if the existence and frequency of directional changes is of biological importance.

CW-PRED: Prediction of C-terminal surface anchoring sorting signals in bacteria and Archaea.

Sorting signals are crucial for the anchoring of proteins to the cell surface in archaea and bacteria. These proteins often feature distinct motifs at their C-terminus, cleaved by sortase or sortase-like enzymes. Gram-positive bacteria exhibit the LPXTGX consensus motif, cleaved by sortases, while Gram-negative bacteria employ exosortases recognizing motifs like PEP. Archaea utilize exosortase homologs known as archaeosortases for signal anchoring. Traditionally identification of such C-terminal sorting signals was performed with profile Hidden Markov Models (pHMMs). The Cell- Wall PREDiction (CW-PRED) method introduced for the first time a custom-made class HMM for proteins in Gram-positive bacteria that contain a cell wall sorting signal which begins with an LPXTG motif, followed by a hydrophobic domain and a tail of positively charged residues. Here we present a new and updated version of CW-PRED for predicting C-terminal sorting signals in Archaea, Gram-positive, and Gram-negative bacteria. We used a large training set and several model enhancements that improve motif identification in order to achieve better discrimination between C-terminal signals and other proteins. Cross-validation demonstrates CW-PRED's superiority in sensitivity and specificity compared to other methods. Application of the method in reference proteomes reveals a large number of potential surface proteins not previously identified. The method is available for academic use at http://195.251.108.230/apps.compgen.org/CW-PRED/ and as standalone software.

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.

Artificial intelligence tools for the identification of antibiotic resistance genes.

The fight against bacterial antibiotic resistance must be given critical attention to avert the current and emerging crisis of treating bacterial infections due to the inefficacy of clinically relevant antibiotics. Intrinsic genetic mutations and transferrable antibiotic resistance genes (ARGs) are at the core of the development of antibiotic resistance. However, traditional alignment methods for detecting ARGs have limitations. Artificial intelligence (AI) methods and approaches can potentially augment the detection of ARGs and identify antibiotic targets and antagonistic bactericidal and bacteriostatic molecules that are or can be developed as antibiotics. This review delves into the literature regarding the various AI methods and approaches for identifying and annotating ARGs, highlighting their potential and limitations. Specifically, we discuss methods for (1) direct identification and classification of ARGs from genome DNA sequences, (2) direct identification and classification from plasmid sequences, and (3) identification of putative ARGs from feature selection.

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.

Proteome-wide neuropeptide identification using NeuroPeptide-HMMer (NP-HMMer).

Neuropeptides are essential neuronal signaling molecules that orchestrate animal behavior and physiology via actions within the nervous system and on peripheral tissues. Due to the small size of biologically active mature peptides, their identification on a proteome-wide scale poses a significant challenge using existing bioinformatics tools like BLAST. To address this, we have developed NeuroPeptide-HMMer (NP-HMMer), a hidden Markov model (HMM)-based tool to facilitate neuropeptide discovery, especially in underexplored invertebrates. NP-HMMer utilizes manually curated HMMs for 46 neuropeptide families, enabling rapid and accurate identification of neuropeptides. Validation of NP-HMMer on Drosophila melanogaster, Daphnia pulex, Tribolium castaneum and Tenebrio molitor demonstrated its effectiveness in identifying known neuropeptides across diverse arthropods. Additionally, we showcase the utility of NP-HMMer by discovering novel neuropeptides in Priapulida and Rotifera, identifying 22 and 19 new peptides, respectively. This tool represents a significant advancement in neuropeptide research, offering a robust method for annotating neuropeptides across diverse proteomes and providing insights into the evolutionary conservation of neuropeptide signaling pathways.

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.

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.

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.

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.

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.

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.

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%.

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.

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.