In Situ Capture of Chromatin Interactions by Biotinylated dCas9.

Cis-regulatory elements (CREs) are commonly recognized by correlative chromatin features, yet the molecular composition of the vast majority of CREs in chromatin remains unknown. Here, we describe a CRISPR affinity purification in situ of regulatory elements (CAPTURE) approach to unbiasedly identify locus-specific chromatin-regulating protein complexes and long-range DNA interactions. Using an in vivo biotinylated nuclease-deficient Cas9 protein and sequence-specific guide RNAs, we show high-resolution and selective isolation of chromatin interactions at a single-copy genomic locus. Purification of human telomeres using CAPTURE identifies known and new telomeric factors. In situ capture of individual constituents of the enhancer cluster controlling human ß-globin genes establishes evidence for composition-based hierarchical organization. Furthermore, unbiased analysis of chromatin interactions at disease-associated cis-elements and developmentally regulated super-enhancers reveals spatial features that causally control gene transcription. Thus, comprehensive and unbiased analysis of locus-specific regulatory composition provides mechanistic insight into genome structure and function in development and disease.

Author Correction: CPS1 maintains pyrimidine pools and DNA synthesis in KRAS/LKB1-mutant lung cancer cells.

Further analysis has revealed that the signal reported in Extended Data Fig. 1c of this Letter is attributed to phosphorylethanolamine, not carbamoyl phosphate. A newly developed derivatization method revealed that the level of carbamoyl phosphate in these NSCLC extracts is below the detection threshold of approximately 10 nanomoles. These findings do not alter the overall conclusions of the Letter; see associated Amendment for full details. The Letter has not been corrected online.

Annotation and evaluation of base editing outcomes in multiple cell types using CRISPRbase.

