Search

Background: The epigenome, the set of modifications to DNA and associated molecules that control gene expression, cellular identity, and function, plays a major role in mediating cellular responses to outside factors. Thus, evaluation of the epigenetic state can provide insights into cellular adaptions occurring over the course of disease., Methods: We performed epigenome-wide association studies of primary sclerosing cholangitis (PSC) and primary biliary cholangitis (PBC) using the Illumina MethylationEPIC Bead Chip., Results: We found evidence of increased epigenetic age acceleration and differences in predicted immune cell composition in patients with PSC and PBC. Epigenetic profiles demonstrated differences in predicted protein levels including increased levels of tumor necrosis factor receptor superfamily member 1B in patients with cirrhotic compared to noncirrhotic PSC and PBC. Epigenome-wide association studies of PSC discovered strongly associated 5'-C-phosphate-G-3' sites in genes including vacuole membrane protein 1 and SOCS3, and epigenome-wide association studies of PBC found strong 5'-C-phosphate-G-3' associations in genes including NOD-like receptor family CARD domain containing 5, human leukocyte antigen-E, and PSMB8. Analyses identified disease-associated canonical pathways and upstream regulators involved with immune signaling and activation of macrophages and T-cells. A comparison of PSC and PBC data found relatively little overlap at the 5'-C-phosphate-G-3' and gene levels with slightly more overlap at the level of pathways and upstream regulators., Conclusions: This study provides insights into methylation profiles of patients that support current concepts of disease mechanisms and provide novel data to inspire future research. Studies to corroborate our findings and expand into other -omics layers will be invaluable to further our understanding of these rare diseases with the goal to improve and individualize prognosis and treatment., (Copyright © 2024 The Author(s). Published by Wolters Kluwer Health, Inc. on behalf of the American Association for the Study of Liver Diseases.)

Our objective was to test whether p53 expression status is associated with survival for women diagnosed with the most common ovarian carcinoma histotypes (high-grade serous carcinoma [HGSC], endometrioid carcinoma [EC], and clear cell carcinoma [CCC]) using a large multi-institutional cohort from the Ovarian Tumor Tissue Analysis (OTTA) consortium. p53 expression was assessed on 6,678 cases represented on tissue microarrays from 25 participating OTTA study sites using a previously validated immunohistochemical (IHC) assay as a surrogate for the presence and functional effect of TP53 mutations. Three abnormal expression patterns (overexpression, complete absence, and cytoplasmic) and the normal (wild type) pattern were recorded. Survival analyses were performed by histotype. The frequency of abnormal p53 expression was 93.4% (4,630/4,957) in HGSC compared to 11.9% (116/973) in EC and 11.5% (86/748) in CCC. In HGSC, there were no differences in overall survival across the abnormal p53 expression patterns. However, in EC and CCC, abnormal p53 expression was associated with an increased risk of death for women diagnosed with EC in multivariate analysis compared to normal p53 as the reference (hazard ratio [HR] = 2.18, 95% confidence interval [CI] 1.36-3.47, p = 0.0011) and with CCC (HR = 1.57, 95% CI 1.11-2.22, p = 0.012). Abnormal p53 was also associated with shorter overall survival in The International Federation of Gynecology and Obstetrics stage I/II EC and CCC. Our study provides further evidence that functional groups of TP53 mutations assessed by abnormal surrogate p53 IHC patterns are not associated with survival in HGSC. In contrast, we validate that abnormal p53 IHC is a strong independent prognostic marker for EC and demonstrate for the first time an independent prognostic association of abnormal p53 IHC with overall survival in patients with CCC., (© 2023 The Authors. The Journal of Pathology: Clinical Research published by The Pathological Society of Great Britain and Ireland and John Wiley & Sons Ltd.)

