Development of Genetic Markers in Eucalyptus Species by Target Enrichment and Exome Sequencing

The advent of next-generation sequencing has facilitated large-scale discovery, validation and assessment of genetic markers for high density genotyping. The present study was undertaken to identify markers in genes supposedly related to wood property traits in three Eucalyptus species. Ninety four genes involved in xylogenesis were selected for hybridization probe based nuclear genomic DNA target enrichment and exome sequencing. Genomic DNA was isolated from the leaf tissues and used for on-array probe hybridization followed by Illumina sequencing. The raw sequence reads were trimmed and high-quality reads were mapped to the E. grandis reference sequence and the presence of single nucleotide variants (SNVs) and insertions/ deletions (InDels) were identified across the three species. The average read coverage was 216X and a total of 2294 SNVs and 479 InDels were discovered in E. camaldulensis, 2383 SNVs and 518 InDels in E. tereticornis, and 1228 SNVs and 409 InDels in E. grandis. Additionally, SNV calling and InDel detection were conducted in pair-wise comparisons of E. tereticornis vs. E. grandis, E. camaldulensis vs. E. tereticornis and E. camaldulensis vs. E. grandis. This study presents an efficient and high throughput method on development of genetic markers for family– based QTL and association analysis in Eucalyptus.     

Comprehensive characterization of somatic variants associated with intronic polyadenylation in human cancers

Somatic single nucleotide variants (SNVs) in cancer genome affect gene expression through various mechanisms depending on their genomic location. While somatic SNVs near canonical splice sites have been reported to cause abnormal splicing of cancer-related genes, whether these SNVs can affect gene expression through other mechanisms remains an open question. Here, we analyzed RNA sequencing and exome data from 4,998 cancer patients covering ten cancer types and identified 152 somatic SNVs near splice sites that were associated with abnormal intronic polyadenylation (IPA). IPA-associated somatic variants favored the localization near the donor splice sites compared to the acceptor splice sites. A proportion of SNV-associated IPA events overlapped with premature cleavage and polyadenylation events triggered by U1 small nuclear ribonucleoproteins (snRNP) inhibition. GC content, intron length and polyadenylation signal were three genomic features that differentiated between SNV-associated IPA and intron retention. Notably, IPA-associated SNVs were enriched in tumor suppressor genes (TSGs), including the well-known TSGs such as PTEN and CDH1 with recurrent SNV-associated IPA events. Minigene assay confirmed that SNVs from PTEN, CDH1, VEGFA, GRHL2, CUL3 and WWC2 could lead to IPA. This work reveals that IPA acts as a novel mechanism explaining the functional consequence of somatic SNVs in human cancer.     

Common Genomic Aberrations in Mouse and Human Breast Cancers with Concurrent P53 Deficiency and Activated PTEN-PI3K-AKT Pathway

Simultaneous P53 loss and activation of the PTEN-restricted PI3K-AKT pathway frequently occur in aggressive breast cancers. P53 loss causes genome instability, while PTEN loss and/or activating mutations of PIK3CA and AKT promote cancer cell proliferation that also increases incidences of genomic aberrations. However, the genomic alterations associated with P53 loss and activated PTEN-PI3K-AKT signaling in breast cancer have not been defined. Spatiotemporally controlled breast cancer models with inactivation of both P53 and Pten in adult mice have not been established for studying genomic alterations. Herein, we deleted both floxed Pten and Tp53 genes in the mammary gland epithelial cells in adult mice using a RCAS virus-mediated Cre-expressing system. These mice developed small tumors in 21 weeks, and poorly differentiated larger tumors in 26 weeks. In these tumors, we identified 360 genes mutated by nonsynonymous point mutations and small insertions and deletions (NSPMs/InDels), 435 genes altered by copy number amplifications (CNAs), and 450 genes inactivated by copy number deletions (CNDs). Importantly, 22.2%, 75.9% and 27.3% of these genes were also altered in human breast tumors with P53 and PTEN losses or P53 loss and activated PI3K-AKT signaling by NSPMs/InDels, CNAs and CNDs, respectively. Therefore, inactivation of P53 and Pten in adult mice causes rapid-growing breast tumors, and these tumors recapitulate a significant number of genetic aberrations in human breast tumors with inactivated P53 and activated PTEN-PI3K-AKT signaling. Further characterization of these commonly altered genes in breast cancer should help to identify novel cancer-driving genes and molecular targets for developing therapeutics.     

