Tepotinib in patients with non-small cell lung cancer with high-level MET amplification detected by liquid biopsy: VISION Cohort B.

High-level MET amplification (METamp) is a primary driver in ∼1%-2% of non-small cell lung cancers (NSCLCs). Cohort B of the phase 2 VISION trial evaluates tepotinib, an oral MET inhibitor, in patients with advanced NSCLC with high-level METamp who were enrolled by liquid biopsy. While the study was halted before the enrollment of the planned 60 patients, the results of 24 enrolled patients are presented here. The objective response rate (ORR) is 41.7% (95% confidence interval [CI], 22.1-63.4), and the median duration of response is 14.3 months (95% CI, 2.8-not estimable). In exploratory biomarker analyses, focal METamp, RB1 wild-type, MYC diploidy, low circulating tumor DNA (ctDNA) burden at baseline, and early molecular response are associated with better outcomes. Adverse events include edema (composite term; any grade: 58.3%; grade 3: 12.5%) and constipation (any grade: 41.7%; grade 3: 4.2%). Tepotinib provides antitumor activity in high-level METamp NSCLC (ClinicalTrials.gov: NCT02864992).

CONET: copy number event tree model of evolutionary tumor history for single-cell data.

Copy number alterations constitute important phenomena in tumor evolution. Whole genome single-cell sequencing gives insight into copy number profiles of individual cells, but is highly noisy. Here, we propose CONET, a probabilistic model for joint inference of the evolutionary tree on copy number events and copy number calling. CONET employs an efficient, regularized MCMC procedure to search the space of possible model structures and parameters. We introduce a range of model priors and penalties for efficient regularization. CONET reveals copy number evolution in two breast cancer samples, and outperforms other methods in tree reconstruction, breakpoint identification and copy number calling.

The genomic and transcriptional landscape of primary central nervous system lymphoma.

Primary lymphomas of the central nervous system (PCNSL) are mainly diffuse large B-cell lymphomas (DLBCLs) confined to the central nervous system (CNS). Molecular drivers of PCNSL have not been fully elucidated. Here, we profile and compare the whole-genome and transcriptome landscape of 51 CNS lymphomas (CNSL) to 39 follicular lymphoma and 36 DLBCL cases outside the CNS. We find recurrent mutations in JAK-STAT, NFkB, and B-cell receptor signaling pathways, including hallmark mutations in MYD88 L265P (67%) and CD79B (63%), and CDKN2A deletions (83%). PCNSLs exhibit significantly more focal deletions of HLA-D (6p21) locus as a potential mechanism of immune evasion. Mutational signatures correlating with DNA replication and mitosis are significantly enriched in PCNSL. TERT gene expression is significantly higher in PCNSL compared to activated B-cell (ABC)-DLBCL. Transcriptome analysis clearly distinguishes PCNSL and systemic DLBCL into distinct molecular subtypes. Epstein-Barr virus (EBV)+ CNSL cases lack recurrent mutational hotspots apart from IG and HLA-DRB loci. We show that PCNSL can be clearly distinguished from DLBCL, having distinct expression profiles, IG expression and translocation patterns, as well as specific combinations of genetic alterations.

Landscape and Clonal Dominance of Co-occurring Genomic Alterations in Non-Small-Cell Lung Cancer Harboring MET Exon 14 Skipping.

MET exon 14 skipping alterations (METex14) comprise a diverse set of actionable oncogene drivers in non-small-cell lung cancer (NSCLC). Recent studies have established the efficacy of tyrosine kinase inhibitors for this patient population. The landscape of co-occurring genetic alterations in METex14 NSCLC and their potential impact to therapeutic sensitivities has not yet been fully described. MATERIALS AND METHODS: METex14 NSCLC cases were collected from three cohorts: the VISION trial, and data sets from Guardant360 and GenePlus. Clinicopathologic characteristics and METex14 mutation sites were analyzed and compared across data sets. Co-occurring genetic alterations and the clonality relationships to METex14 were evaluated. RESULTS: Of 40,824 NSCLCs, 692 METex14 cases (1.7%) were identified, including 332 in Guardant360, 188 in VISION, and 172 in GenePlus. The demographics and mutation type and/or sites were similar in the Asian versus Western cohorts. MET amplification, which were found to be associated with sensitivity to MET kinase inhibitors, co-occurs in 7.6%-13.8% of cases, whereas kinase domain secondary mutation of MET co-occurs in 5%-6%. When co-occurring with METex14, EGFR mutations were often identified as the dominant clone (78%, 7 of 9), whereas when co-occurring, METex14 (39%, 7 of 18) and KRAS (44%, 8 of 18) had similar rates of clonal dominance. PIK3CA and PTEN mutations were almost always subclones (89%, 16 of 18) to METex14. Moreover, RET-CCDC6 fusion and EGFR mutation were detected following crizotinib treatment in two patients, suggesting novel mechanisms of resistance. CONCLUSION: METex14 mutations frequently co-occur with other potential driver oncogenes with differing patterns of clonal dominance observed among the drivers. This cellular context can provide insights into whether METex14 is acting as a primary oncogenic driver or resistance mechanism and help guide treatment choices.