Background: Primary sclerosing cholangitis (PSC) patients have a risk of developing cholangiocarcinoma (CCA). Establishing predictive models for CCA in PSC is important., Methods: In a large cohort of 1,459 PSC patients seen at Mayo Clinic (1993-2020), we quantified the impact of clinical/laboratory variables on CCA development using univariate and multivariate Cox models and predicted CCA using statistical and artificial intelligence (AI) approaches. We explored plasma bile acid (BA) levels' predictive power of CCA (subset of 300 patients, BA cohort)., Results: Eight significant risk factors (false discovery rate: 20%) were identified with univariate analysis; prolonged inflammatory bowel disease (IBD) was the most important one. IBD duration, PSC duration, and total bilirubin remained significant (p < 0.05) with multivariate analysis. Clinical/laboratory variables predicted CCA with cross-validated C-indexes of 0.68-0.71 at different time points of disease, significantly better compared to commonly used PSC risk scores. Lower chenodeoxycholic acid, higher conjugated fraction of lithocholic acid and hyodeoxycholic acid, and higher ratio of cholic acid to chenodeoxycholic acid were predictive of CCA. BAs predicted CCA with a cross-validated C-index of 0.66 (std: 0.11, BA cohort), similar to clinical/laboratory variables (C-index = 0.64, std: 0.11, BA cohort). Combining BAs with clinical/laboratory variables leads to the best average C-index of 0.67 (std: 0.13, BA cohort)., Conclusions: In a large PSC cohort, we identified clinical and laboratory risk factors for CCA development and demonstrated the first AI based predictive models that performed significantly better than commonly used PSC risk scores. More predictive data modalities are needed for clinical adoption of these models., (© 2023. The Author(s).)

Background: Better understanding of prognostic factors in tubo-ovarian high-grade serous carcinoma (HGSC) is critical, as diagnosis confers an aggressive disease course. Variation in tumor DNA methylation shows promise predicting outcome, yet prior studies were largely platform-specific and unable to evaluate multiple molecular features., Methods: We analyzed genome-wide DNA methylation in 1,040 frozen HGSC, including 325 previously reported upon, seeking a multi-platform quantitative methylation signature that we evaluated in relation to clinical features, tumor characteristics, time to recurrence/death, extent of CD8+ tumor-infiltrating lymphocytes (TIL), gene expression molecular subtypes, and gene expression of the ATP-binding cassette transporter TAP1., Results: Methylation signature was associated with shorter time to recurrence, independent of clinical factors (N = 715 new set, hazard ratio (HR), 1.65; 95% confidence interval (CI), 1.10-2.46; P = 0.015; N = 325 published set HR, 2.87; 95% CI, 2.17-3.81; P = 2.2 × 10-13) and remained prognostic after adjustment for gene expression molecular subtype and TAP1 expression (N = 599; HR, 2.22; 95% CI, 1.66-2.95; P = 4.1 × 10-8). Methylation signature was inversely related to CD8+ TIL levels (P = 2.4 × 10-7) and TAP1 expression (P = 0.0011) and was associated with gene expression molecular subtype (P = 5.9 × 10-4) in covariate-adjusted analysis., Conclusions: Multi-center analysis identified a novel quantitative tumor methylation signature of HGSC applicable to numerous commercially available platforms indicative of shorter time to recurrence/death, adjusting for other factors. Along with immune cell composition analysis, these results suggest a role for DNA methylation in the immunosuppressive microenvironment., Impact: This work aids in identification of targetable epigenome processes and stratification of patients for whom tailored treatment may be most beneficial., (©2023 American Association for Cancer Research.)

Aims: In this methylome-wide association study of cholestatic liver diseases (primary sclerosing cholangitis and primary biliary cholangitis), the authors aimed to elucidate changes in methylome and pathway enrichment to identify candidate genes . Patients & methods: Reduced representation bisulfite sequencing was performed on liver tissue from 58 patients with primary sclerosing cholangitis (n = 13), primary biliary cholangitis (n = 20), alcoholic liver disease (n = 21) and live liver donors (n = 4). Pathway enrichment and network analysis were used to explore key genes/pathways. Results: Both cholestatic liver diseases were characterized by global hypomethylation, with pathway enrichment demonstrating distinct genes and pathways associated with the methylome. Conclusions: This novel study demonstrated that differential methylation in cholestatic liver disease was associated with unique pathways, suggesting it may drive disease pathogenesis.

