Prognostic implications of genetic variants in advanced non-small cell lung cancer: a genome-wide association study

New insights into the physiology and pathophysiology of the atypical sodium leak channel NALCN

Cell excitability and its modulation by hormones and neurotransmitters involve the concerted action of a large repertoire of membrane proteins, especially ion channels. Unique complements of coexpressed ion channels are exquisitely balanced against each other in different excitable cell types, establishing distinct electrical properties that are tailored for diverse physiological contributions, and dysfunction of any component may induce a disease state. A crucial parameter controlling cell excitability is the resting membrane potential (RMP) set by extra- and intracellular concentrations of ions, mainly Na+, K+, and Cl-, and their passive permeation across the cell membrane through leak ion channels. Indeed, dysregulation of RMP causes significant effects on cellular excitability. This review describes the molecular and physiological properties of the Na+ leak channel NALCN, which associates with its accessory subunits UNC-79, UNC-80, and NLF-1/FAM155 to conduct depolarizing background Na+ currents in various excitable cell types, especially neurons. Studies of animal models clearly demonstrate that NALCN contributes to fundamental physiological processes in the nervous system including the control of respiratory rhythm, circadian rhythm, sleep, and locomotor behavior. Furthermore, dysfunction of NALCN and its subunits is associated with severe pathological states in humans. The critical involvement of NALCN in physiology is now well established, but its study has been hampered by the lack of specific drugs that can block or agonize NALCN currents in vitro and in vivo. Molecular tools and animal models are now available to accelerate our understanding of how NALCN contributes to key physiological functions and the development of novel therapies for NALCN channelopathies.

Genome-wide mapping of gene-microbe interactions in the murine lung microbiota based on quantitative microbial profiling

BACKGROUND

Mammalian lungs comprise a complex microbial ecosystem that interacts with host physiology. Previous research demonstrates that the environment significantly contributes to bacterial community structure in the upper and lower respiratory tract. However, the influence of host genetics on the makeup of lung microbiota remains ambiguous, largely due to technical difficulties related to sampling, as well as challenges inherent to investigating low biomass communities. Thus, innovative approaches are warranted to clarify host-microbe interactions in the mammalian lung.

RESULTS

Here, we aimed to characterize host genomic regions associated with lung bacterial traits in an advanced intercross mouse line (AIL). By performing quantitative microbial profiling (QMP) using the highly precise method of droplet digital PCR (ddPCR), we refined 16S rRNA gene amplicon-based traits to identify and map candidate lung-resident taxa using a QTL mapping approach. In addition, the two abundant core taxa Lactobacillus and Pelomonas were chosen for independent microbial phenotyping using genus-specific primers. In total, this revealed seven significant loci involving eight bacterial traits. The narrow confidence intervals afforded by the AIL population allowed us to identify several promising candidate genes related to immune and inflammatory responses, cell apoptosis, DNA repair, and lung functioning and disease susceptibility. Interestingly, one genomic region associated with Lactobacillus abundance contains the well-known anti-inflammatory cytokine Il10, which we confirmed through the analysis of Il10 knockout mice.

CONCLUSIONS

Our study provides the first evidence for a role of host genetic variation contributing to variation in the lung microbiota. This was in large part made possible through the careful curation of 16S rRNA gene amplicon data and the incorporation of a QMP-based methods. This approach to evaluating the low biomass lung environment opens new avenues for advancing lung microbiome research using animal models.

Machine learning algorithm-generated and multi-center validated melanoma prognostic signature with inspiration for treatment management

BACKGROUND

Although immunotherapy and targeted treatments have dramatically improved the survival of melanoma patients, the intra- or intertumoral heterogeneity and drug resistance have hindered the further expansion of clinical benefits.

METHODS

The 96 combination frames constructed by ten machine learning algorithms identified a prognostic consensus signature based on 1002 melanoma samples from nine independent cohorts. Clinical features and 26 published signatures were employed to compare the predictive performance of our model.

RESULTS

A machine learning-based prognostic signature (MLPS) with the highest average C-index was developed via 96 algorithm combinations. The MLPS has a stable and excellent predictive performance for overall survival, superior to common clinical traits and 26 collected signatures. The low MLPS group with a better prognosis had significantly enriched immune-related pathways, tending to be an immune-hot phenotype and possessing potential immunotherapeutic responses to anti-PD-1, anti-CTLA-4, and MAGE-A3. On the contrary, the high MLPS group with more complex genomic alterations and poorer prognoses is more sensitive to the BRAF inhibitor dabrafenib, confirmed in patients with BRAF mutations.

CONCLUSION

MLPS could independently and stably predict the prognosis of melanoma, considered a promising biomarker to identify patients suitable for immunotherapy and those with BRAF mutations who would benefit from dabrafenib.

Identification of THSD7B and PRMT9 mutations as risk factors for familial lung adenocarcinoma: A case report

RATIONALE

Lung tumors arise from the unrestrained malignant growth of pulmonary epithelial cells. Lung cancer cases include both small and non-small cell lung cancers, with lung adenocarcinoma (LUAD) accounting for roughly half of all non-small cell lung cancer cases. Research focused on familial cancers suggests that approximately 8% of lung cancer cases are linked to genetic susceptibility or heritability. The precise genetic factors that underlie the onset of lung cancer, however, remain to be firmly established.

PATIENT CONCERNS

A 43-year-old presented with nodules in the lower left lung lobe. Following initial antibiotic treatment in a local hospital, these nodules remained present and the patient subsequently underwent the resection of the left lower lobe of the lung. The patient also had 4 family members with a history of LUAD.

DIAGNOSIS

Immunohistochemical staining results including cytokeratin 7 (+), TTF-1 (+), new aspartic proteinase A (+), CK5/6 (-), P63 (-), and Ki-67 (5%+) were consistent with a diagnosis of LUAD.

INTERVENTION

Whole exome sequencing analyses of 5 patients and 6 healthy family members were performed to explore potential mutations associated with familial LUAD.

OUTCOMES

Whole exome sequencing was conducted, confirming that the proband and their 4 other family members with LUAD harbored heterozygous THSD7B (c.A4000G:p.S1334G) mutations and homozygous PRMT9 (c.G40T:p.G14C) mutations, as further confirmed via Sanger sequencing. These mutations were predicted to be deleterious using the SIFT, PolyPhen2, and MutationTaster algorithms. Protein structure analyses indicated that the mutation of the serine at amino acid position 1334 in THSD7B to a glycine would reduce the minimum free energy from 8.08 kcal/mol to 68.57 kcal/mol. The identified mutation in the PRMT9 mutation was not present in the predicted protein structure. I-Mutant2.0 predictions indicated that both of these mutations (THSD7B:p.S1334G and PRMT9: p.G14C) were predicted to reduce protein stability.

LESSONS

Heterozygous THSD7B (c.A4000G:p.S1334G) and the homozygous PRMT9 (c.G40T:p.G14C) mutations were found to be linked to LUAD incidence in the analyzed family. Early analyses of these genetic loci and timely genetic counseling may provide benefits and aid in the early diagnosis of familial LUAD.

NALCN is a potential biomarker and therapeutic target in human cancers

Background: Sodium leak channel non-selective (NALCN), known as a voltage-independent Na⁺ channel, is increasingly considered to play vital roles in tumorigenesis and metastasis of human cancers. However, no comprehensive pan-cancer analysis of NALCN has been conducted. Our study aims to explore the potential diagnostic, prognostic and therapeutic value of NALCN in human cancers. Methods: Through comprehensive application of datasets from Human Protein Atlas (HPA), The Cancer Genome Atlas (TCGA), Cancer Cell Line Encyclopedia (CCLE), Enhanced Version of Tumor Immune Estimation Resource (TIMER2.0), Tumor and Immune System Interaction Database (TISIDB), The University of Alabama at Birmingham Cancer data analysis Portal (UALCAN), cBioPortal, GeneMANIA and Search Tool for the Retrieval of Interaction Gene/Proteins (STRING) databases, we explored the potential roles of NALCN in different cancers. The differential expression, prognostic implications, pathological stages and grades, molecular and immune subtypes, diagnostic accuracy, tumor mutation burden (TMB), microsatellite instability (MSI), mismatch repair (MMR) genes, immune checkpoint genes, chemokine genes, major histocompatibility complex (MHC)-related genes, tumor-infiltrating immune cells (TIICs), promoter methylation, mutations, copy number alteration (CNA), and functional enrichment related to NALCN were analyzed. Results: Most cancers lowly expressed NALCN. Upregulated NALCN expression was associated with poor or better prognosis in different cancers. Moreover, NALCN was correlated with clinicopathological features in multiple cancers. NALCN showed high diagnostic accuracy in 5 caner types. NALCN is highly linked with immune-related biomarkers, immune-related genes and TIICs. Significant methylation changes and genetic alteration of NALCN can be observed in many cancers. Enrichment analysis showed that NALCN is closely related to multiple tumor-related signaling pathways. Conclusion: Our study revealed the vital involvement of NALCN in cancer. NALCN can be used as a prognostic biomarker for immune infiltration and clinical outcomes, and has potential diagnostic and therapeutic implications.

