Climate-related naturally occurring epimutation and their roles in plant adaptation in A. thaliana.

DNA methylation has been proposed to be an important mechanism that allows plants to respond to their environments sometimes entirely uncoupled from genetic variation. To understand the genetic basis, biological functions and climatic relationships of DNA methylation at a population scale in Arabidopsis thaliana, we performed a genome-wide association analysis with high-quality single nucleotide polymorphisms (SNPs), and found that ~56% on average, especially in the CHH sequence context (71%), of the differentially methylated regions (DMRs) are not tagged by SNPs. Among them, a total of 3235 DMRs are significantly associated with gene expressions and potentially heritable. 655 of the 3235 DMRs are associated with climatic variables, and we experimentally verified one of them, HEI10 (HUMAN ENHANCER OF CELL INVASION NO.10). Such epigenetic loci could be subjected to natural selection thereby affecting plant adaptation, and would be expected to be an indicator of accessions at risk. We therefore incorporated these climate-related DMRs into a gradient forest model, and found that the natural A. thaliana accessions in Southern Europe that may be most at risk under future climate change. Our findings highlight the importance of integrating DNA methylation that is independent of genetic variations, and climatic data to predict plants' vulnerability to future climate change.

Regulation of capsule spine formation in castor.

Castor (Ricinus communis L.) is a dicotyledonous oilseed crop that can have either spineless or spiny capsules. Spines are protuberant structures that differ from thorns or prickles. The developmental regulatory mechanisms governing spine formation in castor or other plants have remained largely unknown. Herein, using map-based cloning in 2 independent F2 populations, F2-LYY5/DL01 and F2-LYY9/DL01, we identified the RcMYB106 (myb domain protein 106) transcription factor as a key regulator of capsule spine development in castor. Haplotype analyses demonstrated that either a 4,353-bp deletion in the promoter or a single nucleotide polymorphism leading to a premature stop codon in the RcMYB106 gene could cause the spineless capsule phenotype in castor. Results of our experiments indicated that RcMYB106 might target the downstream gene RcWIN1 (WAX INDUCER1), which encodes an ethylene response factor known to be involved in trichome formation in Arabidopsis (Arabidopsis thaliana) to control capsule spine development in castor. This hypothesis, however, remains to be further tested. Nevertheless, our study reveals a potential molecular regulatory mechanism underlying the spine capsule trait in a nonmodel plant species.

The double flower variant of yellowhorn is due to a LINE1 transposon-mediated insertion.

As essential organs of reproduction in angiosperms, flowers, and the genetic mechanisms of their development have been well characterized in many plant species but not in the woody tree yellowhorn (Xanthoceras sorbifolium). Here, we focused on the double flower phenotype in yellowhorn, which has high ornamental value. We found a candidate C-class gene, AGAMOUS1 (XsAG1), through bovine serum albumin sequencing and genetics analysis with a Long Interpersed Nuclear Elements 1 (LINE1) transposable element fragment (Xsag1-LINE1-1) inserted into its second intron that caused a loss-of-C-function and therefore the double flower phenotype. In situ hybridization of XsAG1 and analysis of the expression levels of other ABC genes were used to identify differences between single- and double-flower development processes. These findings enrich our understanding of double flower formation in yellowhorn and provide evidence that transposon insertions into genes can reshape plant traits in forest trees.

Semirational Design Strategy To Enhance the Thermostability and Catalytic Activity of Cytochrome P450 105D7 for the Degradation of the Pharmaceutically Active Compounds: Diclofenac.

Pharmaceutically active compounds are an important category of emerging pollutants, and their biological transformation processes in the environment are crucial for understanding and evaluating the migration, transformation, and environmental fate of emerging pollutants. The cytochrome P450 105 enzyme family has been proven to play an important role in the degradation of exogenous environmental pollutants. However, its thermostability and catalytic activity still need to be improved to better adapt to complex environmental conditions. This work elucidates the key mechanisms and important residues of the degradation reaction through multiple computational strategies, establishes a mutation library, and obtains 21 single-point mutation designs. Experimental verification showed that 16 single mutants had enhanced thermostability, with the R89F and L197Y mutants showing the highest increases in thermostability at 135 and 119% relative to the wild-type enzyme, respectively. Additionally, as a result of the higher specific activity of D390Q, it was selected for combination mutagenesis, ultimately resulting in three combination mutants (R89F/L197Y, R89F/D390Q, and R89F/L197Y/D390Q) with enhanced thermostability and catalytic activity. This study provides a modification approach for constructing efficient enzyme variants through semirational design and can contribute to the development of control technologies for emerging pollutants.

