Machine learning-based immune phenotypes correlate with STK11/KEAP1 co-mutations and prognosis in resectable NSCLC: a sub-study of the TNM-I trial.

BACKGROUND: We aim to implement an immune cell score model in routine clinical practice for resected non-small-cell lung cancer (NSCLC) patients (NCT03299478). Molecular and genomic features associated with immune phenotypes in NSCLC have not been explored in detail. PATIENTS AND METHODS: We developed a machine learning (ML)-based model to classify tumors into one of three categories: inflamed, altered, and desert, based on the spatial distribution of CD8+ T cells in two prospective (n = 453; TNM-I trial) and retrospective (n = 481) stage I-IIIA NSCLC surgical cohorts. NanoString assays and targeted gene panel sequencing were used to evaluate the association of gene expression and mutations with immune phenotypes. RESULTS: Among the total of 934 patients, 24.4% of tumors were classified as inflamed, 51.3% as altered, and 24.3% as desert. There were significant associations between ML-derived immune phenotypes and adaptive immunity gene expression signatures. We identified a strong association of the nuclear factor-κB pathway and CD8+ T-cell exclusion through a positive enrichment in the desert phenotype. KEAP1 [odds ratio (OR) 0.27, Q = 0.02] and STK11 (OR 0.39, Q = 0.04) were significantly co-mutated in non-inflamed lung adenocarcinoma (LUAD) compared to the inflamed phenotype. In the retrospective cohort, the inflamed phenotype was an independent prognostic factor for prolonged disease-specific survival and time to recurrence (hazard ratio 0.61, P = 0.01 and 0.65, P = 0.02, respectively). CONCLUSIONS: ML-based immune phenotyping by spatial distribution of T cells in resected NSCLC is able to identify patients at greater risk of disease recurrence after surgical resection. LUADs with concurrent KEAP1 and STK11 mutations are enriched for altered and desert immune phenotypes.

A tissue-bioengineering strategy for modeling rare human kidney diseases in vivo.

The lack of animal models for some human diseases precludes our understanding of disease mechanisms and our ability to test prospective therapies in vivo. Generation of kidney organoids from Tuberous Sclerosis Complex (TSC) patient-derived-hiPSCs allows us to recapitulate a rare kidney tumor called angiomyolipoma (AML). Organoids derived from TSC2-/- hiPSCs but not from isogenic TSC2+/- or TSC2+/+ hiPSCs share a common transcriptional signature and a myomelanocytic cell phenotype with kidney AMLs, and develop epithelial cysts, replicating two major TSC-associated kidney lesions driven by genetic mechanisms that cannot be consistently recapitulated with transgenic mice. Transplantation of multiple TSC2-/- renal organoids into the kidneys of immunodeficient rats allows us to model AML in vivo for the study of tumor mechanisms, and to test the efficacy of rapamycin-loaded nanoparticles as an approach to rapidly ablate AMLs. Collectively, our experimental approaches represent an innovative and scalable tissue-bioengineering strategy for modeling rare kidney disease in vivo.

Inhibition of MAPK pathway is essential for suppressing Rheb-Y35N driven tumor growth.

Rheb is a Ras family GTPase, which binds to and activates mammalian target of rapamycin complex 1 (mTORC1) when GTP loaded. Recently, cancer genome sequencing efforts have identified recurrent Rheb Tyr35Asn mutations in kidney and endometrial carcinoma. Here we show that Rheb-Y35N causes not only constitutive mTORC1 activation, but sustained activation of the MEK-ERK pathway in a TSC1/TSC2/TBC1D7 protein complex and mTORC1-independent manner, contributing to intrinsic resistance to rapamycin. Rheb-Y35N transforms NIH3T3 cells, resulting in aggressive tumor formation in xenograft nude mice, which could be suppressed by combined treatment with rapamycin and an extracellular signal-regulated kinase (ERK) inhibitor. Furthermore, Rheb-Y35N inhibits AMPKα activation in response to nutrient depletion or 5-aminoimidazole-4-carboxamide ribonucleotide (AICAR), leading to attenuated phosphorylation of BRAF-S729 and retained mitogen-activated protein kinase (MAPK) activation. Finally, we demonstrate that Rheb-WT can bind AMPK to facilitate AMPK activation, whereas Rheb-Y35N competitively binds AMPK, impairing AMPK phosphorylation. In summary, our findings indicate that Rheb-Y35N is a dominantly active tumor driver that activates both mTORC1 and MAPK to promote tumor growth, suggesting a combination of mTORC1 and MAPK inhibitors may be of therapeutic value in patients whose cancers sustain this mutation.

