State-of-the-art in transposable element modulation affected by drugs in malignant prostatic cancer cells.

Over recent years, the investigation of transposable elements (TEs) has granted researchers a deeper comprehension of their characteristics and functions, particularly regarding their significance in the mechanisms contributing to cancer development. This manuscript focuses on prostate carcinoma cell lines and offers a comprehensive review intended to scrutinize the associations and interactions between TEs and genes, as well as their response to treatment using various chemical drugs, emphasizing their involvement in cancer progression. We assembled a compendium of articles retrieved from the PubMed database to construct networks demonstrating correlations with genes and pharmaceuticals. In doing so, we linked the transposition of certain TE types to the expression of specific transcripts directly implicated in carcinogenesis. Additionally, we underline that treatment employing different drugs revealed unique patterns of TE reactivation. Our hypothesis gathers the current understanding and guides research toward evidence-based investigations, emphasizing the association between antiviral drugs, chemotherapy, and the reduced expression of TEs in patients affected by prostate cancer.

TumFlow: An AI Model for Predicting New Anticancer Molecules.

Melanoma is the fifth most common cancer in the United States. Conventional drug discovery methods are inherently time-consuming and costly, which imposes significant limitations. However, the advent of Artificial Intelligence (AI) has opened up new possibilities for simulating and evaluating numerous drug candidates, thereby mitigating the requisite time and resources. In this context, normalizing flow models by employing machine learning techniques to create new molecular structures holds promise for accelerating the discovery of effective anticancer therapies. This manuscript introduces TumFlow, a novel AI model designed to generate new molecular entities with potential therapeutic value in cancer treatment. It has been trained on the NCI-60 dataset, encompassing thousands of molecules tested across 60 tumour cell lines, with an emphasis on the melanoma SK-MEL-28 cell line. The model successfully generated new molecules with predicted improved efficacy in inhibiting tumour growth while being synthetically feasible. This represents a significant advancement over conventional generative models, which often produce molecules that are challenging or impossible to synthesize. Furthermore, TumFlow has also been utilized to optimize molecules known for their efficacy in clinical melanoma treatments. This led to the creation of novel molecules with a predicted enhanced likelihood of effectiveness against melanoma, currently undocumented on PubChem.

DNA sequence symmetries from randomness: the origin of the Chargaff's second parity rule.

Most living organisms rely on double-stranded DNA (dsDNA) to store their genetic information and perpetuate themselves. This biological information has been considered as the main target of evolution. However, here we show that symmetries and patterns in the dsDNA sequence can emerge from the physical peculiarities of the dsDNA molecule itself and the maximum entropy principle alone, rather than from biological or environmental evolutionary pressure. The randomness justifies the human codon biases and context-dependent mutation patterns in human populations. Thus, the DNA 'exceptional symmetries,' emerged from the randomness, have to be taken into account when looking for the DNA encoded information. Our results suggest that the double helix energy constraints and, more generally, the physical properties of the dsDNA are the hard drivers of the overall DNA sequence architecture, whereas the selective biological processes act as soft drivers, which only under extraordinary circumstances overtake the overall entropy content of the genome.

Abrogation of esophageal carcinoma development in miR-31 knockout rats.

MicroRNA-31 (miR-31) is overexpressed in esophageal squamous cell carcinoma (ESCC), a deadly disease associated with dietary Zn deficiency and inflammation. In a Zn deficiency-promoted rat ESCC model with miR-31 up-regulation, cancer-associated inflammation, and a high ESCC burden following N-nitrosomethylbenzylamine (NMBA) exposure, systemic antimiR-31 delivery reduced ESCC incidence from 85 to 45% (P = 0.038) and miR-31 gene knockout abrogated development of ESCC (P = 1 × 10-6). Transcriptomics, genome sequencing, and metabolomics analyses in these Zn-deficient rats revealed the molecular basis of ESCC abrogation by miR-31 knockout. Our identification of EGLN3, a known negative regulator of nuclear factor κB (NF-κB), as a direct target of miR-31 establishes a functional link between oncomiR-31, tumor suppressor target EGLN3, and up-regulated NF-κB-controlled inflammation signaling. Interaction among oncogenic miR-31, EGLN3 down-regulation, and inflammation was also documented in human ESCCs. miR-31 deletion resulted in suppression of miR-31-associated EGLN3/NF-κB-controlled inflammatory pathways. ESCC-free, Zn-deficient miR-31-/- rat esophagus displayed no genome instability and limited metabolic activity changes vs. the pronounced mutational burden and ESCC-associated metabolic changes of Zn-deficient wild-type rats. These results provide conclusive evidence that miR-31 expression is necessary for ESCC development.

