RESUMO
In the brain, functional connections form a network whose topological organization can be described by graph-theoretic network diagnostics. These include characterizations of the community structure, such as modularity and participation coefficient, which have been shown to change over the course of childhood and adolescence. To investigate if such changes in the functional network are associated with changes in cognitive performance during development, network studies often rely on an arbitrary choice of preprocessing parameters, in particular the proportional threshold of network edges. Because the choice of parameter can impact the value of the network diagnostic, and therefore downstream conclusions, we propose to circumvent that choice by conceptualizing the network diagnostic as a function of the parameter. As opposed to a single value, a network diagnostic curve describes the connectome topology at multiple scales-from the sparsest group of the strongest edges to the entire edge set. To relate these curves to executive function and other covariates, we use scalar-on-function regression, which is more flexible than previous functional data-based models used in network neuroscience. We then consider how systematic differences between networks can manifest in misalignment of diagnostic curves, and consequently propose a supervised curve alignment method that incorporates auxiliary information from other variables. Our algorithm performs both functional regression and alignment via an iterative, penalized, and nonlinear likelihood optimization. The illustrated method has the potential to improve the interpretability and generalizability of neuroscience studies where the goal is to study heterogeneity among a mixture of function- and scalar-valued measures.
RESUMO
MOTIVATION: Multiplexed immunofluorescence (mIF) is an emerging assay for multichannel protein imaging that can decipher cell-level spatial features in tissues. However, existing automated cell phenotyping methods, such as clustering, face challenges in achieving consistency across experiments and often require subjective evaluation. As a result, mIF analyses often revert to marker gating based on manual thresholding of raw imaging data. RESULTS: To address the need for an evaluable semi-automated algorithm, we developed GammaGateR, an R package for interactive marker gating designed specifically for segmented cell-level data from mIF images. Based on a novel closed-form gamma mixture model, GammaGateR provides estimates of marker-positive cell proportions and soft clustering of marker-positive cells. The model incorporates user-specified constraints that provide a consistent but slide-specific model fit. We compared GammaGateR against the newest unsupervised approach for annotating mIF data, employing two colon datasets and one ovarian cancer dataset for the evaluation. We showed that GammaGateR produces highly similar results to a silver standard established through manual annotation. Furthermore, we demonstrated its effectiveness in identifying biological signals, achieved by mapping known spatial interactions between CD68 and MUC5AC cells in the colon and by accurately predicting survival in ovarian cancer patients using the phenotype probabilities as input for machine learning methods. GammaGateR is a highly efficient tool that can improve the replicability of marker gating results, while reducing the time of manual segmentation. AVAILABILITY AND IMPLEMENTATION: The R package is available at https://github.com/JiangmeiRubyXiong/GammaGateR.
Assuntos
Algoritmos , Análise de Célula Única , Humanos , Análise de Célula Única/métodos , Software , Processamento de Imagem Assistida por Computador/métodos , Feminino , Neoplasias Ovarianas/metabolismo , Imunofluorescência/métodos , Biomarcadores/metabolismoRESUMO
Spatial heterogeneity in the tumor microenvironment (TME) plays a critical role in gaining insights into tumor development and progression. Conventional metrics typically capture the spatial differential between TME cellular patterns by either exploring the cell distributions in a pairwise fashion or aggregating the heterogeneity across multiple cell distributions without considering the spatial contribution. As such, none of the existing approaches has fully accounted for the simultaneous heterogeneity caused by both cellular diversity and spatial configurations of multiple cell categories. In this article, we propose an approach to leverage spatial entropy measures at multiple distance ranges to account for the spatial heterogeneity across different cellular organizations. Functional principal component analysis (FPCA) is applied to estimate FPC scores which are then served as predictors in a Cox regression model to investigate the impact of spatial heterogeneity in the TME on survival outcome, potentially adjusting for other confounders. Using a non-small cell lung cancer dataset (n = 153) as a case study, we found that the spatial heterogeneity in the TME cellular composition of CD14+ cells, CD19+ B cells, CD4+ and CD8+ T cells, and CK+ tumor cells, had a significant non-zero effect on the overall survival (p = 0.027). Furthermore, using a publicly available multiplexed ion beam imaging (MIBI) triple-negative breast cancer dataset (n = 33), our proposed method identified a significant impact of cellular interactions between tumor and immune cells on the overall survival (p = 0.046). In simulation studies under different spatial configurations, the proposed method demonstrated a high predictive power by accounting for both clinical effect and the impact of spatial heterogeneity.
RESUMO
Multiplex imaging is a powerful tool to analyze the structural and functional states of cells in their morphological and pathological contexts. However, hypothesis testing with multiplex imaging data is a challenging task due to the extent and complexity of the information obtained. Various computational pipelines have been developed and validated to extract knowledge from specific imaging platforms. A common problem with customized pipelines is their reduced applicability across different imaging platforms: Every multiplex imaging technique exhibits platform-specific characteristics in terms of signal-to-noise ratio and acquisition artifacts that need to be accounted for to yield reliable and reproducible results. We propose a pixel classifier-based image preprocessing step that aims to minimize platform-dependency for all multiplex image analysis pipelines. Signal detection and noise reduction as well as artifact removal can be posed as a pixel classification problem in which all pixels in multiplex images can be assigned to two general classes of either I) signal of interest or II) artifacts and noise. The resulting feature representation maps contain pixel-scale representations of the input data, but exhibit significantly increased signal-to-noise ratios with normalized pixel values as output data. We demonstrate the validity of our proposed image preprocessing approach by comparing the results of two well-accepted and widely-used image analysis pipelines.
Assuntos
Processamento de Imagem Assistida por Computador , Tomografia Computadorizada por Raios X , Processamento de Imagem Assistida por Computador/métodos , Tomografia Computadorizada por Raios X/métodos , Artefatos , Razão Sinal-Ruído , AlgoritmosRESUMO
MOTIVATION: Multiplexed imaging is a nascent single-cell assay with a complex data structure susceptible to technical variability that disrupts inference. These in situ methods are valuable in understanding cell-cell interactions, but few standardized processing steps or normalization techniques of multiplexed imaging data are available. RESULTS: We implement and compare data transformations and normalization algorithms in multiplexed imaging data. Our methods adapt the ComBat and functional data registration methods to remove slide effects in this domain, and we present an evaluation framework to compare the proposed approaches. We present clear slide-to-slide variation in the raw, unadjusted data and show that many of the proposed normalization methods reduce this variation while preserving and improving the biological signal. Furthermore, we find that dividing multiplexed imaging data by its slide mean, and the functional data registration methods, perform the best under our proposed evaluation framework. In summary, this approach provides a foundation for better data quality and evaluation criteria in multiplexed imaging. AVAILABILITY AND IMPLEMENTATION: Source code is provided at: https://github.com/statimagcoll/MultiplexedNormalization and an R package to implement these methods is available here: https://github.com/ColemanRHarris/mxnorm. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Assuntos
Algoritmos , Software , ImunofluorescênciaRESUMO
The tumor microenvironment (TME), which characterizes the tumor and its surroundings, plays a critical role in understanding cancer development and progression. Recent advances in imaging techniques enable researchers to study spatial structure of the TME at a single-cell level. Investigating spatial patterns and interactions of cell subtypes within the TME provides useful insights into how cells with different biological purposes behave, which may consequentially impact a subject's clinical outcomes. We utilize a class of well-known spatial summary statistics, the K-function and its variants, to explore inter-cell dependence as a function of distances between cells. Using techniques from functional data analysis, we introduce an approach to model the association between these summary spatial functions and subject-level outcomes, while controlling for other clinical scalar predictors such as age and disease stage. In particular, we leverage the additive functional Cox regression model (AFCM) to study the nonlinear impact of spatial interaction between tumor and stromal cells on overall survival in patients with non-small cell lung cancer, using multiplex immunohistochemistry (mIHC) data. The applicability of our approach is further validated using a publicly available multiplexed ion beam imaging (MIBI) triple-negative breast cancer dataset.
Assuntos
Carcinoma Pulmonar de Células não Pequenas , Neoplasias Pulmonares , Ciência de Dados , Humanos , Imuno-Histoquímica , Microambiente TumoralRESUMO
"Smart"-scales are a new tool for frequent monitoring of weight change as well as weigh-in behavior. These scales give researchers the opportunity to discover patterns in the frequency that individuals weigh themselves over time, and how these patterns are associated with overall weight loss. Our motivating data come from an 18-month behavioral weight loss study of 55 adults classified as overweight or obese who were instructed to weigh themselves daily. Adherence to daily weigh-in routines produces a binary times series for each subject, indicating whether a participant weighed in on a given day. To characterize weigh-in by time-invariant patterns rather than overall adherence, we propose using hierarchical clustering with dynamic time warping (DTW). We perform an extensive simulation study to evaluate the performance of DTW compared to Euclidean and Jaccard distances to recover underlying patterns in adherence time series. In addition, we compare cluster performance using cluster validation indices (CVIs) under the single, average, complete, and Ward linkages and evaluate how internal and external CVIs compare for clustering binary time series. We apply conclusions from the simulation to cluster our real data and summarize observed weigh-in patterns. Our analysis finds that the adherence trajectory pattern is significantly associated with weight loss.
Assuntos
Obesidade , Redução de Peso , Adulto , Humanos , Fatores de Tempo , Simulação por Computador , Análise por ConglomeradosRESUMO
OBJECTIVE: Computing phenotypes that provide high-fidelity, time-dependent characterizations and yield personalized interpretations is challenging, especially given the complexity of physiological and healthcare systems and clinical data quality. This paper develops a methodological pipeline to estimate unmeasured physiological parameters and produce high-fidelity, personalized phenotypes anchored to physiological mechanics from electronic health record (EHR). METHODS: A methodological phenotyping pipeline is developed that computes new phenotypes defined with unmeasurable computational biomarkers quantifying specific physiological properties in real time. Working within the inverse problem framework, this pipeline is applied to the glucose-insulin system for ICU patients using data assimilation to estimate an established mathematical physiological model with stochastic optimization. This produces physiological model parameter vectors of clinically unmeasured endocrine properties, here insulin secretion, clearance, and resistance, estimated for individual patient. These physiological parameter vectors are used as inputs to unsupervised machine learning methods to produce phenotypic labels and discrete physiological phenotypes. These phenotypes are inherently interpretable because they are based on parametric physiological descriptors. To establish potential clinical utility, the computed phenotypes are evaluated with external EHR data for consistency and reliability and with clinician face validation. RESULTS: The phenotype computation was performed on a cohort of 109 ICU patients who received no or short-acting insulin therapy, rendering continuous and discrete physiological phenotypes as specific computational biomarkers of unmeasured insulin secretion, clearance, and resistance on time windows of three days. Six, six, and five discrete phenotypes were found in the first, middle, and last three-day periods of ICU stays, respectively. Computed phenotypic labels were predictive with an average accuracy of 89%. External validation of discrete phenotypes showed coherence and consistency in clinically observable differences based on laboratory measurements and ICD 9/10 codes and clinical concordance from face validity. A particularly clinically impactful parameter, insulin secretion, had a concordance accuracy of 83%±27%. CONCLUSION: The new physiological phenotypes computed with individual patient ICU data and defined by estimates of mechanistic model parameters have high physiological fidelity, are continuous, time-specific, personalized, interpretable, and predictive. This methodology is generalizable to other clinical and physiological settings and opens the door for discovering deeper physiological information to personalize medical care.
Assuntos
Algoritmos , Registros Eletrônicos de Saúde , Humanos , Reprodutibilidade dos Testes , Fenótipo , Biomarcadores , Unidades de Terapia IntensivaRESUMO
Uniform actin filament length is required for synchronized contraction of skeletal muscle. In myopathies linked to mutations in tropomyosin (Tpm) genes, irregular thin filaments are a common feature, which may result from defects in length maintenance mechanisms. The current work investigated the effects of the myopathy-causing p.R91C variant in Tpm3.12, a tropomyosin isoform expressed in slow-twitch muscle fibers, on the regulation of actin severing and depolymerization by cofilin-2. The affinity of cofilin-2 for F-actin was not significantly changed by either Tpm3.12 or Tpm3.12-R91C, though it increased two-fold in the presence of troponin (without Ca2+). Saturation of the filament with cofilin-2 removed both Tpm variants from the filament, although Tpm3.12-R91C was more resistant. In the presence of troponin (±Ca2+), Tpm remained on the filament, even at high cofilin-2 concentrations. Both Tpm3.12 variants inhibited filament severing and depolymerization by cofilin-2. However, the inhibition was more efficient in the presence of Tpm3.12-R91C, indicating that the pathogenic variant impaired cofilin-2-dependent actin filament turnover. Troponin (±Ca2+) further inhibited but did not completely stop cofilin-2-dependent actin severing and depolymerization.
Assuntos
Doenças Musculares , Tropomiosina , Humanos , Citoesqueleto de Actina , Actinas/genética , Cofilina 2/genética , Doenças Musculares/genética , Mutação , Tropomiosina/genética , Troponina/genéticaRESUMO
By collecting data continuously over 24 hours, accelerometers and other wearable devices can provide novel insights into circadian rhythms and their relationship to human health. Existing approaches for analyzing diurnal patterns using these data, including the cosinor model and functional principal components analysis, have revealed and quantified population-level diurnal patterns, but considerable subject-level variability remained uncaptured in features such as wake/sleep times and activity intensity. This remaining informative variability could provide a better understanding of chronotypes, or behavioral manifestations of one's underlying 24-hour rhythm. Curve registration, or alignment, is a technique in functional data analysis that separates "vertical" variability in activity intensity from "horizontal" variability in time-dependent markers like wake and sleep times; this data-driven approach is well-suited to studying chronotypes using accelerometer data. We develop a parametric registration framework for 24-hour accelerometric rest-activity profiles represented as dichotomized into epoch-level states of activity or rest. Specifically, we estimate subject-specific piecewise linear time-warping functions parametrized with a small set of parameters. We apply this method to data from the Baltimore Longitudinal Study of Aging and illustrate how estimated parameters give a more flexible quantification of chronotypes compared to traditional approaches.
RESUMO
The ability of individuals to engage in physical activity is a critical component of overall health and quality of life. However, there is a natural decline in physical activity associated with the aging process. Establishing normative trends of physical activity in aging populations is essential to developing public health guidelines and informing clinical perspectives regarding individuals' levels of physical activity. Beyond overall quantity of physical activity, patterns regarding the timing of activity provide additional insights into latent health status. Wearable accelerometers, paired with statistical methods from functional data analysis, provide the means to estimate diurnal patterns in physical activity. To date, these methods have been only applied to study aging trends in populations based in the United States. Here, we apply curve registration and functional regression to 24 h activity profiles for 88,793 men (N = 39,255) and women (N = 49,538) ages 42-78 from the UK Biobank accelerometer study to understand how physical activity patterns vary across ages and by gender. Our analysis finds that daily patterns in both the volume of physical activity and probability of being active change with age, and that there are marked gender differences in these trends. This work represents the largest-ever population analyzed using tools of this kind, and suggest that aging trends in physical activity are reproducible in different populations across countries.
Assuntos
Bancos de Espécimes Biológicos , Exercício Físico , Qualidade de Vida , Adulto , Idoso , Feminino , Humanos , Masculino , Pessoa de Meia-Idade , Reino Unido , Articulação do PunhoRESUMO
Reaching movements, as a basic yet complex motor behavior, are a foundational model system in neuroscience. In particular, there has been a significant recent expansion of investigation into the neural circuit mechanisms of reach behavior in mice. Nevertheless, quantification of mouse reach kinematics remains lacking, limiting comparison to the primate literature. In this study, we quantitatively demonstrate the homology of mouse reach kinematics to primate reach and also discover novel late-phase correlational structure that implies online control. Overall, our results highlight the decelerative phase of reach as important in driving successful outcome. Specifically, we develop and implement a novel statistical machine-learning algorithm to identify kinematic features associated with successful reaches and find that late-phase kinematics are most predictive of outcome, signifying online reach control as opposed to preplanning. Moreover, we identify and characterize late-phase kinematic adjustments that are yoked to midflight position and velocity of the limb, allowing for dynamic correction of initial variability, with head-fixed reaches being less dependent on position in comparison to freely behaving reaches. Furthermore, consecutive reaches exhibit positional error correction but not hot-handedness, implying opponent regulation of motor variability. Overall, our results establish foundational mouse reach kinematics in the context of neuroscientific investigation, characterizing mouse reach production as an active process that relies on dynamic online control mechanisms.NEW & NOTEWORTHY Mice use reaching movements to grasp and manipulate objects in their environment, similar to primates. To better establish mouse reach as a model for motor control, we implement several analytical frameworks, from basic kinematic relationships to statistical machine learning, to quantify mouse reach, finding many canonical features of primate reaches are conserved in mice, as well as evidence for midflight course corrections, expanding the utility of mouse reach paradigms for motor control studies.
Assuntos
Movimento , Animais , Fenômenos Biomecânicos , Feminino , Aprendizado de Máquina , Masculino , Camundongos Endogâmicos C57BLRESUMO
We introduce a novel method for separating amplitude and phase variability in exponential family functional data. Our method alternates between two steps: the first uses generalized functional principal components analysis to calculate template functions, and the second estimates smooth warping functions that map observed curves to templates. Existing approaches to registration have primarily focused on continuous functional observations, and the few approaches for discrete functional data require a pre-smoothing step; these methods are frequently computationally intensive. In contrast, we focus on the likelihood of the observed data and avoid the need for preprocessing, and we implement both steps of our algorithm in a computationally efficient way. Our motivation comes from the Baltimore Longitudinal Study on Aging, in which accelerometer data provides valuable insights into the timing of sedentary behavior. We analyze binary functional data with observations each minute over 24 hours for 592 participants, where values represent activity and inactivity. Diurnal patterns of activity are obscured due to misalignment in the original data but are clear after curves are aligned. Simulations designed to mimic the application indicate that the proposed methods outperform competing approaches in terms of estimation accuracy and computational efficiency. Code for our method and simulations is publicly available.
Assuntos
Interpretação Estatística de Dados , Análise de Componente Principal/métodos , Tempo , Algoritmos , Gráficos por Computador/estatística & dados numéricos , Simulação por Computador/estatística & dados numéricos , Humanos , Estudos Longitudinais , Atividade Motora , Tamanho da AmostraRESUMO
OBJECTIVES: To assess the association between medical marijuana laws (MMLs) and the odds of a positive opioid test, an indicator for prior use. METHODS: We analyzed 1999-2013 Fatality Analysis Reporting System (FARS) data from 18 states that tested for alcohol and other drugs in at least 80% of drivers who died within 1 hour of crashing (n = 68 394). Within-state and between-state comparisons assessed opioid positivity among drivers crashing in states with an operational MML (i.e., allowances for home cultivation or active dispensaries) versus drivers crashing in states before a future MML was operational. RESULTS: State-specific estimates indicated a reduction in opioid positivity for most states after implementation of an operational MML, although none of these estimates were significant. When we combined states, we observed no significant overall association (odds ratio [OR] = 0.79; 95% confidence interval [CI] = 0.61, 1.03). However, age-stratified analyses indicated a significant reduction in opioid positivity for drivers aged 21 to 40 years (OR = 0.50; 95% CI = 0.37, 0.67; interaction P < .001). CONCLUSIONS: Operational MMLs are associated with reductions in opioid positivity among 21- to 40-year-old fatally injured drivers and may reduce opioid use and overdose.
Assuntos
Acidentes de Trânsito/mortalidade , Acidentes de Trânsito/estatística & dados numéricos , Analgésicos Opioides/sangue , Condução de Veículo/estatística & dados numéricos , Maconha Medicinal , Adolescente , Adulto , Distribuição por Idade , Concentração Alcoólica no Sangue , Feminino , Humanos , Masculino , Prevalência , Adulto JovemRESUMO
OBJECTIVE: Macrophage activation syndrome (MAS) is a devastating cytokine storm syndrome complicating many inflammatory diseases and characterized by fever, pancytopenia, and systemic inflammation. It is clinically similar to hemophagocytic lymphohistiocytosis (HLH), which is caused by viral infection of a host with impaired cellular cytotoxicity. Murine models of MAS and HLH illustrate that interferon-γ (IFNγ) is the driving stimulus for hemophagocytosis and immunopathology. This study was undertaken to investigate the inflammatory contributors to a murine model of Toll-like receptor 9 (TLR-9)-induced fulminant MAS. METHODS: Wild-type, transgenic, and cytokine-inhibited mice were treated with an IL-10 receptor blocking antibody and a TLR-9 agonist, and parameters of MAS were evaluated. RESULTS: Fulminant MAS was characterized by dramatic elevations in IFNγ, IL-12, and IL-6 levels. Increased serum IFNγ levels were associated with enhanced IFNγ production within some hepatic cell populations but also with decreased numbers of IFNγ-positive cells. Surprisingly, IFNγ-knockout mice developed immunopathology and hemophagocytosis comparable to that seen in wild-type mice. However, IFNγ-knockout mice did not become anemic and had greater numbers of splenic erythroid precursors. IL-12 neutralization phenocopied disease in IFNγ-knockout mice. Interestingly, type I IFNs contributed to the severity of hypercytokinemia and weight loss, but their absence did not otherwise affect MAS manifestations. CONCLUSION: These data demonstrate that both fulminant MAS and hemophagocytosis can arise independently of IFNγ, IL-12, or type I IFNs. They also suggest that IFNγ-mediated dyserythropoiesis, not hemophagocytosis, is the dominant cause of anemia in fulminant TLR-9-induced MAS. Thus, our data establish a novel mechanism for the acute anemia of inflammation, but suggest that a variety of triggers can result in hemophagocytic disease.
Assuntos
Anemia/fisiopatologia , Eritropoese/fisiologia , Interferon gama/fisiologia , Linfo-Histiocitose Hemofagocítica/fisiopatologia , Síndrome de Ativação Macrofágica/fisiopatologia , Anemia/etiologia , Animais , Modelos Animais de Doenças , Mediadores da Inflamação/imunologia , Interferon Tipo I/fisiologia , Interferon gama/genética , Interleucina-12/fisiologia , Linfo-Histiocitose Hemofagocítica/complicações , Linfo-Histiocitose Hemofagocítica/imunologia , Síndrome de Ativação Macrofágica/complicações , Síndrome de Ativação Macrofágica/imunologia , Camundongos , Camundongos Knockout , Transdução de Sinais/fisiologia , Receptor Toll-Like 9/agonistasRESUMO
Immune modulation is considered a hallmark of cancer initiation and progression, with immune cell density being consistently associated with clinical outcomes of individuals with cancer. Multiplex immunofluorescence (mIF) microscopy combined with automated image analysis is a novel and increasingly used technique that allows for the assessment and visualization of the tumor microenvironment (TME). Recently, application of this new technology to tissue microarrays (TMAs) or whole tissue sections from large cancer studies has been used to characterize different cell populations in the TME with enhanced reproducibility and accuracy. Generally, mIF data has been used to examine the presence and abundance of immune cells in the tumor and stroma compartments; however, this aggregate measure assumes uniform patterns of immune cells throughout the TME and overlooks spatial heterogeneity. Recently, the spatial contexture of the TME has been explored with a variety of statistical methods. In this PSB workshop, speakers will present some of the state-of-the-art statistical methods for assessing the TIME from mIF data.
Assuntos
Biologia Computacional , Neoplasias , Humanos , Reprodutibilidade dos Testes , Microambiente TumoralRESUMO
Emotional experience is central to a fulfilling life. Although exposure to negative experiences is inevitable, an individual's emotion regulation response may buffer against psychopathology. Identification of neural activation patterns associated with emotion regulation via an fMRI task is a promising and non-invasive means of furthering our understanding of the how the brain engages with negative experiences. Prior work has applied multivariate pattern analysis to identify signatures of response to negative emotion-inducing images; we adapt these techniques to establish novel neural signatures associated with conscious efforts to modulate emotional response. We model voxel-level activation via LASSO principal components regression and linear discriminant analysis to predict if a subject was engaged in emotion regulation and to identify brain regions which define this emotion regulation signature. We train our models using 82 participants and evaluate them on a holdout sample of 40 participants, demonstrating an accuracy up to 82.5% across three classes. Our results suggest that emotion regulation produces a unique signature that is differentiable from passive viewing of negative and neutral imagery.
Assuntos
Regulação Emocional , Humanos , Emoções/fisiologia , Encéfalo/diagnóstico por imagem , Encéfalo/fisiologia , Mapeamento Encefálico , Imageamento por Ressonância MagnéticaRESUMO
Introduction: Given the traffic safety and occupational injury prevention implications associated with cannabis impairment, there is a need for objective and validated measures of recent cannabis use. Pupillary light response may offer an approach for detection. Method: Eighty-four participants (mean age: 32, 42% female) with daily, occasional, and no-use cannabis use histories participated in pupillary light response tests before and after smoking cannabis ad libitum or relaxing for 15 min (no use). The impact of recent cannabis consumption on trajectories of the pupillary light response was modeled using functional data analysis tools. Logistic regression models for detecting recent cannabis use were compared, and average pupil trajectories across cannabis use groups and times since light test administration were estimated. Results: Models revealed small, significant differences in pupil response to light after cannabis use comparing the occasional use group to the no-use control group, and similar statistically significant differences in pupil response patterns comparing the daily use group to the no-use comparison group. Trajectories of pupillary light response estimated using functional data analysis found that acute cannabis smoking was associated with less initial and sustained pupil constriction compared to no cannabis smoking. Conclusion: These analyses show the promise of pairing pupillary light response and functional data analysis methods to assess recent cannabis use.
RESUMO
Introduction: Studies show that acute cannabis use significantly increases heart rate (HR) and mildly raises blood pressure in the minutes following smoked or inhaled use of cannabis. However, less is known about how the THC concentration of the product or an individual's frequency of use (i.e., tolerance) may affect the magnitude of the change in HR. It is also relatively unexamined how the physical effects of increased HR after acute cannabis use relate to self-reported drug effects or blood THC levels. Aims: To describe the relationship between THC concentration of product used, self-reported subjective intoxication, THC blood levels, and frequency of cannabis use with the change in HR after acute cannabis use. Materials and Methods: Participants (n = 140) were given 15 min to smoke self-supplied cannabis ad libitum, HR was measured at baseline and an average of 2 min post-cannabis smoking. The ARCI-Marijuana scale and Visual Analog Scales (VAS) were administered, and blood samples were taken at both time points. Participants were asked about their frequency of use. Information about the product used was recorded from the package. Linear regression was used to analyze the relationship between changes in HR (post-pre cannabis use) and post-cannabis use HR, blood THC concentration, THC product concentration, frequency of use, and self-reported drug effect. Results: There was a significantly higher HR among those who smoked cannabis compared to the controls (p < 0.001), which did not significantly differ by frequency of use (p = 0.18). Higher concentration THC (extract) products did not produce a significantly different HR than lower concentration (flower) products (p = 0.096). VAS score was associated with an increase in HR (p < 0.05). Overall, blood THC levels were not significantly related to the change in HR (p = 0.69); however, when probed, there was a slight positive association among the occasional use group only. Discussion: Cardiovascular effects of cannabis consumption may not be as subject to tolerance with daily cannabis use and do not directly increase with THC concentration of the product. This is a departure from other effects (i.e., cognitive, subjective drug effects) where tolerance is well established. These findings also suggest that, at least among those with daily use, higher concentration THC products (>60%) do not necessarily produce cardiovascular physiological effects that are significantly more robust than lower concentration (<20%) products.
RESUMO
In the brain, functional connections form a network whose topological organization can be described by graph-theoretic network diagnostics. These include characterizations of the community structure, such as modularity and participation coefficient, which have been shown to change over the course of childhood and adolescence. To investigate if such changes in the functional network are associated with changes in cognitive performance during development, network studies often rely on an arbitrary choice of pre-processing parameters, in particular the proportional threshold of network edges. Because the choice of parameter can impact the value of the network diagnostic, and therefore downstream conclusions, we propose to circumvent that choice by conceptualizing the network diagnostic as a function of the parameter. As opposed to a single value, a network diagnostic curve describes the connectome topology at multiple scales-from the sparsest group of the strongest edges to the entire edge set. To relate these curves to executive function and other covariates, we use scalar-on-function regression, which is more flexible than previous functional data-based models used in network neuroscience. We then consider how systematic differences between networks can manifest in misalignment of diagnostic curves, and consequently propose a supervised curve alignment method that incorporates auxiliary information from other variables. Our algorithm performs both functional regression and alignment via an iterative, penalized, and nonlinear likelihood optimization. The illustrated method has the potential to improve the interpretability and generalizability of neuroscience studies where the goal is to study heterogeneity among a mixture of function- and scalar-valued measures.