Background: Ovarian clear cell carcinoma (OCCC) is a rare ovarian cancer histotype that tends to be resistant to standard platinum-based chemotherapeutics. We sought to better understand the role of DNA methylation in clinical and biological subclassification of OCCC., Methods: We interrogated genome-wide methylation using DNA from fresh frozen tumors from 271 cases, applied nonsmooth nonnegative matrix factorization (nsNMF) clustering, and evaluated clinical associations and biological pathways., Results: Two approximately equally sized clusters that associated with several clinical features were identified. Compared with Cluster 2 ( N = 137), Cluster 1 cases ( N = 134) presented at a more advanced stage, were less likely to be of Asian ancestry, and tended to have poorer outcomes including macroscopic residual disease following primary debulking surgery ( P < 0.10). Subset analyses of targeted tumor sequencing and IHC data revealed that Cluster 1 tumors showed TP53 mutation and abnormal p53 expression, and Cluster 2 tumors showed aneuploidy and ARID1A/PIK3CA mutation ( P < 0.05). Cluster-defining CpGs included 1,388 CpGs residing within 200 bp of the transcription start sites of 977 genes; 38% of these genes ( N = 369 genes) were differentially expressed across cluster in transcriptomic subset analysis ( P < 10 -4 ). Differentially expressed genes were enriched for six immune-related pathways, including IFNα and IFNγ responses ( P < 10 -6 )., Conclusions: DNA methylation clusters in OCCC correlate with disease features and gene expression patterns among immune pathways., Impact: This work serves as a foundation for integrative analyses that better understand the complex biology of OCCC in an effort to improve potential for development of targeted therapeutics., (©2021 American Association for Cancer Research.)

Rapid cycling (RC) burdens bipolar disorder (BD) patients further by causing more severe disability and increased suicidality. Because diagnosing RC can be challenging, RC patients are at risk of rapid decline due to delayed suitable treatment. Here, we aimed to identify the differences in the circulating cell-free DNA (cfDNA) methylome between BD patients with and without RC. The cfDNA methylome could potentially be developed as a diagnostic test for BD RC. We extracted cfDNA from plasma samples of BD1 patients (46 RC and 47 non-RC). cfDNA methylation levels were measured by 850K Infinium MethylationEPIC array. Principal component analysis (PCA) was conducted to assess global differences in methylome. cfDNA methylation levels were compared between RC groups using a linear model adjusted for age and sex. PCA suggested differences in methylation profiles between RC groups ( p = 0.039) although no significant differentially methylated probes (DMPs; q > 0.15) were found. The top four CpG sites which differed between groups at p < 1E-05 were located in CGGPB1 , PEX10 , NR0B2 , and TP53I11 . Gene set enrichment analysis (GSEA) on top DMPs ( p < 0.05) showed significant enrichment of gene sets related to nervous system tissues, such as neurons, synapse, and glutamate neurotransmission. Other top notable gene sets were related to parathyroid regulation and calcium signaling. To conclude, our study demonstrated the feasibility of utilizing a microarray method to identify circulating cfDNA methylation sites associated with BD RC and found the top differentially methylated CpG sites were mostly related to the nervous system and the parathyroid., Competing Interests: MF was a consultant (for Mayo Clinic) to Janssen, Mitsubishi Tanabe Pharma Corporation, Myriad, Sunovion, and Teva Pharmaceuticals; none of this funding contributed to any work carried out in this study. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Ho, Winham, McCauley, Kundakovic, Robertson, Sun, Ordog, Webb, Frye and Veldic.)

Background: Primary sclerosing cholangitis (PSC) is a rare, chronic cholestatic liver disease that often progresses to end-stage liver disease and/or the development of hepatobiliary neoplasia. Lack of prognostic tools and treatment options for PSC is driven in part by our poor understanding of its pathogenesis, which is thought to be complex, the interaction of genetic variants, environmental influences and biological response throughout the course of disease. The PSC Scientific Community Resource (PSC-SCR) seeks to overcome previous shortcomings by facilitating novel research in PSC with the ultimate goals of individualizing patient care and improving patient outcomes., Methods: PSC patients who receive their health care at Mayo Clinic or a collaborating site are identified by chart review and invited in person or by mail to participate. Non-Mayo patients are offered enrollment if they provide sufficient access to their medical records to evaluate inclusion/exclusion criteria. Controls without liver disease are identified with assistance of the Mayo Clinic Biobank. Participant consent is obtained at the beginning of the recruitment process by mail-in, electronic or face-to-face protocols. Clinical data is extracted from the medical record by qualified physicians and entered in a custom designed database. Participants fill out a custom-designed, comprehensive questionnaire, which collects scientifically relevant demographic and clinical information. Biospecimens are collected using mail-in kits thar are returned via overnight carrier service and processed by the biospecimen accessioning and processing facility at Mayo Clinic, which coordinates sample transfers and provides required sample preparation services. The resource is currently being utilized to perform omics-scale projects investigating the exposome, metabolome, methylome, immunome and microbiome in PSC. Datasets and residual biospecimens will be shared with researchers proposing scientifically sound PSC-focused research with approval of the appropriate review boards., Discussion: Patient-based studies leveraging the latest technologies for targeted and wide-scale interrogation of multiple omics layers offer promise to accelerate PSC research through discovery of unappreciated aspects of disease pathogenesis. However, the rarity of PSC severely limits such studies. Here we describe our effort to overcome this limitation, the PSC-SCR, a repository of patient biospecimens coupled with clinical and omics data for use by the broader PSC research community., (© 2021. The Author(s).)

