miRNome Profiling Analysis Reveals Novel Hepatocellular Carcinoma Diagnostic, Prognostic and Treatment-Related Candidate Biomarkers: Post hoc Analysis of SORAMIC Trial.

INTRODUCTION: Early diagnosis of hepatocellular carcinoma (HCC) as well as evaluation of prognosis and prediction of treatment efficacy remains challenging due to the missing specific non-invasive biomarkers. The aim of this study was to identify disease-specific microRNA (miRNA) patterns for diagnosis, prediction of prognosis, and treatment response in patients with HCC. METHODS: The study population included 42 HCC patients from SORAMIC clinical trial: 22 patients received sorafenib monotherapy, 20 patients underwent 90Y radioembolization in combination with sorafenib. 20 individuals were included in the control group. HCC patients underwent collection of plasma samples before and 7-9 weeks after the beginning of the treatment. Isolation of circulating miRNAs, preparation of small RNA sequencing libraries and next-generation sequencing were performed. Association analysis for novel diagnostic, prognostic, and treatment-related candidate biomarkers was performed. RESULTS: A total of 42 differentially expressed (16 up-regulated and 26 down-regulated) miRNAs were identified comparing baseline and control group plasma samples. hsa-miR-215-5p and hsa-miR-192-5p were down-regulated, while hsa-miR-483-5p and hsa-miR-23b-3p were up-regulated comparing baseline and 7-9 weeks post-sorafenib monotherapy samples. hsa-miR-215-5p was the sole down-regulated miRNA in the same combination therapy comparison. hsa-miR-183-5p, hsa-miR-28-3p, and hsa-miR-1246 were found to be significantly up-regulated comparing non-responders versus responders to sorafenib. High hsa-miR-215-5p expression was significantly associated with worse HCC patients' prognosis. CONCLUSIONS: Systematic miRNA profiling of highly characterized samples from SORAMIC study revealed a subset of potential miRNA biomarkers for HCC diagnosis and prognosis of sorafenib-treated patients' survival.

5'-Isoforms of miR-1246 Have Distinct Targets and Stronger Functional Impact Compared with Canonical miR-1246 in Colorectal Cancer Cells In Vitro.

Colorectal cancer (CRC) is a multifactorial disease involving genetic and epigenetic factors, such as miRNAs. Sequencing-based studies have revealed that miRNAs have many isoforms (isomiRs) with modifications at the 3'- and 5'-ends or in the middle, resulting in distinct targetomes and, consequently, functions. In the present study, we aimed to evaluate the putative targets and functional role of miR-1246 and its two 5'-isoforms (ISO-miR-1246_a and ISO-miR-1246_G) in vitro. Commercial Caco-2 cells of CRC origin were analyzed for the expression of WT-miR-1246 and its 5'-isoforms using small RNA sequencing data, and the overabundance of the two miR-1246 isoforms was determined in cells. The transcriptome analysis of Caco-2 cells transfected with WT-miR-1246, ISO-miR-1246_G, and ISO-miR-1246_a indicated the minor overlap of the targetomes between the studied miRNA isoforms. Consequently, an enrichment analysis showed the involvement of the potential targets of the miR-1246 isoforms in distinct signaling pathways. Cancer-related pathways were predominantly more enriched in dysregulated genes in ISO-miR-1246_G and ISO-miR-1246_a, whereas cell cycle pathways were more enriched in WT-miR-1246. The functional analysis of WT-miR-1246 and its two 5'-isoforms revealed that the inhibition of any of these molecules had a tumor-suppressive role (reduced cell viability and migration and promotion of early cell apoptosis) in CRC cells. However, the 5'-isoforms had a stronger effect on viability compared with WT-miR-1246. To conclude, this research shows that WT-miR-1246 and its two 5'-isoforms have different targetomes and are involved in distinct signaling pathways but collectively play an important role in CRC pathogenesis.

Diagnostic and Prognostic Value of IL-10, FABP2 and LPS Levels in HCC Patients.

