FastMM: an efficient toolbox for personalized constraint-based metabolic modeling.

BACKGROUND: Constraint-based metabolic modeling has been applied to understand metabolism related disease mechanisms, to predict potential new drug targets and anti-metabolites, and to identify biomarkers of complex diseases. Although the state-of-art modeling toolbox, COBRA 3.0, is powerful, it requires substantial computing time conducting flux balance analysis, knockout analysis, and Markov Chain Monte Carlo (MCMC) sampling, which may limit its application in large scale genome-wide analysis. RESULTS: Here, we rewrote the underlying code of COBRA 3.0 using C/C++, and developed a toolbox, termed FastMM, to effectively conduct constraint-based metabolic modeling. The results showed that FastMM is 2~400 times faster than COBRA 3.0 in performing flux balance analysis and knockout analysis and returns consistent outputs. When applied to MCMC sampling, FastMM is 8 times faster than COBRA 3.0. FastMM is also faster than some efficient metabolic modeling applications, such as Cobrapy and Fast-SL. In addition, we developed a Matlab/Octave interface for fast metabolic modeling. This interface was fully compatible with COBRA 3.0, enabling users to easily perform complex applications for metabolic modeling. For example, users who do not have deep constraint-based metabolic model knowledge can just type one command in Matlab/Octave to perform personalized metabolic modeling. Users can also use the advance and multiple threading parameters for complex metabolic modeling. Thus, we provided an efficient and user-friendly solution to perform large scale genome-wide metabolic modeling. For example, FastMM can be applied to the modeling of individual cancer metabolic profiles of hundreds to thousands of samples in the Cancer Genome Atlas (TCGA). CONCLUSION: FastMM is an efficient and user-friendly toolbox for large-scale personalized constraint-based metabolic modeling. It can serve as a complementary and invaluable improvement to the existing functionalities in COBRA 3.0. FastMM is under GPL license and can be freely available at GitHub site: https://github.com/GonghuaLi/FastMM.

Correction to: FastMM: an efficient toolbox for personalized constraint-based metabolic modeling.

An amendment to this paper has been published and can be accessed via the original article.

Dietary antioxidant vitamin C influences the evolutionary path of insecticide resistance in Drosophila melanogaster.

Herbivorous insects encounter a variety of toxic environmental substances ranging from ingested plant defensive compounds to human-introduced insecticidal agents. Dietary antioxidants are known to reduce the negative physiological impacts of toxins in mammalian systems through amelioration of reactive oxygen-related cellular damage. The analogous impacts to insects caused by multigenerational exposure to pesticides and the effects on adaptive responses within insect populations, however, are currently unknown. To address these research gaps, we used Drosophila as a model system to explore adaptive phenotypic responses to acute dichlorodiphenyltrichloroethane (DDT) exposure in the presence of the dietary antioxidant vitamin C and to examine the structural genomic consequences of this exposure. DDT resistance increased significantly among four replicates exposed to a low concentration of DDT for 10 generations. In contrast, dietary intake of vitamin C significantly reduced DDT resistance after mutigenerational exposure to the same concentration of DDT. As to the genomic consequences, no significant differences were predicted in overall nucleotide substitution rates across the genome between any of the treatments. Despite this, replicates exposed to a low concentration of DDT without vitamin C showed the highest number of synonymous and non-synonymous variants (3196 in total), followed by the DDT plus vitamin C (1174 in total), and vitamin C alone (728 in total) treatments. This study demonstrates the potential role of diet (specifically, antioxidant intake) on adaptive genome responses, and thus on the evolution of pesticide resistance within insect populations.

Gain and loss events in the evolution of the apolipoprotein family in vertebrata.