A Toll-like receptor-1 variant and its characteristic cellular phenotype is associated with severe malaria in Papua New Guinean children.

Genetic factors are likely to contribute to low severe malaria case fatality rates in Melanesian populations, but association studies can be underpowered and may not provide plausible mechanistic explanations if significant associations are detected. In preparation for a genome-wide association study, 29 candidate single-nucleotide polymorphisms (SNPs) with minor allele frequencies >5% were examined in a case-control study of 504 Papua New Guinean children with severe malaria. In parallel, an immunological substudy was performed on convalescent peripheral blood mononuclear cells (PBMCs) from cases and controls. Following stimulation with a Toll-like receptor (TLR) 1/2 agonist, effector cytokines and chemokines were assayed. The only significant genetic association observed involved a nonsynonymous SNP (TLR1rs4833095) in the TLR1 gene. A recessive (TT) genotype was associated with reduced odds of severe malaria of 0.52 (95% confidence interval (0.29-0.90), P=0.006). Concentrations of pro-inflammatory cytokines interleukin-1ß and tumour necrosis factor α were significantly higher in severe malaria cases compared with healthy controls, but lower in children with the protective recessive (TT) genotype. A genetic variant in TLR1 may contribute to the low severe malaria case fatality rates in this region through a reduced pro-inflammatory cellular phenotype.

Tsc1-Tp53 loss induces mesothelioma in mice, and evidence for this mechanism in human mesothelioma.

Mesothelioma is diagnosed in ∼2500 patients in the United States every year, most often arising in the pleural space, but also occurring as primary peritoneal mesothelioma. The vast majority of patients with mesothelioma die of their disease within 3 years. We developed a new mouse model of mesothelioma by bladder or intraperitoneal injection of adenovirus Cre into mice with conditional alleles of each of Tp53 and Tsc1. Such mice began to develop malignant ascites about 6 months after injection, which was due to peritoneal mesothelioma, on the basis of tumor morphology and immunohistochemical staining. Mesothelioma cell lines were established, which showed loss of both Tsc1 and Tp53, with mammalian target of rapamycin complex (mTORC)1 activation. Treatment of mice with malignant ascites due to mesothelioma with rapamycin led to a marked reduction in ascites, extended survival and a 95-99% reduction in the mesothelioma tumor volume, in comparison with vehicle-treated mice. To see whether TSC1/TSC2 loss was a common genetic event in human mesothelioma, we examined nine human mesothelioma cell lines, and found that four of nine showed persistent activation of mTORC1, although none had loss of TSC1 or TSC2. A tissue microarray analysis of 198 human mesothelioma specimens showed that 33% of cases had reduced TSC2 expression and 60% showed activation of mTOR, indicating that mTOR activation is common in human mesothelioma, suggesting that it is a potential therapeutic target.

Pivotal role of augmented αB-crystallin in tumor development induced by deficient TSC1/2 complex.

Tuberous sclerosis complex 1 (TSC1) and TSC2 are suppressors of mechanistic target of rapamycin (mTOR). mTOR is the major component of two protein complexes: mTOR complex 1 (mTORC1) and mTORC2. Inactive mutation of either TSC1 or TSC2 unleashes mTOR signaling and consequently causes TSC, a benign tumor syndrome affecting multiple organs. We report here that expression of αB-crystallin was upregulated in Tsc1-/- or Tsc2-/- mouse embryonic fibroblasts, Eker rat uterine leiomyoma-derived Tsc2-deficient ELT3 cells, mutant Tsc2-associated mouse kidney tumors, and human lung lymphangioleiomyomatosis nodules. αB-crystallin was transcriptionally activated by mTOR complex 2 (mTORC2): nuclear factor-kappa B (NFκB) signaling cascade. The augmented αB-crystallin was critical for the migration, invasion and apoptotic resistance of Tsc2-defective cells. Disruption of αB-crystallin suppressed Tsc2-null cell proliferation and tumorigenesis. Therefore, enhanced αB-crystallin has an essential role in TSC1/2 complex deficiency-mediated tumorigenesis, and inhibition of αB-crystallin may complement the current therapy for TSC.

Variation in human genes encoding adhesion and proinflammatory molecules are associated with severe malaria in the Vietnamese.

