Identification of Novel Candidate Genes for Early-Onset Colorectal Cancer Susceptibility

Approximately 25–30% of colorectal cancer (CRC) cases are expected to result from a genetic predisposition, but in only 5–10% of these cases highly penetrant germline mutations are found. The remaining CRC heritability is still unexplained, and may be caused by a hitherto-undefined set of rare variants with a moderately penetrant risk. Here we aimed to identify novel risk factors for early-onset CRC using whole-exome sequencing, which was performed on a cohort of CRC individuals (n = 55) with a disease onset before 45 years of age. We searched for genes that were recurrently affected by rare variants (minor allele frequency ≤0.001) with potentially damaging effects and, subsequently, re-sequenced the candidate genes in a replication cohort of 174 early-onset or familial CRC individuals. Two functionally relevant genes with low frequency variants with potentially damaging effects, PTPN12 and LRP6, were found in at least three individuals. The protein tyrosine phosphatase PTP-PEST, encoded by PTPN12, is a regulator of cell motility and LRP6 is a component of the WNT-FZD-LRP5-LRP6 complex that triggers WNT signaling. All variants in LRP6 were identified in individuals with an extremely early-onset of the disease (≤30 years of age), and two of the three variants showed increased WNT signaling activity in vitro. In conclusion, we present PTPN12 and LRP6 as novel candidates contributing to the heterogeneous susceptibility to CRC.     

km23‐1/DYNLRB1 regulation of MEK/ERK signaling and R‐Ras in invasive human colorectal cancer cells

We previously found that km23‐1/DYNLRB1 is required for transforming growth factor‐β (TGFβ) production through Ras/ERK pathways in TGFβ‐sensitive epithelial cells and in human colorectal cancer (CRC) cells. Here we demonstrate that km23‐1/DYNLRB1 is required for mitogen‐activated protein kinase kinase (MEK) activation in human CRC cells, detected by km23‐1/DYNLRB1‐siRNA inhibition of phospho‐(p)‐MEK immunostaining in RKO cells. Furthermore, we show that CRISPR‐Cas9 knock‐out (KO) of km23‐1/DYNLRB1 reduced cell migration in two additional CRC models, HCT116 and DLD‐1. Of interest, in contrast to our previous work showing that dynein motor activity was required for TGFβ‐mediated nuclear translocation of Smad2, in the current report, we demonstrate for the first time that disruption of dynein motor activity did not reduce TGFβ‐mediated activation of MEK1/2 or c‐Jun N‐terminal kinase (JNK). Moreover, size exclusion chromatography of RKO cell lysates revealed that B‐Raf, extracellular signal‐regulated kinase (ERK), and p‐ERK were not present in the large molecular weight fractions containing dynein holocomplex components. Furthermore, sucrose gradient fractionation of cell lysates from both HCT116 and CBS CRC cells demonstrated that km23‐1/DYNLRB1 co‐sedimented with Ras, p‐ERK, and ERK in fractions that did not contain components of holo‐dynein. Thus, km23‐1/DYNLRB1 may be associated with activated Ras/ERK signaling complexes in cell compartments that do not contain the dynein holoprotein complex, suggesting dynein‐independent km23‐1/DYNLRB1 functions in Ras/ERK signaling. Finally, of the Ras isoforms, R‐Ras is most often associated with cell migration, adhesion, and protrusive activity. Here, we show that a significant fraction of km23‐1/DYNLRB1 and RRas wase co‐localized at the protruding edges of migrating HCT116 cells, suggesting an important role for the km23‐1/DYNLRB1‐R‐Ras complex in CRC invasion.     

Pathway crosstalk between Ras/Raf and PI3K in promotion of M-CSF-induced MEK/ERK-mediated osteoclast survival

