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.

Serotonin type 3 receptor subunit gene polymorphisms associated with psychosomatic symptoms in irritable bowel syndrome: A multicenter retrospective study.

BACKGROUND: Single-nucleotide polymorphisms (SNPs) of the serotonin type 3 receptor subunit (HTR3) genes have been associated with psychosomatic symptoms, but it is not clear whether these associations exist in irritable bowel syndrome (IBS). AIM: To assess the association of HTR3 polymorphisms with depressive, anxiety, and somatization symptoms in individuals with IBS. METHODS: In this retrospective study, 623 participants with IBS were recruited from five specialty centers in Germany, Sweden, the United States, the United Kingdom, and Ireland. Depressive, anxiety, and somatization symptoms and sociodemographic characteristics were collected. Four functional SNPs - HTR3A c.-42C>T, HTR3B c.386A>C, HTR3C c.489C>A, and HTR3E c.*76G>A - were genotyped and analyzed using the dominant and recessive models. We also performed separate analyses for sex and IBS subtypes. SNP scores were calculated as the number of minor alleles of the SNPs above. The impact of HTR3C c.489C>A was tested by radioligand-binding and calcium influx assays. RESULTS: Depressive and anxiety symptoms significantly worsened with increasing numbers of minor HTR3C c.489C>A alleles in the dominant model (F depressive = 7.475, P depressive = 0.006; F anxiety = 6.535, P anxiety = 0.011). A higher SNP score (range 0-6) was linked to a worsened depressive symptoms score (F = 7.710, P-linear trend = 0.006) in IBS. The potential relevance of the HTR3C SNP was corroborated, showing changes in the expression level of 5-HT3AC variant receptors. CONCLUSION: We have provided the first evidence that HTR3C c.489C>A is involved in depressive and anxiety symptoms in individuals with IBS. The SNP score indicated that an increasing number of minor alleles is linked to the worsening of depressive symptoms in IBS.

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.

Inbreeding, Native American ancestry and child mortality: linking human selection and paediatric medicine.

The children of related parents show increased risk of early mortality. The Native American genome typically exhibits long stretches of homozygosity, and Latin Americans are highly heterogeneous regarding the individual burden of homozygosity, the proportion and the type of Native American ancestry. We analysed nationwide mortality and genome-wide genotype data from admixed Chileans to investigate the relationship between common causes of child mortality, homozygosity and Native American ancestry. Results from two-stage linear-Poisson regression revealed a strong association between the sum length of runs of homozygosity (SROH) above 1.5 Megabases (Mb) in each genome and mortality due to intracranial non-traumatic haemorrhage of foetus and newborn (5% increased risk of death per Mb in SROH, P = 1 × 10-3) and disorders related to short gestation and low birth weight (P = 3 × 10-4). The major indigenous populations in Chile are Aymara-Quechua in the north of the country and the Mapuche-Huilliche in the south. The individual proportion of Aymara-Quechua ancestry was associated with an increased risk of death due to anencephaly and similar malformations (P = 4 × 10-5), and the risk of death due to Edwards and Patau trisomy syndromes decreased 4% per 1% Aymara-Quechua ancestry proportion (P = 4 × 10-4) and 5% per 1% Mapuche-Huilliche ancestry proportion (P = 2 × 10-3). The present results suggest that short gestation, low birth weight and intracranial non-traumatic haemorrhage mediate the negative effect of inbreeding on human selection. Independent validation of the identified associations between common causes of child death, homozygosity and fine-scale ancestry proportions may inform paediatric medicine.

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.

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.

Amerindian Ancestry Influences Genetic Susceptibility to Chronic Obstructive Pulmonary Disease.

The contribution of genetic ancestry on chronic obstructive pulmonary disease (COPD) predisposition remains unclear. To explore this relationship, we analyzed the associations between 754,159 single nucleotide polymorphisms (SNPs) and risk of COPD (n = 214 cases, 193 healthy controls) in Talca, Chile, considering the genetic ancestry and established risk factors. The proportion of Mapuche ancestry (PMA) was based on a panel of 45 Mapuche reference individuals. Five PRDM15 SNPs and two PPP1R12B SNPs were associate with COPD risk (p = 0.05 to 5×10-4) in those individuals with lower PMA. Based on linkage disequilibrium and sliding window analyses, an adjacent PRDM15 SNPs were associated with COPD risk in the lower PMA group (p = 10-3 to 3.77×10-8). Our study is the first to report an association between PPP1R12B and COPD risk, as well as effect modification between ethnicity and PRDM15 SNPs in determining COPD risk. Our results are biologically plausible given that PPP1R12B and PRDM15 are involved in immune dysfunction and autoimmunity, providing mechanistic evidence for COPD pathogenesis and highlighting the importance to conduct more genome wide association studies (GWAS) in admixed populations with Amerindian descent.