The genetic basis for susceptibility to malaria has been studied widely in African populations but less is known of the contribution of specific genetic variants in Asian populations. We genotyped 67 single-nucleotide polymorphisms (SNPs) in 1030 severe malaria cases and 2840 controls from Vietnam. After data quality control, genotyping data of 956 cases and 2350 controls were analysed for 65 SNPs (3 gender confirmation, 62 positioned in/near 42 malarial candidate genes). A total of 14 SNPs were monomorphic and 2 (rs8078340 and rs33950507) were not in Hardy-Weinberg equilibrium in controls (P<0.01). In all, 7/46 SNPs in 6 genes (ICAM1, IL1A, IL17RC, IL13, LTA and TNF) were associated with severe malaria, with 3/7 SNPs in the TNF/LTA region. Genotype-phenotype correlations between SNPs and clinical parameters revealed that genotypes of rs708567 (IL17RC) correlate with parasitemia (P=0.028, r(2)=0.0086), with GG homozygotes having the lowest parasite burden. Additionally, rs708567 GG homozygotes had a decreased risk of severe malaria (P=0.007, OR=0.78 (95% CI; 0.65-0.93)) and death (P=0.028, OR=0.58 (95% CI; 0.37-0.93)) than those with AA and AG genotypes. In summary, variants in six genes encoding adhesion and proinflammatory molecules are associated with severe malaria in the Vietnamese. Further replicative studies in independent populations will be necessary to confirm these findings.

Amplification of EGFR T790M causes resistance to an irreversible EGFR inhibitor.

Epidermal growth factor receptor (EGFR) tyrosine kinase inhibitors, gefitinib and erlotinib are effective therapies against mutant non-small cell lung cancers (NSCLCs). Treatment is limited by the development of resistance in part explained by the gain of a secondary EGFR mutation, T790M, at the gatekeeper residue. Irreversible EGFR inhibitors, including PF00299804, are effective in vitro and in vivo against EGFR mutant tumors that contain EGFR T790M and are currently under clinical development. In this study, we generate models of resistance to PF00299804, using cell lines with EGFR T790M and show that the PF00299804-resistant models develop focal amplification of EGFR that preferentially involves the T790M-containing allele. These PF00299804-resistant cell lines remain dependent on EGFR for growth as downregulation of EGFR by shRNA compromises their viability. We show that resistance to PF00299804 arises, at least in part, through selection of a pre-existing EGFR T790M-amplified clone both in vitro and using a xenograft model in vivo. Our findings show that EGFR T790M is a common resistance mechanism to both reversible, and when amplified, the irreversible EGFR kinase inhibitors further emphasizing the need to develop more potent therapies against EGFR T790M. These findings can be used to guide studies of patient tumor specimens from ongoing clinical trials of irreversible EGFR kinase inhibitors.

TSC1 loss synergizes with KRAS activation in lung cancer development in the mouse and confers rapamycin sensitivity.

Germline TSC1 or TSC2 mutations cause tuberous sclerosis complex (TSC), a hamartoma syndrome with lung involvement. To explore the potential interaction between TSC1 and KRAS activation in lung cancer, mice in which Tsc1 loss and Kras(G12D) expression occur in a small fraction of lung epithelial cells were generated. Mice with a combined Tsc1-Kras(G12D) mutation had dramatically reduced tumor latency (median survival: 11.6-15.6 weeks) in comparison with Kras(G12D) alone mutant mice (median survival: 27.5 weeks). Tsc1-Kras(G12D) tumors showed consistent activation of mTOR (mammalian target of rapamycin)C1 and responded to treatment with rapamycin, leading to significantly improved survival, whereas rapamycin had minor effects on cancers in Kras(G12D) alone mice. Loss of heterozygosity for TSC1 or TSC2 was found in 22% of 86 human lung cancer specimens. However, none of the 80 lung cancer lines studied showed evidence of the lack of expression of either TSC1 or TSC2 or a signaling pattern corresponding to complete loss. These data indicate that Tsc1 loss synergizes with the Kras mutation to enhance lung tumorigenesis in the mouse, but that this is a rare event in human lung cancer. Rapamycin may have unique benefit for patients with lung cancer, for whom the TSC1/TSC2 function is limited.

Distinct clinical characteristics of tuberous sclerosis complex patients with no mutation identified.

Tuberous Sclerosis Complex (TSC) is a multi-system disorder that is highly variable in its clinical presentation. Current molecular diagnostic methods permit identification of mutations in either TSC1 or TSC2 in 75-85% of TSC patients. Here we examine the clinical characteristics of those TSC patients who have no mutation identified (NMI). A retrospective review of our patient population that had comprehensive testing for mutations in TSC1/TSC2 identified 23/157 (15%) that were NMI. NMI patients had a lower incidence of brain findings on imaging studies, neurological features, and renal findings than those with TSC2 mutations. In contrast, NMI patients had a lower incidence of seizures than TSC patients with TSC1 mutations, but had a higher incidence of both renal angiomyolipomas and pulmonary lymphangioleiomyomatosis. This distinct constellation of findings suggest that NMI patients may have a unique molecular pathogenesis, different from that seen in TSC patients with the usual mutations in TSC1 and TSC2. We suggest that the mechanisms of disease in these patients include both mosaicism for a TSC2 mutation, and unusual non-coding region mutations in TSC2.