Pathway conversion enables a double-lock mechanism to maintain DNA methylation and genome stability.

The CMT2 and RNA-directed DNA methylation (RdDM) pathways have been proposed to separately maintain CHH methylation in specific regions of the Arabidopsis thaliana genome. Here, we show that dysfunction of the chromatin remodeler DDM1 causes hundreds of genomic regions to switch from CMT2 dependency to RdDM dependency in DNA methylation. These converted loci are enriched at the edge regions of long transposable elements (TEs). Furthermore, we found that dysfunction in both DDM1 and RdDM causes strong reactivation of TEs and a burst of TE transposition in the first generation of mutant plants, indicating that the DDM1 and RdDM pathways together are critical to maintaining TE repression and protecting genomic stability. Our findings reveal the existence of a pathway conversion-based backup mechanism to guarantee the maintenance of DNA methylation and genome integrity.

Characterization of VOC emissions and health risk assessment in the plastic manufacturing industry.

Volatile organic compounds (VOCs) significantly contribute to ozone pollution formation, and many VOCs are known to be harmful to human health. Plastic has become an indispensable material in various industries and daily use scenarios, yet the VOC emissions and associated health risks in the plastic manufacturing industry have received limited attention. In this study, we conducted sampling in three typical plastic manufacturing factories to analyze the emission characteristics of VOCs, ozone formation potential (OFP), and health risks for workers. Isopropanol was detected at relatively high concentrations in all three factories, with concentrations in organized emissions reaching 322.3 µg/m3, 344.8 µg/m3, and 22.6 µg/m3, respectively. Alkanes are the most emitted category of VOCs in plastic factories. However, alkenes and oxygenated volatile organic compounds (OVOCs) exhibit higher OFP. In organized emissions of different types of VOCs in the three factories, alkenes and OVOCs contributed 22.8%, 67%, and 37.8% to the OFP, respectively, highlighting the necessity of controlling them. The hazard index (HI) for all three factories was less than 1, indicating a low non-carcinogenic toxic risk; however, there is still a possibility of non-cancerous health risks in two of the factories, and a potential lifetime cancer risk in all of the three factories. For workers with job tenures exceeding 5 years, there may be potential health risks, hence wearing masks with protective capabilities is necessary. This study provides evidence for reducing VOC emissions and improving management measures to ensure the health protection of workers in the plastic manufacturing industry.

Impacts of QM region sizes and conformation numbers on modelling enzyme reactions: a case study of polyethylene terephthalate hydrolase.

A quantum mechanics/molecular mechanics (QM/MM) approach is a broadly used tool in computational enzymology. Treating the QM region with a high-level DFT method is one of the important branches. Here, taking leaf-branch compost cutinase-catalyzed polyethylene terephthalate depolymerization as an example, the convergence behavior of energy barriers as well as key structural and charge features with respect to the size of the QM region (up to 1000 atoms) is systematically investigated. BP86/6-31G(d)//CHARMM and M06-2X/6-311G(d,p)//CHARMM level of theories were applied for geometry optimizations and single-point energy calculations, respectively. Six independent enzyme conformations for all the four catalytic steps (steps (i)-(iv)) were considered. Most of the twenty-four cases show that at least 500 QM atoms are needed while only two rare cases show that ∼100 QM atoms are sufficient for convergence when only a single conformation was considered. This explains why most previous studies showed that 500 or more QM atoms are required while a few others showed that ∼100 QM atoms are sufficient for DFT/MM calculations. More importantly, average energy barriers and key structural/charge features from six conformations show an accelerated convergence than that in a single conformation. For instance, to reach energy barrier convergence (within 2.0 kcal mol-1) for step (ii), only ∼100 QM atoms are required if six conformations are considered while 500 or more QM atoms are needed with a single conformation. The convergence is accelerated to be more rapid if hundreds and thousands of conformations were considered, which aligns with previous findings that only several dozens of QM atoms are required for convergence with semi-empirical QM/MM MD simulations.

