AmiCa: Atlas of miRNA-gene correlations in cancer.

The increasing availability of RNA sequencing data has opened up numerous opportunities to analyze various RNA interactions, including microRNA-target interactions (MTIs). In response to the necessity for a specialized tool to study MTIs in cancer and normal tissues, we developed AmiCa (https://amica.omics.si/), a web server designed for comprehensive analysis of mature microRNA (miRNA) and gene expression in 32 cancer types. Data from 9498 tumor samples and 626 normal samples from The Cancer Genome Atlas were obtained through the Genomic Data Commons and used to calculate differential expression and miRNA-target gene (MTI) correlations. AmiCa provides data on differential expression of miRNAs/genes for cancers for which normal tissue samples were available. In addition, the server calculates and presents correlations separately for tumor and normal samples for cancers for which normal samples are available. Furthermore, it enables the exploration of miRNA/gene expression in all cancer types with different miRNA/gene expression. In addition, AmiCa includes a ranking system for genes and miRNAs that can be used to identify those that are particularly highly expressed in certain cancers compared to other cancers, facilitating targeted and cancer-specific research. Finally, the functionality of AmiCa is illustrated by two case studies.

Preservation of 5-Hydroxymethylcytosine Levels in LRIG1 across Genomic DNA and Cell-Free DNA in Glioma Patients.

Cell-free DNA (cfDNA) has recently emerged as a promising minimally invasive diagnostic biomarker for various cancers. In this study, our aim was to identify cfDNA biomarkers by investigating genes that displayed significant differences between glioma patients and their corresponding controls. To accomplish this, we utilized publicly available data from the Gene Expression Omnibus, focusing on 5-hydroxymethylcytosine (5hmC) profiles in both cfDNA and genomic DNA (gDNA) from glioma patients and healthy individuals. The intersection of gene lists derived from these comparative analyses unveiled LRIG1 and ZNF703 as the two genes with elevated 5hmC levels in both the cfDNA of glioma patients and gDNA of glioma tissue compared to their respective controls. The gene expression data revealed both genes were upregulated in glioma tissue compared to normal brain tissue. Integration of 5hmC data revealed a strong positive correlation in the glioma tissue group between 5hmC and the gene expression of the LRIG1 gene. Furthermore, exploration using the AmiCa web tool indicated that LRIG1 gene expression was elevated compared to 17 other cancers included in the database, emphasizing its potential as a distinctive biomarker across multiple cancer types.

Steadfast Toll Like Receptor 4 (TLR4) 5-Hydroxymethylcytosine Levels in Cell-Free DNA: A Promising Consistency Marker for Colorectal Cancer Patients.

Cell-free DNA (cfDNA) from patient blood is emerging as a noninvasive diagnostic avenue for various cancers. We aimed to identify reliable biomarkers in cfDNA by investigating genes exhibiting significant differences between colorectal cancer and control samples. Our objective was to identify genes that showed a positive difference between cancer and control samples. To achieve this, we conducted an in silico analysis to identify genes that exhibit no significant variation in methylation between genomic DNA (gDNA) and cfDNA. We collected experimental data from publicly available repositories, which included 5-hydroxymethylcytosine (5hmC) profiles of gDNA and cfDNA samples from both cancer patients and healthy individuals. By comparing and overlapping these two groups, we identified 187 genes of interest, of which 53 genes had a positive difference among colon cancer patients and healthy individuals. Next, we performed an ANOVA test on these genes, resulting in the identification of 12 genes that showed statistically significant higher levels of 5hmC in cfDNA and gDNA from cancer patients compared to healthy individuals. Additionally, we compared the 5hmC status of these genes between cfDNA and gDNA from cancer patients. Interestingly, we found that the 5hmC of the toll like receptor 4 (TLR4) gene was not statistically different between cfDNA and gDNA from cancer patients, indicating consistency between cfDNA and gDNA. These findings have important implications, not only for experimental validation but also for the development of more sensitive and robust noninvasive methods to improve diagnostic, prognostic, and treatment options for colon cancer.

In Silico Gene Prioritization Highlights the Significance of Bone Morphogenetic Protein 4 (BMP4) Promoter Methylation across All Methylation Clusters in Colorectal Cancer.

The cytosine-phosphate-guanine (CpG) island methylator phenotype (CIMP) represents one of the pathways involved in the development of colorectal cancer, characterized by genome-wide hypermethylation. To identify samples exhibiting hypermethylation, we used unsupervised hierarchical clustering on genome-wide methylation data. This clustering analysis revealed the presence of four distinct subtypes within the tumor samples, namely, CIMP-H, CIMP-L, cluster 3, and cluster 4. These subtypes demonstrated varying levels of methylation, categorized as high, intermediate, and very low. To gain further insights, we mapped significant probes from all clusters to Ensembl Regulatory build 89, with a specific focus on those located within promoter regions or bound regions. By intersecting the methylated promoter and bound regions across all methylation subtypes, we identified a total of 253 genes exhibiting aberrant methylation patterns in the promoter regions across all four subtypes of colorectal cancer. Among these genes, our comprehensive genome-wide analysis highlights bone morphogenic protein 4 (BMP4) as the most prominent candidate. This significant finding was derived through the utilization of various bioinformatics tools, emphasizing the potential role of BMP4 in colorectal cancer development and progression.

Genes associated with somatic cell count index in Brown Swiss cattle.

Subclinical mastitis (SM) is one of the most common diseases of cows in milk production herds caused by contagious and/or environmental pathogens. Since there are no visible abnormalities in the milk or udder, the detection of SM requires special diagnostic tests. Somatic cell count (SCC) is the most common test used to detect changes in milk due to the inflammatory process. Previously, we developed somatic cell count index (SCCI), a new method for the accurate prediction of milk yield losses caused by elevated SCC. The aim of this study was to identify new candidate genetic markers for SCCI in the Slovenian population of Brown Swiss (BS) cattle. For that purpose, we analyzed samples of BS cows, which were genotyped using single-nucleotide polymorphism (SNP) microarray ICBF International Dairy and Beef v3 (ICBF, Ireland) for a total of 53,262 SNP markers. After quality control, the set of 18,136 SNPs was used in association analysis. Our association analysis revealed that 130 SNPs were associated with SCCI, which were used for haplotype and overlap analysis. Haplotypes generated from the genotyped data for those 130 SNPs revealed 10 haplotype blocks among 22 SNPs. Additionally, all 130 SNPs, mastitis-related quantitative trait loci, and protein-coding genes are shown on the bovine genome. Overlap analysis shows that the majority of significantly associated SNPs (70) are intergenic, while 60 SNPs are mapped within, upstream, or downstream of the protein-coding genes. However, those genes can serve as strong candidate genes for the marker-assisted selection programs in our and possibly other populations of cattle.

The integrative knowledge base for miRNA-mRNA expression in colorectal cancer.

"miRNA colorectal cancer" (https://mirna-coadread.omics.si/) is a freely available web application for studying microRNA and mRNA expression and their correlation in colorectal cancer. To the best of our knowledge, "miRNA colorectal cancer" has the largest knowledge base of miRNA-target gene expressions and correlations in colorectal cancer, based on the largest available sample size from the same source of data. Data from high-throughput molecular profiling of 295 colon and rectum adenocarcinoma samples from The Cancer Genome Atlas was analyzed and integrated into our knowledge base. The objective of developing this web application was to help researchers to discover the behavior and role of miRNA-target gene interactions in colorectal cancer. For this purpose, results of differential expression and correlation analyses of miRNA and mRNA data collected in our knowledge base are available through web forms. To validate our knowledge base experimentally, we selected genes FN1, TGFB2, RND3, ZEB1 and ZEB2 and miRNAs hsa-miR-200a/b/c-3p, hsa-miR-141-3p and hsa-miR-429. Both approaches revealed a negative correlation between miRNA hsa-miR-200b/c-3p and its target gene FN1 and between hsa-miR-200a-3p and its target TGFB2, thus supporting the usefulness of the developed knowledge base.

Association between histopathological changes and expression of selected microRNAs in skin of adult patients with IgA vasculitis.

AIMS: IgA vasculitis (IgAV) is a common small-vessel systemic vasculitisthat is histologically characterised by granulocyte infiltration and IgA deposition in vessel walls. Information on microRNA (miRNA) involvement inIgAVis limited. The aim of this study was to analyse the association between histopathological changes and expression profiles of 14 miRNAs in the affected skin of 70 adult patients with IgAV. METHODS AND RESULTS: miRNA expression analysis was performed by quantitative real-time polymerase chain reaction and evaluation of histopathological changes by light and immunofluorescence microscopy on formalin-fixed paraffin-embedded skin excision samples. In IgAV-affected skin, granulocyte infiltration was significantly associated with vessel fibrinoid necrosis. Of the analysed miRNAs, four showed two-fold increased expression (let-7d, let-7f, miR-21-5p, and miR-203-3p), five showed five-fold increased expression (let-7b, miR-17-5p, miR-155-5p, miR-423-5p, and miR-451a), and threeshowed 15-fold increased expression (let-7a, miR-21-3p, miR-223-3p), as compared with controls (all P < 0.001). miR-146a-5p and miR-148b-3p showed three-fold decreased expression (P = 0.981 and P < 0.001). The expression of miR-223-3p also showed a significant positive association with granulocyte infiltration and fibrinoid necrosis. CONCLUSIONS: Altered miRNA expression, especially of miRNA-223-3p, may be associated with the skin inflammatory state in IgAV. The majority of aberrantly expressed miRNAs in IgAV-affected skin are known to influence the nuclear factor-κB signalling pathway, which is crucial for activation of key proinflammatory genes, including those encoding tumour necrosis factor-α, interleukin (IL)-6, and IL-8. Furthermore, miR-146a-5p and miR-148b-3p, which are negative regulators of inflammatory gene expression, showed decreased expression and could contribute to the exaggerated inflammation. Further investigation of miRNA expression in the affected tissues could improve our knowledge of IgAV pathogenesis, and possibly help to identify novel biomarkers in body fluids.

Genes CEP55, FOXD3, FOXF2, GNAO1, GRIA4, and KCNA5 as potential diagnostic biomarkers in colorectal cancer.

BACKGROUND: Colorectal cancer (CRC) is one of the leading causes of death by cancer worldwide and in need of novel potential diagnostic biomarkers for early discovery. METHODS: We conducted a two-step study. We first employed bioinformatics on data from The Cancer Genome Atlas to obtain potential biomarkers and then experimentally validated some of them on our clinical samples. Our aim was to find a methylation alteration common to all clusters, with the potential of becoming a diagnostic biomarker in CRC. RESULTS: Unsupervised clustering of methylation data resulted in four clusters, none of which had a known common genetic or epigenetic event, such as mutations or methylation. The intersect among clusters and regulatory regions resulted in 590 aberrantly methylated probes, belonging to 198 differentially expressed genes. After performing pathway and functional analysis on differentially expressed genes, we selected six genes: CEP55, FOXD3, FOXF2, GNAO1, GRIA4 and KCNA5, for further experimental validation on our own clinical samples. In silico analysis demonstrated that CEP55 was hypomethylated in 98.7% and up-regulated in 95.0% of samples. Genes FOXD3, FOXF2, GNAO1, GRIA4 and KCNA5 were hypermethylated in 97.9, 81.1, 80.3, 98.4 and 94.0%, and down-regulated in 98.3, 98.9, 98.1, 98.1 and 98.6% of samples, respectively. Our experimental data show CEP55 was hypomethylated in 97.3% of samples and down-regulated in all samples, while FOXD3, FOXF2, GNAO1, GRIA4 and KCNA5 were hypermethylated in 100.0, 90.2, 100.0, 99.1 and 100.0%, and down-regulated in 68.0, 76.0, 96.0, 95.2 and 84.0% of samples, respectively. Results of in silico and our experimental analyses showed that more than 97% of samples had at least four methylation markers altered. CONCLUSIONS: Using bioinformatics followed by experimental validation, we identified a set of six genes that were differentially expressed in CRC compared to normal mucosa and whose expression seems to be methylation dependent. Moreover, all of these six genes were common in all methylation clusters and mutation statuses of CRC and as such are believed to be an early event in human CRC carcinogenesis and to represent potential CRC biomarkers.

Comparison of SNP Genotypes Related to Proliferative Vitreoretinopathy (PVR) across Slovenian and European Subpopulations.

The present study investigated the distribution of genotypes within single nucleotide polymorphisms (SNPs) in genes, related to PVR pathogenesis across European subpopulations. Genotype distributions of 42 SNPs among 96 Slovenian healthy controls were investigated and compared to genotype frequencies in 503 European individuals (Ensembl database) and their subpopulations. Furthermore, a case-control status was simulated to evaluate effects of allele frequency changes on statistically significant results in gene-association studies investigating functional polymorphisms. In addition, 96 healthy controls were investigated within 4 SNPs: rs17561 (IL1A), rs2069763 (IL2), rs2229094 (LTA), and rs1800629 (TNF) in comparison to PVR patients. Significant differences (P < 0.05) in distribution of genotypes among 96 Slovenian participants and a European population were found in 10 SNPs: rs3024498 (IL10), rs315952 (IL1RN), rs2256965 (LST1), rs2256974 (LST1), rs909253 (LTA), rs2857602 (LTA), rs3138045 (NFKB1A), rs3138056 (NFKB1A), rs7656613 (PDGFRA), and rs1891467 (TGFB2), which additionally showed significant differences in genotype distribution among European subpopulations. This analysis also showed statistically significant differences in genotype distributions between healthy controls and PVR patients in rs17561 of the IL1A gene (OR, 3.00; 95% CI, 0.77-11.75; P = 0.036) and in rs1800629 of the TNF gene (OR, 0.48; 95% CI, 0.27-0.87; P = 0.014). Furthermore, we have shown that a small change (0.02) in minor allele frequency (MAF) significantly affects the statistical p value in case-control studies. In conclusion, the study showed differences in genotype distributions in healthy populations across different European countries. Differences in distribution of genotypes may have had influenced failed replication results in previous PVR-related SNP-association studies.

Genome-wide and species-wide in silico screening for intragenic MicroRNAs in human, mouse and chicken.

MicroRNAs (miRNAs) are non-coding RNAs (ncRNAs) involved in regulation of gene expression. Intragenic miRNAs, especially those exhibiting a high degree of evolutionary conservation, have been shown to be coordinately regulated and/or expressed with their host genes, either with synergistic or antagonistic correlation patterns. However, the degree of cross-species conservation of miRNA/host gene co-location is not known and co-expression information is incomplete and fragmented among several studies. Using the genomic resources (miRBase and Ensembl) we performed a genome-wide in silico screening (GWISS) for miRNA/host gene pairs in three well-annotated vertebrate species: human, mouse, and chicken. Approximately half of currently annotated miRNA genes resided within host genes: 53.0% (849/1,600) in human, 48.8% (418/855) in mouse, and 42.0% (210/499) in chicken, which we present in a central publicly available Catalog of intragenic miRNAs (http://www.integratomics-time.com/miR-host/catalog). The miRNA genes resided within either protein-coding or ncRNA genes, which include long intergenic ncRNAs (lincRNAs) and small nucleolar RNAs (snoRNAs). Twenty-seven miRNA genes were found to be located within the same host genes in all three species and the data integration from literature and databases showed that most (26/27) have been found to be co-expressed. Particularly interesting are miRNA genes located within genes encoding for miRNA silencing machinery (DGCR8, DICER1, and SND1 in human and Cnot3, Gdcr8, Eif4e, Tnrc6b, and Xpo5 in mouse). We furthermore discuss a potential for phenotype misattribution of miRNA host gene polymorphism or gene modification studies due to possible collateral effects on miRNAs hosted within them. In conclusion, the catalog of intragenic miRNAs and identified 27 miRNA/host gene pairs with cross-species conserved co-location, co-expression, and potential co-regulation, provide excellent candidates for further functional annotation of intragenic miRNAs in health and disease.

Obesity gene atlas in mammals.

Obesity in humans has increased at an alarming rate over the past two decades and has become one of the leading public health problems worldwide. Studies have revealed a large number of genes/markers that are associated with obesity and/or obesity-related phenotypes, indicating an urgent need to develop a central database for helping the community understand the genetic complexity of obesity. In the present study, we collected a total of 1,736 obesity associated loci and created a freely available obesity database, including 1,515 protein-coding genes and 221 microRNAs (miRNAs) collected from four mammalian species: human, cattle, rat, and mouse. These loci were integrated as orthologs on comparative genomic views in human, cattle, and mouse. The database and genomic views are freely available online at: http://www.integratomics-time.com/fat_deposition. Bioinformatics analyses of the collected data revealed some potential novel obesity related molecular markers which represent focal points for testing more targeted hypotheses and designing experiments for further studies. We believe that this centralized database on obesity and adipogenesis will facilitate development of comparative systems biology approaches to address this important health issue in human and their potential applications in animals.