Interpretable deep recommender system model for prediction of kinase inhibitor efficacy across cancer cell lines.

Computational models for drug sensitivity prediction have the potential to significantly improve personalized cancer medicine. Drug sensitivity assays, combined with profiling of cancer cell lines and drugs become increasingly available for training such models. Multiple methods were proposed for predicting drug sensitivity from cancer cell line features, some in a multi-task fashion. So far, no such model leveraged drug inhibition profiles. Importantly, multi-task models require a tailored approach to model interpretability. In this work, we develop DEERS, a neural network recommender system for kinase inhibitor sensitivity prediction. The model utilizes molecular features of the cancer cell lines and kinase inhibition profiles of the drugs. DEERS incorporates two autoencoders to project cell line and drug features into 10-dimensional hidden representations and a feed-forward neural network to combine them into response prediction. We propose a novel interpretability approach, which in addition to the set of modeled features considers also the genes and processes outside of this set. Our approach outperforms simpler matrix factorization models, achieving R [Formula: see text] 0.82 correlation between true and predicted response for the unseen cell lines. The interpretability analysis identifies 67 biological processes that drive the cell line sensitivity to particular compounds. Detailed case studies are shown for PHA-793887, XMD14-99 and Dabrafenib.

Feature selection strategies for drug sensitivity prediction.

Drug sensitivity prediction constitutes one of the main challenges in personalized medicine. Critically, the sensitivity of cancer cells to treatment depends on an unknown subset of a large number of biological features. Here, we compare standard, data-driven feature selection approaches to feature selection driven by prior knowledge of drug targets, target pathways, and gene expression signatures. We asses these methodologies on Genomics of Drug Sensitivity in Cancer (GDSC) dataset, evaluating 2484 unique models. For 23 drugs, better predictive performance is achieved when the features are selected according to prior knowledge of drug targets and pathways. The best correlation of observed and predicted response using the test set is achieved for Linifanib (r = 0.75). Extending the drug-dependent features with gene expression signatures yields the most predictive models for 60 drugs, with the best performing example of Dabrafenib. For many compounds, even a very small subset of drug-related features is highly predictive of drug sensitivity. Small feature sets selected using prior knowledge are more predictive for drugs targeting specific genes and pathways, while models with wider feature sets perform better for drugs affecting general cellular mechanisms. Appropriate feature selection strategies facilitate the development of interpretable models that are indicative for therapy design.

Tepotinib in Non-Small-Cell Lung Cancer with MET Exon 14 Skipping Mutations.

