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.

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.

Locoregional risk assessment after neoadjuvant chemotherapy in patients with primary breast cancer: clinical utility of the CPS + EG score.

PURPOSE: Locoregional control is a prerequisite to cure primary breast cancer but the prediction of locoregional recurrence to guide further local therapy following neoadjuvant chemotherapy remains a challenge. The CPS + EG score was designed to predict distant recurrences. Here we examine its ability to predict both not only distant but also locoregional recurrences with respect to accuracy and clinical applicability. METHODS: Clinical data from 432 patients with primary breast cancer treated with neoadjuvant chemotherapy between 2003 and 2011 were prospectively collected. Using the Kaplan-Meier method we analyzed the risk of local and distant recurrences according to individual CPS + EG scores, stratified by type of surgery. Possible confounding of the relationship between recurrence risk and CPS + EG score by established risk factors was accounted for in multiple survival regression models. Additionally, we analyzed the performance of the CPS + EG score to predict isolated locoregional recurrence by censoring patients with prior or simultaneous distant metastases. RESULTS: 5-year locoregional recurrence-free survival was 90%, and 5-year distant metastases-free survival was 82%. The CPS + EG score stratified patients into six prognostic groups with distinct 5-year locoregional recurrence-free survival, ranging from 100 to 41% (p = 0.02) and 5-year distant metastases-free survival, ranging from 96 to 35% (p < 0.0001). 8 patients (17%) with CPS + EG scores ≥ 4 experienced locoregional recurrence-5 of them presented with simultaneous distant disease. CONCLUSION: The CPS + EG score, originally designed to predict distant relapse, is also valuable for assessing local recurrence risks. Our data demonstrate that distant and locoregional recurrence risks are closely related. As prognosis of patients with high risk of locoregional failure based on CPS + EG is dominated by distant recurrences, escalating local therapies may have limited impact on overall prognosis.

miR-16 and miR-125b are involved in barrier function dysregulation through the modulation of claudin-2 and cingulin expression in the jejunum in IBS with diarrhoea.

OBJECTIVE: Micro-RNAs (miRNAs) play a crucial role in controlling intestinal epithelial barrier function partly by modulating the expression of tight junction (TJ) proteins. We have previously shown differential messenger RNA (mRNA) expression correlated with ultrastructural abnormalities of the epithelial barrier in patients with diarrhoea-predominant IBS (IBS-D). However, the participation of miRNAs in these differential mRNA-associated findings remains to be established. Our aims were (1) to identify miRNAs differentially expressed in the small bowel mucosa of patients with IBS-D and (2) to explore putative target genes specifically involved in epithelial barrier function that are controlled by specific dysregulated IBS-D miRNAs. DESIGN: Healthy controls and patients meeting Rome III IBS-D criteria were studied. Intestinal tissue samples were analysed to identify potential candidates by: (a) miRNA-mRNA profiling; (b) miRNA-mRNA pairing analysis to assess the co-expression profile of miRNA-mRNA pairs; (c) pathway analysis and upstream regulator identification; (d) miRNA and target mRNA validation. Candidate miRNA-mRNA pairs were functionally assessed in intestinal epithelial cells. RESULTS: IBS-D samples showed distinct miRNA and mRNA profiles compared with healthy controls. TJ signalling was associated with the IBS-D transcriptional profile. Further validation of selected genes showed consistent upregulation in 75% of genes involved in epithelial barrier function. Bioinformatic analysis of putative miRNA binding sites identified hsa-miR-125b-5p and hsa-miR-16 as regulating expression of the TJ genes CGN (cingulin) and CLDN2 (claudin-2), respectively. Consistently, protein expression of CGN and CLDN2 was upregulated in IBS-D, while the respective targeting miRNAs were downregulated. In addition, bowel dysfunction, perceived stress and depression and number of mast cells correlated with the expression of hsa-miR-125b-5p and hsa-miR-16 and their respective target proteins. CONCLUSIONS: Modulation of the intestinal epithelial barrier function in IBS-D involves both transcriptional and post-transcriptional mechanisms. These molecular mechanisms include miRNAs as master regulators in controlling the expression of TJ proteins and are associated with major clinical symptoms.

NGS allele counts versus called genotypes for testing genetic association.