Functional and genomic characterization of three novel cell lines derived from a metastatic gallbladder cancer tumor.

BACKGROUND: Gallbladder cancer (GBC) is the most common tumor of the biliary tract. The incidence of GBC shows a large geographic variability, being particularly frequent in Native American populations. In Chile, GBC represents the second cause of cancer-related death among women. We describe here the establishment of three novel cell lines derived from the ascitic fluid of a Chilean GBC patient, who presented 46% European, 36% Mapuche, 12% Aymara and 6% African ancestry. RESULTS: After immunocytochemical staining of the primary cell culture, we isolated and comprehensively characterized three independent clones (PUC-GBC1, PUC-GBC2 and PUC-GBC3) by short tandem repeat DNA profiling and RNA sequencing as well as karyotype, doubling time, chemosensitivity, in vitro migration capability and in vivo tumorigenicity assay. Primary culture cells showed high expression of CK7, CK19, CA 19-9, MUC1 and MUC16, and negative expression of mesothelial markers. The three isolated clones displayed an epithelial phenotype and an abnormal structure and number of chromosomes. RNA sequencing confirmed the increased expression of cytokeratin and mucin genes, and also of TP53 and ERBB2 with some differences among the three cells lines, and revealed a novel exonic mutation in NF1. The PUC-GBC3 clone was the most aggressive according to histopathological features and the tumorigenic capacity in NSG mice. CONCLUSIONS: The first cell lines established from a Chilean GBC patient represent a new model for studying GBC in patients of Native American descent.

PRDM15 Is Associated with Risk of Chronic Obstructive Pulmonary Disease in a Rural Population in Chile.

BACKGROUND: Genome-wide association studies (GWAS) have accelerated our understanding of the genetic underpinnings of chronic obstructive pulmonary disease (COPD); however, GWAS populations have typically consisted of European descent, with ∼1% of Latin American ancestry. OBJECTIVE: To overcome this limitation, we conducted a GWAS in a rural Chilean population with increased COPD risk to investigate genetic variation of COPD risk in this understudied minority population. METHOD: We carried out a case-control study of 214 COPD patients (defined by the GOLD criteria) and 193 healthy controls in Talca, Chile. DNA was extracted from venous blood and genotyped on the Illumina Global Screening Array (n = 754,159 markers). After exclusion based on Hardy-Weinberg equilibrium (p ≤ 0.001), call rates (<95%), and minor allele frequencies (<0.5%) in controls, 455,564 markers were available for logistic regression. RESULTS: PRDM15 rs1054761 C allele (p = 2.22 × 10-7) was associated with decreased COPD risk. Three PRDM15 SNPs located on chromosome 21 were significantly associated with COPD risk (p < 10-6). Two of these SNPs, rs1054761 and rs4075967, were located on a noncoding transcript variant region of the gene. CONCLUSION: PRDM15 overexpression may play a role in the B-cell dysregulation in COPD pathogenesis. To the best of our knowledge, the association between PRDM15 and COPD risk was not previously found in GWAS studies in largely European populations, highlighting the importance of investigating novel variants associated with COPD risk among ethnically diverse populations.

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.

Functional and genomic characterization of three novel cell lines derived from a metastatic gallbladder cancer tumor

BACKGROUND: Gallbladder cancer (GBC) is the most common tumor of the biliary tract. The incidence of GBC shows a large geographic variability, being particularly frequent in Native American populations. In Chile, GBC represents the second cause of cancer-related death among women. We describe here the establishment of three novel cell lines derived from the ascitic fluid of a Chilean GBC patient, who presented 46% European, 36% Mapuche, 12% Aymara and 6% African ancestry. RESULTS: After immunocytochemical staining of the primary cell culture, we isolated and comprehensively characterized three independent clones (PUC-GBC1, PUC-GBC2 and PUC-GBC3) by short tandem repeat DNA profiling and RNA sequencing as well as karyotype, doubling time, chemosensitivity, in vitro migration capability and in vivo tumorigenicity assay. Primary culture cells showed high expression of CK7, CK19, CA 19-9, MUC1 and MUC16, and negative expression of mesothelial markers. The three isolated clones displayed an epithelial phenotype and an abnormal structure and number of chromosomes. RNA sequencing confirmed the increased expression of cytokeratin and mucin genes, and also of TP53 and ERBB2 with some differences among the three cells lines, and revealed a novel exonic mutation in NF1. The PUC-GBC3 clone was the most aggressive according to histopathological features and the tumorigenic capacity in NSG mice. CONCLUSIONS: The first cell lines established from a Chilean GBC patient represent a new model for studying GBC in patients of Native American descent.