The effects of the gas-liquid interface and gas phase on Cl/ClO radical interaction with water molecules.

In the marine boundary layer (MBL), chlorine (Cl) and chlorine monoxide (ClO) are powerful oxidants with high concentrations. The gas-liquid interface is also ubiquitous in the MBL as a favorable site for atmospheric reactions. Understanding the role of water in Cl/ClO radical chemistry is essential for predicting their behavior in the atmosphere and developing effective strategies for mitigating their harmful effects. However, the research studies on the system of Cl/ClO radicals on the surface of water droplets are still insufficient. In previous studies, we have found unique results related to the hydroxyl radical at the interface using ab initio molecular dynamics (AIMD). In this work, we have used AIMD to investigate interactions between Cl/ClO radicals and water molecules at the gas-liquid interface. Radical mobility, radial distribution functions, coordination, and population analyses were conducted to investigate the surface preference, bonding pattern, and track Cl/ClO radicals in the water droplets. In addition, density functional theory (DFT) analysis was conducted to compare the results at the gas-liquid interface with those in the gas phase. We found that Cl/ClO radicals tend to remain near the gas-liquid interface in water droplet systems and outside of water clusters in gas phase systems. The ClO radical can form O*-H and Cl-O bonds with water molecules; however, neither the O*-O hemibond nor the Cl-H bond was detected in all systems. Different dominant structures were obtained for ClO in the interface and gas phase. The ClO radical can be bonded to one water molecule from its oxygen side, (H2O)0-Cl-O*-(H2O)1 at the interface, or to two water molecules from the chlorine and oxygen sides, (H2O)1-Cl-O*-(H2O)1 in the gas phase. Meanwhile, the Cl radical can only form a dominant structure like Cl*-(H2O)1 at the gas-liquid interface by making a Cl*-O hemibond. Providing a thorough explanation of the Cl/ClO radical behavior at the gas-liquid interface, this study will improve our understanding of the MBL's oxidizing capacity and pollution causes.

Regioselective oxidation of heterocyclic aromatic hydrocarbons catalyzed by cytochrome P450: A case study of carbazole.

Recently there are increasing interests in accurately evaluating the health effects of heterocyclic PAHs. However, the activation mechanism and possible metabolites of heterocyclic PAHs catalyzed by human CYP1A1 is still elusive to a great extent. Here, leveraged to high level QM/MM calculations, the corresponding activation pathways of a representative heterocyclic PAHs, carbazole, were systematically explored. The first stage is electrophilic addition or hydrogen abstraction from N-H group. Electrophilic addition was evidenced to be more feasible and regioselectivity at C3 and C4 sites were identified. Correlations between energy barriers and key structural/electrostatic parameters reveal that O-Cα distance and Fe-O-Cα angle are the main origin for the catalytic regioselectivity. Electrophilic addition was determined as the rate-determining step and the subsequent possible reactions include epoxidation, NIH shift (the hydrogen migration from the site of hydroxylation to the adjacent carbon) and proton shuttle. The corresponding products are epoxides, ketones and hydroxylated carbazoles, respectively. The main metabolites (hydroxylated carbazoles) are estimated to be more toxic than carbazole. The regioselectivity of carbazole activated by CYP1A1 is different from the environmental processes (gas and aqueous phase). Collectively, these results will inform the in-depth understanding the metabolic processes of heterocyclic PAHs and aid the accurate evaluation of their health effects.

Theoretical Study on the Gas Phase and Gas-Liquid Interface Reaction Mechanism of Criegee Intermediates with Glycolic Acid Sulfate.