Human-like hyperplastic prostate with low ZIP1 induced solely by Zn deficiency in rats.

Prostate cancer is a leading cause of cancer death in men over 50 years of age, and there is a characteristic marked decrease in Zn content in the malignant prostate cells. The cause and consequences of this loss have thus far been unknown. We found that in middle-aged rats a Zn-deficient diet reduces prostatic Zn levels (P = 0.025), increases cellular proliferation, and induces an inflammatory phenotype with COX-2 overexpression. This hyperplastic/inflammatory prostate has a human prostate cancer-like microRNA profile, with up-regulation of the Zn-homeostasis-regulating miR-183-96-182 cluster (fold change = 1.41-2.38; P = 0.029-0.0003) and down-regulation of the Zn importer ZIP1 (target of miR-182), leading to a reduction of prostatic Zn. This inverse relationship between miR-182 and ZIP1 also occurs in human prostate cancer tissue, which is known for Zn loss. The discovery that the Zn-depleted middle-aged rat prostate has a metabolic phenotype resembling that of human prostate cancer, with a 10-fold down-regulation of citric acid (P = 0.0003), links citrate reduction directly to prostatic Zn loss, providing the underlying mechanism linking dietary Zn deficiency with miR-183-96-182 overexpression, ZIP1 down-regulation, prostatic Zn loss, and the resultant citrate down-regulation, changes mimicking features of human prostate cancer. Thus, dietary Zn deficiency during rat middle age produces changes that mimic those of human prostate carcinoma and may increase the risk for prostate cancer, supporting the need for assessment of Zn supplementation in its prevention.

GMIEC: a shiny application for the identification of gene-targeted drugs for precision medicine.

BACKGROUND: Precision medicine is a medical approach that takes into account individual genetic variability and often requires Next Generation Sequencing data in order to predict new treatments. Here we present GMIEC, Genomic Modules Identification et Characterization for genomics medicine, an application that is able to identify specific drugs at the level of single patient integrating multi-omics data such as RNA-sequencing, copy-number variation, methylation, Chromatin Immuno-Precipitation and Exome/Whole Genome sequencing. It is also possible to include clinical data related to each patient. GMIEC has been developed as a web-based R-Shiny platform and gives as output a table easy to use and explore. RESULTS: We present GMIEC, a Shiny application for genomics medicine. The tool allows the users the integration of two or more multiple omics datasets (e.g. gene-expression, copy-number), at sample level, to identify groups of genes that share common genomic and corresponding drugs. We demonstrate the characteristics of our application by using it to analyze a prostate cancer data set. CONCLUSIONS: GMIEC provides a simple interface for genomics medicine. GMIEC was develop with Shiny to provide an application that does not require advanced programming skills. GMIEC consists of three sub-application for the analysis (GMIEC-AN), the visualization (GMIEC-VIS) and the exploration of results (GMIEC-RES). GMIEC is an open source software and is available at https://github.com/guidmt/GMIEC-shiny.

Correction to: Transposable Elements Activity is Positively Related to Rate of Speciation in Mammals.

The original version of the article unfortunately contained tagging error in Given and Surname of all the authors.

Transposable Elements Activity is Positively Related to Rate of Speciation in Mammals.

Transposable elements (TEs) play an essential role in shaping eukaryotic genomes and generating variability. Speciation and TE activity bursts could be strongly related in mammals, in which simple gradualistic models of differentiation do not account for the currently observed species variability. In order to test this hypothesis, we designed two parameters: the Density of insertion (DI) and the Relative rate of speciation (RRS). DI is the ratio between the number of TE insertions in a genome and its size, whereas the RRS is a conditional parameter designed to identify potential speciation bursts. Thus, by analyzing TE insertions in mammals, we defined the genomes as "hot" (high DI) and "cold" (low DI). Then, comparing TE activity among 29 taxonomical families of the whole Mammalia class, 16 intra-order pairs of mammalian species, and four superorders of Eutheria, we showed that taxa with high rates of speciation are associated with "hot" genomes, whereas taxa with low ones are associated with "cold" genomes. These results suggest a remarkable correlation between TE activity and speciation, also being consistent with patterns describing variable rates of differentiation and accounting for the different time frames of the speciation bursts.

