Prostate cancer patient stratification by molecular signatures in the Veterans Precision Oncology Data Commons.

Veterans are at an increased risk for prostate cancer, a disease with extraordinary clinical and molecular heterogeneity, compared with the general population. However, little is known about the underlying molecular heterogeneity within the veteran population and its impact on patient management and treatment. Using clinical and targeted tumor sequencing data from the National Veterans Affairs health system, we conducted a retrospective cohort study on 45 patients with advanced prostate cancer in the Veterans Precision Oncology Data Commons (VPODC), most of whom were metastatic castration-resistant. We characterized the mutational burden in this cohort and conducted unsupervised clustering analysis to stratify patients by molecular alterations. Veterans with prostate cancer exhibited a mutational landscape broadly similar to prior studies, including KMT2A and NOTCH1 mutations associated with neuroendocrine prostate cancer phenotype, previously reported to be enriched in veterans. We also identified several potential novel mutations in PTEN, MSH6, VHL, SMO, and ABL1 Hierarchical clustering analysis revealed two subgroups containing therapeutically targetable molecular features with novel mutational signatures distinct from those reported in the Catalogue of Somatic Mutations in Cancer database. The clustering approach presented in this study can potentially be used to clinically stratify patients based on their distinct mutational profiles and identify actionable somatic mutations for precision oncology.

Immune gene variation associated with chromosome-scale differences among individual zebrafish genomes.

Immune genes have evolved to maintain exceptional diversity, offering robust defense against pathogens. We performed genomic assembly to examine immune gene variation in zebrafish. Gene pathway analysis identified immune genes as significantly enriched among genes with evidence of positive selection. A large subset of genes was absent from analysis of coding sequences due to apparent lack of reads, prompting us to examine genes overlapping zero coverage regions (ZCRs), defined as 2 kb stretches without mapped reads. Immune genes were identified as highly enriched within ZCRs, including over 60% of major histocompatibility complex (MHC) genes and NOD-like receptor (NLR) genes, mediators of direct and indirect pathogen recognition. This variation was most highly concentrated throughout one arm of chromosome 4 carrying a large cluster of NLR genes, associated with large-scale structural variation covering more than half of the chromosome. Our genomic assemblies uncovered alternative haplotypes and distinct complements of immune genes among individual zebrafish, including the MHC Class II locus on chromosome 8 and the NLR gene cluster on chromosome 4. While previous studies have shown marked variation in NLR genes between vertebrate species, our study highlights extensive variation in NLR gene regions between individuals of the same species. Taken together, these findings provide evidence of immune gene variation on a scale previously unknown in other vertebrate species and raise questions about potential impact on immune function.

Contrived Materials and a Data Set for the Evaluation of Liquid Biopsy Tests: A Blood Profiling Atlas in Cancer (BLOODPAC) Community Study.

The Blood Profiling Atlas in Cancer (BLOODPAC) Consortium is a collaborative effort involving stakeholders from the public, industry, academia, and regulatory agencies focused on developing shared best practices on liquid biopsy. This report describes the results from the JFDI (Just Freaking Do It) study, a BLOODPAC initiative to develop standards on the use of contrived materials mimicking cell-free circulating tumor DNA, to comparatively evaluate clinical laboratory testing procedures. Nine independent laboratories tested the concordance, sensitivity, and specificity of commercially available contrived materials with known variant-allele frequencies (VAFs) ranging from 0.1% to 5.0%. Each participating laboratory utilized its own proprietary evaluation procedures. The results demonstrated high levels of concordance and sensitivity at VAFs of >0.1%, but reduced concordance and sensitivity at a VAF of 0.1%; these findings were similar to those from previous studies, suggesting that commercially available contrived materials can support the evaluation of testing procedures across multiple technologies. Such materials may enable more objective comparisons of results on materials formulated in-house at each center in multicenter trials. A unique goal of the collaborative effort was to develop a data resource, the BLOODPAC Data Commons, now available to the liquid-biopsy community for further study. This resource can be used to support independent evaluations of results, data extension through data integration and new studies, and retrospective evaluation of data collection.

A highly diverse set of novel immunoglobulin-like transcript (NILT) genes in zebrafish indicates a wide range of functions with complex relationships to mammalian receptors.

