Genome-wide Analysis Identifies Novel Gallstone-susceptibility Loci Including Genes Regulating Gastrointestinal Motility ======================================================================================================================== * Cameron J Fairfield * Thomas M Drake * Riinu Pius * Andrew D Bretherick * Archie Campbell * David W Clark * Jonathan A Fallowfield * Caroline Hayward * Neil C Henderson * Andrii Iakovliev * Peter K Joshi * Nicholas L Mills * David J Porteous * Prakash Ramachandran * Robert K Semple * Catherine A Shaw * Cathie LM Sudlow * Paul RHJ Timmers * James F Wilson * Stephen J Wigmore * Athina Spiliopoulou * Ewen M Harrison ## Abstract **Objective** Genome-wide association studies (GWAS) have identified several risk loci for gallstone disease. As with most polygenic traits, it is likely many genetic determinants are undiscovered. The aim of this study was to identify novel genetic variants that represent new targets for gallstone research and treatment. **Design** We performed a GWAS of 28,627 gallstone cases and 348,373 controls in the UK Biobank and a GWA meta-analysis (43,639 cases and 506,798 controls) with the FinnGen cohort. We assessed pathway enrichment using gene-based then gene-set analysis and tissue expression of implicated genes in Genotype-Tissue Expression project data. We constructed a polygenic risk score (PRS) and evaluated phenotypic traits associated with the score. **Results** Seventy-five risk loci were identified (P<5*10−8) of which forty-six were novel. Pathway enrichment revealed associations with lipid homeostasis, glucuronidation, phospholipid metabolism and gastrointestinal motility. *ANO1* and *TMEM147*, both in novel loci, are strongly expressed in the gallbladder and gastrointestinal tract. Both regulate gastrointestinal motility. The gallstone risk allele rs7599-A leads to suppression of hepatic *TMEM147* expression suggesting the protein protects against gallstone formation. Individuals in the highest decile of the PRS demonstrated a 6-fold increased risk of gallstones compared to the lowest risk category. The PRS was strongly associated with increased body mass index, serum liver enzyme and C-reactive protein concentrations and decreased lipoprotein cholesterol concentrations. **Conclusion** This GWAS demonstrates the polygenic nature of gallstone risk and identifies 46 novel susceptibility loci. For the first time, we implicate genes influencing gastrointestinal motility in the pathogenesis of gallstones. **What is already known on this subject?** * Genome-wide association studies (GWAS) have identified 29 genetic variants within independent loci which increase the risk of gallstone disease. * Most of these variants lie within or near to genes that regulate lipid or bile acid metabolism. Gallstones are known to have a significant genetic component with 25-50% of gallstone disease due to genetic risk factors. * Much of this risk is not accounted for by the known gallstone-susceptibility loci. **What are the new findings?** * We performed a GWAS in the UK Biobank (28,627 gallstone cases, 348,373 controls) and a GWA meta-analysis (43,639 cases and 506,798 controls) with the FinnGen cohort. * We identified a total of 75 gallstone-susceptibility loci with 46 of these being new and the remaining 29 being those already identified. * We annotated the variants based on their position within or near to genes and assessed pathway enrichment through gene-set analysis. * We identify two novel gallstone-susceptibility loci in which the lead variants lie within genes governing gastrointestinal motility which highly expressed in gallbladder (*ANO1* and *TMEM147*). * We demonstrate further loci involved in primary cilia function. * We report significant association of a polygenic risk score with gallstones using independent subsets of the study population. **How might it impact on clinical practice in the foreseeable future?** * The genes and pathways identified represent novel targets for development of medication targeting primary or secondary prevention of gallstones. This may be of particular benefit to those unable to undergo cholecystectomy. * The individual variants or polygenic risk score identified in this GWAS could form the basis for identification of individuals with high risk of gallstones to support screening or treatment of gallstone disease. This screening may be of particular benefit in populations with elevated risk of gallstones such as haemolytic disease or bariatric surgery. ## Introduction Gallstone disease is one of the most common reasons for acute presentation to hospitals worldwide1 and generates the greatest healthcare expenditure of any gastrointestinal condition in the United States2. Prevalence of gallstones reaches 20-40% in high-risk populations3,4 and gallstones are highly heritable5. Identification of novel pathways influencing gallstone development may improve risk stratification in high risk populations and support development of medications targeting gallstone formation and other cholestatic disease for use in primary and secondary prevention. A total of seven genome-wide association studies (GWAS) have been published revealing 29 gallstone-promoting loci6–12. Most gallstones are composed of cholesterol which crystallises in bile at high concentration and/or in bile lacking protective bile acids13. Many of the identified variants relate to cholesterol and bile acid metabolism. Despite evidence of impaired gallbladder and gastrointestinal motility contributing to the formation of gallstones, there have been no gallstone-susceptibility variants identified that influence these characteristics14. Using larger and more richly phenotyped cohorts may identify such variants and contribute to the knowledge of gallstone formation. The most recent GWAS and GWAS meta-analysis used publicly-available summary GWAS data from the Global Biobank Engine (GBE)15 based on UK Biobank (UKB) data16. GBE is based on three-level ICD-10 and ICD-9 codes and considered gallstones to be code “K80” (“cholelithiasis”) or self-reported gallstones at study enrolment. The UKB captures several other informative variables including OPCS codes, other ICD-9 and 10 codes, self-reported procedures relating to gallstones and primary care data. By incorporating these additional data, a greater proportion of patients with gallstone disease can be classified as cases rather than misclassified as controls. The aim of this study was to perform an updated GWAS using more richly phenotyped participants to identify novel pathways contributing to gallstone formation. The genetic associations identified were replicated in the Generation Scotland: Scottish Family Health Study (GS-SFHS) cohort17 and meta-analysed with FinnGen summary association statistics18. A polygenic risk score was derived and analysis for gene-environment interaction undertaken. ## Materials and Methods Methods for identification of participants, determination of gallstone status and genotyping for the UKB, GS-SFHS and the FinnGen cohort are provided in Supplementary Materials 1. Diagnostic codes used for the identification of gallstones are shown in Supplementary Tables 1-3. The UKB received ethical approval (Research Ethics Committee [REC] reference number: 11/NW/0382). UKB data access was approved under projects 30439 (phenotype data) and 19655 (genotype data). Ethical approval for GS-SFHS was granted by NHS Tayside REC (reference number: 05/S1401/89). The FinnGen study protocol was approved by the Ethical Review Board of the Hospital District of Helsinki and Uusimaa (Nr HUS/990/2017). ### Genome-wide Association Study #### Initial Analysis Genotyped and imputed SNPs were analysed. Participants of self-reported European ancestry were included and outliers for heterozygosity and unexpected runs of homozygosity were excluded. One randomly-selected participant from each pair of related individuals in the UKB (Kinship > 0.088419) was excluded. Association with gallstones was analysed using logistic regression adjusted for age, sex, the first 20 genetic principal components and batch, with batch included as a random effect. Imputation dosage was used for imputed SNPs. ![Formula][1] SNPs were excluded based on minor allele frequency (MAF < 0.001), imputation quality (INFO < 0.3) and departure from Hardy-Weinberg equilibrium (HWE: P < 5*10−6). The MAF was selected on the basis of ensuring 25 or more minor alleles in the cases (*MAF* = 25/(2 * *n_cases)*) which corresponded to a MAF of 0.0004. This was rounded up to 0.001. Genome-wide significance was determined as P<5*10−8 in the UKB and Bonferroni-corrected *P**C* = 0.05 in the replication cohorts (0.05 divided by the number of loci tested). Replication was undertaken by analysing the lead SNP within each locus or a linkage disequilibrium (LD) proxy. P values smaller than 4.94 * 10−324 (the limit of numerical precision on a Linux 64-bit operating system) were rounded to 4.94 * 10−324 (Supplementary Materials 1). #### Genome-Wide Association Meta-analysis A genome-wide association meta-analysis was undertaken using UKB data and FinnGen data in an inverse variance-weighted fixed effects model. SNPs shared between the datasets (10,044,759) were analysed. RsID was converted to the format “rs1234:A:T” to avoid ambiguity at multi-allelic SNPs. #### Linkage Disequilibrium Clumping and Conditional Analyses LD clumping was performed using the Functional Mapping and Annotation of Genome-Wide Association Studies (FUMA) web application version 1.3.5d20. Loci were established for each lead SNP with a minimum distance of 250Kb between loci and using an *r*2 < 0.25 to indicate separate independent SNPs within the same locus. Details of the FUMA parameters are available in Supplementary Materials 2A. Five loci with rare lead SNPs not present in the 1000 Genomes Phase 3 reference panel (MAFUKB⍰0.01-0.001) were generated manually by creating a 500 KB window centred on the lead SNP. Each locus was re-analysed whilst conditioning on the lead SNP and further signals with genome-wide significance were identified. This process was repeated with additional lead SNPs until no remaining SNPs reached genome-wide significance. #### Population Stratification The genomic inflation factor *λ**gc* was calculated using the summary GWAS statistics. Evidence of test statistic inflation in *λ**gc* was investigated with linkage disequilibrium score regression21. ### Sensitivity Analysis As the definition of gallstones was based on clinical detection of gallstones it is possible that detection bias may influence the results. Gallstones detected during investigations for other conditions may be recorded despite causing no symptoms. Genetic variants causing other hepatic pathology may therefore influence the likelihood of detection, rather than formation, of gallstones. Individuals requiring surgical management of gallstones are likely to be symptomatic, minimising detection bias. A sensitivity analysis in which cases underwent a cholecystectomy, or any gallstone-specific procedure as defined in Supplementary Table 2, and controls were free from gallstone disease was undertaken at each of the identified lead SNPs. Those with gallstones who did not undergo surgery were excluded. ### Association of Lithogenic Loci with Phenotypic Traits Associations with serum biochemistry values for each identified variant were assessed in an age- and sex-adjusted linear regression model. Serum bilirubin (total, conjugated and unconjugated), alanine aminotransferase (ALT), aspartate aminotransferase (AST), gamma-glutamyl transferase (GGT), alkaline phosphatase (ALP), total cholesterol (TC), low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), triglycerides, apolipoprotein A and B (ApoA and ApoB), lipoprotein A (LipoA), haemoglobin (Hb), red blood cells (RBC), reticulocyte percentage, C-reactive protein (CRP) and Cystatin C were assessed. Each trait was assessed visually via histogram and log-transformed in the case of skewed distribution. ### Gene-Environment Interaction The *UGT1A1* polymorphism rs4148325 associates with bilirubin and causes Gilbert’s syndrome but is a comparatively modest driver of gallstones in GWAS12,22. In patients with haemolytic disease or cystic fibrosis, the variant is a significant risk factor for gallstones23–25. This raises the possibility of a gene-environment interaction in which defects of bilirubin metabolism only lead to gallstone formation in the context of increased bilirubin production or decreased bilirubin excretion. Each identified locus was re-analysed for a gene-environment interaction with serum unconjugated bilirubin as a covariate and as an interaction term with genotype. A Wald test was performed to determine the presence of a gene-environment interaction and significance was determined using a Bonferroni correction. This was not conducted at the known *UGT1A1* locus as unconjugated bilirubin is likely to be a mediator of the relationship between the variant and gallstones. ### Polygenic Risk Score A polygenic risk score (PRS) was derived using a random 80%/20% split of the UKB cohort. A separate GWAS was conducted in the 80% and the summary associations statistics used to derive the PRS which was validated in the 20% split using PRSice26. This was conducted as an additive model with default PRSice parameters on directly genotyped SNPs. The PRS was evaluated for association with the same phenotypic traits as individual risk loci in a simple, unadjusted linear or logistic regression. ### Pathway Enrichment We conducted a gene-based analysis in MAGMA to identify genes and gene-sets showing significant association with gallstones. Methods for gene-based and gene-set analysis have been described previously27. In brief, summary association statistics, following FUMA-based annotation, were grouped by gene and the average χ2 statistic across all SNPs in each gene computed and converted to gene-based P values. Genes demonstrating significant association (after Bonferroni correction) were taken forward for a competitive gene-set analysis in which gene-sets were tested for association with the phenotype. This tests whether the mean association of the phenotype with genes in the gene set is greater than that of genes not in the gene set. The mean association was adjusted for gene size, gene density (based on genetic distance using between-SNP LD), the inverse of the mean minor allele frequency in the gene, and the log value of these three variables. This adjustment is necessary to prevent spurious associations driven by the number and frequency of possible mutations within a given gene rather than a true effect. 4,761 curated gene sets and 5,917 gene ontology (GO) terms were available for analysis from Molecular Signatures Database v6.2. Gene tissue expression analysis was undertaken in MAGMA using the identified genes. Loci without clear relation to gallstones were assessed by literature search based on the FUMA gene annotations for genes nearest to the lead SNP to identify further pathways not present or identified within the curated gene-sets. Lead SNPs were analysed in the Genotype-Tissue expression project (GTEx) to assess their effect on expression of putative causal genes. ### Plotting and Statistical Analysis of Results Genomic analyses were performed using SNPtest (version 2.5.4 for gene-environment interaction model, version 2.5.2 for all other analyses) and QCTOOLS (version 2.0.6)28 within a Linux high-performance compute cluster. GWA meta-analysis was undertaken in METAL29. Post-GWAS analysis, regression analyses and plotting were performed using R version 3.6.330. Methods and results are reported in accordance with the STREGA guidelines31 (Supplementary Materials 3). ### Patient Involvement During the design stage of the project, the NHS Lothian Public and Patient Involvement and Engagement was involved in drafting lay summaries of the proposed research. Published results will be returned to the UK Biobank team who inform participants of research outputs. ## Results ### Gallstones Cohort A total of 502,616 participants entered the UKB study of which 377,000 were taken forward for GWAS. A history of gallstones was present in 28,627 with 348,373 controls. The median follow-up available was 9.0 years (range 7.3-11.9). Baseline demographics of the UKB cohort are shown in Supplementary Table 4. A total of 24,096 participants entered GS-SFHS of which 6,317 unrelated individuals were taken forward for GWAS. A history of gallstones was present in 1,089 with 5,228 controls. The mean follow-up available was 10.7 years (range 9.6-14.7) (Figure 1). ![Figure1](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/07/19/2021.07.16.21260637/F1.medium.gif) [Figure1](http://medrxiv.org/content/early/2021/07/19/2021.07.16.21260637/F1) In FinnGen 173,437 participants were analysed for gallstones (15,012 cases, 158,425 controls). ### Genome-wide Association Study #### Initial UKB Analysis A total of 4,065 SNPs were identified with genome-wide significant P values after removal of low-quality SNPs (see Supplementary Materials 4). An additional 11,539 SNPs demonstrated borderline significance (*P* < 5 * 10−5). LD-clumping revealed 50 loci with at least one significant signal (see Supplementary Materials 2B). Five additional loci contained variants with a significant association with gallstones for which the lead SNP was not present in the 1000 Genomes reference panel (MAF 0.01 to 0.001) and were therefore not annotated by the FUMA application. LD-clumping for those loci was performed using boundaries of 250KB from the lead SNP. Of the 55 loci identified in the UKB, 27 loci are newly identified variants, 7 are known from previous GWAS using other datasets and 21 variants were previously identified only in studies using GBE data10,12 (Table 1). Four variants were replicated in GS-SFHS (rs11887534, rs1993453, rs1260326 and rs62129966; P < 0.00091, corrected for 55 variants). 33 variants from the UKB GWAS reached the replication threshold in the FinnGen cohort (P < 0.00012, corrected for 42 variants available) and 38 reached nominal significance (*P* < 0.05). 10 of the 27 novel loci identified in the UKB GWAS reached the replication threshold in FinnGen. Results for replication in GS-SFHS and FinnGen at each of the 55 loci identified in the UKB are shown in Supplementary Table 5. Loci with replication in either GS-SFHS or FinnGen are highlighted in bold in Table 1. View this table: [Table 1:](http://medrxiv.org/content/early/2021/07/19/2021.07.16.21260637/T1) Table 1: Gallstone-Susceptibility Loci in the UK Biobank GWAS The D19H variant of *ABCG8* (rs11887534) demonstrated a highly significant result (*P* < 4.94 * 10−324). 23 adjacent variants demonstrated identical P values and almost identical *β*-coefficients. As rs11887534 is a known protein-altering variant (D19H) and previously identified as the lead SNP by other GWAS, this was considered as the lead SNP. #### GWA Meta-analysis GWA meta-analysis with FinnGen data revealed a total of 63 loci associated with gallstone disease (*P* < 5 * 10−8). Of these, 43 were identified in the UKB, one was previously identified with GBE data10 and a further 19 were novel loci bringing the total number of novel gallstone-susceptibility loci to 46 and the overall total of gallstone-susceptibility loci to 75. Of the 29 previously reported loci, 28 were identified in the UKB GWAS and all 29 were identified in the meta-analysis. Twelve loci identified in the UKB GWAS did not retain significance on meta-analysis. The 63 lead SNPs from the meta-analysis are shown in Table 2. View this table: [Table 2:](http://medrxiv.org/content/early/2021/07/19/2021.07.16.21260637/T2) Table 2: Gallstone-Susceptibility Loci in GWAS Meta-analysis A combined circular Manhattan plot for association of variants with gallstones in the meta-analysis as well as the 12 loci only identified in the UKB is shown in Figure 2. Rectangular Manhattan plots for the meta-analysis and separately for the initial UKB analysis are available in Supplementary Figures 1 to 4. ![Figure2](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/07/19/2021.07.16.21260637/F2.medium.gif) [Figure2](http://medrxiv.org/content/early/2021/07/19/2021.07.16.21260637/F2) ### Conditional Analyses Conditional analyses, adjusting for the lead SNP, were undertaken for the 55 loci identified in the UKB GWAS. After conditional analyses, a total of 6 loci were found to have additional independent signal suggesting that more than one variant within these loci has a causal relationship with gallstones. The remaining 49 loci did not have any additional SNPs reaching genome-wide significance (*P* < 5 * 10−8). The *ABCG8* locus (rs11887534) contained 10 variants with independent and significant association with gallstones within *ABCG8, ABCG5* and *DYNC2LI1*. The *TM4SF4* locus (rs4681515) contained 2 variants, the *ABCB4* locus (rs7802555) contained 3 variants, the *SLC10A2* locus (rs144846334) contained 3 variants, the *FOXA3* locus (rs34255979) contained 2 variants and the *TNRC6B* locus (rs11089985) contained 2 variants with independent and significant association with gallstones. Full results of the conditional analyses are provided in Supplementary Materials 4. This study did not detect a second independent signal previously identified at the locus containing *TTC39B*9 although a more stringent genome-wide significance threshold was used in this study (*P =* 5 * 10−8 versus *P* = 6.49 * 10−4). ### Population Stratification After adjusting for the first 20 genetic principal components, there was evidence of test statistic inflation (*λ**gc* = 1.103). Test statistic inflation is seen in traits caused by multiple genes (polygenicity) but can also be a feature of unmeasured population substructure in which genetically-similar individuals also share similar environmental exposures32. Linkage disequilibrium score regression determines the extent to which polygenicity drives the test statistic inflation. The linkage disequilibrium score regression intercept was 1.042 and the proportion of test statistic inflation ascribed to causes other than polygenicity was estimated to be 10.9-17.2% confirming that polygenicity is the main driver of test statistic inflation. *λ**gc* was similar to the GBE GWAS12. The quantile-quantile plot for the UKB is shown in Supplementary Figure 5. ### Sensitivity Analysis A sensitivity analysis was conducted in each of the 55 lead SNPs using 23,422 individuals undergoing surgery for gallstones with 348,373 controls with no record of gallstones. 11 loci did not retain significance at the genome-wide threshold. Ten retained a suggestive association (*P* < 5 * 10−5; *ABO, ADAMTS20, FOXP1, GMDS-DT, GOT1*, chr18:rs542917327, *ITCH, NPAS2, SH3BP4, TBC1D32*) whilst one did not (*P* = 2.75 * 10−4; *LMNB1*). Direction and size of effect was very similar for all 55 loci suggesting detection bias had minimal influence on results. Full results at each locus are available in Supplementary Materials 4. A forest plot comparing effect sizes in shown in Supplementary Figure 6. ### Association of Lithogenic Loci with Phenotypic Traits Assessment of serum lipids was undertaken for each of the 55 lithogenic variants from the UKB in all 377,000 individuals taken forward for GWAS. The lithogenic alleles were heterogeneous in their influence on serum lipids with some lithogenic alleles promoting hyperlipidaemia whilst others demonstrated a protective effect (see Figure 3). The lithogenic variants of *ABCG8* and *HNF4A* (rs11887534-C and rs1800961-T respectively) both resulted in reductions of cholesterol of 0.10mmol/L (95% CI 0.09-0.12). The lithogenic variants of *GCKR* and *ABO* demonstrated highly significant increases in serum cholesterol (rs1260326-C: 0.05mmol/L [95% CI 0.05-0.06]; and rs115478735-T: 0.07 [0.06-0.07]). Figures for the other serum biochemistry markers are available in Supplementary Figures 7-12. Numeric results are shown in Supplementary Materials 5. ![Figure3](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/07/19/2021.07.16.21260637/F3.medium.gif) [Figure3](http://medrxiv.org/content/early/2021/07/19/2021.07.16.21260637/F3) ### Gene-Environment Interaction Gene-environment interaction of SNPs with serum unconjugated bilirubin was undertaken at 54 of the lead SNPs from the main GWAS (excluding *UGT1A1*) in the 377,000 individuals. After Bonferroni correction for 54 analyses, a significance level of 0.00093 was determined. No variants reached statistical significance although nominal significance was seen in *LIN28B* (P=0.012), *ABCB4* (P=0.023), *JMJD1C*(P=0.008), *ANPEP*(P=0.041), *CLDN7*(P=0.011), *ATP8B1* (P=0.019) and *HNF4A* (P=0.002). ### Polygenic Risk Score The PRS was significantly associated with gallstone disease in the validation cohort (maximum variance explained on the liability scale *R*2 = 2.5%, *P* = 2.62 * 10−181, when using a GWAS P value threshold of 0.0001). There was a 6-fold increase of risk of gallstones in those with the highest decile of genetic risk compared to those in the lowest with an approximately linear increase in risk with each increasing PRS decile (Figure 4). The PRS was strongly associated with increased BMI, waist circumference, hip circumference, liver enzymes (other than ALP), cystatin C and CRP (greater genetically-predicted risk of gallstones resulted in higher values for those traits). The PRS was associated with reductions in cholesterol, HDL, LDL, ApoB and ALP. Generally, for all traits with which the PRS was significantly associated there was an approximately linear increase (or decrease) with each increase in PRS decile and no apparent threshold effect. There were no obvious differences in levels of trigs, ApoA, LipoA, blood count markers, HbA1c or glucose (Supplementary Figures 13-15). ![Figure4](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/07/19/2021.07.16.21260637/F4.medium.gif) [Figure4](http://medrxiv.org/content/early/2021/07/19/2021.07.16.21260637/F4) ### Pathway Enrichment 172 genes showed significant association with gallstones (*P* < 2.59 * 10−6; Supplementary Materials 4) and were highly expressed in liver (*P* = 4.09 * 10−6; gallbladder expression was not available in MAGMA). 13 gene-sets were identified as significantly associated with gallstones (Table 3). Many of these gene-sets relate to glucoronidation, lipid metabolism and phospholipid transport. Additionally, Waldenström’s macroglobulinaemia and curated gene-sets of cancer drug resistance were identified which both contained *ABCB1* and *ABCB4*, both regulate biliary composition but are highly pleiotropic. View this table: [Table 3:](http://medrxiv.org/content/early/2021/07/19/2021.07.16.21260637/T3) Table 3: MAGMA Pathway Analysis - Significant Gene-Sets Literature search identified two genes governing gastrointestinal motility (*ANO1* and *TMEM147*). Rs7599 lies within the 3’UTR region and the lithogenic allele (A) is associated with significantly reduced expression of *TMEM147* in the liver (*P* = 3.0 * 10−47). Data on effects of rs56363382 and all other significant SNPs in the *ANO1* loci were not available from GTEx. Both *ANO1* and *TMEM147* were found to be highly expressed in gallbladder mucosa in the Human Protein Atlas. 5 identified genes were involved in primary cilia function (*DYNC2LI1, TBC1D32, ADAMTS20, PKD2L1* and *POC1B*). ## Discussion We report the largest and most detailed GWA of gallstones conducted worldwide with functional annotation of variants and gene-set analysis. We analysed 377,000 individuals (28,627 cases and 348,373 controls) from the UKB and replicated findings in a Scottish cohort (GS-SFHS) and a Finnish Cohort (FinnGen). We performed a genome-wide meta-analysis with UKB and FinnGen data. We identified 46 new lithogenic loci and replicated all 29 previously identified loci bringing the total number of identified gallstone-susceptibility loci to 75. Of the 75 loci, 12 were identified only in the UKB GWAS and not replicated or significant in the GWA meta-analysis. Through pathway enrichment we identified 13 gene-sets strongly associated with gallstones. Through FUMA annotation we identified two novel pathways that may influence gallstone development, gastrointestinal motility and ciliogenesis. Our results highlight important pathways that may be targeted by novel therapies and form the basis for further research into gallstones and other cholestatic disease. Two newly identified loci are highly expressed in the gallbladder and govern gastrointestinal motility. *ANO1* (also known as *TMEM16A*) is a calcium-activated chloride channel widely expressed in the gastrointestinal tract with roles in automaticity and contraction of the smooth muscle cells via the interstitial cells of Cajal33. Administration of *ANO1* inhibitors (TMinh-23) in murine models impedes gastric emptying with improved oral but not parenteral glucose tolerance supporting the role of *ANO1* as a gastric motility promoter33,34. *ANO1* is also expressed in the biliary epithelium where it has another role in secretion of fluid from the apical cholangiocyte membrane in response to bile acids suggesting it may promote increased bile flow35. *ANO1* may influence gallstone formation via gallbladder contractility and biliary flow. Given the possible impact on gallbladder motility and biliary composition, we suggest that any trials of TMinh-23 assess for development of gallstones as an important safety consideration. *TMEM147* is known to inhibit transfer of M3 muscarinic receptors to the cell surface36. M3 receptors are highly expressed in smooth muscle cells of the gallbladder and may influence gallbladder contractility and biliary stasis37. The variants may also influence gallstone formation through increased colonic transit time which is also associated with gallstones and believed to alter enterohepatic bile acid circulation14. Five loci relate to function of primary cilia. *TBC1D32* controls ciliogenesis within the neural tube and interacts with the sonic hedgehog protein which has wider roles in organogenesis and axonal guidance38. The metalloproteinase encoded by *ADAMTS20* also plays a pleiotropic role in ciliogenesis and hedgehog signaling with mutant variants causing ciliopathies39. *POC1B* is essential for formation of the centriole and ciliogenesis during cell division40. *PKD2L1* is involved in regulating calcium concentration within the primary cilia41 and disorders of calcium signaling of the primary cilia are thought to result in cholestasis42. Three independent signals within the *ABCG8* locus were intronic or 3’-UTR variants of *DYNC2LI1. DYNC2LI1* has a role in ciliogenesis with mutations causing severe ciliopathies43. Variants within the sonic hedgehog gene and *DYNC2LI1* have previously been associated with gallbladder cancer and gallstone formation44 whilst *TBC1D32* and *ADAMTS20* (novel loci, not replicated) and *PKD2L1* and *POC1B* (identified in the meta-analysis) are novel loci. We suggest that abnormalities in cholangiocyte ciliogenesis and ciliary function may contribute to biliary stasis with gallstone formation. Abnormal expression or function of cholangiocyte cilia causes other cholestatic disease and malignancies45. At present, therapies targeting primary cilia are limited. Ricolinostat, a selective inhibitor of histone deacetylase 6, is under investigation in several malignancies but its effect on gallstones remains unknown45. MAGMA analysis identified thirteen gene-sets. Most were related to lipid homeostasis and glucoronidation. Gene-sets for Waldenström’s macroglobulinaemia and multidrug-resistance pathways were enriched with enrichment driven by *ABCB4* and *ABCB1. ABCB4* encodes a phosphatidylcholine flippase which effluxes of phosphatidylcholine into bile46. Phosphatidylcholine acts as a solvent for biliary cholesterol. *ABCB4* mutations are a cause of early-onset gallstones and low-phospholipid-associated cholelithiasis47. *ABCB1* is a neighbouring gene involved in phospholipid and drug excretion with broad specificity46. The pathway enrichments are most likely a feature of pleiotropy as these genes have several distinct functions. However, serum IgG and IgM extracted from patients with Waldenström’s macroglobulinaemia results in nucleation of biliary cholesterol, possibly representing an alternate mechanism directly linking gallstones with Waldenström’s macroglobulinaemia48. Based on our results we suggest the following priorities for further research. 1. Direct genotyping at all 75 lead SNPs in individuals with corresponding cholecystectomy specimens for biochemical analysis of biliary content, gallbladder epithelial protein expression and cholangiocyte protein expression. 2. Assessment of gallbladder emptying with hepatobiliary scintigraphy in patients genotyped at *ANO1* and *TMEM147* to assess the influence of the lithogenic variants on gallbladder motility. 3. Assessment of gallstone risk in murine models with pharmacological modification of gallbladder motility (TMinh-23) or ciliotherapy (ricolinostat). 4. Validation of our findings in a large cohort in which radiological evaluation of gallstone disease has been undertaken. The strengths of this study include a more complete definition of gallstones than in previous studies which reduces the rate of case misclassification. Notably, previous studies have identified up to 18,417 cases whilst we identify 28,627 from the UKB. We use three large cohorts with replication of the majority of novel findings. Many of the loci we identify have plausible relationships to gallstones which we explore in pathway analysis, tissue expression analyses and through assessing the impact of risk alleles on serum biomarkers. The limitations of this study include misclassification of individuals with asymptomatic gallstones that have not been detected clinically. Relying on clinical disease as the definition of exposure has been common to all earlier GWAS of gallstones although our selection of relevant codes has improved our sensitivity relative to other studies. The phenotypic characterisation of participants in the UKB is limited to less invasive procedures such as venepuncture and no data on biliary composition or pathological records of gallstone composition following cholecystectomy were available. This impedes the interpretation of cholesterol-influencing SNPs as SNPs may lead to reduced serum cholesterol through biliary excretion, increased serum cholesterol through intestinal absorption or a balanced combination of both with no net change in serum cholesterol. All three scenarios may lead to similar changes in biliary cholesterol yet demonstrate opposing effects on serum cholesterol. Finally, the definition of gallstones in the FinnGen cohort was not identical to the UKB as we used the published summary association statistics from the FinnGen GWAS although the majority of UKB cases would have met the definition used in FinnGen. In summary, we have performed a GWAS of gallstone disease using UKB data, performed a meta-analysis with the FinnGen cohort and identified 46 novel variants associated with gallstone disease. We identify the first known gallstone-susceptibility variants that influence gastrointestinal motility, and we identify key priorities for research into gallstone formation and preventative treatments. ## Supporting information Supplementary Materials 1 [[supplements/260637_file07.docx]](pending:yes) Figure [[supplements/260637_file08.docx]](pending:yes) Supplementary Materials 1 [[supplements/260637_file09.docx]](pending:yes) Supplementary Materials 2A [[supplements/260637_file10.docx]](pending:yes) Supplementary Materials 3 [[supplements/260637_file11.docx]](pending:yes) ## Data Availability Data from the individuals studies is available on request to verified authors. ## Funding This work was funded by the Medical Research Council (MR/T008008/1). ## Competing Interests No authors declare any competing interests. ## Ethical Approval The UKB received ethical approval (Research Ethics Committee reference number: 11/NW/0382). Ethical approval was granted for Generation Scotland: Scottish Family Health Study by NHS Tayside Research Ethics Committee (reference number: 05/S1401/89). ## Author Contributions CJF conducted the genome-wide association study and drafted the manuscript. CJF, ADB, PKJ, DWC, PRHJT and JFW developed the code required for data analysis. CJF, TMD and JFW assisted with data access requests. RP, EMH, AI and JFW assisted with data storage and server maintenance. DJP, CH and AC contributed data from Generation Scotland. CJF, EMH, AS, SJW, CAS, TMD, RKS, NLM, CLMS and RP were involved in study design and conception. JAF, NCH, PR, NLM, RKS, SJW and EMH were involved in critical appraisal of the manuscript and placing the research in the biological context. ADB, AC, CH, AI, JFW and AS were involved in critical appraisal of the manuscript in relation to GWAS methodology. CJF, EMH, RKS, AS and AI were involved in pathway enrichment analyses. ## Acknowledgements Generation Scotland received core support from the Chief Scientist Office of the Scottish Government Health Directorates [CZD/16/6] and the Scottish Funding Council [HR03006] and is currently supported by the Wellcome Trust [216767/Z/19/Z]. Genotyping of the GS:SFHS samples was carried out by the Genetics Core Laboratory at the Edinburgh Clinical Research Facility, University of Edinburgh, Scotland and was funded by the Medical Research Council UK and the Wellcome Trust (Wellcome Trust Strategic Award “STratifying Resilience and Depression Longitudinally” (STRADL) Reference 104036/Z/14/Z). This work was funded by a Medical Research Council Clinical Research Training Fellowship awarded to CJF (MR/T008008/1). ADB acknowledges funding from the Wellcome PhD training fellowship for clinicians (204979/Z/16/Z), the Edinburgh Clinical Academic Track (ECAT) programme. CH and JFW acknowledges support from the MRC Human Genetics Unit programme grant, “Quantitative traits in health and disease” (U. MC_UU_00007/10). NCH is supported by a Wellcome Trust Senior Research Fellowship in Clinical Science (ref. 219542/Z/19/Z). NLM is supported by the British Heart Foundation through a Senior Clinical Research Fellowship (FS/16/14/32023) and a Research Excellence Award (RE/18/5/34216). RKS is supported by the Wellcome Trust (grant 210752/Z/18/Z). AS acknowledges support from the Academy of Medical Sciences/the Wellcome Trust/ the Government Department of Business, Energy and Industrial Strategy/the British Heart Foundation/Diabetes UK Springboard Award [SBF006\1109]. * Received July 16, 2021. * Revision received July 16, 2021. * Accepted July 19, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. 1.Harrison, E. M. et al. Hospital volume and patient outcomes after cholecystectomy in Scotland: Retrospective, national population based study. British Medical Journal 344, e3330s (2012). 2. 2.Everhart, J. E. & Ruhl, C. E. Burden of digestive diseases in the United States Part III: Liver, biliary tract, and pancreas. Gastroenterology 136, 1134–1144 (2009). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1053/j.gastro.2009.02.038&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19245868&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F19%2F2021.07.16.21260637.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000264716500005&link_type=ISI) 3. 3.Jørgensen, T. Prevalence of Gallstones in a Danish Population. American Journal of Epidemiology 126, 912–921 (1987). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=3310613&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F19%2F2021.07.16.21260637.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1987K494200014&link_type=ISI) 4. 4.Everhart, J. E., Khare, M., Hill, M. & Maurer, K. R. Prevalence and ethnic differences in gallbladder disease in the United States. Gastroenterology 117, 632–639 (1999). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0016-5085(99)70456-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10464139&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F19%2F2021.07.16.21260637.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000082261300021&link_type=ISI) 5. 5.Katsika, D. et al. Genetic and environmental influences on symptomatic gallstone disease: A Swedish study of 43,141 twin pairs. Hepatology 41, 1138–1143 (2005). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/hep.20654&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15747383&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F19%2F2021.07.16.21260637.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000228758800023&link_type=ISI) 6. 6.Buch, S. et al. A genome-wide association scan identifies the hepatic cholesterol transporter ABCG8 as a susceptibility factor for human gallstone disease. Nature Genetics 39, 995–999 (2007). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng2101&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17632509&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F19%2F2021.07.16.21260637.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000248446900017&link_type=ISI) 7. 7.Gudbjartsson, D. F. et al. Large-scale whole-genome sequencing of the Icelandic population. Nature Genetics 47, 435–444 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3247&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25807286&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F19%2F2021.07.16.21260637.atom) 8. 8.Joshi, A. D. et al. Four Susceptibility Loci for Gallstone Disease Identified in a Meta-analysis of Genome-Wide Association Studies. Gastroenterology 151, 351–363.e28 (2016). 9. 9.Rodriguez, S. et al. Lipids, obesity and gallbladder disease in women: Insights from genetic studies using the cardiovascular gene-centric 50K SNP array. European Journal of Human Genetics 24, 106–112 (2016). 10. 10.Ferkingstad, E. et al. Genome-wide association meta-analysis yields 20 loci associated with gallstone disease. Nature Communications 9, (2018). 11. 11.Bustos, B. I. et al. Variants in ABCG8 and TRAF3 genes confer risk for gallstone disease in admixed Latinos with Mapuche Native American ancestry. Scientific Reports 9, 1–12 (2019). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-019-40666-8&link_type=DOI) 12. 12.Gellert-Kristensen, H. et al. Identification and Replication of Six Loci Associated With Gallstone Disease. Hepatology 70, 597–609 (2019). 13. 13.Portincasa, P., Moschetta, A. & Palasciano, G. Cholesterol gallstone disease. The Lancet 368, 230–239 (2006). 14. 14.van Erpecum, K. J. & Berge-Henegouwen, G. P. V. Gallstones: An intestinal disease? Gut 44, 435–438 (1999). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ3V0am5sIjtzOjU6InJlc2lkIjtzOjg6IjQ0LzMvNDM1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDcvMTkvMjAyMS4wNy4xNi4yMTI2MDYzNy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 15. 15.McInnes, G. et al. Global Biobank Engine: Enabling genotype-phenotype browsing for biobank summary statistics. Bioinformatics 35, 2495–2497 (2019). 16. 16.Sudlow, C. et al. UK biobank: An open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLoS medicine 12, e1001779 (2015). 17. 17.Smith, B. H. et al. Cohort Profile: Generation Scotland: Scottish Family Health Study (GS:SFHS). The study, its participants and their potential for genetic research on health and illness. International Journal of Epidemiology 42, 689–700 (2013). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dys084&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22786799&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F19%2F2021.07.16.21260637.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000322955900012&link_type=ISI) 18. 18.FinnGen. FinnGen research project is an expedition to the frontier of genomics and medicine. FinnGen (2020). 19. 19.Manichaikul, A. et al. Robust relationship inference in genome-wide association studies. Bioinformatics 26, 2867–2873 (2010). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btq559&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20926424&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F19%2F2021.07.16.21260637.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000283919800010&link_type=ISI) 20. 20.Watanabe, K., Taskesen, E., van Bochoven, A. & Posthuma, D. Functional mapping and annotation of genetic associations with FUMA. Nature Communications 8, 1826 (2017). 21. 21.Bulik-Sullivan, B. K. et al. LD Score regression distinguishes confounding from polygenicity in genome-wide association studies. Nature Genetics 47, 291–295 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3211&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25642630&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F19%2F2021.07.16.21260637.atom) 22. 22.Buch, S. et al. Loci From a Genome-Wide Analysis of Bilirubin Levels Are Associated With Gallstone Risk and Composition. Gastroenterology 139, 1942–1951.e2 (2010). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1053/j.gastro.2010.09.003&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20837016&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F19%2F2021.07.16.21260637.atom) 23. 23.Fertrin, K. Y., Melo, M. B., Assis, Â.M., Saad, S. T. O. & Costa, F. F. UDP-glucuronosyltransferase 1 gene promoter polymorphism is associated with increased serum bilirubin levels and cholecystectomy in patients with sickle cell anemia. Clinical Genetics 64, 160–162 (2003). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1034/j.1399-0004.2003.00113.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12859413&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F19%2F2021.07.16.21260637.atom) 24. 24.Chaar, V. et al. Association of UGT1A1 polymorphism with prevalence and age at onset of cholelithiasis in sickle cell anemia. Haematologica 90, 188–199 (2005). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiaGFlbWF0b2wiO3M6NToicmVzaWQiO3M6ODoiOTAvMi8xODgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wNy8xOS8yMDIxLjA3LjE2LjIxMjYwNjM3LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 25. 25.Wasmuth, H. E. et al. Coinheritance of Gilbert syndromeAssociated UGT1A1 mutation increases gallstone risk in cystic fibrosis. Hepatology 43, 738–741 (2006). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16557566&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F19%2F2021.07.16.21260637.atom) 26. 26.Euesden, J., Lewis, C. M. & O’Reilly, P. F. PRSice: Polygenic Risk Score software. Bioinformatics 31, 1466–1468 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btu848&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25550326&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F19%2F2021.07.16.21260637.atom) 27. 27.Leeuw, C. A. de, Mooij, J. M., Heskes, T. & Posthuma, D. MAGMA: Generalized Gene-Set Analysis of GWAS Data. PLOS Computational Biology 11, e1004219 (2015). 28. 28.Marchini, J. & Howie, B. Genotype imputation for genome-wide association studies. Nature Reviews Genetics 11, 499–511 (2010). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrg2796&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20517342&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F19%2F2021.07.16.21260637.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000278998500012&link_type=ISI) 29. 29.Willer, C. J., Li, Y. & Abecasis, G. R. METAL: Fast and efficient meta-analysis of genomewide association scans. Bioinformatics 26, 2190–2191 (2010). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btq340&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20616382&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F19%2F2021.07.16.21260637.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000281738900017&link_type=ISI) 30. 30.R Core Team. R: A Language and Environment for Statistical Computing. (2017). 31. 31.Little, J. et al. STrengthening the REporting of Genetic Association Studies (STREGA) of the STROBE Statement. PLOS Medicine 6, e1000022 (2009). 32. 32.Yang, J. et al. Genomic inflation factors under polygenic inheritance. European Journal of Human Genetics 19, 807–812 (2011). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ejhg.2011.39&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21407268&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F19%2F2021.07.16.21260637.atom) 33. 33.Mazzone, A. et al. Altered Expression of Ano1 Variants in Human Diabetic Gastroparesis. The Journal of Biological Chemistry 286, 13393–13403 (2011). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEyOiIyODYvMTUvMTMzOTMiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wNy8xOS8yMDIxLjA3LjE2LjIxMjYwNjM3LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 34. 34.Cil, O. et al. Slowed gastric emptying and improved oral glucose tolerance produced by a nanomolar-potency inhibitor of calcium-activated chloride channel TMEM16A. The FASEB Journal 33, 11247–11257 (2019). 35. 35.Li, Q., Dutta, A., Kresge, C., Bugde, A. & Feranchak, A. P. Bile acids stimulate cholangiocyte fluid secretion by activation of transmembrane member 16A Cl-channels. Hepatology (Baltimore, Md.) 68, 187–199 (2018). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/hep.29804&link_type=DOI) 36. 36.Rosemond, E. et al. Regulation of M3 muscarinic receptor expression and function by transmembrane protein 147. Molecular Pharmacology 79, 251–261 (2011). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoibW9scGhhcm0iO3M6NToicmVzaWQiO3M6ODoiNzkvMi8yNTEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wNy8xOS8yMDIxLjA3LjE2LjIxMjYwNjM3LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 37. 37.Parkman, H. P., Pagano, A. P. & Ryan, J. P. Subtypes of muscarinic receptors regulating gallbladder cholinergic contractions. American Journal of Physiology-Gastrointestinal and Liver Physiology 276, G1243–G1250 (1999). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10330016&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F19%2F2021.07.16.21260637.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000080344400022&link_type=ISI) 38. 38.Ko, H. W. et al. Broad-minded links cell cycle-related kinase to cilia assembly and Hedgehog signal transduction. Developmental cell 18, 237–247 (2010). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.devcel.2009.12.014&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20159594&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F19%2F2021.07.16.21260637.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000275100400010&link_type=ISI) 39. 39.Nandadasa, S. et al. Secreted metalloproteases ADAMTS9 and ADAMTS20 have a non-canonical role in ciliary vesicle growth during ciliogenesis. Nature Communications 10, 953 (2019). 40. 40.Atorino, E. S., Hata, S., Funaya, C., Neuner, A. & Schiebel, E. CEP44 ensures the formation of bona fide centriole wall, a requirement for the centriole-to-centrosome conversion. Nature Communications 11, 903 (2020). 41. 41.DeCaen, P. G., Delling, M., Vien, T. N. & Clapham, D. E. Direct recording and molecular identification of the calcium channel of primary cilia. Nature 504, 315–318 (2013). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature12832&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24336289&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F19%2F2021.07.16.21260637.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000328121500044&link_type=ISI) 42. 42.Amaya, M. J. & Nathanson, M. H. Calcium signaling and the secretory activity of bile duct epithelia. Cell Calcium 55, 317–324 (2014). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ceca.2014.02.003&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24612866&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F19%2F2021.07.16.21260637.atom) 43. 43.Taylor, S. P. et al. Mutations in DYNC2LI1 disrupt cilia function and cause short rib polydactyly syndrome. Nature Communications 6, 7092 (2015). 44. 44.Dixit, R., Pandey, M., Tripathi, S. K., Dwivedi, A. N. D. & Shukla, V. K. Comparative Analysis of Mutational Profile of Sonic hedgehog Gene in Gallbladder Cancer. Digestive Diseases and Sciences 62, 708–714 (2017). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10620-016-4438-1&link_type=DOI) 45. 45.Gradilone, S. A. et al. HDAC6 Is Overexpressed in Cystic Cholangiocytes and Its Inhibition Reduces Cystogenesis. The American Journal of Pathology 184, 600–608 (2014). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajpath.2013.11.027&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24434010&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F19%2F2021.07.16.21260637.atom) 46. 46.van Helvoort, A. et al. MDR1 P-Glycoprotein Is a Lipid Translocase of Broad Specificity, While MDR3 P-Glycoprotein Specifically Translocates Phosphatidylcholine. Cell 87, 507– 517 (1996). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0092-8674(00)81370-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8898203&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F19%2F2021.07.16.21260637.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1996VQ46600016&link_type=ISI) 47. 47.Rosmorduc, O. et al. ABCB4 gene mutation-associated cholelithiasis in adults. Gastroenterology 125, 452–459 (2003). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0016-5085(03)00898-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12891548&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F19%2F2021.07.16.21260637.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000184531200025&link_type=ISI) 48. 48.Harvey, P. R., Upadhya, G. A. & Strasberg, S. M. Immunoglobulins as nucleating proteins in the gallbladder bile of patients with cholesterol gallstones. Journal of Biological Chemistry 266, 13996–14003 (1991). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEyOiIyNjYvMjEvMTM5OTYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wNy8xOS8yMDIxLjA3LjE2LjIxMjYwNjM3LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) [1]: /embed/graphic-1.gif