Estimation of Personal Symptom Networks Using the Ising Model for Adult Survivors of Childhood Cancer: A Simulation Study with Real-World Data Application.

Purpose: Childhood cancer survivors experience interconnected symptoms, patterns of which can be elucidated by network analysis. However, current symptom networks are constructed based on the average survivors without considering individual heterogeneities. We propose to evaluate personal symptom network estimation using the Ising model with covariates through simulations and estimate personal symptom network for adult childhood cancer survivors. Patients and Methods: We adopted the Ising model with covariates to construct networks by employing logistic regressions for estimating associations between binary symptoms. Simulation experiments assessed the robustness of this method in constructing personal symptom network. Real-world data illustration included 1708 adult childhood cancer survivors from the St. Jude Lifetime Cohort Study (SJLIFE), a retrospective cohort study with prospective follow-up to characterize the etiology and late effects for childhood cancer survivors. Patients' baseline symptoms in 10 domains (cardiac, pulmonary, sensation, nausea, movement, pain, memory, fatigue, anxiety, depression) and individual characteristics (age, sex, race/ethnicity, attained education, personal income, and marital status) were self-reported using survey. Treatment variables (any chemo or radiation therapy) were obtained from medical records. Personal symptom network of 10 domains was estimated using the Ising model, incorporating individual characteristics and treatment data. Results: Simulations confirmed the robustness of the Ising model with covariates in constructing personal symptom networks. Real-world data analysis identified age, sex, race/ethnicity, education, marital status, and treatment (any chemo and radiation therapy) as major factors influencing symptom co-occurrence. Older childhood cancer survivors showed stronger cardiac-fatigue associations. Survivors of racial/ethnic minorities had stronger pain-fatigue associations. Female survivors with above-college education demonstrated stronger pain-anxiety associations. Unmarried survivors who received radiation had stronger association between movement and memory problems. Conclusion: The Ising model with covariates accurately estimates personal symptom networks. Individual heterogeneities exist in symptom co-occurrence patterns for childhood cancer survivors. The estimated personal symptom network offers insights into interconnected symptom experiences.

Distinguishable neural circuit mechanisms associated with the clinical efficacy of rTMS in aMCI patients.

Repetitive transcranial magnetic stimulation is used in early-stage Alzheimer's disease to slow progression, but heterogeneity in response results in different treatment outcomes. The mechanisms underlying this heterogeneity are unclear. This study used resting-state neuroimaging to investigate the variability in episodic memory improvement from angular gyrus repetitive transcranial magnetic stimulation and tracked the neural circuits involved. Thirty-four amnestic mild cognitive impairment patients underwent angular gyrus repetitive transcranial magnetic stimulation (4 weeks, 20 Hz, 100% resting motor threshold) and were divided into high-response and low-response groups based on minimal clinically important differences in auditory verbal learning test scores. Baseline and pre/post-treatment neural circuit activities were compared. Results indicated that the orbital middle frontal gyrus in the orbitofrontal cortex network and the precuneus in the default mode network had higher local activity in the low-response group. After treatment, changes in local and remote connectivity within brain regions of the orbitofrontal cortex, default mode network, visual network, and sensorimotor network showed opposite trends and were related to treatment effects. This suggests that the activity states of brain regions within the orbitofrontal cortex and default mode network could serve as imaging markers for early cognitive compensation in amnestic mild cognitive impairment patients and predict the aftereffects of repetitive transcranial magnetic stimulation response.

Mapping socio-geographical disparities in the occurrence of teenage maternity in Colombia using multilevel analysis of individual heterogeneity and discriminatory accuracy (MAIHDA).

