A multi-tissue, splicing-based joint transcriptome-wide association study identifies susceptibility genes for breast cancer.

Splicing-based transcriptome-wide association studies (splicing-TWASs) of breast cancer have the potential to identify susceptibility genes. However, existing splicing-TWASs test the association of individual excised introns in breast tissue only and thus have limited power to detect susceptibility genes. In this study, we performed a multi-tissue joint splicing-TWAS that integrated splicing-TWAS signals of multiple excised introns in each gene across 11 tissues that are potentially relevant to breast cancer risk. We utilized summary statistics from a meta-analysis that combined genome-wide association study (GWAS) results of 424,650 women of European ancestry. Splicing-level prediction models were trained in GTEx (v.8) data. We identified 240 genes by the multi-tissue joint splicing-TWAS at the Bonferroni-corrected significance level; in the tissue-specific splicing-TWAS that combined TWAS signals of excised introns in genes in breast tissue only, we identified nine additional significant genes. Of these 249 genes, 88 genes in 62 loci have not been reported by previous TWASs, and 17 genes in seven loci are at least 1 Mb away from published GWAS index variants. By comparing the results of our splicing-TWASs with previous gene-expression-based TWASs that used the same summary statistics and expression prediction models trained in the same reference panel, we found that 110 genes in 70 loci that are identified only by the splicing-TWASs. Our results showed that for many genes, expression quantitative trait loci (eQTL) did not show a significant impact on breast cancer risk, whereas splicing quantitative trait loci (sQTL) showed a strong impact through intron excision events.

Predicted Proteome Association Studies of Breast, Prostate, Ovarian, and Endometrial Cancers Implicate Plasma Protein Regulation in Cancer Susceptibility.

BACKGROUND: Predicting protein levels from genotypes for proteome-wide association studies (PWAS) may provide insight into the mechanisms underlying cancer susceptibility. METHODS: We performed PWAS of breast, endometrial, ovarian, and prostate cancers and their subtypes in several large European-ancestry discovery consortia (effective sample size: 237,483 cases/317,006 controls) and tested the results for replication in an independent European-ancestry GWAS (31,969 cases/410,350 controls). We performed PWAS using the cancer GWAS summary statistics and two sets of plasma protein prediction models, followed by colocalization analysis. RESULTS: Using Atherosclerosis Risk in Communities (ARIC) models, we identified 93 protein-cancer associations [false discovery rate (FDR) < 0.05]. We then performed a meta-analysis of the discovery and replication PWAS, resulting in 61 significant protein-cancer associations (FDR < 0.05). Ten of 15 protein-cancer pairs that could be tested using Trans-Omics for Precision Medicine (TOPMed) protein prediction models replicated with the same directions of effect in both cancer GWAS (P < 0.05). To further support our results, we applied Bayesian colocalization analysis and found colocalized SNPs for SERPINA3 protein levels and prostate cancer (posterior probability, PP = 0.65) and SNUPN protein levels and breast cancer (PP = 0.62). CONCLUSIONS: We used PWAS to identify potential biomarkers of hormone-related cancer risk. SNPs in SERPINA3 and SNUPN did not reach genome-wide significance for cancer in the original GWAS, highlighting the power of PWAS for novel locus discovery, with the added advantage of providing directions of protein effect. IMPACT: PWAS and colocalization are promising methods to identify potential molecular mechanisms underlying complex traits.

A joint transcriptome-wide association study across multiple tissues identifies candidate breast cancer susceptibility genes.