Hepatocellular carcinoma (HCC) still lacks valuable diagnostic and prognostic tools. This study aimed to investigate the potential diagnostic and prognostic value of baseline interleukin (IL)-10, fatty acid-binding protein 2 (FABP2) and lipopolysaccharide (LPS) levels in patients with HCC. Serum levels of IL-10, FABP2 and LPS in 47 newly diagnosed HCC patients and 50 healthy individuals were estimated and compared. The best cut-off points for baseline IL-10, FABP2 and LPS levels predicting overall survival (OS) were evaluated. Both levels of FABP2 and IL-10 were significantly higher in HCC patients vs. control group (median 2095 vs. 1772 pg/mL, p = 0.026; 9.94 vs. 4.89 pg/mL, p < 0.001) and may serve as potential biomarkers in complex HCC diagnostic tools. The cut-off value of 2479 pg/mL for FABP2 was determined to have the highest sensitivity (66.7%) and specificity (55.6%) to distinguish patients with a median OS longer than 17 months. However, the median OS of patients with high and low levels of FABP2 were not significantly different (p = 0.896). The prognostic value of LPS as well as FABP2 and IL-10 for HCC patients appears to be limited.

The microRNA expression in crypt-top and crypt-bottom colonic epithelial cell populations demonstrates cell-type specificity and correlates with endoscopic activity in ulcerative colitis.

BACKGROUND AND AIMS: Colonic epithelial barrier dysfunction is one of the early events in ulcerative colitis (UC) and microRNAs (miRNAs) participate in its regulation. However, cell type-specific miRNome during UC is still unknown. Thus, we aimed to explore miRNA expression patterns in colon tissue and epithelial cells at active and quiescent UC. METHODS: Small RNA-sequencing in colon tissue, crypt-bottom (CD44+), and crypt-top (CD66a+) colonic epithelial cells from two cohorts of UC patients (n=74) and healthy individuals (n=50) was performed. Data analysis encompassed differential expression, weighted gene co-expression network, correlation, gene-set enrichment analyses. RESULTS: Differentially expressed colonic tissue miRNAs showed potential involvement in regulation of interleukin-4 and interleukin-13 signalling during UC. As this pathway plays role in intestinal barrier regulation, consecutive analysis of spatially distinct colonic epithelial cell populations was performed. Cell-type (crypt-top and crypt-bottom) specific miRNA expression patterns were identified in both active and quiescent UC. Target genes of differentially expressed epithelial miRNAs at different disease activity were overrepresented in epithelial cell migration and therefore intestinal barrier integrity regulation. The pro-inflammatory miRNA co-expression module M1 correlated with endoscopic disease activity and successfully distinguished active and quiescent UC not only in both epithelial cell populations, but also in the colon tissue. The anti-inflammatory module M2 was specific to crypt-bottom cells and significantly enriched in the quiescent UC patients. CONCLUSIONS: miRNA expression was specific to colonic epithelial cell populations and UC state, reflecting endoscopic disease activity. Irrespective of the UC state, deregulated epithelial miRNAs were associated with regulation of intestinal barrier integrity.

Circulating microbiome in patients with portal hypertension.

Portal hypertension (PH) in liver cirrhosis leads to increased gut permeability and the translocation of bacteria across the gut-liver axis. Microbial DNA has recently been detected in different blood compartments; however, this phenomenon has not been thoroughly analyzed in PH. This study aimed to explore circulating bacterial DNA signatures, inflammatory cytokines, and gut permeability markers in different blood compartments (peripheral and hepatic veins) of patients with cirrhosis and PH. The 16S rRNA blood microbiome profiles were determined in 58 patients with liver cirrhosis and 46 control patients. Taxonomic differences were analyzed in relation to PH, liver function, inflammatory cytokines, and gut permeability markers. Circulating plasma microbiome profiles in patients with cirrhosis were distinct from those of the controls and were characterized by enrichment of Comamonas, Cnuella, Dialister, Escherichia/Shigella, and Prevotella and the depletion of Bradyrhizobium, Curvibacter, Diaphorobacter, Pseudarcicella, and Pseudomonas. Comparison of peripheral and hepatic vein blood compartments of patients with cirrhosis did not reveal differentially abundant taxa. Enrichment of the genera Bacteroides, Escherichia/Shigella, and Prevotella was associated with severe PH (SPH) in both blood compartments; however, circulating microbiome profiles could not predict PH severity. Escherichia/Shigella and Prevotella abundance was correlated with IL-8 levels in the hepatic vein. In conclusion, we demonstrated a distinct circulating blood microbiome profile in patients with cirrhosis, showing that specific bacterial genera in blood are marginally associated with SPH, Model for End-Stage Liver Disease score, and inflammation biomarkers; however, circulating microbial composition failed to predict PH severity.

High-throughput method for the hybridisation-based targeted enrichment of long genomic fragments for PacBio third-generation sequencing.