BACKGROUND: The prevalence of teenage pregnancy in Colombia is higher than the worldwide average. The identification of socio-geographical disparities might help to prioritize public health interventions. AIM: To describe variation in the probability of teenage maternity across geopolitical departments and socio-geographical intersectional strata in Colombia. METHODS: A cross-sectional study based on live birth certificates in Colombia. Teenage maternity was defined as a woman giving birth aged 19 or younger. Multilevel analysis of individual heterogeneity and discriminatory accuracy (MAIHDA) was applied using multilevel Poisson and logistic regression. Two different approaches were used: (1) intersectional: using strata defined by the combination of health insurance, region, area of residency, and ethnicity as the second level (2) geographical: using geopolitical departments as the second level. Null, partial, and full models were obtained. General contextual effect (GCE) based on the variance partition coefficient (VPC) was considered as the measure of disparity. Proportional change in variance (PCV) was used to identify the contribution of each variable to the between-strata variation and to identify whether this variation, if any, was due to additive or interaction effects. Residuals were used to identify strata with potential higher-order interactions. RESULTS: The prevalence of teenage mothers in Colombia was 18.30% (95% CI 18.20-18.40). The highest prevalence was observed in Vichada, 25.65% (95% CI: 23.71-27.78), and in the stratum containing mothers with Subsidized/Unaffiliated healthcare insurance, Mestizo, Rural area in the Caribbean region, 29.08% (95% CI 28.55-29.61). The VPC from the null model was 1.70% and 9.16% using the geographical and socio-geographical intersectional approaches, respectively. The higher PCV for the intersectional model was attributed to health insurance. Positive and negative interactions of effects were observed. CONCLUSION: Disparities were observed between intersectional socio-geographical strata but not between geo-political departments. Our results indicate that if resources for prevention are limited, using an intersectional socio-geographical approach would be more effective than focusing on geopolitical departments especially when focusing resources on those groups which show the highest prevalence. MAIHDA could potentially be applied to many other health outcomes where resource decisions must be made.

A quantitative approach to the intersectional study of mental health inequalities during the COVID-19 pandemic in UK young adults.

PURPOSE: Mental health inequalities across social identities/positions during the COVID-19 pandemic have been mostly reported independently from each other or in a limited way (e.g., at the intersection between age and sex or gender). We aim to provide an inclusive socio-demographic mapping of different mental health measures in the population using quantitative methods that are consistent with an intersectional perspective. METHODS: Data included 8,588 participants from two British cohorts (born in 1990 and 2000-2002, respectively), collected in February/March 2021 (during the third UK nationwide lockdown). Measures of anxiety and depressive symptomatology, loneliness, and life satisfaction were analysed using Multilevel Analysis of Individual Heterogeneity and Discriminatory Accuracy (MAIHDA) models. RESULTS: We found evidence of large mental health inequalities across intersectional strata. Large proportions of those inequalities were accounted for by the additive effects of the variables used to define the intersections, with some of the largest gaps associated with sexual orientation (with sexual minority groups showing substantially worse outcomes). Additional inequalities were found by cohort/generation, birth sex, racial/ethnic groups, and socioeconomic position. Intersectional effects were observed mostly in intersections defined by combinations of privileged and marginalised social identities/positions (e.g., lower-than-expected life satisfaction in South Asian men in their thirties from a sexual minority and a disadvantaged childhood social class). CONCLUSION: We found substantial inequalities largely cutting across intersectional strata defined by multiple co-constituting social identities/positions. The large gaps found by sexual orientation extend the existing evidence that sexual minority groups were disproportionately affected by the pandemic. Study implications and limitations are discussed.

Evaluating the importance of individual heterogeneity in reproduction to Weddell seal population dynamics using integral projection models.

