Cas12n nucleases, early evolutionary intermediates of type V CRISPR, comprise a distinct family of miniature genome editors.

Type V CRISPR-associated systems (Cas)12 family nucleases are considered to have evolved from transposon-associated TnpB, and several of these nucleases have been engineered as versatile genome editors. Despite the conserved RNA-guided DNA-cleaving functionality, these Cas12 nucleases differ markedly from the currently identified ancestor TnpB in aspects such as guide RNA origination, effector complex composition, and protospacer adjacent motif (PAM) specificity, suggesting the presence of earlier evolutionary intermediates that could be mined to develop advanced genome manipulation biotechnologies. Using evolutionary and biochemical analyses, we identify that the miniature type V-U4 nuclease (referred to as Cas12n, 400-700 amino acids) is likely the earliest evolutionary intermediate between TnpB and large type V CRISPR systems. We demonstrate that with the exception of CRISPR array emergence, CRISPR-Cas12n shares several similar characteristics with TnpB-ωRNA, including a miniature and likely monomeric nuclease for DNA targeting, origination of guide RNA from nuclease coding sequence, and generation of a small sticky end following DNA cleavage. Cas12n nucleases recognize a unique 5'-AAN PAM sequence, of which the A nucleotide at the -2 position is also required for TnpB. Moreover, we demonstrate the robust genome-editing capacity of Cas12n in bacteria and engineer a highly efficient CRISPR-Cas12n (termed Cas12Pro) with up to 80% indel efficiency in human cells. The engineered Cas12Pro enables base editing in human cells. Our results further expand the understanding regarding type V CRISPR evolutionary mechanisms and enrich the miniature CRISPR toolbox for therapeutic applications.

A KDPG sensor RccR governs Pseudomonas aeruginosa carbon metabolism and aminoglycoside antibiotic tolerance.

Pseudomonas aeruginosa harbors sophisticated transcription factor (TF) networks to coordinately regulate cellular metabolic states for rapidly adapting to changing environments. The extraordinary capacity in fine-tuning the metabolic states enables its success in tolerance to antibiotics and evading host immune defenses. However, the linkage among transcriptional regulation, metabolic states and antibiotic tolerance in P. aeruginosa remains largely unclear. By screening the P. aeruginosa TF mutant library constructed by CRISPR/Cas12k-guided transposase, we identify that rccR (PA5438) is a major genetic determinant in aminoglycoside antibiotic tolerance, the deletion of which substantially enhances bacterial tolerance. We further reveal the inhibitory roles of RccR in pyruvate metabolism (aceE/F) and glyoxylate shunt pathway (aceA and glcB), and overexpression of aceA or glcB enhances bacterial tolerance. Moreover, we identify that 2-keto-3-deoxy-6-phosphogluconate (KDPG) is a signal molecule that directly binds to RccR. Structural analysis of the RccR/KDPG complex reveals the detailed interactions. Substitution of the key residue R152, K270 or R277 with alanine abolishes KDPG sensing by RccR and impairs bacterial growth with glycerol or glucose as the sole carbon source. Collectively, our study unveils the connection between aminoglycoside antibiotic tolerance and RccR-mediated central carbon metabolism regulation in P. aeruginosa, and elucidates the KDPG-sensing mechanism by RccR.

Delta-5 elongase knockout reduces docosahexaenoic acid and lipid synthesis and increases heat sensitivity in a diatom.

Recent global marine lipidomic analysis reveals a strong relationship between ocean temperature and phytoplanktonic abundance of omega-3 long-chain polyunsaturated fatty acids (LC-PUFAs), especially eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA), which are essential for human nutrition and primarily sourced from phytoplankton in marine food webs. In phytoplanktonic organisms, EPA may play a major role in regulating the phase transition temperature of membranes, while the function of DHA remains unexplored. In the oleaginous diatom Phaeodactylum tricornutum, DHA is distributed mainly on extraplastidial phospholipids, which is very different from the EPA enriched in thylakoid lipids. Here, CRISPR/Cas9-mediated knockout of delta-5 elongase (ptELO5a), which encodes a delta-5 elongase (ELO5) catalyzing the elongation of EPA to synthesize DHA, led to a substantial interruption of DHA synthesis in P. tricornutum. The ptELO5a mutants showed some alterations in transcriptome and glycerolipidomes, including membrane lipids and triacylglycerols under normal temperature (22°C), and were more sensitive to elevated temperature (28°C) than wild type. We conclude that PtELO5a-mediated synthesis of small amounts of DHA has indispensable functions in regulating membrane lipids, indirectly contributing to storage lipid accumulation, and maintaining thermomorphogenesis in P. tricornutum. This study also highlights the significance of DHA synthesis and lipid composition for environmental adaptation of P. tricornutum.

