Differential alternative polyadenylation response to high-fat diet between polygenic obese and healthy lean mice.

Obesity's complex etiology due to the interplay of environment and genetics makes it a more challenging research and health problem. Some of the contributing genetic factors that have not yet been examined in detail entail mRNA polyadenylation (PA). Genes with multiple PA sites express mRNA isoforms differing in coding sequence or 3'UTR through alternative polyadenylation (APA). Alterations in PA have been associated with various diseases; however, its contribution to obesity is not well-researched. Following an 11-week high-fat diet, the APA sites in the hypothalamus of two unique mouse models for polygenic obesity (Fat line) and healthy leanness (Lean line) were determined using whole transcriptome termini site sequencing (WTTS-seq). We found 17 genes of interest with differentially expressed APA (DE-APA) isoforms, among which seven were previously associated with obesity or obesity-related traits (Pdxdc1, Smyd3, Rpl14, Copg1, Pcna, Ric3, Stx3) but have not yet been studied in the context of APA. The remaining ten genes (Ccdc25, Dtd2, Gm14403, Hlf, Lyrm7, Mrpl3, Pisd-ps3, Sbsn, Slx1b, Spon1) represent novel candidates associated with obesity/adiposity due to variability brought about by differential usage of APA sites. Our results provide insights into the relationship between PA and the hypothalamus in the context of obesity, by being the first study of DE-APA sites and DE-APA isoforms in these mouse models. Future studies are needed further to explore the role of APA isoforms in polygenic obesity by expanding the scope of research to other metabolically important tissues (such as liver and adipose tissues) and investigating the potential for targeting PA as a therapeutic strategy for obesity management.

Genome-wide screening for genetic variants in polyadenylation signal (PAS) sites in mouse selection lines for fatness and leanness.

Alternative polyadenylation (APA) determines mRNA stability, localisation, translation and protein function. Several diseases, including obesity, have been linked to APA. Studies have shown that single nucleotide polymorphisms in polyadenylation signals (PAS-SNPs) can influence APA and affect phenotype and disease susceptibility. However, these studies focussed on associations between single PAS-SNP alleles with very large effects and phenotype. Therefore, we performed a genome-wide screening for PAS-SNPs in the polygenic mouse selection lines for fatness and leanness by whole-genome sequencing. The genetic variants identified in the two lines were overlapped with locations of PAS sites obtained from the PolyASite 2.0 database. Expression data for selected genes were extracted from the microarray expression experiment performed on multiple tissue samples. In total, 682 PAS-SNPs were identified within 583 genes involved in various biological processes, including transport, protein modifications and degradation, cell adhesion and immune response. Moreover, 63 of the 583 orthologous genes in human have been previously associated with human diseases, such as nervous system and physical disorders, and immune, endocrine, and metabolic diseases. In both lines, PAS-SNPs have also been identified in genes broadly involved in APA, such as Polr2c, Eif3e and Ints11. Five PAS-SNPs within 5 genes (Car, Col4a1, Itga7, Lat, Nmnat1) were prioritised as potential functional variants and could contribute to the phenotypic disparity between the two selection lines. The developed PAS-SNPs catalogue presents a key resource for planning functional studies to uncover the role of PAS-SNPs in APA, disease susceptibility and fat deposition.

CaNDis: a web server for investigation of causal relationships between diseases, drugs and drug targets.

MOTIVATION: Causal biological interaction networks represent cellular regulatory pathways. Their fusion with other biological data enables insights into disease mechanisms and novel opportunities for drug discovery. RESULTS: We developed Causal Network of Diseases (CaNDis), a web server for the exploration of a human causal interaction network, which we expanded with data on diseases and FDA-approved drugs, on the basis of which we constructed a disease-disease network in which the links represent the similarity between diseases. We show how CaNDis can be used to identify candidate genes with known and novel roles in disease co-occurrence and drug-drug interactions. AVAILABILITYAND IMPLEMENTATION: CaNDis is freely available to academic users at http://candis.ijs.si and http://candis.insilab.org. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Transcriptomics of receptive endometrium in women with sonographic features of adenomyosis.