BACKGROUND: A splice-site mutation that results in a loss of transcription of exon 14 in the oncogenic driver MET occurs in 3 to 4% of patients with non-small-cell lung cancer (NSCLC). We evaluated the efficacy and safety of tepotinib, a highly selective MET inhibitor, in this patient population. METHODS: In this open-label, phase 2 study, we administered tepotinib (at a dose of 500 mg) once daily in patients with advanced or metastatic NSCLC with a confirmed MET exon 14 skipping mutation. The primary end point was the objective response by independent review among patients who had undergone at least 9 months of follow-up. The response was also analyzed according to whether the presence of a MET exon 14 skipping mutation was detected on liquid biopsy or tissue biopsy. RESULTS: As of January 1, 2020, a total of 152 patients had received tepotinib, and 99 patients had been followed for at least 9 months. The response rate by independent review was 46% (95% confidence interval [CI], 36 to 57), with a median duration of response of 11.1 months (95% CI, 7.2 to could not be estimated) in the combined-biopsy group. The response rate was 48% (95% CI, 36 to 61) among 66 patients in the liquid-biopsy group and 50% (95% CI, 37 to 63) among 60 patients in the tissue-biopsy group; 27 patients had positive results according to both methods. The investigator-assessed response rate was 56% (95% CI, 45 to 66) and was similar regardless of the previous therapy received for advanced or metastatic disease. Adverse events of grade 3 or higher that were considered by investigators to be related to tepotinib therapy were reported in 28% of the patients, including peripheral edema in 7%. Adverse events led to permanent discontinuation of tepotinib in 11% of the patients. A molecular response, as measured in circulating free DNA, was observed in 67% of the patients with matched liquid-biopsy samples at baseline and during treatment. CONCLUSIONS: Among patients with advanced NSCLC with a confirmed MET exon 14 skipping mutation, the use of tepotinib was associated with a partial response in approximately half the patients. Peripheral edema was the main toxic effect of grade 3 or higher. (Funded by Merck [Darmstadt, Germany]; VISION ClinicalTrials.gov number, NCT02864992.).

Distinct human circulating NKp30+FcεRIγ+CD8+ T cell population exhibiting high natural killer-like antitumor potential.

CD8+ T cells are considered prototypical cells of adaptive immunity. Here, we uncovered a distinct CD8+ T cell population expressing the activating natural killer (NK) receptor NKp30 in the peripheral blood of healthy individuals. We revealed that IL-15 could de novo induce NKp30 expression in a population of CD8+ T cells and drive their differentiation toward a broad innate transcriptional landscape. The adaptor FcεRIγ was concomitantly induced and was shown to be crucial to enable NKp30 cell-surface expression and function in CD8+ T cells. FcεRIγ de novo expression required promoter demethylation and was accompanied by acquisition of the signaling molecule Syk and the "innate" transcription factor PLZF. IL-15-induced NKp30+CD8+ T cells exhibited high NK-like antitumor activity in vitro and were able to synergize with T cell receptor signaling. Importantly, this population potently controlled tumor growth in a preclinical xenograft mouse model. Our study, while blurring the borders between innate and adaptive immunity, reveals a unique NKp30+FcεRIγ+CD8+ T cell population with high antitumor therapeutic potential.

Molecular Classification Substitutes for the Prognostic Variables Stage, Age, and MYCN Status in Neuroblastoma Risk Assessment.

BACKGROUND: Current risk stratification systems for neuroblastoma patients consider clinical, histopathological, and genetic variables, and additional prognostic markers have been proposed in recent years. We here sought to select highly informative covariates in a multistep strategy based on consecutive Cox regression models, resulting in a risk score that integrates hazard ratios of prognostic variables. METHODS: A cohort of 695 neuroblastoma patients was divided into a discovery set (n=75) for multigene predictor generation, a training set (n=411) for risk score development, and a validation set (n=209). Relevant prognostic variables were identified by stepwise multivariable L1-penalized least absolute shrinkage and selection operator (LASSO) Cox regression, followed by backward selection in multivariable Cox regression, and then integrated into a novel risk score. RESULTS: The variables stage, age, MYCN status, and two multigene predictors, NB-th24 and NB-th44, were selected as independent prognostic markers by LASSO Cox regression analysis. Following backward selection, only the multigene predictors were retained in the final model. Integration of these classifiers in a risk scoring system distinguished three patient subgroups that differed substantially in their outcome. The scoring system discriminated patients with diverging outcome in the validation cohort (5-year event-free survival, 84.9±3.4 vs 63.6±14.5 vs 31.0±5.4; P<.001), and its prognostic value was validated by multivariable analysis. CONCLUSION: We here propose a translational strategy for developing risk assessment systems based on hazard ratios of relevant prognostic variables. Our final neuroblastoma risk score comprised two multigene predictors only, supporting the notion that molecular properties of the tumor cells strongly impact clinical courses of neuroblastoma patients.

Genetic Contribution to Alcohol Dependence: Investigation of a Heterogeneous German Sample of Individuals with Alcohol Dependence, Chronic Alcoholic Pancreatitis, and Alcohol-Related Cirrhosis.

