Colon Crypts of Subjects With Familial Adenomatous Polyposis Show an Increased Number of LGR5+ Ectopic Stem Cells

INTRODUCTION

Familial adenomatous polyposis (FAP) is a hereditary colorectal cancer (CRC) syndrome characterized by accelerated adenoma development due to inherited (or de novo) mutations in the APC regulator of WNT signaling pathway (APC) gene. The mechanism underlying this accelerated polyp development in subjects with FAP has not been defined. Given that LGR5+ stem cells drive crypt cell proliferation, we hypothesized that FAP crypts would demonstrate aberrant leucine-rich repeat-containing G-protein-coupled receptor 5 (LGR5) staining patterns.

METHODS

Biopsies were taken from 11 healthy subjects, 7 subjects with Lynch syndrome, 4 subjects with FAP, and 1 subject with MUTYH-associated polyposis syndrome during routine screening or surveillance colonoscopy. Crypt staining was evaluated by immunohistochemistry of paraffin-embedded tissue sections. Stem cell numbers were estimated by immunofluorescence staining of isolated crypts using antibodies against LGR5 and other proteins.

RESULTS

Subjects with FAP exhibited a greater number of LGR5+ stem cells in their crypts than healthy subjects and subjects with Lynch syndrome and MUTYH-associated polyposis syndrome. Most crypts of subjects with FAP harbored LGR5+ cells located above the lower third of the crypts.

DISCUSSION

These findings support a model in which inactivation of one copy of APC leads to increased numbers of LGR5+ stem cells, many of which are ectopic, in colon crypts of subjects with FAP. Overabundant and ectopic LGR5+ stem cells could lead to an expanded proliferative zone of dividing cells more likely to develop mutations that would contribute to the accelerated adenoma development observed in FAP.

SNPs associated with colorectal cancer at 15q13.3 affect risk enhancers that modulate GREM1 gene expression

Several genome wide association studies of colorectal cancer (CRC) have identified single nucleotide polymorphisms (SNPs) on chromosome 15q13.3 associated with CRC risk. To identify functional variant(s) underlying this association, we investigated SNPs in linkage disequilibrium with the risk‐associated SNP rs4779584 that overlapped regulatory regions/enhancer elements characterized in colon‐related tissues and cells. We identified several SNP‐containing regulatory regions that exhibited enhancer activity in vitro, including one SNP (rs1406389) that correlated with allele‐specific effects on enhancer activity. Deletion of either this enhancer or another enhancer that had previously been reported in this region correlated with decreased expression of GREM1 following CRISPR/Cas9 genome editing. That GREM1 is one target of these enhancers was further supported by an expression quantitative trait loci correlation between rs1406389 and GREM1 expression in the transverse but not sigmoid colon in the Genotype‐Tissue Expression dataset. Taken together, we conclude that the 15q13.3 region contains at least two functional variants that map to distinct enhancers and impact CRC risk through modulation of GREM1 expression.

Novel Common Genetic Susceptibility Loci for Colorectal Cancer

BACKGROUND

Previous genome-wide association studies (GWAS) have identified 42 loci (P < 5 × 10-8) associated with risk of colorectal cancer (CRC). Expanded consortium efforts facilitating the discovery of additional susceptibility loci may capture unexplained familial risk.

METHODS

We conducted a GWAS in European descent CRC cases and control subjects using a discovery-replication design, followed by examination of novel findings in a multiethnic sample (cumulative n = 163 315). In the discovery stage (36 948 case subjects/30 864 control subjects), we identified genetic variants with a minor allele frequency of 1% or greater associated with risk of CRC using logistic regression followed by a fixed-effects inverse variance weighted meta-analysis. All novel independent variants reaching genome-wide statistical significance (two-sided P < 5 × 10-8) were tested for replication in separate European ancestry samples (12 952 case subjects/48 383 control subjects). Next, we examined the generalizability of discovered variants in East Asians, African Americans, and Hispanics (12 085 case subjects/22 083 control subjects). Finally, we examined the contributions of novel risk variants to familial relative risk and examined the prediction capabilities of a polygenic risk score. All statistical tests were two-sided.

RESULTS

The discovery GWAS identified 11 variants associated with CRC at P < 5 × 10-8, of which nine (at 4q22.2/5p15.33/5p13.1/6p21.31/6p12.1/10q11.23/12q24.21/16q24.1/20q13.13) independently replicated at a P value of less than .05. Multiethnic follow-up supported the generalizability of discovery findings. These results demonstrated a 14.7% increase in familial relative risk explained by common risk alleles from 10.3% (95% confidence interval [CI] = 7.9% to 13.7%; known variants) to 11.9% (95% CI = 9.2% to 15.5%; known and novel variants). A polygenic risk score identified 4.3% of the population at an odds ratio for developing CRC of at least 2.0.

CONCLUSIONS

This study provides insight into the architecture of common genetic variation contributing to CRC etiology and improves risk prediction for individualized screening.