The Role of the Adapter Protein Anks1a in the Regulation of Breast Cancer Cell Motility

Epithelial-mesenchymal transition (EMT) is a critical step in tumor progression that leads to the acquisition by cancer cells the capacity for migration using the mesenchymal motility mode regulated by the Rac→WAVE→Arp2/3 signaling pathway. Earlier it was shown that proteins interacting with Rac can regulate mesenchymal migration and thus determine the metastatic potential of the cells. The search for new regulators of cell migration is an important theoretical and practical task. The adaptor protein Anks1a is one of the proteins interacting with Rac, whose expression is altered in many types of tumors. The aim of this study was to find whether Anks1a affects the migration of cancer cells and to identify the mechanism underlying this effect. It was suggested that Anks1a can influence cancer cell migration either as a Rac1 effector or by activating human epidermal growth factor receptor 2 (HER2) exchange. We investigated how upregulation and inhibition of Anks1a expression affected migration of breast cancer cells with different HER2 status. Anks1a was shown to interact with the activated form of Rac1. In the MDA-MB-231 cells (triple negative cancer), which lack HER2, Anks1a accumulated at the active cell edge, which is characterized by enrichment with active Rac1, whereas no such accumulation was observed in the HER2-overexpressing SK-BR-3 cells. Downregulation of the ANKS1a expression with esiRNA had almost no effect on the cancer cell motility, except a slight increase in the average migration rate of MDA-MB-231 cells. Among three cell lines tested, overexpression of Anks1a increased the migration rate of HER2-overexpressng SK-BR-3 cells only. We showed that Anks1a is an effector of activated Rac1, but its influence on the cell migration in this capacity was minimal, at least in the studied breast cancer cells. Anks1a affected the motility of breast cancer cells due to its involvement in the EGF receptor exchange.

Genetic predisposition to lung adenocarcinoma outcome is a feature already present in patients' noninvolved lung tissue

Emerging evidence suggests that the prognosis of patients with lung adenocarcinoma can be determined from germline variants and transcript levels in nontumoral lung tissue. Gene expression data from noninvolved lung tissue of 483 lung adenocarcinoma patients were tested for correlation with overall survival using multivariable Cox proportional hazard and multivariate machine learning models. For genes whose transcript levels are associated with survival, we used genotype data from 414 patients to identify germline variants acting as cis-expression quantitative trait loci (eQTLs). Associations of eQTL variant genotypes with gene expression and survival were tested. Levels of four transcripts were inversely associated with survival by Cox analysis (CLCF1, hazard ratio [HR] = 1.53; CNTNAP1, HR = 2.17; DUSP14, HR = 1.78; and MT1F: HR = 1.40). Machine learning analysis identified a signature of transcripts associated with lung adenocarcinoma outcome that was largely overlapping with the transcripts identified by Cox analysis, including the three most significant genes (CLCF1, CNTNAP1, and DUSP14). Pathway analysis indicated that the signature is enriched for ECM components. We identified 32 cis-eQTLs for CNTNAP1, including 6 with an inverse correlation and 26 with a direct correlation between the number of minor alleles and transcript levels. Of these, all but one were prognostic: the six with an inverse correlation were associated with better prognosis (HR < 1) while the others were associated with worse prognosis. Our findings provide supportive evidence that genetic predisposition to lung adenocarcinoma outcome is a feature already present in patients' noninvolved lung tissue.

Exploring the Genomic Patterns in Human and Mouse Cerebellums Via Single-Cell Sequencing and Machine Learning Method

In mammals, the cerebellum plays an important role in movement control. Cellular research reveals that the cerebellum involves a variety of sub-cell types, including Golgi, granule, interneuron, and unipolar brush cells. The functional characteristics of cerebellar cells exhibit considerable differences among diverse mammalian species, reflecting a potential development and evolution of nervous system. In this study, we aimed to recognize the transcriptional differences between human and mouse cerebellum in four cerebellar sub-cell types by using single-cell sequencing data and machine learning methods. A total of 321,387 single-cell sequencing data were used. The 321,387 cells included 4 cell types, i.e., Golgi (5,048, 1.57%), granule (250,307, 77.88%), interneuron (60,526, 18.83%), and unipolar brush (5,506, 1.72%) cells. Our results showed that by using gene expression profiles as features, the optimal classification model could achieve very high even perfect performance for Golgi, granule, interneuron, and unipolar brush cells, respectively, suggesting a remarkable difference between the genomic profiles of human and mouse. Furthermore, a group of related genes and rules contributing to the classification was identified, which might provide helpful information for deepening the understanding of cerebellar cell heterogeneity and evolution.

Genome-wide association study of hospitalized COVID-19 patients in the United Arab Emirates

BACKGROUND

The heterogeneity in symptomatology and phenotypic profile attributable to COVID-19 is widely unknown. The objective of this manuscript is to conduct a trans-ancestry genome wide association study (GWAS) meta-analysis of COVID-19 severity to improve the understanding of potentially causal targets for SARS-CoV-2.

METHODS

This cross-sectional study recruited 646 participants in the UAE that were divided into two phenotypic groups based on the severity of COVID-19 phenotypes, hospitalized (n=482) and non-hospitalized (n=164) participants. Hospitalized participants were COVID-19 patients that developed acute respiratory distress syndrome (ARDS), pneumonia or progression to respiratory failure that required supplemental oxygen therapy or mechanical ventilation support or had severe complications such as septic shock or multi-organ failure. We conducted a trans-ancestry meta-analysis GWAS of European (n=302), American (n=102), South Asian (n=99), and East Asian (n=107) ancestry populations. We also carried out comprehensive post-GWAS analysis, including enrichment of SNP associations in tissues and cell-types, expression quantitative trait loci and differential expression analysis.

FINDINGS

Eight genes demonstrated a strong association signal: VWA8 gene in locus 13p14·11 (SNP rs10507497; p=9·54 x10-7), PDE8B gene in locus 5q13·3 (SNP rs7715119; p=2·19 x10-6), CTSC gene in locus 11q14·2 (rs72953026; p=2·38 x10-6), THSD7B gene in locus 2q22·1 (rs7605851; p=3·07x10-6), STK39 gene in locus 2q24·3 (rs7595310; p=4·55 x10-6), FBXO34 gene in locus 14q22·3 (rs10140801; p=8·26 x10-6), RPL6P27 gene in locus 18p11·31 (rs11659676; p=8·88 x10-6), and METTL21C gene in locus 13q33·1 (rs599976; p=8·95 x10-6). The genes are expressed in the lung, associated to tumour progression, emphysema, airway obstruction, and surface tension within the lung, as well as an association to T-cell-mediated inflammation and the production of inflammatory cytokines.

INTERPRETATION

We have discovered eight highly plausible genetic association with hospitalized cases in COVID-19. Further studies must be conducted on worldwide population genetics to facilitate the development of population specific therapeutics to mitigate this worldwide challenge.

FUNDING

This review was commissioned as part of a project to study the host cell receptors of coronaviruses funded by Khalifa University's CPRA grant (Reference number 2020-004).

Functional annotation of lung cancer‒associated genetic variants by cell type‒specific epigenome and long-range chromatin interactome

Functional interpretation of noncoding genetic variants associated with complex human diseases and traits remains a challenge. In an effort to enhance our understanding of common germline variants associated with lung cancer, we categorize regulatory elements based on eight major cell types of human lung tissue. Our results show that 21.68% of lung cancer‒associated risk variants are linked to noncoding regulatory elements, nearly half of which are cell type‒specific. Integrative analysis of high-resolution long-range chromatin interactome maps and single-cell RNA-sequencing data of lung tumors uncovers number of putative target genes of these variants and functionally relevant cell types, which display a potential biological link to cancer susceptibility. The present study greatly expands the scope of functional annotation of lung cancer‒associated genetic risk factors and dictates probable cell types involved in lung carcinogenesis.

Genetic Determinants of Lung Cancer Prognosis in Never Smokers: A Pooled Analysis in the International Lung Cancer Consortium

BACKGROUND

Lung cancer remains the leading cause of cancer death worldwide, with 15% to 20% occurring in never smokers. To assess genetic determinants for prognosis among never smokers, we conducted a genome-wide investigation in the International Lung Cancer Consortium (ILCCO).

METHODS

Genomic and clinical data from 1,569 never-smoking patients with lung cancer of European ancestry from 10 ILCCO studies were included. HRs and 95% confidence intervals of overall survival were estimated. We assessed whether the associations were mediated through mRNA expression-based 1,553 normal lung tissues from the lung expression quantitative trait loci (eQTL) dataset and Genotype-Tissue Expression (GTEx). For cross-ethnicity generalization, we assessed the associations in a Japanese study (N = 887).

RESULTS

One locus at 13q22.2 was associated with lung adenocarcinoma survival at genome-wide level, with carriers of rs12875562-T allele exhibiting poor prognosis [HR = 1.71 (1.41-2.07), P = 3.60 × 10-8], and altered mRNA expression of LMO7DN in lung tissue (GTEx, P = 9.40 × 10-7; Lung eQTL dataset, P = 0.003). Furthermore, 2 of 11 independent loci that reached the suggestive significance level (P < 10-6) were significant eQTL affecting mRNA expression of nearby genes in lung tissues, including CAPZB at 1p36.13 and UBAC1 at 9q34.3. One locus encoding NWD2/KIAA1239 at 4p14 showed associations in both European [HR = 0.50 (0.38-0.66), P = 6.92 × 10-7] and Japanese populations [HR = 0.79 (0.67-0.94), P = 0.007].

