BMPR2 as a Novel Predisposition Gene for Hereditary Colorectal Polyposis.

BACKGROUND & AIMS: Colorectal cancer (CRC) is one of the most prevalent tumors worldwide, with incidence quickly increasing (particularly in the context of early-onset cases), despite important prevention efforts, mainly in the form of population-wide screening programs. Although many cases present a clear familial component, the current list of hereditary CRC genes leaves a considerable proportion of the cases unexplained. METHODS: In this work, we used whole-exome sequencing approaches on 19 unrelated patients with unexplained colonic polyposis to identify candidate CRC predisposition genes. The candidate genes were then validated in an additional series of 365 patients. CRISPR-Cas9 models were used to validate BMPR2 as a potential candidate for CRC risk. RESULTS: We found 8 individuals carrying 6 different variants in the BMPR2 gene (approximately 2% of our cohort of patients with unexplained colonic polyposis). CRISPR-Cas9 models of 3 of these variants showed that the p.(Asn442Thrfs∗32) truncating variant completely abrogated BMP pathway function in a similar way to the BMPR2 knockout. Missense variants p.(Asn565Ser), p.(Ser967Pro) had varying effects on cell proliferation levels, with the former impairing cell control inhibition via noncanonical pathways. CONCLUSIONS: Collectively, these results support loss-of-function BMPR2 variants as candidates to be involved in CRC germline predisposition.

Characterization of additive gene-environment interactions for colorectal cancer risk.

BACKGROUND: Colorectal cancer (CRC) is a common, fatal cancer. Identifying subgroups who may benefit more from intervention is of critical public health importance. Previous studies have assessed multiplicative interaction between genetic risk scores and environmental factors, but few have assessed additive interaction, the relevant public health measure. METHODS: Using resources from colorectal cancer consortia including 45,247 CRC cases and 52,671 controls, we assessed multiplicative and additive interaction (relative excess risk due to interaction, RERI) using logistic regression between 13 harmonized environmental factors and genetic risk score including 141 variants associated with CRC risk. RESULTS: There was no evidence of multiplicative interaction between environmental factors and genetic risk score. There was additive interaction where, for individuals with high genetic susceptibility, either heavy drinking [RERI = 0.24, 95% confidence interval, CI, (0.13, 0.36)], ever smoking [0.11 (0.05, 0.16)], high BMI [female 0.09 (0.05, 0.13), male 0.10 (0.05, 0.14)], or high red meat intake [highest versus lowest quartile 0.18 (0.09, 0.27)] was associated with excess CRC risk greater than that for individuals with average genetic susceptibility. Conversely, we estimate those with high genetic susceptibility may benefit more from reducing CRC risk with aspirin/NSAID use [-0.16 (-0.20, -0.11)] or higher intake of fruit, fiber, or calcium [highest quartile versus lowest quartile -0.12 (-0.18, -0.050); -0.16 (-0.23, -0.09); -0.11 (-0.18, -0.05), respectively] than those with average genetic susceptibility. CONCLUSIONS: Additive interaction is important to assess for identifying subgroups who may benefit from intervention. The subgroups identified in this study may help inform precision CRC prevention.

CTNND1 is involved in germline predisposition to early-onset gastric cancer by affecting cell-to-cell interactions.