Criegee intermediates (CIs) are important zwitterionic oxidants in the atmosphere, which affect the budget of OH radicals, amines, alcohols, organic/inorganic acids, etc. In this study, quantum chemical calculation and Born-Oppenheimer molecular dynamic (BOMD) simulation were performed to show the reaction mechanisms of C2 CIs with glycolic acid sulfate (GAS) at the gas-phase and gas-liquid interface, respectively. The results indicate that CIs can react with COOH and OSO3H groups of GAS and generate hydroperoxide products. Intramolecular proton transfer reactions occurred in the simulations. Moreover, GAS acts as a proton donor and participates in the hydration of CIs, during which the intramolecular proton transfer also occurs. As GAS widely exists in atmospheric particulate matter, the reaction with GAS is one of the sink pathways of CIs in areas polluted by particulate matter.

Theoretical Study on the Gas-Phase and Aqueous Interface Reaction Mechanism of Criegee Intermediates with 2-Methylglyceric Acid and the Nucleation of Products.

Criegee intermediates (CIs) are important in the sink of many atmospheric substances, including alcohols, organic acids, amines, etc. In this work, the density functional theory (DFT) method was used to calculate the energy barriers for the reactions of CH3CHOO with 2-methyl glyceric acid (MGA) and to evaluate the interaction of the three functional groups of MGA. The results show that the reactions involving the COOH group of MGA are negligibly affected, and that hydrogen bonding can affect the reactions involving α-OH and ß-OH groups. The water molecule has a negative effect on the reactions of the COOH group. It decreases the energy barriers of reactions involving the α-OH and ß-OH groups as a catalyst. The Born-Oppenheimer molecular dynamic (BOMD) was applied to simulate the reactions of CH3CHOO with MGA at the gas-liquid interface. Water molecule plays the role of proton transfer in the reaction. Gas-phase calculations and gas-liquid interface simulations demonstrate that the reaction of CH3CHOO with the COOH group is the main pathway in the atmosphere. The molecular dynamic (MD) simulations suggest that the reaction products can form clusters in the atmosphere to participate in the formation of particles.

Unexpected enhancement of sulfuric acid-driven new particle formation by alcoholic amines: The role of ion-induced nucleation.

New particle formation (NPF) contributes more than half of the global aerosol. Diethanolamine (DEA) and methyldiethanolamine (MDEA) are the most common amines used to remove CO2 and H2S, which are lost to the atmosphere from CO2 chemical absorbers, livestock and consumer products and are involved in sulfuric acid (SA)-driven NPF. Ion-induced nucleation (IIN) is an important nucleation pathway for NPF. We investigated the role of IIN on DEA and MDEA enhancing SA-driven NPF using density functional method (DFT), molecular dynamics (MD) simulation and atmospheric cluster dynamics code (ACDC). The effects of SO42-, H3O+, NH4+, HSO4-, NO3-, ammonia, methylamine, dimethylamine, trimethylamine and water (W) on the nucleation of SA-DEA were further investigated. The enhancement ability of DEA is greater than that of dimethylamine (DMA) and MDEA. Participation in SA-based NPF is a removal pathway for DEA and MDEA. DEA-SA clusters are generated that not only aggregate DEA and SA molecules, but also increase further growth of atmospheric ions. The very low Gibbs formation free energy highlights the importance of ion-induced nucleation for SA-based NPF. The order of the ability of common atmospheric ions to increase the (SA)(DEA) cluster nucleation is SO42- > H3O+ > NH4+ > HSO4- > NO3-. The addition of 20 water molecules increases the (SA)(DEA)9 cluster from 1.882 nm to 2.053 nm, promoting SA-based NPF. The atmospheric ions accelerate the aggregation rate of the (SA)5(DEA)5 cluster within 15 ns?

Mitotically heritable epigenetic modifications of CmMYB6 control anthocyanin biosynthesis in chrysanthemum.