CONCLUSIONS

Based on the largest genomic investigation on the lung cancer prognosis of never smokers to date, we observed that lung cancer prognosis is affected by inherited genetic variants.

IMPACT

We identified one locus near LMO7DN at genome-wide level and several potential prognostic genes with cis-effect on mRNA expression. Further functional genomics work is required to understand their role in tumor progression.

Integrated genomic analysis identifies a genetic mutation model predicting response to immune checkpoint inhibitors in melanoma

Several biomarkers such as tumor mutation burden (TMB), neoantigen load (NAL), programmed cell‐death receptor 1 ligand (PD‐L1) expression, and lactate dehydrogenase (LDH) have been developed for predicting response to immune checkpoint inhibitors (ICIs) in melanoma. However, some limitations including the undefined cut‐off value, poor uniformity of test platform, and weak reliability of prediction have restricted the broad application in clinical practice. In order to identify a clinically actionable biomarker and explore an effective strategy for prediction, we developed a genetic mutation model named as immunotherapy score (ITS) for predicting response to ICIs therapy in melanoma, based on whole‐exome sequencing data from previous studies. We observed that patients with high ITS had better durable clinical benefit and survival outcomes than patients with low ITS in three independent cohorts, as well as in the meta‐cohort. Notably, the prediction capability of ITS was more robust than that of TMB. Remarkably, ITS was not only an independent predictor of ICIs therapy, but also combined with TMB or LDH to better predict response to ICIs than any single biomarker. Moreover, patients with high ITS harbored the immunotherapy‐sensitive characteristics including high TMB and NAL, ultraviolet light damage, impaired DNA damage repair pathway, arrested cell cycle signaling, and frequent mutations in NF1 and SERPINB3/4. Overall, these findings deserve prospective investigation in the future and may help guide clinical decisions on ICIs therapy for patients with melanoma.

Association of PTEN Gene SNPs rs2299939 With PFS in Patients With Small Cell Lung Cancer Treated With Early Radiotherapy

Lung cancer has higher morbidity and mortality than most cancers. It is common that there are some phenomenons of secondary drug resistance, radiotherapy resistance and poor prognosis during the treatment of small cell lung cancer (SCLC). Recent studies revealed that the single-nucleotide polymorphisms (SNPs) are associated with the curative effect among patients with the same pathological type and stage. Our study analyzed the start time of radiotherapy and the relationship between PTEN gene rs2299939 polymorphisms and survival time among 116 SCLC patients. The results showed that early radiotherapy significantly improved the time of survival in patients compared with late radiotherapy (P = 0.029). Simultaneously, the study found that patients with the rs2299939 AA genotype showed significant sensitivity to both early and late radiotherapy, but early radiotherapy is better. The median survival time of CC genotype patients was 12 months in the early radiotherapy group while it was 9 months in the late radiotherapy group, thus recommending early radiotherapy among these patients. In addition, it was found that rs2299939 could regulate the expression of related genes in peripheral blood and lung tissues by eQTL analysis. This study revealed that the early radiotherapy could prolong the PFS of SCLC and shall be performed in SCLC treatment.

Hyponatremia and Cancer Progression: Possible Association with Sodium-Transporting Proteins

Hyponatremia, the phenomenon of serum sodium level falling below 135 mmol/L, is seen frequently in cancer patients and has been correlated with poor prognosis. Hyponatremia has classically been attributed to the "syndrome of inappropriate antidiuretic hormone secretion," leading to prolonged fluid retention. However, this is unlikely to be the only mechanism. In this study, we advance the hypothesis that upregulation of various sodium-transporting proteins during the cancer process makes a significant contribution to the pathophysiology of cancer-associated hyponatremia. Such sodium-transporting proteins include voltage-gated sodium channels, especially its hypoxia-promoted persistent current, epithelial sodium channels, and transient receptor potential channels. Thus, hyponatremia follows cancer, whereby drop in blood serum level occurs as a result of uptake of sodium from extracellular fluid by cancer cells. Indeed, the sodium content of cancer cells/tissues is higher than normal. In turn, the rise in the intracellular sodium concentration brings about a range of cellular effects, including extracellular acidification that promotes invasiveness and thus leads to poor prognosis. This perspective offers novel therapies for cancer and the associated hyponatremia.

Gene-environment and gene-gene interactions between CHRNA3 rs1051730, XRCC1 rs25487, and ERCC1 rs735482 variants highly elevate the risk of lung cancer

Background

Gene-gene and gene-environment interactions play an important role in cancer susceptibility. In this work, we studied the association of XRCC1 rs25487, ERCC1 rs735482, and CHRNA3 rs1051730 variants with lung cancer and assessed the modulatory effect of potential interaction between these variants on disease risk.

Results

In this study, 86 primary lung cancer patients and 64 control subjects were genotyped for CHRNA3 rs1051730, XRCC1 rs25487, and ERCC1 rs735482 by real-time PCR. The frequency of the three studied variants was higher among lung cancer patients than in control subjects, but with no statistical significance. ERCC1 rs735482 variant was associated with 6.9-fold increased risk to develop lung cancer among smokers (p = 0.03). Concomitant presence of CHRNA3 and ERCC1 wild alleles was associated with 2.7-fold elevated risk of lung cancer (p < 0.0001), while concomitant presence of CHRNA3 rs1051730 variant allele with ERCC1 wild allele was associated with 20-fold elevated risk (p < 0.000). Concomitant presence of both variants, ERCC1 rs735482 and CHRNA3 rs1051730, was associated with 9.9-fold elevated risk (p < 0.0001). Meanwhile, the concomitant presence of XRCC1 rs25487 with either ERCC1 rs735482 or CHRNA3 rs1051730 or both was not associated with increased risk of the disease.

Conclusion

Our results emphasize the role of gene-gene interaction in the pathogenesis of lung cancer. Large-scale further studies to clarify the underlying mechanisms are needed.

THUMP domain containing 2 protein possibly induces resistance to cisplatin and 5-fluorouracil in in vitro human esophageal squamous cell carcinoma cells as revealed by transposon activation mutagenesis

BACKGROUND

Although chemotherapy is a core treatment for esophageal cancer, some patients develop drug resistance. Gene screening with transposons (i.e. mobile genetic elements) is a novel procedure for identifying chemotherapy-resistant genes. Transposon insertion can randomly affect nearby gene expression. By identifying the affected genes, candidate genes can be found. The present study aimed to identify cisplatin (CDDP)/5-fluorouracil (5-FU)-resistant genes in in vitro human esophageal squamous cell carcinoma with transposons.

METHODS

After establishing transposon-tagged cells, we obtained CDDP/5-FU-resistant colonies. A polymerase chain reaction and sequencing were used to identify the transposon inserted site and candidate CDDP/5-FU resistant genes. Focusing on one candidate gene, we confirmed CDDP/5-FU resistance by comparing the IC₅₀ between drug-resistant and wild-type cells. Furthermore, we investigated gene expression by a real-time polymerase chain reaction. Finally, we mediated the candidate gene level with small interfering RNA to confirm the resistance.

RESULTS

Thirty-nine candidate genes for CDDP/5-FU resistance were identified. Nineteen were for CDDP resistance and 27 were for 5-FU resistance. Seven genes, THUMP domain-containing protein 2 (THUMPD2), nuclear factor interleukin-3-regulated protein (NFIL3), tyrosine-protein kinase transmembrane receptor 2 (ROR2), C-X-C chemokine receptor type 4 (CXCR4), thrombospondin type-1 domain-containing protein 2 (THSD7B) alpha-parvin (PARVA) and TEA domain transcription factor 1 (TEAD1), were detected as candidate genes in both colonies. Regarding THUMPD2, its expression was downregulated and knocking down THUMPD2 suggested drug resistance in both drugs.

CONCLUSIONS

Thirty-nine candidate genes were identified with transposons. The downregulation of THUMPD2 was suggested to play a role in multidrug resistance in in vitro esophageal squamous cell carcinoma.

Pipeline design to identify key features and classify the chemotherapy response on lung cancer patients using large-scale genetic data

BACKGROUND

During the last decade, the interest to apply machine learning algorithms to genomic data has increased in many bioinformatics applications. Analyzing this type of data entails difficulties for managing high-dimensional data, class imbalance for knowledge extraction, identifying important features and classifying individuals. In this study, we propose a general framework to tackle these challenges with different machine learning algorithms and techniques. We apply the configuration of this framework on lung cancer patients, identifying genetic signatures for classifying response to drug treatment response. We intersect these relevant SNPs with the GWAS Catalog of the National Human Genome Research Institute and explore the Regulomedb, GTEx databases for functional analysis purposes.

RESULTS