Molecular basis for cell-wall recycling regulation by transcriptional repressor MurR in Escherichia coli.

The cell-wall recycling process is important for bacterial survival in nutrient-limited conditions and, in certain cases, is directly involved in antibiotic resistance. In the sophisticated cell-wall recycling process in Escherichia coli, the transcriptional repressor MurR controls the expression of murP and murQ, which are involved in transporting and metabolizing N-acetylmuramic acid (MurNAc), generating N-acetylmuramic acid-6-phosphate (MurNAc-6-P) and N-acetylglucosamine-6-phosphate (GlcNAc-6-P). Here, we report that both MurNAc-6-P and GlcNAc-6-P can bind to MurR and weaken the DNA binding ability of MurR. Structural characterizations of MurR in complex with MurNAc-6-P or GlcNAc-6-P as well as in the apo form revealed the detailed ligand recognition chemistries. Further studies showed that only MurNAc-6-P, but not GlcNAc-6-P, is capable of derepressing the expression of murQP controlled by MurR in cells and clarified the substrate specificity through the identification of key residues responsible for ligand binding in the complex structures. In summary, this study deciphered the molecular mechanism of the cell wall recycling process regulated by MurR in E. coli.

Programmed genome editing by a miniature CRISPR-Cas12f nuclease.

The RNA-guided CRISPR-associated (Cas) nucleases are versatile tools for genome editing in various organisms. The large sizes of the commonly used Cas9 and Cas12a nucleases restrict their flexibility in therapeutic applications that use the cargo-size-limited adeno-associated virus delivery vehicle. More compact systems would thus offer more therapeutic options and functionality for this field. Here, we report a miniature class 2 type V-F CRISPR-Cas genome-editing system from Acidibacillus sulfuroxidans (AsCas12f1, 422 amino acids). AsCas12f1 is an RNA-guided endonuclease that recognizes 5' T-rich protospacer adjacent motifs and creates staggered double-stranded breaks to target DNA. We show that AsCas12f1 functions as an effective genome-editing tool in both bacteria and human cells using various delivery methods, including plasmid, ribonucleoprotein and adeno-associated virus. The small size of AsCas12f1 offers advantages for cellular delivery, and characterizations of AsCas12f1 may facilitate engineering more compact genome-manipulation technologies.

Molecular basis for the PAM expansion and fidelity enhancement of an evolved Cas9 nuclease.

Clustered regularly interspaced short palindromic repeats (CRISPR)-Cas systems have been harnessed as powerful genome editing tools in diverse organisms. However, the off-target effects and the protospacer adjacent motif (PAM) compatibility restrict the therapeutic applications of these systems. Recently, a Streptococcus pyogenes Cas9 (SpCas9) variant, xCas9, was evolved to possess both broad PAM compatibility and high DNA fidelity. Through determination of multiple xCas9 structures, which are all in complex with single-guide RNA (sgRNA) and double-stranded DNA containing different PAM sequences (TGG, CGG, TGA, and TGC), we decipher the molecular mechanisms of the PAM expansion and fidelity enhancement of xCas9. xCas9 follows a unique two-mode PAM recognition mechanism. For non-NGG PAM recognition, xCas9 triggers a notable structural rearrangement in the DNA recognition domains and a rotation in the key PAM-interacting residue R1335; such mechanism has not been observed in the wild-type (WT) SpCas9. For NGG PAM recognition, xCas9 applies a strategy similar to WT SpCas9. Moreover, biochemical and cell-based genome editing experiments pinpointed the critical roles of the E1219V mutation for PAM expansion and the R324L, S409I, and M694I mutations for fidelity enhancement. The molecular-level characterizations of the xCas9 nuclease provide critical insights into the mechanisms of the PAM expansion and fidelity enhancement of xCas9 and could further facilitate the engineering of SpCas9 and other Cas9 orthologs.

Selective cadmium regulation mediated by a cooperative binding mechanism in CadR.

