Gallbladder Cancer Risk and Indigenous South American Mapuche Ancestry: Instrumental Variable Analysis Using Ancestry-Informative Markers.

A strong association between the proportion of indigenous South American Mapuche ancestry and the risk of gallbladder cancer (GBC) has been reported in observational studies. Chileans show the highest incidence of GBC worldwide, and the Mapuche are the largest indigenous people in Chile. We set out to assess the confounding-free effect of the individual proportion of Mapuche ancestry on GBC risk and to investigate the mediating effects of gallstone disease and body mass index (BMI) on this association. Genetic markers of Mapuche ancestry were selected based on the informativeness for assignment measure, and then used as instrumental variables in two-sample Mendelian randomization analyses and complementary sensitivity analyses. Results suggested a putatively causal effect of Mapuche ancestry on GBC risk (inverse variance-weighted (IVW) risk increase of 0.8% per 1% increase in Mapuche ancestry proportion, 95% CI 0.4% to 1.2%, p = 6.7 × 10-5) and also on gallstone disease (3.6% IVW risk increase, 95% CI 3.1% to 4.0%), pointing to a mediating effect of gallstones on the association between Mapuche ancestry and GBC. In contrast, the proportion of Mapuche ancestry showed a negative effect on BMI (IVW estimate -0.006 kg/m2, 95% CI -0.009 to -0.003). The results presented here may have significant implications for GBC prevention and are important for future admixture mapping studies. Given that the association between the individual proportion of Mapuche ancestry and GBC risk previously noted in observational studies appears to be free of confounding, primary and secondary prevention strategies that consider genetic ancestry could be particularly efficient.

Development and internal validation of a multifactorial risk prediction model for gallbladder cancer in a high-incidence country.

Since 2006, Chile has been implementing a gallbladder cancer (GBC) prevention program based on prophylactic cholecystectomy for gallstone patients aged 35 to 49 years. The effectiveness of this prevention program has not yet been comprehensively evaluated. We conducted a retrospective study of 473 Chilean GBC patients and 2137 population-based controls to develop and internally validate three GBC risk prediction models. The Baseline Model accounted for gallstones while adjusting for sex and birth year. Enhanced Model I also included the non-genetic risk factors: body mass index, educational level, Mapuche surnames, number of children and family history of GBC. Enhanced Model II further included Mapuche ancestry and the genotype for rs17209837. Multiple Cox regression was applied to assess the predictive performance, quantified by the area under the precision-recall curve (AUC-PRC) and the number of cholecystectomies needed (NCN) to prevent one case of GBC at age 70 years. The AUC-PRC for the Baseline Model (0.44%, 95%CI 0.42-0.46) increased by 0.22 (95%CI 0.15-0.29) when non-genetic factors were included, and by 0.25 (95%CI 0.20-0.30) when incorporating non-genetic and genetic factors. The overall NCN for Chileans with gallstones (115, 95%CI 104-131) decreased to 92 (95%CI 60-128) for Chileans with a higher risk than the median according to Enhanced Model I, and to 80 (95%CI 59-110) according to Enhanced Model II. In conclusion, age, sex and gallstones are strong risk factors for GBC, but consideration of other non-genetic factors and individual genotype data improves risk prediction and may optimize allocation of financial resources and surgical capacity.

Gallstones, Cholecystectomy, and Kidney Cancer: Observational and Mendelian Randomization Results Based on Large Cohorts.