CRISPR-Cas base editing (BE) system is a powerful tool to expand the scope and efficiency of genome editing with single-nucleotide resolution. The editing efficiency, product purity, and off-target effect differ among various BE systems. Herein, we developed CRISPRbase (http://crisprbase.maolab.org), by integrating 1 252 935 records of base editing outcomes in more than 50 cell types from 17 species. CRISPRbase helps to evaluate the putative editing precision of different BE systems by integrating multiple annotations, functional predictions and a blasting system for single-guide RNA sequences. We systematically assessed the editing window, editing efficiency and product purity of various BE systems. Intensive efforts were focused on increasing the editing efficiency and product purity of base editors since the byproduct could be detrimental in certain applications. Remarkably, more than half of cancer-related off-target mutations were non-synonymous and extremely damaging to protein functions in most common tumor types. Luckily, most of these cancer-related mutations were passenger mutations (4840/5703, 84.87%) rather than cancer driver mutations (863/5703, 15.13%), indicating a weak effect of off-target mutations on carcinogenesis. In summary, CRISPRbase is a powerful and convenient tool to study the outcomes of different base editors and help researchers choose appropriate BE designs for functional studies.

In vivo proximity proteomics uncovers palmdelphin (PALMD) as a Z-disc-associated mitigator of isoproterenol-induced cardiac injury.

Z-discs are core ultrastructural organizers of cardiomyocytes that modulate many facets of cardiac pathogenesis. Yet a comprehensive proteomic atlas of Z-disc-associated components remain incomplete. Here, we established an adeno-associated virus (AAV)-delivered, cardiomyocyte-specific, proximity-labeling approach to characterize the Z-disc proteome in vivo. We found palmdelphin (PALMD) as a novel Z-disc-associated protein in both adult murine cardiomyocytes and human pluripotent stem cell-derived cardiomyocytes. Germline and cardiomyocyte-specific Palmd knockout mice were grossly normal at baseline but exhibited compromised cardiac hypertrophy and aggravated cardiac injury upon long-term isoproterenol treatment. By contrast, cardiomyocyte-specific PALMD overexpression was sufficient to mitigate isoproterenol-induced cardiac injury. PALMD ablation perturbed the transverse tubule (T-tubule)-sarcoplasmic reticulum (SR) ultrastructures, which formed the Z-disc-associated junctional membrane complex (JMC) essential for calcium handling and cardiac function. These phenotypes were associated with the reduction of nexilin (NEXN), a crucial Z-disc-associated protein that is essential for both Z-disc and JMC structures and functions. PALMD interacted with NEXN and enhanced its protein stability while the Nexn mRNA level was not affected. AAV-based NEXN addback rescued the exacerbated cardiac injury in isoproterenol-treated PALMD-depleted mice. Together, this study discovered PALMD as a potential target for myocardial protection and highlighted in vivo proximity proteomics as a powerful approach to nominate novel players regulating cardiac pathogenesis.

Identifying modifications on DNA-bound histones with joint deep learning of multiple binding sites in DNA sequence.

MOTIVATION: Histone modifications are epigenetic markers that impact gene expression by altering the chromatin structure or recruiting histone modifiers. Their accurate identification is key to unraveling the mechanisms by which they regulate gene expression. However, the solutions for this task can be improved by exploiting multiple relationships from dataset and exploring designs of learning models, for example jointly learning technology. RESULTS: This article proposes a deep learning-based multi-objective computational approach, iHMnBS, to identify which of the seven typical histone modifications a DNA sequence may choose to bind, and which parts of the DNA sequence bind to them. iHMnBS employs a customized dataset that allows the marking of modifications contained in histones that may bind to any position in the DNA sequence. iHMnBS tries to mine the information implicit in this richer data by means of deep neural networks. In comprehensive comparisons, iHMnBS outperforms a baseline method, and the probability of binding to modified histones assigned to a representative nucleotide of a DNA sequence can serve as a reference for biological experiments. Since the interaction between transcription factors and histone modifications has an important role in gene expression, we extracted a number of sequence patterns that may bind to transcription factors, and explored their possible impact on disease. AVAILABILITY AND IMPLEMENTATION: The source code is available at https://github.com/lennylv/iHMnBS. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The emergence of influenza B as a major respiratory pathogen in the absence of COVID-19 during the 2021-2022 flu season in China.

BACKGROUND: The emergence of COVID-19 and the implementation of preventive measures and behavioral changes have led to a significant decrease in the prevalence of other respiratory viruses. However, the manner in which seasonal viruses will reemerge in the absence of COVID-19-related restrictions remains unknown. METHODS: Patients presenting with influenza-like illness in two hospitals in Beijing were subjected to testing for COVID-19, influenza A, and influenza B to determine the causative agent for viral infections. The prevalence of influenza B across China was confirmed using data from the Centers for Disease Control, China (China CDC). Clinical characteristics, laboratory findings, imaging results, and mortality data were collected for a cohort of 70 hospitalized patients with confirmed influenza B from 9 hospitals across China. RESULTS: Starting from October 2021, a substantial increase in the number of patients visiting the designated fever clinics in Beijing was observed, with this trend continuing until January 2022. COVID-19 tests conducted on these patients yielded negative results, while the positivity rate for influenza rose from approximately 8% in October 2021 to over 40% by late January 2022. The cases started to decline after this peak. Data from China CDC confirmed that influenza B is a major pathogen during the season. Sequencing of the viral strain revealed the presence of the Victoria-like lineage of the influenza B strain, with minor variations from the Florida/39/2018 strain. Analysis of the hospitalized patients' characteristics indicated that severe cases were relatively more prevalent among younger individuals, with an average age of 40.9 ± 24.1 years. Among the seven patients who succumbed to influenza, the average age was 30 ± 30.1 years. These patients exhibited secondary infections involving either bacterial or fungal pathogens and displayed elevated levels of cell death markers (such as LDH) and coagulation pathway markers (D-dimer). CONCLUSION: Influenza B represents a significant infection threat and can lead to substantial morbidity and mortality, particularly among young patients. To mitigate morbidity and mortality rates, it is imperative to implement appropriate vaccination and other preventive strategies.

Simultaneous Prediction of Interaction Sites on the Protein and Peptide Sides of Complexes through Multilayer Graph Convolutional Networks.

Identifying the binding residues of protein-peptide complexes is essential for understanding protein function mechanisms and exploring drug discovery. Recently, many computational methods have been developed to predict the interaction sites of either protein or peptide. However, to our knowledge, no prediction method can simultaneously identify the interaction sites on both the protein and peptide sides. Here, we propose a deep graph convolutional network (GCN)-based method called GraphPPepIS to predict the interaction sites of protein-peptide complexes using protein and peptide structural information. We also propose a companion method, SeqPPepIS, for assisting with the lack of structural information and the flexibility of peptides. SepPPepIS replaces the peptide structural features in GraphPPepIS by learning features from peptide sequences. We performed a comprehensive evaluation of the benchmark data sets, and the results show that our two methods outperform state-of-the-art methods on the accurate interaction sites of both protein and peptide sides. We show that our methods can help improve protein-peptide docking. For docking data sets, our methods maintain robust performance in identifying binding sites, thereby enhancing the prediction of peptide binding poses. Finally, we visualized the analysis of protein and peptide graph embedding to demonstrate the learning ability of graph convolution in predicting interaction sites, which was mainly obtained through the shared parameters of a protein graph and peptide graph.

CPS1 maintains pyrimidine pools and DNA synthesis in KRAS/LKB1-mutant lung cancer cells.

Metabolic reprogramming by oncogenic signals promotes cancer initiation and progression. The oncogene KRAS and tumour suppressor STK11, which encodes the kinase LKB1, regulate metabolism and are frequently mutated in non-small-cell lung cancer (NSCLC). Concurrent occurrence of oncogenic KRAS and loss of LKB1 (KL) in cells specifies aggressive oncological behaviour. Here we show that human KL cells and tumours share metabolomic signatures of perturbed nitrogen handling. KL cells express the urea cycle enzyme carbamoyl phosphate synthetase-1 (CPS1), which produces carbamoyl phosphate in the mitochondria from ammonia and bicarbonate, initiating nitrogen disposal. Transcription of CPS1 is suppressed by LKB1 through AMPK, and CPS1 expression correlates inversely with LKB1 in human NSCLC. Silencing CPS1 in KL cells induces cell death and reduces tumour growth. Notably, cell death results from pyrimidine depletion rather than ammonia toxicity, as CPS1 enables an unconventional pathway of nitrogen flow from ammonia into pyrimidines. CPS1 loss reduces the pyrimidine to purine ratio, compromises S-phase progression and induces DNA-polymerase stalling and DNA damage. Exogenous pyrimidines reverse DNA damage and rescue growth. The data indicate that the KL oncological genotype imposes a metabolic vulnerability related to a dependence on a cross-compartmental pathway of pyrimidine metabolism in an aggressive subset of NSCLC.

Development of a multicriteria decision-making model for evaluating hybrid offspring in the sweetpotato (Ipomoea batatas L.) breeding process.

Sweetpotato variety breeding is always a long process. Screening of hybrid offspring is dominated by empirical judgment in this process. Data analysis and decision fatigue have been troubling breeders. In recent years, the low-efficiency screening mode has been unable to meet the requirements of sweetpotato germplasm innovation. Therefore, it is necessary to construct a high-efficiency method that can screen germplasms for different usages, for mining elite genotypes, and to create dedicated sweetpotato varieties. In this article, the multicriteria decision-making (MCDM) model was constructed based on six agronomic traits, including fresh root yield, vine length, vine diameter, branch number, root number and the spatial distribution of storage roots, and five quality traits, including dry matter content, marketable root yield, uniformity of roots, starch content and the edible quality score. Among these, the edible quality score was calculated by using fuzzy comprehensive evaluation to integrate the sensory scores of color, odor, sweetness, stickiness and fibrous taste. The MCDM model was compared with the traditional screening method via an evaluation in 25 sweetpotato materials. The interference of subjective factors on the evaluation results was significantly reduced. The MCDM model is more overall, more accurate and faster than the traditional screening method in the selection of elite sweetpotato materials. It could be programmed to serve the breeders in combination with the traditional screening method.

Sulfur-decorated Fe/C composite synthesized from MIL-88A(Fe) for peroxymonosulfate activation towards tetracycline degradation: Multiple active sites and non-radical pathway dominated mechanism.

Peroxymonosulfate (PMS)-mediated advanced oxidation processes gain growing attention in degrading antibiotics (e.g., tetracycline (TC)) in wastewater for their high capacity and relatively low cost, while designing efficient catalysts for PMS activation remains a challenge. In this study, a sulfur-doped Fe/C catalyst (Fe@C-S) synthesized from iron metal-organic frameworks (Fe-MOFs) was developed for PMS activation towards TC removal. Under optimal conditions, the TC removal efficiency of Fe@C-S150/PMS system within 40 min was 91.2%. Meanwhile, the k value for Fe@C-S150/PMS system (0.2038 min-1) was 3.36-fold as high as the S-free Fe@C-based PMS system. Also, Fe@C-S150/PMS system showed high robustness in different water matrices. Further studies found that the TC degradation mechanism was mainly ascribed to the non-radical pathway (1O2 and electron transfer). Fe nanoparticles, S and CO groups on the catalyst all participated in the generation of reactive oxygen species (ROS). Besides, S species could enhance the Fe2+/Fe3+ redox cycle and accelerate the electron transfer process. This work highlights the critical role of S in enhancing the catalytic performance of Fe/C-based catalysts for PMS activation, which would provide meaningful insights into the design of high-performance PMS activators for the sustainable remediation of emerging contaminants-polluted water bodies.

Evolution law of comet-shaped defects in magnetorheological finishing.

Ultraprecision optical systems such as high-energy laser systems have very strict requirements for surface defects of optical components. Comet-shaped defects are common in magnetorheological finishing, which is one of the main ultraprecision manufacturing technologies. Computational fluid dynamics is used to simulate the pressure and velocity of the magnetorheological fluid in pits during magnetorheological finishing for the first time, to the best of our knowledge. The simulation model analyzes the distribution of velocity and pressure in the process of the formation of comet-shaped defects. After the simulation, magnetorheological polishing experiments were carried out on the fused quartz silica glass. The polishing experiments use the rotation speed and immersion depth of the polishing wheel to verify the velocity and pressure distribution in the simulation experiment. The polishing experiment results are consistent with the simulation results, which show that the high polishing speed will increase the depth of the defect, and the uneven velocity and pressure distribution at the pits make the pits evolve into comet-shaped defects. We explain the cause of comet-shaped defects from the perspective of material removal and provide an optimization method of process parameters for high-efficiency removal of comet-shaped defects.

Wi-Fi Fingerprint-Based Indoor Localization Method via Standard Particle Swarm Optimization.

With the continuous development and improvement in Internet-of-Things (IoT) technology, indoor localization has received considerable attention. Particularly, owing to its unique advantages, the Wi-Fi fingerprint-based indoor-localization method has been widely investigated. However, achieving high-accuracy localization remains a challenge. This study proposes an application of the standard particle swarm optimization algorithm to Wi-Fi fingerprint-based indoor localization, wherein a new two-panel fingerprint homogeneity model is adopted to characterize fingerprint similarity to achieve better performance. In addition, the performance of the localization method is experimentally verified. The proposed localization method outperforms conventional algorithms, with improvements in the localization accuracy of 15.32%, 15.91%, 32.38%, and 36.64%, compared to those of KNN, SVM, LR, and RF, respectively.

Research on the Necessity of Lie Group Strapdown Inertial Integrated Navigation Error Model Based on Euler Angle.

In response to the lack of specific demonstration and analysis of the research on the necessity of the Lie group strapdown inertial integrated navigation error model based on the Euler angle, two common integrated navigation systems, strapdown inertial navigation system/global navigation satellite system (SINS/GNSS) and strapdown inertial navigation system/doppler velocity log (SINS/DVL), are used as subjects, and the piecewise constant system (PWCS) matrix, based on the Lie group error model, is established. From three aspects of variance estimation, the observability and performance of the system with large misalignment angles for low, medium, and high accuracy levels, traditional error model, Lie group left error model, and right error model are compared. The necessity of research on Lie group error model is analyzed quantitatively and qualitatively. The experimental results show that Lie group error model has better stability of variance estimation, estimation accuracy, and observability than traditional error model, as well as higher practical value.

An Improved Online Fast Self-Calibration Method for Dual-Axis RINS Based on Backtracking Scheme.

In the field of high accuracy dual-axis rotational inertial navigation system (RINS), the calibration accuracy of the gyroscopes and accelerometers is of great importance. Although rotation modulation can suppress the navigation error caused by scale factor error and bias error in a static condition, it cannot suppress the scale factor errors thoroughly during the maneuvering process of the vehicle due to the two degrees of rotation freedom. The self-calibration method has been studied by many researchers. However, traditional calibration methods need several hours to converge, which is unable to meet the demand for quick response to positioning and orientation. To solve the above problems, we do the following work in this study: (1) we propose a 39-dimensional online calibration Kalman filtering (KF) model to estimate all calibration parameters; (2) Error relationship between calibration parameters error and navigation error are derived; (3) A backtracking filtering scheme is proposed to shorten the calibration process. Experimental results indicate that the proposed method can shorten the calibration process and improve the calibration accuracy simultaneously.

Physiological and Transcriptome Analysis on Diploid and Polyploid Populus ussuriensis Kom. under Salt Stress.

Populus ussuriensis Kom. is a valuable forest regeneration tree species in the eastern mountainous region of Northeast China. It is known that diploid P. ussuriensis (CK) performed barely satisfactorily under salt stress, but the salt stress tolerance of polyploid (i.e., triploid (T12) and tetraploid (F20)) P. ussuriensis is still unknown. In order to compare the salt stress tolerance and salt stress response mechanism between diploid and polyploid P. ussuriensis, phenotypic observation, biological and biochemistry index detections, and transcriptome sequencing (RNA-seq) were performed on CK, T12, and F20. Phenotypic observation and leaf salt injury index analysis indicated CK suffered more severe salt injury than T12 and F20. SOD and POD activity detections indicated the salt stress response capacity of T12 was stronger than that of CK and F20. MDA content, proline content and relative electric conductivity detections indicated CK suffered the most severe cell-membrane damage, and T12 exhibited the strongest osmoprotective capacity under salt stress. Transcriptome analysis indicated the DEGs of CK, T12, and F20 under salt stress were different in category and change trend, and there were abundant WRKY, NAM, MYB and AP2/ERF genes among the DEGs in CK, T12, and F20 under salt stress. GO term enrichment indicated the basic growth progresses of CK, and F20 was obviously influenced, while T12 immediately launched more salt stress response processes in 36 h after salt stress. KEGG enrichment indicated the DEGs of CK mainly involved in plant−pathogen interaction, ribosome biogenesis in eukaryotes, protein processing in endoplasmic reticulum, degradation of aromatic compounds, plant hormone signal transduction, photosynthesis, and carbon metabolism pathways. The DEGs of T12 were mainly involved in plant−pathogen interaction, cysteine and methionine metabolism, phagosomes, biosynthesis of amino acids, phenylalanine, tyrosine and tryptophan biosynthesis, plant hormone signal transduction, and starch and sucrose metabolism pathways. The DEGs of F20 were mainly involved in plant hormone signal transduction, plant−pathogen interaction, zeatin biosynthesis, and glutathione metabolism pathways. In conclusion, triploid exhibited stronger salt stress tolerance than tetraploid and diploid P. ussuriensis (i.e., T12 > F20 > CK). The differences between the DEGs of CK, T12, and F20 probably are the key clues for discovering the salt stress response signal transduction network in P. Ussuriensis.

Monitoring residue levels and dietary risk assessment of thiamethoxam and its metabolite clothianidin for Chinese consumption of Chinese kale.

BACKGROUND: Thiamethoxam is widely used to control pests in Chinese kale, popularly consumed leafy vegetables. The potential risk to the environment and human health has aroused much public concern. Therefore, it is important to investigate the degradation behavior, residue distribution and dietary risk assessment of thiamethoxam in Chinese kale. RESULTS: A sensitive analytical method for determination of thiamethoxam and its metabolite clothianidin residue in Chinese kale was established and validated through a quick, easy, cheap, effective, rugged, and safe (QuEChERS) technique with ultra-performance liquid chromatography tandem mass spectrometry (UPLC-MS/MS). The recoveries were 85.4-101.2% for thiamethoxam and 79.5-108.1% for clothianidin, with the relative standard deviations (RSDs) of 0.9-10.2% and 1.8-6.0%, respectively. For the dissipation kinetics, the data showed that thiamethoxam in Chinese kale was degraded with the half-lives of 4.1 to 4.5 days. In the terminal residue experiments, the residues of thiamethoxam were 0.017-0.357 mg kg-1 after application 2-3 times with a preharvest interval (PHI) of 7 days under the designed dosages. The chronic and acute dietary exposure assessment risk quotient (RQ) values of thiamethoxam in Chinese kale for different Chinese consumers were 0.08-0.19% and 0.05-0.12%, respectively, and those of clothianidin were 0.01-0.04% and 0.02-0.04%, respectively, all of the RQ values were lower than 100%. CONCLUSION: Thiamethoxam in Chinese kale was rapidly degraded following first-order kinetics models. The dietary risk of thiamethoxam and clothianidin through Chinese kale was negligible to consumers. The results from this study are important reference for Chinese governments to developing criteria for the safe and rational use of thiamethoxam, setting maximum residue levels (MRLs), monitoring the quality safety of agricultural products and protecting consumer health. © 2021 Society of Chemical Industry.

Residue analysis and dietary risk assessment of tebuconazole in loquat and sugarcane after open-field application in China.

Tebuconazole is an effective but potentially carcinogenic fungicide. Tebuconazole residues in loquat and sugarcane were investigated, together with an analysis of their dietary risk. Solid-phase extraction and high-performance liquid chromatography-tandem mass spectrometry showed fortified recoveries of 87% to 98%, with standard deviations ≤ 6.4%. The tebuconazole residues were below 6.34 mg kg-1 for whole loquat and 0.80 mg kg-1 for sugarcane. Based on the final residue test, the chronic hazard index was 90.74% in whole loquat and sugarcane, while the acute hazard index ranged from 66.32% to 109.49% for loquat and 33.37% to 76.93% for sugarcane. The assessment of potential dietary risk showed that while loquats may represent an acute risk to human health, whereas sugarcane consumption was safe. A 21-day pre-harvest interval is proposed to comply with maximum residue limit regulations and the observed risks when using 430 g/L tebuconazole for treating sugarcane. These findings offer guidance for tebuconazole formulation and use.

The Role of NLRP3 Inflammasome in IgA Nephropathy.

Immunoglobulin A nephropathy (IgAN) is the most common primary glomerular disease worldwide today. The NLRP3 inflammasome is a polyprotein complex and an important participant in inflammation. Accumulating studies have shown that the NLRP3 inflammasome participates in a variety of kidney diseases, including IgAN. This review focuses on the role of the NLRP3 inflammasome in IgAN and summarizes multiple involved pathways, which may provide novel treatments for IgAN treatment.

A genome wide transcriptional study of Populus alba x P. tremula var. glandulosa in response to nitrogen deficiency stress.

Poplar 84 K (Populus alba x P. tremula var. glandulosa) is a good resource for genetic engineering due to its rapid growth and wide adaptability, and it is also an excellent ornamental tree species. In this study, we used 84 K plantlets grown in the nitrogen-limited medium as experimental materials to explore the molecular mechanism in 84 K leaves under nitrogen deficiency. A total of 5,868 differentially expressed genes (DEGs) were identified using the transcriptional information from RNA-seq data. GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment results revealed that the DEGs were mainly involved in energy metabolism and anthocyanin biosynthesis. We then identified differentially expressed transcription factors (TFs) and constructed TF centered gene co-expression networks for chlorophyll and anthocyanin biosynthesis pathway genes. Twenty potential regulators were finally identified. We speculated the transcription factors that control the pigmentation in leaves with the MYB-bHLH-WD40 (MBW) pigment regulatory model. Such identification will clarify the genetic basis of the secondary metabolism in 84 K, and being a source of candidate genes for future plant genetic engineering. Our work broadens the researchers' understanding of the regulation of anthocyanin synthesis in trees and provides new perspectives for ornamental 84 K poplar breeding. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s12298-021-01012-3.

Cumulative risk assessment of phthalates in edible vegetable oil consumed by Chinese residents.

BACKGROUND: Phthalates have been widely used as plasticizers in various industries and are widely focused on in the international community as a result of their reproductive toxicity. Exposure of Chinese residents to phthalates via edible vegetable oil occurs often. In the present study, gas chromatography-mass spectrometry was used to detect the two main phthalates bis(2-ethylhexyl) phthalate (DEHP) and dibutyl phthalate (DBP) in four major edible vegetable oil sources: an edible oil blend, soybean oil, peanut oil and rapeseed oil (a total of 1016 samples), as collected throughout China. Furthermore, cumulative risk assessment was used to estimate the reproductive health risk to Chinese residents caused by the phthalates that come from edible vegetable oils. RESULTS: Both phthalates were detected in four major edible vegetable oil sources. The phthalate with the highest detection rate was DBP (13.48%), followed by DEHP (7.78%). The results of the cumulative risk assessment showed that the hazard indices of these two phthalates in edible vegetable oils were less than 1, except in soybean oil. Nevertheless, the two phthalates had the lowest detection rates in soybean oil, which were 1.94% (DEHP) and 5.16% (DBP). In China, contamination levels of phthalates in the soils where oil crops are cultivated have a great influence on the phthalate concentrations in edible vegetable oils. CONCLUSION: It is recommended that Chinese residents who are consuming soybean oil choose well-known brands and regularly change their brand of consumption. The phthalates in edible vegetable oils pose a relatively small reproductive health risk to Chinese residents. © 2019 Society of Chemical Industry.