Background: Many loci have been found to be associated with risk of epithelial ovarian cancer (EOC). However, although there is considerable variation in progression-free survival (PFS), no loci have been found to be associated with outcome at genome-wide levels of significance., Methods: We carried out a genome-wide association study (GWAS) of PFS in 2,352 women with EOC who had undergone cytoreductive surgery and standard carboplatin/paclitaxel chemotherapy., Results: We found seven SNPs at 12q24.33 associated with PFS ( P < 5 × 10 -8 ), the top SNP being rs10794418 (HR = 1.24; 95% CI, 1.15-1.34; P = 1.47 × 10 -8 ). High expression of a nearby gene, ULK1 , is associated with shorter PFS in EOC, and with poor prognosis in other cancers. SNP rs10794418 is also associated with expression of ULK1 in ovarian tumors, with the allele associated with shorter PFS being associated with higher expression, and chromatin interactions were detected between the ULK1 promoter and associated SNPs in serous and endometrioid EOC cell lines. ULK1 knockout ovarian cancer cell lines showed significantly increased sensitivity to carboplatin in vitro ., Conclusions: The locus at 12q24.33 represents one of the first genome-wide significant loci for survival for any cancer. ULK1 is a plausible candidate for the target of this association., Impact: This finding provides insight into genetic markers associated with EOC outcome and potential treatment options. See related commentary by Peres and Monteiro, p. 1604 ., (©2021 American Association for Cancer Research.)

Background: PTEN loss is a putative driver in histotypes of ovarian cancer (high-grade serous (HGSOC), endometrioid (ENOC), clear cell (CCOC), mucinous (MOC), low-grade serous (LGSOC)). We aimed to characterise PTEN expression as a biomarker in epithelial ovarian cancer in a large population-based study., Methods: Tumours from 5400 patients from a multicentre observational, prospective cohort study of the Ovarian Tumour Tissue Analysis Consortium were used to evaluate associations between immunohistochemical PTEN patterns and overall survival time, age, stage, grade, residual tumour, CD8+ tumour-infiltrating lymphocytes (TIL) counts, expression of oestrogen receptor (ER), progesterone receptor (PR) and androgen receptor (AR) by means of Cox proportional hazard models and generalised Cochran-Mantel-Haenszel tests., Results: Downregulation of cytoplasmic PTEN expression was most frequent in ENOC (most frequently in younger patients; p value = 0.0001) and CCOC and was associated with longer overall survival in HGSOC (hazard ratio: 0.78, 95% CI: 0.65-0.94, p value = 0.022). PTEN expression was associated with ER, PR and AR expression (p values: 0.0008, 0.062 and 0.0002, respectively) in HGSOC and with lower CD8 counts in CCOC (p value < 0.0001). Heterogeneous expression of PTEN was more prevalent in advanced HGSOC (p value = 0.019) and associated with higher CD8 counts (p value = 0.0016)., Conclusions: PTEN loss is a frequent driver in ovarian carcinoma associating distinctly with expression of hormonal receptors and CD8+ TIL counts in HGSOC and CCOC histotypes.

Venous thromboembolism (VTE) is a significant contributor to morbidity and mortality. To advance our understanding of the biology contributing to VTE, we conducted a genome-wide association study (GWAS) of VTE and a transcriptome-wide association study (TWAS) based on imputed gene expression from whole blood and liver. We meta-analyzed GWAS data from 18 studies for 30 234 VTE cases and 172 122 controls and assessed the association between 12 923 718 genetic variants and VTE. We generated variant prediction scores of gene expression from whole blood and liver tissue and assessed them for association with VTE. Mendelian randomization analyses were conducted for traits genetically associated with novel VTE loci. We identified 34 independent genetic signals for VTE risk from GWAS meta-analysis, of which 14 are newly reported associations. This included 11 newly associated genetic loci (C1orf198, PLEK, OSMR-AS1, NUGGC/SCARA5, GRK5, MPHOSPH9, ARID4A, PLCG2, SMG6, EIF5A, and STX10) of which 6 replicated, and 3 new independent signals in 3 known genes. Further, TWAS identified 5 additional genetic loci with imputed gene expression levels differing between cases and controls in whole blood (SH2B3, SPSB1, RP11-747H7.3, RP4-737E23.2) and in liver (ERAP1). At some GWAS loci, we found suggestive evidence that the VTE association signal for novel and previously known regions colocalized with expression quantitative trait locus signals. Mendelian randomization analyses suggested that blood traits may contribute to the underlying risk of VTE. To conclude, we identified 16 novel susceptibility loci for VTE; for some loci, the association signals are likely mediated through gene expression of nearby genes.