The present study investigated the genetic contribution to alcohol dependence (AD) using genome-wide association data from three German samples. These comprised patients with: (i) AD; (ii) chronic alcoholic pancreatitis (ACP); and (iii) alcohol-related liver cirrhosis (ALC). Single marker, gene-based, and pathway analyses were conducted. A significant association was detected for the ADH1B locus in a gene-based approach (puncorrected = 1.2 × 10-6; pcorrected = 0.020). This was driven by the AD subsample. No association with ADH1B was found in the combined ACP + ALC sample. On first inspection, this seems surprising, since ADH1B is a robustly replicated risk gene for AD and may therefore be expected to be associated also with subgroups of AD patients. The negative finding in the ACP + ALC sample, however, may reflect genetic stratification as well as random fluctuation of allele frequencies in the cases and controls, demonstrating the importance of large samples in which the phenotype is well assessed.

A Systematic Approach to Defining the microRNA Landscape in Metastasis.

The microRNA (miRNA) landscape changes during the progression of cancer. We defined a metastasis-associated miRNA landscape using a systematic approach. We profiled and validated miRNA and mRNA expression in a unique series of human colorectal metastasis tissues together with their matched primary tumors and corresponding normal tissues. We identified an exclusive miRNA signature that is differentially expressed in metastases. Three of these miRNAs were identified as key drivers of an EMT-regulating network acting though a number of novel targets. These targets include SIAH1, SETD2, ZEB2, and especially FOXN3, which we demonstrated for the first time as a direct transcriptional suppressor of N-cadherin. The modulation of N-cadherin expression had significant impact on migration, invasion, and metastasis in two different in vivo models. The significant deregulation of the miRNAs defining the network was confirmed in an independent patient set as well as in a database of diverse malignancies derived from more than 6,000 patients. Our data define a novel metastasis-orchestrating network based on systematic hypothesis generation from metastasis tissues.

MYCN amplification confers enhanced folate dependence and methotrexate sensitivity in neuroblastoma.

MYCN amplification occurs in 20% of neuroblastomas and is strongly related to poor clinical outcome. We have identified folate-mediated one-carbon metabolism as highly upregulated in neuroblastoma tumors with MYCN amplification and have validated this finding experimentally by showing that MYCN amplified neuroblastoma cell lines have a higher requirement for folate and are significantly more sensitive to the antifolate methotrexate than cell lines without MYCN amplification. We have demonstrated that methotrexate uptake in neuroblastoma cells is mediated principally by the reduced folate carrier (RFC; SLC19A1), that SLC19A1 and MYCN expression are highly correlated in both patient tumors and cell lines, and that SLC19A1 is a direct transcriptional target of N-Myc. Finally, we assessed the relationship between SLC19A1 expression and patient survival in two independent primary tumor cohorts and found that SLC19A1 expression was associated with increased risk of relapse or death, and that SLC19A1 expression retained prognostic significance independent of age, disease stage and MYCN amplification. This study adds upregulation of folate-mediated one-carbon metabolism to the known consequences of MYCN amplification, and suggests that this pathway might be targeted in poor outcome tumors with MYCN amplification and high SLC19A1 expression.

Revised risk estimation and treatment stratification of low- and intermediate-risk neuroblastoma patients by integrating clinical and molecular prognostic markers.

PURPOSE: To optimize neuroblastoma treatment stratification, we aimed at developing a novel risk estimation system by integrating gene expression-based classification and established prognostic markers. EXPERIMENTAL DESIGN: Gene expression profiles were generated from 709 neuroblastoma specimens using customized 4 × 44 K microarrays. Classification models were built using 75 tumors with contrasting courses of disease. Validation was performed in an independent test set (n = 634) by Kaplan-Meier estimates and Cox regression analyses. RESULTS: The best-performing classifier predicted patient outcome with an accuracy of 0.95 (sensitivity, 0.93; specificity, 0.97) in the validation cohort. The highest potential clinical value of this predictor was observed for current low-risk patients [5-year event-free survival (EFS), 0.84 ± 0.02 vs. 0.29 ± 0.10; 5-year overall survival (OS), 0.99 ± 0.01 vs. 0.76 ± 0.11; both P < 0.001] and intermediate-risk patients (5-year EFS, 0.88 ± 0.06 vs. 0.41 ± 0.10; 5-year OS, 1.0 vs. 0.70 ± 0.09; both P < 0.001). In multivariate Cox regression models for low-risk/intermediate-risk patients, the classifier outperformed risk assessment of the current German trial NB2004 [EFS: hazard ratio (HR), 5.07; 95% confidence interval (CI), 3.20-8.02; OS: HR, 25.54; 95% CI, 8.40-77.66; both P < 0.001]. On the basis of these findings, we propose to integrate the classifier into a revised risk stratification system for low-risk/intermediate-risk patients. According to this system, we identified novel subgroups with poor outcome (5-year EFS, 0.19 ± 0.08; 5-year OS, 0.59 ± 0.1), for whom we propose intensified treatment, and with beneficial outcome (5-year EFS, 0.87 ± 0.05; 5-year OS, 1.0), who may benefit from treatment de-escalation. CONCLUSIONS: Combination of gene expression-based classification and established prognostic markers improves risk estimation of patients with low-risk/intermediate-risk neuroblastoma. We propose to implement our revised treatment stratification system in a prospective clinical trial.

