Single-cell and spatial transcriptomics analysis of non-small cell lung cancer.

Lung cancer is the second most frequently diagnosed cancer and the leading cause of cancer-related mortality worldwide. Tumour ecosystems feature diverse immune cell types. Myeloid cells, in particular, are prevalent and have a well-established role in promoting the disease. In our study, we profile approximately 900,000 cells from 25 treatment-naive patients with adenocarcinoma and squamous-cell carcinoma by single-cell and spatial transcriptomics. We note an inverse relationship between anti-inflammatory macrophages and NK cells/T cells, and with reduced NK cell cytotoxicity within the tumour. While we observe a similar cell type composition in both adenocarcinoma and squamous-cell carcinoma, we detect significant differences in the co-expression of various immune checkpoint inhibitors. Moreover, we reveal evidence of a transcriptional "reprogramming" of macrophages in tumours, shifting them towards cholesterol export and adopting a foetal-like transcriptional signature which promotes iron efflux. Our multi-omic resource offers a high-resolution molecular map of tumour-associated macrophages, enhancing our understanding of their role within the tumour microenvironment.

Tracing the First Days of the SARS-CoV-2 Pandemic in Greece and the Role of the First Imported Group of Travelers.

The first SARS-CoV-2 case in Greece was confirmed on February 26, 2020, and since then, multiple strains have circulated the country, leading to regional and country-wide outbreaks. Our aim is to enlighten the events that took place during the first days of the SARS-CoV-2 pandemic in Greece, focusing on the role of the first imported group of travelers. We used whole-genome SARS-CoV-2 sequences obtained from the infected travelers of the group as well as Greece-derived and globally subsampled sequences and applied dedicated phylogenetics and phylodynamics tools as well as in-house-developed bioinformatics pipelines. Our analyses reveal the genetic variants circulating in Greece during the first days of the pandemic and the role of the group's imported strains in the course of the first pandemic wave in Greece. The strain that dominated in Greece throughout the first wave, bearing the D614G mutation, was primarily imported from a certain group of travelers, while molecular and clinical data suggest that the infection of the travelers occurred in Egypt. Founder effects early in the pandemic are important for the success of certain strains, as those arriving early, several times, and to diverse locations lead to the formation of large transmission clusters that can be estimated using molecular epidemiology approaches and can be a useful surveillance tool for the prioritization of nonpharmaceutical interventions and combating present and future outbreaks. IMPORTANCE The strain that dominated in Greece during the first pandemic wave was primarily imported from a group of returning travelers in February 2020, while molecular and clinical data suggest that the origin of the transmission was Egypt. The observed molecular transmission clusters reflect the transmission dynamics of this particular strain bearing the D614G mutation while highlighting the necessity of their use as a surveillance tool for the prioritization of nonpharmaceutical interventions and combating present and future outbreaks.

Radiomics/Radiogenomics in Lung Cancer: Basic Principles and Initial Clinical Results.

Lung cancer is the leading cause of cancer-related deaths worldwide, and elucidation of its complicated pathobiology has been traditionally targeted by studies incorporating genomic as well other high-throughput approaches. Recently, a collection of methods used for cancer imaging, supplemented by quantitative aspects leading towards imaging biomarker assessment termed "radiomics", has introduced a novel dimension in cancer research. Integration of genomics and radiomics approaches, where identifying the biological basis of imaging phenotypes is feasible due to the establishment of associations between molecular features at the genomic-transcriptomic-proteomic level and radiological features, has recently emerged termed radiogenomics. This review article aims to briefly describe the main aspects of radiogenomics, while discussing its basic limitations related to lung cancer clinical applications for clinicians, researchers and patients.

Cell adhesion tunes inflammatory TPL2 kinase signal transduction.

Signaling through adhesion-related molecules is important for cancer growth and metastasis and cancer cells are resistant to anoikis, a form of cell death ensued by cell detachment from the extracellular matrix. Herein, we report that detached carcinoma cells and immortalized fibroblasts display defects in TNF and CD40 ligand (CD40L)-induced MEK-ERK signaling. Cell detachment results in reduced basal levels of the MEK kinase TPL2, compromises TPL2 activation and sensitizes carcinoma cells to death-inducing receptor ligands, mimicking the synthetic lethal interactions between TPL2 inactivation and TNF or CD40L stimulation. Focal Adhesion Kinase (FAK), which is activated in focal adhesions and mediates anchorage-dependent survival signaling, was found to sustain steady state TPL2 protein levels and to be required for TNF-induced TPL2 signal transduction. We show that when FAK levels are reduced, as seen in certain types of malignancy or malignant cell populations, the formation of cIAP2:RIPK1 complexes increases, leading to reduced TPL2 expression levels by a dual mechanism: first, by the reduction in the levels of NF-κΒ1 which is required for TPL2 stability; second, by the engagement of an RelA NF-κΒ pathway that elevates interleukin-6 production, leading to activation of STAT3 and its transcriptional target SKP2 which functions as a TPL2 E3 ubiquitin ligase. These data underscore a new mode of regulation of TNF family signal transduction on the TPL2-MEK-ERK branch by adhesion-related molecules that may have important ramifications for cancer therapy.

