To understand the genetic drivers of immune recognition and evasion in colorectal cancer, we analyzed 1,211 colorectal cancer primary tumor samples, including 179 classified as microsatellite ...instability-high (MSI-high). This set includes The Cancer Genome Atlas colorectal cancer cohort of 592 samples, completed and analyzed here. MSI-high, a hypermutated, immunogenic subtype of colorectal cancer, had a high rate of significantly mutated genes in important immune-modulating pathways and in the antigen presentation machinery, including biallelic losses of
and
genes due to copy-number alterations and copy-neutral loss of heterozygosity. WNT/β-catenin signaling genes were significantly mutated in all colorectal cancer subtypes, and activated WNT/β-catenin signaling was correlated with the absence of T-cell infiltration. This large-scale genomic analysis of colorectal cancer demonstrates that MSI-high cases frequently undergo an immunoediting process that provides them with genetic events allowing immune escape despite high mutational load and frequent lymphocytic infiltration and, furthermore, that colorectal cancer tumors have genetic and methylation events associated with activated WNT signaling and T-cell exclusion.
This multi-omic analysis of 1,211 colorectal cancer primary tumors reveals that it should be possible to better monitor resistance in the 15% of cases that respond to immune blockade therapy and also to use WNT signaling inhibitors to reverse immune exclusion in the 85% of cases that currently do not.
.
Guidelines for initiating colorectal cancer (CRC) screening are based on family history but do not consider lifestyle, environmental, or genetic risk factors. We developed models to determine risk of ...CRC, based on lifestyle and environmental factors and genetic variants, and to identify an optimal age to begin screening.
We collected data from 9748 CRC cases and 10,590 controls in the Genetics and Epidemiology of Colorectal Cancer Consortium and the Colorectal Transdisciplinary study, from 1992 through 2005. Half of the participants were used to develop the risk determination model and the other half were used to evaluate the discriminatory accuracy (validation set). Models of CRC risk were created based on family history, 19 lifestyle and environmental factors (E-score), and 63 CRC-associated single-nucleotide polymorphisms identified in genome-wide association studies (G-score). We evaluated the discriminatory accuracy of the models by calculating area under the receiver operating characteristic curve values, adjusting for study, age, and endoscopy history for the validation set. We used the models to project the 10-year absolute risk of CRC for a given risk profile and recommend ages to begin screening in comparison to CRC risk for an average individual at 50 years of age, using external population incidence rates for non-Hispanic whites from the Surveillance, Epidemiology, and End Results program registry.
In our models, E-score and G-score each determined risk of CRC with greater accuracy than family history. A model that combined both scores and family history estimated CRC risk with an area under the receiver operating characteristic curve value of 0.63 (95% confidence interval, 0.62–0.64) for men and 0.62 (95% confidence interval, 0.61–0.63) for women; area under the receiver operating characteristic curve values based on only family history ranged from 0.53 to 0.54 and those based only E-score or G-score ranged from 0.59 to 0.60. Although screening is recommended to begin at age 50 years for individuals with no family history of CRC, starting ages calculated based on combined E-score and G-score differed by 12 years for men and 14 for women, for individuals with the highest vs the lowest 10% of risk.
We used data from 2 large international consortia to develop CRC risk calculation models that included genetic and environmental factors along with family history. These determine risk of CRC and starting ages for screening with greater accuracy than the family history only model, which is based on the current screening guideline. These scoring systems might serve as a first step toward developing individualized CRC prevention strategies.
Display omitted
Genome-wide association studies (GWAS) have successfully identified tens of thousands of genetic variants associated with various phenotypes, but together they explain only a fraction of ...heritability, suggesting many variants have yet to be discovered. Recently it has been recognized that incorporating functional information of genetic variants can improve power for identifying novel loci. For example, S-PrediXcan and TWAS tested the association of predicted gene expression with phenotypes based on GWAS summary statistics by leveraging the information on genetic regulation of gene expression and found many novel loci. However, as genetic variants may have effects on more than one gene and through different mechanisms, these methods likely only capture part of the total effects of these variants. In this paper, we propose a summary statistics-based mixed effects score test (sMiST) that tests for the total effect of both the effect of the mediator by imputing genetically predicted gene expression, like S-PrediXcan and TWAS, and the direct effects of individual variants. It allows for multiple functional annotations and multiple genetically predicted mediators. It can also perform conditional association analysis while adjusting for other genetic variants (e.g., known loci for the phenotype). Extensive simulation and real data analyses demonstrate that sMiST yields p-values that agree well with those obtained from individual level data but with substantively improved computational speed. Importantly, a broad application of sMiST to GWAS is possible, as only summary statistics of genetic variant associations are required. We apply sMiST to a large-scale GWAS of colorectal cancer using summary statistics from ∼120, 000 study participants and gene expression data from the Genotype-Tissue Expression (GTEx) project. We identify several novel and secondary independent genetic loci.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Genome-wide association studies (GWASs) have identified more than 70 loci associated with type 2 diabetes (T2D), but for most, the underlying causal variants, associated genes, and functional ...mechanisms remain unknown. At a T2D- and fasting-proinsulin-associated locus on 11q13.4, we have identified a functional regulatory DNA variant, a candidate target gene, and a plausible underlying molecular mechanism. Fine mapping, conditional analyses, and exome array genotyping in 8,635 individuals from the Metabolic Syndrome in Men study confirmed a single major association signal between fasting proinsulin and noncoding variants (p = 7.4 × 10−50). Measurement of allele-specific mRNA levels in human pancreatic islet samples heterozygous for rs11603334 showed that the T2D-risk and proinsulin-decreasing allele (C) is associated with increased ARAP1 expression (p < 0.02). We evaluated four candidate functional SNPs for allelic effects on transcriptional activity by performing reporter assays in rodent pancreatic beta cell lines. The C allele of rs11603334, located near one of the ARAP1 promoters, exhibited 2-fold higher transcriptional activity than did the T allele (p < 0.0001); three other candidate SNPs showed no allelic differences. Electrophoretic mobility shift assays demonstrated decreased binding of pancreatic beta cell transcriptional regulators PAX6 and PAX4 to the rs11603334 C allele. Collectively, these data suggest that the T2D-risk allele of rs11603334 could abrogate binding of a complex containing PAX6 and PAX4 and thus lead to increased promoter activity and ARAP1 expression in human pancreatic islets. This work suggests that increased ARAP1 expression might contribute to T2D susceptibility at this GWAS locus.
Lipid and lipoprotein subclasses are associated with metabolic and cardiovascular diseases, yet the genetic contributions to variability in subclass traits are not fully understood. We conducted ...single-variant and gene-based association tests between 15.1M variants from genome-wide and exome array and imputed genotypes and 72 lipid and lipoprotein traits in 8,372 Finns. After accounting for 885 variants at 157 previously identified lipid loci, we identified five novel signals near established loci at HIF3A, ADAMTS3, PLTP, LCAT, and LIPG. Four of the signals were identified with a low-frequency (0.005<minor allele frequency MAF<0.05) or rare (MAF<0.005) variant, including Arg123His in LCAT. Gene-based associations (P<10-10) support a role for coding variants in LIPC and LIPG with lipoprotein subclass traits. 30 established lipid-associated loci had a stronger association for a subclass trait than any conventional trait. These novel association signals provide further insight into the molecular basis of dyslipidemia and the etiology of metabolic disorders.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Alcohol is a consistently identified risk factor for colon cancer. However, the molecular mechanism underlying its effect on normal colon crypt cells remains poorly understood. We employed ...RNA-sequencing to asses transcriptomic response to ethanol exposure (0.2% vol:vol) in 3D organoid lines derived from healthy colon (n = 34). Paired regression analysis identified 2,162 differentially expressed genes in response to ethanol. When stratified by colon location, a far greater number of differentially expressed genes were identified in organoids derived from the left versus right colon, many of which corresponded to cell-type specific markers. To test the hypothesis that the effects of ethanol treatment on colon organoid populations were in part due to differential cell composition, we incorporated external single cell RNA-sequencing data from normal colon biopsies to estimate cellular proportions following single cell deconvolution. We inferred cell-type-specific changes, and observed an increase in transit amplifying cells following ethanol exposure that was greater in organoids from the left than right colon, with a concomitant decrease in more differentiated cells. If this occurs in the colon following alcohol consumption, this would lead to an increased zone of cells in the lower crypt where conditions are optimal for cell division and the potential to develop mutations.
We used an unbiased genome-wide approach to identify exonic variants segregating with diabetes in a multigenerational Finnish family. At least eight members of this family presented with diabetes ...with age of diagnosis ranging from 18 to 51 years and a pattern suggesting autosomal dominant inheritance. We sequenced the exomes of four affected members of this family and performed follow-up genotyping of additional affected and unaffected family members. We uncovered a novel nonsynonymous variant (p.Trp314Arg) in the Wolfram syndrome 1 (WFS1) gene that segregates completely with the diabetic phenotype. Multipoint parametric linkage analysis with 13 members of this family identified a single linkage signal with maximum logarithm of odds score 3.01 at 4p16.2-p16.1, corresponding to a region harboring the WFS1 locus. Functional studies demonstrate a role for this variant in endoplasmic reticulum stress, which is consistent with the β-cell failure phenotype seen in mutation carriers. This represents the first compelling report of a mutation in WFS1 associated with dominantly inherited nonsyndromic adult-onset diabetes.
Risk variants identified so far for colorectal cancer explain only a small proportion of familial risk of this cancer, particularly in Asians.
We performed a genome-wide association study (GWAS) of ...colorectal cancer in East Asians, including 23,572 colorectal cancer cases and 48,700 controls. To identify novel risk loci, we selected 60 promising risk variants for replication using data from 58,131 colorectal cancer cases and 67,347 controls of European descent. To identify additional risk variants in known colorectal cancer loci, we performed conditional analyses in East Asians.
An indel variant, rs67052019 at 1p13.3, was found to be associated with colorectal cancer risk at
× 10
in Asians (
per allele deletion = 1.13, 95% confidence interval = 1.08-1.18). This association was replicated in European descendants using a variant (rs2938616) in complete linkage disequilibrium with rs67052019 (
= 7.7 × 10
). Of the remaining 59 variants, 12 showed an association at
< 0.05 in the European-ancestry study, including rs11108175 and rs9634162 at
< 5 × 10
and two variants with an association near the genome-wide significance level (rs60911071,
= 5.8 × 10
; rs62558833,
= 7.5 × 10
) in the combined analyses of Asian- and European-ancestry data. In addition, using data from East Asians, we identified 13 new risk variants at 11 loci reported from previous GWAS.
In this large GWAS, we identified three novel risk loci and two highly suggestive loci for colorectal cancer risk and provided evidence for potential roles of multiple genes and pathways in the etiology of colorectal cancer. In addition, we showed that additional risk variants exist in many colorectal cancer risk loci identified previously.
Our study provides novel data to improve the understanding of the genetic basis for colorectal cancer risk.
Disturbances in one-carbon metabolism, intracellular reactions involved in nucleotide synthesis and methylation, likely increase the risk of colorectal cancer (CRC). However, results have been ...inconsistent. To explore whether this inconsistency could be explained by intertumoral heterogeneity, we evaluated a comprehensive panel of one-carbon metabolism biomarkers and some single nucleotide polymorphisms (SNPs) in relation to the risk of molecular subtypes of CRC defined by mutations in the KRAS and BRAF oncogenes. This nested case-control study included 488 CRC cases and 947 matched controls from two population-based cohorts in the Northern Sweden Health and Disease Study. We analyzed 14 biomarkers and 17 SNPs in prediagnostic blood and determined KRAS and BRAF mutation status in tumor tissue. In a multivariate network analysis, no variable displayed a strong association with the risk of specific CRC subtypes. A non-synonymous SNP in the CTH gene, rs1021737, had a stronger association compared with other variables. In subsequent univariate analyses, participants with variant rs1021737 genotype had a decreased risk of KRAS-mutated CRC (OR per allele = 0.72, 95% CI = 0.50, 1.05), and an increased risk of BRAF-mutated CRC (OR per allele = 1.56, 95% CI = 1.07, 2.30), with weak evidence for heterogeneity (Pheterogeneity = 0.01). This subtype-specific SNP association was not replicated in a case-case analysis of 533 CRC cases from The Cancer Genome Atlas (P = 0.85). In conclusion, we found no support for clear subtype-specific roles of one-carbon metabolism biomarkers and SNPs in CRC development, making differences in CRC molecular subtype distributions an unlikely explanation for the varying results on the role of one-carbon metabolism in CRC development across previous studies. Further investigation of the CTH gene in colorectal carcinogenesis with regards to KRAS and BRAF mutations or other molecular characteristics of the tumor may be warranted.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Genome-wide association studies (GWASs) have identified more than 150 loci associated with blood lipid and cholesterol levels; however, the functional and molecular mechanisms for many associations ...are unknown. We examined the functional regulatory effects of candidate variants at the GALNT2 locus associated with high-density lipoprotein cholesterol (HDL-C). Fine-mapping and conditional analyses in the METSIM study identified a single locus harboring 25 noncoding variants (r2 > 0.7 with the lead GWAS variants) strongly associated with total cholesterol in medium-sized HDL (e.g., rs17315646, p = 3.5 × 10−12). We used luciferase reporter assays in HepG2 cells to test all 25 variants for allelic differences in regulatory enhancer activity. rs2281721 showed allelic differences in transcriptional activity (75-fold T versus 27-fold C more than the empty-vector control), as did a separate 780-bp segment containing rs4846913, rs2144300, and rs6143660 (49-fold AT– haplotype versus 16-fold CC+ haplotype more). Using electrophoretic mobility shift assays, we observed differential CEBPB binding to rs4846913, and we confirmed this binding in a native chromatin context by performing chromatin-immunoprecipitation (ChIP) assays in HepG2 and Huh-7 cell lines of differing genotypes. Additionally, sequence reads in HepG2 DNase-I-hypersensitivity and CEBPB ChIP-seq signals spanning rs4846913 showed significant allelic imbalance. Allelic-expression-imbalance assays performed with RNA from primary human hepatocyte samples and expression-quantitative-trait-locus (eQTL) data in human subcutaneous adipose tissue samples confirmed that alleles associated with increased HDL-C are associated with a modest increase in GALNT2 expression. Together, these data suggest that at least rs4846913 and rs2281721 play key roles in influencing GALNT2 expression at this HDL-C locus.