Detoxification of the highly toxic cadmium element is essential for the survival of living organisms. Pseudomonas putida CadR, a MerR family transcriptional regulator, has been reported to exhibit an ultraspecific response to the cadmium ion. Our crystallographic and spectroscopic studies reveal that the extra cadmium selectivity of CadR is mediated by the unexpected cooperation of thiolate-rich site I and histidine-rich site II. Cadmium binding in site I mediates the reorientation of protein domains and facilitates the assembly of site II. Subsequently, site II bridge-links 2 DNA binding domains through ligands His140/His145 in the C-terminal histidine-rich tail. With dynamic transit between 2 conformational states, this bridge could stabilize the regulator into an optimal conformation that is critical for enhancing the transcriptional activity of the cadmium detoxification system. Our results provide dynamic insight into how nature utilizes the unique cooperative binding mechanism in multisite proteins to recognize cadmium ions specifically.

Genome-wide colocalization of RNA-DNA interactions and fusion RNA pairs.

Fusion transcripts are used as biomarkers in companion diagnoses. Although more than 15,000 fusion RNAs have been identified from diverse cancer types, few common features have been reported. Here, we compared 16,410 fusion transcripts detected in cancer (from a published cohort of 9,966 tumor samples of 33 cancer types) with genome-wide RNA-DNA interactions mapped in two normal, noncancerous cell types [using iMARGI, an enhanced version of the mapping of RNA-genome interactions (MARGI) assay]. Among the top 10 most significant RNA-DNA interactions in normal cells, 5 colocalized with the gene pairs that formed fusion RNAs in cancer. Furthermore, throughout the genome, the frequency of a gene pair to exhibit RNA-DNA interactions is positively correlated with the probability of this gene pair to present documented fusion transcripts in cancer. To test whether RNA-DNA interactions in normal cells are predictive of fusion RNAs, we analyzed these in a validation cohort of 96 lung cancer samples using RNA sequencing (RNA-seq). Thirty-seven of 42 fusion transcripts in the validation cohort were found to exhibit RNA-DNA interactions in normal cells. Finally, by combining RNA-seq, single-molecule RNA FISH, and DNA FISH, we detected a cancer sample with EML4-ALK fusion RNA without forming the EML4-ALK fusion gene. Collectively, these data suggest an RNA-poise model, where spatial proximity of RNA and DNA could poise for the creation of fusion transcripts.

Pulsation, translation and P1 deformation of two aspherical bubbles in liquid.

In this work, the interactions between the axial translational motions and aspherical oscillations of two gas bubbles in an incompressible liquid are considered. Representing the surface function by the Legendre polynomial of first order, we derive a dynamic model to describe the motions of two aspherical bubbles in Lagrangian mechanics. An apple-shaped bubble from simulations based on the model can be well consistent with known experimental observation. The bubble appears as the shape of a sphere at maximum expansion. The maximum asymmetry of the bubbles occurs during collapse. The surface tension is a key factor to stable oscillatory deformation. It is also found that the aspherical amplitudes of two bubbles decrease with increasing distance or decreasing driving pressure.

The left-right symmetrical and asymmetrical deformations in a three-bubble system.

This paper studies the simplest system that can possess left-right symmetrical and asymmetrical surroundings, three bubbles in a line. Assuming that the deformations are small, the surfaces of bubbles are described by a combination of the first three Legendre polynomials, that is, spherical symmetrical mode P0, L-R antisymmetrical mode P1, and symmetrical mode P2. A dynamical model is built to describe aspherical oscillations of central and two side bubbles. It is found that when three identical bubbles are separated uniformly, the central bubble only has a P2 component and P1 component tends to zero, while two side bubbles have both P1 and P2 components. When three identical bubbles are separated by different distances, they can be degenerated into a two-bubble system and a free bubble. The bubble deformations contain both P1 and P2 components in the two-bubble system, while both aspherical components P1 and P2 of the free bubble tend to zero. If side bubbles are different in ambient radii but located symmetrically on the left and right of the central bubble, the side bubble pulsated more strongly plays an important role on the deformation of the central one.

Auxilin-like protein MoSwa2 promotes effector secretion and virulence as a clathrin uncoating factor in the rice blast fungus Magnaporthe oryzae.