XRCC5 as a risk gene for alcohol dependence: evidence from a genome-wide gene-set-based analysis and follow-up studies in Drosophila and humans.

Genetic factors have as large role as environmental factors in the etiology of alcohol dependence (AD). Although genome-wide association studies (GWAS) enable systematic searches for loci not hitherto implicated in the etiology of AD, many true findings may be missed owing to correction for multiple testing. The aim of the present study was to circumvent this limitation by searching for biological system-level differences, and then following up these findings in humans and animals. Gene-set-based analysis of GWAS data from 1333 cases and 2168 controls identified 19 significantly associated gene-sets, of which 5 could be replicated in an independent sample. Clustered in these gene-sets were novel and previously identified susceptibility genes. The most frequently present gene, ie in 6 out of 19 gene-sets, was X-ray repair complementing defective repair in Chinese hamster cells 5 (XRCC5). Previous human and animal studies have implicated XRCC5 in alcohol sensitivity. This phenotype is inversely correlated with the development of AD, presumably as more alcohol is required to achieve the desired effects. In the present study, the functional role of XRCC5 in AD was further validated in animals and humans. Drosophila mutants with reduced function of Ku80-the homolog of mammalian XRCC5-due to RNAi silencing showed reduced sensitivity to ethanol. In humans with free access to intravenous ethanol self-administration in the laboratory, the maximum achieved blood alcohol concentration was influenced in an allele-dose-dependent manner by genetic variation in XRCC5. In conclusion, our convergent approach identified new candidates and generated independent evidence for the involvement of XRCC5 in alcohol dependence.

Integrated pathway-based approach identifies association between genomic regions at CTCF and CACNB2 and schizophrenia.

In the present study, an integrated hierarchical approach was applied to: (1) identify pathways associated with susceptibility to schizophrenia; (2) detect genes that may be potentially affected in these pathways since they contain an associated polymorphism; and (3) annotate the functional consequences of such single-nucleotide polymorphisms (SNPs) in the affected genes or their regulatory regions. The Global Test was applied to detect schizophrenia-associated pathways using discovery and replication datasets comprising 5,040 and 5,082 individuals of European ancestry, respectively. Information concerning functional gene-sets was retrieved from the Kyoto Encyclopedia of Genes and Genomes, Gene Ontology, and the Molecular Signatures Database. Fourteen of the gene-sets or pathways identified in the discovery dataset were confirmed in the replication dataset. These include functional processes involved in transcriptional regulation and gene expression, synapse organization, cell adhesion, and apoptosis. For two genes, i.e. CTCF and CACNB2, evidence for association with schizophrenia was available (at the gene-level) in both the discovery study and published data from the Psychiatric Genomics Consortium schizophrenia study. Furthermore, these genes mapped to four of the 14 presently identified pathways. Several of the SNPs assigned to CTCF and CACNB2 have potential functional consequences, and a gene in close proximity to CACNB2, i.e. ARL5B, was identified as a potential gene of interest. Application of the present hierarchical approach thus allowed: (1) identification of novel biological gene-sets or pathways with potential involvement in the etiology of schizophrenia, as well as replication of these findings in an independent cohort; (2) detection of genes of interest for future follow-up studies; and (3) the highlighting of novel genes in previously reported candidate regions for schizophrenia.

Cell competition is a tumour suppressor mechanism in the thymus.