BACKGROUND & AIMS: Gallstones (cholelithiasis) constitute a major health burden with high costs related to surgical removal of the gallbladder (cholecystectomy), generally indicated for symptomatic gallstones. The association between gallstones and cholecystectomy and kidney cancer is controversial. We comprehensively investigated this association, considering age at cholecystectomy and time from cholecystectomy to kidney cancer diagnosis, and assessed the causal effect of gallstones on kidney cancer risk by Mendelian randomization (MR). METHODS: We compared the risk of kidney cancer in cholecystectomized and noncholecystectomized patients (16.6 million in total) from the Swedish nationwide cancer, census, patient, and death registries using hazard ratios (HRs). For 2-sample and multivariable MR, we used summary statistics based on 408,567 UK Biobank participants. RESULTS: During a median follow-up of 13 years, 2627 of 627,870 cholecystectomized Swedish patients developed kidney cancer (HR, 1.17; 95% CI, 1.12-1.22). Kidney cancer risk was particularly increased in the first 6 months after cholecystectomy (HR, 3.79; 95% CI, 3.18-4.52) and in patients cholecystectomized before age 40 years (HR, 1.55; 95% CI, 1.39-1.72). MR results based on 18,417 patients with gallstones and 1788 patients with kidney cancer from the United Kingdom revealed a causal effect of gallstones on kidney cancer risk (9.6% risk increase per doubling in gallstone prevalence; 95% CI, 1.2%-18.8%). CONCLUSIONS: Both observational and causal MR estimates based on large prospective cohorts support an increased risk of kidney cancer in patients with gallstones. Our findings provide solid evidence for the compelling need to diagnostically rule out kidney cancer before and during gallbladder removal, to prioritize kidney cancer screening in patients undergoing cholecystectomy in their 30s, and to investigate the underlying mechanisms linking gallstones and kidney cancer in future studies.

Identification of Circulating lncRNAs Associated with Gallbladder Cancer Risk by Tissue-Based Preselection, Cis-eQTL Validation, and Analysis of Association with Genotype-Based Expression.

Long noncoding RNAs (lncRNAs) play key roles in cell processes and are good candidates for cancer risk prediction. Few studies have investigated the association between individual genotypes and lncRNA expression. Here we integrate three separate datasets with information on lncRNA expression only, both lncRNA expression and genotype, and genotype information only to identify circulating lncRNAs associated with the risk of gallbladder cancer (GBC) using robust linear and logistic regression techniques. In the first dataset, we preselect lncRNAs based on expression changes along the sequence "gallstones → dysplasia → GBC". In the second dataset, we validate associations between genetic variants and serum expression levels of the preselected lncRNAs (cis-lncRNA-eQTLs) and build lncRNA expression prediction models. In the third dataset, we predict serum lncRNA expression based on individual genotypes and assess the association between genotype-based expression and GBC risk. AC084082.3 and LINC00662 showed increasing expression levels (p-value = 0.009), while C22orf34 expression decreased in the sequence from gallstones to GBC (p-value = 0.04). We identified and validated two cis-LINC00662-eQTLs (r2 = 0.26) and three cis-C22orf34-eQTLs (r2 = 0.24). Only LINC00662 showed a genotyped-based serum expression associated with GBC risk (OR = 1.25 per log2 expression unit, 95% CI 1.04-1.52, p-value = 0.02). Our results suggest that preselection of lncRNAs based on tissue samples and exploitation of cis-lncRNA-eQTLs may facilitate the identification of circulating noncoding RNAs linked to cancer risk.

A New Pipeline for the Normalization and Pooling of Metabolomics Data.

Pooling metabolomics data across studies is often desirable to increase the statistical power of the analysis. However, this can raise methodological challenges as several preanalytical and analytical factors could introduce differences in measured concentrations and variability between datasets. Specifically, different studies may use variable sample types (e.g., serum versus plasma) collected, treated, and stored according to different protocols, and assayed in different laboratories using different instruments. To address these issues, a new pipeline was developed to normalize and pool metabolomics data through a set of sequential steps: (i) exclusions of the least informative observations and metabolites and removal of outliers; imputation of missing data; (ii) identification of the main sources of variability through principal component partial R-square (PC-PR2) analysis; (iii) application of linear mixed models to remove unwanted variability, including samples' originating study and batch, and preserve biological variations while accounting for potential differences in the residual variances across studies. This pipeline was applied to targeted metabolomics data acquired using Biocrates AbsoluteIDQ kits in eight case-control studies nested within the European Prospective Investigation into Cancer and Nutrition (EPIC) cohort. Comprehensive examination of metabolomics measurements indicated that the pipeline improved the comparability of data across the studies. Our pipeline can be adapted to normalize other molecular data, including biomarkers as well as proteomics data, and could be used for pooling molecular datasets, for example in international consortia, to limit biases introduced by inter-study variability. This versatility of the pipeline makes our work of potential interest to molecular epidemiologists.