Plant pathogens exploit the extracellular matrix (ECM) to inhibit host immunity during their interactions with the host. The formation of ECM involves a series of continuous steps of vesicular transport events. To understand how such vesicle trafficking impacts ECM and virulence in the rice blast fungus Magnaporthe oryzae, we characterised MoSwa2, a previously identified actin-regulating kinase MoArk1 interacting protein, as an orthologue of the auxilin-like clathrin uncoating factor Swa2 of the budding yeast Saccharomyces cerevisiae. We found that MoSwa2 functions as an uncoating factor of the coat protein complex II (COPII) via an interaction with the COPII subunit MoSec24-2. Loss of MoSwa2 led to a deficiency in the secretion of extracellular proteins, resulting in both restricted growth of invasive hyphae and reduced inhibition of host immunity. Additionally, extracellular fluid (ECF) proteome analysis revealed that MoSwa2-regulated extracellular proteins include many redox proteins such as the berberine bridge enzyme-like (BBE-like) protein MoSef1. We further found that MoSef1 functions as an apoplastic virulent factor that inhibits the host immune response. Our studies revealed a novel function of a COPII uncoating factor in vesicular transport that is critical in the suppression of host immunity and pathogenicity of M. oryzae.

Differential nucleosome occupancy modulates alternative splicing in Arabidopsis thaliana.

Alternative splicing (AS) is a major gene regulatory mechanism in plants. Recent evidence supports co-transcriptional splicing in plants, hence the chromatin state can impact AS. However, how dynamic changes in the chromatin state such as nucleosome occupancy influence the cold-induced AS remains poorly understood. Here, we generated transcriptome (RNA-Seq) and nucleosome positioning (MNase-Seq) data for Arabidopsis thaliana to understand how nucleosome positioning modulates cold-induced AS. Our results show that characteristic nucleosome occupancy levels are strongly associated with the type and abundance of various AS events under normal and cold temperature conditions in Arabidopsis. Intriguingly, exitrons, alternatively spliced internal regions of protein-coding exons, exhibit distinctive nucleosome positioning pattern compared to other alternatively spliced regions. Likewise, nucleosome patterns differ between exitrons and retained introns, pointing to their distinct regulation. Collectively, our data show that characteristic changes in nucleosome positioning modulate AS in plants in response to cold.

Electronic characteristics of BRCA1 mutations in DNA.

The efficient and low-cost way for gene mutation detection and identification are conducive for the detection of disease. Here, we report the electronic characteristics of the gene of breast cancer 1 in four common mutation types: duplication, single nucleotide variant, deletion, and indel. The electronic characteristics are investigated by the combination of density functional theory and non-equilibrium Green's function formulation with decoherence. The magnitude of conductance of these DNA molecules and mutational changes are found to be detectable experimentally. In this study, we also find the significant mutation type dependent on the change of conductance. Hence these mutations are expected to be identifiable. We find deletion type mutation shows the largest change in relative conductance (~97%), whereas the indel mutation shows the smallest change in relative conductance (~27%). Therefore, this work presents a possibility of electronic detection and identification of mutations in DNA, which could be an efficient method as compared to the conventional methods.

Effect of cytosine hydroxymethylation on DNA charge transport.

DNA hydroxymethylation plays a very important role in some biological processes, such as DNA methylation process. In addition, its presence can also cause some diseases. In this paper, the electrical properties of cytosine hydroxymethylated (Chm) DNA sequences are studied. The density functional theory (DFT) and Landauer-Büttiker framework are used to study the decoherence conductance and transmission of the Chm strands in different configurations, which provides a theoretical basis for the detection of Chm. The results show that the conductance of the hydroxymethylated DNA strand is smaller than that of the native and methylated strands. The length dependence of the Chm strands is also studied. With the length increasing, the conductance becomes larger. This study shows that DNA methylation can be detected electrically.

Genetic diversity and natural selection on the thrombospondin-related adhesive protein (TRAP) gene of Plasmodium falciparum on Bioko Island, Equatorial Guinea and global comparative analysis.