BACKGROUND: Various apolipoproteins widely distributed among vertebrata play key roles in lipid metabolism and have a direct correlation with human diseases as diagnostic markers. However, the evolutionary progress of apolipoproteins in species remains unclear. Nine human apolipoproteins and well-annotated genome data of 30 species were used to identify 210 apolipoprotein family members distributed among species from fish to humans. Our study focused on the evolution of nine exchangeable apolipoproteins (ApoA-I/II/IV/V, ApoC-I~IV and ApoE) from Chondrichthyes, Holostei, Teleostei, Amphibia, Sauria (including Aves), Prototheria, Marsupialia and Eutheria. RESULTS: In this study, we reported the overall distribution and the frequent gain and loss evolutionary events of apolipoprotein family members in vertebrata. Phylogenetic trees of orthologous apolipoproteins indicated evident divergence between species evolution and apolipoprotein phylogeny. Successive gain and loss events were found by evaluating the presence and absence of apolipoproteins in the context of species evolution. For example, only ApoA-I and ApoA-IV occurred in cartilaginous fish as ancient apolipoproteins. ApoA-II, ApoE, and ApoC-I/ApoC-II were found in Holostei, Coelacanthiformes, and Teleostei, respectively, but the latter three apolipoproteins were absent from Aves. ApoC-I was also absent from Cetartiodactyla. The apolipoprotein ApoC-III emerged in terrestrial animals, and ApoC-IV first arose in Eutheria. The results indicate that the order of the emergence of apolipoproteins is most likely ApoA-I/ApoA-IV, ApoE, ApoA-II, ApoC-I/ApoC-II, ApoA-V, ApoC-III, and ApoC-IV. CONCLUSIONS: This study reveals not only the phylogeny of apolipoprotein family members in species from Chondrichthyes to Eutheria but also the occurrence and origin of new apolipoproteins. The broad perspective of gain and loss events and the evolutionary scenario of apolipoproteins across vertebrata provide a significant reference for the research of apolipoprotein function and related diseases.

Regulatory network reconstruction of five essential microRNAs for survival analysis in breast cancer by integrating miRNA and mRNA expression datasets.

Although many of the genetic loci associated with breast cancer risk have been reported, there is a lack of systematic analysis of regulatory networks composed of different miRNAs and mRNAs on survival analysis in breast cancer. To reconstruct the microRNAs-genes regulatory network in breast cancer, we employed the expression data from The Cancer Genome Atlas (TCGA) related to five essential miRNAs including miR-21, miR-22, miR-210, miR-221, and miR-222, and their associated functional genomics data from the GEO database. Then, we performed an integration analysis to identify the essential target factors and interactions for the next survival analysis in breast cancer. Based on the results of our integrated analysis, we have identified significant common regulatory signatures including differentially expressed genes, enriched pathways, and transcriptional regulation such as interferon regulatory factors (IRFs) and signal transducer and activator of transcription 1 (STAT1). Finally, a reconstructed regulatory network of five miRNAs and 34 target factors was established and then applied to survival analysis in breast cancer. When we used expression data for individual miRNAs, only miR-21 and miR-22 were significantly associated with a survival change. However, we identified 45 significant miRNA-gene pairs that predict overall survival in breast cancer out of 170 one-on-one interactions in our reconstructed network covering all of five miRNAs, and several essential factors such as PSMB9, HLA-C, RARRES3, UBE2L6, and NMI. In our study, we reconstructed regulatory network of five essential microRNAs for survival analysis in breast cancer by integrating miRNA and mRNA expression datasets. These results may provide new insights into regulatory network-based precision medicine for breast cancer.

Analysis of false-positive results of diatom test in the diagnosis of drowning-would not be an impediment.

Diatom test has been a significant tool for the diagnosis of drowning. The reliability of the diatom test is still in strong dispute in the field of forensic science because of the false-positive results. This study was designed to quantitatively compare the numbers of the diatoms in false-positive cases and true drowning cases. Diatom samples from 64 victims were used in this study: 32 cases are confirmed drowning victims and other 32 cases died from non-drowned death. Samples were subject for the diatom test those were analyzed by the microwave digestion-vacuum filtration-automated scanning electron microscopy method (MD-VF-Auto SEM method) that we developed before. The results did show that there are false-positive diatoms detected in the liver and kidney tissues of non-drowned bodies: 6/20 in liver tissues and 7/20 in kidney tissues. However, the quantitative studies showed that there are statistical differences with the numbers of diatoms in the false-positive cases and in the true drowning cases. Diatom test of single organ is difficult for us to distinguish the sources of the diatoms detected. Therefore, comprehensive analysis of multiple organs would be more useful for the diagnosis of drowning.