BACKGROUND: CDH1 and CTNNA1 remain as the main genes for hereditary gastric cancer. However, they only explain a small fraction of gastric cancer cases with suspected inherited basis. In this study, we aimed to identify new hereditary genes for early-onset gastric cancer patients (EOGC; < 50 years old). METHODS: After germline exome sequencing in 20 EOGC patients and replication of relevant findings by gene-panel sequencing in an independent cohort of 152 patients, CTNND1 stood out as an interesting candidate gene, since its protein product (p120ctn) directly interacts with E-cadherin. We proceeded with functional characterization by generating two knockout CTNND1 cellular models by gene editing and introducing the detected genetic variants using a lentiviral delivery system. We assessed ß-catenin and E-cadherin levels, cell detachment, as well as E-cadherin localization and cell-to-cell interaction by spheroid modeling. RESULTS: Three CTNND1 germline variants [c.28_29delinsCT, p.(Ala10Leu); c.1105C > T, p.(Pro369Ser); c.1537A > G, p.(Asn513Asp)] were identified in our EOGC cohorts. Cells encoding CTNND1 variants displayed altered E-cadherin levels and intercellular interactions. In addition, the p.(Pro369Ser) variant, located in a key region in the E-cadherin/p120ctn binding domain, showed E-cadherin mislocalization. CONCLUSIONS: Defects in CTNND1 could be involved in germline predisposition to gastric cancer by altering E-cadherin and, consequently, cell-to-cell interactions. In the present study, CTNND1 germline variants explained 2% (3/172) of the cases, although further studies in larger external cohorts are needed.

Germline mutations in WNK2 could be associated with serrated polyposis syndrome.

BACKGROUND: Patients with serrated polyposis syndrome (SPS) have multiple and/or large serrated colonic polyps and higher risk for colorectal cancer. SPS inherited genetic basis is mostly unknown. We aimed to identify new germline predisposition factors for SPS by functionally evaluating a candidate gene and replicating it in additional SPS cohorts. METHODS: After a previous whole-exome sequencing in 39 SPS patients from 16 families (discovery cohort), we sequenced specific genes in an independent validation cohort of 211 unrelated SPS cases. Additional external replication was also available in 297 SPS cases. The WNK2 gene was disrupted in HT-29 cells by gene editing, and WNK2 variants were transfected using a lentiviral delivery system. Cells were analysed by immunoblots, real-time PCR and functional assays monitoring the mitogen-activated protein kinase (MAPK) pathway, cell cycle progression, survival and adhesion. RESULTS: We identified 2 rare germline variants in the WNK2 gene in the discovery cohort, 3 additional variants in the validation cohort and 10 other variants in the external cohorts. Variants c.2105C>T (p.Pro702Leu), c.4820C>T (p.Ala1607Val) and c.6157G>A (p.Val2053Ile) were functionally characterised, displaying higher levels of phospho-PAK1/2, phospho-ERK1/2, CCND1, clonogenic capacity and MMP2. CONCLUSION: After whole-exome sequencing in SPS cases with familial aggregation and replication of results in additional cohorts, we identified rare germline variants in the WNK2 gene. Functional studies suggested germline WNK2 variants affect protein function in the context of the MAPK pathway, a molecular hallmark in this disease.

Data-Independent Acquisition Mass Spectrometry Analysis of FFPE Rectal Cancer Samples Offers In-Depth Proteomics Characterization of the Response to Neoadjuvant Chemoradiotherapy.

Locally advanced rectal cancer (LARC) presents a challenge in identifying molecular markers linked to the response to neoadjuvant chemoradiotherapy (nCRT). This study aimed to utilize a sensitive proteomic method, data-independent mass spectrometry (DIA-MS), to extensively analyze the LARC proteome, seeking individuals with favorable initial responses suitable for a watch-and-wait approach. This research addresses the unmet need to understand the response to treatment, potentially guiding personalized strategies for LARC patients. Post-treatment assessment included MRI scans and proctoscopy. This research involved 97 LARC patients treated with intense chemoradiotherapy, comprising radiation and chemotherapy. Out of 97 LARC included in this study, we selected 20 samples with the most different responses to nCRT for proteome profiling (responders vs. non-responders). This proteomic approach shows extensive proteome coverage in LARC samples. The analysis identified a significant number of proteins compared to a prior study. A total of 915 proteins exhibited differential expression between the two groups, with certain signaling pathways associated with response mechanisms, while top candidates had good predictive potential. Proteins encoded by genes SMPDL3A, PCTP, LGMN, SYNJ2, NHLRC3, GLB1, and RAB43 showed high predictive potential of unfavorable treatment outcome, while RPA2, SARNP, PCBP2, SF3B2, HNRNPF, RBBP4, MAGOHB, DUT, ERG28, and BUB3 were good predictive biomarkers of favorable treatment outcome. The identified proteins and related biological processes provide promising insights that could enhance the management and care of LARC patients.

