We use cookies to distinguish you from other users and to provide you with a better experience on our websites. Close this message to accept cookies or find out how to manage your cookie settings.
To save content items to your account,
please confirm that you agree to abide by our usage policies.
If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your account.
Find out more about saving content to .
To save content items to your Kindle, first ensure no-reply@cambridge.org
is added to your Approved Personal Document E-mail List under your Personal Document Settings
on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part
of your Kindle email address below.
Find out more about saving to your Kindle.
Note you can select to save to either the @free.kindle.com or @kindle.com variations.
‘@free.kindle.com’ emails are free but can only be saved to your device when it is connected to wi-fi.
‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.
We sought to assess the degree to which environmental risk factors affect CHD prevalence using a case–control study.
Methods:
A hospital-based study was conducted by collecting data from outpatients between January 2016 and January 2021, which included 31 CHD cases and 72 controls from eastern China. Risk ratios were estimated using univariate and multivariate logistic regression models and mediating effect analysis.
Results:
Residential characteristics (usage of cement flooring, odds ratio = 17.04[1.954–148.574], P = 0.01; musty smell, odds ratio = 3.105[1.198–8.051], P = 0.02) and indoor total volatile organic compound levels of participants’ room (odds ratio = 31.846[8.187–123.872, P < 0.001), benzene level (odds ratio = 7.370[2.289–23.726], P = 0.001) increased the risk of CHDs in offspring. And folic acid plays a masking effect, which mitigates the affection of the total volatile organic compound (indirect effect = -0.072[−0.138,-0.033]) and formaldehyde (indirect effect = −0.109[-0.381,-0.006]) levels on the incidence of CHDs. While food intake including milk (odds ratio = 0.396[0.16–0.977], P = 0.044), sea fish (odds ratio = 0.273[0.086–0.867], P = 0.028), and wheat (odds ratio = 0.390[0.154–0.990], P = 0.048) were all protective factors for the occurrence of CHDs. Factors including women reproductive history (history of conception control, odds ratio = 2.648[1.062–6.603], P = 0.037; history of threatened abortion, odds ratio = 2.632[1.005–6.894], P = 0.049; history of dysmenorrhoea (odds ratio = 2.720[1.075–6.878], P = 0.035); sleep status (napping habit during daytime, odds ratio = 0.856[0.355–2.063], P = 0.047; poor sleep quality, odds ratio = 3.180[1.037–9.754], P = 0.043); and work status (working time > 40h weekly, odds ratio = 2.882[1.172–7.086], P = 0.021) also influenced the CHDs incidence to differing degrees.
Conclusion:
Diet habits, nutrients intake, psychological status of pregnant women, and residential air quality were associated with fetal CHDs. Indoor total volatile organic compound content was significantly correlated with CHDs risk, and folic acid may serve as a masking factor that reduce the harmful effects of air pollutants.
The relationship between emotional symptoms and cognitive impairments in major depressive disorder (MDD) is key to understanding cognitive dysfunction and optimizing recovery strategies. This study investigates the relationship between subjective and objective cognitive functions and emotional symptoms in MDD and evaluates their contributions to social functioning recovery.
Methods
The Prospective Cohort Study of Depression in China (PROUD) involved 1,376 MDD patients, who underwent 8 weeks of antidepressant monotherapy with assessments at baseline, week 8, and week 52. Measures included the Hamilton Depression Rating Scale (HAMD-17), Quick Inventory of Depressive Symptomatology-Self Report (QIDS-SR16), Chinese Brief Cognitive Test (C-BCT), Perceived Deficits Questionnaire for Depression-5 (PDQ-D5), and Sheehan Disability Scale (SDS). Cross-lagged panel modeling (CLPM) was used to analyze temporal relationships.
Results
Depressive symptoms and cognitive measures demonstrated significant improvement over 8 weeks (p < 0.001). Baseline subjective cognitive dysfunction predicted depressive symptoms at week 8 (HAMD-17: β = 0.190, 95% CI: 0.108–0.271; QIDS-SR16: β = 0.217, 95% CI: 0.126–0.308). Meanwhile, baseline depressive symptoms (QIDS-SR16) also predicted subsequent subjective cognitive dysfunction (β = 0.090, 95% CI: 0.003-0.177). Recovery of social functioning was driven by improvements in depressive symptoms (β = 0.384, p < 0.0001) and subjective cognition (β = 0.551, p < 0.0001), with subjective cognition contributing more substantially (R2 = 0.196 vs. 0.075).
Conclusions
Subjective cognitive dysfunction is more strongly associated with depressive symptoms and plays a significant role in social functioning recovery, highlighting the need for targeted interventions addressing subjective cognitive deficits in MDD.
This study proposes two novel time-varying model-averaging methods for time-varying parameter regression models. When the number of predictors is small, we propose a novel time-varying complete subset-averaging (TVCSA) procedure, where the optimal time-varying subset size is obtained by minimizing the local leave-h-out cross-validation criterion. The TVCSA method is asymptotically optimal for achieving the lowest possible local mean squared error. When the number of predictors is relatively large, we propose a factor TVCSA method to reduce the computational burden by first reducing the dimension of predictors by extracting a few factors using principal component analysis and then obtaining the TVCSA forecasts from time-varying models with the generated factors. We show that the TVCSA estimator remains asymptotically optimal in the presence of generated factors. Monte Carlo simulation studies have provided favorable evidence for the TVCSA methods relative to the popular model-averaging methods in the literature. Empirical applications to equity premiums and inflation forecasting highlight the practical merits of the proposed methods.
Special economic zones (SEZs) and investment facilitation are among the most discussed policy topics in recent years. They are important investment policy tools in promoting countries’ economic growth. They are closely linked but also have a number of important distinctions. There is a lack of discussion over the nexus between SEZs and investment facilitation. This chapter takes a closer look at the interconnections between SEZs and investment facilitation. Through a brief overview of the development of these two policy instruments, it analyzes their linkages and discusses their areas of divergence. It shows that SEZs and investment facilitation can be complementary and mutually supportive. It advocates enhanced coherence and possible coordination between investment facilitation and SEZ policy schemes.
The whitefly, Bemisia tabaci is a cryptic species complex in which one member, Middle East-Asia Minor 1 (MEAM1) has invaded globally. After invading large countries like Australia, China, and the USA, MEAM1 spread rapidly across each country. In contrast, our analysis of MEAM1 in India showed a very different pattern. Despite the detection of MEAM1 being contemporaneous with invasions in Australia, the USA, and China, MEAM1 has not spread widely and instead remains restricted to the southern regions. An assessment of Indian MEAM1 genetic diversity showed a level of diversity equivalent to that found in its presumed home range and significantly higher than that expected across the invaded range. The high level of diversity and restricted distribution raises the prospect that its home range extends into India. Similarly, while the levels of diversity in Australia and the USA conformed to that expected for the invaded range, China did not. It suggests that China may also be part of its home range. We also observed that diversity across the invaded range was primarily accounted for by a single haplotype, Hap1, which accounted for 79.8% of all records. It was only the invasion of Hap1 that enabled outbreaks to occur and MEAM1’s discovery.
The emotion regulation network (ERN) in the brain provides a framework for understanding the neuropathology of affective disorders. Although previous neuroimaging studies have investigated the neurobiological correlates of the ERN in major depressive disorder (MDD), whether patients with MDD exhibit abnormal functional connectivity (FC) patterns in the ERN and whether the abnormal FC in the ERN can serve as a therapeutic response signature remain unclear.
Methods
A large functional magnetic resonance imaging dataset comprising 709 patients with MDD and 725 healthy controls (HCs) recruited across five sites was analyzed. Using a seed-based FC approach, we first investigated the group differences in whole-brain resting-state FC of the 14 ERN seeds between participants with and without MDD. Furthermore, an independent sample (45 MDD patients) was used to evaluate the relationship between the aforementioned abnormal FC in the ERN and symptom improvement after 8 weeks of antidepressant monotherapy.
Results
Compared to the HCs, patients with MDD exhibited aberrant FC between 7 ERN seeds and several cortical and subcortical areas, including the bilateral middle temporal gyrus, bilateral occipital gyrus, right thalamus, calcarine cortex, middle frontal gyrus, and the bilateral superior temporal gyrus. In an independent sample, these aberrant FCs in the ERN were negatively correlated with the reduction rate of the HAMD17 score among MDD patients.
Conclusions
These results might extend our understanding of the neurobiological underpinnings underlying unadaptable or inflexible emotional processing in MDD patients and help to elucidate the mechanisms of therapeutic response.
The multi-colour complete light curves and low-resolution spectra of two short period eclipsing Am binaries V404 Aur and GW Gem are presented. The stellar atmospheric parameters of the primary stars were derived through the spectra fitting. The observed and TESS-based light curves of them were analysed by using the Wilson-Devinney code. The photometric solutions suggest that both V404 Aur and GW Gem are semi-detached systems with the secondary component filling its critical Roche Lobe, while the former should be a marginal contact binary. The $O-C$ analysis found that the period of V404 Aur is decreasing at a rate of $dP/dt=-1.06(\pm0.01)\times 10^{-7}\,\mathrm{d}\,\mathrm{ yr}^{-1}$, while the period of GW Gem is increasing at $dP/dt=+2.41(\pm0.01)\times 10^{-8} \mathrm{d}\,\mathrm{yr}^{-1}$. The period decrease of V404 Aur may mainly be caused by the combined effects of the angular momentum loss (AML) via an enhanced stellar wind of the more evolved secondary star and mass transfer between two components. The period increase of GW Gem supports the mass transfer from the secondary to the primary. Both targets may be in the broken contact stage predicted by the thermal relaxation oscillations theory and will eventually evolve to the contact stage. We have collected about 54 well-known eclipsing Am binaries with absolute parameters from the literature. The relations of these parameters are summarised. There are some components that have a higher degree of evolution. The majority of their hydrogen shell may have been stripped away and the stellar internal layer exposed. The accretion processes from such evolved components may be very important for the formation of Am peculiarity in binaries.
Establishing the invariance property of an instrument (e.g., a questionnaire or test) is a key step for establishing its measurement validity. Measurement invariance is typically assessed by differential item functioning (DIF) analysis, i.e., detecting DIF items whose response distribution depends not only on the latent trait measured by the instrument but also on the group membership. DIF analysis is confounded by the group difference in the latent trait distributions. Many DIF analyses require knowing several anchor items that are DIF-free in order to draw inferences on whether each of the rest is a DIF item, where the anchor items are used to identify the latent trait distributions. When no prior information on anchor items is available, or some anchor items are misspecified, item purification methods and regularized estimation methods can be used. The former iteratively purifies the anchor set by a stepwise model selection procedure, and the latter selects the DIF-free items by a LASSO-type regularization approach. Unfortunately, unlike the methods based on a correctly specified anchor set, these methods are not guaranteed to provide valid statistical inference (e.g., confidence intervals and p-values). In this paper, we propose a new method for DIF analysis under a multiple indicators and multiple causes (MIMIC) model for DIF. This method adopts a minimal \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$L_1$$\end{document} norm condition for identifying the latent trait distributions. Without requiring prior knowledge about an anchor set, it can accurately estimate the DIF effects of individual items and further draw valid statistical inferences for quantifying the uncertainty. Specifically, the inference results allow us to control the type-I error for DIF detection, which may not be possible with item purification and regularized estimation methods. We conduct simulation studies to evaluate the performance of the proposed method and compare it with the anchor-set-based likelihood ratio test approach and the LASSO approach. The proposed method is applied to analysing the three personality scales of the Eysenck personality questionnaire-revised (EPQ-R).
Although dietary factors have been examined as potential risk factors for liver cancer, the evidence is still inconclusive. Using a diet-wide association analysis, our research evaluated the associations of 126 foods and nutrients on the risk of liver cancer in a Chinese population. We obtained the diet consumption of 72,680 women in the Shanghai Women’s Health Study using baseline dietary questionnaires. The association between each food and nutrient and liver cancer risk was quantified by Cox regression model. A false discovery rate of 0.05 was used to determine the foods and nutrients which need to be verified. Totally 256 incident liver cancer cases were identified in 1,267,391 person-years during the follow-up duration. At the statistical significance level (P ≤ 0.05), higher intakes of cooked wheaten foods, pear, grape and copper were inversely associated with liver cancer risk, while spinach, leafy vegetables, eggplant and carrots showed the positive associations. After considering multiple comparisons, no dietary variable was associated with liver cancer risk. Similar findings were seen in the stratification, secondary and sensitivity analyses. Our findings observed no significant association between dietary factors and liver cancer risk after considering multiple comparisons in Chinese women. More evidence is needed to explore the associations between diet and female liver cancer occurrence.
Soft robots show an advantage when conducting tasks in complex environments due to their enormous flexibility and adaptability. However, soft robots suffer interactions and nonlinear deformation when interacting with soft and fluid materials. The reason behind is the free boundary interactions, which refers to undetermined contact between soft materials, specifically containing nonlinear deformation in air and nonlinear interactions in fluid for soft robot simulation. Therefore, we propose a new approach using material point method (MPM), which can solve the free boundary interactions problem, to simulate soft robots under such environments. The proposed approach can autonomously predict the flexible and versatile behaviors of soft robots. Our approach entails incorporating automatic differentiation into the algorithm of MPM to simplify the computation and implement an efficient implicit time integration algorithm. We perform two groups of experiments with an ordinary pneumatic soft finger in different free boundary interactions. The results indicate that it is possible to simulate soft robots with nonlinear interactions and deformation, and such environmental effects on soft robots can be restored.
The breaking and energy distribution of mode-1 depression internal solitary wave interactions with Gaussian ridges are examined through laboratory experiments. A series of processes, such as shoaling, breaking, transmission and reflection, are captured completely by measuring the velocity field in a large region. It is found that the maximum interface descent ($a_{max}$) during wave shoaling is an important parameter for diagnosing the type of wave–ridge interaction and energy distribution. The wave breaking on the ridge depends on the modified blockage parameter $\zeta _m$, the ratio of the sum of the upper layer depth and $a_{max}$ to the water depth at the top of the ridge. As $\zeta _m$ increases, the interaction type transitions from no breaking to plunging and mixed plunging–collapsing breaking. Within the scope of this experiment, the energy distribution can be characterized solely by $\zeta _m$. The transmission energy decreases monotonically with increasing $\zeta _m$, and there is a linear relationship between $\zeta _m^2$ and the reflection coefficient. The value of $a_{max}$ can be determined from the basic initial parameters of the experiment. Based on the incident wave parameters, the depth of the upper and lower layers, and the topographic parameters, two new simple methods for predicting $a_{max}$ on the ridge are proposed.
This study aimed to assess the relationship between COVID-19 infection-related conditions and depressive symptoms among medical staff after easing the zero-COVID policy in China, and to further examine the mediating role of professional burnout.
Methods
A total of 1716 medical staff from all levels of health care institutions in 16 administrative districts of Beijing, China, were recruited to participate at the end of 2022 in this cross-sectional study. Several multiple linear regressions and mediating effects tests were performed to analyze the data.
Results
At the beginning of the end of the zero-COVID policy in China, 91.84% of respondents reported infection with COVID-19. After adjusting for potential confounding variables, the severity of infection symptoms was significantly positively associated with high levels of depressive symptoms (β = 0.06, P < 0.001), and this association was partially mediated by professional burnout. Specifically, emotional exhaustion (95% CI, 0.131, 0.251) and depersonalization (95% CI, 0.009, 0.043) significantly mediated the association between the severity of infection symptoms and depressive symptoms.
Conclusions
The mental health of medical staff with more severe symptoms of COVID-19 infection should be closely monitored. Also, interventions aimed at reducing emotional exhaustion and depersonalization may effectively reduce their risk of developing depressive symptoms.
In hypersonic flight the shock wave and turbulent boundary layer interaction (STBLI) sharply increases wall heat transfer that intensifies the aerodynamic heating problems. In this work the STBLI is modelled by compression ramp flow with a Mach number of 5, a Reynolds number based on momentum thickness of 4652 and a wall to recovery temperature ratio of 0.5. The aerodynamic heat generation and transport mechanisms are investigated in the interaction based on theoretical analysis and direct numerical simulation (DNS) that agrees with previous studies. A prediction correlation of wall heat flux in STBLI is deduced theoretically and validated by some representative data including the present DNS, which improves the prediction accuracy and can be applied to a wider $Ma$ range compared with the canonical Q-P theory. The correlation indicates that the sharp increase of wall heat transfer in the STBLI can be explained by the boundary layer compression and the convection transport enhancement. Based on the DNS results, the aerodynamic heat generation and transport mechanisms are revealed in the separation, recirculation and reattachment zones in the STBLI. From this perspective, the peak heat flux can be further explained by the enhancement of near-wall turbulent energy dissipation, compression aerodynamic heat generation and the near-wall turbulent transport. The generation and transport of compression aerodynamic heat reveal the underlying mechanism of the strong correlation between the peak heat flux ratios and the pressure ratios in STBLIs.
First-episode schizophrenia (FES) is a progressive psychiatric disorder influenced by genetics, environmental factors, and brain function. The functional gradient deficits of drug-naïve FES and its relationship to gene expression profiles and treatment outcomes are unknown.
Methods
In this study, we engaged a cohort of 116 FES and 100 healthy controls (HC), aged 7 to 30 years, including 15 FES over an 8-week antipsychotic medication regimen. Our examination focused on primary-to-transmodal alterations in voxel-based connection gradients in FES. Then, we employed network topology, Neurosynth, postmortem gene expression, and support vector regression to evaluate integration and segregation functions, meta-analytic cognitive terms, transcriptional patterns, and treatment predictions.
Results
FES displayed diminished global connectome gradients (Cohen's d = 0.32–0.57) correlated with compensatory integration and segregation functions (Cohen's d = 0.31–0.36). Predominant alterations were observed in the default (67.6%) and sensorimotor (21.9%) network, related to high-order cognitive functions. Furthermore, we identified notable overlaps between partial least squares (PLS1) weighted genes and dysregulated genes in other psychiatric conditions. Genes linked with gradient alterations were enriched in synaptic signaling, neurodevelopment process, specific astrocytes, cortical layers (layer II and IV), and developmental phases from late/mid fetal to young adulthood. Additionally, the onset age influenced the severity of FES, with discernible differences in connection gradients between minor- and adult-FES. Moreover, the connectivity gradients of FES at baseline significantly predicted treatment outcomes.
Conclusions
These results offer significant theoretical foundations for elucidating the intricate interplay between macroscopic functional connection gradient changes and microscopic transcriptional patterns during the onset and progression of FES.
Folate metabolism is involved in the development and progression of various cancers. We investigated the association of single nucleotide polymorphisms (SNP) in folate-metabolising genes and their interactions with serum folate concentrations with overall survival (OS) and liver cancer-specific survival (LCSS) of newly diagnosed hepatocellular carcinoma (HCC) patients. We detected the genotypes of six SNP in three genes related to folate metabolism: methylenetetrahydrofolate reductase (MTHFR), 5-methyltetrahydrofolate-homocysteine methyltransferase reductase (MTRR) and 5-methyltetrahydrofolate-homocysteine methyltransferase (MTR). Cox proportional hazard models were used to calculate multivariable-adjusted hazard ratios (HR) and 95 % CI. This analysis included 970 HCC patients with genotypes of six SNP, and 864 of them had serum folate measurements. During a median follow-up of 722 d, 393 deaths occurred, with 360 attributed to HCC. In the fully-adjusted models, the MTRR rs1801394 polymorphism was significantly associated with OS in additive (per G allele: HR = 0·84, 95 % CI: 0·71, 0·99), co-dominant (AG v. AA: HR = 0·77; 95 % CI: 0·62, 0·96) and dominant (AG + GG v. AA: HR = 0·78; 95 % CI: 0·63, 0·96) models. Carrying increasing numbers of protective alleles was linked to better LCSS (HR10–12 v. 2–6 = 0·70; 95 % CI: 0·49, 1·00) and OS (HR10–12 v. 2–6 = 0·67; 95 % CI: 0·47, 0·95). Furthermore, we observed significant interactions on both multiplicative and additive scales between serum folate levels and MTRR rs1801394 polymorphism. Carrying the variant G allele of the MTRR rs1801394 is associated with better HCC prognosis and may enhance the favourable association between higher serum folate levels and improved survival among HCC patients.
The associations between obesity and liver diseases are complex and diverse. To explore the causal relationships between obesity and liver diseases, we applied two-sample Mendelian randomisation (MR) and multivariable MR analysis. The data of exposures (BMI and WHRadjBMI) and outcomes (liver diseases and liver function biomarker) were obtained from the open genome-wide association study database. A two-sample MR study revealed that the genetically predicted BMI and WHRadjBMI were associated with non-alcoholic fatty liver disease, liver fibrosis and autoimmune hepatitis. Obesity was not associated with primary biliary cholangitis, liver failure, liver cell carcinoma, viral hepatitis and secondary malignant neoplasm of liver. A higher WHRadjBMI was associated with higher levels of biomarkers of lipid accumulation and metabolic disorders. These findings indicated independent causal roles of obesity in non-alcoholic fatty liver disease, liver fibrosis and impaired liver metabolic function rather than in viral or autoimmune liver disease.
To evaluate the variations in COVID-19 case fatality rates (CFRs) across different regions and waves, and the impact of public health interventions, social and economic characteristics, and demographic factors on COVID-19 CFRs, we collected data from 30 countries with the highest incidence rate in three waves. We summarized the CFRs of different countries and continents in each wave through meta-analysis. Spearman’s correlation and multiple linear regression were employed to estimate the correlation between influencing factors and reduction rates of CFRs. Significant differences in CFRs were observed among different regions during the three waves (P < 0.001). An association was found between the changes in fully vaccinated rates (rs = 0.41), population density (rs = 0.43), the proportion of individuals over 65 years old (rs = 0.43), and the reduction rates of case fatality rate. Compared to Wave 1, the reduction rates in Wave 2 were associated with population density (β = 0.19, 95%CI: 0.05–0.33) and smoking rates (β = −4.66, 95%CI: −8.98 – −0.33), while in Wave 3 it was associated with booster vaccine rates (β = 0.60, 95%CI: 0.11–1.09) and hospital beds per thousand people (β = 4.15, 95%CI: 1.41–6.89). These findings suggest that the COVID-19 CFRs varied across different countries and waves, and promoting booster vaccinations, increasing hospital bed capacity, and implementing tobacco control measures can help reduce CFRs.
High-risk Human Papillomavirus (HPV) infections are a leading cause of cervical diseases among Han Chinese women of reproductive age. Despite studies like Mai et al. (2021) addressing HPV prevalence in Southern China, awareness remains low, especially in Southwest China. Our study addresses this gap.
Objective:
This hospital-based, retrospective study analyzes the prevalence of high-risk HPV and its association with cervical intraepithelial neoplasia (CIN) among Han Chinese women of reproductive age in Southwest China.
Methods:
Data were collected from 724 women undergoing routine health exams from December 2022 to April 2023. A total of 102 women with high-risk HPV infections were identified. A survey assessed HPV awareness, CIN incidence, and socio-demographic factors influencing awareness.
Results:
Of the 724 women, 102 (14.1%) were diagnosed with high-risk HPV, with HPV-16 being the most common subtype (22.5%). Awareness was significantly lower among unmarried women (OR: 6.632, p = 0.047), those with high school education or less (OR: 20.571, p = 0.003), and rural residents (OR: 19.483, p = 0.020). HPV-16 was detected in 54.55% of women with high-grade CIN.
Conclusion:
There is an urgent need for targeted education and HPV vaccination in Southwest China, particularly for women with lower education, rural residents, and older individuals. Subtype-specific strategies are essential for preventing and managing CIN.
Pharmaceutical distribution routing problem is a key problem for pharmaceutical enterprises, since efficient schedules can enhance resource utilization and reduce operating costs. Meanwhile, it is a complicated combinatorial optimization problem. Existing research mainly focused on delivery route lengths or distribution costs minimization, while seldom considered customer priority and carbon emissions simultaneously. However, considering the customer priority and carbon emissions simultaneously will not only help to enhance customer satisfaction, but also help to reduce the carbon emissions. In this article, we consider the customer priority and carbon emission minimization simultaneously in the pharmaceutical distribution routing problem, the corresponding problem is named pharmaceutical distribution routing problem considering customer priority and carbon emissions. A corresponding mathematical model is formulated, the objectives of which are minimizing fixed cost, refrigeration cost, fuel consumption cost, carbon emission cost, and penalty cost for violating time windows. Moreover, a hybrid genetic algorithm (HGA) is proposed to solve the problem. The framework of the proposed HGA is genetic algorithm (GA), where an effective local search based on variable neighborhood search (VNS) is specially designed and incorporated to improve the intensification abilities. In the proposed HGA, crossover with adaptive probability and mutation with adaptive probability are utilized to enhance the algorithm performance. Finally, the proposed HGA is compared with four optimization algorithms, and experimental results have demonstrated the effectiveness of the HGA.