Genome-wide association studies (GWASs) have identified more than 200 genomic loci for breast cancer risk, but specific causal genes in most of these loci have not been identified. In fact, transcriptome-wide association studies (TWASs) of breast cancer performed using gene expression prediction models trained in breast tissue have yet to clearly identify most target genes. To identify candidate genes, we performed a GWAS analysis in a breast cancer dataset from UK Biobank (UKB) and combined the results with the GWAS results of the Breast Cancer Association Consortium (BCAC) by a meta-analysis. Using the summary statistics from the meta-analysis, we performed a joint TWAS analysis that combined TWAS signals from multiple tissues. We used expression prediction models trained in 11 tissues that are potentially relevant to breast cancer from the Genotype-Tissue Expression (GTEx) data. In the GWAS analysis, we identified eight loci distinct from those reported previously. In the TWAS analysis, we identified 309 genes at 108 genomic loci to be significantly associated with breast cancer at the Bonferroni threshold. Of these, 17 genes were located in eight regions that were at least 1 Mb away from published GWAS hits. The remaining TWAS-significant genes were located in 100 known genomic loci from previous GWASs of breast cancer. We found that 21 genes located in known GWAS loci remained statistically significant after conditioning on previous GWAS index variants. Our study provides insights into breast cancer genetics through mapping candidate target genes in a large proportion of known GWAS loci and discovering multiple new loci.

Polygenic transcriptome risk scores for COPD and lung function improve cross-ethnic portability of prediction in the NHLBI TOPMed program.

While polygenic risk scores (PRSs) enable early identification of genetic risk for chronic obstructive pulmonary disease (COPD), predictive performance is limited when the discovery and target populations are not well matched. Hypothesizing that the biological mechanisms of disease are shared across ancestry groups, we introduce a PrediXcan-derived polygenic transcriptome risk score (PTRS) to improve cross-ethnic portability of risk prediction. We constructed the PTRS using summary statistics from application of PrediXcan on large-scale GWASs of lung function (forced expiratory volume in 1 s [FEV1] and its ratio to forced vital capacity [FEV1/FVC]) in the UK Biobank. We examined prediction performance and cross-ethnic portability of PTRS through smoking-stratified analyses both on 29,381 multi-ethnic participants from TOPMed population/family-based cohorts and on 11,771 multi-ethnic participants from TOPMed COPD-enriched studies. Analyses were carried out for two dichotomous COPD traits (moderate-to-severe and severe COPD) and two quantitative lung function traits (FEV1 and FEV1/FVC). While the proposed PTRS showed weaker associations with disease than PRS for European ancestry, the PTRS showed stronger association with COPD than PRS for African Americans (e.g., odds ratio [OR] = 1.24 [95% confidence interval [CI]: 1.08-1.43] for PTRS versus 1.10 [0.96-1.26] for PRS among heavy smokers with ≥ 40 pack-years of smoking) for moderate-to-severe COPD. Cross-ethnic portability of the PTRS was significantly higher than the PRS (paired t test p < 2.2 × 10-16 with portability gains ranging from 5% to 28%) for both dichotomous COPD traits and across all smoking strata. Our study demonstrates the value of PTRS for improved cross-ethnic portability compared to PRS in predicting COPD risk.

A Trans-Ethnic Genome-Wide Association Study of Uterine Fibroids.

