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.

Kinetics of moisture loss applied to the baking of snacks with pregelatinized cassava starch

Moisture loss kinetics is a complex process defined as the liquid removal from a solid by thermal application. The purposes of the study were to obtain kinetic curves of moisture loss during the baking of cassava snacks and establish which processes govern the moisture loss, recognize which mathematical models describe the moisture loss curves more precisely, and determine activation energy (Ea) and effective diffusivity (D_eff ). Experimental data were obtained through baking at four temperatures formulations for snacks with different dehydrated cassava puree (DCP) proportions. Page's and Chávez Méndez's models showed the best fits. We calculated D_eff and Ea employing the analytical solution of Fick's Second Law for the geometry of plane plates. D_eff values increase with DCP but did not show a trend. The range found was from 5.22E-06 to 2.93E-05 m² /s. The results of Ea showed that the mixture of flours produced an increase in the energy necessary to initiate the effective diffusion (24.84 kJ/mol) compared to the samples without mixing (15.54 kJ/mol). Moisture loss curves show that the diffusion process governs a large part of the process. PRACTICAL APPLICATION: Given the need to increase research for the development of the cassava industry, which currently has low competitiveness compared to less expensive products such as corn, various efforts are being made to generate new products that can replace wheat flour, at least in part. However, it is necessary to research how this substitution affects the various steps of the production system, including baking. During baking, one of the most significant processes is moisture loss. In this sense, the kinetic modeling of the moisture loss process parameters is mainly helpful in the food industry. The mathematical models of moisture loss processes are used to design new or improved baking systems or even control the process.

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.

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.

Depressive symptoms evaluated by the Calgary Depression Scale for Schizophrenia (CDSS): genetic vulnerability and sex effects

The present study compares the occurrence of depressive symptoms evaluated by the Calgary Depression Scale for Schizophrenia (CDSS) in patients of Multiplex (MS) and Simplex Schizophrenia families (SS). The Positive and Negative Syndrome Scale (PANSS) was used to evaluate psychopathology. A total of 206 paranoid schizophrenia patients were studied according DSM-IV criteria. The Family Interview for Genetic Studies (FIGS) was used to study the families. A result in the FIGS for a positive family history of schizophrenia was referred as MS (patients); its lack as SS (patients). CDSS scores were compared among MS and SS patients and possible sex differences intra- and inter-groups were explored. In the analysis of our sample (30) 19% of the total persons with schizophrenia group was depressed. The depressive symptoms measured by the CDSS were higher in females and the MS males group. Males from MS group showed more depressive symptoms than males from SS group. No differences with females from both groups were found. Findings in this study underscore the importance of gender and family history in understanding the heterogeneity of schizophrenia. This study suggests that sex and familiar history are important to consider in studying depressive symptoms.

Statistical properties of the variation at linked microsatellite loci: implications for the history of human Y chromosomes

It has recently been suggested that observed levels of variation at microsatellite loci can be used to infer patterns of selection in genomes and to assess demographic history. In order to evaluate the feasibility of these suggestions it is necessary to know something about how levels of variation at microsatellite loci are expected to fluctuate due simply to stochasticity in the processes of mutation and inheritance (genetic sampling). Here we use recently derived properties of the stepwise mutation model to place confidence intervals around the variance in repeat score that is expected at mutation-drift equilibrium and outline a statistical test for whether an observed value differs significantly from expectation. We also develop confidence intervals for the time course of the buildup of variation following a complete elimination of variation, such as might be caused by a selective sweep or an extreme population bottleneck. We apply these methods to the variation observed at human Y-specific microsatellites. Although a number of authors have suggested the possibility of a very recent sweep, our analyses suggest that a sweep or extreme bottleneck is unlikely to have occurred anytime during the last approximately 74,000 years. To generate this result we use a recently estimated mutation rate for microsatellite loci of 5.6 x 10(-4) along with the variation observed at autosomal microsatellite loci to estimate the human effective population size. This estimate is 18,000, implying an effective number of 4,500 Y chromosomes. One important general conclusion to emerge from this study is that in order to reject mutation-drift equilibrium at a set of linked microsatellite loci it is necessary to have an unreasonably large number of loci unless the observed variance is far below that expected at mutation-drift equilibrium.

Aspects of nonrandom turnover involved in the concerted evolution of intergenic spacers within the ribosomal DNA of Drosophila melanogaster

Polymerase chain reaction (PCR)-amplified, sequenced, and digitally typed intergenic spacers (IGSs) of the ribosomal (r)DNA in D. melanogaster reveal unexpected features of the mechanisms of turnover involved with the concerted evolution of the gene family. Characterization of the structure of three isolated IGS length variants reveals breakage "hot spots" within the 330-base-pair (bp) subrepeat array found in the spacers. Internal mapping of variant repeats within the 240-bp subrepeat array using a novel digital DNA typing procedure (minisatellite variant repeat [MVR]-PCR) shows an unexpected pattern of clustering of variant repeats. Each 240-bp subrepeat array consists of essentially two halves with the repeats in each half identified by specific mutations. This bipartite structure, observed in a cloned IGS unit, in the majority of genomic DNA of laboratory and wild flies and in PCR-amplified products, has been widely homogenized yet is not predicted by a model of unequal crossing over with randomly placed recombination breakpoints. Furthermore, wild populations contain large numbers of length variants in contrast to uniformly shared length variants in laboratory stocks. High numbers of length variants coupled to the observation of a homogenized bipartite structure of the 240-bp subrepeat array suggest that the unit of turnover and homogenization is smaller than the IGS and might involve gene conversion. The use of PCR for the structural analysis of members of the rDNA gene family coupled to digital DNA typing provides powerful new inroads into the mechanisms of DNA turnover affecting the course of molecular evolution in this family.