The machine learning based solution proposed in this study is a scalable and flexible alternative to the classical uni-variate regression approach to analyze large-scale data. From 36 experiments executed using the machine learning framework design, we obtain good classification performance from the top 5 models with the highest cross-validation score and the smallest standard deviation. One thousand two hundred twenty four SNPs corresponding to the key features from the top 20 models (cross validation F1 mean >= 0.65) were compared with the GWAS Catalog finding no intersection with genome-wide significant reported hits. From these, new genetic signatures in MAE, CEP104, PRKCZ and ADRB2 show relevant biological regulatory functionality related to lung physiology.

CONCLUSIONS

We have defined a machine learning framework using data with an unbalanced large data-set of SNP-arrays and imputed genotyping data from a pharmacogenomics study in lung cancer patients subjected to first-line platinum-based treatment. This approach found genome signals with no genome-wide significance in the uni-variate regression approach (GWAS Catalog) that are valuable for classifying patients, only few of them with related biological function. The effect results of these variants can be explained by the recently proposed omnigenic model hypothesis, which states that complex traits can be influenced mostly by genes outside not only by the "core genes", mainly found by the genome-wide significant SNPs, but also by the rest of genes outside of the "core pathways" with apparent unrelated biological functionality.

Ion Channels in Cancer: Are Cancer Hallmarks Oncochannelopathies?

Genomic instability is a primary cause and fundamental feature of human cancer. However, all cancer cell genotypes generally translate into several common pathophysiological features, often referred to as cancer hallmarks. Although nowadays the catalog of cancer hallmarks is quite broad, the most common and obvious of them are 1) uncontrolled proliferation, 2) resistance to programmed cell death (apoptosis), 3) tissue invasion and metastasis, and 4) sustained angiogenesis. Among the genes affected by cancer, those encoding ion channels are present. Membrane proteins responsible for signaling within cell and among cells, for coupling of extracellular events with intracellular responses, and for maintaining intracellular ionic homeostasis ion channels contribute to various extents to pathophysiological features of each cancer hallmark. Moreover, tight association of these hallmarks with ion channel dysfunction gives a good reason to classify them as special type of channelopathies, namely oncochannelopathies. Although the relation of cancer hallmarks to ion channel dysfunction differs from classical definition of channelopathies, as disease states causally linked with inherited mutations of ion channel genes that alter channel's biophysical properties, in a broader context of the disease state, to which pathogenesis ion channels essentially contribute, such classification seems absolutely appropriate. In this review the authors provide arguments to substantiate such point of view.

Cancer- and behavior-related genes are targeted by selection in the Tasmanian devil (Sarcophilus harrisii)

Devil Facial Tumor Disease (DFTD) is an aggressive cancer notorious for its rare etiology and its impact on Tasmanian devil populations. Two regions underlying an evolutionary response to this cancer were recently identified using genomic time-series pre- and post-DTFD arrival. Here, we support that DFTD shaped the genome of the Tasmanian devil in an even more extensive way than previously reported. We detected 97 signatures of selection, including 148 protein coding genes having a human orthologue, linked to DFTD. Most candidate genes are associated with cancer progression, and an important subset of candidate genes has additional influence on social behavior. This confirms the influence of cancer on the ecology and evolution of the Tasmanian devil. Our work also demonstrates the possibility to detect highly polygenic footprints of short-term selection in very small populations.

Associations between genetic polymorphisms of membrane transporter genes and prognosis after chemotherapy: meta-analysis and finding from Seoul Breast Cancer Study (SEBCS)

Membrane transporters can be major determinants of the pharmacokinetic profiles of anticancer drugs. The associations between genetic variations of ATP-binding cassette (ABC) and solute carrier (SLC) genes and cancer survival were investigated through a meta-analysis and an association study in the Seoul Breast Cancer Study (SEBCS). Including the SEBCS, the meta-analysis was conducted among 38 studies of genetic variations of transporters on various cancer survivors. The population of SEBCS consisted of 1338 breast cancer patients who had been treated with adjuvant chemotherapy. A total of 7750 SNPs were selected from 453 ABC and/or SLC genes typed by an Affymetrix 6.0 chip. ABCB1 rs1045642 was associated with poor progression-free survival in a meta-analysis (HR = 1.33, 95% CI: 1.07-1.64). ABCB1, SLC8A1, and SLC12A8 were associated with breast cancer survival in SEBCS (P_gene < 0.05). ABCB1 rs1202172 was differentially associated with survival depending on the chemotherapy (P_interaction = 0.035). Our finding provides suggestive associations of membrane transporters on cancer survival.

A Decade of GWAS Results in Lung Cancer

Genome-wide association studies (GWAS) were successful to identify genetic factors robustly associated with lung cancer. This review aims to synthesize the literature in this field and accelerate the translation of GWAS discoveries into results that are closer to clinical applications. A chronologic presentation of published GWAS on lung cancer susceptibility, survival, and response to treatment is presented. The most important results are tabulated to provide a concise overview in one read. GWAS have reported 45 lung cancer susceptibility loci with varying strength of evidence and highlighted suspected causal genes at each locus. Some genetic risk loci have been refined to more homogeneous subgroups of lung cancer patients in terms of histologic subtypes, smoking status, gender, and ethnicity. Overall, these discoveries are an important step for future development of new therapeutic targets and biomarkers to personalize and improve the quality of care for patients. GWAS results are on the edge of offering new tools for targeted screening in high-risk individuals, but more research is needed if GWAS are to pay off the investment. Complementary genomic datasets and functional studies are needed to refine the underlying molecular mechanisms of lung cancer preliminarily revealed by GWAS and reach results that are medically actionable. Cancer Epidemiol Biomarkers Prev; 27(4); 363-79. ©2018 AACRSee all articles in this CEBP Focus section, "Genome-Wide Association Studies in Cancer."

Genomic profiling in advanced stage non-small-cell lung cancer patients with platinum-based chemotherapy identifies germline variants with prognostic value in SMYD2

OBJECTIVE

The aim of the study was to investigate the relationship between germline variations as a prognosis biomarker in patients with advanced Non-Small-Cell-Lung-Cancer (NSCLC) subjected to first-line platinum-based treatment.

MATERIALS AND METHODS

We carried out a two-stage genome-wide-association study in non-small-cell lung cancer patients with platinum-based chemotherapy in an exploratory sample of 181 NSCLC patients from Caucasian origin, followed by a validation on 356 NSCLC patients from the same ancestry (Valencia, Spain).

RESULTS

We identified germline variants in SMYD2 as a prognostic factor for survival in patients with advanced NSCLC receiving chemotherapy. SMYD2 alleles are associated to a decreased overall survival and with a reduced Time to Progression. In addition, enrichment pathway analysis identified 361 variants in 40 genes to be involved in poorer outcome in advanced-stage NSCLC patients.

CONCLUSION

Germline SMYD2 alleles are associated with bad clinical outcome of first-line platinum-based treatment in advanced NSCLC patients. This result supports the role of SMYD2 in the carcinogenic process, and might be used as prognostic signature directing patient stratification and the choice of therapy.

MICROABSTRACT

A two-Stage Genome wide association study in Caucasian population reveals germline genetic variation in SMYD2 associated to progression disease in first-line platinum-based treatment in advanced NSCLC patients. SMYD2 profiling might have prognostic / predictive value directing choice of therapy and enlighten current knowledge on pathways involved in human carcinogenesis as well in resistance to chemotherapy.

Regulatory variants in cancer-related pathway genes predict survival of patients with surgically resected non-small cell lung cancer

BACKGROUND

We conducted this study to identify genetic variants in cancer-related pathway genes which can predict prognosis of NSCLC patients after surgery, using a comprehensive list of regulatory single nucleotide polymorphisms (SNPs) prioritized by RegulomeDB.

METHOD

A total of 509 potentially functional SNPs in cancer-related pathway genes selected from RegulomeDB were evaluated. These SNPs were analyzed in a discovery set (n=354), and a replication study was performed in an independent set (n=772). The association of the SNPs with overall survival (OS) and disease-free survival (DFS) were analyzed.

RESULTS

In the discovery set, 76 SNPs were significantly associated with OS or DFS. Among the 76 SNPs, the association was consistently observed for 5 SNPs (ERCC1 rs2298881C>A, BRCA2 rs3092989G>A, NELFE rs440454C>T, PPP2R4 rs2541164G>A, and LTBP4 rs3786527G>A) in the validation set. In combined analysis, ERCC1 rs2298881C>A, BRCA2 rs3092989, NELFE rs440454C>T, and PPP2R4 rs2541164G>A were significantly associated with OS and DFS (adjusted HR ·aHR· for OS=1.46, 0.62, 078, and 0.76, respectively; P=0.003, 0.002, 0.007, and 0.003 respectively; and aHR for DFS=1.27, 0.69, 0.86, and 0.82, respectively; P=0.02, 0.002, 0.03, and 0.008, respectively). The LTBP4 rs3786527G>A was significantly associated with better OS (aHR=0.75; P=0.003).

CONCLUSION

Our results suggest that five SNPs in the cancer-related pathway genes may be useful for the prediction of the prognosis in patients with surgically resected NSCLC.

KCa3.1 channel inhibition leads to an ICAM-1 dependent increase of cell-cell adhesion between A549 lung cancer and HMEC-1 endothelial cells