Generation of human memory stem T cells after haploidentical T-replete hematopoietic stem cell transplantation.

Memory stem T cells (TSCM) have been proposed as key determinants of immunologic memory. However, their exact contribution to a mounting immune response, as well as the mechanisms and timing of their in vivo generation, are poorly understood. We longitudinally tracked TSCM dynamics in patients undergoing haploidentical hematopoietic stem cell transplantation (HSCT), thereby providing novel hints on the contribution of this subset to posttransplant immune reconstitution in humans. We found that donor-derived TSCM are highly enriched early after HSCT. We showed at the antigen-specific and clonal level that TSCM lymphocytes can differentiate directly from naive precursors infused within the graft and that the extent of TSCM generation might correlate with interleukin 7 serum levels. In vivo fate mapping through T-cell receptor sequencing allowed defining the in vivo differentiation landscapes of human naive T cells, supporting the notion that progenies of single naive cells embrace disparate fates in vivo and highlighting TSCM as relevant novel players in the diversification of immunological memory after allogeneic HSCT.

Estrogen mediated-activation of miR-191/425 cluster modulates tumorigenicity of breast cancer cells depending on estrogen receptor status.

MicroRNAs (miRNAs), single-stranded non-coding RNAs, influence myriad biological processes that can contribute to cancer. Although tumor-suppressive and oncogenic functions have been characterized for some miRNAs, the majority of microRNAs have not been investigated for their ability to promote and modulate tumorigenesis. Here, we established that the miR-191/425 cluster is transcriptionally dependent on the host gene, DALRD3, and that the hormone 17ß-estradiol (estrogen or E2) controls expression of both miR-191/425 and DALRD3. MiR-191/425 locus characterization revealed that the recruitment of estrogen receptor α (ERα) to the regulatory region of the miR-191/425-DALRD3 unit resulted in the accumulation of miR-191 and miR-425 and subsequent decrease in DALRD3 expression levels. We demonstrated that miR-191 protects ERα positive breast cancer cells from hormone starvation-induced apoptosis through the suppression of tumor-suppressor EGR1. Furthermore, enforced expression of the miR-191/425 cluster in aggressive breast cancer cells altered global gene expression profiles and enabled us to identify important tumor promoting genes, including SATB1, CCND2, and FSCN1, as targets of miR-191 and miR-425. Finally, in vitro and in vivo experiments demonstrated that miR-191 and miR-425 reduced proliferation, impaired tumorigenesis and metastasis, and increased expression of epithelial markers in aggressive breast cancer cells. Our data provide compelling evidence for the transcriptional regulation of the miR-191/425 cluster and for its context-specific biological determinants in breast cancers. Importantly, we demonstrated that the miR-191/425 cluster, by reducing the expression of an extensive network of genes, has a fundamental impact on cancer initiation and progression of breast cancer cells.

The role of p19 and p21 H-Ras proteins and mutants in miRNA expression in cancer and a Costello syndrome cell model.

BACKGROUND: P19 H-Ras, a second product derived from the H-Ras gene by alternative splicing, induces a G1/S phase delay, thereby maintaining cells in a reversible quiescence state. When P21 H-Ras is mutated in tumour cells, the alternative protein P19 H-Ras is also mutated. The H-Ras mutation Q61L is frequently detected in different tumours, which acts as constitutive activator of Ras functions and is considered to be a strong activating mutant. Additionally, a rare congenital disorder named Costello Syndrome, is described as a H-Ras disorder in children, mainly due to mutation G12S in p19 and p21 H-Ras proteins, which is present in 90 % of the Costello Syndrome patients. Our aim is to better understand the role of p19 and p21 H-Ras proteins in the cancer and Costello Syndrome development, concerning the miRNAs expression. METHODS: Total miRNAs expression regulated by H-Ras proteins were first analyzed in human miRNA microarrays assays. Previously selected miRNAs, were further analyzed in developed cell lines containing H-Ras protein mutants, that included the G12S Costello Syndrome mutant, with PCR Real-Time Taq Man miRNA Assays primers. RESULTS: This study describes how p19 affects the RNA world and shows that: i) miR-342, miR-206, miR-330, miR-138 and miR-99b are upregulated by p19 but not by p19W164A mutant; ii) anti-miR-206 can restore the G2 phase in the presence of p19; iii) p19 and p21Q61L regulate their own alternative splicing; iv) miR-206 and miR-138 are differentially regulated by p19 and p21 H-Ras and v) P19G12S Costello mutants show a clear upregulation of miR-374, miR-126, miR-342, miR-330, miR-335 and let-7. CONCLUSIONS: These results allow us to conclude that the H-Ras G12S mutation plays an important role in miRNA expression and open up a new line of study to understand the consequences of this mutation on Costello syndrome. Furthermore, they suggest that oncogenes may have a sufficiently important impact on miRNA expression to promote the development of numerous cancers.