Generation of Non-Small Cell Lung Cancer Patient-Derived Xenografts to Study Intratumor Heterogeneity.

Recent advances in sequencing technologies have allowed the in-depth molecular study of tumors, even at the single cell level. Sequencing efforts have uncovered a previously unappreciated heterogeneity among tumor cells, which has been postulated to be the driving force of tumor evolution and to facilitate recurrence, metastasis, and drug resistance. In the current study, focused on early-stage operable non-small cell lung cancer, we used tumor growth in patient-derived xenograft (PDX) models in mice as a fast-forward tumor evolution process to investigate the molecular characteristics of tumor cells that grow in mice, as well as the parameters that affect the grafting efficiency. We found that squamous cell carcinomas grafted significantly more efficiently compared with adenocarcinomas. Advanced stage, patient age and primary tumor size were positively correlated with grafting. Additionally, we isolated and characterized circulating tumor cells (CTC) from patients' peripheral blood and found that the presence of CTCs expressing epithelial-to-mesenchymal (EMT) markers correlated with the grafting potential. Interestingly, exome sequencing of the PDX tumor identified genetic alterations in DNA repair and genome integrity genes that were under-represented in the human primary counterpart. In conclusion, through the generation of a PDX biobank of NSCLC, we identified the clinical and molecular properties of tumors that affected growth in mice.

Single-cell transcriptional analysis reveals ILC-like cells in zebrafish.

Innate lymphoid cells (ILCs) are important mediators of the immune response and homeostasis in barrier tissues of mammals. However, the existence and function of ILCs in other vertebrates are poorly understood. Here, we use single-cell RNA sequencing to generate a comprehensive atlas of zebrafish lymphocytes during tissue homeostasis and after immune challenge. We profiled 14,080 individual cells from the gut of wild-type zebrafish, as well as of rag1-deficient zebrafish that lack T and B cells, and discovered populations of ILC-like cells. We uncovered a rorc-positive subset of ILCs that could express cytokines associated with type 1, 2, and 3 responses upon immune challenge. Specifically, these ILC-like cells expressed il22 and tnfa after exposure to inactivated bacteria or il13 after exposure to helminth extract. Cytokine-producing ILC-like cells express a specific repertoire of novel immune-type receptors, likely involved in recognition of environmental cues. We identified additional novel markers of zebrafish ILCs and generated a cloud repository for their in-depth exploration.

Single-cell RNA-sequencing uncovers transcriptional states and fate decisions in haematopoiesis.

The success of marker-based approaches for dissecting haematopoiesis in mouse and human is reliant on the presence of well-defined cell surface markers specific for diverse progenitor populations. An inherent problem with this approach is that the presence of specific cell surface markers does not directly reflect the transcriptional state of a cell. Here, we used a marker-free approach to computationally reconstruct the blood lineage tree in zebrafish and order cells along their differentiation trajectory, based on their global transcriptional differences. Within the population of transcriptionally similar stem and progenitor cells, our analysis reveals considerable cell-to-cell differences in their probability to transition to another committed state. Once fate decision is executed, the suppression of transcription of ribosomal genes and upregulation of lineage-specific factors coordinately controls lineage differentiation. Evolutionary analysis further demonstrates that this haematopoietic programme is highly conserved between zebrafish and higher vertebrates.

CD4-Transgenic Zebrafish Reveal Tissue-Resident Th2- and Regulatory T Cell-like Populations and Diverse Mononuclear Phagocytes.

CD4+ T cells are at the nexus of the innate and adaptive arms of the immune system. However, little is known about the evolutionary history of CD4+ T cells, and it is unclear whether their differentiation into specialized subsets is conserved in early vertebrates. In this study, we have created transgenic zebrafish with vibrantly labeled CD4+ cells allowing us to scrutinize the development and specialization of teleost CD4+ leukocytes in vivo. We provide further evidence that CD4+ macrophages have an ancient origin and had already emerged in bony fish. We demonstrate the utility of this zebrafish resource for interrogating the complex behavior of immune cells at cellular resolution by the imaging of intimate contacts between teleost CD4+ T cells and mononuclear phagocytes. Most importantly, we reveal the conserved subspecialization of teleost CD4+ T cells in vivo. We demonstrate that the ancient and specialized tissues of the gills contain a resident population of il-4/13b-expressing Th2-like cells, which do not coexpress il-4/13a Additionally, we identify a contrasting population of regulatory T cell-like cells resident in the zebrafish gut mucosa, in marked similarity to that found in the intestine of mammals. Finally, we show that, as in mammals, zebrafish CD4+ T cells will infiltrate melanoma tumors and obtain a phenotype consistent with a type 2 immune microenvironment. We anticipate that this unique resource will prove invaluable for future investigation of T cell function in biomedical research, the development of vaccination and health management in aquaculture, and for further research into the evolution of adaptive immunity.