Uterine fibroids affect up to 77% of women by menopause and account for up to $34 billion in healthcare costs each year. Although fibroid risk is heritable, genetic risk for fibroids is not well understood. We conducted a two-stage case-control meta-analysis of genetic variants in European and African ancestry women with and without fibroids classified by a previously published algorithm requiring pelvic imaging or confirmed diagnosis. Women from seven electronic Medical Records and Genomics (eMERGE) network sites (3,704 imaging-confirmed cases and 5,591 imaging-confirmed controls) and women of African and European ancestry from UK Biobank (UKB, 5,772 cases and 61,457 controls) were included in the discovery genome-wide association study (GWAS) meta-analysis. Variants showing evidence of association in Stage I GWAS (P < 1 × 10-5) were targeted in an independent replication sample of African and European ancestry individuals from the UKB (Stage II) (12,358 cases and 138,477 controls). Logistic regression models were fit with genetic markers imputed to a 1000 Genomes reference and adjusted for principal components for each race- and site-specific dataset, followed by fixed-effects meta-analysis. Final analysis with 21,804 cases and 205,525 controls identified 326 genome-wide significant variants in 11 loci, with three novel loci at chromosome 1q24 (sentinel-SNP rs14361789; P = 4.7 × 10-8), chromosome 16q12.1 (sentinel-SNP rs4785384; P = 1.5 × 10-9) and chromosome 20q13.1 (sentinel-SNP rs6094982; P = 2.6 × 10-8). Our statistically significant findings further support previously reported loci including SNPs near WT1, TNRC6B, SYNE1, BET1L, and CDC42/WNT4. We report evidence of ancestry-specific findings for sentinel-SNP rs10917151 in the CDC42/WNT4 locus (P = 1.76 × 10-24). Ancestry-specific effect-estimates for rs10917151 were in opposite directions (P-Het-between-groups = 0.04) for predominantly African (OR = 0.84) and predominantly European women (OR = 1.16). Genetically-predicted gene expression of several genes including LUZP1 in vagina (P = 4.6 × 10-8), OBFC1 in esophageal mucosa (P = 8.7 × 10-8), NUDT13 in multiple tissues including subcutaneous adipose tissue (P = 3.3 × 10-6), and HEATR3 in skeletal muscle tissue (P = 5.8 × 10-6) were associated with fibroids. The finding for HEATR3 was supported by SNP-based summary Mendelian randomization analysis. Our study suggests that fibroid risk variants act through regulatory mechanisms affecting gene expression and are comprised of alleles that are both ancestry-specific and shared across continental ancestries.

Opportunities and challenges for transcriptome-wide association studies.

Transcriptome-wide association studies (TWAS) integrate genome-wide association studies (GWAS) and gene expression datasets to identify gene-trait associations. In this Perspective, we explore properties of TWAS as a potential approach to prioritize causal genes at GWAS loci, by using simulations and case studies of literature-curated candidate causal genes for schizophrenia, low-density-lipoprotein cholesterol and Crohn's disease. We explore risk loci where TWAS accurately prioritizes the likely causal gene as well as loci where TWAS prioritizes multiple genes, some likely to be non-causal, owing to sharing of expression quantitative trait loci (eQTL). TWAS is especially prone to spurious prioritization with expression data from non-trait-related tissues or cell types, owing to substantial cross-cell-type variation in expression levels and eQTL strengths. Nonetheless, TWAS prioritizes candidate causal genes more accurately than simple baselines. We suggest best practices for causal-gene prioritization with TWAS and discuss future opportunities for improvement. Our results showcase the strengths and limitations of using eQTL datasets to determine causal genes at GWAS loci.

Genetic architecture of gene expression traits across diverse populations.

For many complex traits, gene regulation is likely to play a crucial mechanistic role. How the genetic architectures of complex traits vary between populations and subsequent effects on genetic prediction are not well understood, in part due to the historical paucity of GWAS in populations of non-European ancestry. We used data from the MESA (Multi-Ethnic Study of Atherosclerosis) cohort to characterize the genetic architecture of gene expression within and between diverse populations. Genotype and monocyte gene expression were available in individuals with African American (AFA, n = 233), Hispanic (HIS, n = 352), and European (CAU, n = 578) ancestry. We performed expression quantitative trait loci (eQTL) mapping in each population and show genetic correlation of gene expression depends on shared ancestry proportions. Using elastic net modeling with cross validation to optimize genotypic predictors of gene expression in each population, we show the genetic architecture of gene expression for most predictable genes is sparse. We found the best predicted gene in each population, TACSTD2 in AFA and CHURC1 in CAU and HIS, had similar prediction performance across populations with R2 > 0.8 in each population. However, we identified a subset of genes that are well-predicted in one population, but poorly predicted in another. We show these differences in predictive performance are due to allele frequency differences between populations. Using genotype weights trained in MESA to predict gene expression in independent populations showed that a training set with ancestry similar to the test set is better at predicting gene expression in test populations, demonstrating an urgent need for diverse population sampling in genomics. Our predictive models and performance statistics in diverse cohorts are made publicly available for use in transcriptome mapping methods at https://github.com/WheelerLab/DivPop.