BACKGROUND: Women with uterine adenomyosis seeking assisted reproduction have been associated with compromised endometrial receptivity to embryo implantation. To understand the mechanisms involved in this process, we aimed to compare endometrial transcriptome profiles during the window of implantation (WOI) between women with and without adenomyosis. METHODS: We obtained endometrial biopsies LH-timed to the WOI from women with sonographic features of adenomyosis (n=10) and controls (n=10). Isolated RNA samples were subjected to RNA sequencing (RNA-seq) by the Illumina NovaSeq 6000 platform and endometrial receptivity classification with a molecular tool for menstrual cycle phase dating (beREADY®, CCHT). The program language R and Bioconductor packages were applied to analyse RNA-seq data in the setting of the result of accurate endometrial dating. To suggest robust candidate pathways, the identified differentially expressed genes (DEGs) associated with the adenomyosis group in the receptive phase were further integrated with 151, 173 and 42 extracted genes from published studies that were related to endometrial receptivity in healthy uterus, endometriosis and adenomyosis, respectively. Enrichment analyses were performed using Cytoscape ClueGO and CluePedia apps. RESULTS: Out of 20 endometrial samples, 2 were dated to the early receptive phase, 13 to the receptive phase and 5 to the late receptive phase. Comparison of the transcriptomics data from all 20 samples provided 909 DEGs (p<0.05; nonsignificant after adjusted p value) in the adenomyosis group but only 4 enriched pathways (Bonferroni p value < 0.05). The analysis of 13 samples only dated to the receptive phase provided suggestive 382 DEGs (p<0.05; nonsignificant after adjusted p value) in the adenomyosis group, leading to 33 enriched pathways (Bonferroni p value < 0.05). These included pathways were already associated with endometrial biology, such as "Expression of interferon (IFN)-induced genes" and "Response to IFN-alpha". Data integration revealed pathways indicating a unique effect of adenomyosis on endometrial molecular organization (e.g., "Expression of IFN-induced genes") and its interference with endometrial receptivity establishment (e.g., "Extracellular matrix organization" and "Tumour necrosis factor production"). CONCLUSIONS: Accurate endometrial dating and RNA-seq analysis resulted in the identification of altered response to IFN signalling as the most promising candidate of impaired uterine receptivity in adenomyosis.

Genetic variants of the hypoxia-inducible factor 3 alpha subunit (Hif3a) gene in the Fat and Lean mouse selection lines.

BACKGROUND: Adipose tissue hypoxia and members of the hypoxia-inducible factor alpha (HIFA) are involved in development of obesity. However, the mechanism and functions of HIF3A, one of three HIFA paralogs, in fat deposition have not been sufficiently studied. METHODS AND RESULTS: In the present study, we investigated whether Hif3a sequence variants are associated with divergent fat deposition in mouse selection lines for fatness and leanness. Sequencing and RFLP were used to analyse sequence variants within Hif3a. To identify candidate regulatory variants, we performed literature screening and used databases and bioinformatics tools like Ensembl, MethPrimer, TargetScanMouse, miRDB, PolyAsite, RISE, LncRRIsearch, RNAfold, PredictProtein, CAIcal, and switches.ELM Resource. There are 90 sequence variants in Hif3a between the two mouse lines. While most Fat line variants locate within intronic regions, Lean line variants are mainly in 3' UTR. We constructed a map of Hif3a potential regulatory regions and identified 39 regulatory variants by integrating data on constrained and regulatory elements, CpGs, and miRNAs and lncRNAs binding sites. Moreover, 3' UTR and two exonic variants may influence mRNA stability, translation rate and protein functionality. We propose as priority candidates for further functional studies a missense (rs37398126) and synonymous (rs37739792) variants, and intronic (rs47471302) variant that overlap conserved element in promoter region and predicted lncRNAs binding site. CONCLUSION: The results indicate a potential involvement of Hif3a in fat deposition. Additionally, approach used in the present study may serve as a general guideline for constructing an integrative gene map for prioritizing candidate gene variants with phenotypic effects.