Early metastasis leads to poor prognosis of lung cancer patients, whose 5-year survival rate is only 15%. We could recently show that the Ca²⁺ sensitive K⁺ channel K_Ca3.1 promotes aggressive behavior of non-small cell lung cancer (NSCLC) cells and that it can serve as a prognostic marker in NSCLC. Since NSCLC patients die of metastases, we investigated whether K_Ca3.1 channels contribute to poor patient prognosis by regulating distinct steps of the metastatic cascade. We investigated the extravasation of NSCLC cells and focused on their adhesion to endothelial cells and on transendothelial migration. We quantified the adhesion forces between NSCLC cells and endothelial cells by applying single cell force spectroscopy, and we monitored transendothelial migration using live-cell imaging. Inhibition of K_Ca3.1 channels with senicapoc or K_Ca3.1 silencing increases the adhesion force of A549 lung cancer cells to human microvascular endothelial cells (HMEC-1). Western blotting, immunofluorescence staining and biotinylation assays indicate that the elevated adhesion force is due to increased expression of ICAM-1 in both cell lines when K_Ca3.1 channels are downregulated. Consistent with this interpretation, an anti-ICAM-1 blocking antibody abolishes the K_Ca3.1-dependent increase in adhesion. Senicapoc inhibits transendothelial migration of A549 cells by 50%. Selectively silencing K_Ca3.1 channels in either NSCLC or endothelial cells reveals that transendothelial migration depends predominantly on endothelial K_Ca3.1 channels.

In conclusion, our findings disclose a novel function of K_Ca3.1 channels in cancer. K_Ca3.1 channels regulate ICAM-1 dependent cell-cell adhesion between endothelial and cancer cells that affects the transmigration step of the metastatic cascade.

Germline Genetic Variants and Lung Cancer Survival in African Americans

Background: African Americans have the highest lung cancer mortality in the United States. Genome-wide association studies (GWASs) of germline variants influencing lung cancer survival have not yet been conducted with African Americans. We examined five previously reported GWAS catalog variants and explored additional genome-wide associations among African American lung cancer cases.Methods: Incident non-small cell lung cancer cases (N = 286) in the Southern Community Cohort Study were genotyped on the Illumina HumanExome BeadChip. We used Cox proportional hazards models to estimate HRs and 95% confidence intervals (CIs) for overall mortality. Two independent African American studies (N = 316 and 298) were used for replication.Results: One previously reported variant, rs1878022 on 12q23.3, was significantly associated with mortality (HR = 0.70; 95% CI: 0.54-0.92). Replication findings were in the same direction, although attenuated (HR = 0.87 and 0.94). Meta-analysis had a HR of 0.83 (95% CI, 0.71-0.97). Analysis of common variants identified an association between chromosome 6q21.33 and mortality (HR = 0.46; 95% CI, 0.33-0.66).Conclusions: We identified an association between rs1878022 in CMKLR1 and lung cancer survival. However, our results in African Americans have a different direction of effect compared with a prior study in European Americans, suggesting a different genetic architecture or presence of gene-environment interactions. We also identified variants on chromosome 6 within the gene-rich HLA region, which has been previously implicated in lung cancer risk and survival.Impact: We found evidence that inherited genetic risk factors influence lung cancer survival in African Americans. Replication in additional populations is necessary to confirm potential genetic differences in lung cancer survival across populations. Cancer Epidemiol Biomarkers Prev; 26(8); 1288-95. ©2017 AACR.

Confirmation of five novel susceptibility loci for systemic lupus erythematosus (SLE) and integrated network analysis of 82 SLE susceptibility loci

We recently identified ten novel SLE susceptibility loci in Asians and uncovered several additional suggestive loci requiring further validation. This study aimed to replicate five of these suggestive loci in a Han Chinese cohort from Hong Kong, followed by meta-analysis (11,656 cases and 23,968 controls) on previously reported Asian and European populations, and to perform bioinformatic analyses on all 82 reported SLE loci to identify shared regulatory signatures. We performed a battery of analyses for these five loci, as well as joint analyses on all 82 SLE loci. All five loci passed genome-wide significance: MYNN (rs10936599, Pmeta = 1.92 × 10-13, OR = 1.14), ATG16L2 (rs11235604, Pmeta = 8.87 × 10 -12, OR = 0.78), CCL22 (rs223881, Pmeta = 5.87 × 10-16, OR = 0.87), ANKS1A (rs2762340, Pmeta = 4.93 × 10-15, OR = 0.87) and RNASEH2C (rs1308020, Pmeta = 2.96 × 10-19, OR = 0.84) and co-located with annotated gene regulatory elements. The novel loci share genetic signatures with other reported SLE loci, including effects on gene expression, transcription factor binding, and epigenetic characteristics. Most (56%) of the correlated (r2 > 0.8) SNPs from the 82 SLE loci were implicated in differential expression (9.81 × 10-198 < P < 5 × 10-3) of cis-genes. Transcription factor binding sites for p53, MEF2A and E2F1 were significantly (P < 0.05) over-represented in SLE loci, consistent with apoptosis playing a critical role in SLE. Enrichment analysis revealed common pathways, gene ontology, protein domains, and cell type-specific expression. In summary, we provide evidence of five novel SLE susceptibility loci. Integrated bioinformatics using all 82 loci revealed that SLE susceptibility loci share many gene regulatory features, suggestive of conserved mechanisms of SLE etiopathogenesis.

Genetic susceptibility variants for lung cancer: replication study and assessment as expression quantitative trait loci

Many single nucleotide polymorphisms (SNPs) have been associated with lung cancer but lack confirmation and functional characterization. We retested the association of 56 candidate SNPs with lung adenocarcinoma risk and overall survival in a cohort of 823 Italian patients and 779 healthy controls, and assessed their function as expression quantitative trait loci (eQTLs). In the replication study, eight SNPs (rs401681, rs3019885, rs732765, rs2568494, rs16969968, rs6495309, rs11634351, and rs4105144) associated with lung adenocarcinoma risk and three (rs9557635, rs4105144, and rs735482) associated with survival. Five of these SNPs acted as cis-eQTLs, being associated with the transcription of IREB2 (rs2568494, rs16969968, rs11634351, rs6495309), PSMA4 (rs6495309) and ERCC1 (rs735482), out of 10,821 genes analyzed in lung. For these three genes, we obtained experimental evidence of differential allelic expression in lung tissue, pointing to the existence of in-cis genomic variants that regulate their transcription. These results suggest that these SNPs exert their effects on cancer risk/outcome through the modulation of mRNA levels of their target genes.

Anks1a regulates COPII-mediated anterograde transport of receptor tyrosine kinases critical for tumorigenesis

ErbB2 signalling, which is amplified by EphA2 binding, is an important therapeutic target for breast cancer. Despite the importance of the EphA2/ErbB2 complex in promoting breast tumorigenesis, the mechanism by which these receptor tyrosine kinases (RTKs) are exported from the endoplasmic reticulum (ER) remains poorly understood. Here we report that the PTB adaptor Anks1a is specifically localized to the ER on its own serine phosphorylation. Once there, Anks1a acts as an important regulator of COPII-mediated EphA2 ER export. The Anks1a ankyrin repeat domain binds EphA2 and causes it to accumulate at sites of ER exit. Simultaneously, the Anks1a PTB domain binds Sec23. This induces internalization of EphA2 via COPII vesicles, while Anks1a remains behind on the ER membrane. EphA2 also binds ErbB2 in the ER and seems to load ErbB2 into growing COPII carriers. Together, our study reveals a novel mechanism that regulates the loading of RTKs into COPII vesicles.

EphA2/ErbB2 complex is important in promoting breast cancer but the mechanism by which these receptor tyrosine kinases are exported from the endoplasmic reticulum is unknown. Here the authors show that Anks1a acts as a cargo adaptor in sorting EphA2 into COPII vesicles, thus modulating the surface level of EphA2.

The study of the relation of DNA repair pathway genes SNPs and the sensitivity to radiotherapy and chemotherapy of NSCLC

To analyze the relation between SNPs in DNA repair pathway-related genes and sensitivity of tumor radio-chemotherapy, 26 SNPs in 20 DNA repair genes were genotyped on 176 patients of NSCLC undertaking radio-chemotherapy treatment. In squamous cell carcinoma (SCC), as the rs2228000, rs2228001 (XPC), rs2273953 (TP73), rs2279744 (MDM2), rs2299939 (PTEN) and rs8178085, rs12334811 (DNA-PKcs) affected the sensitivity to chemotherapy, so did the rs8178085, rs12334811 to radiotherapy. Moreover rs344781, rs2273953 and rs12334811 were related with the survival time of SCC. In general, the “good” genotype GG (rs12334811) showed greater efficacy of radio-chemotherapy and MSF (24 months) on SCC. In adenocarcinoma, as the rs2699887 (PIK3), rs12334811 (DNA-PKcs) influenced the sensitivity to chemotherapy, so did the rs2299939, rs2735343 (PTEN) to radiotherapy. And rs402710, rs80270, rs2279744 and rs2909430 impacted the survival time of the adenocarcinoma patients. Both GG (rs2279744) and AG (rs2909430) showed a shorter survival time (MFS = 6). Additionally, some SNPs such as rs2228000, rs2228001 and rs344781 were found to regulate the expression of DNA repair pathway genes through eQTLs dataset analysis. These results indicate that SNPs in DNA repair pathway genes might regulate the expression and affect the DNA damage repair, and thereby impact the efficacy of radio-chemotherapy and the survival time of NSCLC.