Identifying and accounting for unobserved individual heterogeneity in vital rates in demographic models is important for estimating population-level vital rates and identifying diverse life-history strategies, but much less is known about how this individual heterogeneity influences population dynamics. We aimed to understand how the distribution of individual heterogeneity in reproductive and survival rates influenced population dynamics using vital rates from a Weddell seal population by altering the distribution of individual heterogeneity in reproduction, which also altered the distribution of individual survival rates through the incorporation of our estimate of the correlation between the two rates and assessing resulting changes in population growth. We constructed an integral projection model (IPM) structured by age and reproductive state using estimates of vital rates for a long-lived mammal that has recently been shown to exhibit large individual heterogeneity in reproduction. Using output from the IPM, we evaluated how population dynamics changed with different underlying distributions of unobserved individual heterogeneity in reproduction. Results indicate that the changes to the underlying distribution of individual heterogeneity in reproduction cause very small changes in the population growth rate and other population metrics. The largest difference in the estimated population growth rate resulting from changes to the underlying distribution of individual heterogeneity was less than 1%. Our work highlights the differing importance of individual heterogeneity at the population level compared to the individual level. Although individual heterogeneity in reproduction may result in large differences in the lifetime fitness of individuals, changing the proportion of above- or below-average breeders in the population results in much smaller differences in annual population growth rate. For a long-lived mammal with stable and high adult-survival that gives birth to a single offspring, individual heterogeneity in reproduction has a limited effect on population dynamics. We posit that the limited effect of individual heterogeneity on population dynamics may be due to canalization of life-history traits.

Individual life histories: neither slow nor fast, just diverse.

The slow-fast continuum is a commonly used framework to describe variation in life-history strategies across species. Individual life histories have also been assumed to follow a similar pattern, especially in the pace-of-life syndrome literature. However, whether a slow-fast continuum commonly explains life-history variation among individuals within a population remains unclear. Here, we formally tested for the presence of a slow-fast continuum of life histories both within populations and across species using detailed long-term individual-based demographic data for 17 bird and mammal species with markedly different life histories. We estimated adult lifespan, age at first reproduction, annual breeding frequency, and annual fecundity, and identified the main axes of life-history variation using principal component analyses. Across species, we retrieved the slow-fast continuum as the main axis of life-history variation. However, within populations, the patterns of individual life-history variation did not align with a slow-fast continuum in any species. Thus, a continuum ranking individuals from slow to fast living is unlikely to shape individual differences in life histories within populations. Rather, individual life-history variation is likely idiosyncratic across species, potentially because of processes such as stochasticity, density dependence, and individual differences in resource acquisition that affect species differently and generate non-generalizable patterns across species.

Unpacking the lifelong secrets of little penguins: individual quality, energy allocation, and stochasticity in defining fitness.

While the heterogeneity among individuals of a population is more and more documented, questions on the paths through which it arises, particularly whether it is linked to fixed heterogeneity or chance alone, are still widely debated. Here, we tested how individual quality, energy allocation trade-offs, and environmental stochasticity define individual fitness. To do so, we simultaneously investigated the contribution of 18 life-history traits to the fitness of breeding little penguins (Eudyptula minor), using a structural equation model. Fitness was highly variable amongst the 162 birds monitored over their entire lifespan. It increased with the individual penguin's ability to increase (a) the number of breeding events (i.e., living longer, breeding younger, breeding more often, and producing more second clutches) and (b) the breeding success per event through increased foraging performances (i.e., mass gained at sea). While all three processes (stochasticity, individual quality, and allocation trade-offs) affected fitness, interindividual variability in fitness was mainly driven by individual quality, birds consistently breeding earlier in the season and displaying higher foraging efficiency exhibiting higher fitness. Why some birds consistently can perform better at sea and breed earlier remains a question to investigate to understand how selection applies to these traits.

Hurricanes affect diversification among individual life courses of a primate population.

Extreme climatic events may influence individual-level variability in phenotypes, survival and reproduction, and thereby drive the pace of evolution. Climate models predict increases in the frequency of intense hurricanes, but no study has measured their impact on individual life courses within animal populations. We used 45 years of demographic data of rhesus macaques to quantify the influence of major hurricanes on reproductive life courses using multiple metrics of dynamic heterogeneity accounting for life course variability and life-history trait variances. To reduce intraspecific competition, individuals may explore new reproductive stages during years of major hurricanes, resulting in higher temporal variation in reproductive trajectories. Alternatively, individuals may opt for a single optimal life-history strategy due to trade-offs between survival and reproduction. Our results show that heterogeneity in reproductive life courses increased by 4% during years of major hurricanes, despite a 2% reduction in the asymptotic growth rate due to an average decrease in mean fertility and survival by that is, shortened life courses and reduced reproductive output. In agreement with this, the population is expected to achieve stable population dynamics faster after being perturbed by a hurricane ( ρ = 1.512 ; 95% CI: 1.488, 1.538), relative to ordinary years ρ = 1.482 ; 1.475 , 1.490 . Our work suggests that natural disasters force individuals into new demographic roles to potentially reduce competition during unfavourable environments where mean reproduction and survival are compromised. Variance in lifetime reproductive success and longevity are differently affected by hurricanes, and such variability is mostly driven by survival.

