Malignant pleural mesothelioma is a disease primarily associated with exposure to the carcinogen asbestos. Whereas other carcinogen-related tumors are associated with a high tumor mutation burden, mesothelioma is not. We sought to resolve this discrepancy. We used mate-pair (n = 22), RNA (n = 28), and T cell receptor sequencing along with in silico predictions and immunologic assays to understand how structural variants of chromosomes affect the transcriptome. We observed that inter- or intrachromosomal rearrangements were present in every specimen and were frequently in a pattern of chromoanagenesis such as chromoplexy or chromothripsis. Transcription of rearrangement-related junctions was predicted to result in many potential neoantigens, some of which were proven to bind patient-specific major histocompatibility complex molecules and to expand intratumoral T cell clones. T cells responsive to these predicted neoantigens were also present in a patient's circulating T cell repertoire. Analysis of genomic array data from the mesothelioma cohort in The Cancer Genome Atlas suggested that multiple chromothriptic-like events negatively impact survival. Our findings represent the discovery of potential neoantigen expression driven by structural chromosomal rearrangements. These results may have implications for the development of novel immunotherapeutic strategies and the selection of patients to receive immunotherapies.

The COVID-19 pandemic demands assimilation of all biomedical knowledge to decode mechanisms of pathogenesis. Despite the recent renaissance in neural networks, a platform for the real-time synthesis of the exponentially growing biomedical literature and deep omics insights is unavailable. Here, we present the nferX platform for dynamic inference from over 45 quadrillion possible conceptual associations from unstructured text, and triangulation with insights from single-cell RNA-sequencing, bulk RNA-seq and proteomics from diverse tissue types. A hypothesis-free profiling of ACE2 suggests tongue keratinocytes, olfactory epithelial cells, airway club cells and respiratory ciliated cells as potential reservoirs of the SARS-CoV-2 receptor. We find the gut as the putative hotspot of COVID-19, where a maturation correlated transcriptional signature is shared in small intestine enterocytes among coronavirus receptors (ACE2, DPP4, ANPEP). A holistic data science platform triangulating insights from structured and unstructured data holds potential for accelerating the generation of impactful biological insights and hypotheses.

SUMMARY A long-standing concept in vision science has held that a single photoreceptor expresses a single type of opsin, the protein component of visual pigment. However, the number of examples in the literature of photoreceptors from vertebrates and invertebrates that break this rule is increasing. Here, we describe a newly discovered Limulus opsin, Limulus opsin5, which is significantly different from previously characterized Limulus opsins, opsins1 and 2. We show that opsin5 is co-expressed with opsins1 and 2 in Limulus lateral and ventral eye photoreceptors and provide the first evidence that the expression of co-expressed opsins can be differentially regulated. We show that the relative levels of opsin5 and opsin1 and 2 in the rhabdom change with a diurnal rhythm and that their relative levels are also influenced by the animal9s central circadian clock. An analysis of the sequence of opsin5 suggests it is sensitive to visible light (400–700 nm) but that its spectral properties may be different from that of opsins1 and 2. Changes in the relative levels of these opsins may underlie some of the dramatic day–night changes in Limulus photoreceptor function and may produce a diurnal change in their spectral sensitivity.

Technology to generate single cell RNA-sequencing (scRNA-seq) datasets and tools to annotate them have advanced rapidly in the past several years. Such tools generally rely on existing transcriptomic datasets or curated databases of cell type defining genes, while the application of scalable natural language processing (NLP) methods to enhance analysis workflows has not been adequately explored. Here we deployed an NLP framework to objectively quantify associations between a comprehensive set of over 20,000 human protein-coding genes and over 500 cell type terms across over 26 million biomedical documents. The resultant gene-cell type associations (GCAs) are significantly stronger between a curated set of matched cell type-marker pairs than the complementary set of mismatched pairs (Mann Whitney p = 6.15 × 10⁻⁷⁶, r = 0.24; cohen’s D = 2.6). Building on this, we developed an augmented annotation algorithm (single cell Annotation via Literature Encoding, or scALE) that leverages GCAs to categorize cell clusters identified in scRNA-seq datasets, and we tested its ability to predict the cellular identity of 133 clusters from nine datasets of human breast, colon, heart, joint, ovary, prostate, skin, and small intestine tissues. With the optimized settings, the true cellular identity matched the top prediction in 59% of tested clusters and was present among the top five predictions for 91% of clusters. scALE slightly outperformed an existing method for reference data driven automated cluster annotation, and we demonstrate that integration of scALE can meaningfully improve the annotations derived from such methods. Further, contextualization of differential expression analyses with these GCAs highlights poorly characterized markers of well-studied cell types, such as CLIC6 and DNASE1L3 in retinal pigment epithelial cells and endothelial cells, respectively. Taken together, this study illustrates for the first time how the systematic application of a literature-derived knowledge graph can expedite and enhance the annotation and interpretation of scRNA-seq data.