Multiple functional risk variants in a SMAD7 enhancer implicate a colorectal cancer risk haplotype

Genome-wide association studies (GWAS) of colorectal cancer (CRC) have led to the identification of a number of common variants associated with modest risk. Several risk variants map within the vicinity of TGFβ/BMP signaling pathway genes, including rs4939827 within an intron of SMAD7 at 18q21.1. A previous study implicated a novel SNP (novel 1 or rs58920878) as a functional variant within an enhancer element in SMAD7 intron 4. In this study, we show that four SNPs including novel 1 (rs6507874, rs6507875, rs8085824, and rs58920878) in linkage disequilibrium (LD) with the index SNP rs4939827 demonstrate allele-specific enhancer effects in a large, multi-component enhancer of SMAD7. All four SNPs demonstrate allele-specific protein binding to nuclear extracts of CRC cell lines. Furthermore, some of the risk-associated alleles correlate with increased expression of SMAD7 in normal colon tissues. Finally, we show that the enhancer is responsive to BMP4 stimulation. Taken together, we propose that the associated CRC risk at 18q21.1 is due to four functional variants that regulate SMAD7 expression and potentially perturb a BMP negative feedback loop in TGFβ/BMP signaling pathways.

Identification and characterization of functional risk variants for colorectal cancer mapping to chromosome 11q23.1

Genome-wide association studies of colorectal cancer (CRC) have identified a number of common variants associated with modest risk, including rs3802842 at chromosome 11q23.1. Several genes map to this region but rs3802842 does not map to any known transcribed or regulatory sequences. We reasoned, therefore, that rs3802842 is not the functional single-nucleotide polymorphism (SNP), but is in linkage disequilibrium (LD) with a functional SNP(s). We performed ChIP-seq for histone modifications in SW480 and HCT-116 CRC cells, and incorporated ChIP-seq and DNase I hypersensitivity data available through ENCODE within a 137-kb genomic region containing rs3802842 on 11q23.1. We identified SNP rs10891246 in LD with rs3802842 that mapped within a bidirectional promoter region of genes C11orf92 and C11orf93. Following mutagenesis to the risk allele, the promoter demonstrated lower levels of reporter gene expression. A second SNP rs7130173 was identified in LD with rs3802842 that mapped to a candidate enhancer region, which showed strong unidirectional activity in both HCT-116 and SW480 CRC cells. The risk allele of rs7130173 demonstrated reduced enhancer activity compared with the common allele, and reduced nuclear protein binding affinity in electromobility shift assays compared with the common allele suggesting differential transcription factor (TF) binding. SNPs rs10891246 and rs7130173 are on the same haplotype, and expression quantitative trait loci (eQTL) analyses of neighboring genes implicate C11orf53, C11orf92 and C11orf93 as candidate target genes. These data imply that rs10891246 and rs7130173 are functional SNPs mapping to 11q23.1 and that C11orf53, C11orf92 and C11orf93 represent novel candidate target genes involved in CRC etiology.

RNA splicing is responsive to MBNL1 dose

Myotonic dystrophy (DM1) is a highly variable, multi-system disorder resulting from the expansion of an untranslated CTG tract in DMPK. In DM1 expanded CUG repeat RNAs form hairpin secondary structures that bind and aberrantly sequester the RNA splice regulator, MBNL1. RNA splice defects resulting as a consequence of MBNL1 depletion have been shown to play a key role in the development of DM1 pathology. In patient populations, both the number and severity of DM1 symptoms increase broadly as a function of CTG tract length. However significant variability in the DM1 phenotype is observed in patients encoding similar CTG repeat numbers. Here we demonstrate that a gradual decrease in MBNL1 levels results both in the expansion of the repertoire of splice defects and an increase in the severity of the splice alterations. Thus, MBNL1 loss does not have an all or none outcome but rather shows a graded effect on the number and severity of the ensuing splice defects. Our results suggest that once a critical threshold is reached, relatively small dose variations of free MBNL1 levels, which may reflect modest changes in the size of the CUG tract or the extent of hairpin secondary structure formation, can significantly alter the number and severity of splice abnormalities and thus contribute to the phenotype variability observed in DM1 patients.

RNA steady-state defects in myotonic dystrophy are linked to nuclear exclusion of SHARP

We describe a new mechanism by which CTG tract expansion affects myotonic dystrophy (DM1). Changes to the levels of a panel of RNAs involved in muscle development and function that are downregulated in DM1 are due to aberrant localization of the transcription factor SHARP (SMART/HDAC1-associated repressor protein). Mislocalization of SHARP in DM1 is consistent with increased CRM1-mediated export of SHARP to the cytoplasm. A direct link between CTG repeat expression and SHARP mislocalization is demonstrated as expression of expanded CTG repeats in normal cells recapitulates cytoplasmic SHARP localization. These results demonstrate a role for the inactivation of SHARP transcription in DM1 biology.