Categories
Nevin Manimala Statistics

Interictal Connectivity Revealed by Granger Analysis of Stereoelectroencephalography: Association With Ictal Onset Zone, Resection, and Outcome

Neurosurgery. 2022 Aug 22. doi: 10.1227/neu.0000000000002079. Online ahead of print.

ABSTRACT

BACKGROUND: Stereoelectroencephalography (sEEG) facilitates electrical sampling and evaluation of complex deep-seated, dispersed, and multifocal locations. Granger causality (GC), previously used to study seizure networks using interictal data from subdural grids, may help identify the seizure-onset zone from interictal sEEG recordings.

OBJECTIVE: To examine whether statistical analysis of interictal sEEG helps identify surgical target sites and whether surgical resection of highly ranked nodes correspond to favorable outcomes.

METHODS: Ten minutes of extraoperative recordings from sequential patients who underwent sEEG evaluation were analyzed (n = 20). GC maps were compared with clinically defined surgical targets using rank order statistics. Outcomes of patients with focal resection/ablation with median follow-up of 3.6 years were classified as favorable (Engel 1, 2) or poor (Engel 3, 4) to assess their relationship with the removal of highly ranked nodes using the Wilcoxon rank-sum test.

RESULTS: In 12 of 20 cases, the rankings of contacts (based on the sum of outward connection weights) mapped to the seizure-onset zone showed higher causal node connectivity than predicted by chance (P ≤ .02). A very low aggregate probability (P < 10-18, n = 20) suggests that causal node connectivity predicts seizure networks. In 8 of 16 with outcome data, causal connectivity in the resection was significantly greater than in the remaining contacts (P ≤ .05). We found a significant association between favorable outcome and the presence of highly ranked nodes in the resection (P < .05).

CONCLUSION: Granger analysis can identify seizure foci from interictal sEEG and correlates highly ranked nodes with favorable outcome, potentially informing surgical decision-making without reliance on ictal recordings.

PMID:36084171 | DOI:10.1227/neu.0000000000002079

Categories
Nevin Manimala Statistics

ceRNAR: An R package for identification and analysis of ceRNA-miRNA triplets

PLoS Comput Biol. 2022 Sep 9;18(9):e1010497. doi: 10.1371/journal.pcbi.1010497. Online ahead of print.

ABSTRACT