Contribution of Multiparameter Flow Cytometry Immunophenotyping to the Diagnostic Screening and Classification of Pediatric Cancer

Pediatric cancer is a relatively rare and heterogeneous group of hematological and non-hematological malignancies which require multiple procedures for its diagnostic screening and classification. Until now, flow cytometry (FC) has not been systematically applied to the diagnostic work-up of such malignancies, particularly for solid tumors. Here we evaluated a FC panel of markers for the diagnostic screening of pediatric cancer and further classification of pediatric solid tumors. The proposed strategy aims at the differential diagnosis between tumoral vs. reactive samples, and hematological vs. non-hematological malignancies, and the subclassification of solid tumors. In total, 52 samples from 40 patients suspicious of containing tumor cells were analyzed by FC in parallel to conventional diagnostic procedures. The overall concordance rate between both approaches was of 96% (50/52 diagnostic samples), with 100% agreement for all reactive/inflammatory and non-infiltrated samples as well as for those corresponding to solid tumors (n = 35), with only two false negative cases diagnosed with Hodgkin lymphoma and anaplastic lymphoma, respectively. Moreover, clear discrimination between samples infiltrated by hematopoietic vs. non-hematopoietic tumor cells was systematically achieved. Distinct subtypes of solid tumors showed different protein expression profiles, allowing for the differential diagnosis of neuroblastoma (CD56^hi/GD2⁺/CD81^hi), primitive neuroectodermal tumors (CD271^hi/CD99⁺), Wilms tumors (>1 cell population), rhabdomyosarcoma (_nuMYOD1⁺/_numyogenin⁺), carcinomas (CD45⁻/EpCAM⁺), germ cell tumors (CD56⁺/CD45⁻/NG2⁺/CD10⁺) and eventually also hemangiopericytomas (CD45⁻/CD34⁺). In summary, our results show that multiparameter FC provides fast and useful complementary data to routine histopathology for the diagnostic screening and classification of pediatric cancer.     

Comparative Oncogenomics Implicates the Neurofibromin 1 Gene (NF1) as a Breast Cancer Driver

Identifying genomic alterations driving breast cancer is complicated by tumor diversity and genetic heterogeneity. Relevant mouse models are powerful for untangling this problem because such heterogeneity can be controlled. Inbred Chaos3 mice exhibit high levels of genomic instability leading to mammary tumors that have tumor gene expression profiles closely resembling mature human mammary luminal cell signatures. We genomically characterized mammary adenocarcinomas from these mice to identify cancer-causing genomic events that overlap common alterations in human breast cancer. Chaos3 tumors underwent recurrent copy number alterations (CNAs), particularly deletion of the RAS inhibitor Neurofibromin 1 (Nf1) in nearly all cases. These overlap with human CNAs including NF1, which is deleted or mutated in 27.7% of all breast carcinomas. Chaos3 mammary tumor cells exhibit RAS hyperactivation and increased sensitivity to RAS pathway inhibitors. These results indicate that spontaneous NF1 loss can drive breast cancer. This should be informative for treatment of the significant fraction of patients whose tumors bear NF1 mutations.     

Cancer progression by non-clonal chromosome aberrations

The establishment of the correct conceptual framework is vital to any scientific discipline including cancer research. Influenced by hematologic cancer studies, the current cancer concept focuses on the stepwise patterns of progression as defined by specific recurrent genetic aberrations. This concept has faced a tough challenge as the majority of cancer cases follow non-linear patterns and display stochastic progression. In light of the recent discovery that genomic instability is directly linked to stochastic non-clonal chromosome aberrations (NCCAs), and that cancer progression can be characterized as a dynamic relationship between NCCAs and recurrent clonal chromosome aberrations (CCAs), we propose that the dynamics of NCCAs is a key element for karyotypic evolution in solid tumors. To support this viewpoint, we briefly discuss various basic elements responsible for cancer initiation and progression within an evolutionary context. We argue that even though stochastic changes can be detected at various levels of genetic organization, such as at the gene level and epigenetic level, it is primarily detected at the chromosomal or genome level. Thus, NCCA-mediated genomic variation plays a dominant role in cancer progression. To further illustrate the involvement of NCCA/CCA cycles in the pattern of cancer evolution, four cancer evolutionary models have been proposed based on the comparative analysis of karyotype patterns of various types of cancer.     