Identifying Novel Susceptibility Genes for Colorectal Cancer Risk From a Transcriptome-Wide Association Study of 125,478 Subjects.

BACKGROUND AND AIMS: Susceptibility genes and the underlying mechanisms for the majority of risk loci identified by genome-wide association studies (GWAS) for colorectal cancer (CRC) risk remain largely unknown. We conducted a transcriptome-wide association study (TWAS) to identify putative susceptibility genes. METHODS: Gene-expression prediction models were built using transcriptome and genetic data from the 284 normal transverse colon tissues of European descendants from the Genotype-Tissue Expression (GTEx), and model performance was evaluated using data from The Cancer Genome Atlas (n = 355). We applied the gene-expression prediction models and GWAS data to evaluate associations of genetically predicted gene-expression with CRC risk in 58,131 CRC cases and 67,347 controls of European ancestry. Dual-luciferase reporter assays and knockdown experiments in CRC cells and tumor xenografts were conducted. RESULTS: We identified 25 genes associated with CRC risk at a Bonferroni-corrected threshold of P < 9.1 × 10-6, including genes in 4 novel loci, PYGL (14q22.1), RPL28 (19q13.42), CAPN12 (19q13.2), MYH7B (20q11.22), and MAP1L3CA (20q11.22). In 9 known GWAS-identified loci, we uncovered 9 genes that have not been reported previously, whereas 4 genes remained statistically significant after adjusting for the lead risk variant of the locus. Through colocalization analysis in GWAS loci, we additionally identified 12 putative susceptibility genes that were supported by TWAS analysis at P < .01. We showed that risk allele of the lead risk variant rs1741640 affected the promoter activity of CABLES2. Knockdown experiments confirmed that CABLES2 plays a vital role in colorectal carcinogenesis. CONCLUSIONS: Our study reveals new putative susceptibility genes and provides new insight into the biological mechanisms underlying CRC development.

Germline Mutations in FAF1 Are Associated With Hereditary Colorectal Cancer.

BACKGROUND & AIMS: A significant proportion of colorectal cancer (CRC) cases have familial aggregation but little is known about the genetic factors that contribute to these cases. We performed an exhaustive functional characterization of genetic variants associated with familial CRC. METHODS: We performed whole-exome sequencing analyses of 75 patients from 40 families with a history of CRC (including early-onset cases) of an unknown germline basis (discovery cohort). We also sequenced specific genes in DNA from an external replication cohort of 473 families, including 488 patients with colorectal tumors that had normal expression of mismatch repair proteins (validation cohort). We disrupted the Fas-associated factor 1 gene (FAF1) in DLD-1 CRC cells using CRISPR/Cas9 gene editing; some cells were transfected with plasmids that express FAF1 missense variants. Cells were analyzed by immunoblots, quantitative real-time polymerase chain reaction, and functional assays monitoring apoptosis, proliferation, and assays for Wnt signaling or nuclear factor (NF)-kappa-B activity. RESULTS: We identified predicted pathogenic variant in the FAF1 gene (c.1111G>A; p.Asp371Asn) in the discovery cohort; it was present in 4 patients of the same family. We identified a second variant in FAF1 in the validation cohort (c.254G>C; p.Arg85Pro). Both variants encoded unstable FAF1 proteins. Expression of these variants in CRC cells caused them to become resistant to apoptosis, accumulate beta-catenin in the cytoplasm, and translocate NF-kappa-B to the nucleus. CONCLUSIONS: In whole-exome sequencing analyses of patients from families with a history of CRC, we identified variants in FAF1 that associate with development of CRC. These variants encode unstable forms of FAF1 that increase resistance of CRC cells to apoptosis and increase activity of beta-catenin and NF-kappa-B.