BACKGROUND: Thrombospondin-related adhesive protein (TRAP) is a transmembrane protein that plays a crucial role during the invasion of Plasmodium falciparum into liver cells. As a potential malaria vaccine candidate, the genetic diversity and natural selection of PfTRAP was assessed and the global PfTRAP polymorphism pattern was described. METHODS: 153 blood spot samples from Bioko malaria patients were collected during 2016-2018 and the target TRAP gene was amplified. Together with the sequences from database, nucleotide diversity and natural selection analysis, and the structural prediction were preformed using bioinformatical tools. RESULTS: A total of 119 Bioko PfTRAP sequences were amplified successfully. On Bioko Island, PfTRAP shows its high degree of genetic diversity and heterogeneity, with π value for 0.01046 and Hd for 0.99. The value of dN-dS (6.2231, p < 0.05) hinted at natural selection of PfTRAP on Bioko Island. Globally, the African PfTRAPs showed more diverse than the Asian ones, and significant genetic differentiation was discovered by the fixation index between African and Asian countries (Fst > 0.15, p < 0.05). 667 Asian isolates clustered in 136 haplotypes and 739 African isolates clustered in 528 haplotypes by network analysis. The mutations I116T, L221I, Y128F, G228V and P299S were predicted as probably damaging by PolyPhen online service, while mutations L49V, R285G, R285S, P299S and K421N would lead to a significant increase of free energy difference (ΔΔG > 1) indicated a destabilization of protein structure. CONCLUSIONS: Evidences in the present investigation supported that PfTRAP gene from Bioko Island and other malaria endemic countries is highly polymorphic (especially at T cell epitopes), which provided the genetic information background for developing an PfTRAP-based universal effective vaccine. Moreover, some mutations have been shown to be detrimental to the protein structure or function and deserve further study and continuous monitoring.

Mechanistic insights into staphylopine-mediated metal acquisition.

Metal acquisition is vital to pathogens for successful infection within hosts. Staphylopine (StP), a broad-spectrum metallophore biosynthesized by the major human pathogen, Staphylococcus aureus, plays a central role in transition-metal acquisition and bacterial virulence. The StP-like biosynthesis loci are present in various pathogens, and the proteins responsible for StP/metal transportation have been determined. However, the molecular mechanisms of how StP/metal complexes are recognized and transported remain unknown. We report multiple structures of the extracytoplasmic solute-binding protein CntA from the StP/metal transportation system in apo form and in complex with StP and three different metals. We elucidated a sophisticated metal-bound StP recognition mechanism and determined that StP/metal binding triggers a notable interdomain conformational change in CntA. Furthermore, CRISPR/Cas9-mediated single-base substitution mutations and biochemical analysis highlight the importance of StP/metal recognition for StP/metal acquisition. These discoveries provide critical insights into the study of novel metal-acquisition mechanisms in microbes.

Two-Step Measurement of Water-Energy-Food Symbiotic Coordination and Identification of Key Influencing Factors in the Yangtze River Basin.

With the intensification of people's production and life behaviors, the systemic risks of water, energy and food in the Yangtze River Basin have become increasingly prominent, which has become a bottleneck for sustainable development of social, economic and ecological in the basin. Therefore, studying the symbiotic coordination between water, energy and food is of great significance to promoting regional sustainable development. First, from the perspective of water-energy-food symbiosis, with the water-energy-food ecosystem conceptual model as the nexus, the two-step measurement model of the symbiotic index and the symbiotic level index is used to study the water-energy-food symbiosis of the Yangtze River. Then, we use the BP-DEMATEL-GTCW model to identify the key influencing factors that affect the symbiotic security of the water-energy-food ecosystem. In this research, it is found that the average value of the symbiotic degree of the water-energy-food ecosystem of the 11 provinces or municipalities in the Yangtze River Basin only reached the risk grade. It can also be seen from the identification results of key influencing factors that energy microsystem-related indicators have a greater impact on the symbiotic development of the entire WEF ecosystem. Therefore, special attention needs to be paid to increasing energy sources and reducing expenditure. Relevant departments need to effectively develop primary energy production and expand energy-saving investment through multiple channels to expand energy self-sufficiency and ultimately promote the coordinated and effective development of water, energy and food in the Yangtze River Basin.

Assessment and Prediction of Water Resources Vulnerability Based on a NRS-RF Model: A Case Study of the Song-Liao River Basin, China.