Multiple novel immunoglobulin-like transcripts (NILTs) have been identified from salmon, trout, and carp. NILTs typically encode activating or inhibitory transmembrane receptors with extracellular immunoglobulin (Ig) domains. Although predicted to provide immune recognition in ray-finned fish, we currently lack a definitive framework of NILT diversity, thereby limiting our predictions for their evolutionary origin and function. In order to better understand the diversity of NILTs and their possible roles in immune function, we identified five NILT loci in the Atlantic salmon (Salmo salar) genome, defined 86 NILT Ig domains within a 3-Mbp region of zebrafish (Danio rerio) chromosome 1, and described 41 NILT Ig domains as part of an alternative haplotype for this same genomic region. We then identified transcripts encoded by 43 different NILT genes which reflect an unprecedented diversity of Ig domain sequences and combinations for a family of non-recombining receptors within a single species. Zebrafish NILTs include a sole putative activating receptor but extensive inhibitory and secreted forms as well as membrane-bound forms with no known signaling motifs. These results reveal a higher level of genetic complexity, interindividual variation, and sequence diversity for NILTs than previously described, suggesting that this gene family likely plays multiple roles in host immunity.

Host transcriptional responses in nasal swabs identify potential SARS-CoV-2 infection in PCR negative patients.

We analyzed RNA sequencing data from nasal swabs used for SARS-CoV-2 testing. 13% of 317 PCR-negative samples contained over 100 reads aligned to multiple regions of the SARS-CoV-2 genome. Differential gene expression analysis compares the host gene expression in potential false-negative (FN: PCR negative, sequencing positive) samples to subjects with multiple SARS-CoV-2 viral loads. The host transcriptional response in FN samples was distinct from true negative samples (PCR & sequencing negative) and similar to low viral load samples. Gene Ontology analysis shows viral load-dependent changes in gene expression are functionally distinct; 23 common pathways include responses to viral infections and associated immune responses. GO analysis reveals FN samples had a high overlap with high viral load samples. Deconvolution of RNA-seq data shows similar cell content across viral loads. Hence, transcriptome analysis of nasal swabs provides an additional level of identifying SARS-CoV-2 infection.

Tcf-1 promotes genomic instability and T cell transformation in response to aberrant ß-catenin activation.

Understanding the mechanisms promoting chromosomal translocations of the rearranging receptor loci in leukemia and lymphoma remains incomplete. Here we show that leukemias induced by aberrant activation of ß-catenin in thymocytes, which bear recurrent Tcra/Myc-Pvt1 translocations, depend on Tcf-1. The DNA double strand breaks (DSBs) in the Tcra site of the translocation are Rag-generated, whereas the Myc-Pvt1 DSBs are not. Aberrantly activated ß-catenin redirects Tcf-1 binding to novel DNA sites to alter chromatin accessibility and down-regulate genome-stability pathways. Impaired homologous recombination (HR) DNA repair and replication checkpoints lead to retention of DSBs that promote translocations and transformation of double-positive (DP) thymocytes. The resulting lymphomas, which resemble human T cell acute lymphoblastic leukemia (T-ALL), are sensitive to PARP inhibitors (PARPis). Our findings indicate that aberrant ß-catenin signaling contributes to translocations in thymocytes by guiding Tcf-1 to promote the generation and retention of replication-induced DSBs allowing their coexistence with Rag-generated DSBs. Thus, PARPis could offer therapeutic options in hematologic malignancies with active Wnt/ß-catenin signaling.

The Biomedical Research Hub: a federated platform for patient research data.

OBJECTIVE: The objective was to develop and operate a cloud-based federated system for managing, analyzing, and sharing patient data for research purposes, while allowing each resource sharing patient data to operate their component based upon their own governance rules. The federated system is called the Biomedical Research Hub (BRH). MATERIALS AND METHODS: The BRH is a cloud-based federated system built over a core set of software services called framework services. BRH framework services include authentication and authorization, services for generating and assessing findable, accessible, interoperable, and reusable (FAIR) data, and services for importing and exporting bulk clinical data. The BRH includes data resources providing data operated by different entities and workspaces that can access and analyze data from one or more of the data resources in the BRH. RESULTS: The BRH contains multiple data commons that in aggregate provide access to over 6 PB of research data from over 400 000 research participants. DISCUSSION AND CONCLUSION: With the growing acceptance of using public cloud computing platforms for biomedical research, and the growing use of opaque persistent digital identifiers for datasets, data objects, and other entities, there is now a foundation for systems that federate data from multiple independently operated data resources that expose FAIR application programming interfaces, each using a separate data model. Applications can be built that access data from one or more of the data resources.