A Transcriptome-Wide Association Study Among 97,898 Women to Identify Candidate Susceptibility Genes for Epithelial Ovarian Cancer Risk.

Large-scale genome-wide association studies (GWAS) have identified approximately 35 loci associated with epithelial ovarian cancer (EOC) risk. The majority of GWAS-identified disease susceptibility variants are located in noncoding regions, and causal genes underlying these associations remain largely unknown. Here, we performed a transcriptome-wide association study to search for novel genetic loci and plausible causal genes at known GWAS loci. We used RNA sequencing data (68 normal ovarian tissue samples from 68 individuals and 6,124 cross-tissue samples from 369 individuals) and high-density genotyping data from European descendants of the Genotype-Tissue Expression (GTEx V6) project to build ovarian and cross-tissue models of genetically regulated expression using elastic net methods. We evaluated 17,121 genes for their cis-predicted gene expression in relation to EOC risk using summary statistics data from GWAS of 97,898 women, including 29,396 EOC cases. With a Bonferroni-corrected significance level of P < 2.2 × 10-6, we identified 35 genes, including FZD4 at 11q14.2 (Z = 5.08, P = 3.83 × 10-7, the cross-tissue model; 1 Mb away from any GWAS-identified EOC risk variant), a potential novel locus for EOC risk. All other 34 significantly associated genes were located within 1 Mb of known GWAS-identified loci, including 23 genes at 6 loci not previously linked to EOC risk. Upon conditioning on nearby known EOC GWAS-identified variants, the associations for 31 genes disappeared and three genes remained (P < 1.47 × 10-3). These data identify one novel locus (FZD4) and 34 genes at 13 known EOC risk loci associated with EOC risk, providing new insights into EOC carcinogenesis.Significance: Transcriptomic analysis of a large cohort confirms earlier GWAS loci and reveals FZD4 as a novel locus associated with EOC risk. Cancer Res; 78(18); 5419-30. ©2018 AACR.

Erratum to: A multi-stage genome-wide association study of uterine fibroids in African Americans.

The article "A multi-stage genome-wide association study of uterine fibroids in African Americans", written by Jacklyn N. Hellwege, was originally published Online First without open access. After publication in volume 136, issue 10, page 1363-1373 the author decided to opt for Open Choice and to make the article an open access publication. Therefore, the copyright of the article has been changed to

Trans-ethnic predicted expression genome-wide association analysis identifies a gene for estrogen receptor-negative breast cancer.

Genome-wide association studies (GWAS) have identified more than 90 susceptibility loci for breast cancer, but the underlying biology of those associations needs to be further elucidated. More genetic factors for breast cancer are yet to be identified but sample size constraints preclude the identification of individual genetic variants with weak effects using traditional GWAS methods. To address this challenge, we utilized a gene-level expression-based method, implemented in the MetaXcan software, to predict gene expression levels for 11,536 genes using expression quantitative trait loci and examine the genetically-predicted expression of specific genes for association with overall breast cancer risk and estrogen receptor (ER)-negative breast cancer risk. Using GWAS datasets from a Challenge launched by National Cancer Institute, we identified TP53INP2 (tumor protein p53-inducible nuclear protein 2) at 20q11.22 to be significantly associated with ER-negative breast cancer (Z = -5.013, p = 5.35×10-7, Bonferroni threshold = 4.33×10-6). The association was consistent across four GWAS datasets, representing European, African and Asian ancestry populations. There are 6 single nucleotide polymorphisms (SNPs) included in the prediction of TP53INP2 expression and five of them were associated with estrogen-receptor negative breast cancer, although none of the SNP-level associations reached genome-wide significance. We conducted a replication study using a dataset outside of the Challenge, and found the association between TP53INP2 and ER-negative breast cancer was significant (p = 5.07x10-3). Expression of HP (16q22.2) showed a suggestive association with ER-negative breast cancer in the discovery phase (Z = 4.30, p = 1.70x10-5) although the association was not significant after Bonferroni adjustment. Of the 249 genes that are 250 kb within known breast cancer susceptibility loci identified from previous GWAS, 20 genes (8.0%) were statistically significant associated with ER-negative breast cancer (p<0.05), compared to 582 (5.2%) of 11,287 genes that are not close to previous GWAS loci. This study demonstrated that expression-based gene mapping is a promising approach for identifying cancer susceptibility genes.