Flower color, which is determined by various chemical pigments, is a vital trait for ornamental plants, in which anthocyanin is a major component. However, the epigenetic regulation of anthocyanin biosynthesis remains poorly understood. During chrysanthemum cultivation, we found a heterochromatic chrysanthemum accession (YP) whose progeny generated by asexual reproduction contained both yellow-flowered (YP-Y) and pink-flowered (YP-P) plants. In this study, we aimed to elucidate the epigenetic mechanisms of different flower colors in the YP plant progeny. Metabolome and transcriptome analyses revealed that the difference in flower color between YP-Y and YP-P was caused by expression variation of the anthocyanin biosynthesis gene CmMYB6. Bisulfite sequencing revealed that methylation at the CmMYB6 promoter, especially in the CHH context, was higher in YP-Y than YP-P. After demethylation of the CmMYB6 promoter using the dCas9-TET1cd system, the flower color returned from yellow to pink. Furthermore, the methylation status of the CmMYB6 promoter was higher in YP-Y over three consecutive generations, indicating that this methylation status was heritable mitotically. Finally, investigation of other chrysanthemum cultivars showed that the methylation of CmMYB6 decreased gradually with the increase in anthocyanin content. These results lay an epigenetic foundation for the improvement of flower color in horticultural plants.

Hydrolysis Mechanism of Carbamate Methomyl by a Novel Esterase PestE: A QM/MM Approach.

Methomyl is one of the most important carbamates that has caused potential hazardous effects on both human beings and the environment. Here, we systematically investigated the hydrolysis mechanism of methomyl catalyzed by esterase PestE using molecular dynamics simulations (MD) and quantum mechanics/molecular mechanics (QM/MM) calculations. The hydrolysis mechanism involves two elementary steps: (Ⅰ) serine-initiated nucleophilic attack and (Ⅱ) C-O bond cleavage. Our work elicits the atomic level details of the hydrolysis mechanism and free energy profiles along the reaction pathway. The Boltzmann-weighted average potential barriers are 19.1 kcal/mol and 7.5 kcal/mol for steps Ⅰ and Ⅱ, respectively. We identified serine-initiated nucleophilic attack as the rate determining-step. The deep learning-based kcat prediction model indicated that the barrier of the rate-determining step is 15.4 kcal/mol, which is in good agreement with the calculated results using Boltzmann-weighted average method. We have elucidated the importance of the protein-substrate interactions and the roles of the key active site residues during the hydrolysis process through noncovalent interactions analysis and electrostatic potential (ESP) analysis. The results provide practical value for achieving efficient degradation of carbamates by hydrolases.

Transgenerationally Transmitted DNA Demethylation of a Spontaneous Epialleles Using CRISPR/dCas9-TET1cd Targeted Epigenetic Editing in Arabidopsis.

CRISPR/dCas9 is an important DNA modification tool in which a disarmed Cas9 protein with no nuclease activity is fused with a specific DNA modifying enzyme. A previous study reported that overexpression of the TET1 catalytic domain (TET1cd) reduces genome-wide methylation in Arabidopsis. A spontaneous naturally occurring methylation region (NMR19-4) was identified in the promoter region of the PPH (Pheophytin Pheophorbide Hydrolase) gene, which encodes an enzyme that can degrade chlorophyll and accelerate leaf senescence. The methylation status of NMR19-4 is associated with PPH expression and leaf senescence in Arabidopsis natural accessions. In this study, we show that the CRISPR/dCas9-TET1cd system can be used to target the methylation of hypermethylated NMR19-4 region to reduce the level of methylation, thereby increasing the expression of PPH and accelerating leaf senescence. Furthermore, hybridization between transgenic demethylated plants and hypermethylated ecotypes showed that the demethylation status of edited NMR19-4, along with the enhanced PPH expression and accelerated leaf senescence, showed Mendelian inheritance in F1 and F2 progeny, indicating that spontaneous epialleles are stably transmitted trans-generationally after demethylation editing. Our results provide a rational approach for future editing of spontaneously mutated epialleles and provide insights into the epigenetic mechanisms that control plant leaf senescence.

Interaction Structure and Affinity of Zwitterionic Amino Acids with Important Metal Cations (Cd2+, Cu2+, Fe3+, Hg2+, Mn2+, Ni2+ and Zn2+) in Aqueous Solution: A Theoretical Study.