BCG-malaria co-Infection has paradoxical effects on C57BL/6 and A/J mouse strains.

Bacillus Calmette-Guérin (BCG) infection of the spleen is a potent modifier of splenic function. Prior to malaria infection, we infected two mouse strains of differing susceptibility to Plasmodium chabaudi AS (C57BL/6 and A/J) with this mycobacterium. We then evaluated aspects of spleen cell composition, architecture and cytokine expression, and correlated these with the outcome. BCG preinfection resulted in protection of the A/J mice but paradoxically resulted in mortality of the C57BL/6 mice. The latter developed higher parasitaemias that peaked earlier than the A/J mice rendered resistant by BCG. BCG infection induced remarkable changes to splenic histology examined by H&E staining, but there were no consistent differences between mouse strains. C57BL/6 mice had higher absolute numbers of all immune cell phenotypes than did A/J mice, and higher macrophage and dendritic cell proportions. BCG-induced resistance in A/J mice was associated with an increased CD4+ expression of IFN-gamma whilst induced death in C57BL/6 mice was associated with excessive IFN-gamma expression. A moderate TH1 response in the A/J model may have been responsible for the improved survival, and an excessive TH1 response in the C57BL/6 model may have contributed to their death.

Haplotypic diversity in human CEACAM genes: effects on susceptibility to meningococcal disease.

Adhesion between the opacity-associated adhesin (Opa) proteins of Neisseria meningitidis and human carcino-embryonic antigen cell adhesion molecule (CEACAM) proteins is an important stage in the pathogenesis of meningococcal disease, a globally important bacterial infection. Most disease is caused by a small number of meningococcal genotypes known as hyperinvasive lineages. As these are also carried asymptomatically, acquisition of them alone cannot explain why only some hosts develop meningococcal disease. Our aim was to determine whether genetic diversity in CEACAM is associated with susceptibility to meningococcal disease. Frequency distributions of alleles, genotypes and haplotypes were compared in four CEACAM genes in 384 case samples and 190 controls. Linkage disequilibrium among polymorphic sites, haplotype structures and relationships were also analysed. A number of polymorphisms were observed in CEACAM genes but the diversity of CEACAM1, to which most Opa proteins bind, was lower, and a small number of high-frequency haplotypes were detected. Dose-dependent associations of three CEACAM haplotypes with meningococcal disease were observed, with the effect of carrying these haplotypes amplified in homozygous individuals. Two haplotypes were protective while one haplotype in CEACAM6 was associated with a twofold increase in disease susceptibility. These data imply that human CEACAM may be one determinant of human susceptibility to meningococcal disease.

Genetic variation at the TNF locus and the risk of severe sequelae of ocular Chlamydia trachomatis infection in Gambians.

Tumor necrosis factor (TNF) is thought to be a key mediator of the inflammatory and fibrotic response to Chlamydia trachomatis (Ct) infection. A large matched-pair case-control study investigated putative functional single nucleotide polymorphisms (SNPs) across the major histocompatibility complex (MHC) class III region, including TNF and its immediate neighbors nuclear factor of kappa light polypeptide gene enhancer in B cells (IkappaBL), inhibitor like 1 and lymphotoxin alpha (LTA) in relation to the risk of scarring sequelae of ocular Ct infection. Haplotype and linkage disequilibrium analysis demonstrated two haplotypes, differing at position TNF-308, conferring an increased risk of trichiasis. The TNF-308A allele, and its bearing haplotype, correlated with increased TNF production in lymphocyte cultures stimulated with chlamydial elementary body antigen. Thus TNF-308A may determine directly, or be a marker of a high TNF producer phenotype associated with increased risk of sequelae of chlamydial infection. Multivariate analysis provided evidence for the presence of additional risk-associated variants near the TNF locus.

Genetic determinants of C-reactive protein in COPD.