A multi-stage genome-wide association study of uterine fibroids in African Americans.

Uterine fibroids are benign tumors of the uterus affecting up to 77% of women by menopause. They are the leading indication for hysterectomy, and account for $34 billion annually in the United States. Race/ethnicity and age are the strongest known risk factors. African American (AA) women have higher prevalence, earlier onset, and larger and more numerous fibroids than European American women. We conducted a multi-stage genome-wide association study (GWAS) of fibroid risk among AA women followed by in silico genetically predicted gene expression profiling of top hits. In Stage 1, cases and controls were confirmed by pelvic imaging, genotyped and imputed to 1000 Genomes. Stage 2 used self-reported fibroid and GWAS data from 23andMe, Inc. and the Black Women's Health Study. Associations with fibroid risk were modeled using logistic regression adjusted for principal components, followed by meta-analysis of results. We observed a significant association among 3399 AA cases and 4764 AA controls at rs739187 (risk-allele frequency = 0.27) in CYTH4 (OR (95% confidence interval) = 1.23 (1.16-1.30), p value = 7.82 × 10-9). Evaluation of the genetic association results with MetaXcan identified lower predicted gene expression of CYTH4 in thyroid tissue as significantly associated with fibroid risk (p value = 5.86 × 10-8). In this first multi-stage GWAS for fibroids among AA women, we identified a novel risk locus for fibroids within CYTH4 that impacts gene expression in thyroid and has potential biological relevance for fibroids.

A Low-Frequency Inactivating AKT2 Variant Enriched in the Finnish Population Is Associated With Fasting Insulin Levels and Type 2 Diabetes Risk.

To identify novel coding association signals and facilitate characterization of mechanisms influencing glycemic traits and type 2 diabetes risk, we analyzed 109,215 variants derived from exome array genotyping together with an additional 390,225 variants from exome sequence in up to 39,339 normoglycemic individuals from five ancestry groups. We identified a novel association between the coding variant (p.Pro50Thr) in AKT2 and fasting plasma insulin (FI), a gene in which rare fully penetrant mutations are causal for monogenic glycemic disorders. The low-frequency allele is associated with a 12% increase in FI levels. This variant is present at 1.1% frequency in Finns but virtually absent in individuals from other ancestries. Carriers of the FI-increasing allele had increased 2-h insulin values, decreased insulin sensitivity, and increased risk of type 2 diabetes (odds ratio 1.05). In cellular studies, the AKT2-Thr50 protein exhibited a partial loss of function. We extend the allelic spectrum for coding variants in AKT2 associated with disorders of glucose homeostasis and demonstrate bidirectional effects of variants within the pleckstrin homology domain of AKT2.

C. elegans and mutants with chronic nicotine exposure as a novel model of cancer phenotype.

We previously investigated MET and its oncogenic mutants relevant to lung cancer in C. elegans. The inactive orthlogues of the receptor tyrosine kinase Eph and MET, namely vab-1 and RB2088 respectively, the temperature sensitive constitutively active form of KRAS, SD551 (let-60; GA89) and the inactive c-CBL equivalent mutants in sli-1 (PS2728, PS1258, and MT13032) when subjected to chronic exposure of nicotine resulted in a significant loss in egg-laying capacity and fertility. While the vab-1 mutant revealed increased circular motion in response to nicotine, the other mutant strains failed to show any effect. Overall locomotion speed increased with increasing nicotine concentration in all tested mutant strains except in the vab-1 mutants. Moreover, chronic nicotine exposure, in general, upregulated kinases and phosphatases. Taken together, these studies provide evidence in support of C. elegans as initial in vivo model to study nicotine and its effects on oncogenic mutations identified in humans.

MicroRNA biogenesis and cellular proliferation.