While M-CSF-mediated MEK/ERK activation promotes osteoclast survival, the signaling pathway by which M-CSF activates MEK/ERK is unresolved. Functions for PI3K, Ras, and Raf have been implicated in support of osteoclast survival, although interaction between these signaling components has not been examined. Therefore, the interplay between PI3K, Ras and Raf in M-CSF-promoted MEK/ERK activation and osteoclast survival was investigated. M-CSF activates Ras to coordinate activation of PI3K and Raf/MEK/ERK, since Ras inhibition decreased PI3K activation and PI3K inhibition did not block M-CSF-mediated Ras activation. As further support for Ras-mediated signaling, constitutively active (ca) Ras promoted MEK/ERK activation and osteoclast survival, which was blocked by inhibition of PI3K or Raf. Moreover, PI3K-selective or Raf-selective caRas were only partially able to promote osteoclast survival when compared to parental caRas. We then examined whether PI3K and Raf function linearly or in parallel downstream of Ras. Expression of caPI3K increased MEK/ERK activation and promoted osteoclast survival downstream of M-CSF, supporting this hypothesis. Blocking Raf did not decrease osteoclast survival and MEK/ERK activation promoted by caPI3K. In addition, PI3K-selective Ras-mediated survival was not blocked by Raf inhibition. Taken together, our data support that Raf signaling is separate from Ras/PI3K signaling and PI3K signaling is separate from Ras/Raf signaling. These data therefore support a role for Ras in coordinate activation of PI3K and Raf acting in parallel to mediate MEK/ERK-promoted osteoclast survival induced by M-CSF.     

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.     

Small G Proteins Rac1 and Ras Regulate Serine/Threonine Protein Phosphatase 5 (PP5)·Extracellular Signal-Regulated Kinase (ERK) Complexes Involved in the Feedback Regulation of Raf1

Serine/threonine protein phosphatase 5 (PP5, PPP5C) is known to interact with the chaperonin heat shock protein 90 (HSP90) and is involved in the regulation of multiple cellular signaling cascades that control diverse cellular processes, such as cell growth, differentiation, proliferation, motility, and apoptosis. Here, we identify PP5 in stable complexes with extracellular signal-regulated kinases (ERKs). Studies using mutant proteins reveal that the formation of PP5·ERK1 and PP5·ERK2 complexes partially depends on HSP90 binding to PP5 but does not require PP5 or ERK1/2 activity. However, PP5 and ERK activity regulates the phosphorylation state of Raf1 kinase, an upstream activator of ERK signaling. Whereas expression of constitutively active Rac1 promotes the assembly of PP5·ERK1/2 complexes, acute activation of ERK1/2 fails to influence the phosphatase-kinase interaction. Introduction of oncogenic HRas (HRas^V12) has no effect on PP5-ERK1 binding but selectively decreases the interaction of PP5 with ERK2, in a manner that is independent of PP5 and MAPK/ERK kinase (MEK) activity, yet paradoxically requires ERK2 activity. Additional studies conducted with oncogenic variants of KRas4B reveal that KRas^L61, but not KRas^V12, also decreases the PP5-ERK2 interaction. The expression of wild type HRas or KRas proteins fails to reduce PP5-ERK2 binding, indicating that the effect is specific to HRas^V12 and KRas^L61 gain-of-function mutations. These findings reveal a novel, differential responsiveness of PP5-ERK1 and PP5-ERK2 interactions to select oncogenic Ras variants and also support a role for PP5·ERK complexes in regulating the feedback phosphorylation of PP5-associated Raf1.     

A Combination of CD28 (rs1980422) and IRF5 (rs10488631) Polymorphisms Is Associated with Seropositivity in Rheumatoid Arthritis: A Case Control Study