Intratumor heterogeneity and clonal evolution revealed in castration-resistant prostate cancer by longitudinal genomic analysis

Intratumor heterogeneity is a key driver for local relapse and treatment failure. Thus, using multifocal prostate cancer as a model to investigate tumor inter-clonal relationships and tumor evolution could aid in our understanding of drug resistance. Previous studies discovered genomic alterations by comparing hormone-sensitive prostate cancer (HSPC) with castration-resistant prostate cancer (CRPC) in large cohorts. However, most studies did not sequentially sample tumors from the same patient. In our study, we performed whole-exome sequencing (WES) on 14 specimens from five locally relapsed patients before and after androgen-deprivation therapy. We described the landscape of genomic alterations before and after treatment and identified critical driver events that could have contributed to the evolution of CRPC. In addition to confirming known cancer genes such as TP53 and CDK12, we also identified new candidate genes that may play a role in the progression of prostate cancer, including MYO15A, CHD6 and LZTR1. At copy number alteration (CNA) level, gain of 8q24.13-8q24.3 was observed in 60% of patients and was the most commonly altered locus in both HSPC and CRPC tumors. Finally, utilizing phylogenetic reconstruction, we explored the clonal progression pattern from HSPC to CRPC in each patient. Our findings highlight the complex and heterogeneous mechanisms underlying the development of drug resistance, and underscore the potential value of monitoring tumor clonal architectures during disease progression in a clinical setting.     

Stochastic cancer progression driven by non-clonal chromosome aberrations

Cancer research has previously focused on the identification of specific genes and pathways responsible for cancer initiation and progression based on the prevailing viewpoint that cancer is caused by a stepwise accumulation of genetic aberrations. This viewpoint, however, is not consistent with the clinical finding that tumors display high levels of genetic heterogeneity and distinctive karyotypes. We show that chromosomal instability primarily generates stochastic karyotypic changes leading to the random progression of cancer. This was accomplished by tracing karyotypic patterns of individual cells that contained either defective genes responsible for genome integrity or were challenged by onco-proteins or carcinogens that destabilized the genome. Analysis included the tracing of patterns of karyotypic evolution during different stages of cellular immortalization. This study revealed that non-clonal chromosomal aberrations (NCCAs) (both aneuploidy and structural aberrations) and not recurrent clonal chromosomal aberrations (CCAs) are directly linked to genomic instability and karyotypic evolution. Discovery of "transitional CCAs" during in vitro immortalization clearly demonstrates that karyotypic evolution in solid tumors is not a continuous process. NCCAs and their dynamic interplay with CCAs create infinite genomic combinations leading to clonal diversity necessary for cancer cell evolution. The karyotypic chaos observed within the cell crisis stage prior to establishment of the immortalization further supports the ultimate importance of genetic aberrations at the karyotypic or genome level. Therefore, genomic instability generated NCCAs are a key driving force in cancer progression. The dynamic relationship between NCCAs and CCAs provides a mechanism underlying chromosomal based cancer evolution and could have broad clinical applications.     

Exonic variants in multiple myeloma patients associated with relapsed/ refractory and response to bortezomib regimens

Novel treatment in multiple myeloma represented by proteasome inhibitors, immunomodulatory drugs and monoclonal antibodies have produced a deep response. However, relapses are possible, and all classes of drugs are refractory to patients. Next-generation sequencing has improved our understanding of the multiple myeloma genome related to drug resistance and has discovered many genomic variants. Therefore, this study was conducted to investigate new variants associated with drug resistance in MM patients who relapsed and refractory to bortezomib regimen and daratumumab treatment using next-generation sequencing for whole-exome sequencing. Peripheral blood samples were collected in EDTA tubes from six patients; four were in relapsed and refractory to bortezomib regimens and daratumumab; two patients responded to bortezomib regimens. Whole-exome sequencing was performed by the MGI-DNBSEQ-G400 instrument. We identified 21 variants in multiple myeloma patients. Seventeen variants were found in relapsed and refractory multiple myeloma in 11 genes (GNAQ, PMS1, CREB1, NSUNS2, PIK3CG, ROS1, PMS2, FIT4, KDM5A, STK11 and ZFHX3). And four variants were identified in two patients with response to bortezomib regimens in 4 genes (RAF1, CREB1, ZFHX3 and INSR). We have observed several genetic variants in many genes that may have been associated with the poor prognosis and poor response to treatment in these patients. These values should be further confirmed in large sample studies using the RNA-seq technique to identify genome expression.     