Accurate and rapid prediction of tuberculosis drug resistance from genome sequence data using traditional machine learning algorithms and CNN.

Effective and timely antibiotic treatment depends on accurate and rapid in silico antimicrobial-resistant (AMR) predictions. Existing statistical rule-based Mycobacterium tuberculosis (MTB) drug resistance prediction methods using bacterial genomic sequencing data often achieve varying results: high accuracy on some antibiotics but relatively low accuracy on others. Traditional machine learning (ML) approaches have been applied to classify drug resistance for MTB and have shown more stable performance. However, there is no study that uses deep learning architecture like Convolutional Neural Network (CNN) on a large and diverse cohort of MTB samples for AMR prediction. We developed 24 binary classifiers of MTB drug resistance status across eight anti-MTB drugs and three different ML algorithms: logistic regression, random forest and 1D CNN using a training dataset of 10,575 MTB isolates collected from 16 countries across six continents, where an extended pan-genome reference was used for detecting genetic features. Our 1D CNN architecture was designed to integrate both sequential and non-sequential features. In terms of F1-scores, 1D CNN models are our best classifiers that are also more accurate and stable than the state-of-the-art rule-based tool Mykrobe predictor (81.1 to 93.8%, 93.7 to 96.2%, 93.1 to 94.8%, 95.9 to 97.2% and 97.1 to 98.2% for ethambutol, rifampicin, pyrazinamide, isoniazid and ofloxacin respectively). We applied filter-based feature selection to find AMR relevant features. All selected variant features are AMR-related ones in CARD database. 78.8% of them are also in the catalogue of MTB mutations that were recently identified as drug resistance-associated ones by WHO. To facilitate ML model development for AMR prediction, we packaged every step into an automated pipeline and shared the source code at https://github.com/KuangXY3/MTB-AMR-classification-CNN .

Immunogenomic determinants of tumor microenvironment correlate with superior survival in high-risk neuroblastoma.

BACKGROUND: Tumor-infiltrating CD8+ T cells and neoantigens are predictors of a favorable prognosis and response to immunotherapy with checkpoint inhibitors in many types of adult cancer, but little is known about their role in pediatric malignancies. Here, we analyzed the prognostic strength of T cell-inflamed gene expression and neoantigen load in high-risk neuroblastoma. We also compared transcriptional programs in T cell-inflamed and non-T cell-inflamed high-risk neuroblastomas to investigate possible mechanisms of immune exclusion. METHODS: A defined T cell-inflamed gene expression signature was used to categorize high-risk neuroblastomas in the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) program (n=123), and the Gabriella Miller Kids First (GMKF) program (n=48) into T cell-inflamed, non-T cell-inflamed, and intermediate groups. Associations between the T cell-inflamed and non-T cell-inflamed group, MYCN amplification, and survival were analyzed by Cox proportional hazards models. Additional survival analysis was conducted after integrating neoantigen load predicted from somatic mutations. Pathways activated in non-T cell-inflamed relative to T cell-inflamed tumors were analyzed using causal network analysis. RESULTS: Patients with T cell-inflamed high-risk tumors showed improved overall survival compared with those with non-T cell-inflamed tumors (p<0.05), independent of MYCN amplification status, in both TARGET and GMKF cohorts. Higher neoantigen load was also associated with better event-free and overall survival (p<0.005) and was independent of the T cell-inflamed signature. Activation of MYCN, ASCL1, SOX11, and KMT2A transcriptional programs was inversely correlated with the T cell-inflamed signature in both cohorts. CONCLUSIONS: Our results indicate that tumors from children with high-risk neuroblastoma harboring a strong T cell-inflamed signature have a more favorable clinical outcome, and neoantigen load is a prognosis predictor, independent of T cell inflammation. Strategies to target SOX11 and other signaling pathways associated with non-T cell-inflamed tumors should be pursued as potential immune-potentiating interventions.

African Americans and European Americans exhibit distinct gene expression patterns across tissues and tumors associated with immunologic functions and environmental exposures.