IntroductionThe aim of the study was to analyse genetic architecture of RA by utilizing multiparametric statistical methods such as linear discriminant analysis (LDA) and redundancy analysis (RDA).MethodsA total of 1393 volunteers, 499 patients with RA and 894 healthy controls were included in the study. The presence of shared epitope (SE) in HLA-DRB1 and 11 SNPs (PTPN22 C/T (rs2476601), STAT4 G/T (rs7574865), CTLA4 A/G (rs3087243), TRAF1/C5 A/G (rs3761847), IRF5 T/C (rs10488631), TNFAIP3 C/T (rs5029937), AFF3 A/T (rs11676922), PADI4 C/T (rs2240340), CD28 T/C (rs1980422), CSK G/A (rs34933034) and FCGR3A A/C (rs396991), rheumatoid factor (RF), anti–citrullinated protein antibodies (ACPA) and clinical status was analysed using the LDA and RDA.ResultsHLA-DRB1, PTPN22, STAT4, IRF5 and PADI4 significantly discriminated between RA patients and healthy controls in LDA. The correlation between RA diagnosis and the explanatory variables in the model was 0.328 (Trace = 0.107; F = 13.715; P = 0.0002). The risk variants of IRF5 and CD28 genes were found to be common determinants for seropositivity in RDA, while positivity of RF alone was associated with the CTLA4 risk variant in heterozygous form. The correlation between serologic status and genetic determinants on the 1st ordinal axis was 0.468, and 0.145 on the 2nd one (Trace = 0.179; F = 6.135; P = 0.001). The risk alleles in AFF3 gene together with the presence of ACPA were associated with higher clinical severity of RA.ConclusionsThe association among multiple risk variants related to T cell receptor signalling with seropositivity may play an important role in distinct clinical phenotypes of RA. Our study demonstrates that multiparametric analyses represent a powerful tool for investigation of mutual relationships of potential risk factors in complex diseases such as RA.     

Mitogen/extracellular signal-regulated kinase kinase-5 promoter region polymorphisms affect the risk of sporadic colorectal cancer in a southern Chinese population

Mitogen/extracellular signal-regulated kinase kinase-5 (MEK5), which belongs to a network of mitogen-activated protein kinase pathways, play a pivotal role in carcinogenesis. The purpose of this study was to investigate whether variants in the MEK5 gene promoter were involved in susceptivity of individuals to sporadic colorectal cancer (CRC). In the present hospital-based case-control study of 737 patients with sporadic CRC and 703 healthy control subjects in a southern Chinese population, the two polymorphisms of MEK5 promoter (i.e., rs7172582C>T and rs3743354T>C) were genotyped by TaqMan assay. There were significant differences between cases and controls in the genotype and allele distribution of the MEK5 gene rs3743354T>C polymorphism. The rs3743354 CC genotype was associated with a significantly decreased risk of CRC when compared with the TT genotype (adjusted odds ratios [ORs]=0.43; 95% confidence interval [CI], 0.24-0.77). Compared to the T allele, a significant correlation was detected between the presence of the C allele and decreased risk of CRC (adjusted OR=0.79; 95% CI, 0.61-0.94). The decreased risk of CRC associated with rs3743354 variant genotypes (i.e., CT+CC) was found in the smoker subgroup (adjusted OR=0.63; 95% CI=0.45-0.88). Further, environmental factors, including smoking and drinking, interacted with rs3743354C variant genotypes to reduce CRC risk. Western blot analysis showed that the levels of MEK5 protein in sporadic CRC neoplastic tissues and adjacent normal colorectal epithelium tissues were lower in the carriers of rs3743354 CC genotypes than that in those with rs3743354 TT genotypes or those with rs3743354 TC genotypes. However, no significant association was found between the rs7172582C>T polymorphism and risk of CRC. These data indicate that the rs3743354 polymorphism in the MEK5 promoter may affect the risk of developing CRC.     

A comprehensive investigation on common polymorphisms in the MDR1/ABCB1 transporter gene and susceptibility to colorectal cancer

ATP Binding Cassette B1 (ABCB1) is a transporter with a broad substrate specificity involved in the elimination of several carcinogens from the gut. Several polymorphic variants within the ABCB1 gene have been reported as modulators of ABCB1-mediated transport. We investigated the impact of ABCB1 genetic variants on colorectal cancer (CRC) risk. A hybrid tagging/functional approach was performed to select 28 single nucleotide polymorphisms (SNPs) that were genotyped in 1,321 Czech subjects, 699 CRC cases and 622 controls. In addition, six potentially functional SNPs were genotyped in 3,662 German subjects, 1,809 cases and 1,853 controls from the DACHS study. We found that three functional SNPs (rs1202168, rs1045642 and rs868755) were associated with CRC risk in the German population. Carriers of the rs1202168_T and rs868755_T alleles had an increased risk for CRC (P_trend = 0.016 and 0.029, respectively), while individuals bearing the rs1045642_C allele showed a decreased risk of CRC (P_trend = 0.022). We sought to replicate the most significant results in an independent case-control study of 3,803 subjects, 2,169 cases and 1,634 controls carried out in the North of Germany. None of the SNPs tested were significantly associated with CRC risk in the replication study. In conclusion, in this study of about 8,800 individuals we show that ABCB1 gene polymorphisms play at best a minor role in the susceptibility to CRC.     

Vitamin D receptor polymorphism and colorectal cancer-specific and all-cause mortality

BackgroundThe vitamin D receptor (VDR) gene is present in colorectal cancer (CRC) cells and its genetic variants have been associated with an increased risk of CRC. The association with colorectal cancer prognosis remains widely unexplored.Methods1397 colorectal cancer patients participating in two cancer cohorts (ESTHER II and VERDI) and in a population-based case–control study (DACHS) were followed for 5 years. Unadjusted and adjusted hazard ratios for all-cause mortality (469 events) and CRC-specific mortality (336 events) were estimated for VDR variants rs731236 (TaqI), rs2228570 (FokI), rs11568820 (Cdx2), and rs1989969 (VDR-5132).ResultsNo association was found between VDR polymorphism and CRC specific and all-cause mortality. Adjusted hazard ratios ranged from 0.79 (95% CI 0.57–1.12) to 1.14 (95% CI 0.89–1.46) for CRC-specific mortality and from 0.89 (95% CI 0.67–1.18) to 1.22 (95% CI 0.99–1.50) for all-cause mortality. All 95% confidence intervals included the null value.ConclusionsOur findings do not support the hypothesis that the common VDR gene variants investigated in this study are of clinical relevance with respect to CRC prognosis.     

Phosphorylation of Serine 779 in Fibroblast Growth Factor Receptor 1 and 2 by Protein Kinase C? Regulates Ras/Mitogen-activated Protein Kinase Signaling and Neuronal Differentiation

The FGF receptors (FGFRs) control a multitude of cellular processes both during development and in the adult through the initiation of signaling cascades that regulate proliferation, survival, and differentiation. Although FGFR tyrosine phosphorylation and the recruitment of Src homology 2 domain proteins have been widely described, we have previously shown that FGFR is also phosphorylated on Ser⁷⁷⁹ in response to ligand and binds the 14-3-3 family of phosphoserine/threonine-binding adaptor/scaffold proteins. However, whether this receptor phosphoserine mode of signaling is able to regulate specific signaling pathways and biological responses is unclear. Using PC12 pheochromocytoma cells and primary mouse bone marrow stromal cells as models for growth factor-regulated neuronal differentiation, we show that Ser⁷⁷⁹ in the cytoplasmic domains of FGFR1 and FGFR2 is required for the sustained activation of Ras and ERK but not for other FGFR phosphotyrosine pathways. The regulation of Ras and ERK signaling by Ser⁷⁷⁹ was critical not only for neuronal differentiation but also for cell survival under limiting growth factor concentrations. PKCϵ can phosphorylate Ser⁷⁷⁹ in vitro, whereas overexpression of PKCϵ results in constitutive Ser⁷⁷⁹ phosphorylation and enhanced PC12 cell differentiation. Furthermore, siRNA knockdown of PKCϵ reduces both growth factor-induced Ser⁷⁷⁹ phosphorylation and neuronal differentiation. Our findings show that in addition to FGFR tyrosine phosphorylation, the phosphorylation of a conserved serine residue, Ser⁷⁷⁹, can quantitatively control Ras/MAPK signaling to promote specific cellular responses.     

Multiple Autoimmune-Associated Variants Confer Decreased IL-2R Signaling in CD4+CD25hi T Cells of Type 1 Diabetic and Multiple Sclerosis Patients

IL-2 receptor (IL-2R) signaling is essential for optimal stability and function of CD4⁺CD25^hiFOXP3⁺ regulatory T cells (Treg); a cell type that plays an integral role in maintaining tolerance. Thus, we hypothesized that decreased response to IL-2 may be a common phenotype of subjects who have autoimmune diseases associated with variants in the IL2RA locus, including T1D and MS, particularly in cells expressing the high affinity IL-2R alpha chain (IL-2RA or CD25). To examine this question we used phosphorylation of STAT5 (pSTAT5) as a downstream measure of IL-2R signaling, and found a decreased response to IL-2 in CD4⁺CD25^hi T cells of T1D and MS, but not SLE patients. Since the IL2RArs2104286 haplotype is associated with T1D and MS, we measured pSTAT5 in controls carrying the rs2104286 risk haplotype to test whether this variant contributed to reduced IL-2 responsiveness. Consistent with this, we found decreased pSTAT5 in subjects carrying the rs2104286 risk haplotype. Reduced IL-2R signaling did not result from lower CD25 expression on CD25^hi cells; instead we detected increased CD25 expression on naive Treg from controls carrying the rs2104286 risk haplotype, and subjects with T1D and MS. However the rs2104286 risk haplotype correlated with increased soluble IL-2RA levels, suggesting that shedding of the IL-2R may account in part for the reduced IL-2R signaling associated with the rs2104286 risk haplotype. In addition to risk variants in IL2RA, we found that the T1D-associated risk variant of PTPN2rs1893217 independently contributed to diminished IL-2R signaling. However, even when holding genotype constant at IL2RA and PTPN2, we still observed a significant signaling defect in T1D and MS patients. Together, these data suggest that multiple mechanisms converge in disease leading to decreased response to IL-2, a phenotype that may eventually lead to loss of tolerance and autoimmunity.     

Advanced Glycation End Products Affect Osteoblast Proliferation and Function by Modulating Autophagy Via the Receptor of Advanced Glycation End Products/Raf Protein/Mitogen-activated Protein Kinase/Extracellular Signal-regulated Kinase Kinase/Extracellular Signal-regulated Kinase (RAGE/Raf/MEK/ERK) Pathway

The interaction between advanced glycation end products (AGEs) and receptor of AGEs (RAGE) is associated with the development and progression of diabetes-associated osteoporosis, but the mechanisms involved are still poorly understood. In this study, we found that AGE-modified bovine serum albumin (AGE-BSA) induced a biphasic effect on the viability of hFOB1.19 cells; cell proliferation was stimulated after exposure to low dose AGE-BSA, but cell apoptosis was stimulated after exposure to high dose AGE-BSA. The low dose AGE-BSA facilitates proliferation of hFOB1.19 cells by concomitantly promoting autophagy, RAGE production, and the Raf/MEK/ERK signaling pathway activation. Furthermore, we investigated the effects of AGE-BSA on the function of hFOB1.19 cells. Interestingly, the results suggest that the short term effects of low dose AGE-BSA increase osteogenic function and decrease osteoclastogenic function, which are likely mediated by autophagy and the RAGE/Raf/MEK/ERK signal pathway. In contrast, with increased treatment time, the opposite effects were observed. Collectively, AGE-BSA had a biphasic effect on the viability of hFOB1.19 cells in vitro, which was determined by the concentration of AGE-BSA and treatment time. A low concentration of AGE-BSA activated the Raf/MEK/ERK signal pathway through the interaction with RAGE, induced autophagy, and regulated the proliferation and function of hFOB1.19 cells.     

Lack of association between miR-146a rs2910164 C/G locus and colorectal cancer: from a case–control study to a meta-analysis

Previous studies suggested that miR-146a rs2910164 (C/G) locus was predicted to influence the risk of cancer. However, the relationship of miR-146a rs2910164 locus with colorectal cancer (CRC) susceptibility was controversial. We recruited 1003 CRC patients and 1303 controls, and performed a case–control study to clarify the correlation of miR-146a rs2910164 locus with CRC risk. Subsequently, a comprehensive meta-analysis was conducted to verify our findings. In the case–control study, we suggested that miR-146a rs2910164 variants did not alter CRC risk (CG vs. CC: adjusted P=0.465; GG vs. CC: adjusted P=0.436, CG/GG vs. CC: adjusted P=0.387 and GG vs. CC/CG: adjusted P=0.589), even in subgroup analysis. Next, we conducted a pooled-analysis to identify the correlation of miR-146a rs2910164 locus with CRC risk. In this pooled-analysis, 7947 CRC cases and 12,168 controls were included. We found that miR-146a rs2910164 polymorphism did not influence the risk of CRC (G vs. C: P=0.537; GG vs. CC: P=0.517, CG/GG vs. CC: P=0.520 and GG vs. CC/CG: P=0.167). Our findings suggest that miR-146a rs2910164 C/G polymorphism is not correlated with the susceptibility of CRC. In the future, more case–control studies are needed to confirm our results.     

Caveolin-1 Is Required for Kinase Suppressor of Ras 1 (KSR1)-Mediated Extracellular Signal-Regulated Kinase 1/2 Activation,H-RasV12-Induced Senescence,and Transformation

The molecular scaffold kinase suppressor of Ras 1 (KSR1) regulates the activation of the Raf/MEK/extracellular signal-regulated kinase (ERK) signal transduction pathway. KSR1 disruption in mouse embryo fibroblasts (MEFs) abrogates growth factor-induced ERK activation, H-Ras^V12-induced replicative senescence, and H-Ras^V12-induced transformation. Caveolin-1 has been primarily described as a major component of the coating structure of caveolae, which can serve as a lipid binding adaptor protein and coordinates the assembly of Ras, Raf, MEK, and ERK. In this study, we show that KSR1 interacts with caveolin-1 and is responsible for MEK and ERK redistribution to caveolin-1-rich fractions. The interaction between KSR1 and caveolin-1 is essential for optimal activation of ERK as a KSR1 mutant unable to interact with caveolin-1 does not efficiently mediate growth factor-induced ERK activation at the early stages of pathway activation. Furthermore, abolishing the KSR1–caveolin-1 interaction increases growth factor demands to promote H-Ras^V12-induced proliferation and has adverse effects on H-Ras^V12-induced cellular senescence and transformation. These data show that caveolin-1 is necessary for optimal KSR1-dependent ERK activation by growth factors and oncogenic Ras.     

Oligonol, an oligomerized lychee fruit-derived polyphenol, activates the Ras/Raf-1/MEK1/2 cascade independent of the IL-6 signaling pathway in rat primary adipocytes

Oligonol is a lychee fruit-derived low-molecular form of polyphenol. In this study, the effect of Oligonol on the mitogen activated-protein kinase (MAPK) signaling pathway in primary adipocytes was investigated to examine the mechanism underlying the enhanced levels of phosphorylated extracellular-signaling regulatory kinase1/2 (ERK1/2) that accompany an in vitro increase in lipolysis. Oligonol significantly elevated the levels of activated Ras and the phosphorylation of Raf-1 and MAPK/ERK kinase1/2 (MEK1/2) with no increase in pan-Raf-1 and -MEK1/2 proteins. The increase in phosphorylation of Raf-1 and MEK1/2 with Oligonol was inhibited completely by pretreatment with GW5074, a selective Raf-1 inhibitor, or PD98059, a selective MEK1/2 inhibitor. IL-6 also activated the MAPK signaling pathway in adipocytes through the association with its receptor. IL-6-induced phosphorylation of Raf-1 and MEK1/2 was significantly inhibited by pretreatment with the IL-6 receptor antibody. Under such a condition, however, the levels of phosphorylated Raf-1 and MEK1/2 with Oligonol still remained significantly higher, and there was a significant decrease in secretion of IL-6 from adipocytes, compared with untreated control cells. These results suggest that Oligonol activates the Ras/Raf-1/MEK1/2 signaling pathway, independent of the IL-6 signaling pathway, leading to activation of ERK1/2 proteins in primary adipocytes.     

Systematic Quantification of Negative Feedback Mechanisms in the Extracellular Signal-regulated Kinase (ERK) Signaling Network

Cell responses are actuated by tightly controlled signal transduction pathways. Although the concept of an integrated signaling network replete with interpathway cross-talk and feedback regulation is broadly appreciated, kinetic data of the type needed to characterize such interactions in conjunction with mathematical models are lacking. In mammalian cells, the Ras/ERK pathway controls cell proliferation and other responses stimulated by growth factors, and several cross-talk and feedback mechanisms affecting its activation have been identified. In this work, we take a systematic approach to parse the magnitudes of multiple regulatory mechanisms that attenuate ERK activation through canonical (Ras-dependent) and non-canonical (PI3K-dependent) pathways. In addition to regulation of receptor and ligand levels, we consider three layers of ERK-dependent feedback: desensitization of Ras activation, negative regulation of MEK kinase (e.g. Raf) activities, and up-regulation of dual-specificity ERK phosphatases. Our results establish the second of these as the dominant mode of ERK self-regulation in mouse fibroblasts. We further demonstrate that kinetic models of signaling networks, trained on a sufficient diversity of quantitative data, can be reasonably comprehensive, accurate, and predictive in the dynamical sense.     

Identification of potential functional variants and genes at 18q21.1 associated with the carcinogenesis of colorectal cancer

Genome-wide association studies (GWAS) have identified more than 160 susceptibility loci for colorectal cancer (CRC). The effects of these variants, particularly their mechanisms, however, remain unclear. In this study, a comprehensive functional annotation of CRC-related GWAS signals was firstly conducted to identify the potential causal variants. We found that the SNP rs7229639 in intron 3 of SMAD7 at 18q21.1 might serve as a putative functional variant in CRC. The SNP rs7229639 is located in a region with evidence of regulatory potential. Dual-luciferase reporter assays revealed that three other SNPs (rs77544449, rs60385309 and rs72917785), in strong linkage disequilibrium (LD) with rs7229639, exhibited allele-specific enhancer activity, of which one of the target genes may conceivably be LIPG, as suggested by eQTL association data and Hi-C data. We also verified that LIPG promoted malignancy of CRC cells in vitro, with supporting clinical data indicating that LIPG is upregulated and correlated with a poor prognosis in CRC. Finally, pitavastatin was observed to exhibit an anti-CRC activity and modest inhibition of LIPG mRNA levels. Collectively, our data suggest that these functional variants at 18q21.1 are involved in the pathogenesis of CRC by modulating enhancer activity, and possibly LIPG expression, thus indicating a promising therapeutic target for CRC. The results of functional annotation in our investigation could also serve as an inventory for CRC susceptibility SNPs and offer guides for post-GWAS downstream functional studies.