Identification of differentially expressed genes and splicing events in early-onset colorectal cancer.

Background: The incidence of colorectal cancer (CRC) has been steadily increasing in younger individuals over the past several decades for reasons that are incompletely defined. Identifying differences in gene expression profiles, or transcriptomes, in early-onset colorectal cancer (EOCRC, < 50 years old) patients versus later-onset colorectal cancer (LOCRC, > 50 years old) patients is one approach to understanding molecular and genetic features that distinguish EOCRC. Methods: We performed RNA-sequencing (RNA-seq) to characterize the transcriptomes of patient-matched tumors and adjacent, uninvolved (normal) colonic segments from EOCRC (n=21) and LOCRC (n=22) patients. The EOCRC and LOCRC cohorts were matched for demographic and clinical characteristics. We used The Cancer Genome Atlas Colon Adenocarcinoma (TCGA-COAD) database for validation. We used a series of computational and bioinformatic tools to identify EOCRC-specific differentially expressed genes, molecular pathways, predicted cell populations, differential gene splicing events, and predicted neoantigens. Results: We identified an eight-gene signature in EOCRC comprised of ALDOB, FBXL16, IL1RN, MSLN, RAC3, SLC38A11, WBSCR27 and WNT11, from which we developed a score predictive of overall CRC patient survival. On the entire set of genes identified in normal tissues and tumors, cell type deconvolution analysis predicted a differential abundance of immune and non-immune populations in EOCRC versus LOCRC. Gene set enrichment analysis identified increased expression of splicing machinery in EOCRC. We further found differences in alternative splicing (AS) events, including one within the long non-coding RNA, HOTAIRM1. Additional analysis of AS found seven events specific to EOCRC that encode potential neoantigens. Conclusion: Our transcriptome analyses identified genetic and molecular features specific to EOCRC which may inform future screening, development of prognostic indicators, and novel drug targets.

A pair of RNA binding proteins inhibit ion transporter expression to maintain lifespan.

Regulation of lifespan by transcription factors has been well established. More recently, a role for RNA binding proteins (RBPs) in regulating lifespan has also emerged. In both cases, a major challenge is to determine which regulatory targets are functionally responsible for the observed lifespan phenotype. We recently identified a pair of neuronal RBPs, exc-7/ELAVL and mbl-1/Muscleblind, which in Caenorhabditis elegans display synthetic (nonadditive) lifespan defects: single mutants do not affect lifespan, but exc-7; mbl-1 double mutants have strongly reduced lifespan. Such a strong synthetic phenotype represented an opportunity to use transcriptomics to search for potential causative targets that are synthetically regulated. Focus on such genes would allow us to narrow our target search by ignoring the hundreds of genes altered only in single mutants, and provide a shortlist of synthetically regulated candidate targets that might be responsible for the double mutant phenotype. We identified a small handful of genes synthetically dysregulated in double mutants and systematically tested each candidate gene for functional contribution to the exc-7; mbl-1 lifespan phenotype. We identified 1 such gene, the ion transporter nhx-6, which is highly upregulated in double mutants. Overexpression of nhx-6 causes reduced lifespan, and deletion of nhx-6 in an exc-7; mbl-1 background partially restores both lifespan and healthspan. Together, these results reveal that a pair of RBPs mediate lifespan in part by inhibiting expression of an ion transporter, and provide a template for how synthetic phenotypes (including lifespan) can be dissected at the transcriptomic level to reveal potential causative genes.

TDP-1 and FUST-1 co-inhibit exon inclusion and control fertility together with transcriptional regulation.

Gene expression is a multistep process and crosstalk among regulatory layers plays an important role in coordinating gene expression. To identify functionally relevant gene expression coordination, we performed a systematic reverse-genetic interaction screen in C. elegans, combining RNA binding protein (RBP) and transcription factor (TF) mutants to generate over 100 RBP;TF double mutants. We identified many unexpected double mutant phenotypes, including two strong genetic interactions between the ALS-related RBPs, fust-1 and tdp-1, and the homeodomain TF ceh-14. Losing any one of these genes alone has no effect on the health of the organism. However, fust-1;ceh-14 and tdp-1;ceh-14 double mutants both exhibit strong temperature-sensitive fertility defects. Both double mutants exhibit defects in gonad morphology, sperm function, and oocyte function. RNA-Seq analysis of double mutants identifies ceh-14 as the main controller of transcript levels, while fust-1 and tdp-1 control splicing through a shared role in exon inhibition. A skipped exon in the polyglutamine-repeat protein pqn-41 is aberrantly included in tdp-1 mutants, and genetically forcing this exon to be skipped in tdp-1;ceh-14 double mutants rescues their fertility. Together our findings identify a novel shared physiological role for fust-1 and tdp-1 in promoting C. elegans fertility and a shared molecular role in exon inhibition.