Molecular cloning, expression profiling, and functional analysis of a broad-complex isoform 2/3 (Br-Z2/Z3) transcription factor in the diamondback moth, Plutella xylostella (L.).

The diamondback moth, Plutella xylostella (L.) (Lepidoptera: Plutellidae), is a widespread and destructive pest of cruciferous crops. New strategies for controlling it are needed because it is rapidly developing resistance to conventional pesticides. In insects, transcription factors (TFs) including broad-complex (Br-C) are thought to be useful for insecticide development because they are able to regulate the transcription of functional genes involved in responses to external stimuli including insecticides. In the present study, we cloned and sequenced the open reading frames (ORFs) of three BTB-ZF encoding genes from the diamondback moth deposited in the National Center for Biotechnology Information (NCBI) database under accessions MG753773, MG288674, and MG753772. The lengths of these ORFs were 1,680, 1,428, and 1,647 bp, respectively. The phylogenetic analysis based on the predicted amino acid sequences of ZF domains showed that MG753773 and MG288674 belonged to Z2/Z3 and Z7 of Br-C while MG753772 belonged to Ttk types. In the agreement, the highest expression level of MG753773 occurred during the prepupal stage, MG288674 and MG753772 were expressed during all stages and peaked in the adult and egg stages, respectively. RNA interference silencing of MG753773 in the late third instar larvae significantly decreased survival and pupation of the insects. With precocene II, transcription of MG753773 increased (4×) in the fourth instar larva 24 hr later; 48 hr later the rate of prepupation and pupation was significantly higher. These findings will contribute to the development of new regulators of the growth and development for diamondback moth control.

Folate Deficiency and Gene Polymorphisms of MTHFR, MTR and MTRR Elevate the Hyperhomocysteinemia Risk.

BACKGROUND: Hyperhomocysteinemia (HHcy) is an independent risk factor for cardiovascular diseases (CVDs). We aimed to investigate the joint effect of homocysteine metabolism gene polymorphisms, as well as the folate deficiency on the risk of HHcy in a Chinese hypertensive population. METHODS: This study enrolled 480 hypertensive patients aged 28 - 75 from six hospitals in different Chinese regions from 9/2005 - 12/2005. Known genotypes of methylenetetrahydrofolate reductase (MTHFR) C677T and A1298C, methionine synthase (MTR) A2756G, and methionine synthase reductase (MTRR) A66G were detected by PCRRFLP methods. Serum Hcy was measured by high-performance liquid chromatography and serum folate was measured by chemiluminescent immunoassay. RESULTS: MTHFR C677T and MTR A2756G can independently elevate the risk of HHcy (TT vs. CC + CT, p < 0.001 and AG + GG vs. AA, p = 0.026, respectively), whereas MTHFR A1298C decreased HHcy risk (AC + CC vs. AA, p < 0.001) and showed a protective effect against HHcy risk. Importantly, the joint effect of these risk genotypes showed significantly higher odds of HHcy than non-risk genotypes, especially the patients with four risk genotypes. It is noteworthy that this deleterious effect was aggravated by folate deficiency. These findings were verified by generalized multifactor dimensionality reduction model (p = 0.001) and a cumulative effects model (p < 0.001). CONCLUSIONS: We have first demonstrated that the joint effect of homocysteine metabolism gene polymorphisms and folate deficiency lead to dramatic elevations in the HHcy risk.

Inferring therapeutic targets from heterogeneous data: HKDC1 is a novel potential therapeutic target for cancer.