A Phase I Trial of DFMO Targeting Polyamine Addiction in Patients with Relapsed/Refractory Neuroblastoma

BackgroundNeuroblastoma (NB) is the most common cancer in infancy and most frequent cause of death from extracranial solid tumors in children. Ornithine decarboxylase (ODC) expression is an independent indicator of poor prognosis in NB patients. This study investigated safety, response, pharmacokinetics, genetic and metabolic factors associated with ODC in a clinical trial of the ODC inhibitor difluoromethylornithine (DFMO) ± etoposide for patients with relapsed or refractory NB.ConclusionsDFMO doses of 500-1500mg/m2/day are safe and well tolerated in children with relapsed NB. Children with the minor T allele at rs2302616 of the ODC gene with relapsed or refractory NB had higher levels of urinary polyamine markers and responded better to therapy containing DFMO, compared to those with the major G allele at this locus. These findings suggest that this patient subset may display dependence on polyamines and be uniquely susceptible to therapies targeting this pathway.

Trial Registration

Clinicaltrials.gov NCT#01059071     

From humble beginnings to success in the clinic: Chimeric antigen receptor-modified T-cells and implications for immunotherapy

In the past 50 years, disease burden has steadily shifted from infectious disease to cancer. Standard chemotherapy has long been the mainstay of cancer medical management, and despite vast efforts towards more targeted and personalized drug therapy, many cancers remain refractory to treatment, with high rates of relapse and poor prognosis. Recent dramatic immunotherapy clinical trials have demonstrated that engineering T-cells with chimeric antigen receptors (CARs) to target CD19 can lead to complete remission in relapsed or refractory B-cell malignancies, generating a great deal of enthusiasm in the field. Here we provide a comprehensive overview of the history of adoptive T-cell therapy, including CARs, in solid tumors as well as hematologic malignancies. CAR therapy has the potential to fundamentally transform cancer treatment with specific and even personalized targeting of tissue- and tumor-specific antigens. However, before CARs become standard first-line treatment modalities, critical issues regarding efficacy, combinatorial regimens, and mechanisms of treatment failure and toxicity will need to be addressed.     

Recurrent Chromosome 22 Deletions in Osteoblastoma Affect Inhibitors of the Wnt/Beta-Catenin Signaling Pathway

Osteoblastoma is a bone forming tumor with histological features highly similar to osteoid osteoma; the discrimination between the tumor types is based on size and growth pattern. The vast majority of osteoblastomas are benign but there is a group of so-called aggressive osteoblastomas that can be diagnostically challenging at the histopathological level. The genetic aberrations required for osteoblastoma development are not known and no genetic difference between conventional and aggressive osteoblastoma has been reported. In order to identify recurrent genomic aberrations of importance for tumor development we applied cytogenetic and/or SNP array analyses on nine conventional and two aggressive osteoblastomas. The conventional osteoblastomas showed few or no acquired genetic aberrations while the aggressive tumors displayed heavily rearranged genomes. In one of the aggressive osteoblastomas, three neighboring regions in chromosome band 22q12 were homozygously deleted. Hemizygous deletions of these regions were found in two additional cases, one aggressive and one conventional. In total, 10 genes were recurrently and homozygously lost in osteoblastoma. Four of them are functionally involved in regulating osteogenesis and/or tumorigenesis. MN1 and NF2 have previously been implicated in the development of leukemia and solid tumors, and ZNRF3 and KREMEN1 are inhibitors of the Wnt/beta-catenin signaling pathway. In line with deletions of the latter two genes, high beta-catenin protein expression has previously been reported in osteoblastoma and aberrations affecting the Wnt/beta-catenin pathway have been found in other bone lesions, including osteoma and osteosarcoma.     

Clinical relevance of cancer stem cell chemotherapeutic assay for recurrent ovarian cancer

