Antagonism of B cell enhancer networks by STAT5 drives leukemia and poor patient survival.

The transcription factor STAT5 has a critical role in B cell acute lymphoblastic leukemia (B-ALL). How STAT5 mediates this effect is unclear. Here we found that activation of STAT5 worked together with defects in signaling components of the precursor to the B cell antigen receptor (pre-BCR), including defects in BLNK, BTK, PKCß, NF-κB1 and IKAROS, to initiate B-ALL. STAT5 antagonized the transcription factors NF-κB and IKAROS by opposing regulation of shared target genes. Super-enhancers showed enrichment for STAT5 binding and were associated with an opposing network of transcription factors, including PAX5, EBF1, PU.1, IRF4 and IKAROS. Patients with a high ratio of active STAT5 to NF-κB or IKAROS had more-aggressive disease. Our studies indicate that an imbalance of two opposing transcriptional programs drives B-ALL and suggest that restoring the balance of these pathways might inhibit B-ALL.

Scientific selection: A century of increasing crop varietal diversity in US wheat.

A prevalent and persistent biodiversity concern is that modern cropping systems lead to an erosion in crop genetic diversity. Although certain trait uniformity provides advantages in crop management and marketing, farmers facing risks from change in climate, pests, and markets are also incentivized to adopt new varieties to address complex and spatially variable genetics, environment, and crop management interactions to optimize crop performance. In this study, we applied phylogenetically blind and phylogenetically informed diversity metrics to reveal significant increases in both the spatial and temporal diversity of the US wheat crop over the past century. Contrary to commonly held perceptions on the negative impact of modern cropping systems on crop genetic diversity, our results demonstrated a win-win outcome where the widespread uptake of scientifically selected varieties increased both crop production and crop diversity.

Nodule-specific PLAT domain proteins are expanded in the Medicago lineage and required for nodulation.

Symbiotic nitrogen fixation in legumes is mediated by an interplay of signaling processes between plant hosts and rhizobial symbionts. In legumes, several secreted protein families have undergone expansions and play key roles in nodulation. Thus, identifying lineage-specific expansions (LSEs) of nodulation-associated genes can be a strategy to discover candidate gene families. Using bioinformatic tools, we identified 13 LSEs of nodulation-related secreted protein families, each unique to either Glycine, Arachis or Medicago lineages. In the Medicago lineage, nodule-specific Polycystin-1, Lipoxygenase, Alpha Toxin (PLAT) domain proteins (NPDs) expanded to five members. We examined NPD function using CRISPR/Cas9 multiplex genome editing to create Medicago truncatula NPD knockout lines, targeting one to five NPD genes. Mutant lines with differing combinations of NPD gene inactivations had progressively smaller nodules, earlier onset of nodule senescence, or ineffective nodules compared to the wild-type control. Double- and triple-knockout lines showed dissimilar nodulation phenotypes but coincided in upregulation of a DHHC-type zinc finger and an aspartyl protease gene, possible candidates for the observed disturbance of proper nodule function. By postulating that gene family expansions can be used to detect candidate genes, we identified a family of nodule-specific PLAT domain proteins and confirmed that they play a role in successful nodule formation.

The deubiquitinase USP9X suppresses pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDA) remains a lethal malignancy despite much progress concerning its molecular characterization. PDA tumours harbour four signature somatic mutations in addition to numerous lower frequency genetic events of uncertain significance. Here we use Sleeping Beauty (SB) transposon-mediated insertional mutagenesis in a mouse model of pancreatic ductal preneoplasia to identify genes that cooperate with oncogenic Kras(G12D) to accelerate tumorigenesis and promote progression. Our screen revealed new candidate genes for PDA and confirmed the importance of many genes and pathways previously implicated in human PDA. The most commonly mutated gene was the X-linked deubiquitinase Usp9x, which was inactivated in over 50% of the tumours. Although previous work had attributed a pro-survival role to USP9X in human neoplasia, we found instead that loss of Usp9x enhances transformation and protects pancreatic cancer cells from anoikis. Clinically, low USP9X protein and messenger RNA expression in PDA correlates with poor survival after surgery, and USP9X levels are inversely associated with metastatic burden in advanced disease. Furthermore, chromatin modulation with trichostatin A or 5-aza-2'-deoxycytidine elevates USP9X expression in human PDA cell lines, indicating a clinical approach for certain patients. The conditional deletion of Usp9x cooperated with Kras(G12D) to accelerate pancreatic tumorigenesis in mice, validating their genetic interaction. We propose that USP9X is a major tumour suppressor gene with prognostic and therapeutic relevance in PDA.