MOTIVATION: The discovery of therapeutic targets is important for cancer treatment. Although dozens of targets have been used in cancer therapies, cancer remains a serious disease with a high mortality rate. Owing to the expansion of cancer-related data, we now have the opportunity to infer therapeutic targets using computational biology methods. RESULTS: Here, we describe a method, termed anticancer activity enrichment analysis, used to determine genes that could be used as therapeutic targets. The results show that these genes have high likelihoods of being developed into clinical targets (>60%). Combined with gene expression data, we predicted 50 candidate targets for lung cancer, of which 19 of the top 20 genes are targeted by approved drugs or drugs used in clinical trials. A hexokinase family member, hexokinase domain-containing protein 1 (HKDC1), is the only one of the top 20 genes that has not been targeted by either an approved drug or one being used in clinical trials. Further investigations indicate that HKDC1 is a novel potential therapeutic target for lung cancer. CONCLUSION: We developed a protocol to identify potential therapeutic targets from heterogeneous data. We suggest that HKDC1 is a novel potential therapeutic target for lung cancer. CONTACT: huangjf@mail.kiz.ac.cn SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A COMPREHENSIVE STUDY ON APOPTOSIS INDUCTION BY AZADIRACHTIN IN Spodoptera frugiperda CULTURED CELL LINE Sf9.

The induction of apoptosis by azadirachtin, a well-known botanical tetranortriterpenoid isolated from the neem tree (Azadirachta indica A. Juss) and other members of the Meliaceae, was investigated in Spodoptera frugiperda cultured cell line (Sf9). Morphological changes in Sf9 cells treated by various concentrations of azadirachtin were observed at different times under light microscopy. Morphological and biochemical analysis indicated that Sf9 cells treated by 1.5 µg/mL azadirachtin showed typical morphological changes, which were indicative of apoptosis and a clear DNA ladder. The flow cytometry analysis showed the apoptosis rate reached a maximum value of 32.66% at 24 h with 1.5 µg/mL azadirachtin in Sf9 cells. The inhibition of Sf9 cell proliferation suggested that the effect of azadirachtin was dose dependent and the EC50 at 48 and 72 h was 2.727 × 10(-6) and 6.348 × 10(-9) µg/mL, respectively. The treatment of azadirachtin in Sf9 cells could significantly increase the activity of Sf caspase-1, but showed no effect on the activity of Topo I, suggesting that the apoptosis induced by azadirachtinin Sf9 cells is through caspase-dependent pathway. These results provided not only a series of morphological, biochemical, and toxicological comprehensive evidences for induction of apoptosis by azadirachtin, but also a reference model for screening insect cell apoptosis inducers from natural compounds.

Joint associations of folate, homocysteine and MTHFR, MTR and MTRR gene polymorphisms with dyslipidemia in a Chinese hypertensive population: a cross-sectional study.

