Adenocarcinoma in situ and minimally invasive adenocarcinoma are the pre-invasive forms of lung adenocarcinoma. The genomic and immune profiles of these lesions are poorly understood. Here we report exome and transcriptome sequencing of 98 lung adenocarcinoma precursor lesions and 99 invasive adenocarcinomas. We have identified EGFR, RBM10, BRAF, ERBB2, TP53, KRAS, MAP2K1 and MET as significantly mutated genes in the pre/minimally invasive group. Classes of genome alterations that increase in frequency during the progression to malignancy are revealed. These include mutations in TP53, arm-level copy number alterations, and HLA loss of heterozygosity. Immune infiltration is correlated with copy number alterations of chromosome arm 6p, suggesting a link between arm-level events and the tumor immune environment.

Lung cancer has become the most common cancer type and caused the most cancer deaths. Lung adenocarcinoma (LUAD) is one of the major type of lung cancer. This study aimed to establish a signature based on immune related genes that can predict patients’ OS for LUAD. The expression data of 976 LUAD patients from The Cancer Genome Atlas database (training set) and the Gene Expression Omnibus database (four testing sets) and 1534 immune related genes from the ImmPort database were used for generation and validation of the signature. The glmnet Cox proportional hazards model was used to find the best gene model and construct the signature. To assess the independently prognostic ability of the signature, the Kaplan–Meier survival analysis and Cox’s proportional hazards model were performed. A gene model consisting of 30 immune related genes with the highest frequency after 1000 iterations was used as our signature. The signature demonstrated robust prognostic ability in both training set and testing set and could serve as an independent predictor for LUAD patients in all datasets except GSE31210. Besides, the signature could predict the overall survival (OS) of LUAD patients in different subgroups. And this signature was strongly associated with important clinicopathological factors like recurrence and TNM stage. More importantly, patients with high risk score presented high tumor mutation burden. This signature could predict prognosis and reflect the tumor immune microenvironment of LUAD patients, which can promote individualized treatment and provide potential novel targets for immunotherapy.

Characterization and integration of the genome, epigenome, transcriptome, proteome and metabolome of different datasets is difficult owing to a lack of ground truth. Here we develop and characterize suites of publicly available multi-omics reference materials of matched DNA, RNA, protein and metabolites derived from immortalized cell lines from a family quartet of parents and monozygotic twin daughters. These references provide built-in truth defined by relationships among the family members and the information flow from DNA to RNA to protein. We demonstrate how using a ratio-based profiling approach that scales the absolute feature values of a study sample relative to those of a concurrently measured common reference sample produces reproducible and comparable data suitable for integration across batches, labs, platforms and omics types. Our study identifies reference-free ‘absolute’ feature quantification as the root cause of irreproducibility in multi-omics measurement and data integration and establishes the advantages of ratio-based multi-omics profiling with common reference materials.

The kidney renal papillary cell carcinoma (KIRP) is a relatively rare type of kidney cancer. There has been no investigation to find a robust signature to predict the survival outcome of KIRP patients in the aspect of tumor immunology. In this study, 285 KIRP samples from The Cancer Genome Atlas (TCGA) were randomly divided into training and testing set. A total of 1534 immune‐related genes from The Immunology Database and Analysis Portal (ImmPort) were used as candidates to construct the signature. Using univariate Cox analysis, we evaluated the relationship between overall survival and immune‐related genes expression and found 272 immune‐related genes with predicting prognostic ability. In order to construct an efficient predictive model, the Cox proportional hazards model with an elastic‐net penalty was used and identified 23 groups after 1000 iterations. As a result, 15‐genes model showing more stable than other gene groups was chosen to construct our immune‐related risk signature. In line with our expectations, the signature can independently predict the survival outcome of KIRP patients. Patients with high‐immune risk were found correlated with advanced stage. We also found that the high‐immune risk patients with higher PBRM1 and SETD2 mutations, increasing chromosomal instability, together with the gene set enrichment analysis (GSEA) results showing intensive connection of our signature with immune pathways. In conclusion, our study constructs a robust 15‐gene signature for predicting KIRP patients’ survival outcome on the basis of tumor immune environment and may provide possible relationship between prognosis and immune‐related biological function.

Certified RNA reference materials are indispensable for assessing the reliability of RNA sequencing to detect intrinsically small biological differences in clinical settings, such as molecular subtyping of diseases. As part of the Quartet Project for quality control and data integration of multi-omics profiling, we established four RNA reference materials derived from immortalized B-lymphoblastoid cell lines from four members of a monozygotic twin family. Additionally, we constructed ratio-based transcriptome-wide reference datasets between two samples, providing cross-platform and cross-laboratory ‘ground truth’. Investigation of the intrinsically subtle biological differences among the Quartet samples enables sensitive assessment of cross-batch integration of transcriptomic measurements at the ratio level. The Quartet RNA reference materials, combined with the ratio-based reference datasets, can serve as unique resources for assessing and improving the quality of transcriptomic data in clinical and biological settings.