IntroductionDisease recurrence and progression of ovarian cancer is common with the development of platinum-resistant or refractory disease. This is due in large part to the presence of chemo-resistant cancer stem cells (CSCs) that contribute to tumor propagation, maintenance, and treatment resistance. We developed a CSCs drug cytotoxicity assay (ChemoID) to identify the most effective chemotherapy treatment from a panel of FDA approved chemotherapies.MethodsAscites and pleural fluid samples were collected under physician order from 45 consecutive patients affected by 3rd-5th relapsed ovarian cancer. Test results from the assay were used to treat patients with the highest cell kill drugs, taking into consideration their health status and using dose reductions, as needed. A retrospective chart review of CT and PET scans was used to determine patients' outcomes for tumor response, time to recurrence, progression-free survival (PFS), and overall survival (OS).ResultsWe observed that recurrent ovarian cancer patients treated with high-cell kill chemotherapy agents guided by the CSCs drug response assay had an improvement in the median PFS corresponding to 5.4 months (3rd relapse), 3.6 months (4th relapse), and 3.9 months (5th relapse) when compared to historical data. Additionally, we observed that ovarian cancer patients identified as non-responders by the CSC drug response assay had 30 times the hazard of death compared to those women that were identified as responders with respective median survivals of 6 months vs. 13 months. We also found that ChemoID treated patients on average had an incremental cost-effectiveness ratio (ICER) between -$18,421 and $7,241 per life-year saved (LYS).ConclusionsThis study demonstrated improved PFS and OS for recurrent ovarian cancer patients treated with assay-guided chemotherapies while decreasing the cost of treatment.     

A transcriptomic computational analysis of mastic oil-treated Lewis lung carcinomas reveals molecular mechanisms targeting tumor cell growth and survival

Background

DNA copy number alterations are frequently observed in ovarian cancer, but it remains a challenge to identify the most relevant alterations and the specific causal genes in those regions.

Methods

We obtained high-resolution 500K SNP array data for 52 ovarian tumors and identified the most statistically significant minimal genomic regions with the most prevalent and highest-level copy number alterations (recurrent CNAs). Within a region of recurrent CNA, comparison of expression levels in tumors with a given CNA to tumors lacking that CNA and to whole normal ovary samples was used to select genes with CNA-specific expression patterns. A public expression array data set of laser capture micro-dissected (LCM) non-malignant fallopian tube epithelia and LCM ovarian serous adenocarcinoma was used to evaluate the effect of cell-type mixture biases.

Results

Fourteen recurrent deletions were detected on chromosomes 4, 6, 9, 12, 13, 15, 16, 17, 18, 22 and most prevalently on X and 8. Copy number and expression data suggest several apoptosis mediators as candidate drivers of the 8p deletions. Sixteen recurrent gains were identified on chromosomes 1, 2, 3, 5, 8, 10, 12, 15, 17, 19, and 20, with the most prevalent gains localized to 8q and 3q. Within the 8q amplicon, PVT1, but not MYC, was strongly over-expressed relative to tumors lacking this CNA and showed over-expression relative to normal ovary. Likewise, the cell polarity regulators PRKCI and ECT2 were identified as putative drivers of two distinct amplicons on 3q. Co-occurrence analyses suggested potential synergistic or antagonistic relationships between recurrent CNAs. Genes within regions of recurrent CNA showed an enrichment of Cancer Census genes, particularly when filtered for CNA-specific expression.

Conclusion

These analyses provide detailed views of ovarian cancer genomic changes and highlight the benefits of using multiple reference sample types for the evaluation of CNA-specific expression changes.     

Uncovering Genomic Features and Maternal Origin of Korean Native Chicken by Whole Genome Sequencing

The Korean Native Chicken (KNC) is an important endemic biological resource in Korea. While numerous studies have been conducted exploring this breed, none have used next-generation sequencing to identify its specific genomic features. We sequenced five strains of KNC and identified 10.9 million SNVs and 1.3 million InDels. Through the analysis, we found that the highly variable region common to all 5 strains had genes like PCHD15, CISD1, PIK3C2A, and NUCB2 that might be related to the phenotypic traits of the chicken such as auditory sense, growth rate and egg traits. In addition, we assembled unaligned reads that could not be mapped to the reference genome. By assembling the unaligned reads, we were able to present genomic sequences characteristic to the KNC. Based on this, we also identified genes related to the olfactory receptors and antigen that are common to all 5 strains. Finally, through the reconstructed mitochondrial genome sequences, we performed phylogenomic analysis and elucidated the maternal origin of the artificially restored KNC. Our results revealed that the KNC has multiple maternal origins which are in agreement with Korea''s history of chicken breed imports. The results presented here provide a valuable basis for future research on genomic features of KNC and further understanding of KNC''s origin.     

Covariation of viral recombination with single nucleotide variants during virus evolution revealed by CoVaMa