The COVID-19 pandemic has affected African American populations disproportionately with respect to prevalence, and mortality. Expression profiles represent snapshots of combined genetic, socio-environmental (including socioeconomic and environmental factors), and physiological effects on the molecular phenotype. As such, they have potential to improve biological understanding of differences among populations, and provide therapeutic biomarkers and environmental mitigation strategies. Here, we undertook a large-scale assessment of patterns of gene expression between African Americans and European Americans, mining RNA-Seq data from 25 non-diseased and diseased (tumor) tissue-types. We observed the widespread enrichment of pathways implicated in COVID-19 and integral to inflammation and reactive oxygen stress. Chemokine CCL3L3 expression is up-regulated in African Americans. GSTM1, encoding a glutathione S-transferase that metabolizes reactive oxygen species and xenobiotics, is upregulated. The little-studied F8A2 gene is up to 40-fold more highly expressed in African Americans; F8A2 encodes HAP40 protein, which mediates endosome movement, potentially altering the cellular response to SARS-CoV-2. African American expression signatures, superimposed on single cell-RNA reference data, reveal increased number or activity of esophageal glandular cells and lung ACE2-positive basal keratinocytes. Our findings establish basal prognostic signatures that can be used to refine approaches to minimize risk of severe infection and improve precision treatment of COVID-19 for African Americans. To enable dissection of causes of divergent molecular phenotypes, we advocate routine inclusion of metadata on genomic and socio-environmental factors for human RNA-sequencing studies.

Uniform genomic data analysis in the NCI Genomic Data Commons.