Interaction between genetic ancestry and common breast cancer susceptibility variants in Colombian women.

Latino women show lower incidences of breast cancer (BC) than non-Hispanic whites. Large-scale genetic association studies have identified variants robustly associated with BC risk in European women. We examine here the relevance of these variants to Colombian BC and possible interactions with genetic ancestry. Native American, European and African proportions were estimated for 1022 Colombian BC cases and 1023 controls. Logistic regression was applied to assess the association between 78 variants and BC risk and interactions between the variants and ancestry proportions. We constructed a multifactorial risk score combining established BC risk factors, associated risk variants and individual ancestry proportions. Each 1% increase in the Native American proportion translated into a 2.2% lower BC risk (95% CI: 1.4-2.9). Thirteen variants were associated with BC in Colombian women, with allele frequencies and risk effects partially different from European women. Ancestry proportions moderated the risk effects of two variants. The ability of Native American proportions to separate Colombian cases and controls (area-under-the-curve (AUC) = 0.61) was similar to the discriminative ability of family history of BC in first-degree female relatives (AUC = 0.58) or the combined effect of all 13 associated risk variants (AUC = 0.57). Our findings demonstrate ample potential for individualized BC prevention in Hispanic women taking advantage of individual Native American proportions, information on established susceptibility factors and recently identified common risk variants.

Methods and results from the genome-wide association group at GAW20.

BACKGROUND: This paper summarizes the contributions from the Genome-wide Association Study group (GWAS group) of the GAW20. The GWAS group contributions focused on topics such as association tests, phenotype imputation, and application of empirical kinships. The goals of the GWAS group contributions were varied. A real or a simulated data set based on the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN) study was employed by different methods. Different outcomes and covariates were considered, and quality control procedures varied throughout the contributions. RESULTS: The consideration of heritability and family structure played a major role in some contributions. The inclusion of family information and adaptive weights based on data were found to improve power in genome-wide association studies. It was proven that gene-level approaches are more powerful than single-marker analysis. Other contributions focused on the comparison between pedigree-based kinship and empirical kinship matrices, and investigated similar results in heritability estimation, association mapping, and genomic prediction. A new approach for linkage mapping of triglyceride levels was able to identify a novel linkage signal. CONCLUSIONS: This summary paper reports on promising statistical approaches and findings of the members of the GWAS group applied on real and simulated data which encompass the current topics of epigenetic and pharmacogenomics.

Subtypes of Native American ancestry and leading causes of death: Mapuche ancestry-specific associations with gallbladder cancer risk in Chile.

Latin Americans are highly heterogeneous regarding the type of Native American ancestry. Consideration of specific associations with common diseases may lead to substantial advances in unraveling of disease etiology and disease prevention. Here we investigate possible associations between the type of Native American ancestry and leading causes of death. After an aggregate-data study based on genome-wide genotype data from 1805 admixed Chileans and 639,789 deaths, we validate an identified association with gallbladder cancer relying on individual data from 64 gallbladder cancer patients, with and without a family history, and 170 healthy controls. Native American proportions were markedly underestimated when the two main types of Native American ancestry in Chile, originated from the Mapuche and Aymara indigenous peoples, were combined together. Consideration of the type of Native American ancestry was crucial to identify disease associations. Native American ancestry showed no association with gallbladder cancer mortality (P = 0.26). By contrast, each 1% increase in the Mapuche proportion represented a 3.7% increased mortality risk by gallbladder cancer (95%CI 3.1-4.3%, P = 6×10-27). Individual-data results and extensive sensitivity analyses confirmed the association between Mapuche ancestry and gallbladder cancer. Increasing Mapuche proportions were also associated with an increased mortality due to asthma and, interestingly, with a decreased mortality by diabetes. The mortality due to skin, bladder, larynx, bronchus and lung cancers increased with increasing Aymara proportions. Described methods should be considered in future studies on human population genetics and human health. Complementary individual-based studies are needed to apportion the genetic and non-genetic components of associations identified relying on aggregate-data.