Hybridisation-based targeted enrichment is a widely used and well-established technique in high-throughput second-generation short-read sequencing. Despite the high potential to genetically resolve highly repetitive and variable genomic sequences by, for example PacBio third-generation sequencing, targeted enrichment for long fragments has not yet established the same high-throughput due to currently existing complex workflows and technological dependencies. We here describe a scalable targeted enrichment protocol for fragment sizes of >7 kb. For demonstration purposes we developed a custom blood group panel of challenging loci. Test results achieved > 65% on-target rate, good coverage (142.7×) and sufficient coverage evenness for both non-paralogous and paralogous targets, and sufficient non-duplicate read counts (83.5%) per sample for a highly multiplexed enrichment pool of 16 samples. We genotyped the blood groups of nine patients employing highly accurate phased assemblies at an allelic resolution that match reference blood group allele calls determined by SNP array and NGS genotyping. Seven Genome-in-a-Bottle reference samples achieved high recall (96%) and precision (99%) rates. Mendelian error rates were 0.04% and 0.13% for the included Ashkenazim and Han Chinese trios, respectively. In summary, we provide a protocol and first example for accurate targeted long-read sequencing that can be used in a high-throughput fashion.

Gut microbial similarity in twins is driven by shared environment and aging.

BACKGROUND: Human gut microbiome composition is influenced by genetics, diet and environmental factors. We investigated the microbial composition in several gastrointestinal (GI) compartments to evaluate the impact of genetics, delivery mode, diet, household sharing and aging on microbial similarity in monozygotic and dizygotic twins. METHODS: Fecal, biopsy and saliva samples were obtained from total 108 twins. DNA and/or RNA was extracted and the region V1-V2 of the 16S rRNA gene was amplified and sequenced. Bray-Curtis similarity was used for further microbiome comparisons, Mann-Whitney test was applied to evaluate the significant differences between groups and Spearman test was applied to reveal potential correlations between data. FINDINGS: The global bacterial profiles were grouped into two clusters separating the upper and lower GI. The upper GI microbiome composition was strictly dependent on the Helicobacter pylori status. With a positivity rate of 55%, H. pylori completely colonized the stomach and separated infected twins from non-infected twins irrespective of zygosity status. Lower GI microbiome similarity between the twins was defined mainly by household-sharing and aging; whereas delivery mode and host genetics had no influence. There was a progredient decrease in the bacterial similarity with aging. Shared vs. non-shared phylotypes analysis showed that in both siblings the shared phylotypes progressively diminished with aging, while the non-shared phylotypes increased. INTERPRETATION: Our findings strongly highlight the aging and shared household as they key determinants in gut microbial similarity and drift in twins irrespective of their zygotic state. FUNDING: This work was supported by the grant of the Research Council of Lithuania (Project no. APP-2/2016) and also partially supported by the funds of European Commission through the "European funds for regional development" (EFRE) as well as by the regional Ministry of Economy, Science and Digitalization as part of the "LiLife" Project as part of the "Autonomy in old Age" research group (Project ID: ZS/2018/11/95324).

Atrophic gastritis and gastric cancer tissue miRNome analysis reveals hsa-miR-129-1 and hsa-miR-196a as potential early diagnostic biomarkers.