Effects of population density and environmental conditions on life-history prevalence in a migratory fish.

Individual variation in life-history traits can have important implications for the ability of populations to respond to environmental variability and change. In migratory animals, flexibility in the timing of life-history events, such as juvenile emigration from natal areas, can influence the effects of population density and environmental conditions on habitat use and population dynamics. We evaluated the functional relationships between population density and environmental covariates and the abundance of juveniles expressing different life-history pathways in a migratory fish, Chinook salmon (Oncorhynchus tshawytscha), in the Wenatchee River basin in Washington State, USA. We found that the abundance of younger emigrants from natal streams was best described by an accelerating or near-linear function of spawners, whereas the abundance of older emigrants was best described by a decelerating function of spawners. This supports the hypothesis that emigration timing varies in response to density in natal areas, with younger-emigrating life-history pathways comprising a larger proportion of emigrants when densities of conspecifics are high. We also observed positive relationships between winter stream discharge and abundance of younger emigrants, supporting the hypothesis that habitat conditions can also influence the prevalence of different life-history pathways. Our results suggest that early emigration, and a resultant increase in the use of downstream rearing habitats, may increase at higher population densities and with greater winter precipitation. Winter precipitation is projected to increase in this system due to climate warming. Characterizing relationships between life-history prevalence and environmental conditions may improve our understanding of species habitat requirements and is a first step in understanding the dynamics of species with diverse life-history strategies. As environmental conditions change-due to climate change, management, or other factors-resultant life-history changes are likely to have important demographic implications that will be challenging to predict when life-history diversity is not accounted for in population models.

LRBmat: A novel gut microbial interaction and individual heterogeneity inference method for colorectal cancer.

The gut microbial community has been shown to play a significant role in various diseases, including colorectal cancer (CRC), which is a major public health concern worldwide. The accurate diagnosis and etiological analysis of CRC are crucial issues. Numerous methods have utilized gut microbiota to address these challenges; however, few have considered the complex interactions and individual heterogeneity of the gut microbiota, which are important issues in genetics and intestinal microbiology, particularly in high-dimensional cases. This paper presents a novel method called Binary matrix based on Logistic Regression (LRBmat) to address these concerns. The binary matrix in LRBmat can directly mitigate or eliminate the influence of heterogeneity, while also capturing information on gut microbial interactions with any order. LRBmat is highly adaptable and can be combined with any machine learning method to enhance its capabilities. The proposed method was evaluated using real CRC data and demonstrated superior classification performance compared to state-of-the-art methods. Furthermore, the association rules extracted from the binary matrix of the real data align well with biological properties and existing literature, thereby aiding in the etiological analysis of CRC.

Adaption, neutrality and life-course diversity.

Heterogeneity among individuals in fitness components is what selection acts upon. Evolutionary theories predict that selection in constant environments acts against such heterogeneity. But observations reveal substantial non-genetic and also non-environmental variability in phenotypes. Here, we examine whether there is a relationship between selection pressure and phenotypic variability by analysing structured population models based on data from a large and diverse set of species. Our findings suggest that non-genetic, non-environmental variation is in general neither truly neutral, selected for, nor selected against. We find much variations among species and populations within species, with mean patterns suggesting nearly neutral evolution of life-course variability. Populations that show greater diversity of life courses do not show, in general, increased or decreased population growth rates. Our analysis suggests we are only at the beginning of understanding the evolution and maintenance of non-genetic non-environmental variation.