Cytochrome P450 Allele CYP3A7*1C Associates with Adverse Outcomes in Chronic Lymphocytic Leukemia, Breast, and Lung Cancer

CYP3A enzymes metabolize endogenous hormones and chemotherapeutic agents used to treat cancer, thereby potentially affecting drug effectiveness. Here, we refined the genetic basis underlying the functional effects of a CYP3A haplotype on urinary estrone glucuronide (E1G) levels and tested for an association between CYP3A genotype and outcome in patients with chronic lymphocytic leukemia (CLL), breast, or lung cancers. The most significantly associated SNP was rs45446698, an SNP that tags the CYP3A7*1C allele; this SNP was associated with a 54% decrease in urinary E1G levels. Genotyping this SNP in 1,008 breast cancer, 1,128 lung cancer, and 347 CLL patients, we found that rs45446698 was associated with breast cancer mortality (HR, 1.74; P = 0.03), all-cause mortality in lung cancer patients (HR, 1.43; P = 0.009), and CLL progression (HR, 1.62; P = 0.03). We also found borderline evidence of a statistical interaction between the CYP3A7*1C allele, treatment of patients with a cytotoxic agent that is a CYP3A substrate, and clinical outcome (Pinteraction = 0.06). The CYP3A7*1C allele, which results in adult expression of the fetal CYP3A7 gene, is likely to be the functional allele influencing levels of circulating endogenous sex hormones and outcome in these various malignancies. Further studies confirming these associations and determining the mechanism by which CYP3A7*1C influences outcome are required. One possibility is that standard chemotherapy regimens that include CYP3A substrates may not be optimal for the approximately 8% of cancer patients who are CYP3A7*1C carriers.

The Relevance of Genomic Signatures at Adhesion GPCR Loci in Humans

Adhesion G protein-coupled receptors (aGPCRs) have a long evolutionary history dating back to very basal unicellular eukaryotes. Almost every vertebrate is equipped with a set of different aGPCRs. Genomic sequence data of several hundred extinct and extant species allows for reconstruction of aGPCR phylogeny in vertebrates and non-vertebrates in general but also provides a detailed view into the recent evolutionary history of human aGPCRs. Mining these sequence sources with bioinformatic tools can unveil many facets of formerly unappreciated aGPCR functions. In this review, we extracted such information from the literature and open public sources and provide insights into the history of aGPCR in humans. This includes comprehensive analyses of signatures of selection, variability of human aGPCR genes, and quantitative traits at human aGPCR loci. As indicated by a large number of genome-wide genotype-phenotype association studies, variations in aGPCR contribute to specific human phenotypes. Our survey demonstrates that aGPCRs are significantly involved in adaptation processes, phenotype variations, and diseases in humans.

PTEN polymorphisms contribute to clinical outcomes of advanced lung adenocarcinoma patients treated with platinum-based chemotherapy

This study aimed to elucidate the impact of PTEN single nucleotide polymorphism (SNP) on clinical outcomes for advanced lung adenocarcinoma (LAC) patients treated with platinum-based chemotherapy. Three functional SNPs (rs11202607 G>A, rs701848 A>G, and rs11202592 G>C) of PTEN gene were genotyped by using DNA from blood samples of 618 advanced LAC patients, and their relationships with clinical outcomes were analyzed. The carriers of homozygous mutant of rs701848 and rs11202592 polymorphisms revealed significantly worse overall survival (OS) than those with heterozygote or wild-type homozygote (18.83 vs. 21.47 vs. 24.37 months, P = 0.034 and 13.40 vs. 19.03 vs. 21.90 months, P = 0.025, respectively). Subgroup analysis revealed that this association was particularly significant in tumor-lymph-node metastasis (TNM) stage III patients. The objective response rates (ORR) and disease control rates (DCR) of patients with genotype AA, AG, and GG in PTEN rs701848 polymorphism were statistically different (24.1 vs 16.6 vs 12.2 %, P = 0.017 and 82.7 vs 76.0 vs 70.2 %, P = 0.029, respectively). Haplotype analysis revealed a protective effect of the haplotype G-A-A (in the order of rs11202592, rs701848, and rs11202607) on chemotherapy efficacy and survival. Taken together, PTEN polymorphisms may contribute to survival and chemotherapy efficacy of advanced LAC patients treated with platinum-based agents.

Integration of Genome-Wide SNP Data and Gene-Expression Profiles Reveals Six Novel Loci and Regulatory Mechanisms for Amino Acids and Acylcarnitines in Whole Blood

Profiling amino acids and acylcarnitines in whole blood spots is a powerful tool in the laboratory diagnosis of several inborn errors of metabolism. Emerging data suggests that altered blood levels of amino acids and acylcarnitines are also associated with common metabolic diseases in adults. Thus, the identification of common genetic determinants for blood metabolites might shed light on pathways contributing to human physiology and common diseases. We applied a targeted mass-spectrometry-based method to analyze whole blood concentrations of 96 amino acids, acylcarnitines and pathway associated metabolite ratios in a Central European cohort of 2,107 adults and performed genome-wide association (GWA) to identify genetic modifiers of metabolite concentrations. We discovered and replicated six novel loci associated with blood levels of total acylcarnitine, arginine (both on chromosome 6; rs12210538, rs17657775), propionylcarnitine (chromosome 10; rs12779637), 2-hydroxyisovalerylcarnitine (chromosome 21; rs1571700), stearoylcarnitine (chromosome 1; rs3811444), and aspartic acid traits (chromosome 8; rs750472). Based on an integrative analysis of expression quantitative trait loci in blood mononuclear cells and correlations between gene expressions and metabolite levels, we provide evidence for putative causative genes: SLC22A16 for total acylcarnitines, ARG1 for arginine, HLCS for 2-hydroxyisovalerylcarnitine, JAM3 for stearoylcarnitine via a trans-effect at chromosome 1, and PPP1R16A for aspartic acid traits. Further, we report replication and provide additional functional evidence for ten loci that have previously been published for metabolites measured in plasma, serum or urine. In conclusion, our integrative analysis of SNP, gene-expression and metabolite data points to novel genetic factors that may be involved in the regulation of human metabolism. At several loci, we provide evidence for metabolite regulation via gene-expression and observed overlaps with GWAS loci for common diseases. These results form a strong rationale for subsequent functional and disease-related studies.

A penalized robust semiparametric approach for gene-environment interactions

In genetic and genomic studies, gene-environment (G×E) interactions have important implications. Some of the existing G×E interaction methods are limited by analyzing a small number of G factors at a time, by assuming linear effects of E factors, by assuming no data contamination, and by adopting ineffective selection techniques. In this study, we propose a new approach for identifying important G×E interactions. It jointly models the effects of all E and G factors and their interactions. A partially linear varying coefficient model is adopted to accommodate possible nonlinear effects of E factors. A rank-based loss function is used to accommodate possible data contamination. Penalization, which has been extensively used with high-dimensional data, is adopted for selection. The proposed penalized estimation approach can automatically determine if a G factor has an interaction with an E factor, main effect but not interaction, or no effect at all. The proposed approach can be effectively realized using a coordinate descent algorithm. Simulation shows that it has satisfactory performance and outperforms several competing alternatives. The proposed approach is used to analyze a lung cancer study with gene expression measurements and clinical variables. Copyright © 2015 John Wiley & Sons, Ltd.

Using Drosophila melanogaster to identify chemotherapy toxicity genes

The severity of the toxic side effects of chemotherapy shows a great deal of interindividual variability, and much of this variation is likely genetically based. Simple DNA tests predictive of toxic side effects could revolutionize the way chemotherapy is carried out. Due to the challenges in identifying polymorphisms that affect toxicity in humans, we use Drosophila fecundity following oral exposure to carboplatin, gemcitabine and mitomycin C as a model system to identify naturally occurring DNA variants predictive of toxicity. We use the Drosophila Synthetic Population Resource (DSPR), a panel of recombinant inbred lines derived from a multiparent advanced intercross, to map quantitative trait loci affecting chemotoxicity. We identify two QTL each for carboplatin and gemcitabine toxicity and none for mitomycin. One QTL is associated with fly orthologs of a priori human carboplatin candidate genes ABCC2 and MSH2, and a second QTL is associated with fly orthologs of human gemcitabine candidate genes RRM2 and RRM2B. The third, a carboplatin QTL, is associated with a posteriori human orthologs from solute carrier family 7A, INPP4A&B, and NALCN. The fourth, a gemcitabine QTL that also affects methotrexate toxicity, is associated with human ortholog GPx4. Mapped QTL each explain a significant fraction of variation in toxicity, yet individual SNPs and transposable elements in the candidate gene regions fail to singly explain QTL peaks. Furthermore, estimates of founder haplotype effects are consistent with genes harboring several segregating functional alleles. We find little evidence for nonsynonymous SNPs explaining mapped QTL; thus it seems likely that standing variation in toxicity is due to regulatory alleles.