BACKGROUND: Dyslipidemia is a well-established risk factor for cardiovascular disease. Serum lipids were affected by several gene polymorphisms, folate, homocysteine and other metabolite levels. We aim to investigate the effects of homocysteine metabolism enzyme polymorphisms (MTHTR C677T, MTHFR A1298C, MTR A2756G and MTRR A66G) and their interactions with folate, homocysteine on serum lipid levels in Chinese patients with hypertension. METHODS: Participants were 480 hypertensive adults that enrolled in September to December 2005 from six different Chinese hospitals (Harbin, Shanghai, Shenyang, Beijing, Xi'an, and Nanjing). Known MTHFR C677T, MTHFR A1298C, MTR A2756G and MTRR A66G genotypes were determined by PCR-RFLP methods. Serum folate was measured by chemiluminescent immunoassay, homocysteine was measured by high-performance liquid chromatography, serum lipids parameters were determined by an automatic biochemistry analyzer, low-density lipoprotein was calculated by Friedewald's equation. Unitary linear regression model was used to assess the associations of gene polymorphisms, folate and homocysteine on serum lipid profiles. Unconditional logistic regression model was applied to test the interactions of folate, homocysteine and gene polymorphisms on dyslipidemia. RESULTS: No correlations between gene polymorphisms and homocysteine on serum lipid profiles. Compared with normal folate patients, patients with low folate showed higher odds of hypertriglyceridemia (OR = 2.02, 95 % CI: 1.25-3.25, P = 0.004) and low levels of high-density lipoprotein cholesterol (OR = 1.88, 95 % CI: 1.07-3.28, P = 0.027). Each of four gene polymorphisms (MTHTR C677T, MTHFR A1298C, MTR A2756G and MTRR A66G) combined with low folate showed higher odds of hypertriglyceridemia (P for trend: 0.049, 0.004, 0.007 and 0.005, respectively). MTHFR C677T and A1298C with low folate showed higher odds of low levels of high-density lipoprotein cholesterol (P for trend: 0.008 and 0.031). CONCLUSIONS: Low folate status and homocysteine metabolism gene polymorphisms (MTHTR C677T, MTHFR A1298C, MTR A2756G and MTRR A66G) may have a synergistic effect increased the incidence of dyslipidemia in Chinese hypertensive population.

The characterization and comparison of amyloidogenic segments and non-amyloidogenic segments shed light on amyloid formation.

Amyloid fibrillar aggregates of proteins or peptides are involved in the etiology of several neurodegenerative diseases and represent a major problem in healthcare. Short regions in the protein trigger this aggregation. It is important to understand the basis of such short regions aggregation and amyloidosis for therapeutic intervention. In this study, we describe specific physico-chemical properties of amyloidogenic segments and compare them with non-amyloidogenic segments. First, amyloidogenic segments are characterized by lower values for average net charge, electrostatic potential, solvent accessible surface area and B-factor when compared to the non-amyloidogenic segments of the same proteins. Second, they are enriched in hydrophobic residues and have a tendency to form hydrogen bonds. Thus, amyloidogenic segments have distinct physico-chemical properties that are different from those of non-amyloidogenic segments. Third, and quite unexpectedly, our dynamic simulation studies support the hypothesis that amyloidogenic segments have lower average flexibility than non-amyloidogenic segments. Furthermore, the presence of amyloidogenic segments in disordered proteins does not contradict the observation that amyloidogenic segments are less flexible.

Repeated functional convergent effects of NaV1.7 on acid insensitivity in hibernating mammals.

Hibernating mammals need to be insensitive to acid in order to cope with conditions of high CO2; however, the molecular basis of acid tolerance remains largely unknown. The African naked mole-rat (Heterocephalus glaber) and hibernating mammals share similar environments and physiological features. In the naked mole-rat, acid insensitivity has been shown to be conferred by the functional motif of the sodium ion channel NaV1.7. There is now an opportunity to evaluate acid insensitivity in other taxa. In this study, we tested for functional convergence of NaV1.7 in 71 species of mammals, including 22 species that hibernate. Our analyses revealed a functional convergence of amino acid sequences, which occurred at least six times independently in mammals that hibernate. Evolutionary analyses determined that the convergence results from both parallel and divergent evolution of residues in the functional motif. Our findings not only identify the functional molecules responsible for acid insensitivity in hibernating mammals, but also open new avenues to elucidate the molecular underpinnings of acid insensitivity in mammals.

Bio-degraded of sulfamethoxazole by microbial consortia without addition nutrients: Mineralization, nitrogen removal, and proteomic characterization.

Sulfamethoxazole (SMX) is widely employed as an antibiotic, while its residue in environment has become a common public concern. Using 100 mg/L SMX as the sole nutrient source, the acclimated sludge obtained by this study displayed an excellent SMX degradation performance. The addition of SMX resulted in significant microbiological differentiation within the acclimated sludge. Microbacterium (6.6%) was identified as the relatively dominant genera in metabolism group that used SMX as sole carbon source. Highly expressed proteins from this strain strongly suggested its essential role in SMX degradation, while the degradation of SMX by other strains (Thaurea 78%) in co-metabolism group appeared to also rely on this strain. The interactions of differentially expressed proteins were primarily involved in metabolic pathways including TCA cycle and nitrogen metabolism. It is concluded that the sulfonamides might serve not only as the carbon source but also as the nitrogen source in the reactor. A total of 24 intermediates were identified, 13 intermediates were newly reported. The constructed pathway suggested the mineralizing and nitrogen conversion ability towards SMX. Batch experiments also proved that the acclimated sludge displayed ability to biodegrade other sulfonamides, including SM2 and SDZ and SMX-N could be removed completely.

Molecular dynamics simulation of PNPLA3 I148M polymorphism reveals reduced substrate access to the catalytic cavity.

A missense mutation I148M in PNPLA3 (patatin-like phospholipase domain-containing 3 protein) is significantly correlated with nonalcoholic fatty liver disease (NAFLD). To glean insights into mutation's effect on enzymatic activity, we performed molecular dynamics simulation and flexible docking studies. Our data show that the size of the substrate-access entry site is significantly reduced in mutants, which limits the access of palmitic acid to the catalytic dyad. Besides, the binding free energy calculations suggest low affinity for substrate to mutant enzyme. The substrate-bound system simulations reveal that the spatial arrangement of palmitic acid is distinct in wild-type from that in mutant. The substrate recognition specificity is lost due to the loop where the I148M mutation was located. Our results provide strong evidence for the mechanism by which I148M affects the enzyme activity and suggest that mediating the dynamics may offer a potential avenue for NAFLD.

Hearing aid for vertebrates via multiple episodic adaptive events on prestin genes.

Auditory detection is essential for survival and reproduction of vertebrates, yet the genetic changes underlying the evolution and diversity of hearing are poorly documented. Recent discoveries concerning prestin, which is responsible for cochlear amplification by electromotility, provide an opportunity to redress this situation. We identify prestin genes from the genomes of 14 vertebrates, including three fishes, one amphibian, one lizard, one bird, and eight mammals. An evolutionary analysis of these sequences and 34 previously known prestin genes reveals for the first time that this hearing gene was under positive selection in the most recent common ancestor (MRCA) of tetrapods. This discovery might document the genetic basis of enhanced high sound sensibility in tetrapods. An investigation of the adaptive gain and evolution of electromotility, an important evolutionary innovation for the highest hearing ability of mammals, detects evidence for positive selections on the MRCA of mammals, therians, and placentals, respectively. It is suggested that electromotility determined by prestin might initially appear in the MRCA of mammals, and its functional improvements might occur in the MRCA of therian and placental mammals. Our patch clamp experiments further support this hypothesis, revealing the functional divergence of voltage-dependent nonlinear capacitance of prestin from platypus, opossum, and gerbil. Moreover, structure-based cdocking analyses detect positively selected amino acids in the MRCA of placental mammals that are key residues in sulfate anion transport. This study provides new insights into the adaptation and functional diversity of hearing sensitivity in vertebrates by evolutionary and functional analysis of the hearing gene prestin.

CDRUG: a web server for predicting anticancer activity of chemical compounds.

Cancer is the leading cause of death worldwide. Screening anticancer candidates from tens of millions of chemical compounds is expensive and time-consuming. A rapid and user-friendly web server, known as CDRUG, is described here to predict the anticancer activity of chemical compounds. In CDRUG, a hybrid score was developed to measure the similarity of different compounds. The performance analysis shows that CDRUG has the area under curve of 0.878, indicating that CDRUG is effective to distinguish active and inactive compounds.

Dimorphism of Candida tropicalis and its effect on nitrogen and phosphorus removal and sludge settleability.

Candida tropicalis PNY, a novel dimorphic strain with the capacity of simultaneous carbon, nitrogen and phosphorus removal in anaerobic and aerobic conditions, was isolated from activated sludge. Dimorphism of C. tropicalis PNY had effect on removing nitrogen and phosphorous and slightly affected COD removal under aerobic condition. Sample with high hypha formation rate (40 ± 5%) had more removal efficiencies of NH4+-N (50 mg/L) and PO43--P (10 mg/L), which could achieve 82.19% and 97.53%, respectively. High hypha cells dosage exhibited good settleability and filamentous overgrowth was not observed. According to label-free quantitative proteomics assays. Up-regulated proteins involved in the mitogen-activated protein kinase (MAPK) pathway indicated the active growth and metabolism process of sample with high hypha formation rate (40 ± 5%). And proteins concerning about glutamate synthetase and SPX domain-contain protein explain for the nutrient removal mechanism including assimilation of ammonia and polyphosphates synthesis.

Efficient degradation of hexabromocyclododecane using montmorillonite supported nano-zero-valent iron and Citrobacter sp. Y3.

The coupling of modified nanoscale zero-valent iron (nZVI) with organohalide-degrading bacteria provides a promising solution for the remediation of hexabromocyclododecane (HBCD)-contaminated environments. However, the interactions between modified nZVI and dehalogenase bacteria are intricate, and the mechanisms of synergistic action and electron transfer are not clear, and requires further specific investigation. In this study, HBCD was used as a model pollutant, and stable isotope analysis revealed that organic montmorillonite (OMt)-supported nZVI coupled with the degrading bacterial strain Citrobacter sp. Y3 (nZVI/OMt-Y3) can use [13C]HBCD as the sole carbon source and degrade or even mineralise it into 13CO2 with a maximum conversion rate of 100% within approximately 5 days. Analysis of the intermediates showed that the degradation of HBCD mainly involves three different pathways: dehydrobromination, hydroxylation, and debromination. The proteomics results showed that nZVI introduction promoted the transport of electrons and debromination. Combining the results from XPS, FTIR, and Raman spectroscopy with the analysis results of proteinomics and biodegradation products, we verified the process of electron transport and proposed a metabolic mechanism of HBCD degradation by the nZVI/OMt-Y3. Moreover, this study provides insightful avenues and models for the further remediation of HBCD and other similar pollutants in the environment.

Exploring the detailed spatiotemporal characteristics of PM2.5: Generating a full-coverage and hourly PM2.5 dataset in the Sichuan Basin, China.

Fine particulate matter (PM2.5) has received worldwide attention due to its threat to public health. In the Sichuan Basin (SCB), PM2.5 is causing heavy health burdens due to its high concentrations and population density. Compared with other heavily polluted areas, less effort has been made to generate a full-coverage PM2.5 dataset of the SCB, in which the detailed PM2.5 spatiotemporal characteristics remain unclear. Considering commonly existing spatiotemporal autocorrelations, the top-of-atmosphere reflectance (TOAR) with a high coverage rate and other auxiliary data were employed to build commonly used random forest (RF) models to generate accurate hourly PM2.5 concentration predictions with a 0.05° × 0.05° spatial resolution in the SCB in 2016. Specifically, with historical concentrations predicted from a spatial RF (S-RF) and observed at stations, an alternative spatiotemporal RF (AST-RF) and spatiotemporal RF (ST-RF) were built in grids with stations (type 1). The predictions from the AST-RF in grids without stations (type 2) and observations in type 1 formed the PM2.5 dataset. The LOOCV R2, RMSE and MAE were 0.94/0.94, 8.71/8.62 µg∕m3 and 5.58/5.57 µg∕m3 in the AST-RF/ST-RF, respectively. Using the produced dataset, spatiotemporal analysis was conducted for a detailed understanding of the spatiotemporal characteristics of PM2.5 in the SCB. The PM2.5 concentrations gradually increased from the edge to the center of the SCB in spatial distribution. Two high-concentration areas centered on Chengdu and Zigong were observed throughout the year, while another high-concentration area centered on Dazhou was only observed in winter. The diurnal variation had double peaks and double valleys in the SCB. The concentrations were high at night and low in daytime, which suggests that characterizing the relationship between PM2.5 and adverse health outcomes by daily means might be inaccurate with most human activities conducted in daytime.