Circulating Levels of Insulin-like Growth Factor 1 and Insulin-like Growth Factor Binding Protein 3 Associate With Risk of Colorectal Cancer Based on Serologic and Mendelian Randomization Analyses.

BACKGROUND & AIMS: Human studies examining associations between circulating levels of insulin-like growth factor 1 (IGF1) and insulin-like growth factor binding protein 3 (IGFBP3) and colorectal cancer risk have reported inconsistent results. We conducted complementary serologic and Mendelian randomization (MR) analyses to determine whether alterations in circulating levels of IGF1 or IGFBP3 are associated with colorectal cancer development. METHODS: Serum levels of IGF1 were measured in blood samples collected from 397,380 participants from the UK Biobank, from 2006 through 2010. Incident cancer cases and cancer cases recorded first in death certificates were identified through linkage to national cancer and death registries. Complete follow-up was available through March 31, 2016. For the MR analyses, we identified genetic variants associated with circulating levels of IGF1 and IGFBP3. The association of these genetic variants with colorectal cancer was examined with 2-sample MR methods using genome-wide association study consortia data (52,865 cases with colorectal cancer and 46,287 individuals without [controls]) RESULTS: After a median follow-up period of 7.1 years, 2665 cases of colorectal cancer were recorded. In a multivariable-adjusted model, circulating level of IGF1 associated with colorectal cancer risk (hazard ratio per 1 standard deviation increment of IGF1, 1.11; 95% confidence interval [CI] 1.05-1.17). Similar associations were found by sex, follow-up time, and tumor subsite. In the MR analyses, a 1 standard deviation increment in IGF1 level, predicted based on genetic factors, was associated with a higher risk of colorectal cancer risk (odds ratio 1.08; 95% CI 1.03-1.12; P = 3.3 × 10-4). Level of IGFBP3, predicted based on genetic factors, was associated with colorectal cancer risk (odds ratio per 1 standard deviation increment, 1.12; 95% CI 1.06-1.18; P = 4.2 × 10-5). Colorectal cancer risk was associated with only 1 variant in the IGFBP3 gene region (rs11977526), which also associated with anthropometric traits and circulating level of IGF2. CONCLUSIONS: In an analysis of blood samples from almost 400,000 participants in the UK Biobank, we found an association between circulating level of IGF1 and colorectal cancer. Using genetic data from 52,865 cases with colorectal cancer and 46,287 controls, a higher level of IGF1, determined by genetic factors, was associated with colorectal cancer. Further studies are needed to determine how this signaling pathway might contribute to colorectal carcinogenesis.

Colorectal cancer genetic variants are also associated with serrated polyposis syndrome susceptibility.

BACKGROUND: Serrated polyposis syndrome (SPS) is a clinical entity characterised by large and/ormultiple serrated polyps throughout the colon and increased risk for colorectal cancer (CRC). The basis for SPS genetic predisposition is largely unknown. Common, low-penetrance genetic variants have been consistently associated with CRC susceptibility, however, their role in SPS genetic predisposition has not been yet explored. OBJECTIVE: The aim of this study was to evaluate if common, low-penetrance genetic variants for CRC risk are also implicated in SPS genetic susceptibility. METHODS: A case-control study was performed in 219 SPS patients and 548 asymptomatic controls analysing 65 CRC susceptibility variants. A risk prediction model for SPS predisposition was developed. RESULTS: Statistically significant associations with SPS were found for seven genetic variants (rs4779584-GREM1, rs16892766-EIF3H, rs3217810-CCND2, rs992157-PNKD1/TMBIM1, rs704017-ZMIZ1, rs11196172-TCF7L2, rs6061231-LAMA5). The GREM1 risk allele was remarkably over-represented in SPS cases compared with controls (OR=1.573, 1.21-2.04, p value=0.0006). A fourfold increase in SPS risk was observed when comparing subjects within the highest decile of variants (≥65) with those in the first decile (≤50). CONCLUSIONS: Genetic variants for CRC risk are also involved in SPS susceptibility, being the most relevant ones rs4779584-GREM1, rs16892766-EIF3H and rs3217810-CCND2.