Cell competition is an emerging principle underlying selection for cellular fitness during development and disease. Competition may be relevant for cancer, but an experimental link between defects in competition and tumorigenesis is elusive. In the thymus, T lymphocytes develop from precursors that are constantly replaced by bone-marrow-derived progenitors. Here we show that in mice this turnover is regulated by natural cell competition between 'young' bone-marrow-derived and 'old' thymus-resident progenitors that, although genetically identical, execute differential gene expression programs. Disruption of cell competition leads to progenitor self-renewal, upregulation of Hmga1, transformation, and T-cell acute lymphoblastic leukaemia (T-ALL) resembling the human disease in pathology, genomic lesions, leukaemia-associated transcripts, and activating mutations in Notch1. Hence, cell competition is a tumour suppressor mechanism in the thymus. Failure to select fit progenitors through cell competition may explain leukaemia in X-linked severe combined immune deficiency patients who showed thymus-autonomous T-cell development after therapy with gene-corrected autologous progenitors.

Genome-wide association-, replication-, and neuroimaging study implicates HOMER1 in the etiology of major depression.

BACKGROUND: Genome-wide association studies are a powerful tool for unravelling the genetic background of complex disorders such as major depression. METHODS: We conducted a genome-wide association study of 604 patients with major depression and 1364 population based control subjects. The top hundred findings were followed up in a replication sample of 409 patients and 541 control subjects. RESULTS: Two SNPs showed nominally significant association in both the genome-wide association study and the replication samples: 1) rs9943849 (p(combined) = 3.24E-6) located upstream of the carboxypeptidase M (CPM) gene and 2) rs7713917 (p(combined) = 1.48E-6), located in a putative regulatory region of HOMER1. Further evidence for HOMER1 was obtained through gene-wide analysis while conditioning on the genotypes of rs7713917 (p(combined) = 4.12E-3). Homer1 knockout mice display behavioral traits that are paradigmatic of depression, and transcriptional variants of Homer1 result in the dysregulation of cortical-limbic circuitry. This is consistent with the findings of our subsequent human imaging genetics study, which revealed that variation in single nucleotide polymorphism rs7713917 had a significant influence on prefrontal activity during executive cognition and anticipation of reward. CONCLUSION: Our findings, combined with evidence from preclinical and animal studies, suggest that HOMER1 plays a role in the etiology of major depression and that the genetic variation affects depression via the dysregulation of cognitive and motivational processes.

Prognostic impact of gene expression-based classification for neuroblastoma.

PURPOSE: To evaluate the impact of a predefined gene expression-based classifier for clinical risk estimation and cytotoxic treatment decision making in neuroblastoma patients. PATIENTS AND METHODS: Gene expression profiles of 440 internationally collected neuroblastoma specimens were investigated by microarray analysis, 125 of which were examined prospectively. Patients were classified as either favorable or unfavorable by a 144-gene prediction analysis for microarrays (PAM) classifier established previously on a separate set of 77 patients. PAM classification results were compared with those of current prognostic markers and risk estimation strategies. RESULTS: The PAM classifier reliably distinguished patients with contrasting clinical courses (favorable [n = 249] and unfavorable [n = 191]; 5-year event free survival [EFS] 0.84 +/- 0.03 v 0.38 +/- 0.04; 5-year overall survival [OS] 0.98 +/- 0.01 v 0.56 +/- 0.05, respectively; both P < .001). Moreover, patients with divergent outcome were robustly discriminated in both German and international cohorts and in prospectively analyzed samples (P

Detection and quantification of the nifH gene in shoot and root of cucumber plants.

A real-time polymerase chain reaction (PCR) method was applied to quantify the nifH gene pool in cucumber shoot and root and to evaluate how nitrogen (N) supply and plant age affect the nifH gene pool. In shoots, the relative abundance of the nifH gene was affected neither by different stages of plant growth nor by N supply. In roots, higher numbers of diazotrophic bacteria were found compared with that in the shoot. The nifH gene pool in roots significantly increased with plant age, and unexpectedly, the pool size was positively correlated with N supply. The relative abundance of nifH gene copy numbers in roots was also positively correlated (r = 0.96) with total N uptake of the plant. The data suggest that real-time PCR-based nifH gene quantification in combination with N-content analysis can be used as an efficient way to perform further studies to evaluate the direct contribution of the N2-fixing plant-colonizing plant growth promoting bacteria to plant N nutrition.