Germline polymorphisms and survival of lung adenocarcinoma patients: a genome-wide study in two European patient series

In lung cancer, the survival of patients with the same clinical stage varies widely for unknown reasons. In this two-phase study, we examined the hypothesis that germline variations influence the survival of patients with lung adenocarcinoma. First, we analyzed existing genotype and clinical data from 289 UK-resident patients with lung adenocarcinoma, identifying 86 single nucleotide polymorphisms (SNPs) that associated with survival (p < 0.01). We then genotyped these candidate SNPs in a validation series of 748 patients from Italy that resulted genetically compatible with the UK series based on principal component analysis. In a Cox proportional hazard model adjusted for age, sex and clinical stage, four SNPs were confirmed on the basis of their having a hazard ratio (HR) indicating the same direction of effect in the two series and p < 0.05. The strongest association was provided by rs2107561, an intronic SNP of PTPRG, protein tyrosine phosphatase, receptor type, G; the C allele was associated with poorer survival in both patient series (pooled analysis loge HR = 0.31; 95% CI: 0.15-0.46, p = 8.5 × 10(-5) ). PTPRG mRNA levels in 43 samples of lung adenocarcinoma were 40% of those observed in noninvolved lung tissue from the same patients. PTPRG overexpression significantly inhibited the clonogenicity of A549 lung carcinoma cells and the anchorage-independent growth of the NCI-H460 large cell lung cancer line. These four germline variants represent promising candidates that, with further study, may help predict clinical outcome. In addition, the PTPRG locus may have a role in tumor progression.

Insight in genome-wide association of metabolite quantitative traits by exome sequence analyses

Metabolite quantitative traits carry great promise for epidemiological studies, and their genetic background has been addressed using Genome-Wide Association Studies (GWAS). Thus far, the role of less common variants has not been exhaustively studied. Here, we set out a GWAS for metabolite quantitative traits in serum, followed by exome sequence analysis to zoom in on putative causal variants in the associated genes. ¹H Nuclear Magnetic Resonance (¹H-NMR) spectroscopy experiments yielded successful quantification of 42 unique metabolites in 2,482 individuals from The Erasmus Rucphen Family (ERF) study. Heritability of metabolites were estimated by SOLAR. GWAS was performed by linear mixed models, using HapMap imputations. Based on physical vicinity and pathway analyses, candidate genes were screened for coding region variation using exome sequence data. Heritability estimates for metabolites ranged between 10% and 52%. GWAS replicated three known loci in the metabolome wide significance: CPS1 with glycine (P-value  = 1.27×10⁻³²), PRODH with proline (P-value  = 1.11×10⁻¹⁹), SLC16A9 with carnitine level (P-value  = 4.81×10⁻¹⁴) and uncovered a novel association between DMGDH and dimethyl-glycine (P-value  = 1.65×10⁻¹⁹) level. In addition, we found three novel, suggestively significant loci: TNP1 with pyruvate (P-value  = 1.26×10⁻⁸), KCNJ16 with 3-hydroxybutyrate (P-value  = 1.65×10⁻⁸) and 2p12 locus with valine (P-value  = 3.49×10⁻⁸). Exome sequence analysis identified potentially causal coding and regulatory variants located in the genes CPS1, KCNJ2 and PRODH, and revealed allelic heterogeneity for CPS1 and PRODH. Combined GWAS and exome analyses of metabolites detected by high-resolution ¹H-NMR is a robust approach to uncover metabolite quantitative trait loci (mQTL), and the likely causative variants in these loci. It is anticipated that insight in the genetics of intermediate phenotypes will provide additional insight into the genetics of complex traits.

Human metabolic individuality is under strict control of genetic and environmental factors. In our study, we aimed to find the genetic determinants of circulating molecules in sera of large set of individuals representing the general population. First, we performed a hypothesis-free genome wide screen in this population to identify genetic regions of interest. Our study confirmed four known gene metabolite connections, but also pointed to four novel ones. Genome-wide screens enriched for common intergenic variants may miss causal genetic variations directly changing the protein sequence. To investigate this further, we zoomed into regions of interest and tested whether the association signals obtained in the first stage were direct, or whether they represent causal variations, which were not captured in the initial panel. These subsequent tests showed that protein coding and regulatory variations are involved in metabolite levels. For two genomic regions we also found that genes harbour more than one causal variant influencing metabolite levels independent of each other. We also observed strong connection between markers of cardio-metabolic health and metabolites. Taken together, our novel loci are of interest for further research to investigate the causal relation to for instance type 2 diabetes and cardiovascular disease.

Acetylcholine receptor pathway in lung cancer: New twists to an old story

Genome-wide association studies revealed that allelic variation in the α5-α3-β4 nicotine acetylcholine receptor (nAChR) cluster on chromosome 15q24-15q25.1 was associated with lung cancer risk. nAChRs are membrane ligand-gated cation channels whose activation is triggered by the binding of the endogenous neurotransmitter acetylcholine (ACh) or other biologic compounds including nicotine. nAChRs have been found on lung cancer cells, underscoring the idea that the non-neuronal nAChR pathway has important implications for lung cancer. Several studies focusing on the treatment with nAChR antagonists with improved selectivity might trigger novel strategies for the intervention and prevention of lung cancer. Here we review the genetic risk factors for lung cancer in the nAChR gene cluster, the roles of nicotine receptors, and the molecular mechanisms of acetylcholine receptor pathways to lead to more opportunities for intervention and prevention of lung cancer.

Genomic rearrangements define lineage relationships between adjacent lepidic and invasive components in lung adenocarcinoma

The development of adenocarcinoma of the lung is believed to proceed from in situ disease (adenocarcinoma in situ, AIS) to minimally invasive disease with prominent lepidic growth (minimally invasive adenocarcinoma, MIA), then to fully invasive adenocarcinoma (AD), but direct evidence for this model has been lacking. Because some lung adenocarcinomas show prominent lepidic growth (AD-L), we designed a study to address the lineage relationship between the lepidic (noninvasive) component (L) and the adjacent nonlepidic growth component representing invasive disease within individual tumors. Lineage relationships were evaluated by next-generation DNA sequencing to define large genomic rearrangements in microdissected tissue specimens collected by laser capture. We found a strong lineage relationship between the majority of adjacent lepidic and invasive components, supporting a putative AIS-AD transition. Notably, many rearrangements were detected in the less aggressive lepidic component, although the invasive component exhibited an overall higher rate of genomic rearrangement. Furthermore, a significant number of genomic rearrangements were present in histologically normal lung adjacent to tumor, but not in host germline DNA, suggesting field defects restricted to zonal regions near a tumor. Our results offer a perspective on the genetic pathogenesis underlying adenocarcinoma development and its clinical management.

Genome-wide association study of survival in early-stage non-small cell lung cancer

BACKGROUND

Lung cancer, especially non-small cell lung cancer (NSCLC), is the leading cause of cancer-related deaths all over the world. Studies have indicated that molecular biomarkers, including genetic variants, may provide additional values for the targeted treatments and clinical outcomes of NSCLC patients. To better understand the effects of molecular biomarkers on the treatment of NSCLC, we conducted a genome-wide analysis to investigate the prognostic implications of genetic variants in early-stage NSCLC patients with surgery.

METHODS

A genome wide scan of 906,703 single-nucleotide polymorphisms (SNPs) was conducted in a cohort with 365 early-stage NSCLC patients with surgery, followed by a fast-track replication in another independent cohort of 327 NSCLC patients from Nanjing, China. Cox models were used to screen and validate significant SNPs associated with the overall survival of early-stage NSCLC patients.

RESULTS

We found that rs10023113 in calcium/calmodulin-dependent protein kinase II delta (CAMK2D) was consistently associated with survival of early-stage NSCLC in the GWAS scan and the replication cohort [GWAS scan: hazard ratio (HR) 2.84; 95 % confidence interval (CI) 1.90-4.23, P = 1.29 × 10(-6); replication cohort: HR 2.19, 95 % CI 1.15-4.21, P = 1.80 × 10(-2)]. When combining all the patients, the results showed that the variant allele of rs10023113 was significantly associated with poor prognosis of early-stage NSCLC with P value of 3.40 × 10(-7) (HR 2.30, 95 % CI 1.67-3.17).

CONCLUSIONS

These findings suggest that CAMK2D rs10023113 may be a potentially prognostic marker for overall survival of early-stage NSCLC patients in Chinese population.

Genetic susceptibility to lung cancer and co-morbidities