Gallstones, Body Mass Index, C-Reactive Protein, and Gallbladder Cancer: Mendelian Randomization Analysis of Chilean and European Genotype Data.

BACKGROUND AND AIMS: Gallbladder cancer (GBC) is a neglected disease with substantial geographical variability: Chile shows the highest incidence worldwide, while GBC is relatively rare in Europe. Here, we investigate the causal effects of risk factors considered in current GBC prevention programs as well as C-reactive protein (CRP) level as a marker of chronic inflammation. APPROACH AND RESULTS: We applied two-sample Mendelian randomization (MR) using publicly available data and our own data from a retrospective Chilean and a prospective European study. Causality was assessed by inverse variance weighted (IVW), MR-Egger regression, and weighted median estimates complemented with sensitivity analyses on potential heterogeneity and pleiotropy, two-step MR, and mediation analysis. We found evidence for a causal effect of gallstone disease on GBC risk in Chileans (P = 9 × 10-5 ) and Europeans (P = 9 × 10-5 ). A genetically elevated body mass index (BMI) increased GBC risk in Chileans (P = 0.03), while higher CRP concentrations increased GBC risk in Europeans (P = 4.1 × 10-6 ). European results suggest causal effects of BMI on gallstone disease (P = 0.008); public Chilean data were not, however, available to enable assessment of the mediation effects among causal GBC risk factors. CONCLUSIONS: Two risk factors considered in the current Chilean program for GBC prevention are causally linked to GBC risk: gallstones and BMI. For Europeans, BMI showed a causal effect on gallstone risk, which was itself causally linked to GBC risk.

Robust Huber-LASSO for improved prediction of protein, metabolite and gene expression levels relying on individual genotype data.

Least absolute shrinkage and selection operator (LASSO) regression is often applied to select the most promising set of single nucleotide polymorphisms (SNPs) associated with a molecular phenotype of interest. While the penalization parameter λ restricts the number of selected SNPs and the potential model overfitting, the least-squares loss function of standard LASSO regression translates into a strong dependence of statistical results on a small number of individuals with phenotypes or genotypes divergent from the majority of the study population-typically comprised of outliers and high-leverage observations. Robust methods have been developed to constrain the influence of divergent observations and generate statistical results that apply to the bulk of study data, but they have rarely been applied to genetic association studies. In this article, we review, for newcomers to the field of robust statistics, a novel version of standard LASSO that utilizes the Huber loss function. We conduct comprehensive simulations and analyze real protein, metabolite, mRNA expression and genotype data to compare the stability of penalization, the cross-iteration concordance of the model, the false-positive and true-positive rates and the prediction accuracy of standard and robust Huber-LASSO. Although the two methods showed controlled false-positive rates ≤2.1% and similar true-positive rates, robust Huber-LASSO outperformed standard LASSO in the accuracy of predicted protein, metabolite and gene expression levels using individual SNP data. The conducted simulations and real-data analyses show that robust Huber-LASSO represents a valuable alternative to standard LASSO in genetic studies of molecular phenotypes.

Epigenome-Wide Analysis of Methylation Changes in the Sequence of Gallstone Disease, Dysplasia, and Gallbladder Cancer.