Given the fundamental roles of microRNAs (miRNAs) in physiological, developmental, and pathologic processes, we hypothesized that genes involved in miRNA biogenesis contribute to human complex traits. For 13 such genes, we evaluated the relationship between transcription and 2 classes of complex traits, namely cellular growth and sensitivity to various chemotherapeutic agents in a set of lymphoblastoid cell lines. We found a highly significant correlation between argonaute RNA-induced silencing complex catalytic component 2 (AGO2) expression and cellular growth rate (Bonferroni-adjusted P < 0.05), and report additional miRNA biogenesis genes with suggestive associations with either cellular growth rate or chemotherapeutic sensitivity. AGO2 expression was found to be correlated with multiple drug sensitivity phenotypes. Furthermore, small interfering RNA knockdown of AGO2 resulted in cellular growth inhibition in an ovarian cancer cell line (OVCAR-3), supporting the role of this miRNA biogenesis gene in cell proliferation in cancer cells. Expression quantitative trait loci mapping indicated that genetic variation (in the form of both single-nucleotide polymorphisms and copy number variations) that may regulate the expression of AGO2 can have downstream effects on cellular growth-dependent complex phenotypes.

Integrative analyses of genetic variation, epigenetic regulation, and the transcriptome to elucidate the biology of platinum sensitivity.

BACKGROUND: Using genome-wide genetic, gene expression, and microRNA expression (miRNA) data, we developed an integrative approach to investigate the genetic and epigenetic basis of chemotherapeutic sensitivity. RESULTS: Through a sequential multi-stage framework, we identified genes and miRNAs whose expression correlated with platinum sensitivity, mapped these to genomic loci as quantitative trait loci (QTLs), and evaluated the associations between these QTLs and platinum sensitivity. A permutation analysis showed that top findings from our approach have a much lower false discovery rate compared to those from a traditional GWAS of drug sensitivity. Our approach identified five SNPs associated with 10 miRNAs and the expression level of 15 genes, all of which were associated with carboplatin sensitivity. Of particular interest was one SNP (rs11138019), which was associated with the expression of both miR-30d and the gene ABCD2, which were themselves correlated with both carboplatin and cisplatin drug-specific phenotype in the HapMap samples. Functional study found that knocking down ABCD2 in vitro led to increased apoptosis in ovarian cancer cell line SKOV3 after cisplatin treatment. Over-expression of miR-30d in vitro caused a decrease in ABCD2 expression, suggesting a functional relationship between the two. CONCLUSIONS: We developed an integrative approach to the investigation of the genetic and epigenetic basis of human complex traits. Our approach outperformed standard GWAS and provided hints at potential biological function. The relationships between ABCD2 and miR-30d, and ABCD2 and platin sensitivity were experimentally validated, suggesting a functional role of ABCD2 and miR-30d in sensitivity to platinating agents.

The impact of microRNA expression on cellular proliferation.

As an important class of non-coding regulatory RNAs, microRNAs (miRNAs) play a key role in a range of biological processes. These molecules serve as post-transcriptional regulators of gene expression and their regulatory activity has been implicated in disease pathophysiology and pharmacological traits. We sought to investigate the impact of miRNAs on cellular proliferation to gain insight into the molecular basis of complex traits that depend on cellular growth, including, most prominently, cancer. We examined the relationship between miRNA expression and intrinsic cellular growth (iGrowth) in the HapMap lymphoblastoid cell lines derived from individuals of different ethnic backgrounds. We found a substantial enrichment for miRNAs (53 miRNAs, FDR < 0.05) correlated with cellular proliferation in pooled CEU (Caucasian of northern and western European descent) and YRI (individuals from Ibadan, Nigeria) samples. Specifically, 119 miRNAs (59 %) were significantly correlated with iGrowth in YRI; of these miRNAs, 18 were correlated with iGrowth in CEU. To gain further insight into the effect of miRNAs on cellular proliferation in cancer, we showed that over-expression of miR-22, one of the top iGrowth-associated miRNAs, leads to growth inhibition in an ovarian cancer cell line (SKOV3). Furthermore, over-expression of miR-22 down-regulates the expression of its target genes (MXI1 and SLC25A37) in this ovarian cancer cell line, highlighting an miRNA-mediated regulatory network potentially important for cellular proliferation. Importantly, our study identified miRNAs that can be used as molecular targets in cancer therapy.