Lung cancer is a leading cause of cancer death and disease burden in many countries. Understanding of the biological pathways involved in lung cancer aetiology is required to identify key biomolecules that could be of significant clinical value, either as predictive, prognostic or diagnostic markers, or as targets for the development of novel therapies to treat this disease, in addition to smoking avoidance strategies. Genome-wide association studies (GWAS) have enabled significant progress in the past 5 years in investigating genetic susceptibility to lung cancer. Large scale, multi-cohort GWAS of mainly Caucasian, smoking, populations have identified strong associations for lung cancer mapped to chromosomal regions 15q [nicotinic acetylcholine receptor (nAChR) subunits: CHRNA3, CHRNA5], 5p (TERT-CLPTM1L locus) and 6p (BAT3-MSH5). Some studies in Asian populations of smokers have found similar risk loci, whereas GWAS in never smoking Asian females have identified associations in other chromosomal regions, e.g., 3q (TP63), that are distinct from smoking-related lung cancer risk loci. GWAS of smoking behaviour have identified risk loci for smoking quantity at 15q (similar genes to lung cancer susceptibility: CHRNA3, CHRNA5) and 19q (CYP2A6). Other genes have been mapped for smoking initiation and smoking cessation. In chronic obstructive pulmonary disease (COPD), which is a known risk factor for lung cancer, GWAS in large cohorts have also found CHRNA3 and CHRNA5 single nucleotide polymorphisms (SNPs) mapping at 15q as risk loci, as well as other regions at 4q31 (HHIP), 4q24 (FAM13A) and 5q (HTR4). The overlap in risk loci between lung cancer, smoking behaviour and COPD may be due to the effects of nicotine addiction; however, more work needs to be undertaken to explore the potential direct effects of nicotine and its metabolites in gene-environment interaction in these phenotypes. Goals of future genetic susceptibility studies of lung cancer should focus on refining the strongest risk loci in a wide range of populations with lung cancer, and integrating other clinical and biomarker information, in order to achieve the aim of personalised therapy for lung cancer.

The sodium leak channel, NALCN, in health and disease

Ion channels are crucial components of cellular excitability and are involved in many neurological diseases. This review focuses on the sodium leak, G protein-coupled receptors (GPCRs)-activated NALCN channel that is predominantly expressed in neurons where it regulates the resting membrane potential and neuronal excitability. NALCN is part of a complex that includes not only GPCRs, but also UNC-79, UNC-80, NLF-1 and src family of Tyrosine kinases (SFKs). There is growing evidence that the NALCN channelosome critically regulates its ion conduction. Both in mammals and invertebrates, animal models revealed an involvement in many processes such as locomotor behaviors, sensitivity to volatile anesthetics, and respiratory rhythms. There is also evidence that alteration in this NALCN channelosome can cause a wide variety of diseases. Indeed, mutations in the NALCN gene were identified in Infantile Neuroaxonal Dystrophy (INAD) patients, as well as in patients with an Autosomal Recessive Syndrome with severe hypotonia, speech impairment, and cognitive delay. Deletions in NALCN gene were also reported in diseases such as 13q syndrome. In addition, genes encoding NALCN, NLF- 1, UNC-79, and UNC-80 proteins may be susceptibility loci for several diseases including bipolar disorder, schizophrenia, Alzheimer's disease, autism, epilepsy, alcoholism, cardiac diseases and cancer. Although the physiological role of the NALCN channelosome is poorly understood, its involvement in human diseases should foster interest for drug development in the near future. Toward this goal, we review here the current knowledge on the NALCN channelosome in physiology and diseases.

Replication of results of a genome-wide association study on lung cancer survival in a Korean population

Recently, a genome-wide association study (GWAS) identified single nucleotide polymorphisms (SNPs) that may influence the prognosis of early-stage non-small cell lung cancer (NSCLC) in Caucasians. We attempted to replicate the impact of genetic variants identified in the GWAS on lung cancer survival in a Korean population. A total of 363 patients with surgically resected NSCLCs were enrolled, and 12 SNPs were genotyped using the SEQUENOM MassARRAY iPLEX assay, TaqMan assay, or a polymerase chain reaction-restriction fragment length polymorphism analysis. The association between genotypes and overall survival (OS) was analyzed. Among the 12 SNPs, the rs6034368T>C was associated with OS. Patients with the rs6034368C allele showed a better OS than the patients with the rs6034368T allele (adjusted hazard ratio = 0.72, confidence interval = 0.56-0.93, P = 0.01). The rs12446308A>G had an effect on OS, but it was marginally significant (under a codominant model, adjusted hazard ratio = 1.85, confidence interval = 0.98-3.47, P = 0.06). We identified that the rs6034368T>C was associated with survival in early-stage NSCLC in a Korean population.

B-cell malignancies in microRNA Eμ-miR-17~92 transgenic mice

miR-17∼92 is a polycistronic microRNA (miR) cluster (consisting of miR-17, miR-18a, miR-19a, miR-19b, miR-20a, and miR-92a) which frequently is overexpressed in several solid and lymphoid malignancies. Loss- and gain-of-function studies have revealed the role of miR-17∼92 in heart, lung, and B-cell development and in Myc-induced B-cell lymphomas, respectively. Recent studies indicate that overexpression of this locus leads to lymphoproliferation, but no experimental proof that dysregulation of this cluster causes B-cell lymphomas or leukemias is available. To determine whether miR-17∼92- overexpression induces lymphomagenesis/leukemogenesis, we generated a B-cell-specific transgenic mouse model with targeted overexpression of this cluster in B cells. The miR-17∼92 overexpression was driven by the Eµ-enhancer and Ig heavy-chain promoter, and a 3' GFP tag was added to the transgene to track the miR expression. Expression analysis using Northern Blot and quantitative RT-PCR confirmed 2.5- to 25-fold overexpression of all six miRs in the transgenic mice spleens as compared with spleens from wild-type mice. Eµ-miR-17∼92 mice developed B-cell malignancy by the age of 12-18 mo with a penetrance of ∼80% (49% splenic B-cell lymphoproliferative disease, 28% lymphoma). At this stage mice exhibited severe splenomegaly with abnormal B-cell-derived white pulp expansion and enlarged lymph nodes. Interestingly, we found three classes of B-cell lymphomas/leukemias at varying grades of differentiation. These included expansion of CD19(+) and CD5(+) double-positive B cells similar to the aggressive form of human B-cell chronic lymphocytic leukemia, B220(+) CD43(+) B1-cell proliferation, and a CD19(+) aggressive diffuse large B-cell lymphoma-like disease, as assessed by flow cytometry and histopathological analysis.

Gene expression signature of non-involved lung tissue associated with survival in lung adenocarcinoma patients

Lung adenocarcinoma patients of similar clinical stage and undergoing the same treatments often have marked interindividual variations in prognosis. These clinical discrepancies may be due to the genetic background modulating an individual's predisposition to fighting cancer. Herein, we hypothesized that the lung microenvironment, as reflected by its expression profile, may affect lung adenocarcinoma patients' survival. The transcriptome of non-involved lung tissue, excised from a discovery series of 204 lung adenocarcinoma patients, was evaluated using whole-genome expression microarrays (with probes corresponding to 28 688 well-annotated coding sequences). Genes associated with survival status at 60 months were identified by Cox regression analysis (adjusted for gender, age and clinical stage) and retested in a validation series of 78 additional cases. RNA-Seq analysis from non-involved lung tissue of 12 patients was performed to characterize the different isoforms of candidate genes. Ten genes for which the loge-transformed hazard ratios expressed the same direction of effect in the discovery (P < 1.0 × 10(-3)) and validation series comprised the gene expression signature associated with survival: CNTNAP1, PKNOX1, FAM156A, FRMD8, GALNTL1, TXNDC12, SNTB1, PPP3R1, SNX10 and SERPINH1. RNA sequencing highlighted the complex expression pattern of these genes in non-involved lung tissue from different patients and permitted the detection of a read-through gene fusion between PPP3R1 and the flanking gene (CNRIP1) as well as a novel isoform of CNTNAP1. Our findings support the hypothesis that individual genetic characteristics, evidenced by the expression pattern of non-involved tissue, influence the outcome of lung adenocarcinoma patients.

Recessive truncating NALCN mutation in infantile neuroaxonal dystrophy with facial dysmorphism

BACKGROUND

Infantile neuroaxonal dystrophy (INAD) is a recessive disease that results in total neurological degeneration and death in childhood. PLA2G6 mutation is the underlying genetic defect, but rare genetic heterogeneity has been demonstrated. One of the five families we studied did not link to PLA2G6 locus, and in the family one of the two affected siblings additionally had atypical features including facial dysmorphism, pectus carinatum, scoliosis, pes varus, zygodactyly and bilateral cryptorchidism as well as cerebellar atrophy, as previously reported.

METHODS

Sural biopsy was investigated by electron microscopy. PLA2G6 was screened for mutations by Sanger sequencing. In the mutation-free family, candidate disease loci were found via linkage analysis using data from single nucleotide polymorphism genome scans. Exome sequencing was applied to find the variants at the loci.

RESULTS

PLA2G6 mutations were identified in four families including the one with an unusually severe phenotype that led to death within the first 2 years of life. In the remaining family, seven candidate loci totalling 15.2 Mb were found and a homozygous truncating mutation p.Q642X was identified in NALCN at 13q32.3. The patients are around 20-years-old.

CONCLUSIONS

NALCN is the gene responsible for INAD with facial dysmorphism. The patients have lived to adulthood despite severe growth and neuromotor retardation. NALCN forms a voltage-independent ion channel with a role in the regulation of neuronal excitability. Our findings broaden the spectrum of genes associated with neuroaxonal dystrophy. Testing infants with idiopathic severe growth retardation and neurodegeneration for NALCN mutations could benefit families.