Discovering gene re-ranking efficiency and conserved gene-gene relationships derived from gene co-expression network analysis on breast cancer data.

Systemic approaches are essential in the discovery of disease-specific genes, offering a different perspective and new tools on the analysis of several types of molecular relationships, such as gene co-expression or protein-protein interactions. However, due to lack of experimental information, this analysis is not fully applicable. The aim of this study is to reveal the multi-potent contribution of statistical network inference methods in highlighting significant genes and interactions. We have investigated the ability of statistical co-expression networks to highlight and prioritize genes for breast cancer subtypes and stages in terms of: (i) classification efficiency, (ii) gene network pattern conservation, (iii) indication of involved molecular mechanisms and (iv) systems level momentum to drug repurposing pipelines. We have found that statistical network inference methods are advantageous in gene prioritization, are capable to contribute to meaningful network signature discovery, give insights regarding the disease-related mechanisms and boost drug discovery pipelines from a systems point of view.

Bioinformatics methods in drug repurposing for Alzheimer's disease.

Alarming epidemiological features of Alzheimer's disease impose curative treatment rather than symptomatic relief. Drug repurposing, that is reappraisal of a substance's indications against other diseases, offers time, cost and efficiency benefits in drug development, especially when in silico techniques are used. In this study, we have used gene signatures, where up- and down-regulated gene lists summarize a cell's gene expression perturbation from a drug or disease. To cope with the inherent biological and computational noise, we used an integrative approach on five disease-related microarray data sets of hippocampal origin with three different methods of evaluating differential gene expression and four drug repurposing tools. We found a list of 27 potential anti-Alzheimer agents that were additionally processed with regard to molecular similarity, pathway/ontology enrichment and network analysis. Protein kinase C, histone deacetylase, glycogen synthase kinase 3 and arginase inhibitors appear consistently in the resultant drug list and may exert their pharmacologic action in an epidermal growth factor receptor-mediated subpathway of Alzheimer's disease.

ZoomOut: Analyzing Multiple Networks as Single Nodes.

We have developed ZoomOut web server in order to provide the research community with a tool for analysis, visualization and clustering of networks as a super network, based on their calculated feature properties. Networks can be analysed and be further treated as single nodes in a super network that describe their relations. Specifically, the user interface is divided into three main sections: the Workspace, the Networks Feature Calculations and the Clustering Networks section. In the Workspace section, users are able to upload and manage multiple networks for further processing. In the Networks Feature Calculations section, a variety of network properties are calculated as features for each uploaded network. In the Clustering Networks section, users are able to apply clustering by selecting from the list of previously calculated features. All processed networks can also be visualized as a super interactive network, were interconnections among networks are based on the calculated clustering distances. To the best of our knowledge, this is the first available web-service that allows users to manage, quantify and visualize multiple networks at the same time, handling them as parts of a larger network with properties calculated in an upper scale. The ZoomOut web-application is available at http://bioserver-3.bioacademy.gr/Bioserver/ZoomOut.

GeneStoryTeller: a mobile app for quick and comprehensive information retrieval of human genes.

In the last few years, mobile devices such as smartphones and tablets have become an integral part of everyday life, due to their software/hardware rapid development, as well as the increased portability they offer. Nevertheless, up to now, only few Apps have been developed in the field of bioinformatics, capable to perform fast and robust access to services. We have developed the GeneStoryTeller, a mobile application for Android platforms, where users are able to instantly retrieve information regarding any recorded human gene, derived from eight publicly available databases, as a summary story. Complementary information regarding gene-drugs interactions, functional annotation and disease associations for each selected gene is also provided in the gene story. The most challenging part during the development of the GeneStoryTeller was to keep balance between storing data locally within the app and obtaining the updated content dynamically via a network connection. This was accomplished with the implementation of an administrative site where data are curated and synchronized with the application requiring a minimum human intervention.

A Web-based database of genetic association studies in cutaneous melanoma enhanced with network-driven data exploration tools.