Identification of New Genes Involved in Germline Predisposition to Early-Onset Gastric Cancer.

The genetic cause for several families with gastric cancer (GC) aggregation is unclear, with marked relevance in early-onset patients. We aimed to identify new candidate genes involved in GC germline predisposition. Whole-exome sequencing (WES) of germline samples was performed in 20 early-onset GC patients without previous germline mutation identified. WES was also performed in nine tumor samples to analyze the somatic profile using SigProfilerExtractor tool. Sequencing germline data were filtered to select those variants with plausible pathogenicity, rare frequency and previously involved in cancer. Then, a manual filtering was performed to prioritize genes according to current knowledge and function. These genetic variants were prevalidated with Integrative Genomics Viewer 2.8.2 (IGV). Subsequently, a further selection step was carried out according to function and information obtained from tumor samples. After IGV and selection step, 58 genetic variants in 52 different candidate genes were validated by Sanger sequencing. Among them, APC, FAT4, CTNND1 and TLR2 seem to be the most promising genes because of their role in hereditary cancer syndromes, tumor suppression, cell adhesion and Helicobacter pylori recognition, respectively. These encouraging results represent the open door to the identification of new genes involved in GC germline predisposition.

Prognostic Factors for Post-Recurrence Survival in Stage II and III Colorectal Carcinoma Patients.

Background and objectives: This study aimed to evaluate prognostic factors for post-recurrence survival in local and locally advanced colorectal cancer patients. Materials and Methods: A total of 273 patients with stage III and high-risk stage II colorectal cancer were prospectively enrolled. All patients underwent operative treatment of the primary tumor and adjuvant fluorouracil-based chemotherapy. Results: Over the three-year period (2008-2010), a cohort of 273 patients with stage III and high-risk stage II colorectal cancer had been screened. During follow up, 105 (38.5%) patients had disease recurrence. Survival rates 1-, 3- and 5-year after recurrence were 53.9, 18.2 and 6.5%, respectively, and the median post-recurrence survival time was 13 months. Survival analysis showed that age at diagnosis (p < 0.01), gender (p < 0.05), elevated postoperative Ca19-9 (p < 0.01), tumor histology (adenocarcinoma vs. mucinous vs. signet ring tumors, p < 0.01) and tumor stage (II vs. III, p < 0.05) had a significant influence on post-recurrence survival. Recurrence interval and metastatic site were not related to survival following recurrence. Multivariate analysis showed that older age (HR 2.43), mucinous tumors (HR 1.51) and tumors expressing Ca19-9 at baseline (HR 3.51) were independently associated with survival following recurrence. Conclusions: Baseline patient and tumor characteristics largely predicted patient outcomes after disease recurrence. Recurrence intervals in local and locally advanced colorectal cancer were not found to be prognostic factors for post-recurrence survival. Older age, male gender, stage III and mucinous histology were poor prognostic factors after the disease had recurred. Stage II patients had remarkable post-recurrence survival compared to stage III patients.

Systematic meta-analyses, field synopsis and global assessment of the evidence of genetic association studies in colorectal cancer.