The goal of the National Cancer Institute's (NCI's) Genomic Data Commons (GDC) is to provide the cancer research community with a data repository of uniformly processed genomic and associated clinical data that enables data sharing and collaborative analysis in the support of precision medicine. The initial GDC dataset include genomic, epigenomic, proteomic, clinical and other data from the NCI TCGA and TARGET programs. Data production for the GDC started in June, 2015 using an OpenStack-based private cloud. By June of 2016, the GDC had analyzed more than 50,000 raw sequencing data inputs, as well as multiple other data types. Using the latest human genome reference build GRCh38, the GDC generated a variety of data types from aligned reads to somatic mutations, gene expression, miRNA expression, DNA methylation status, and copy number variation. In this paper, we describe the pipelines and workflows used to process and harmonize the data in the GDC. The generated data, as well as the original input files from TCGA and TARGET, are available for download and exploratory analysis at the GDC Data Portal and Legacy Archive ( https://gdc.cancer.gov/ ).

SARS-CoV-2 genomic diversity and the implications for qRT-PCR diagnostics and transmission.

The COVID-19 pandemic has sparked an urgent need to uncover the underlying biology of this devastating disease. Though RNA viruses mutate more rapidly than DNA viruses, there are a relatively small number of single nucleotide polymorphisms (SNPs) that differentiate the main SARS-CoV-2 lineages that have spread throughout the world. In this study, we investigated 129 RNA-seq data sets and 6928 consensus genomes to contrast the intra-host and inter-host diversity of SARS-CoV-2. Our analyses yielded three major observations. First, the mutational profile of SARS-CoV-2 highlights intra-host single nucleotide variant (iSNV) and SNP similarity, albeit with differences in C > U changes. Second, iSNV and SNP patterns in SARS-CoV-2 are more similar to MERS-CoV than SARS-CoV-1. Third, a significant fraction of insertions and deletions contribute to the genetic diversity of SARS-CoV-2. Altogether, our findings provide insight into SARS-CoV-2 genomic diversity, inform the design of detection tests, and highlight the potential of iSNVs for tracking the transmission of SARS-CoV-2.

Hidden genomic diversity of SARS-CoV-2: implications for qRT-PCR diagnostics and transmission.

The COVID-19 pandemic has sparked an urgent need to uncover the underlying biology of this devastating disease. Though RNA viruses mutate more rapidly than DNA viruses, there are a relatively small number of single nucleotide polymorphisms (SNPs) that differentiate the main SARS-CoV-2 clades that have spread throughout the world. In this study, we investigated over 7,000 SARS-CoV-2 datasets to unveil both intrahost and interhost diversity. Our intrahost and interhost diversity analyses yielded three major observations. First, the mutational profile of SARS-CoV-2 highlights iSNV and SNP similarity, albeit with high variability in C>T changes. Second, iSNV and SNP patterns in SARS-CoV-2 are more similar to MERS-CoV than SARS-CoV-1. Third, a significant fraction of small indels fuel the genetic diversity of SARS-CoV-2. Altogether, our findings provide insight into SARS-CoV-2 genomic diversity, inform the design of detection tests, and highlight the potential of iSNVs for tracking the transmission of SARS-CoV-2.

ACE2 and TMPRSS2 expression by clinical, HLA, immune, and microbial correlates across 34 human cancers and matched normal tissues: implications for SARS-CoV-2 COVID-19.

BACKGROUND: Pandemic COVID-19 by severe acute respiratory syndrome (SARS) coronavirus 2 (SARS-CoV-2) infection is facilitated by the ACE2 receptor and protease TMPRSS2. Modestly sized case series have described clinical factors associated with COVID-19, while ACE2 and TMPRSS2 expression analyses have been described in some cell types. Patients with cancer may have worse outcomes to COVID-19. METHODS: We performed an integrated study of ACE2 and TMPRSS2 gene expression across and within organ systems, by normal versus tumor, across several existing databases (The Cancer Genome Atlas, Census of Immune Single Cell Expression Atlas, The Human Cell Landscape, and more). We correlated gene expression with clinical factors (including but not limited to age, gender, race, body mass index, and smoking history), HLA genotype, immune gene expression patterns, cell subsets, and single-cell sequencing as well as commensal microbiome. RESULTS: Matched normal tissues generally display higher ACE2 and TMPRSS2 expression compared with cancer, with normal and tumor from digestive organs expressing the highest levels. No clinical factors were consistently identified to be significantly associated with gene expression levels though outlier organ systems were observed for some factors. Similarly, no HLA genotypes were consistently associated with gene expression levels. Strong correlations were observed between ACE2 expression levels and multiple immune gene signatures including interferon-stimulated genes and the T cell-inflamed phenotype as well as inverse associations with angiogenesis and transforming growth factor-ß signatures. ACE2 positively correlated with macrophage subsets across tumor types. TMPRSS2 was less associated with immune gene expression but was strongly associated with epithelial cell abundance. Single-cell sequencing analysis across nine independent studies demonstrated little to no ACE2 or TMPRSS2 expression in lymphocytes or macrophages. ACE2 and TMPRSS2 gene expression associated with commensal microbiota in matched normal tissues particularly from colorectal cancers, with distinct bacterial populations showing strong associations. CONCLUSIONS: We performed a large-scale integration of ACE2 and TMPRSS2 gene expression across clinical, genetic, and microbiome domains. We identify novel associations with the microbiota and confirm host immunity associations with gene expression. We suggest caution in interpretation regarding genetic associations with ACE2 expression suggested from smaller case series.

Before and After: Comparison of Legacy and Harmonized TCGA Genomic Data Commons' Data.

We present a systematic analysis of the effects of synchronizing a large-scale, deeply characterized, multi-omic dataset to the current human reference genome, using updated software, pipelines, and annotations. For each of 5 molecular data platforms in The Cancer Genome Atlas (TCGA)-mRNA and miRNA expression, single nucleotide variants, DNA methylation and copy number alterations-comprehensive sample, gene, and probe-level studies were performed, towards quantifying the degree of similarity between the 'legacy' GRCh37 (hg19) TCGA data and its GRCh38 (hg38) version as 'harmonized' by the Genomic Data Commons. We offer gene lists to elucidate differences that remained after controlling for confounders, and strategies to mitigate their impact on biological interpretation. Our results demonstrate that the hg19 and hg38 TCGA datasets are very highly concordant, promote informed use of either legacy or harmonized omics data, and provide a rubric that encourages similar comparisons as new data emerge and reference data evolve.

Glucocorticoid receptor modulation decreases ER-positive breast cancer cell proliferation and suppresses wild-type and mutant ER chromatin association.

BACKGROUND: Non-ER nuclear receptor activity can alter estrogen receptor (ER) chromatin association and resultant ER-mediated transcription. Consistent with GR modulation of ER activity, high tumor glucocorticoid receptor (GR) expression correlates with improved relapse-free survival in ER+ breast cancer (BC) patients. METHODS: In vitro cell proliferation assays were used to assess ER-mediated BC cell proliferation following GR modulation. ER chromatin association following ER/GR co-liganding was measured using global ChIP sequencing and directed ChIP analysis of proliferative gene enhancers. RESULTS: We found that GR liganding with either a pure agonist or a selective GR modulator (SGRM) slowed estradiol (E2)-mediated proliferation in ER+ BC models. SGRMs that antagonized transcription of GR-unique genes both promoted GR chromatin association and inhibited ER chromatin localization at common DNA enhancer sites. Gene expression analysis revealed that ER and GR co-activation decreased proliferative gene activation (compared to ER activation alone), specifically reducing CCND1, CDK2, and CDK6 gene expression. We also found that ligand-dependent GR occupancy of common ER-bound enhancer regions suppressed both wild-type and mutant ER chromatin association and decreased corresponding gene expression. In vivo, treatment with structurally diverse SGRMs also reduced MCF-7 Y537S ER-expressing BC xenograft growth. CONCLUSION: These studies demonstrate that liganded GR can suppress ER chromatin occupancy at shared ER-regulated enhancers, including CCND1 (Cyclin D1), regardless of whether the ligand is a classic GR agonist or antagonist. Resulting GR-mediated suppression of ER+ BC proliferative gene expression and cell division suggests that SGRMs could decrease ER-driven gene expression.

Proteomics reveals NNMT as a master metabolic regulator of cancer-associated fibroblasts.

High-grade serous carcinoma has a poor prognosis, owing primarily to its early dissemination throughout the abdominal cavity. Genomic and proteomic approaches have provided snapshots of the proteogenomics of ovarian cancer1,2, but a systematic examination of both the tumour and stromal compartments is critical in understanding ovarian cancer metastasis. Here we develop a label-free proteomic workflow to analyse as few as 5,000 formalin-fixed, paraffin-embedded cells microdissected from each compartment. The tumour proteome was stable during progression from in situ lesions to metastatic disease; however, the metastasis-associated stroma was characterized by a highly conserved proteomic signature, prominently including the methyltransferase nicotinamide N-methyltransferase (NNMT) and several of the proteins that it regulates. Stromal NNMT expression was necessary and sufficient for functional aspects of the cancer-associated fibroblast (CAF) phenotype, including the expression of CAF markers and the secretion of cytokines and oncogenic extracellular matrix. Stromal NNMT expression supported ovarian cancer migration, proliferation and in vivo growth and metastasis. Expression of NNMT in CAFs led to depletion of S-adenosyl methionine and reduction in histone methylation associated with widespread gene expression changes in the tumour stroma. This work supports the use of ultra-low-input proteomics to identify candidate drivers of disease phenotypes. NNMT is a central, metabolic regulator of CAF differentiation and cancer progression in the stroma that may be therapeutically targeted.

IBD-associated Colon Cancers Differ in DNA Methylation and Gene Expression Profiles Compared With Sporadic Colon Cancers.

BACKGROUND AND AIMS: As ulcerative colitis [UC]-associated colorectal cancer [CRC] and sporadic CRC differ in presentation and molecular features, we sought to evaluate differences in the impact of DNA methylation on gene expression. METHODS: DNA methylation was assessed in 11 UC-CRCs and adjacent tissue and 11 sporadic CRCs and adjacent tissue, using Illumina arrays. RNA sequencing was performed on 10 UC-CRCs and adjacent tissue and eight sporadic CRCs and adjacent tissues. Differences in DNA methylation and transcript expression, as well as their correlation in the same tissues, were assessed. Immunohistochemistry was performed for three proteins, ANPEP, FAM92A1, and STK31, all of which exhibited an inverse correlation between DNA methylation and transcript expression in UC. RESULTS: Thirty three loci demonstrated differences in DNA methylation between UC-CRC and adjacent tissue. In contrast, there were 4204 differentially methylated loci between sporadic colon cancer and adjacent tissue. Eight hundred eighty six genes as well as 10 long non-coding RNAs [lncRNA] were differentially expressed between UC-CRC and adjacent tissues. Although there were no differentially methylated loci between UC and sporadic CRC, 997 genes and 38 lncRNAs were differentially expressed between UC-CRC and sporadic CRC. In UC, 18 genes demonstrated a negative correlation between DNA methylation and transcript expression. Evaluation of protein expression related to three genes, ANPEP, FAM92A1, and STK31, confirmed down-regulation of ANPEP and up-regulation of STK31 in UC-CRC. CONCLUSIONS: Regulation of transcript expression by DNA methylation involves genes key to colon carcinogenesis and may account for differences in presentation and outcomes between inflammatory bowel disease and sporadic colon cancer.