Heavy metals are non-biodegradable and carcinogenic pollutants with great bio-accumulation potential. Their ubiquitous occurrence in water and soils has caused serious environmental concerns. Effective strategies that can eliminate the heavy metal pollution are urgently needed. Here the adsorption potential of seven heavy metal cations (Cd2+, Cu2+, Fe3+, Hg2+, Mn2+, Ni2+ and Zn2+) with 20 amino acids was systematically investigated with Density Functional Theory method. The binding energies calculated at B3LYP-D3/def2TZVP level showed that the contribution order of amino acid side chains to the binding affinity was carboxyl > benzene ring > hydroxyl > sulfhydryl > amino group. The affinity order was inversely proportional to the radius and charge transfer of heavy metal cations, approximately following the order of: Ni2+ > Fe3+ > Cu2+ > Hg2+ > Zn2+ > Cd2+ > Mn2+. Compared to the gas-phase in other researches, the water environment has a significant influence on structures and binding energies of the heavy metal and amino acid binary complexes. Collectively, the present results will provide a basis for the design of a chelating agent (e.g., adding carboxyl or a benzene ring) to effectively remove heavy metals from the environment.

Genome-wide identification, classification, and expression analysis of the JmjC domain-containing histone demethylase gene family in birch.

BACKGROUND: Histone methylation occurs primarily on lysine residues and requires a set of enzymes capable of reading, writing, and erasing to control its establishment and deletion, which is essential for maintaining chromatin structure and gene expression. Histone methylation and demethylation are contributed to plant growth and development, and are involved in adapting to environmental stresses. The JmjC domain-containing proteins are extensively studied for their function in histone lysine demethylation in plants, and play a critical role in sustaining histone methylation homeostasis. RESULTS: In this study, a total of 21 JmjC domain-containing histone demethylase proteins (JHDMs) in birch were identified and classified into five subfamilies based on structural characteristics and phylogenetic relationships among Arabidopsis, rice, maize, and birch. Although the BpJMJ genes displayed significant schematic variation, their distribution on the chromosomes is relatively uniform. Additionally, the BpJMJ genes in birch have never experienced a tandem-duplication event proved by WGD analysis and were remaining underwent purifying selection (Ka/Ks < < 1). A typical JmjC domain was found in all BpJMJ genes, some of which have other essential domains for their functions. In the promoter regions of BpJMJ genes, cis-acting elements associated with hormone and abiotic stress responses were overrepresented. Under abiotic stresses, the transcriptome profile reveals two contrasting expression patterns within 21 BpJMJ genes. Furthermore, it was established that most BpJMJ genes had higher expression in young tissues under normal conditions, with BpJMJ06/16 having the highest expression in germinating seeds and participating in the regulation of BpGA3ox1/2 gene expression. Eventually, BpJMJ genes were found to directly interact with genes involved in the "intracellular membrane" in respond to cold stress. CONCLUSIONS: The present study will provide a foundation for future experiments on histone demethylases in birch and a theoretical basis for epigenetic research on growth and development in response to abiotic stresses.

Genome-wide identification and expression analysis of the AT-hook Motif Nuclear Localized gene family in soybean.

BACKGROUND: Soybean is an important legume crop and has significant agricultural and economic value. Previous research has shown that the AT-Hook Motif Nuclear Localized (AHL) gene family is highly conserved in land plants, playing crucial roles in plant growth and development. To date, however, the AHL gene family has not been studied in soybean. RESULTS: To investigate the roles played by the AHL gene family in soybean, genome-wide identification, expression patterns and gene structures were performed to analyze. We identified a total of 63 AT-hook motif genes, which were characterized by the presence of the AT-hook motif and PPC domain in soybean. The AT-hook motif genes were distributed on 18 chromosomes and formed two distinct clades (A and B), as shown by phylogenetic analysis. All the AHL proteins were further classified into three types (I, II and III) based on the AT-hook motif. Type-I was belonged to Clade-A, while Type-II and Type-III were belonged to Clade-B. Our results also showed that the main type of duplication in the soybean AHL gene family was segmented duplication event. To discern whether the AHL gene family was involved in stress response in soybean, we performed cis-acting elements analysis and found that AHL genes were associated with light responsiveness, anaerobic induction, MYB and gibberellin-responsiveness elements. This suggest that AHL genes may participate in plant development and mediate stress response. Moreover, a co-expression network analysis showed that the AHL genes were also involved in energy transduction, and the associated with the gibberellin pathway and nuclear entry signal pathways in soybean. Transcription analysis revealed that AHL genes in Jack and Williams82 have a common expression pattern and are mostly expressed in roots, showing greater sensitivity under drought and submergence stress. Hence, the AHL gene family mainly reacts on mediating stress responses in the roots and provide comprehensive information for further understanding of the AT-hook motif gene family-mediated stress response in soybean. CONCLUSION: Sixty-three AT-hook motif genes were identified in the soybean genome. These genes formed into two distinct phylogenetic clades and belonged to three different types. Cis-acting elements and co-expression network analyses suggested that AHL genes participated in significant biological processes. This work provides important theoretical basis for the understanding of AHLs biological functions in soybean.