ODG: Omics database generator - a tool for generating, querying, and analyzing multi-omics comparative databases to facilitate biological understanding.

BACKGROUND: Rapid generation of omics data in recent years have resulted in vast amounts of disconnected datasets without systemic integration and knowledge building, while individual groups have made customized, annotated datasets available on the web with few ways to link them to in-lab datasets. With so many research groups generating their own data, the ability to relate it to the larger genomic and comparative genomic context is becoming increasingly crucial to make full use of the data. RESULTS: The Omics Database Generator (ODG) allows users to create customized databases that utilize published genomics data integrated with experimental data which can be queried using a flexible graph database. When provided with omics and experimental data, ODG will create a comparative, multi-dimensional graph database. ODG can import definitions and annotations from other sources such as InterProScan, the Gene Ontology, ENZYME, UniPathway, and others. This annotation data can be especially useful for studying new or understudied species for which transcripts have only been predicted, and rapidly give additional layers of annotation to predicted genes. In better studied species, ODG can perform syntenic annotation translations or rapidly identify characteristics of a set of genes or nucleotide locations, such as hits from an association study. ODG provides a web-based user-interface for configuring the data import and for querying the database. Queries can also be run from the command-line and the database can be queried directly through programming language hooks available for most languages. ODG supports most common genomic formats as well as generic, easy to use tab-separated value format for user-provided annotations. CONCLUSIONS: ODG is a user-friendly database generation and query tool that adapts to the supplied data to produce a comparative genomic database or multi-layered annotation database. ODG provides rapid comparative genomic annotation and is therefore particularly useful for non-model or understudied species. For species for which more data are available, ODG can be used to conduct complex multi-omics, pattern-matching queries.

Strategies for optimizing BioNano and Dovetail explored through a second reference quality assembly for the legume model, Medicago truncatula.

BACKGROUND: Third generation sequencing technologies, with sequencing reads in the tens- of kilo-bases, facilitate genome assembly by spanning ambiguous regions and improving continuity. This has been critical for plant genomes, which are difficult to assemble due to high repeat content, gene family expansions, segmental and tandem duplications, and polyploidy. Recently, high-throughput mapping and scaffolding strategies have further improved continuity. Together, these long-range technologies enable quality draft assemblies of complex genomes in a cost-effective and timely manner. RESULTS: Here, we present high quality genome assemblies of the model legume plant, Medicago truncatula (R108) using PacBio, Dovetail Chicago (hereafter, Dovetail) and BioNano technologies. To test these technologies for plant genome assembly, we generated five assemblies using all possible combinations and ordering of these three technologies in the R108 assembly. While the BioNano and Dovetail joins overlapped, they also showed complementary gains in continuity and join numbers. Both technologies spanned repetitive regions that PacBio alone was unable to bridge. Combining technologies, particularly Dovetail followed by BioNano, resulted in notable improvements compared to Dovetail or BioNano alone. A combination of PacBio, Dovetail, and BioNano was used to generate a high quality draft assembly of R108, a M. truncatula accession widely used in studies of functional genomics. As a test for the usefulness of the resulting genome sequence, the new R108 assembly was used to pinpoint breakpoints and characterize flanking sequence of a previously identified translocation between chromosomes 4 and 8, identifying more than 22.7 Mb of novel sequence not present in the earlier A17 reference assembly. CONCLUSIONS: Adding Dovetail followed by BioNano data yielded complementary improvements in continuity over the original PacBio assembly. This strategy proved efficient and cost-effective for developing a quality draft assembly compared to traditional reference assemblies.

Exploring structural variation and gene family architecture with De Novo assemblies of 15 Medicago genomes.