BACKGROUND: Gastric cancer (GC) is one of the most frequently diagnosed tumor globally. In most cases, GC develops in a stepwise manner from chronic gastritis or atrophic gastritis (AG) to cancer. One of the major issues in clinical settings of GC is diagnosis at advanced disease stages resulting in poor prognosis. MicroRNAs (miRNAs) are small noncoding molecules that play an essential role in a variety of fundamental biological processes. However, clinical potential of miRNA profiling in the gastric cancerogenesis, especially in premalignant GC cases, remains unclear. AIM: To evaluate the AG and GC tissue miRNomes and identify specific miRNAs' potential for clinical applications (e.g., non-invasive diagnostics). METHODS: Study included a total of 125 subjects: Controls (CON), AG, and GC patients. All study subjects were recruited at the Departments of Surgery or Gastroenterology, Hospital of Lithuanian University of Health Sciences and divided into the profiling (n = 60) and validation (n = 65) cohorts. Total RNA isolated from tissue samples was used for preparation of small RNA sequencing libraries and profiled using next-generation sequencing (NGS). Based on NGS data, deregulated miRNAs hsa-miR-129-1-3p and hsa-miR-196a-5p were analyzed in plasma samples of independent cohort consisting of CON, AG, and GC patients. Expression level of hsa-miR-129-1-3p and hsa-miR-196a-5p was determined using the quantitative real-time polymerase chain reaction and 2-ΔΔCt method. RESULTS: Results of tissue analysis revealed 20 differentially expressed miRNAs in AG group compared to CON group, 129 deregulated miRNAs in GC compared to CON, and 99 altered miRNAs comparing GC and AG groups. Only 2 miRNAs (hsa-miR-129-1-3p and hsa-miR-196a-5p) were identified to be step-wise deregulated in healthy-premalignant-malignant sequence. Area under the curve (AUC)-receiver operating characteristic analysis revealed that expression level of hsa-miR-196a-5p is significant for discrimination of CON vs AG, CON vs GC and AG vs GC and resulted in AUCs: 88.0%, 93.1% and 66.3%, respectively. Compar-ing results in tissue and plasma samples, hsa-miR-129-1-3p was significantly down-regulated in GC compared to AG (P = 0.0021 and P = 0.024, tissue and plasma, respectively). Moreover, analysis revealed that hsa-miR-215-3p/5p and hsa-miR-934 were significantly deregulated in GC based on Helicobacter pylori (H. pylori) infection status [log2 fold change (FC) = -4.52, P-adjusted = 0.02; log2FC = -4.00, P-adjusted = 0.02; log2FC = 6.09, P-adjusted = 0.02, respectively]. CONCLUSION: Comprehensive miRNome study provides evidence for gradual deregulation of hsa-miR-196a-5p and hsa-miR-129-1-3p in gastric carcinogenesis and found hsa-miR-215-3p/5p and hsa-miR-934 to be significantly deregulated in H. pylori carrying GC patients.

Liquid Biopsy in Gastric Cancer: Analysis of Somatic Cancer Tissue Mutations in Plasma Cell-Free DNA for Predicting Disease State and Patient Survival.

INTRODUCTION: Gastric cancer (GC) diagnosis in late stages and high mortality rates are the main issues that require new noninvasive molecular tools. We aimed to assess somatic mutational profiles in GC tissue and plasma cell-free DNA (cfDNA), evaluate their concordance rate, and analyze the role of multilayer molecular profiling to predict disease state and prognosis. METHODS: Treatment-naive GC patient group (n = 29) was selected. Whole exome sequencing (WES) of GC tissue was performed, and a unique 38-gene panel for deep targeted sequencing of plasma cfDNA was developed. Oncoproteins were measured by enzyme-linked immunosorbent assay, and other variables such as tumor mutational burden and microsatellite instability were evaluated using WES data. RESULTS: The yield of cfDNA was increased 43.6-fold; the integrity of fragments was decreased in GC compared with controls. WES analysis of cancerous tissue and plasma cfDNA (targeted sequencing) mutational profiles revealed 47.8% concordance. The increased quantity of GC tissue-derived alterations detected in cfDNA was associated with worse patients' survival. Analysis of importance of multilayer variables and receiver operating characteristic curve showed that combination of 2 analytes: (i) quantity of tissue matching alterations and (ii) presence of any somatic alteration in plasma cfDNA resulted in area under curve 0.744 when discriminating patients with or without distant metastasis. Furthermore, cfDNA sequence alterations derived from tumor tissue were detected in patients who had even relatively small GC tumors (T1-T2). DISCUSSION: Our results indicate that quantitative and qualitative cfDNA mutational profile analysis is a promising tool for evaluating GC disease status or poorer prognosis.

PNPLA3 and RNF7 Gene Variants are Associated with the Risk of Developing Liver Fibrosis and Cirrhosis in an Eastern European Population.