BACKGROUND AND AIMS: Gallbladder cancer (GBC) is a highly aggressive malignancy of the biliary tract. Most cases of GBC are diagnosed in low-income and middle-income countries, and research into this disease has long been limited. In this study we therefore investigate the epigenetic changes along the model of GBC carcinogenesis represented by the sequence gallstone disease → dysplasia → GBC in Chile, the country with the highest incidence of GBC worldwide. APPROACH AND RESULTS: To perform epigenome-wide methylation profiling, genomic DNA extracted from sections of formalin-fixed, paraffin-embedded gallbladder tissue was analyzed using Illumina Infinium MethylationEPIC BeadChips. Preprocessed, quality-controlled data from 82 samples (gallstones n = 32, low-grade dysplasia n = 13, high-grade dysplasia n = 9, GBC n = 28) were available to identify differentially methylated markers, regions, and pathways as well as changes in copy number variations (CNVs). The number and magnitude of epigenetic changes increased with disease development and predominantly involved the hypermethylation of cytosine-guanine dinucleotide islands and gene promoter regions. The methylation of genes implicated in Wnt signaling, Hedgehog signaling, and tumor suppression increased with tumor grade. CNVs also increased with GBC development and affected cyclin-dependent kinase inhibitor 2A, MDM2 proto-oncogene, tumor protein P53, and cyclin D1 genes. Gains in the targetable Erb-B2 receptor tyrosine kinase 2 gene were detected in 14% of GBC samples. CONCLUSIONS: Our results indicate that GBC carcinogenesis comprises three main methylation stages: early (gallstone disease and low-grade dysplasia), intermediate (high-grade dysplasia), and late (GBC). The identified gradual changes in methylation and CNVs may help to enhance our understanding of the mechanisms underlying this aggressive disease and eventually lead to improved treatment and early diagnosis of GBC.

Genotype-Based Gene Expression in Colon Tissue-Prediction Accuracy and Relationship with the Prognosis of Colorectal Cancer Patients.

Colorectal cancer (CRC) survival has environmental and inherited components. The expression of specific genes can be inferred based on individual genotypes-so called expression quantitative trait loci. In this study, we used the PrediXcan method to predict gene expression in normal colon tissue using individual genotype data from 91 CRC patients and examined the correlation ρ between predicted and measured gene expression levels. Out of 5434 predicted genes, 58% showed a negative ρ value and only 16% presented a ρ higher than 0.10. We subsequently investigated the association between genotype-based gene expression in colon tissue for genes with ρ > 0.10 and survival of 4436 CRC patients. We identified an inverse association between the predicted expression of ARID3B and CRC-specific survival for patients with a body mass index greater than or equal to 30 kg/m2 (HR (hazard ratio) = 0.66 for an expression higher vs. lower than the median, p = 0.005). This association was validated using genotype and clinical data from the UK Biobank (HR = 0.74, p = 0.04). In addition to the identification of ARID3B expression in normal colon tissue as a candidate prognostic biomarker for obese CRC patients, our study illustrates the challenges of genotype-based prediction of gene expression, and the advantage of reassessing the prediction accuracy in a subset of the study population using measured gene expression data.

Genetic Variants in the Regulatory T cell-Related Pathway and Colorectal Cancer Prognosis.

BACKGROUND: High numbers of lymphocytes in tumor tissue, including T regulatory cells (Treg), have been associated with better colorectal cancer survival. Tregs, a subset of CD4+ T lymphocytes, are mediators of immunosuppression in cancer, and therefore variants in genes related to Treg differentiation and function could be associated with colorectal cancer prognosis. METHODS: In a prospective German cohort of 3,593 colorectal cancer patients, we assessed the association of 771 single-nucleotide polymorphisms (SNP) in 58 Treg-related genes with overall and colorectal cancer-specific survival using Cox regression models. Effect modification by microsatellite instability (MSI) status was also investigated because tumors with MSI show greater lymphocytic infiltration and have been associated with better prognosis. Replication of significant results was attempted in 2,047 colorectal cancer patients of the International Survival Analysis in Colorectal Cancer Consortium (ISACC). RESULTS: A significant association of the TGFBR3 SNP rs7524066 with more favorable colorectal cancer-specific survival [hazard ratio (HR) per minor allele: 0.83; 95% confidence interval (CI), 0.74-0.94; P value: 0.0033] was replicated in ISACC (HR: 0.82; 95% CI, 0.68-0.98; P value: 0.03). Suggestive evidence for association was found with two IL7 SNPs, rs16906568 and rs7845577. Thirteen SNPs with differential associations with overall survival according to MSI in the discovery analysis were not confirmed. CONCLUSIONS: Common genetic variation in the Treg pathway implicating genes such as TGFBR3 and IL7 was shown to be associated with prognosis of colorectal cancer patients. IMPACT: The implicated genes warrant further investigation.

Adipose Tissue Properties in Tumor-Bearing Breasts.