Competitive endogenous RNA (ceRNA) represents a novel mechanism of gene regulation that controls several biological and pathological processes. Recently, an increasing number of in silico methods have been developed to accelerate the identification of such regulatory events. However, there is still a need for a tool supporting the hypothesis that ceRNA regulatory events only occur at specific miRNA expression levels. To this end, we present an R package, ceRNAR, which allows identification and analysis of ceRNA-miRNA triplets via integration of miRNA and RNA expression data. The ceRNAR package integrates three main steps: (i) identification of ceRNA pairs based on a rank-based correlation between pairs that considers the impact of miRNA and a running sum correlation statistic, (ii) sample clustering based on gene-gene correlation by circular binary segmentation, and (iii) peak merging to identify the most relevant sample patterns. In addition, ceRNAR also provides downstream analyses of identified ceRNA-miRNA triplets, including network analysis, functional annotation, survival analysis, external validation, and integration of different tools. The performance of our proposed approach was validated through simulation studies of different scenarios. Compared with several published tools, ceRNAR was able to identify true ceRNA triplets with high sensitivity, low false-positive rates, and acceptable running time. In real data applications, the ceRNAs common to two lung cancer datasets were identified in both datasets. The bridging miRNA for one of these, the ceRNA for MAP4K3, was identified by ceRNAR as hsa-let-7c-5p. Since similar cancer subtypes do share some biological patterns, these results demonstrated that our proposed algorithm was able to identify potential ceRNA targets in real patients. In summary, ceRNAR offers a novel algorithm and a comprehensive pipeline to identify and analyze ceRNA regulation. The package is implemented in R and is available on GitHub (https://github.com/ywhsiao/ceRNAR).

PMID:36084156 | DOI:10.1371/journal.pcbi.1010497

Categories
Nevin Manimala Statistics

Conditional GWAS of non-CG transposon methylation in Arabidopsis thaliana reveals major polymorphisms in five genes

PLoS Genet. 2022 Sep 9;18(9):e1010345. doi: 10.1371/journal.pgen.1010345. Online ahead of print.

ABSTRACT

Genome-wide association studies (GWAS) have revealed that the striking natural variation for DNA CHH-methylation (mCHH; H is A, T, or C) of transposons has oligogenic architecture involving major alleles at a handful of known methylation regulators. Here we use a conditional GWAS approach to show that CHG-methylation (mCHG) has a similar genetic architecture-once mCHH is statistically controlled for. We identify five key trans-regulators that appear to modulate mCHG levels, and show that they interact with a previously identified modifier of mCHH in regulating natural transposon mobilization.

PMID:36084135 | DOI:10.1371/journal.pgen.1010345

Categories
Nevin Manimala Statistics

Methodology to estimate natural- and vaccine-induced antibodies to SARS-CoV-2 in a large geographic region

PLoS One. 2022 Sep 9;17(9):e0273694. doi: 10.1371/journal.pone.0273694. eCollection 2022.

ABSTRACT

Accurate estimates of natural and/or vaccine-induced antibodies to SARS-CoV-2 are difficult to obtain. Although model-based estimates of seroprevalence have been proposed, they require inputting unknown parameters including viral reproduction number, longevity of immune response, and other dynamic factors. In contrast to a model-based approach, the current study presents a data-driven detailed statistical procedure for estimating total seroprevalence (defined as antibodies from natural infection or from full vaccination) in a region using prospectively collected serological data and state-level vaccination data. Specifically, we conducted a longitudinal statewide serological survey with 88,605 participants 5 years or older with 3 prospective blood draws beginning September 30, 2020. Along with state vaccination data, as of October 31, 2021, the estimated percentage of those 5 years or older with naturally occurring antibodies to SARS-CoV-2 in Texas is 35.0% (95% CI = (33.1%, 36.9%)). This is 3× higher than, state-confirmed COVID-19 cases (11.83%) for all ages. The percentage with naturally occurring or vaccine-induced antibodies (total seroprevalence) is 77.42%. This methodology is integral to pandemic preparedness as accurate estimates of seroprevalence can inform policy-making decisions relevant to SARS-CoV-2.

PMID:36084125 | DOI:10.1371/journal.pone.0273694

Categories
Nevin Manimala Statistics

Prediction model for an early revision for dislocation after primary total hip arthroplasty

PLoS One. 2022 Sep 9;17(9):e0274384. doi: 10.1371/journal.pone.0274384. eCollection 2022.

ABSTRACT

Dislocation is one of the most common complications after primary total hip arthroplasty (THA). Several patient-related risk factors for dislocation have been reported in the previous literature, but only few prediction models for dislocation have been made. Our aim was to build a prediction model for an early (within the first 2 years) revision for dislocation after primary THA using two different statistical methods. The study data constituted of 37 pre- or perioperative variables and postoperative follow-up data of 16 454 primary THAs performed at our institution in 2008-2021. Model I was a traditional logistic regression model and Model II was based on the elastic net method that utilizes machine learning. The models’ overall performance was measured using the pseudo R2 values. The discrimination of the models was measured using C-index in Model I and Area Under the Curve (AUC) in Model II. Calibration curves were made for both models. At 2 years postoperatively, 95 hips (0.6% prevalence) had been revised for dislocation. The pseudo R2 values were 0.04 in Model I and 0.02 in Model II indicating low predictive capability in both models. The C-index in Model I was 0.67 and the AUC in Model II was 0.73 indicating modest discrimination. The prediction of an early revision for dislocation after primary THA is difficult even in a large cohort of patients with detailed data available because of the reasonably low prevalence and multifactorial nature of dislocation. Therefore, the risk of dislocation should be kept in mind in every primary THA, whether the patient has predisposing factors for dislocation or not. Further, when conducting a prediction model, sophisticated methods that utilize machine learning may not necessarily offer significant advantage over traditional statistical methods in clinical setup.

PMID:36084121 | DOI:10.1371/journal.pone.0274384

Categories
Nevin Manimala Statistics

Geolocation of multiple sociolinguistic markers in Buenos Aires

PLoS One. 2022 Sep 9;17(9):e0274114. doi: 10.1371/journal.pone.0274114. eCollection 2022.

ABSTRACT

Analysis of language geography is increasingly being used for studying spatial patterns of social dynamics. This trend is fueled by social media platforms such as Twitter which provide access to large amounts of natural language data combined with geolocation and user metadata enabling reconstruction of detailed spatial patterns of language use. Most studies are performed on large spatial scales associated with countries and regions, where language dynamics are often dominated by the effects of geographic and administrative borders. Extending to smaller, urban scales, however, allows visualization of spatial patterns of language use determined by social dynamics within the city, providing valuable information for a range of social topics from demographic studies to urban planning. So far, few studies have been made in this domain, due, in part, to the challenges in developing algorithms that accurately classify linguistic features. Here we extend urban-scale geographical analysis of language use beyond lexical meaning to include other sociolinguistic markers that identify language style, dialect and social groups. Some features, which have not been explored with social-media data on the urban scale, can be used to target a range of social phenomena. Our study focuses on Twitter use in Buenos Aires and our approach classifies tweets based on contrasting sets of tokens manually selected to target precise linguistic features. We perform statistical analyses of eleven categories of language use to quantify the presence of spatial patterns and the extent to which they are socially driven. We then perform the first comparative analysis assessing how the patterns and strength of social drivers vary with category. Finally, we derive plausible explanations for the patterns by comparing them with independently generated maps of geosocial context. Identifying these connections is a key aspect of the social-dynamics analysis which has so far received insufficient attention.

PMID:36084118 | DOI:10.1371/journal.pone.0274114

Categories
Nevin Manimala Statistics

Fibroblast growth factor 23 and kidney function in patients with type 1 diabetes

PLoS One. 2022 Sep 9;17(9):e0274182. doi: 10.1371/journal.pone.0274182. eCollection 2022.

ABSTRACT

Diabetic kidney disease (DKD) is a key determinant of morbidity and mortality in patients with type 1 diabetes (T1D). Identifying factors associated with early glomerular filtration rate (GFR) decline in T1D is important in prevention or early intervention for DKD. This study investigated whether phosphate metabolism, including fibroblast growth factor 23 (FGF23) is associated with the kidney function of patients with T1D. We randomly recruited 118 patients with T1D with a normal or mildly impaired kidney function [chronic kidney disease (CKD) stages of G1/G2, A1/A2], and measured their serum FGF23 levels. Serum FGF23 was significantly negatively associated with the estimated GFR (eGFR) (r = -0.292, P = 0.0016), but not urinary albumin creatinine ratio (UACR), and positively associated with serum phosphate (Pi; r = 0.273, P = 0.0027). Serum FGF23 increased with decreasing eGFR quartiles (P for linear trend = 0.0371), while FGF23 was modestly higher in the higher quartiles of UACR (not statistically significant). The multiple linear regression analysis also showed a significant inverse association between FGF23 and eGFR (Model 1: β = -0.149, P = 0.0429; Model 2: β = -0.141, P = 0.0370). The association remained significant after adjustment for Pi. We identified that FGF23 was inversely associated with the eGFR in T1D patients with a normal or mildly impaired kidney function.

PMID:36084108 | DOI:10.1371/journal.pone.0274182

Categories
Nevin Manimala Statistics

Colloidal Stochastic Resonance in Confined Geometries

Phys Rev Lett. 2022 Aug 26;129(9):098001. doi: 10.1103/PhysRevLett.129.098001.

ABSTRACT

We investigate the dynamical properties of a colloidal particle in a double cavity. Without external driving, the particle hops between two free-energy minima with transition mean time depending on the system’s entropic and energetic barriers. We then drive the particle with a periodic force. When the forcing period is set at twice the transition mean time, a statistical synchronization between particle motion and forcing phase marks the onset of a stochastic resonance mechanism. Comparisons between experimental results and predictions from the Fick-Jacobs theory and Brownian dynamics simulation reveal significant hydrodynamic effects, which change both resonant amplification and noise level. We further show that hydrodynamic effects can be incorporated into existing theory and simulation by using an experimentally measured particle diffusivity.

PMID:36083679 | DOI:10.1103/PhysRevLett.129.098001

Categories
Nevin Manimala Statistics

Stretching of a Fractal Polymer around a Disc Reveals Kardar-Parisi-Zhang Scaling

Phys Rev Lett. 2022 Aug 26;129(9):097801. doi: 10.1103/PhysRevLett.129.097801.

ABSTRACT

While stretching of a polymer along a flat surface is hardly different from the classical Pincus problem of pulling chain ends in free space, the role of curved geometry in conformational statistics of the stretched chain is an exciting open question. We use scaling analysis and computer simulations to examine stretching of a fractal polymer chain around a disc in 2D (or a cylinder in 3D) of radius R. We reveal that the typical excursions of the polymer away from the surface and curvature-induced correlation length scale as Δ∼R^{β} and S^{*}∼R^{1/z}, respectively, with the Kardar-Parisi-Zhang (KPZ) growth β=1/3 and dynamic exponents z=3/2. Although probability distribution of excursions does not belong to KPZ universality class, the KPZ scaling is independent of the fractal dimension of the polymer and, thus, is universal across classical polymer models, e.g., SAW, randomly branching polymers, crumpled unknotted rings. Additionally, our Letter establishes a mapping between stretched polymers in curved geometry and the Balagurov-Vaks model of random walks among traps.

PMID:36083665 | DOI:10.1103/PhysRevLett.129.097801

Categories
Nevin Manimala Statistics

Emergence of Hilbert Space Fragmentation in Ising Models with a Weak Transverse Field

Phys Rev Lett. 2022 Aug 26;129(9):090602. doi: 10.1103/PhysRevLett.129.090602.

ABSTRACT

The transverse-field Ising model is one of the fundamental models in quantum many-body systems, yet a full understanding of its dynamics remains elusive in higher than one dimension. Here, we show for the first time the breakdown of ergodicity in d-dimensional Ising models with a weak transverse field in a prethermal regime. We demonstrate that novel Hilbert-space fragmentation occurs in the effective nonintegrable model with d≥2 as a consequence of only one emergent global conservation law of the domain wall number. Our results indicate nontrivial initial-state dependence for nonequilibrium dynamics of the Ising models with a weak transverse field.

PMID:36083664 | DOI:10.1103/PhysRevLett.129.090602