Ultraconserved regions encoding ncRNAs are altered in human leukemias and carcinomas.

Noncoding RNA (ncRNA) transcripts are thought to be involved in human tumorigenesis. We report that a large fraction of genomic ultraconserved regions (UCRs) encode a particular set of ncRNAs whose expression is altered in human cancers. Genome-wide profiling revealed that UCRs have distinct signatures in human leukemias and carcinomas. UCRs are frequently located at fragile sites and genomic regions involved in cancers. We identified certain UCRs whose expression may be regulated by microRNAs abnormally expressed in human chronic lymphocytic leukemia, and we proved that the inhibition of an overexpressed UCR induces apoptosis in colon cancer cells. Our findings argue that ncRNAs and interaction between noncoding genes are involved in tumorigenesis to a greater extent than previously thought.

Transcriptomics Studies Reveal Functions of Transglutaminase 2 in Breast Cancer Cells Using Membrane Permeable and Impermeable Inhibitors.

Transglutaminase 2 (TG2) performs many functions both under physiological and pathological conditions. In cancer, its expression is associated with aggressiveness, propensity to epithelial-mesenchymal transition, and metastasis. Since TG2 performs key functions both outside and inside the cell, using inhibitors with different membrane permeability we analyzed the changes in the transcriptome induced in two triple-negative cell lines (MDA-MB-436 and MDA-MB-231) with aggressive features. By characterizing pathways and gene networks, we were able to define the effects of TG2 inhibitors (AA9, membrane-permeable, and NCEG2, impermeable) in relation to the roles of the enzyme in the intra- and extracellular space within the context of breast cancer. The deregulated genes revealed p53 and integrin signaling to be the common pathways with some genes showing opposite changes in expression. In MDA-MB-436, AA9 induced apoptosis, modulated cadherin, Wnt, gastrin and cholecystokinin receptors (CCKR) mediated signaling, with RHOB and GNG2 playing significant roles, and affected the Warburg effect by decreasing glycolytic enzymes. In MDA-MB-231 cells, AA9 strongly impacted HIF-mediated hypoxia, including AKT and mTOR pathway. These effects suggest an anti-tumor activity by blocking intracellular TG2 functions. Conversely, the use of NCEG2 stimulated the expression of ATP synthase and proteins involved in DNA replication, indicating a potential promotion of cell proliferation through inhibition of extracellular TG2. To effectively utilize these molecules as an anti-tumor strategy, an appropriate delivery system should be evaluated to target specific functions and avoid adverse effects. Additionally, considering combinations with other pathway modulators is crucial.

Reprogramming of miRNA networks in cancer and leukemia.

We studied miRNA profiles in 4419 human samples (3312 neoplastic, 1107 nonmalignant), corresponding to 50 normal tissues and 51 cancer types. The complexity of our database enabled us to perform a detailed analysis of microRNA (miRNA) activities. We inferred genetic networks from miRNA expression in normal tissues and cancer. We also built, for the first time, specialized miRNA networks for solid tumors and leukemias. Nonmalignant tissues and cancer networks displayed a change in hubs, the most connected miRNAs. hsa-miR-103/106 were downgraded in cancer, whereas hsa-miR-30 became most prominent. Cancer networks appeared as built from disjointed subnetworks, as opposed to normal tissues. A comparison of these nets allowed us to identify key miRNA cliques in cancer. We also investigated miRNA copy number alterations in 744 cancer samples, at a resolution of 150 kb. Members of miRNA families should be similarly deleted or amplified, since they repress the same cellular targets and are thus expected to have similar impacts on oncogenesis. We correctly identified hsa-miR-17/92 family as amplified and the hsa-miR-143/145 cluster as deleted. Other miRNAs, such as hsa-miR-30 and hsa-miR-204, were found to be physically altered at the DNA copy number level as well. By combining differential expression, genetic networks, and DNA copy number alterations, we confirmed, or discovered, miRNAs with comprehensive roles in cancer. Finally, we experimentally validated the miRNA network with acute lymphocytic leukemia originated in Mir155 transgenic mice. Most of miRNAs deregulated in these transgenic mice were located close to hsa-miR-155 in the cancer network.