BACKGROUND: Previous studies exploring sequence variation in the model legume, Medicago truncatula, relied on mapping short reads to a single reference. However, read-mapping approaches are inadequate to examine large, diverse gene families or to probe variation in repeat-rich or highly divergent genome regions. De novo sequencing and assembly of M. truncatula genomes enables near-comprehensive discovery of structural variants (SVs), analysis of rapidly evolving gene families, and ultimately, construction of a pan-genome. RESULTS: Genome-wide synteny based on 15 de novo M. truncatula assemblies effectively detected different types of SVs indicating that as much as 22% of the genome is involved in large structural changes, altogether affecting 28% of gene models. A total of 63 million base pairs (Mbp) of novel sequence was discovered, expanding the reference genome space for Medicago by 16%. Pan-genome analysis revealed that 42% (180 Mbp) of genomic sequences is missing in one or more accession, while examination of de novo annotated genes identified 67% (50,700) of all ortholog groups as dispensable - estimates comparable to recent studies in rice, maize and soybean. Rapidly evolving gene families typically associated with biotic interactions and stress response were found to be enriched in the accession-specific gene pool. The nucleotide-binding site leucine-rich repeat (NBS-LRR) family, in particular, harbors the highest level of nucleotide diversity, large effect single nucleotide change, protein diversity, and presence/absence variation. However, the leucine-rich repeat (LRR) and heat shock gene families are disproportionately affected by large effect single nucleotide changes and even higher levels of copy number variation. CONCLUSIONS: Analysis of multiple M. truncatula genomes illustrates the value of de novo assemblies to discover and describe structural variation, something that is often under-estimated when using read-mapping approaches. Comparisons among the de novo assemblies also indicate that different large gene families differ in the architecture of their structural variation.

Hybrid assembly with long and short reads improves discovery of gene family expansions.

BACKGROUND: Long-read and short-read sequencing technologies offer competing advantages for eukaryotic genome sequencing projects. Combinations of both may be appropriate for surveys of within-species genomic variation. METHODS: We developed a hybrid assembly pipeline called "Alpaca" that can operate on 20X long-read coverage plus about 50X short-insert and 50X long-insert short-read coverage. To preclude collapse of tandem repeats, Alpaca relies on base-call-corrected long reads for contig formation. RESULTS: Compared to two other assembly protocols, Alpaca demonstrated the most reference agreement and repeat capture on the rice genome. On three accessions of the model legume Medicago truncatula, Alpaca generated the most agreement to a conspecific reference and predicted tandemly repeated genes absent from the other assemblies. CONCLUSION: Our results suggest Alpaca is a useful tool for investigating structural and copy number variation within de novo assemblies of sampled populations.

The Medicago genome provides insight into the evolution of rhizobial symbioses.

Legumes (Fabaceae or Leguminosae) are unique among cultivated plants for their ability to carry out endosymbiotic nitrogen fixation with rhizobial bacteria, a process that takes place in a specialized structure known as the nodule. Legumes belong to one of the two main groups of eurosids, the Fabidae, which includes most species capable of endosymbiotic nitrogen fixation. Legumes comprise several evolutionary lineages derived from a common ancestor 60 million years ago (Myr ago). Papilionoids are the largest clade, dating nearly to the origin of legumes and containing most cultivated species. Medicago truncatula is a long-established model for the study of legume biology. Here we describe the draft sequence of the M. truncatula euchromatin based on a recently completed BAC assembly supplemented with Illumina shotgun sequence, together capturing ∼94% of all M. truncatula genes. A whole-genome duplication (WGD) approximately 58 Myr ago had a major role in shaping the M. truncatula genome and thereby contributed to the evolution of endosymbiotic nitrogen fixation. Subsequent to the WGD, the M. truncatula genome experienced higher levels of rearrangement than two other sequenced legumes, Glycine max and Lotus japonicus. M. truncatula is a close relative of alfalfa (Medicago sativa), a widely cultivated crop with limited genomics tools and complex autotetraploid genetics. As such, the M. truncatula genome sequence provides significant opportunities to expand alfalfa's genomic toolbox.

TALEN-engineered AR gene rearrangements reveal endocrine uncoupling of androgen receptor in prostate cancer.