X chromosome inactivation (XCI) is a key epigenetic gene expression regulatory process, which may play a role in women's cancer. In particular tissues, some genes are known to escape XCI, yet patterns of XCI in ovarian cancer (OC) and their clinical associations are largely unknown. To examine XCI in OC, we integrated germline genotype with tumor copy number, gene expression and DNA methylation information from 99 OC patients. Approximately 10% of genes showed different XCI status (either escaping or being subject to XCI) compared with the studies of other tissues. Many of these genes are known oncogenes or tumor suppressors (e.g. DDX3X, TRAPPC2 and TCEANC). We also observed strong association between cis promoter DNA methylation and allele-specific expression imbalance (P = 2.0 × 10-10). Cluster analyses of the integrated data identified two molecular subgroups of OC patients representing those with regulated (N = 47) and dysregulated (N = 52) XCI. This XCI cluster membership was associated with expression of X inactive specific transcript (P = 0.002), a known driver of XCI, as well as age, grade, stage, tumor histology and extent of residual disease following surgical debulking. Patients with dysregulated XCI (N = 52) had shorter time to recurrence (HR = 2.34, P = 0.001) and overall survival time (HR = 1.87, P = 0.02) than those with regulated XCI, although results were attenuated after covariate adjustment. Similar findings were observed when restricted to high-grade serous tumors. We found evidence of a unique OC XCI profile, suggesting that XCI may play an important role in OC biology. Additional studies to examine somatic changes with paired tumor-normal tissue are needed., (© The Author(s) 2018. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

The United Kingdom-Primary Biliary Cholangitis (UK-PBC) risk scores are a set of prognostic models that estimate the risk of end-stage liver disease in patients with PBC at 5-, 10- and 15-year intervals. They have not been externally validated outside the United Kingdom. In this retrospective, external validation study, data were abstracted from outpatient charts and discrimination and calibration of the UK-PBC risk scores were assessed. A total of 464 patients with PBC treated with ursodeoxycholic acid were included. The median diagnosis age was 52.4 years, and 88% were female patients. The cumulative incidence of events was 6%, 9%, and 15% at 5, 10, and 15 years, respectively. Concordance (c-statistic) was 0.88, 0.85, and 0.84 using the 5-, 10- and 15-year risk scores, respectively, which was slightly lower than values observed in the United Kingdom validation cohort. Using the 5-year risk score, more events were observed than predicted (25 versus 16.8; P = 0.046); using the 10-year risk score, there was no difference between the observed and predicted number of events (35 versus 44.9; P = 0.14); conversely, using the 15-year risk score, fewer events were observed than predicted (46 versus 67.5; P = 0.009). Limiting evaluation by the 15-year UK-PBC risk score to those with >10 years of follow-up demonstrated no difference between observed and predicted events. Using the 5-year risk score, patients within the highest quartile had statistically significant worse event-free survival compared to the rest of the cohort: 82% versus 98% at 5 years, 73% versus 97% at 10 years, and 58% versus 93% at 15 years. Conclusion: In patients assessed at a North American tertiary medical center, the UK-PBC risk score had excellent discrimination and was reasonably calibrated both in the short and long term. ( Hepatology Communications 2018;2:676-682).