OBJECTIVE: To provide an understanding of the role of common genetic variations in colorectal cancer (CRC) risk, we report an updated field synopsis and comprehensive assessment of evidence to catalogue all genetic markers for CRC (CRCgene2). DESIGN: We included 869 publications after parallel literature review and extracted data for 1063 polymorphisms in 303 different genes. Meta-analyses were performed for 308 single nucleotide polymorphisms (SNPs) in 158 different genes with at least three independent studies available for analysis. Scottish, Canadian and Spanish data from genome-wide association studies (GWASs) were incorporated for the meta-analyses of 132 SNPs. To assess and classify the credibility of the associations, we applied the Venice criteria and Bayesian False-Discovery Probability (BFDP). Genetic associations classified as 'positive' and 'less-credible positive' were further validated in three large GWAS consortia conducted in populations of European origin. RESULTS: We initially identified 18 independent variants at 16 loci that were classified as 'positive' polymorphisms for their highly credible associations with CRC risk and 59 variants at 49 loci that were classified as 'less-credible positive' SNPs; 72.2% of the 'positive' SNPs were successfully replicated in three large GWASs and the ones that were not replicated were downgraded to 'less-credible' positive (reducing the 'positive' variants to 14 at 11 loci). For the remaining 231 variants, which were previously reported, our meta-analyses found no evidence to support their associations with CRC risk. CONCLUSION: The CRCgene2 database provides an updated list of genetic variants related to CRC risk by using harmonised methods to assess their credibility.

Using linkage studies combined with whole-exome sequencing to identify novel candidate genes for familial colorectal cancer.

Colorectal cancer (CRC) is a complex disorder for which the majority of the underlying germline predisposition factors remain still unidentified. Here, we combined whole-exome sequencing (WES) and linkage analysis in families with multiple relatives affected by CRC to identify candidate genes harboring rare variants with potential high-penetrance effects. Forty-seven affected subjects from 18 extended CRC families underwent WES. Genome-wide linkage analysis was performed under linear and exponential models. Suggestive linkage peaks were identified on chromosomes 1q22-q24.2 (maxSNP = rs2134095; LODlinear = 2.38, LODexp = 2.196), 7q31.2-q34 (maxSNP = rs6953296; LODlinear = 2.197, LODexp = 2.149) and 10q21.2-q23.1 (maxSNP = rs1904589; LODlinear = 1.445, LODexp = 2.195). These linkage signals were replicated in 10 independent sets of random markers from each of these regions. To assess the contribution of rare variants predicted to be pathogenic, we performed a family-based segregation test with 89 rare variants predicted to be deleterious from 78 genes under the linkage intervals. This analysis showed significant segregation of rare variants with CRC in 18 genes (weighted p-value > 0.0028). Protein network analysis and functional evaluation were used to suggest a plausible candidate gene for germline CRC predisposition. Etiologic rare variants implicated in cancer germline predisposition may be identified by combining traditional linkage with WES data. This approach can be used with already available NGS data from families with several sequenced members to further identify candidate genes involved germline predisposition to disease. This approach resulted in one candidate gene associated with increased risk of CRC but needs evidence from further studies.

Adiposity, metabolites, and colorectal cancer risk: Mendelian randomization study.

BACKGROUND: Higher adiposity increases the risk of colorectal cancer (CRC), but whether this relationship varies by anatomical sub-site or by sex is unclear. Further, the metabolic alterations mediating the effects of adiposity on CRC are not fully understood. METHODS: We examined sex- and site-specific associations of adiposity with CRC risk and whether adiposity-associated metabolites explain the associations of adiposity with CRC. Genetic variants from genome-wide association studies of body mass index (BMI) and waist-to-hip ratio (WHR, unadjusted for BMI; N = 806,810), and 123 metabolites from targeted nuclear magnetic resonance metabolomics (N = 24,925), were used as instruments. Sex-combined and sex-specific Mendelian randomization (MR) was conducted for BMI and WHR with CRC risk (58,221 cases and 67,694 controls in the Genetics and Epidemiology of Colorectal Cancer Consortium, Colorectal Cancer Transdisciplinary Study, and Colon Cancer Family Registry). Sex-combined MR was conducted for BMI and WHR with metabolites, for metabolites with CRC, and for BMI and WHR with CRC adjusted for metabolite classes in multivariable models. RESULTS: In sex-specific MR analyses, higher BMI (per 4.2 kg/m2) was associated with 1.23 (95% confidence interval (CI) = 1.08, 1.38) times higher CRC odds among men (inverse-variance-weighted (IVW) model); among women, higher BMI (per 5.2 kg/m2) was associated with 1.09 (95% CI = 0.97, 1.22) times higher CRC odds. WHR (per 0.07 higher) was more strongly associated with CRC risk among women (IVW OR = 1.25, 95% CI = 1.08, 1.43) than men (IVW OR = 1.05, 95% CI = 0.81, 1.36). BMI or WHR was associated with 104/123 metabolites at false discovery rate-corrected P ≤ 0.05; several metabolites were associated with CRC, but not in directions that were consistent with the mediation of positive adiposity-CRC relations. In multivariable MR analyses, associations of BMI and WHR with CRC were not attenuated following adjustment for representative metabolite classes, e.g., the univariable IVW OR for BMI with CRC was 1.12 (95% CI = 1.00, 1.26), and this became 1.11 (95% CI = 0.99, 1.26) when adjusting for cholesterol in low-density lipoprotein particles. CONCLUSIONS: Our results suggest that higher BMI more greatly raises CRC risk among men, whereas higher WHR more greatly raises CRC risk among women. Adiposity was associated with numerous metabolic alterations, but none of these explained associations between adiposity and CRC. More detailed metabolomic measures are likely needed to clarify the mechanistic pathways.