Androgen receptor (AR) target genes direct development and survival of the prostate epithelial lineage, including prostate cancer (PCa). Thus, endocrine therapies that inhibit the AR ligand-binding domain (LBD) are effective in treating PCa. AR transcriptional reactivation is central to resistance, as evidenced by the efficacy of AR retargeting in castration-resistant PCa (CRPC) with next-generation endocrine therapies abiraterone and enzalutamide. However, resistance to abiraterone and enzalutamide limits this efficacy in most men, and PCa remains the second-leading cause of male cancer deaths. Here we show that AR gene rearrangements in CRPC tissues underlie a completely androgen-independent, yet AR-dependent, resistance mechanism. We discovered intragenic AR gene rearrangements in CRPC tissues, which we modeled using transcription activator-like effector nuclease (TALEN)-mediated genome engineering. This modeling revealed that these AR gene rearrangements blocked full-length AR synthesis, but promoted expression of truncated AR variant proteins lacking the AR ligand-binding domain. Furthermore, these AR variant proteins maintained the constitutive activity of the AR transcriptional program and a CRPC growth phenotype independent of full-length AR or androgens. These findings demonstrate that AR gene rearrangements are a unique resistance mechanism by which AR transcriptional activity can be uncoupled from endocrine regulation in CRPC.

Patterns of divergence of a large family of nodule cysteine-rich peptides in accessions of Medicago truncatula.

The nodule cysteine-rich (NCR) groups of defensin-like (DEFL) genes are one of the largest gene families expressed in the nodules of some legume plants. They have only been observed in the inverted repeat loss clade (IRLC) of legumes, which includes the model legume Medicago truncatula. NCRs are reported to play an important role in plant-microbe interactions. To understand their diversity we analyzed their expression and sequence polymorphisms among four accessions of M. truncatula. A significant expression and nucleotide variation was observed among the genes. We then used 26 accessions to estimate the selection pressures shaping evolution among the accessions by calculating the nucleotide diversity at non-synonymous and synonymous sites in the coding region. The mature peptides of the orthologous NCRs had signatures of both purifying and diversifying selection pressures, unlike the seed DEFLs, which predominantly exhibited purifying selection. The expression, sequence variation and apparent diversifying selection in NCRs within the Medicago species indicates rapid and recent evolution, and suggests that this family of genes is actively evolving to adapt to different environments and is acquiring new functions.

Clonal Dissemination of Enterobacter cloacae Harboring blaKPC-3 in the Upper Midwestern United States.