Reproductive senescence and mating tactic interact and conflict to drive reproductive success in a passerine.

An understanding of the drivers of individual fitness is a fundamental component of evolutionary ecology and life-history theory. Reproductive senescence, mate and mating tactic choice and latent heterogeneity in individual quality interact to affect individual fitness. We sought to disentangle the effects of these fitness drivers, where longitudinal data are required to understand their respective impacts. We used reproductive allocation and success data from a long-term (1989-2018) study of white-throated dippers Cinclus cinclus in Switzerland to simultaneously examine the effects of female and male age, mating tactic, nest initiation date and individual heterogeneity on reproductive performance. We modelled quadratic and categorical effects of age on reproductive parameters. The probability of polygyny increased with age in both sexes before declining in older age classes. Similarly, hatching probability in monogamous pairs and the number of nestlings hatched in both monogamous and polygynous pairs increased with female age before declining later in life. As predicted, offspring survival in monogamous pairs increased with male age before declining in older age classes, but male age had no effect on offspring survival in polygynous nesting attempts. Our results demonstrate that parental age, mating tactic and individual heterogeneity all affect reproductive success, and that the impacts of senescent decline are expressed across different demographic components as a function of sex-specific senescent decline and mating tactic.

Zero-inflated count distributions for capture-mark-reencounter data.

The estimation of demographic parameters is a key component of evolutionary demography and conservation biology. Capture-mark-recapture methods have served as a fundamental tool for estimating demographic parameters. The accurate estimation of demographic parameters in capture-mark-recapture studies depends on accurate modeling of the observation process. Classic capture-mark-recapture models typically model the observation process as a Bernoulli or categorical trial with detection probability conditional on a marked individual's availability for detection (e.g., alive, or alive and present in a study area). Alternatives to this approach are underused, but may have great utility in capture-recapture studies. In this paper, we explore a simple concept: in the same way that counts contain more information about abundance than simple detection/non-detection data, the number of encounters of individuals during observation occasions contains more information about the observation process than detection/non-detection data for individuals during the same occasion. Rather than using Bernoulli or categorical distributions to estimate detection probability, we demonstrate the application of zero-inflated Poisson and gamma-Poisson distributions. The use of count distributions allows for inference on availability for encounter, as well as a wide variety of parameterizations for heterogeneity in the observation process. We demonstrate that this approach can accurately recover demographic and observation parameters in the presence of individual heterogeneity in detection probability and discuss some potential future extensions of this method.

Causes and consequences of variation in early-life telomere length in a bird metapopulation.

Environmental conditions during early-life development can have lasting effects shaping individual heterogeneity in fitness and fitness-related traits. The length of telomeres, the DNA sequences protecting chromosome ends, may be affected by early-life conditions, and telomere length (TL) has been associated with individual performance within some wild animal populations. Thus, knowledge of the mechanisms that generate variation in TL, and the relationship between TL and fitness, is important in understanding the role of telomeres in ecology and life-history evolution. Here, we investigate how environmental conditions and morphological traits are associated with early-life blood TL and if TL predicts natal dispersal probability or components of fitness in 2746 wild house sparrow (Passer domesticus) nestlings from two populations sampled across 20 years (1994-2013). We retrieved weather data and we monitored population fluctuations, individual survival, and reproductive output using field observations and genetic pedigrees. We found a negative effect of population density on TL, but only in one of the populations. There was a curvilinear association between TL and the maximum daily North Atlantic Oscillation index during incubation, suggesting that there are optimal weather conditions that result in the longest TL. Dispersers tended to have shorter telomeres than non-dispersers. TL did not predict survival, but we found a tendency for individuals with short telomeres to have higher annual reproductive success. Our study showed how early-life TL is shaped by effects of growth, weather conditions, and population density, supporting that environmental stressors negatively affect TL in wild populations. In addition, shorter telomeres may be associated with a faster pace-of-life, as individuals with higher dispersal rates and annual reproduction tended to have shorter early-life TL.