The publicly available online database MelGene provides a comprehensive, regularly updated, collection of data from genetic association studies in cutaneous melanoma (CM), including random-effects meta-analysis results of all eligible polymorphisms. The updated database version includes data from 192 publications with information on 1114 significantly associated polymorphisms across 280 genes, along with new front-end and back-end capabilities. Various types of relationships between data are calculated and visualized as networks. We constructed 13 different networks containing the polymorphisms and the genes included in MelGene. We explored the derived network representations under the following questions: (i) are there nodes that deserve consideration regarding their network connectivity characteristics? (ii) What is the relation of either the genome-wide or nominally significant CM polymorphisms/genes with the ones highlighted by the network representation? We show that our network approach using the MelGene data reveals connections between statistically significant genes/ polymorphisms and other genes/polymorphisms acting as 'hubs' in the reconstructed networks. To the best of our knowledge, this is the first database containing data from a comprehensive field synopsis and systematic meta-analyses of genetic polymorphisms in CM that provides user-friendly tools for in-depth molecular network visualization and exploration. The proposed network connections highlight potentially new loci requiring further investigation of their relation to melanoma risk. Database URL: http://www.melgene.org.

Multimodality GPU-based computer-assisted diagnosis of breast cancer using ultrasound and digital mammography images.

PURPOSE: To improve the computer-aided diagnosis of breast lesions, by designing a pattern recognition system (PR-system) on commercial graphics processing unit (GPU) cards using parallel programming and textural information from multimodality imaging. MATERIAL AND METHODS: Patients with histologically verified breast lesions underwent both ultrasound (US) and digital mammography (DM), lesions were outlined on the images by an experienced radiologist, and textural features were calculated. The PR-system was designed to provide highest possible precision by programming in parallel the multiprocessors of the NVIDIA's GPU cards, GeForce 8800GT or 580GTX, and using the CUDA programming framework and C++. The PR-system was built around the probabilistic neural network classifier, and its performance was evaluated by a re-substitution method, for estimating the system's highest accuracy, and by the external cross-validation method, for assessing the PR-system's unbiased accuracy to new, "unseen" by the system, data. RESULTS: Classification accuracies for discriminating malignant from benign lesions were as follows: 85.5 % using US-features alone, 82.3 % employing DM features alone, and 93.5 % combining US and DM features. Mean accuracy to new "unseen" data for the combined US and DM features was 81 %. Those classification accuracies were about 10 % higher than accuracies achieved on a single CPU, using sequential programming methods, and 150-fold faster. CONCLUSION: The proposed PR-system improves breast-lesion discrimination accuracy, it may be redesigned on site when new verified data are incorporated in its depository, and it may serve as a second opinion tool in a clinical environment.

Segmentation of complementary DNA microarray images by wavelet-based Markov random field model.

A wavelet-based modification of the Markov random field (WMRF) model is proposed for segmenting complementary DNA (cDNA) microarray images. For evaluation purposes, five simulated and a set of five real microarray images were used. The one-level stationary wavelet transform (SWT) of each microarray image was used to form two images, a denoised image, using hard thresholding filter, and a magnitude image, from the amplitudes of the horizontal and vertical components of SWT. Elements from these two images were suitably combined to form the WMRF model for segmenting spots from their background. The WMRF was compared against the conventional MRF and the Fuzzy C means (FCM) algorithms on simulated and real microarray images and their performances were evaluated by means of the segmentation matching factor (SMF) and the coefficient of determination (r2). Additionally, the WMRF was compared against the SPOT and SCANALYZE, and performances were evaluated by the mean absolute error (MAE) and the coefficient of variation (CV). The WMRF performed more accurately than the MRF and FCM (SMF: 92.66, 92.15, and 89.22, r2 : 0.92, 0.90, and 0.84, respectively) and achieved higher reproducibility than the MRF, SPOT, and SCANALYZE (MAE: 497, 1215, 1180, and 503, CV: 0.88, 1.15, 0.93, and 0.90, respectively).

Complementary DNA microarray image processing based on the fuzzy gaussian mixture model.

The objective of this paper was to investigate the segmentation ability of the fuzzy Gaussian mixture model (FGMM) clustering algorithm, applied on complementary DNA (cDNA) images. Following a standard established procedure, a simulated microarray image of 1600 cells, each containing one spot, was produced. For further evaluation of the algorithm, three real microarray images were also used, each containing 6400 spots. For the task of locating spot borders and surrounding background (BG) in each cell, an automatic gridding process was developed and applied on microarray images. The FGMM and the Gaussian mixture model (GMM) algorithms were applied to each cell with the purpose of discriminating foreground (FG) from BG. The segmentation abilities of both algorithms were evaluated by means of the segmentation matching factor, coefficient of determination, and concordance correlation, in respect to the actual classes (FG-BG pixels) of the simulated spots. Pairwise correlation and mean absolute error of the real images among replicates were also calculated. The FGMM was found to perform better and with equal processing time, as compared to the GMM, rendering the FGMM algorithm an efficient alternative for segmenting cDNA microarray images.