Carbapenemase-producing, carbapenem-resistant Enterobacteriaceae, or CP-CRE, are an emerging threat to human and animal health, because they are resistant to many of the last-line antimicrobials available for disease treatment. Carbapenemase-producing Enterobacter cloacae harboring blaKPC-3 recently was reported in the upper midwestern United States and implicated in a hospital outbreak in Fargo, North Dakota (L. M. Kiedrowski, D. M. Guerrero, F. Perez, R. A. Viau, L. J. Rojas, M. F. Mojica, S. D. Rudin, A. M. Hujer, S. H. Marshall, and R. A. Bonomo, Emerg Infect Dis 20:1583-1585, 2014, http://dx.doi.org/10.3201/eid2009.140344). In early 2009, the Minnesota Department of Health began collecting and screening CP-CRE from patients throughout Minnesota. Here, we analyzed a retrospective group of CP-E. cloacae isolates (n = 34) collected between 2009 and 2013. Whole-genome sequencing and analysis revealed that 32 of the strains were clonal, belonging to the ST171 clonal complex and differing collectively by 211 single-nucleotide polymorphisms, and it revealed a dynamic clone under positive selection. The phylogeography of these strains suggests that this clone existed in eastern North Dakota and western Minnesota prior to 2009 and subsequently was identified in the Minneapolis and St. Paul metropolitan area. All strains harbored identical IncFIA-like plasmids conferring a CP-CRE phenotype and an additional IncX3 plasmid. In a single patient with multiple isolates submitted over several months, we found evidence that these plasmids had transferred from the E. cloacae clone to an Escherichia coli ST131 bacterium, rendering it as a CP-CRE. The spread of this clone throughout the upper midwestern United States is unprecedented for E. cloacae and highlights the importance of continued surveillance to identify such threats to human health.

Cross platform analysis of methylation, miRNA and stem cell gene expression data in germ cell tumors highlights characteristic differences by tumor histology.

BACKGROUND: Alterations in methylation patterns, miRNA expression, and stem cell protein expression occur in germ cell tumors (GCTs). Our goal is to integrate molecular data across platforms to identify molecular signatures in the three main histologic subtypes of Type I and Type II GCTs (yolk sac tumor (YST), germinoma, and teratoma). METHODS: We included 39 GCTs and 7 paired adjacent tissue samples in the current analysis. Molecular data available for analysis include DNA methylation data (Illumina GoldenGate Cancer Methylation Panel I), miRNA expression (NanoString nCounter miRNA platform), and stem cell factor expression (SABiosciences Human Embryonic Stem Cell Array). We evaluated the cross platform correlations of the data features using the Maximum Information Coefficient (MIC). RESULTS: In analyses of individual datasets, differences were observed by tumor histology. Germinomas had higher expression of transcription factors maintaining stemness, while YSTs had higher expression of cytokines, endoderm and endothelial markers. We also observed differences in miRNA expression, with miR-371-5p, miR-122, miR-302a, miR-302d, and miR-373 showing elevated expression in one or more histologic subtypes. Using the MIC, we identified correlations across the data features, including six major hubs with higher expression in YST (LEFTY1, LEFTY2, miR302b, miR302a, miR 126, and miR 122) compared with other GCT. CONCLUSIONS: While prognosis for GCTs is overall favorable, many patients experience resistance to chemotherapy, relapse and/or long term adverse health effects following treatment. Targeted therapies, based on integrated analyses of molecular tumor data such as that presented here, may provide a way to secure high cure rates while reducing unintended health consequences.

Evolution of pathogenicity and sexual reproduction in eight Candida genomes.

Candida species are the most common cause of opportunistic fungal infection worldwide. Here we report the genome sequences of six Candida species and compare these and related pathogens and non-pathogens. There are significant expansions of cell wall, secreted and transporter gene families in pathogenic species, suggesting adaptations associated with virulence. Large genomic tracts are homozygous in three diploid species, possibly resulting from recent recombination events. Surprisingly, key components of the mating and meiosis pathways are missing from several species. These include major differences at the mating-type loci (MTL); Lodderomyces elongisporus lacks MTL, and components of the a1/2 cell identity determinant were lost in other species, raising questions about how mating and cell types are controlled. Analysis of the CUG leucine-to-serine genetic-code change reveals that 99% of ancestral CUG codons were erased and new ones arose elsewhere. Lastly, we revise the Candida albicans gene catalogue, identifying many new genes.

A Sleeping Beauty mutagenesis screen reveals a tumor suppressor role for Ncoa2/Src-2 in liver cancer.

The Sleeping Beauty (SB) transposon mutagenesis system is a powerful tool that facilitates the discovery of mutations that accelerate tumorigenesis. In this study, we sought to identify mutations that cooperate with MYC, one of the most commonly dysregulated genes in human malignancy. We performed a forward genetic screen with a mouse model of MYC-induced liver cancer using SB-mediated mutagenesis. We sequenced insertions in 63 liver tumor nodules and identified at least 16 genes/loci that contribute to accelerated tumor development. RNAi-mediated knockdown in a liver progenitor cell line further validate three of these genes, Ncoa2/Src-2, Zfx, and Dtnb, as tumor suppressors in liver cancer. Moreover, deletion of Ncoa2/Src-2 in mice predisposes to diethylnitrosamine-induced liver tumorigenesis. These findings reveal genes and pathways that functionally restrain MYC-mediated liver tumorigenesis and therefore may provide targets for cancer therapy.

SHEAR: sample heterogeneity estimation and assembly by reference.

BACKGROUND: Personal genome assembly is a critical process when studying tumor genomes and other highly divergent sequences. The accuracy of downstream analyses, such as RNA-seq and ChIP-seq, can be greatly enhanced by using personal genomic sequences rather than standard references. Unfortunately, reads sequenced from these types of samples often have a heterogeneous mix of various subpopulations with different variants, making assembly extremely difficult using existing assembly tools. To address these challenges, we developed SHEAR (Sample Heterogeneity Estimation and Assembly by Reference; http://vk.cs.umn.edu/SHEAR), a tool that predicts SVs, accounts for heterogeneous variants by estimating their representative percentages, and generates personal genomic sequences to be used for downstream analysis. RESULTS: By making use of structural variant detection algorithms, SHEAR offers improved performance in the form of a stronger ability to handle difficult structural variant types and better computational efficiency. We compare against the lead competing approach using a variety of simulated scenarios as well as real tumor cell line data with known heterogeneous variants. SHEAR is shown to successfully estimate heterogeneity percentages in both cases, and demonstrates an improved efficiency and better ability to handle tandem duplications. CONCLUSION: SHEAR allows for accurate and efficient SV detection and personal genomic sequence generation. It is also able to account for heterogeneous sequencing samples, such as from tumor tissue, by estimating the subpopulation percentage for each heterogeneous variant.

MMuFLR: missense mutation and frameshift location reporter.

MOTIVATION: Cancer researchers seeking immunotherapy targets in cancer cells need tools to locate highly expressed proteins unique to cancer cells. Missense mutation and frameshift location reporter (MMuFLR), a Galaxy-based workflow, analyzes next-generation sequencing paired read RNA-seq output to reliably identify small frameshift mutations and missense mutations in highly expressed protein-coding genes. MMuFLR ignores known SNPs, low quality reads and poly-A/T sequences. For each frameshift and missense mutation identified, MMuFLR provides the location and sequence of the amino acid substitutions in the novel protein candidates for direct input into epitope evaluation tools. AVAILABILITY: http://toolshed.g2.bx.psu.edu/ CONTACT: rath0096@umn.edu or johns198@umn.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Sex bias occurrence of hepatocellular carcinoma in Poly7 molecular subclass is associated with EGFR.

UNLABELLED: Hepatocellular carcinoma (HCC) is one of the deadliest solid cancers and is the third leading cause of cancer-related death. There is a universal estimated male/female ratio of 2.5, but the reason for this is not well understood. The Sleeping Beauty (SB) transposon system was used to elucidate candidate oncogenic drivers of HCC in a forward genetics screening approach. Sex bias occurrence was conserved in our model, with male experimental mice developing liver tumors at reduced latency and higher tumor penetrance. In parallel, we explored sex differences regarding genomic aberrations in 235 HCC patients. Liver cancer candidate genes were identified from both sexes and genotypes. Interestingly, transposon insertions in the epidermal growth factor receptor (Egfr) gene were common in SB-induced liver tumors from male mice (10/10, 100%) but infrequent in female mice (2/9, 22%). Human single-nucleotide polymorphism data confirmed that polysomy of chromosome 7, locus of EGFR, was more frequent in males (26/62, 41%) than females (2/27, 7%) (P = 0.001). Gene expression-based Poly7 subclass patients were predominantly male (9/9) compared with 67% males (55/82) in other HCC subclasses (P = 0.02), and this subclass was accompanied by EGFR overexpression (P < 0.001). CONCLUSION: Sex bias occurrence of HCC associated with EGFR was confirmed in experimental animals using the SB transposon system in a reverse genetic approach. This study provides evidence for the role of EGFR in sex bias occurrences of liver cancer and as the driver mutational gene in the Poly7 molecular subclass of human HCC.

An insertional mutagenesis screen identifies genes that cooperate with Mll-AF9 in a murine leukemogenesis model.

Patients with a t(9;11) translocation (MLL-AF9) develop acute myeloid leukemia (AML), and while in mice the expression of this fusion oncogene also results in the development of myeloid leukemia, it is with long latency. To identify mutations that cooperate with Mll-AF9, we infected neonatal wild-type (WT) or Mll-AF9 mice with a murine leukemia virus (MuLV). MuLV-infected Mll-AF9 mice succumbed to disease significantly faster than controls presenting predominantly with myeloid leukemia while infected WT animals developed predominantly lymphoid leukemia. We identified 88 candidate cancer genes near common sites of proviral insertion. Analysis of transcript levels revealed significantly elevated expression of Mn1, and a trend toward increased expression of Bcl11a and Fosb in Mll-AF9 murine leukemia samples with proviral insertions proximal to these genes. Accordingly, FOSB and BCL11A were also overexpressed in human AML harboring MLL gene translocations. FOSB was revealed to be essential for growth in mouse and human myeloid leukemia cells using shRNA lentiviral vectors in vitro. Importantly, MN1 cooperated with Mll-AF9 in leukemogenesis in an in vivo BM viral transduction and transplantation assay. Together, our data identified genes that define transcription factor networks and important genetic pathways acting during progression of leukemia induced by MLL fusion oncogenes.