Objective: To evaluate myeloid differentiation primary response gene 88 (MyD88) and Toll-like receptor 4 (TLR4) expression in relation to clinical features of epithelial ovarian cancer, histologic subtypes, and overall survival., Patients and Methods: We conducted centralized immunohistochemical staining, semi-quantitative scoring, and survival analysis in 5263 patients participating in the Ovarian Tumor Tissue Analysis consortium. Patients were diagnosed between January 1, 1978, and December 31, 2014, including 2865 high-grade serous ovarian carcinomas (HGSOCs), with more than 12,000 person-years of follow-up time. Tissue microarrays were stained for MyD88 and TLR4, and staining intensity was classified using a 2-tiered system for each marker (weak vs strong)., Results: Expression of MyD88 and TLR4 was similar in all histotypes except clear cell ovarian cancer, which showed reduced expression compared with other histotypes (P<.001 for both). In HGSOC, strong MyD88 expression was modestly associated with shortened overall survival (hazard ratio [HR], 1.13; 95% CI, 1.01-1.26; P=.04) but was also associated with advanced stage (P<.001). The expression of TLR4 was not associated with survival. In low-grade serous ovarian cancer (LGSOC), strong expression of both MyD88 and TLR4 was associated with favorable survival (HR [95% CI], 0.49 [0.29-0.84] and 0.44 [0.21-0.89], respectively; P=.009 and P=.02, respectively)., Conclusion: Results are consistent with an association between strong MyD88 staining and advanced stage and poorer survival in HGSOC and demonstrate correlation between strong MyD88 and TLR4 staining and improved survival in LGSOC, highlighting the biological differences between the 2 serous histotypes., (Copyright © 2017 Mayo Foundation for Medical Education and Research. Published by Elsevier Inc. All rights reserved.)

Background: Certain diseases can occur with and without a trigger. We use Venous Thromboembolism (VTE) as our example to identify genetic interaction with pregnancy in women with VTE during pre- or postpartum. Pregnancy is one of the major risk factors for VTE as it accounts for 10% of maternal deaths., Methods: We performed a whole genome association analysis using the Cox Proportional Hazard (CoxPH) model adjusted for covariates to identify genetic variants associated with the time-to-event of VTE related to pre- or postpartum during the childbearing age of 18-45 years using a case-only design in a cohort of women with VTE. Women with a VTE event after 45 years of age were censored and contributed only follow-up time., Results: We identified two intragenic single nucleotide polymorphisms (SNPs) at genome-wide significance in the PURB gene located on chromosome 7, and two additional intragenic SNPs, one in the LINGO2 gene on chromosome 9 and one in RDXP2 on chromosome X., Conclusions: We showed that the time-to-event model is a useful approach for identifying potential hazard-modification of the genetic variants when the event of interest (VTE) occurs due to a risk factor (pre- or post-partum)., Competing Interests: The authors declare no conflict of interest.

We previously identified associations with ovarian cancer outcome at five genetic loci. To identify putatively causal genetic variants and target genes, we prioritized two ovarian outcome loci (1q22 and 19p12) for further study. Bioinformatic and functional genetic analyses indicated that MEF2D and ZNF100 are targets of candidate outcome variants at 1q22 and 19p12, respectively. At 19p12, the chromatin interaction of a putative regulatory element with the ZNF100 promoter region correlated with candidate outcome variants. At 1q22, putative regulatory elements enhanced MEF2D promoter activity and haplotypes containing candidate outcome variants modulated these effects. In a public dataset, MEF2D and ZNF100 expression were both associated with ovarian cancer progression-free or overall survival time. In an extended set of 6,162 epithelial ovarian cancer patients, we found that functional candidates at the 1q22 and 19p12 loci, as well as other regional variants, were nominally associated with patient outcome; however, no associations reached our threshold for statistical significance ( p <1×10 -5 ). Larger patient numbers will be needed to convincingly identify any true associations at these loci., Competing Interests: CONFLICTS OF INTEREST No authors have financial or other conflicts of interest.

Sequencing of RNAs (RNA-Seq) has revolutionized the field of transcriptomics, but the reads obtained often contain errors. Read error correction can have a large impact on our ability to accurately assemble transcripts. This is especially true for de novo transcriptome analysis, where a reference genome is not available. Current read error correction methods, developed for DNA sequence data, cannot handle the overlapping effects of non-uniform abundance, polymorphisms and alternative splicing. Here we present SEquencing Error CorrEction in Rna-seq data (SEECER), a hidden Markov Model (HMM)-based method, which is the first to successfully address these problems. SEECER efficiently learns hundreds of thousands of HMMs and uses these to correct sequencing errors. Using human RNA-Seq data, we show that SEECER greatly improves on previous methods in terms of quality of read alignment to the genome and assembly accuracy. To illustrate the usefulness of SEECER for de novo transcriptome studies, we generated new RNA-Seq data to study the development of the sea cucumber Parastichopus parvimensis. Our corrected assembled transcripts shed new light on two important stages in sea cucumber development. Comparison of the assembled transcripts to known transcripts in other species has also revealed novel transcripts that are unique to sea cucumber, some of which we have experimentally validated. Supporting website: http://sb.cs.cmu.edu/seecer/.