The vulnerability of water resources is an important criterion for evaluating the carrying capacity of water resources systems under the influence of climate change and human activities. Moreover, assessment and prediction of river basins' water resources vulnerability are important means to assess the water resources security state of river basins and identify possible problems in future water resources systems. Based on the constructed indicator system of water resources vulnerability assessment in Song-Liao River Basin, this paper uses the neighborhood rough set (abbreviated as NRS) method to reduce the dimensionality of the original indicator system to remove redundant attributes. Then, assessment indicators' standard values after dimensionality reduction are taken as the evaluation sample, and the random forest regression (abbreviated as RF) model is used to assess the water resources vulnerability of the river basin. Finally, based on data under three different future climate and socio-economic scenarios, scenario predictions are made on the vulnerability of future water resources. The results show that the overall water resources vulnerability of the Song-Liao River Basin has not improved significantly in the past 18 years, and the overall vulnerability of the Song-Liao River Basin is in the level V of moderate to high vulnerability. In the future scenario 1, the overall water resources vulnerability of the river basin will improve, and it is expected to achieve an improvement to the level III of moderate to low vulnerability. At the same time, the natural vulnerability and vulnerability of carrying capacity will increase significantly in the future, and the man-made vulnerability will increase slowly, which will deteriorate to the level V of moderate to high vulnerability under Scenario 3. Therefore, taking active measures can significantly reduce the vulnerability of nature and carrying capacity, but man-made vulnerability will become a bottleneck restricting the fragility of the overall water resources of the river basin in the future.

Structural Basis of Staphylococcus aureus Surface Protein SdrC.

Staphylococcus aureus surface proteins play important roles in host tissue colonization, biofilm formation, and bacterial virulence and are thus essential for successful host infections. The surface protein SdrC from S. aureus induces bacterial biofilm formation via an intermolecular homophilic interaction of its N2 domains. However, the molecular mechanism of how the homophilic interaction is achieved is unknown. Here, we report two crystal structures of SdrC N2N3 domains, revealing two possible homophilic interaction mechanisms: Ca2+-mediated intermolecular metal chelation of N2 domains and intermolecular interaction of N2 and N3 domains. Given the unnecessary role of the N3 domain in the induction of biofilm formation, the N2 domain-mediated metal chelation mechanism is likely the mechanism that facilitates SdrC homophilic interaction. Mutation of key Ca2+-chelating residues differentially reduced the level of protein dimer formation, further supporting the key role of metal chelation in the N2 domain interaction. Together, these results reveal the possible mechanism of the homophilic interaction of SdrC N2 domains and pave the way for the rational development of new strategies against this mechanism.

Metabolism and Toxicity of Emodin: Genome-Wide Association Studies Reveal Hepatocyte Nuclear Factor 4α Regulates UGT2B7 and Emodin Glucuronidation.

Emodin is the main toxic component in Chinese medicinal herbs such as rhubarb. Our previous studies demonstrated that genetic polymorphisms of UDP-glucuronosyltransferase 2B7 (UGT2B7) had an effect on the glucuronidation and detoxification of emodin. This study aimed to reveal the transcriptional regulation mechanism of UGT2B7 on emodin glucuronidation and its effect on toxicity. Emodin glucuronic activity and genome and transcriptome data were obtained from 36 clinical human kidney tissues. The genome-wide association studies (GWAS) identified that four single nucleotide polymorphisms (SNPs) (rs6093966, rs2868094, rs2071197, and rs6073433), which were located on the hepatocyte nuclear factor 4α (HNF4A) gene, were significantly associated with the emodin glucuronidation (p < 0.05). Notably, rs2071197 was significantly associated with the gene expression of HNF4A and UGT2B7 and the glucuronidation of emodin. The gene expression of HNF4A showed a high correlation with UGT2B7 (R2 = 0.721, p = 5.83 × 10-11). The luciferase activity was increased 7.68-fold in 293T cells and 2.03-fold in HepG2 cells, confirming a significant transcriptional activation of UGT2B7 promoter by HNF4A. The knockdown of HNF4A in HepG2 cells (36.6%) led to a significant decrease of UGT2B7 (19.8%) and higher cytotoxicity (p < 0.05). The overexpression of HNF4A in HepG2 cells (31.2%) led to a significant increase of UGT2B7 (24.4%) and improved cell viability (p < 0.05). Besides, HNF4A and UGT2B7 were both decreased in HepG2 cells and rats after treatment with emodin. In conclusion, emodin used long term or in high doses could inhibit the expression of HNF4A, thereby reducing the expression of UGT2B7 and causing hepatotoxicity.