Adaptation of viruses to their environments occurs through the acquisition of both novel single-nucleotide variants (SNV) and recombination events including insertions, deletions, and duplications. The co-occurrence of SNVs in individual viral genomes during their evolution has been well-described. However, unlike covariation of SNVs, studying the correlation between recombination events with each other or with SNVs has been hampered by their inherent genetic complexity and a lack of bioinformatic tools. Here, we expanded our previously reported CoVaMa pipeline (v0.1) to measure linkage disequilibrium between recombination events and SNVs within both short-read and long-read sequencing datasets. We demonstrate this approach using long-read nanopore sequencing data acquired from Flock House virus (FHV) serially passaged in vitro. We found SNVs that were either correlated or anti-correlated with large genomic deletions generated by nonhomologous recombination that give rise to Defective-RNAs. We also analyzed NGS data from longitudinal HIV samples derived from a patient undergoing antiretroviral therapy who proceeded to virological failure. We found correlations between insertions in the p6^Gag and mutations in Gag cleavage sites. This report confirms previous findings and provides insights on novel associations between SNVs and specific recombination events within the viral genome and their role in viral evolution.     

Germline Variation Controls the Architecture of Somatic Alterations in Tumors

Studies have suggested that somatic events in tumors can depend on an individual''s constitutional genotype. We used squamous cell carcinomas (SCC) of the skin, which arise in high multiplicity in organ transplant recipients, as a model to compare the pattern of somatic alterations within and across individuals. Specifically, we performed array comparative genomic hybridization on 104 tumors from 25 unrelated individuals who each had three or more independently arisen SCCs and compared the profiles occurring within patients to profiles of tumors across a larger set of 135 patients. In general, chromosomal aberrations in SCCs were more similar within than across individuals (two-sided exact-test p-value ), consistent with the notion that the genetic background was affecting the pattern of somatic changes. To further test this possibility, we performed allele-specific imbalance studies using microsatellite markers mapping to 14 frequently aberrant regions of multiple independent tumors from 65 patients. We identified nine loci which show evidence of preferential allelic imbalance. One of these loci, 8q24, corresponded to a region in which multiple single nucleotide polymorphisms have been associated with increased cancer risk in genome-wide association studies (GWAS). We tested three implicated variants and identified one, rs13281615, with evidence of allele-specific imbalance (p-value = 0.012). The finding of an independently identified cancer susceptibility allele with allele-specific imbalance in a genomic region affected by recurrent DNA copy number changes suggest that it may also harbor risk alleles for SCC. Together these data provide strong evidence that the genetic background is a key driver of somatic events in cancer, opening an opportunity to expand this approach to identify cancer risk alleles.     

Utility of Genomic Analysis in Differentiating Synchronous and Metachronous Lung Adenocarcinomas from Primary Adenocarcinomas with Intrapulmonary Metastasis

Distinguishing synchronous and metachronous primary lung adenocarcinomas from adenocarcinomas with intrapulmonary metastasis is essential for optimal patient management. In this study, multiple lung adenocarcinomas occurring in the same patient were evaluated using comprehensive histopathologic evaluation supplemented with molecular analysis. The cohort included 18 patients with a total of 52 lung adenocarcinomas. Eleven patients had a new diagnosis of multiple adenocarcinomas in the same lobe (n = 5) or different lobe (n = 6). Seven patients had a history of lung cancer and developed multiple new tumors. The final diagnosis was made in resection specimens (n = 49), fine needle aspiration (n = 2), and biopsy (n = 1). Adenocarcinomas were non‐mucinous, and histopathologic comparison of tumors was performed. All tumors save for one were subjected to ALK gene rearrangement testing and targeted Next Generation Sequencing (NGS). Using clinical, radiologic, and morphologic features, a confident conclusion favoring synchronous/metachronous or metastatic disease was made in 65% of patients. Cases that proved challenging included ones with more than three tumors showing overlapping growth patterns and lacking a predominant lepidic component. Genomic signatures unique to each tumor were helpful in determining the relationship of multiple carcinomas in 72% of patients. Collectively, morphologic and genomic data proved to be of greater value and achieved a conclusive diagnosis in 94% of patients. Assessment of the genomic profiles of multiple lung adenocarcinomas complements the histological findings, enabling a more comprehensive assessment of synchronous, metachronous, and metastatic lesions in most patients, thereby improving staging accuracy. Targeted NGS can identify genetic alterations with therapeutic implications.