MicroRNAs in Leukemias: A Clinically Annotated Compendium.

Leukemias are a group of malignancies of the blood and bone marrow. Multiple types of leukemia are known, however reliable treatments have not been developed for most leukemia types. Furthermore, even relatively reliable treatments can result in relapses. MicroRNAs (miRNAs) are a class of short, noncoding RNAs responsible for epigenetic regulation of gene expression and have been proposed as a source of potential novel therapeutic targets for leukemias. In order to identify central miRNAs for leukemia, we conducted data synthesis using two databases: miRTarBase and DISNOR. A total of 137 unique miRNAs associated with 16 types of leukemia were retrieved from miRTarBase and 86 protein-coding genes associated with leukemia were retrieved from the DISNOR database. Based on these data, we formed a visual network of 248 miRNA-target interactions (MTI) between leukemia-associated genes and miRNAs associated with ≥4 leukemia types. We then manually reviewed the literature describing these 248 MTIs for interactions identified in leukemia studies. This manually curated data was then used to visualize a network of 64 MTIs identified in leukemia patients, cell lines and animal models. We also formed a visual network of miRNA-leukemia associations. Finally, we compiled leukemia clinical trials from the ClinicalTrials database. miRNAs with the highest number of MTIs were miR-125b-5p, miR-155-5p, miR-181a-5p and miR-19a-3p, while target genes with the highest number of MTIs were TP53, BCL2, KIT, ATM, RUNX1 and ABL1. The analysis of 248 MTIs revealed a large, highly interconnected network. Additionally, a large MTI subnetwork was present in the network visualized from manually reviewed data. The interconnectedness of the MTI subnetwork suggests that certain miRNAs represent central disease molecules for multiple leukemia types. Additional studies on miRNAs, their target genes and associated biological pathways are required to elucidate the therapeutic potential of miRNAs in leukemia.

Genetic analysis of 39 erythrocytosis and hereditary hemochromatosis-associated genes in the Slovenian family with idiopathic erythrocytosis.

BACKGROUND: Erythrocytosis is a condition with an excessive number of erythrocytes, accompanied by an elevated haemoglobin and/or haematocrit value. Congenital erythrocytosis has a diverse genetic background with several genes involved in erythropoiesis. In clinical practice, nine genes are usually examined, but in approximately 70% of patients, no causative mutation can be identified. In this study, we screened 39 genes, aiming to identify potential disease-driving variants in the family with erythrocytosis of unknown cause. PATIENTS AND METHODS: Two affected family members with elevated haemoglobin and/or haematocrit and negative for acquired causes and one healthy relative from the same family were selected for molecular-genetic analysis of 24 erythrocytosis and 15 hereditary haemochromatosis-associated genes with targeted NGS. The identified variants were further analysed for pathogenicity using various bioinformatic tools and review of the literature. RESULTS: Of the 12 identified variants, two heterozygous variants, the missense variant c.471G>C (NM_022051.2) (p.(Gln157His)) in the EGLN1 gene and the intron variant c.2572-13A>G (NM_004972.3) in the JAK2 gene, were classified as low-frequency variants in European population. None of the two variants were present in a healthy family member. Variant c.2572-13A>G has potential impact on splicing by one prediction tool. CONCLUSION: For the first time, we included 39 genes in the erythrocytosis clinical panel and identified two potential disease-driving variants in the Slovene family studied. Based on the reported functional in vitro studies combined with our bioinformatics analysis, we suggest further functional analysis of variant in the JAK2 gene and evaluation of a cumulative effect of both variants.

Identification of Novel RNA Binding Proteins Influencing Circular RNA Expression in Hepatocellular Carcinoma.