The tissue stroma plays a major role in tumors' natural history. Most programs for tumor progression are not activated as cell-autonomous processes but under the conditions of cross-talks between tumor and stroma. Adipose tissue is a major component of breast stroma. This study compares adipose tissues in tumor-bearing breasts to those in tumor-free breasts with the intention of defining a signature that could translate into markers of cancer risk. In tumor-bearing breasts, we sampled adipose tissues adjacent to, or distant from the tumor. Parameters studied included: adipocytes size and density, immune cell infiltration, vascularization, secretome and gene expression. Adipose tissues from tumor-bearing breasts, whether adjacent to or distant from the tumor, do not differ from each other by any of these parameters. By contrast, adipose tissues from tumor-bearing breasts have the capacity to secrete twice as much interleukin 8 (IL-8) than those from tumor-free breasts and differentially express a set of 137 genes of which a significant fraction belongs to inflammation, integrin and wnt signaling pathways. These observations show that adipose tissues from tumor-bearing breasts have a distinct physiological status from those from tumor-free breasts. We propose that this constitutive status contributes as a non-cell autonomous process to determine permissiveness for tumor growth.

RNA Sequencing of Hepatobiliary Cancer Cell Lines: Data and Applications to Mutational and Transcriptomic Profiling.

Cancer cell lines allow the identification of clinically relevant alterations and the prediction of drug response. However, sequencing data for hepatobiliary cancer cell lines in general, and particularly gallbladder cancer (GBC), are sparse. Here, we apply RNA sequencing to characterize 10 GBC, eight hepatocellular carcinoma, and five cholangiocarcinoma (CCA) cell lines. RNA extraction, quality control, library preparation, sequencing, and pre-processing of sequencing data were implemented using state-of-the-art techniques. Public data from the MSK-IMPACT database and a large cohort of Japanese biliary tract cancer patients were used to illustrate the usage of the released data. The total number of exonic mutations varied from 7207 for the cell line NOZ to 9760 for HuCCT1. Researchers planning experiments that require TP53 mutations could use the cell lines NOZ, OCUG-1, SNU308, or YoMi. Mz-Cha-1 showed mutations in ATM, SNU308 presented SMAD4 mutations, and the only investigated cell line that showed ARID1A mutations was GB-d1. SNU478 was the cell line with the global gene expression pattern most similar to GBC, intrahepatic CCA, and extrahepatic CCA. EGFR, KMT2D, and KMT2C generally presented a higher expression in the investigated cell lines than in Japanese primary GBC tumors. We provide the scientific community with detailed mutation and gene expression data, together with three showcase applications, with the aim of facilitating the design of future in vitro cell culture assays for research on hepatobiliary cancer.

Polymorphisms in the Angiogenesis-Related Genes EFNB2, MMP2 and JAG1 Are Associated with Survival of Colorectal Cancer Patients.

An individual's inherited genetic variation may contribute to the 'angiogenic switch', which is essential for blood supply and tumor growth of microscopic and macroscopic tumors. Polymorphisms in angiogenesis-related genes potentially predispose to colorectal cancer (CRC) or affect the survival of CRC patients. We investigated the association of 392 single nucleotide polymorphisms (SNPs) in 33 angiogenesis-related genes with CRC risk and survival of CRC patients in 1754 CRC cases and 1781 healthy controls within DACHS (Darmkrebs: Chancen der Verhütung durch Screening), a German population-based case-control study. Odds ratios and 95% confidence intervals (CI) were estimated from unconditional logistic regression to test for genetic associations with CRC risk. The Cox proportional hazard model was used to estimate hazard ratios (HR) and 95% CIs for survival. Multiple testing was adjusted for by a false discovery rate. No variant was associated with CRC risk. Variants in EFNB2, MMP2 and JAG1 were significantly associated with overall survival. The association of the EFNB2 tagging SNP rs9520090 (p < 0.0001) was confirmed in two validation datasets (p-values: 0.01 and 0.05). The associations of the tagging SNPs rs6040062 in JAG1 (p-value 0.0003) and rs2241145 in MMP2 (p-value 0.0005) showed the same direction of association with overall survival in the first and second validation sets, respectively, although they did not reach significance (p-values: 0.09 and 0.25, respectively). EFNB2, MMP2 and JAG1 are known for their functional role in angiogenesis and the present study points to novel evidence for the impact of angiogenesis-related genetic variants on the CRC outcome.