Snared in an Evil Time: How Age-Dependent Environmental and Demographic Variability Contribute to Variance in Lifetime Outcomes.

AbstractTo what degree is lifetime success determined by innate individual quality versus external events and random chance, whether success is measured by lifetime reproductive output, life span, years that a tree spends in the canopy, or some other measure? And how do external events and chance interact with development (survival and growth) to drive success? To answer these questions, we extend our earlier age partitioning of luck in lifetime outcomes in two ways: we incorporate effects of external environmental variation, and we subdivide demographic luck into contributions from survival and growth. Applying our methods to four case studies, we find that luck in survival, in growth, or in environmental variation can all be the dominant driver of success, depending on life history, but variation in individual quality remains a lesser driver. Luck in its various forms is most important at very early ages and again close to the time when individuals typically first begin to be successful (e.g., entering the canopy, reaching reproductive maturity), but different forms of luck peak at different times. For example, a favorable year can boost a tree into the canopy, while luck in survival is required to take full advantage of that fortunate event.

When the going gets tough, the tough get going: Effect of extreme climate on an Antarctic seabird's life history.

Individuals differ in many ways. Most produce few offspring; a handful produce many. Some die early; others live to old age. It is tempting to attribute these differences in outcomes to differences in individual traits, and thus in the demographic rates experienced. However, there is more to individual variation than meets the eye of the biologist. Even among individuals sharing identical traits, life history outcomes (life expectancy and lifetime reproduction) will vary due to individual stochasticity, that is to chance. Quantifying the contributions of heterogeneity and chance is essential to understand natural variability. Interindividual differences vary across environmental conditions, hence heterogeneity and stochasticity depend on environmental conditions. We show that favourable conditions increase the contributions of individual stochasticity, and reduce the contributions of heterogeneity, to variance in demographic outcomes in a seabird population. The opposite is true under poor conditions. This result has important consequence for understanding the ecology and evolution of life history strategies.

Complex Disease Individual Molecular Characterization Using Infinite Sparse Graphical Independent Component Analysis.

Identifying individual mechanisms involved in complex diseases, such as cancer, is essential for precision medicine. Their characterization is particularly challenging due to the unknown relationships of high-dimensional omics data and their inter-patient heterogeneity. We propose to model individual gene expression as a combination of unobserved molecular mechanisms (molecular components) that may differ between the individuals. Considering a baseline molecular profile common to all individuals, these molecular components may represent molecular pathways differing from the population background. We defined an infinite sparse graphical independent component analysis (isgICA) to identify these molecular components. This model relies on double sparseness: the source matrix sparseness defines the subset of genes involved in each molecular component, whereas the weight matrix sparseness identifies the subset of molecular components associated with each patient. As the number of molecular components is unknown but likely high, we simultaneously inferred it and the weight matrix sparseness using the beta-Bernoulli process (BBP). We simulated data from a double sparse ICA with 10/30 components with specific sparseness structures for 100/500 individuals and 500/1000/5000 genes with different noise variance levels to evaluate the reconstruction of the latent structures by our model. For all simulations, the isgICA was able to reconstruct with higher accuracy than 2 state-of-the-art methods (ica and fastICA) the number of components, the weight and source matrix sparsenesses (correlation simulated/estimated >.8). Applying our model to the expression of 1063 genes of 614 breast cancer patients, the isgICA identified 22 components. According to the source matrix, 7 of these 22 components seemed to be specifically related to 3 known molecular pathways with a prognostic effect in early breast cancer (immune system, proliferation, and stroma invasion). This proposed algorithm provides an insight into individual molecular heterogeneity to better understand complex disease mechanisms.

A Single-Cell Raman Spectroscopy Analysis of Bone Marrow Mesenchymal Stem/Stromal Cells to Identify Inter-Individual Diversity.