Circular RNAs (circRNAs) are increasingly recognized as having a role in cancer development. Their expression is modified in numerous cancers, including hepatocellular carcinoma (HCC); however, little is known about the mechanisms of their regulation. The aim of this study was to identify regulators of circRNAome expression in HCC. Using publicly available datasets, we identified RNA binding proteins (RBPs) with enriched motifs around the splice sites of differentially expressed circRNAs in HCC. We confirmed the binding of some of the candidate RBPs using ChIP-seq and eCLIP datasets in the ENCODE database. Several of the identified RBPs were found to be differentially expressed in HCC and/or correlated with the overall survival of HCC patients. According to our bioinformatics analyses and published evidence, we propose that NONO, PCPB2, PCPB1, ESRP2, and HNRNPK are candidate regulators of circRNA expression in HCC. We confirmed that the knocking down the epithelial splicing regulatory protein 2 (ESRP2), known to be involved in the maintenance of the adult liver phenotype, significantly changed the expression of candidate circRNAs in a model HCC cell line. By understanding the systemic changes in transcriptome splicing, we can identify new proteins involved in the molecular pathways leading to HCC development and progression.

Molecular signature of eutopic endometrium in endometriosis based on the multi-omics integrative synthesis.

PURPOSE: To synthesise data from genome-wide studies reporting molecular signature of eutopic endometrium through the phases of the menstrual cycle in endometriosis. METHODS: Extraction of data from publications reporting genetic signatures characterising endometrium associated with endometriosis. The nomenclature of extracted differentially expressed transcripts and proteins was adopted according to the HUGO Gene Nomenclature Committee (HGNC). Loci were further sorted according to the different phases of the menstrual cycle, i.e. menstrual (M), proliferative (P), secretory (S), early-secretory (ES), mid-secretory (MS), late-secretory (LS), and not specified (N/S) if the endometrial dating was not available. Enrichment analysis was performed using the DAVID bioinformatics tool. RESULTS: Altered molecular changes were reported by 21 studies, including 13 performed at the transcriptomic, 6 at proteomic, and 2 at epigenomic level. Extracted data resulted in a catalogue of total 670 genetic causes with available 591 official gene symbols, i.e. M = 3, P = 188, S = 81, ES = 82, MS = 173, LS = 36, and N/S = 28. Enriched pathways included oestrogen signalling pathway, extracellular matrix organization, and endothelial cell chemotaxis. Our study revealed that knowledge of endometrium biology in endometriosis is fragmented due to heterogeneity of published data. However, 15 genes reported as dysregulated by at least two studies within the same phase and 33 significantly enriched GO-BP terms/KEGG pathways associated with different phases of the menstrual cycle were identified. CONCLUSIONS: A multi-omics insight into molecular patterns underlying endometriosis could contribute towards identification of endometrial pathological mechanisms that impact fertility capacities of women with endometriosis.

Genetic variability of serotonin pathway associated with schizophrenia onset, progression, and treatment.

Schizophrenia (SZ) onset and treatment outcome have important genetic components, however individual genes do not have strong effects on SZ phenotype. Therefore, it is important to use the pathway-based approach and study metabolic and signaling pathways, such as dopaminergic and serotonergic. Serotonin pathway has an important role in brain signaling, nevertheless, its role in SZ is not as thoroughly examined as that of dopamine pathway. In this study, we reviewed serotonin pathway genes and genetic variations associated with SZ, including variations at DNA, RNA, and epigenetic level. We obtained 30 serotonin pathway genes from Kyoto encyclopedia of genes and genomes and used these genes for the literature review. We extracted 20 protein coding serotonin pathway genes with genetic variations associated with SZ onset, development, and treatment from 31 research papers. Genes associated with SZ are present on all levels of serotonin pathway: serotonin synthesis, transport, receptor binding, intracellular signaling, and reuptake; however, regulatory genes are poorly researched. We summarized common challenges of genetic association studies and presented some solutions. The analysis of reported serotonin pathway-SZ associations revealed lack of information about certain serotonin pathway genes potentially associated with SZ. Furthermore, it is becoming clear that interactions among serotonin pathway genes and their regulators may bring further knowledge about their involvement in SZ.