miR-181b is a biomarker of disease progression in chronic lymphocytic leukemia.

MicroRNAs play a crucial role in chronic lymphocytic leukemia. We investigated whether microRNAs can discriminate patients with a progressive disease from patients with a stable disease. We analyzed microRNA expression on leukemic cells isolated from 358 sequential samples of 114 patients with either stable or progressive disease. We found that during the course of the disease the expression values of miR-181b, the most dysregulated microRNA, decreased in samples of patients with a progressive (P < .001, training and validation sets) but not in samples of patients with a stable disease (P = .3, training set; P = .2, validation set) over time. A drop of ≥ 50% between sequential samples and/or a miR-181b value ≤ 0.005 at the starting time point were significant to differentiate progressive from stable disease (P = .004, training set; P < .001, validation set). These parameters were associated with high risk of requiring treatment (risk ratio, 5.8; 95% confidence interval, 2.5-14.9). We also observed that miR-181b targets Mcl-1 protein and that the decrease of its expression inversely correlated with increased protein levels of MCL1 and BCL2 target genes. We conclude that parameters defined on the basis of the miR-181b expression values specify disease progression in chronic lymphocytic leukemia and are associated with clinical outcome.

Comparative analysis of bats and rodents' genomes suggests a relation between non-LTR retrotransposons, cancer incidence, and ageing.

The presence in nature of species showing drastic differences in lifespan and cancer incidence has recently increased the interest of the scientific community. In particular, the adaptations and the genomic features underlying the evolution of cancer-resistant and long-lived organisms have recently focused on transposable elements (TEs). In this study, we compared the content and dynamics of TE activity in the genomes of four rodent and six bat species exhibiting different lifespans and cancer susceptibility. Mouse, rat, and guinea pig genomes (short-lived and cancer-prone organisms) were compared with that of naked mole rat (Heterocephalus glaber) which is a cancer-resistant organism and the rodent with the longest lifespan. The long-lived bats of the genera Myotis, Rhinolophus, Pteropus and Rousettus were instead compared with Molossus molossus, which is one of the organisms with the shortest lifespan among the order Chiroptera. Despite previous hypotheses stating a substantial tolerance of TEs in bats, we found that long-lived bats and the naked mole rat share a marked decrease of non-LTR retrotransposons (LINEs and SINEs) accumulation in recent evolutionary times.

Uncovering Patterns in Dairy Cow Behaviour: A Deep Learning Approach with Tri-Axial Accelerometer Data.

The accurate detection of behavioural changes represents a promising method of detecting the early onset of disease in dairy cows. This study assessed the performance of deep learning (DL) in classifying dairy cows' behaviour from accelerometry data acquired by single sensors on the cows' left flanks and compared the results with those obtained through classical machine learning (ML) from the same raw data. Twelve cows with a tri-axial accelerometer were observed for 136 ± 29 min each to detect five main behaviours: standing still, moving, feeding, ruminating and resting. For each 8 s time interval, 15 metrics were calculated, obtaining a dataset of 211,720 observation units and 15 columns. The entire dataset was randomly split into training (80%) and testing (20%) datasets. The DL accuracy, precision and sensitivity/recall were calculated and compared with the performance of classical ML models. The best predictive model was an 8-layer convolutional neural network (CNN) with an overall accuracy and F1 score equal to 0.96. The precision, sensitivity/recall and F1 score of single behaviours had the following ranges: 0.93-0.99. The CNN outperformed all the classical ML algorithms. The CNN used to monitor the cows' conditions showed an overall high performance in successfully predicting multiple behaviours using a single accelerometer.

Potential of Full-Fat Silkworm-Based Diets for Laying Quails: Performance and Egg Physical Quality.