Quantitative analysis of somatically acquired and constitutive uniparental disomy in gastrointestinal cancers.

Somatically acquired uniparental disomies (aUPDs) are frequent events in solid tumors and have been associated with cancer-related genes. Studies assessing their functional consequences across several cancer types are therefore necessary. Here, we aimed at integrating aUPD profiles with the mutational status of cancer-related genes in a tumor-type specific manner. Using TCGA datasets for 1,032 gastrointestinal cancers, including colon (COAD), rectum (READ), stomach (STAD), esophageal adenocarcinoma (EAC) and esophageal squamous cell carcinoma (ESCC), we show a non-random distribution of aUPD, suggesting the existence of a cancer-specific landscape of aUPD events. Our analysis indicates that aUPD acts as a "second hit" in Knudson's model in order to achieve biallelic inactivation of tumor suppressor genes. In particular, APC, ARID1A and NOTCH1 were recurrently inactivated by the presence of homozygous mutation as a consequence of aUPD in COAD and READ, STAD and ESCC, respectively. Furthermore, while TP53 showed inactivation caused by aUPD at chromosome arm 17p across all tumor types, copy number losses at this genomic position were also frequent. By experimental and computationally inferring genome ploidy, we demonstrate that an increased number of aUPD events, both affecting the whole chromosome or segments of it, were present in highly aneuploid genomes compared to near-diploid tumors. Finally, the presence of mosaic UPD was detected at a higher frequency in DNA extracted from peripheral blood lymphocytes of patients with colorectal cancer compared to healthy individuals. In summary, our study defines specific profiles of aUPD in gastrointestinal cancers and provides unequivocal evidence of their relevance in cancer.

GeMSTONE: orchestrated prioritization of human germline mutations in the cloud.