Genetic variability of hypoxia-inducible factor alpha (HIFA) genes in familial erythrocytosis: Analysis of the literature and genome databases.

Familial erythrocytosis (FE) is a congenital disorder, defined by elevated red blood cell number, hemoglobin, and hematocrit. Among eight types of FE, type 4 is caused by variants in the EPAS1 gene. Two other hypoxia-inducible factor alpha (HIFA) subunits, HIF1A and HIF3A, have not yet been associated with medical history of FE, but have potential role in the development of erythrocytosis. To improve diagnosis, it is crucial to identify new variants in genes involved in erythrocyte production. Published literature and data from genome browsers were used to obtain HIFA sequence variants associated with erythrocytosis and to locate them on protein sequence and regulatory sites. We retrieved 24 variants from the literature: 2 in HIF1A, 20 in EPAS1 and 2 in HIF3A gene. Sixteen out of 20 variants in the EPAS1 gene are positioned in a conserved region of 13 amino acids within exon 12, next to regulatory post-translational modification and binding sites, suggesting that EPAS1 has an important role in erythropoiesis. The role of HIF1A and HIF3A in the development of erythrocytosis should be further investigated.

GenProBiS: web server for mapping of sequence variants to protein binding sites.

Discovery of potentially deleterious sequence variants is important and has wide implications for research and generation of new hypotheses in human and veterinary medicine, and drug discovery. The GenProBiS web server maps sequence variants to protein structures from the Protein Data Bank (PDB), and further to protein-protein, protein-nucleic acid, protein-compound, and protein-metal ion binding sites. The concept of a protein-compound binding site is understood in the broadest sense, which includes glycosylation and other post-translational modification sites. Binding sites were defined by local structural comparisons of whole protein structures using the Protein Binding Sites (ProBiS) algorithm and transposition of ligands from the similar binding sites found to the query protein using the ProBiS-ligands approach with new improvements introduced in GenProBiS. Binding site surfaces were generated as three-dimensional grids encompassing the space occupied by predicted ligands. The server allows intuitive visual exploration of comprehensively mapped variants, such as human somatic mis-sense mutations related to cancer and non-synonymous single nucleotide polymorphisms from 21 species, within the predicted binding sites regions for about 80 000 PDB protein structures using fast WebGL graphics. The GenProBiS web server is open and free to all users at http://genprobis.insilab.org.

Genome-wide screening for smallest regions of overlaps in cryptorchidism.

Cryptorchidism is a urogenital abnormality associated with increased rates of testicular neoplasia and impaired spermatogenesis. The field is facing expansion of genomics data; however, it lacks protocols for biomarker prioritization. Identification of smallest regions of overlap (SRO) presents an approach for candidate gene identification but has not yet been systematically conducted in cryptorchidism. The aim of this study was to conduct a genome-wide screening for SRO (GW-SRO) associated with cryptorchidism development. We updated the Cryptorchidism Gene Database to version 3, remapped genomic coordinates of loci from older assemblies to the GRCh38 and performed genome-wide screening for overlapping regions associated with cryptorchidism risk. A total of 73 chromosomal loci (68 involved in chromosomal mutations and five copy number variations) described in 37 studies associated with cryptorchidism risk in humans were used for SRO identification. Analysis resulted in 18 SRO, based on deletions, duplications, inversions, derivations and copy number variations. Screening for SRO was challenging owing to heterogeneous reporting of genomic locations. To our knowledge, this is the first GW-SRO study for cryptorchidism and it presents the basis for further narrowing of critical regions for cryptorchidism and planning functional experiments. The developed protocol could also be applied to other multifactorial diseases.

In silico screening of the chicken genome for overlaps between genomic regions: microRNA genes, coding and non-coding transcriptional units, QTL, and genetic variations.