Loss of chromatin remodeler DDM1 causes segregation distortion in Arabidopsis thaliana.

MAIN CONCLUSION: In ddm1 mutants, the DNA methylation is primarily affected in the heterochromatic region of the chromosomes, which is associated with the segregation distortion of SNPs in the F2 progenies. Segregation distortion (SD) is common in most genetic mapping experiments and a valuable resource to determine how gene loci induce deviation. Meiotic DNA crossing over and SD are under the control of several types of epigenetic modifications. DNA methylation is an important regulatory epigenetic modification that is inherited across generations. In the present study, we investigated the relationship between SD and DNA methylation. The ecotypes Col-0/C24 and chromatin remodeler mutants ddm1-10/Col and ddm1-15/C24 were reciprocally crossed to obtain F2 generations. A total of 300 plants for each reciprocally crossed plant in the F2 generations were subjected to next-generation sequencing to detect the single-nucleotide polymorphisms (SNPs) as DNA markers. All SNPs were analyzed using the Chi-square test method to determine their segregation ratio in F2 generations. Through the segregation ratio, whole-genome SNPs were classified into 16 classes. In class 10, the SNPs in the reciprocal crosses of wild type showed the expected Mendelian ratio of 1:2:1, while those in the reciprocal crosses of ddm1 mutants showed distortion. In contrast, all SNPs in class 16 displayed a normal 1:2:1 ratio, and class 1 showed SD, regardless of wild type or mutants, as assessed using CAPS (cleaved amplified polymorphic sequences) marker analysis to confirm the next-generation sequencing. In ddm1 mutants, the DNA methylation is highly reduced throughout the whole genome and more significantly in the heterochromatic regions of chromosomes. Our results showed that the ddm1 mutants exhibit low levels of DNA methylation, which facilitates the SD of SNPs primarily located in the heterochromatic region of chromosomes by reducing the heterozygous ratio. The present study will provide a strong base for future research focusing on the impact of DNA methylation on trait segregation and plant evolution.

Comprehensive Understanding of Fluoroacetate Dehalogenase-Catalyzed Degradation of Fluorocarboxylic Acids: A QM/MM Approach.

Fluorochemicals are persistent, bioaccumulative, and toxic compounds that are widely tributed in the environment. Developing efficient biodegradation strategies to decompose the fluorochemicals via breaking the inert C-F bonds presents a holistic challenge. As a promising biodegradation enzyme candidate, fluoroacetate dehalogenase (FAcD) has been reported as the only non-metallic enzyme to catalyze the cleavage of the strong C-F bond. Here, we systematically investigated the catalytic actions of FAcD toward its natural substrate fluoroacetate using molecular dynamics simulations and quantum mechanism/molecular mechanism calculations. We propose that the enzymatic transformation involves four elementary steps, (I) C-F bond activation, (II) nucleophilic attack, (III) C-O bond cleavage, and (IV) proton transfer. Our results show that nucleophilic attack is the rate-determining step. However, for difluoroacetate and trifluoroacetate, C-F bond activation, instead of nucleophilic attack, becomes the rate-determining step. We show that FAcD, originally recognized as α-fluorocarboxylic acid degradation enzyme, can catalyze the defluorination of difluoroacetate to glyoxylate, which is captured by our high-resolution mass spectrometry experiments. In addition, we employed amino acid electrostatic analysis method to screen potential mutation hotspots for tuning FAcD's electrostatic environment to favor substrate conversion. The comprehensive understanding of catalytic mechanism will inform a rational enzyme engineering strategy to degrade fluorochemicals for benefits of environmental sustainability.