Multi-omics Analysis Reveals Adipose-tumor Crosstalk in Patients with Colorectal Cancer.

Obesity and obesity-driven cancer rates are continuing to rise worldwide. We hypothesize that adipocyte-colonocyte interactions are a key driver of obesity-associated cancers. To understand the clinical relevance of visceral adipose tissue in advancing tumor growth, we analyzed paired tumor-adjacent visceral adipose, normal mucosa, and colorectal tumor tissues as well as presurgery blood samples from patients with sporadic colorectal cancer. We report that high peroxisome proliferator-activated receptor gamma (PPARG) visceral adipose tissue expression is associated with glycoprotein VI (GPVI) signaling-the major signaling receptor for collagen-as well as fibrosis and adipogenesis pathway signaling in colorectal tumors. These associations were supported by correlations between PPARG visceral adipose tissue expression and circulating levels of plasma 4-hydroxyproline and serum intercellular adhesion molecule 1 (ICAM1), as well as gene set enrichment analysis and joint gene-metabolite pathway results integration that yielded significant enrichment of genes defining epithelial-to-mesenchymal transition-as in fibrosis and metastasis-and genes involved in glycolytic metabolism, confirmed this association. We also reveal that elevated prostaglandin-endoperoxide synthase 2 (PTGS2) colorectal tumor expression is associated with a fibrotic signature in adipose-tumor crosstalk via GPVI signaling and dendritic cell maturation in visceral adipose tissue. Systemic metabolite and biomarker profiling confirmed that high PTGS2 expression in colorectal tumors is significantly associated with higher concentrations of serum amyloid A and glycine, and lower concentrations of sphingomyelin, in patients with colorectal cancer. This multi-omics study suggests that adipose-tumor crosstalk in patients with colorectal cancer is a critical microenvironment interaction that could be therapeutically targeted.See related spotlight by Colacino et al., p. 803.

ABCB1/4 gallbladder cancer risk variants identified in India also show strong effects in Chileans.

BACKGROUND: The first large-scale genome-wide association study of gallbladder cancer (GBC) recently identified and validated three susceptibility variants in the ABCB1 and ABCB4 genes for individuals of Indian descent. We investigated whether these variants were also associated with GBC risk in Chileans, who show the highest incidence of GBC worldwide, and in Europeans with a low GBC incidence. METHODS: This population-based study analysed genotype data from retrospective Chilean case-control (255 cases, 2042 controls) and prospective European cohort (108 cases, 181 controls) samples consistently with the original publication. RESULTS: Our results confirmed the reported associations for Chileans with similar risk effects. Particularly strong associations (per-allele odds ratios close to 2) were observed for Chileans with high Native American (=Mapuche) ancestry. No associations were noticed for Europeans, but the statistical power was low. CONCLUSION: Taking full advantage of genetic and ethnic differences in GBC risk may improve the efficiency of current prevention programs.

Low Back Pain in the Emergency Department: Prevalence of Serious Spinal Pathologies and Diagnostic Accuracy of Red Flags.

BACKGROUND: Very little evidence is available on the prevalence of serious spinal pathologies and the diagnostic accuracy of red flags in patients presenting to the emergency department (ED). This systematic review aims to investigate the prevalence of serious spinal pathologies and the diagnostic accuracy of red flags in patients presenting with low back pain to the ED. METHODS: We systematically searched MEDLINE, PUBMED, EMBASE, Cochrane Library, and SCOPUS from inception to January 2019. Two reviewers independently reviewed the references and evaluated methodological quality. RESULTS: We analyzed 22 studies with a total of 41,320 patients. The prevalence of any requiring immediate/urgent treatment was 2.5%-5.1% in prospective and 0.7%-7.4% in retrospective studies (0.0%-7.2% for vertebral fractures, 0.0%-2.1% for spinal cancer, 0.0%-1.9% for infectious disorders, 0.1%-1.9% for pathologies with spinal cord/cauda equina compression, 0.0%-0.9% for vascular pathologies). Examples of red flags which increased the likelihood for a serious condition were suspicion or history of cancer (spinal cancer); intravenous drug use, indwelling vascular catheter, and other infection site (epidural abscess). CONCLUSION: We found a higher prevalence of serious spinal pathologies in the ED compared to the reported prevalence in primary care settings. As the diagnostic accuracy of most red flags was reported only by a single study, further validation in high-quality prospective studies is needed.