RNA sequence data are commonly summarized as read counts. By contrast, so far there is no alternative to genotype calling for investigating the relationship between genetic variants determined by next-generation sequencing (NGS) and a phenotype of interest. Here we propose and evaluate the direct analysis of allele counts for genetic association tests. Specifically, we assess the potential advantage of the ratio of alternative allele counts to the total number of reads aligned at a specific position of the genome (coverage) over called genotypes. We simulated association studies based on NGS data from HapMap individuals. Genotype quality scores and allele counts were simulated using NGS data from the Personal Genome Project. Real data from the 1000 Genomes Project was also used to compare the two competing approaches. The average proportions of probability values lower or equal to 0.05 amounted to 0.0496 for called genotypes and 0.0485 for the ratio of alternative allele counts to coverage in the null scenario, and to 0.69 for called genotypes and 0.75 for the ratio of alternative allele counts to coverage in the alternative scenario (9% power increase). The advantage in statistical power of the novel approach increased with decreasing coverage, with decreasing genotype quality and with decreasing allele frequency - 124% power increase for variants with a minor allele frequency lower than 0.05. We provide computer code in R to implement the novel approach, which does not preclude the use of complementary data quality filters before or after identification of the most promising association signals. Author summary: Genetic association tests usually rely on called genotypes. We postulate here that the direct analysis of allele counts from sequence data improves the quality of statistical inference. To evaluate this hypothesis, we investigate simulated and real data using distinct statistical approaches. We demonstrate that association tests based on allele counts rather than called genotypes achieve higher statistical power with controlled type I error rates.

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

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

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.

Prediction of local recurrence risk after neoadjuvant chemotherapy in patients with primary breast cancer: Clinical utility of the MD Anderson Prognostic Index.

BACKGROUND: Locoregional recurrence after neoadjuvant chemotherapy for primary breast cancer is associated with poor prognosis. It is essential to identify patients at high risk of locoregional recurrence who may benefit from extended local therapy. Here, we examined the prediction accuracy and clinical applicability of the MD Anderson Prognostic Index (MDAPI). METHODS: Prospective clinical data from 456 patients treated between 2003 and 2011 was analyzed. The Kaplan-Meier method was used to examine the probabilities of locoregional recurrence, local recurrence and distant metastases according to individual prognosis score, stratified by type of surgery (breast conserving therapy or mastectomy). The possible confounding of the relationship between recurrence risk and MDAPI by established risk factors was accounted for in multiple survival regression models. To define the clinical utility of the MDAPI we analyzed its performance to predict locoregional recurrence censoring patients with prior or simultaneous distant metastases. RESULTS: Mastectomized patients (42% of the patients) presented with more advanced tumor stage, lower tumor grade, hormone-receptor positive disease and consequently lower pathological complete response rates. Only a few patients presented with high-risk scores (2,7% MDAPI≥3). All patients with high-risk MDAPI score (MDAPI ≥3) who developed locoregional recurrence were simultaneously affected by distant metastases. CONCLUSION: Our data do not support a clinical utility of the MDAPI to guide local therapy.

Using next-generation DNA sequence data for genetic association tests based on allele counts with and without consideration of zero inflation.

The relationship between genetic variability and individual phenotypes is usually investigated by testing for association relying on called genotypes. Allele counts obtained from next-generation sequence data could be used for this purpose too. Genetic association can be examined by treating alternative allele counts (AACs) as the response variable in negative binomial regression. AACs from sequence data often contain an excess of zeros, thus motivating the use of Hurdle and zero-inflated models. Here we examine rough type I error rates and the ability to pick out variants with small probability values for 7 different testing approaches that incorporate AACs as an explanatory or as a response variable. Model comparisons relied on chromosome 3 DNA sequence data from 407 Hispanic participants in the Type 2 Diabetes Genetic Exploration by Next-generation sequencing in Ethnic Samples (T2D-GENES) project 1 with complete information on diastolic blood pressure and related medication. Our results suggest that in the investigation of the relationship between AAC as response variable and individual phenotypes as explanatory variable, Hurdle-negative binomial regression has some advantages. This model showed a good ability to discriminate strongly associated variants and controlled overall type I error rates. However, probability values from Hurdle-negative binomial regression were not obtained for approximately 25 % of the investigated variants because of convergence problems, and the mass of the probability value distribution was concentrated around 1.