MicroRNAs (miRNAs) are a class of non-coding RNAs involved in posttranscriptional regulation of target genes. Regulation requires complementarity between target mRNA and the mature miRNA seed region, responsible for their recognition and binding. It has been estimated that each miRNA targets approximately 200 genes, and genetic variability of miRNA genes has been reported to affect phenotypic variability and disease susceptibility in humans, livestock species, and model organisms. Polymorphisms in miRNA genes could therefore represent biomarkers for phenotypic traits in livestock animals. In our previous study, we collected polymorphisms within miRNA genes in chicken. In the present study, we identified miRNA-related genomic overlaps to prioritize genomic regions of interest for further functional studies and biomarker discovery. Overlapping genomic regions in chicken were analyzed using the following bioinformatics tools and databases: miRNA SNiPer, Ensembl, miRBase, NCBI Blast, and QTLdb. Out of 740 known pre-miRNA genes, 263 (35.5 %) contain polymorphisms; among them, 35 contain more than three polymorphisms The most polymorphic miRNA genes in chicken are gga-miR-6662, containing 23 single nucleotide polymorphisms (SNPs) within the pre-miRNA region, including five consecutive SNPs, and gga-miR-6688, containing ten polymorphisms including three consecutive polymorphisms. Several miRNA-related genomic hotspots have been revealed in chicken genome; polymorphic miRNA genes are located within protein-coding and/or non-coding transcription units and quantitative trait loci (QTL) associated with production traits. The present study includes the first description of an exonic miRNA in a chicken genome, an overlap between the miRNA gene and the exon of the protein-coding gene (gga-miR-6578/HADHB), and the first report of a missense polymorphism located within a mature miRNA seed region. Identified miRNA-related genomic hotspots in chicken can serve researchers as a starting point for further functional studies and association studies with poultry production and health traits and the basis for systematic screening of exonic miRNAs and missense/miRNA seed polymorphisms in other genomes.

Transcription factor HIF1A: downstream targets, associated pathways, polymorphic hypoxia response element (HRE) sites, and initiative for standardization of reporting in scientific literature.