Optimal selection of genetic variants for adjustment of population stratification in European association studies.

Population stratification is usually corrected relying on principal component analysis (PCA) of genome-wide genotype data, even in populations considered genetically homogeneous, such as Europeans. The need to genotype only a small number of genetic variants that show large differences in allele frequency among subpopulations-so-called ancestry-informative markers (AIMs)-instead of the whole genome for stratification adjustment could represent an advantage for replication studies and candidate gene/pathway studies. Here we compare the correction performance of classical and robust principal components (PCs) with the use of AIMs selected according to four different methods: the informativeness for assignment measure ($IN$-AIMs), the combination of PCA and F-statistics, PCA-correlated measurement and the PCA weighted loadings for each genetic variant. We used real genotype data from the Population Reference Sample and The Cancer Genome Atlas to simulate European genetic association studies and to quantify type I error rate and statistical power in different case-control settings. In studies with the same numbers of cases and controls per country and control-to-case ratios reflecting actual rates of disease prevalence, no adjustment for population stratification was required. The unnecessary inclusion of the country of origin, PCs or AIMs as covariates in the regression models translated into increasing type I error rates. In studies with cases and controls from separate countries, no investigated method was able to adequately correct for population stratification. The first classical and the first two robust PCs achieved the lowest (although inflated) type I error, followed at some distance by the first eight $IN$-AIMs.

Expression Patterns of Xenobiotic-Metabolizing Enzymes in Tumor and Adjacent Normal Mucosa Tissues among Patients with Colorectal Cancer: The ColoCare Study.

BACKGROUND: Xenobiotic-metabolizing enzymes (XME) play a critical role in the activation and detoxification of several carcinogens. However, the role of XMEs in colorectal carcinogenesis is unclear. METHODS: We investigated the expression of XMEs in human colorectal tissues among patients with stage I-IV colorectal cancer (n = 71) from the ColoCare Study. Transcriptomic profiling using paired colorectal tumor and adjacent normal mucosa tissues of XMEs (GSTM1, GSTA1, UGT1A8, UGT1A10, CYP3A4, CYP2C9, GSTP1, and CYP2W1) by RNA microarray was compared using Wilcoxon rank-sum tests. We assessed associations between clinicopathologic, dietary, and lifestyle factors and XME expression with linear regression models. RESULTS: GSTM1, GSTA1, UGT1A8, UGT1A10, and CYP3A4 were all statistically significantly downregulated in colorectal tumor relative to normal mucosa tissues (all P ≤ 0.03). Women had significantly higher expression of GSTM1 in normal tissues compared with men (ß = 0.37, P = 0.02). By tumor site, CYP2C9 expression was lower in normal mucosa among patients with rectal cancer versus colon cancer cases (ß = -0.21, P = 0.0005). Smokers demonstrated higher CYP2C9 expression levels in normal mucosa (ß = 0.17, P = 0.02) when compared with nonsmokers. Individuals who used NSAIDs had higher GSTP1 tumor expression compared with non-NSAID users (ß = 0.17, P = 0.03). Higher consumption of cooked vegetables (>1×/week) was associated with higher CYP3A4 expression in colorectal tumor tissues (ß = 0.14, P = 0.007). CONCLUSIONS: XMEs have lower expression in colorectal tumor relative to normal mucosa tissues and may modify colorectal carcinogenesis via associations with clinicopathologic, lifestyle, and dietary factors. IMPACT: Better understanding into the role of drug-metabolizing enzymes in colorectal cancer may reveal biological differences that contribute to cancer development, as well as treatment response, leading to clinical implications in colorectal cancer prevention and management.