BACKGROUND AND AIMS: Genome-wide association studies have revealed an association between the risk of developing liver fibrosis or cirrhosis and the single nucleotide polymorphisms (SNPs) of the PNPLA3, RNF7, MERTK and PCSK7 genes. We aimed to validate these results in an Eastern European population. METHODS: We evaluated the associations between the PNPLA3 (rs738409), RNF7 (rs16851720), MERTK (rs4374383) and PCSK7 (rs236918) variants and liver fibrosis and cirrhosis in a series of consecutive patients recruited at the Department of Gastroenterology, Lithuanian University of Health Sciences Hospital, during the period 2012-2015. The study included 317 individuals with liver cirrhosis, 154 individuals with liver fibrosis, and 498 controls. The studied SNPs were determined using RT-PCR TaqMan assays. RESULTS: MERTK and PCSK7 SNPs were not associated with liver fibrosis or cirrhosis. The PNPLA3 SNP rs738409 was associated with a higher risk of developing liver fibrosis (aOR: 1.65, P=0.001) and cirrhosis (aOR: 1.92, P=5.57*10-7). PNPLA3 genotypes were also associated with higher risk of developing liver fibrosis and cirrhosis in dominant (aOR: 1.98, P=2.20*10-5; aOR: 1.67, P=0.008, respectively) and recessive (aOR: 3.94, P=5.16*10-5; aOR: 3.02, P=0.003, respectively) models. RNF7 rs16851720 was associated with liver cirrhosis comparing CC vs. AA + CA genotypes (aOR: 0.26, P=0.020). CONCLUSION: Our study showed that PNPLA3 rs738409 and RNF7 rs16851720 confer an increased risk of developing liver fibrosis and cirrhosis in this Eastern European population, while the MERTK and PCSK7 SNPs are not associated with these conditions.

Thrombosis Related ABO, F5, MTHFR, and FGG Gene Polymorphisms in Morbidly Obese Patients.

Objective. Obesity is a well-known risk factor for thrombotic complications. The aim of the present study was to determine the frequency of thrombosis related ABO, F5, MTHFR, and FGG gene polymorphisms in morbidly obese patients and compare them with the group of nonobese individuals. Methods. Gene polymorphisms were analyzed in 320 morbidly obese patients (BMI > 40 kg/m2) and 303 control individuals (BMI < 30 kg/m2) of European descent. ABO C>T (rs505922), F5 C>G (rs6427196), MTHFR C>T (rs1801133), and FGG C>T (rs6536024) SNPs were genotyped by RT-PCR. Results. We observed a tendency for MTHFR rs1801133 TT genotype to be linked with morbid obesity when compared to CC genotype; however, the difference did not reach the significant P value (OR 1.84, 95% CI 0.83-4.05, P = 0.129). Overall, the genotypes and alleles of rs505922, rs6427196, rs1801133, and rs6536024 SNPs had similar distribution between morbidly obese and nonobese control individuals. Distribution of height and weight means among individuals carrying different rs505922, rs6427196, rs1801133, and rs6536024 genotypes did not differ significantly. Conclusions. Gene polymorphisms ABO C>T (rs505922), F5 C>G (rs6427196), MTHFR C>T (rs1801133), and FGG C>T (rs6536024) were not associated with height, weight, or morbid obesity among European subjects.

TPMT and ITPA genetic variants in Lithuanian inflammatory bowel disease patients: Prevalence and azathioprine-related side effects.

PURPOSE: Inter-individual thiopurine metabolism variability can influence treatment outcomes in inflammatory bowel disease (IBD) patients. Genetic polymorphisms in thiopurine methyltransferase (TPMT) and inosine triphosphate pyrophosphatase (ITPA) were linked with toxicity of azathioprine (AZA). The aim of the study was to investigate frequencies of TPMT and ITPA polymorphisms in Lithuanian IBD patients and analyze their association with AZA-related adverse events. MATERIALS/METHODS: Polymorphisms in TPMT (TPMT*2,*3B,*3C,*3A) and ITPA (rs1127354, rs7270101) genes were determined using PCR-RFLP and TaqMan(®) genotyping assays. 551 consecutive Lithuanian IBD patients were genotyped. The use of AZA and its side effects were assessed retrospectively according to hospital medical records. RESULTS: Frequencies of TPMT*3A, TPMT*3B and TPMT*3C alleles were 3.1%, 0.5% and 0.1%, respectively. TPMT*2 genetic variant was not detected in the study group. The distribution of minor alleles for ITPA rs1127354 and rs7270101 polymorphisms was 9.9% and 10.5%, respectively. AZA was prescribed in 82 patients and it provoked myelotoxicity in 11%, hepatotoxicity in 6.1%, dyspepsia in 6.1%, and pancreatitis in 3.6% of cases. Among patients who had AZA-related myelotoxicity, 11.1% were TPMT compound heterozygous, 44.4% had heterozygous genotype (P<0.01). Frequencies of ITPA minor alleles were similar among the patients with and without AZA-related side effects. CONCLUSION: Frequencies of TPMT and ITPA variant alleles in Lithuanian IBD group were similar to those observed in the Northern-Eastern Europe Caucasian populations. Polymorphisms in TPMT might be associated with myelotoxicity and leukopenia in AZA treated patients, while ITPA variant alleles appear not to be linked with treatment-related side effects.