The heterogeneity of stem cells represents the main challenge in regenerative medicine development. This issue is particularly pronounced when it comes to the use of primary mesenchymal stem/stromal cells (MSCs) due to a lack of identification markers. Considering the need for additional approaches in MSCs characterization, we applied Raman spectroscopy to investigate inter-individual differences between bone marrow MSCs (BM-MSCs). Based on standard biological tests, BM-MSCs of analyzed donors fulfill all conditions for their characterization, while no donor-related specifics were observed in terms of BM-MSCs morphology, phenotype, multilineage differentiation potential, colony-forming capacity, expression of pluripotency-associated markers or proliferative capacity. However, examination of BM-MSCs at a single-cell level by Raman spectroscopy revealed that despite similar biochemical background, fine differences in the Raman spectra of BM-MSCs of each donor can be detected. After extensive principal component analysis (PCA) of Raman spectra, our study revealed the possibility of this method to diversify BM-MSCs populations, whereby the grouping of cell populations was most prominent when cell populations were analyzed in pairs. These results indicate that Raman spectroscopy, as a label-free assay, could have a huge potential in understanding stem cell heterogeneity and sorting cell populations with a similar biochemical background that can be significant for the development of personalized therapy approaches.

Tobacco Taxes as the Unsung Hero: Impact of a Tax Increase on Advancing Sustainable Development in Colombia.

Objective: Tobacco taxes are a well-established cost-effective policy to prevent Noncommunicable Diseases. This paper evaluates the expected effects of a tobacco tax increase on the Sustainable Development Goals in Colombia. Methods: We use microsimulation to build an artificial society that mimics the observed characteristics of Colombia's population, and from there we simulate the behavioral response to a tax increase of COP$4,750 (an increase that has been discussed by policy makers and legislators) and the subsequent effects in all SDGs. Results: The tobacco tax hike reduces the number of smokers (from 4.51 to 3.45 MM smokers) and smoking intensity, resulting in a drop in the number of cigarettes smoked in Colombia (from 332.3 to 215.5 MM of 20-stick packs). Such reduction is expected to decrease premature mortality, healthcare costs, poverty and people facing catastrophic expenditure on healthcare, to increase health, income and gender equity, and to strengthen domestic resource mobilization even in the presence of illicit cigarettes. Conclusion: Tobacco taxes are an effective intervention for public health and a powerful instrument to advance on the 2030 Sustainable Development Agenda. Relevance: A comprehensive analysis of the impact of tobacco taxes on all areas of Sustainable Development is missing in the empirical literature. Such perspective is needed to break the barriers for further tobacco tax increases by gathering wider societal support, especially from stakeholders and key decision makers from development areas other than health. SDG Nr: SDG3 (health), SDG 1 (no poverty), SDG 4 (education), SDG 5 (gender equality), SDG6 (water), SDG10 (inequality), SDG12 (responsible production and consumption), SDG17 (partnerships).

Modeling Rumor Diffusion Process With the Consideration of Individual Heterogeneity: Take the Imported Food Safety Issue as an Example During the COVID-19 Pandemic.

At present, rumors appear frequently in social platforms. The rumor diffusion will cause a great impact on the network order and the stability of the society. So it's necessary to study the diffusion process and develop the rumor control strategies. This article integrates three heterogeneous factors into the SEIR model and designs an individual state transition mode at first. Secondly, based on the influencing factors such as the trust degree among individuals, an individual information interaction mode is constructed. Finally, an improved SEIR model named SEIR-OM model is established, and the diffusion process of rumors are simulated and analyzed. The results show that: (1) when the average value of the interest correlation is greater, the information content deviation is lower, but the rumor diffusion range will be wider. (2) The increase of the average network degree intensifies influence of rumors, but its impact on the diffusion has a peak. (3) Adopting strategies in advance can effectively reduce the influence of rumors. In addition, the government should enforce rumor-refuting strategies right after the event. Also, the number of rumor-refuting individuals must be paid attention to. Finally, the article verifies the rationality and effectiveness of the SEIR-OM model through the real case.