The present research was conducted to determine the optimal inclusion level of full-fat silkworm chrysalis meal (SWM) into laying quails' diets, focusing on performance traits and egg physical quality. A total of 240 31-day-old female Japanese quails were randomly assigned to four dietary groups (12 replicates/treatment; 5 quails/replicate); quails were initially fed a standard commercial diet for pullets until 63 days of age. When oviposition started, the experimental groups received the following diets: a conventional corn and soybean-based diet (control diet-C) and three other diets, including 4%, 8%, or 12% of full-fat SWM (SWM4, SWM8, SWM12, respectively). Experimental diets were provided until quails reached 119 days of age. Birds displayed satisfactory productive performance throughout the trial. SWM12 and SWM8 had higher (p < 0.001) egg production but also a higher feed conversion ratio compared to C. At the end of the trial, the eggs edible portion increased, and shell weight decreased with increasing the SWM dietary inclusion level (p < 0.001). At the same time, SWM12 displayed an increase in albumen pH (p < 0.05), even though in the normal range for quail egg. Overall, full-fat SWM (up to 12%) can be considered a promising feed ingredient for laying quails, although higher inclusion levels (>8%) require special attention because SWM also contains anti-nutritional factors.

Metabolic characterisation of transglutaminase 2 inhibitor effects in breast cancer cell lines.

Transglutaminase 2 (TG2), which mediates post-translational modifications of multiple intracellular enzymes, is involved in the pathogenesis and progression of cancer. We used 1 H-NMR metabolomics to study the effects of AA9, a novel TG2 inhibitor, on two breast cancer cell lines with distinct phenotypes, MCF-7 and MDA-MB-231. AA9 can promote apoptosis in both cell lines, but it is particularly effective in MD-MB-231, inhibiting transamidation reactions and decreasing cell migration and invasiveness. This metabolomics study provides evidence of a major effect of AA9 on MDA-MB-231 cells, impacting glutamate and aspartate metabolism, rather than on MCF-7 cells, characterised by choline and O-phosphocholine decrease. Interestingly, AA9 treatment induces myo-inositol alteration in both cell lines, indicating action on phosphatidylinositol metabolism, likely modulated by the G protein activity of TG2 on phospholipase C. Considering the metabolic deregulations that characterise various breast cancer subtypes, the existence of a metabolic pathway affected by AA9 further points to TG2 as a promising hot spot. The metabolomics approach provides a powerful tool to monitor the effectiveness of inhibitors and better understand the role of TG2 in cancer.

Dysregulation of miR-31 and miR-21 induced by zinc deficiency promotes esophageal cancer.

Zinc deficiency (ZD) increases the risk of esophageal squamous cell carcinoma (ESCC). In a rat model, chronic ZD induces an inflammatory gene signature that fuels ESCC development. microRNAs regulate gene expression and are aberrantly expressed in cancers. Here we investigated whether chronic ZD (23 weeks) also induces a protumorigenic microRNA signature. Using the nanoString technology, we evaluated microRNA profiles in ZD esophagus and six additional tissues (skin, lung, pancreas, liver, prostate and peripheral blood mononuclear cells [PBMC]). ZD caused overexpression of inflammation genes and altered microRNA expression across all tissues analyzed, predictive of disease development. Importantly, the inflammatory ZD esophagus had a distinct microRNA signature resembling human ESCC or tongue SCC miRNAomes with miR-31 and miR-21 as the top-up-regulated species. Circulating miR-31 was also the top-up-regulated species in PBMCs. In ZD esophagus and tongue, oncogenic miR-31 and miR-21 overexpression was accompanied by down-regulation of their respective tumor-suppressor targets PPP2R2A and PDCD4. Importantly, esophageal miR-31 and miR-21 levels were directly associated with the appearance of ESCC in ZD rats, as compared with their cancer-free Zn-sufficient or Zn-replenished counterparts. In situ hybridization analysis in rat and human tongue SCCs localized miR-31 to tumor cells and miR-21 to stromal cells. In regressing tongue SCCs from Zn-supplemented rats, miR-31 and miR-21 expression was concomitantly reduced, establishing their responsiveness to Zn therapy. A search for putative microRNA targets revealed a bias toward genes in inflammatory pathways. Our finding that ZD causes miR-31 and miR-21 dysregulation associated with inflammation provides insight into mechanisms whereby ZD promotes ESCC.