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.

Combined genome-wide association study of 136 quantitative ear morphology traits in multiple populations reveal 8 novel loci

Human ear morphology, a complex anatomical structure represented by a multidimensional set of correlated and heritable phenotypes, has a poorly understood genetic architecture. In this study, we quantitatively assessed 136 ear morphology traits using deep learning analysis of digital face images in 14,921 individuals from five different cohorts in Europe, Asia, and Latin America. Through GWAS meta-analysis and C-GWASs, a recently introduced method to effectively combine GWASs of many traits, we identified 16 genetic loci involved in various ear phenotypes, eight of which have not been previously associated with human ear features. Our findings suggest that ear morphology shares genetic determinants with other surface ectoderm-derived traits such as facial variation, mono eyebrow, and male pattern baldness. Our results enhance the genetic understanding of human ear morphology and shed light on the shared genetic contributors of different surface ectoderm-derived phenotypes. Additionally, gene editing experiments in mice have demonstrated that knocking out the newly ear-associated gene (Intu) and a previously ear-associated gene (Tbx15) causes deviating mouse ear morphology.

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.

Bulsarapp: Interactive Visual Analysis for Surname Trend Exploration

The study of surnames for a given population, together with their distribution and spatial patterns identification, has been a long-standing problem in the fields of human biology, public health, and social sciences. The ancestry inferred from surname information can be a useful means to understand the dynamics of human populations. This knowledge allows us to characterize geographically the ethnicity of populations, and to understand the complex relationships between identity, migration, and health issues in a demographic view. However, in most cases, a detailed geolocalization of this data can be a daunting task. We propose a visual analytic tool that summarizes the heterogeneous surname and geographic information collected from Argentinean electoral rolls. This tool allows a massive data analysis, and facilitates interdisciplinary studies about population dynamics related to ancestry, migration, and health. It also offers an easy-to-use interface that allows interactive exploration of isonymy and surname origins, their distribution, and spatial trends in a high population density context.

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 (r² = 0.26) and three cis-C22orf34-eQTLs (r² = 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 new-born (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.

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.

Novel genetic loci affecting facial shape variation in humans

The human face represents a combined set of highly heritable phenotypes, but knowledge on its genetic architecture remains limited, despite the relevance for various fields. A series of genome-wide association studies on 78 facial shape phenotypes quantified from 3-dimensional facial images of 10,115 Europeans identified 24 genetic loci reaching study-wide suggestive association (p < 5 × 10-8), among which 17 were previously unreported. A follow-up multi-ethnic study in additional 7917 individuals confirmed 10 loci including six unreported ones (padjusted < 2.1 × 10-3). A global map of derived polygenic face scores assembled facial features in major continental groups consistent with anthropological knowledge. Analyses of epigenomic datasets from cranial neural crest cells revealed abundant cis-regulatory activities at the face-associated genetic loci. Luciferase reporter assays in neural crest progenitor cells highlighted enhancer activities of several face-associated DNA variants. These results substantially advance our understanding of the genetic basis underlying human facial variation and provide candidates for future in-vivo functional studies.

Meta-analysis of genome-wide association studies identifies 8 novel loci involved in shape variation of human head hair

Shape variation of human head hair shows striking variation within and between human populations, while its genetic basis is far from being understood. We performed a series of genome-wide association studies (GWASs) and replication studies in a total of 28 964 subjects from 9 cohorts from multiple geographic origins. A meta-analysis of three European GWASs identified 8 novel loci (1p36.23 ERRFI1/SLC45A1, 1p36.22 PEX14, 1p36.13 PADI3, 2p13.3 TGFA, 11p14.1 LGR4, 12q13.13 HOXC13, 17q21.2 KRTAP, and 20q13.33 PTK6), and confirmed 4 previously known ones (1q21.3 TCHH/TCHHL1/LCE3E, 2q35 WNT10A, 4q21.21 FRAS1, and 10p14 LINC00708/GATA3), all showing genome-wide significant association with hair shape (P < 5e-8). All except one (1p36.22 PEX14) were replicated with nominal significance in at least one of the 6 additional cohorts of European, Native American and East Asian origins. Three additional previously known genes (EDAR, OFCC1, and PRSS53) were confirmed at the nominal significance level. A multivariable regression model revealed that 14 SNPs from different genes significantly and independently contribute to hair shape variation, reaching a cross-validated AUC value of 0.66 (95% CI: 0.62–0.70) and an AUC value of 0.64 in an independent validation cohort, providing an improved accuracy compared with a previous model. Prediction outcomes of 2504 individuals from a multiethnic sample were largely consistent with general knowledge on the global distribution of hair shape variation. Our study thus delivers target genes and DNA variants for future functional studies to further evaluate the molecular basis of hair shape in humans.

Developmental pathways inferred from modularity, morphological integration and fluctuating asymmetry patterns in the human face

Facial asymmetries are usually measured and interpreted as proxies to developmental noise. However, analyses focused on its developmental and genetic architecture are scarce. To advance on this topic, studies based on a comprehensive and simultaneous analysis of modularity, morphological integration and facial asymmetries including both phenotypic and genomic information are needed. Here we explore several modularity hypotheses on a sample of Latin American mestizos, in order to test if modularity and integration patterns differ across several genomic ancestry backgrounds. To do so, 4104 individuals were analyzed using 3D photogrammetry reconstructions and a set of 34 facial landmarks placed on each individual. We found a pattern of modularity and integration that is conserved across sub-samples differing in their genomic ancestry background. Specifically, a signal of modularity based on functional demands and organization of the face is regularly observed across the whole sample. Our results shed more light on previous evidence obtained from Genome Wide Association Studies performed on the same samples, indicating the action of different genomic regions contributing to the expression of the nose and mouth facial phenotypes. Our results also indicate that large samples including phenotypic and genomic metadata enable a better understanding of the developmental and genetic architecture of craniofacial phenotypes.