TDP-1 and FUST-1 co-inhibit exon inclusion and control fertility together with transcriptional regulation.

Gene expression is a multistep, carefully controlled process, and crosstalk between regulatory layers plays an important role in coordinating gene expression. To identify functionally relevant coordination between transcriptional and post-transcriptional gene regulation, we performed a systematic reverse-genetic interaction screen in C. elegans . We combined RNA binding protein (RBP) and transcription factor (TF) mutants, creating over 100 RBP; TF double mutants. This screen identified a variety of unexpected double mutant phenotypes, including two strong genetic interactions between the ALS-related RBPs, fust-1 and tdp-1 , and the homeodomain TF ceh-14 . Losing any one of these genes alone has no significant effect on the health of the organism. However, fust-1; ceh-14 and tdp-1; ceh-14 double mutants both exhibit strong temperature-sensitive fertility defects. Both double mutants exhibit defects in gonad morphology, sperm function, and oocyte function. RNA-seq analysis of double mutants identifies ceh-14 as the main controller of transcript levels, while fust-1 and tdp-1 control splicing through a shared role in exon inhibition. We identify a cassette exon in the polyglutamine-repeat protein pqn-41 which tdp-1 inhibits. Loss of tdp-1 causes the pqn-41 exon to be aberrantly included, and forced skipping of this exon in tdp-1; ceh-14 double mutants rescues fertility. Together our findings identify a novel shared physiological role for fust-1 and tdp-1 in promoting C. elegans fertility in a ceh-14 mutant background and reveal a shared molecular function of fust-1 and tdp-1 in exon inhibition.

A pair of RNA binding proteins inhibit ion transporter expression to maintain lifespan.

Regulation of lifespan by transcription factors has been well established. More recently a role for RNA binding proteins (RBPs) in regulating lifespan has also emerged. In both cases, a major challenge is to determine which regulatory targets are functionally responsible for the observed lifespan phenotype. We recently identified a pair of RBPs, exc-7/ELAVL and mbl-1/Muscleblind, which display synthetic (non-additive) lifespan defects: single mutants do not affect lifespan, but exc-7; mbl-1 double mutants have strongly reduced lifespan. Such a strong synthetic phenotype represented an opportunity to use transcriptomics to search for potential causative targets that are synthetically regulated. Focus on such genes would allow us to narrow our target search by ignoring the hundreds of genes altered only in single mutants, and provide a shortlist of synthetically-regulated candidate targets that might be responsible for the double mutant phenotype. We identified a small handful of genes synthetically dysregulated in double mutants and systematically tested each candidate gene for functional contribution to the exc-7; mbl-1 lifespan phenotype. We identified one such gene, the ion transporter nhx-6, which is highly upregulated in double mutants. Overexpression of nhx-6 causes reduced lifespan, and deletion of nhx-6 in an exc-7; mbl-1 background partially restores both lifespan and healthspan. Together, these results reveal that a pair of RBPs mediate lifespan in part by inhibiting expression of an ion transporter, and provide a template for how synthetic phenotypes (including lifespan) can be dissected at the transcriptomic level to reveal potential causative genes.

Molecular genetics of early-onset colorectal cancer.

Early-onset colorectal cancer (EOCRC) has been rising in global prevalence and incidence over the past several decades. Environmental influences, including generational lifestyle changes and rising obesity, contribute to these increased rates. While the rise in EOCRC is best documented in western countries, it is seen throughout the world, although EOCRC may have distinct genetic mutations in patients of different ethnic backgrounds. Pathological and molecular characterizations show that EOCRC has a distinct presentation compared with later-onset colorectal cancer (LOCRC). Recent studies have identified DNA, RNA, and protein-level alterations unique to EOCRC, revealing much-needed biomarkers and potential novel therapeutic targets. Many molecular EOCRC studies have been performed with Caucasian and Asian EOCRC cohorts, however, studies of other ethnic backgrounds are limited. In addition, certain molecular characterizations that have been conducted for LOCRC have not yet been repeated in EOCRC, including high-throughput analyses of histone modifications, mRNA splicing, and proteomics on large cohorts. We propose that the complex relationship between cancer and aging should be considered when studying the molecular underpinnings of EOCRC. In this review, we summarize current EOCRC literature, focusing on sporadic molecular alterations in tumors, and their clinical implications. We conclude by discussing current challenges and future directions of EOCRC research efforts.