Hypoxia-inducible factor-1α (HIF-1α) has crucial role in adapting cells to hypoxia through expression regulation of many genes. Identification of HIF-1α target genes (HIF-1α-TGs) is important for understanding the adapting mechanism. The aim of the present study was to collect known HIF-1α-TGs and identify their associated pathways. Targets and associated genomics data were retrieved using PubMed, WoS ( http://apps.webofknowledge.com/ ), HGNC ( http://www.genenames.org/ ), NCBI ( http://www.ncbi.nlm.nih.gov/ ), Ensemblv.84 ( http://www.ensembl.org/index.html ), DAVID Bioinformatics Resources ( https://david.ncifcrf.gov /), and Disease Ontology database ( http://disease-ontology.org/ ). From 51 papers, we collected 98 HIF-1α TGs found to be associated with 20 pathways, including metabolism of carbohydrates and pathways in cancer. Reanalysis of genomic coordinates of published HREs (hypoxia response elements) revealed six polymorphisms within HRE sites (HRE-SNPs): ABCG2, ACE, CA9, and CP. Due to large heterogeneity of results presentation in scientific literature, we also propose a first step towards reporting standardization of HIF-1α-target interactions consisting of ten relevant data types. Suggested minimal checklist for reporting will enable faster development of a complete catalog of HIF-1α-TGs, data sharing, bioinformatics analyses, and setting novel more targeted hypotheses. The proposed format for data standardization is not yet complete but presents a baseline for further optimization of the protocol with additional details, for example, regarding the experimental validation.

Construction of an integrative regulatory element and variation map of the murine Tst locus.

BACKGROUND: Given the abundance of new genomic projects and gene annotations, researchers trying to pinpoint causal genetic variants are faced with a challenging task of how to efficiently integrate all current genomic information. The objective of the study was to develop an approach to integrate various genomic annotations for a recently positionally-cloned Tst gene (Thiosulfate Sulfur Transferase, synonym Rhodanese) responsible for the Fob3b2 QTL effect on leanness and improved metabolic parameters. The second aim was to identify and prioritize Tst genetic variants that may be causal for the phenotypic effects. RESULTS: A bioinformatics approach was developed to integrate existing knowledge of regulatory elements of the Tst gene. The entire Tst locus along with flanking segments was sequenced between our unique polygenic mouse Fat and Lean strains that were generated by divergent selection on adiposity for over 60 generations. The bioinformatics-generated regulatory element map of the Tst locus was then combined with genetic variants between the Fat and Lean mice and with comparative analyses of polymorphisms across 17 mouse strains in order to prioritise likely causal polymorphisms. Two candidate regulatory variants were identified, one overlapping an evolutionary constrained Tst intronic element and the other residing in the seed region of a predicted 3'UTR miRNA binding site. CONCLUSIONS: This study developed a map of regulatory elements for the Tst locus in mice and identified candidate genetic variants with increased causal likelihood. This map provides a basis for experimental validation and functional analyses of this novel candidate leanness and antidiabetic gene. Our methodological approach is of general utility for analyzing regulation of loci that have limited annotations and experimental evidence and for identifying candidate causal regulatory genetic variants in post-GWAS or post-QTL- cloning studies.

MicroRNA epigenetic signatures in human disease.

MicroRNAs (miRNAs) are short non-coding RNAs that act as important regulators of gene expression as part of the epigenetic machinery. In addition to posttranscriptional gene silencing by miRNAs, the epigenetic mechanisms also include DNA methylation, histone modifications and their crosstalk. Epigenetic modifications were reported to play an important role in many disease onsets and progressions and can be used to explain several features of complex diseases, such as late onset and fluctuation of symptoms. However, miRNAs not only function as a part of epigenetic machinery, but are also epigenetically modified by DNA methylation and histone modification like any other protein-coding gene. There is a strong connection between epigenome and miRNome, and any dysregulation of this complex system can result in various physiological and pathological conditions. In addition, miRNAs play an important role in toxicogenomics and may explain the relationship between toxicant exposure and tumorigenesis. The present review provides information on 63 miRNA genes shown to be epigenetically regulated in association with 21 diseases, including 11 cancer types: cardiac fibrosis, cardiovascular disease, preeclampsia, Hirschsprung's disease, rheumatoid arthritis, systemic sclerosis, systemic lupus erythematosus, temporal lobe epilepsy, autism, pulmonary fibrosis, melanoma, acute myeloid leukemia, chronic lymphocytic leukemia, colorectal, gastric, cervical, ovarian, prostate, lung, breast, and bladder cancer. The review revealed that hsa-miR-34a, hsa-miR-34b, and hsa-miR-34c are the most frequently reported epigenetically dysregulated miRNAs. There is a need to further study molecular mechanisms of various diseases to better understand the crosstalk between epigenetics and gene expression and to develop new therapeutic options and biomarkers.

Catalog of genetic variants within mature microRNA seed regions in chicken.

MicroRNA (miRNA) is a class of noncoding RNA important in posttranscriptional regulation of target genes. The regulation mechanism requires complementarity between target mRNA and the miRNA region responsible for their recognition and binding, also called the seed region. It has been estimated that each miRNA targets approximately 200 genes and genetic variability of miRNA genes has been associated with phenotypic variation and disease susceptibility in humans, livestock species, and model organisms. Polymorphisms in miRNA genes especially within the seed region could therefore represent biomarkers for phenotypic traits important in livestock animals. Using the updated Version 5.0 of our previously developed bioinformatics tool miRNA SNiPer we assembled polymorphic miRNA genes in chicken. Out of 740 miRNA genes 263 were polymorphic, among them 77 had SNPs located within the mature region, and 29 of them within the miRNA seed region. Because several polymorphisms in databases result from sequencing errors, we performed experimental validation of polymorphisms located within 4 selected miRNA genes in chicken (gga-mir-1614, -1644, -1648, and -1657). We confirmed the presence of nine polymorphisms and identified 3 additional novel polymorphisms within primary miRNA regions in chicken representing 3 layer-type breeds, one layer-type hybrid, and one meat-type intercrossed population. The developed catalog of mir-SNPs in chicken can serve researchers as a starting point for association studies dealing with poultry production traits and designing functional experiments.

The decalog of long non-coding RNA involvement in cancer diagnosis and monitoring.

Long non-coding RNAs (lncRNAs) are transcripts without protein-coding capacity; initially regarded as "transcriptional noise", lately they have emerged as essential factors in both cell biology and mechanisms of disease. In this article, we present basic knowledge of lncRNA molecular mechanisms, associated physiological processes and cancer association, as well as their diagnostic and therapeutic value in the form of a decalog: (1) Non-coding RNAs (ncRNAs) are transcripts without protein-coding capacity divided by size (short and long ncRNAs), function (housekeeping RNA and regulatory RNA) and direction of transcription (sense/antisense, bidirectional, intronic and intergenic), containing a broad range of molecules with diverse properties and functions, such as messenger RNA, transfer RNA, microRNA and long non-coding RNAs. (2) Long non-coding RNAs are implicated in many molecular mechanisms, such as transcriptional regulation, post-transcriptional regulation and processing of other short ncRNAs. (3) Long non-coding RNAs play an important role in many physiological processes such as X-chromosome inactivation, cell differentiation, immune response and apoptosis. (4) Long non-coding RNAs have been linked to hallmarks of cancer: (a) sustaining proliferative signaling; (b) evading growth suppressors; (c) enabling replicative immortality; (d) activating invasion and metastasis; (e) inducing angiogenesis; (f) resisting cell death; and (g) reprogramming energy metabolism. (5) Regarding their impact on cancer cells, lncRNAs are divided into two groups: oncogenic and tumor-suppressor lncRNAs. (6) Studies of lncRNA involvement in cancer usually analyze deregulated expression patterns at the RNA level as well as the effects of single nucleotide polymorphisms and copy number variations at the DNA level. (7) Long non-coding RNAs have potential as novel biomarkers due to tissue-specific expression patterns, efficient detection in body fluids and high stability. (8) LncRNAs serve as novel biomarkers for diagnostic, prognostic and monitoring purposes. (9) Tissue specificity of lncRNAs enables the development of selective therapeutic options. (10) Long non-coding RNAs are emerging as commercial biomarkers and therapeutic agents.

Multiple omics levels of chronic lymphocytic leukemia.

Chronic lymphocytic leukemia (CLL) is a lymphoproliferative malignancy characterized by the proliferation of functionally mature but incompetent B cells. It is the most prevalent type of leukemia in Western populations, accounting for approximately 25% of new leukemia cases. While recent advances, such as ibrutinib and venetoclax treatment have improved patient outlook, aggressive forms of CLL such as Richter transformation still pose a significant challenge. This discrepancy may be due to the heterogeneity of factors contributing to CLL development at multiple -omics levels. However, information on the omics of CLL is fragmented, hindering multi-omics-based research into potential treatment options. To address this, we aggregated and presented a selection of important aspects of various omics levels of the disease in this review. The purpose of the present literature analysis is to portray examples of CLL studies from different omics levels, including genomics, epigenomics, transcriptomics, epitranscriptomics, proteomics, epiproteomics, metabolomics, glycomics and lipidomics, as well as those identified by multi-omics approaches. The review includes the list of 102 CLL-associated genes with relevant genomics information. While single-omics studies yield substantial and useful data, they omit a significant level of complex biological interplay present in the disease. As multi-omics studies integrate several different layers of data, they may be better suited for complex diseases such as CLL and have thus far yielded promising results. Future multi-omics studies may assist clinicians in improved treatment choices based on CLL subtypes as well as allow the identification of novel biomarkers and targets for treatments.