Integrative analysis of whole-genome/exome-sequencing data has been challenging, especially for the non-programming research community, as it requires simultaneously managing a large number of computational tools. Even computational biologists find it unexpectedly difficult to reproduce results from others or optimize their strategies in an end-to-end workflow. We introduce Germline Mutation Scoring Tool fOr Next-generation sEquencing data (GeMSTONE), a cloud-based variant prioritization tool with high-level customization and a comprehensive collection of bioinformatics tools and data libraries (http://gemstone.yulab.org/). GeMSTONE generates and readily accepts a shareable 'recipe' file for each run to either replicate previous results or analyze new data with identical parameters and provides a centralized workflow for prioritizing germline mutations in human disease within a streamlined workflow rather than a pool of program executions.

Mutational Signatures in Cancer (MuSiCa): a web application to implement mutational signatures analysis in cancer samples.

BACKGROUND: Mutational signatures have been proved as a valuable pattern in somatic genomics, mainly regarding cancer, with a potential application as a biomarker in clinical practice. Up to now, several bioinformatic packages to address this topic have been developed in different languages/platforms. MutationalPatterns has arisen as the most efficient tool for the comparison with the signatures currently reported in the Catalogue of Somatic Mutations in Cancer (COSMIC) database. However, the analysis of mutational signatures is nowadays restricted to a small community of bioinformatic experts. RESULTS: In this work we present Mutational Signatures in Cancer (MuSiCa), a new web tool based on MutationalPatterns and built using the Shiny framework in R language. By means of a simple interface suited to non-specialized researchers, it provides a comprehensive analysis of the somatic mutational status of the supplied cancer samples. It permits characterizing the profile and burden of mutations, as well as quantifying COSMIC-reported mutational signatures. It also allows classifying samples according to the above signature contributions. CONCLUSIONS: MuSiCa is a helpful web application to characterize mutational signatures in cancer samples. It is accessible online at http://bioinfo.ciberehd.org/GPtoCRC/en/tools.html and source code is freely available at https://github.com/marcos-diazg/musica .

Colon-specific eQTL analysis to inform on functional SNPs.

BACKGROUND: Genome-wide association studies on colorectal cancer have identified more than 60 susceptibility loci, but for most of them there is no clear knowledge of functionality or the underlying gene responsible for the risk modification. Expression quantitative trail loci (eQTL) may provide functional information for such single nucleotide polymorphisms (SNPs). METHODS: We have performed detailed eQTL analysis specific for colon tissue on a series of 97 colon tumours, their paired adjacent normal mucosa and 47 colon mucosa samples donated by healthy individuals. R package MatrixEQTL was used to search for genome-wide cis-eQTL and trans-eQTL fitting linear models adjusted for age, gender and tissue type to rank transformed expression data. RESULTS: The cis-eQTL analyses has revealed 29,073 SNP-gene associations with permutation-adjusted P-values < 0.01. These correspond to 363 unique genes. The trans-eQTL analysis identified 10,665 significant SNP-gene associations, most of them in the same chromosome, further than 1 Mb of the gene. We provide a web tool to search for specific SNPs or genes. The tool calculates Pearson or Spearman correlation, and allows to select tissue type for analysis. Data and plots can be exported. CONCLUSIONS: This resource should be useful to prioritise SNPs for further functional studies and to identify relevant genes behind identified loci.

Fine-mapping of the HNF1B multicancer locus identifies candidate variants that mediate endometrial cancer risk.

Common variants in the hepatocyte nuclear factor 1 homeobox B (HNF1B) gene are associated with the risk of Type II diabetes and multiple cancers. Evidence to date indicates that cancer risk may be mediated via genetic or epigenetic effects on HNF1B gene expression. We previously found single-nucleotide polymorphisms (SNPs) at the HNF1B locus to be associated with endometrial cancer, and now report extensive fine-mapping and in silico and laboratory analyses of this locus. Analysis of 1184 genotyped and imputed SNPs in 6608 Caucasian cases and 37 925 controls, and 895 Asian cases and 1968 controls, revealed the best signal of association for SNP rs11263763 (P = 8.4 × 10(-14), odds ratio = 0.86, 95% confidence interval = 0.82-0.89), located within HNF1B intron 1. Haplotype analysis and conditional analyses provide no evidence of further independent endometrial cancer risk variants at this locus. SNP rs11263763 genotype was associated with HNF1B mRNA expression but not with HNF1B methylation in endometrial tumor samples from The Cancer Genome Atlas. Genetic analyses prioritized rs11263763 and four other SNPs in high-to-moderate linkage disequilibrium as the most likely causal SNPs. Three of these SNPs map to the extended HNF1B promoter based on chromatin marks extending from the minimal promoter region. Reporter assays demonstrated that this extended region reduces activity in combination with the minimal HNF1B promoter, and that the minor alleles of rs11263763 or rs8064454 are associated with decreased HNF1B promoter activity. Our findings provide evidence for a single signal associated with endometrial cancer risk at the HNF1B locus, and that risk is likely mediated via altered HNF1B gene expression.