Chronic obstructive pulmonary disease (COPD) is associated with a systemic inflammatory state, marked by elevations in serum inflammatory markers including C-reactive protein (CRP). The present study sought to determine epidemiological predictors of CRP levels, to estimate the genetic influence on CRP levels, and to identify genetic variants that affect CRP in a family-based study of COPD. CRP was measured by a high-sensitivity assay in participants from the Boston Early-Onset COPD Study. Predictors of CRP level were determined using multilevel linear models. Variance component analysis was used to estimate heritability and to perform genome-wide linkage analysis for CRP levels. Two variants in the surfactant protein B (SFTPB) gene were tested for association with CRP levels. Increased age, female sex, higher body mass index, greater smoking pack-yrs and reduced forced expiratory volume in one second were all associated with increased CRP levels. There was a significant genetic influence on CRP (heritability = 0.25). Genome-wide linkage analysis revealed several potentially interesting chromosomal regions, though no significant evidence for linkage was found. A short tandem repeat marker near SFTPB was significantly associated with CRP levels. There is a genetic influence on C-reactive protein levels in chronic obstructive pulmonary disease patients. Preliminary evidence suggests an association of the surfactant protein B gene with systemic inflammation in chronic obstructive pulmonary disease.

An investigation of transmission ratio distortion in the central region of the human MHC.

Transmission ratio distortion (TRD) describes a significant departure from expected Mendelian inheritance ratios that is fundamental to both the biology of reproduction and statistical genetics. The relatively high fetal wastage in humans, with consequent selection of alleles in utero, makes it likely that TRD is prevalent in the human genome. The central region of the human major histocompatibility complex (MHC) is a strong TRD candidate, as it houses a number of immune and regulatory genes that may be important in pregnancy outcome. We used a nonhaplotype-based method to select 13 tagging SNPs from three central MHC candidate regions, and analysed their transmission in 380 newborns and their parents (1138 individuals). A TRD of 54:46 was noted in favour of the common allele of a promoter SNP in the CLIC1 gene (P = 0.025), with a similar distortion using haplotypes across the same gene region (P = 0.016). We also found evidence that markers in the CLIC1 gene region may have been subject to recent selection (P < 0.001). The study illustrates the potential benefits of screening for TRD and highlights the difficulties encountered therein.

Increased in vivo transcription of an IL-8 haplotype associated with respiratory syncytial virus disease-susceptibility.

Interleukin-8 (IL-8) has been implicated in the pathogenesis of RSV-induced bronchiolitis. Previously, we have described an association between bronchiolitis disease severity and a specific IL-8 haplotype comprising six single-nucleotide polymorphisms (SNPs) (-251A/+396G/+781T/+1238delA/+1633T/+2767T, haplotype 2). Here we investigated the functional basis for this association by measuring haplotype-specific transcription in vivo in human primary cells. We found a significant increase in transcript level derived from the IL-8 haplotype 2 relative to the mirror haplotype 1 (-251T/+396T/+781C/+1238insA/+1633C/+2767A) in respiratory epithelial cells but not in lymphocytes. A promoter polymorphism, -251A, present on the high producer haplotype, had no significant affect on the allele-specific level of transcription when analyzed in reporter gene experiments in human respiratory epithelial A549 cells. We proceeded to systematically screen for allele-specific protein-DNA binding in this functional haplotype, which revealed significant differential binding at the +781T/C polymorphism. C/EBP beta was identified as being part of a transcription factor binding complex that preferentially bound in the presence of the +781 T allele. These results suggest that the mechanism for disease susceptibility to RSV-induced bronchiolitis may occur through a haplotype-specific increase in IL-8 transcription, which may be mediated by functional polymorphisms within that haplotype.

Complex haplotypic structure of the central MHC region flanking TNF in a West African population.

TNF polymorphisms have been associated with susceptibility to malaria and other infectious and inflammatory conditions. We investigated a sample of 150 West African chromosomes to determine linkage disequilibrium (LD) between 25 SNP markers located in an 80 kb segment of the MHC Class III region encompassing TNF and eight neighbouring genes. We observed 45 haplotypes, and 22 of them comprise 80% of the sample. The pattern of LD is remarkably patchy, such that many markers show no LD with adjacent markers but high LD with markers that are much further away. We introduce a method of examining the implications of LD data for disease association studies based on sample size considerations: this shows that certain TNF polymorphisms would be likely to yield positive associations if the true disease allele resided in LTA or BAT1. We conclude that detailed marker maps are needed to resolve the causal origin of disease associations observed at the TNF locus.

A family with seizures and minor features of tuberous sclerosis and a novel TSC2 mutation.

The authors studied nine members of a family that demonstrated a limited form of tuberous sclerosis complex (TSC). Cutaneous findings were limited to hypopigmented macules in four patients. Five family members had recurrent seizures, and three of these had migrational defects of the cerebral mantle. Mutational analysis of TSC2 indicated the presence of the novel missense change 3106T-->C, 1036S-->P in all family members with seizures. The findings suggest that this mild variant form of TSC is due to a novel TSC2 mutation.