Association of body mass index and outcomes after major lung resection.

OBJECTIVES: Obesity has been thought to predispose patients to excess morbidity after lung resection because of decreased diaphragm excursion, reduced lung volumes and relative immobility. We assessed the relationship of body mass index (BMI) to acute outcomes after major lung resection. METHODS: Information from our database of lung resections was evaluated for the period 1980-2011. Univariate analysis for adverse events (pulmonary, cardiovascular, other and overall) was used to select variables for inclusion in multivariate logistic regression analyses. Missing values were imputed. BMI was categorized as underweight (<18.5), normal (18.5-24.9), overweight (25-29.9), obese (30-34.9) and very obese (≥ 35). RESULTS: Among 1369 patients, there were 703 males (51%) and the mean age was 62 ± 11 years. Complications included the following: pulmonary 12%, cardiovascular 15%, other 16%, mortality 5% and any 29%. The incidence of complications decreased during each decade of study (40, 30, 26, 20%; P < 0.0001) and the incidence of obese/very obese increased during the same intervals (11, 22, 30, 25%; P = 0.0007). Adjusting for age, performance status, coronary artery disease, smoking status, diffusing capacity of the lung for carbon monoxide, forced expiratory volume in 1 s and operation year, being overweight/obese/very obese did not increase the risk of postoperative complications in any category. In fact, patients in this group showed a lower rate of cardiovascular complications than those with BMI ≤ 25 (odds ratio (OR): 0.72; 95% confidence interval (CI): 0.51-1.00; P = 0.048). However, being underweight was importantly associated with an increased risk of pulmonary complications (OR: 2.5; 95% CI: 1.3-4.9; P = 0.0087) and of operative mortality (OR: 2.96; 95% CI: 1.28-6.86; P = 0.011). CONCLUSION: Being overweight or obese does not increase the risk of complications after major lung resection. In contrast, patients who are underweight are at significantly increased risk of pulmonary complications and mortality. Knowledge of the relationship of BMI to perioperative risk for major lung resection is essential in proper risk stratification.

Genetic and epigenetic variants contributing to clofarabine cytotoxicity.

2-chloro-2-fluoro-deoxy-9-D-arabinofuranosyladenine (Clofarabine), a purine nucleoside analog, is used in the treatment of hematologic malignancies and as induction therapy for stem cell transplantation. The discovery of pharmacogenomic markers associated with chemotherapeutic efficacy and toxicity would greatly benefit the utility of this drug. Our objective was to identify genetic and epigenetic variants associated with clofarabine toxicity using an unbiased, whole genome approach. To this end, we employed International HapMap lymphoblastoid cell lines (190 LCLs) of European (CEU) or African (YRI) ancestry with known genetic information to evaluate cellular sensitivity to clofarabine. We measured modified cytosine levels to ascertain the contribution of genetic and epigenetic factors influencing clofarabine-mediated cytotoxicity. Association studies revealed 182 single nucleotide polymorphisms (SNPs) and 143 modified cytosines associated with cytotoxicity in both populations at the threshold P ≤ 0.0001. Correlation between cytotoxicity and baseline gene expression revealed 234 genes at P ≤ 3.98 × 10(-6). Six genes were implicated as: (i) their expression was directly correlated to cytotoxicity, (ii) they had a targeting SNP associated with cytotoxicity, and (iii) they had local modified cytosines associated with gene expression and cytotoxicity. We identified a set of three SNPs and three CpG sites targeting these six genes explaining 43.1% of the observed variation in phenotype. siRNA knockdown of the top three genes (SETBP1, BAG3, KLHL6) in LCLs revealed altered susceptibility to clofarabine, confirming relevance. As clofarabine's toxicity profile includes acute kidney injury, we examined the effect of siRNA knockdown in HEK293 cells. siSETBP1 led to a significant change in HEK293 cell susceptibility to clofarabine.