TCF7L1 Regulates LGR5 Expression in Colorectal Cancer Cells.

Mutations in components of the Wnt/ß-catenin signaling pathway drive colorectal cancer (CRC), in part, by deregulating expression of genes controlled by the T-cell factor (TCF) family of transcription factors. TCFs contain a conserved DNA binding domain that mediates association with TCF binding elements (TBEs) within Wnt-responsive DNA elements (WREs). Intestinal stem cell marker, leucine-rich-repeat containing G-protein-coupled receptor 5 (LGR5), is a Wnt target gene that has been implicated in CRC stem cell plasticity. However, the WREs at the LGR5 gene locus and how TCF factors directly regulate LGR5 gene expression in CRC have not been fully defined. Here, we report that TCF family member, TCF7L1, plays a significant role in regulating LGR5 expression in CRC cells. We demonstrate that TCF7L1 binds to a novel promoter-proximal WRE through association with a consensus TBE at the LGR5 locus to repress LGR5 expression. Using CRISPR activation and interference (CRISPRa/i) technologies to direct epigenetic modulation, we demonstrate that this WRE is a critical regulator of LGR5 expression and spheroid formation capacity of CRC cells. Furthermore, we found that restoring LGR5 expression rescues the TCF7L1-mediated reduction in spheroid formation efficiency. These results demonstrate a role for TCF7L1 in repressing LGR5 gene expression to govern the spheroid formation potential of CRC cells.

Transcriptome Analyses Identify Deregulated MYC in Early Onset Colorectal Cancer.

Despite a global decrease in colorectal cancer (CRC) incidence, the prevalence of early-onset colorectal cancer (EOCRC), or those occurring in individuals before the age of 50, has steadily increased over the past several decades. When compared to later onset colorectal cancer (LOCRC) in individuals over 50, our understanding of the genetic and molecular underpinnings of EOCRCs is limited. Here, we conducted transcriptomic analyses of patient-matched normal colonic segments and tumors to identify gene expression programs involved in carcinogenesis. Amongst differentially expressed genes, we found increased expression of the c-MYC proto-oncogene (MYC) and its downstream targets in tumor samples. We identified tumors with high and low differential MYC expression and found patients with high-MYC tumors were older and overweight or obese. We also detected elevated expression of the PVT1 long-non-coding RNA (lncRNA) in most tumors and found gains in copy number for both MYC and PVT1 gene loci in 35% of tumors evaluated. Our transcriptome analyses indicate that EOCRC can be sub-classified into groups based on differential MYC expression and suggest that deregulated MYC contributes to CRCs that develop in younger patients.

Integrating 3D genomic and epigenomic data to enhance target gene discovery and drug repurposing in transcriptome-wide association studies.

Transcriptome-wide association studies (TWAS) are popular approaches to test for association between imputed gene expression levels and traits of interest. Here, we propose an integrative method PUMICE (Prediction Using Models Informed by Chromatin conformations and Epigenomics) to integrate 3D genomic and epigenomic data with expression quantitative trait loci (eQTL) to more accurately predict gene expressions. PUMICE helps define and prioritize regions that harbor cis-regulatory variants, which outperforms competing methods. We further describe an extension to our method PUMICE +, which jointly combines TWAS results from single- and multi-tissue models. Across 79 traits, PUMICE + identifies 22% more independent novel genes and increases median chi-square statistics values at known loci by 35% compared to the second-best method, as well as achieves the narrowest credible interval size. Lastly, we perform computational drug repurposing and confirm that PUMICE + outperforms other TWAS methods.

Spliceosomal component PRP-40 is a central regulator of microexon splicing.

Microexons (≤27 nt) play critical roles in nervous system development and function but create unique challenges for the splicing machinery. The mechanisms of microexon regulation are therefore of great interest. We performed a genetic screen for alternative splicing regulators in the C. elegans nervous system and identify PRP-40, a core component of the U1 snRNP. RNA-seq reveals that PRP-40 is required for inclusion of alternatively spliced, but not constitutively spliced, exons. PRP-40 is particularly required for inclusion of neuronal microexons, and our data indicate that PRP-40 is a central regulator of microexon splicing. Microexons can be relieved from PRP-40 dependence by artificially increasing exon size or reducing flanking intron size, indicating that PRP-40 is specifically required for microexons surrounded by conventionally sized introns. Knockdown of the orthologous PRPF40A in mouse neuroblastoma cells causes widespread dysregulation of microexons but not conventionally sized exons. PRP-40 regulation of neuronal microexons is therefore a widely conserved phenomenon.