Trans-population analysis of genetic mechanisms of ethnic disparities in neuroblastoma survival.

BACKGROUND: Black patients with neuroblastoma have a higher prevalence of high-risk disease and worse outcome than white patients. We sought to investigate the relationship between genetic variation and the disparities in survival observed in neuroblastoma. METHODS: The analytic cohort was composed of 2709 patients. Principal components were used to assign patients to genomic ethnic clusters for survival analyses. Locus-specific ancestry was calculated for use in association analysis. The shorter spans of linkage disequilibrium in African populations may facilitate the fine mapping of causal variants in regions previously implicated by genome-wide association studies conducted primarily in patients of European descent. Thus, we evaluated 13 single nucleotide polymorphisms known to be associated with susceptibility to high-risk neuroblastoma from genome-wide association studies and all variants with highly divergent allele frequencies in reference African and European populations near the known susceptibility loci. All statistical tests were two-sided. RESULTS: African genomic ancestry was associated with high-risk neuroblastoma (P = .007) and lower event-free survival (P = .04, hazard ratio = 1.4, 95% confidence interval = 1.05 to 1.80). rs1033069 within SPAG16 (sperm associated antigen 16) was determined to have higher risk allele frequency in the African reference population and statistically significant association with high-risk disease in patients of European and African ancestry (P = 6.42 × 10(-5), false discovery rate < 0.0015) in the overall cohort. Multivariable analysis using an additive model demonstrated that the SPAG16 single nucleotide polymorphism contributes to the observed ethnic disparities in high-risk disease and survival. CONCLUSIONS: Our study demonstrates that common genetic variation influences neuroblastoma phenotype and contributes to the ethnic disparities in survival observed and illustrates the value of trans-population mapping.

Identification of novel germline polymorphisms governing capecitabine sensitivity.

BACKGROUND: Capecitabine, an oral 5-fluorouracil (5-FU) prodrug, is widely used in the treatment of breast, colorectal, and gastric cancers. To guide the selection of patients with potentially the greatest benefit of experiencing antitumor efficacy, or, alternatively, of developing toxicities, identifying genomic predictors of capecitabine sensitivity could permit its more informed use. METHODS: The objective of this study was to perform capecitabine sensitivity genome-wide association studies (GWAS) using 503 well genotyped human cell lines from individuals representing multiple different world populations. A meta-analysis that included all ethnic populations then enabled the identification of novel germline determinants (single nucleotide polymorphisms [SNPs]) of capecitabine susceptibility. RESULTS: First, an intrapopulation GWAS of Caucasian individuals identified reference SNP 4702484 (rs4702484) (within adenylate cyclase 2 [ADCY2]) at a level reaching genome-wide significance (P = 5.2 × 10(-8) ). This SNP is located upstream of the 5 methyltetrahydrofolate-homocysteine methyltransferase reductase (MTRR) gene, and it is known that the enzyme for MTRR is involved in the methionine-folate biosynthesis and metabolism pathway, which is the primary target of 5-FU-related compounds, although the authors were unable to identify a direct relation between rs4702484 and MTRR expression in a tested subset of cells. In the meta-analysis, 4 SNPs comprised the top hits, which, again, included rs4702484 and 3 additional SNPs (rs8101143, rs576523, and rs361433) that approached genome-wide significance (P values from 1.9 × 10(-7) to 8.8 × 10(-7) ). The meta-analysis also identified 1 missense variant (rs11722476; serine to asparagine) within switch/sucrose nonfermentable-related, matrix-associated, actin-dependent regulator of chromatin (SMARCAD1), a novel gene for association with capecitabine/5-FU susceptibility. CONCLUSIONS: Toward the goal of individualizing cancer chemotherapy, the current study identified novel SNPs and genes associated with capecitabine sensitivity that are potentially informative and testable in any patient regardless of ethnicity.