Peripheral blood immune parameters, response, and adverse events after neoadjuvant chemotherapy plus durvalumab in early-stage triple-negative breast cancer.

PURPOSE: We evaluated T- and B-cell receptor (TCR and BCR) repertoire diversity and 38 serum cytokines in pre- and post-treatment peripheral blood of 66 patients with triple-negative breast cancer (TNBC) who received neoadjuvant chemotherapy plus durvalumab and assessed associations with pathologic response and immune-related adverse events (irAEs) during treatment. METHODS: Genomic DNA was isolated from buffy coat for TCR and BCR clonotype profiling using the Immunoseq platform and diversity was quantified with Pielou's evenness index. MILLIPLEX MAP Human Cytokine/Chemokine Magnetic Bead Panel was used to measure serum cytokine levels, which were compared between groups using moderated t-statistic with Benjamini-Hochberg correction for multiple testing. RESULTS: TCR and BCR diversity was high (Pielou's index > 0.75) in all samples. Baseline receptor diversities and change in diversity pre- and post-treatment were not associated with pathologic response or irAE status, except for BCR diversity that was significantly lower post-treatment in patients who developed irAE (unadjusted p = 0.0321). Five cytokines increased after treatment in patients with pathologic complete response (pCR) but decreased in patients with RD, most prominently IL-8. IFNγ, IL-7, and GM-CSF levels were higher in pre-treatment than in post-treatment samples of patients who developed irAEs but were lower in those without irAEs. CONCLUSION: Baseline peripheral blood cytokine levels may predict irAEs in patients treated with immune checkpoint inhibitors and chemotherapy, and increased post-treatment B-cell clonal expansion might mediate irAEs.

Pre-analytical effects on whole transcriptome and targeted RNA sequencing analysis in cytology: The effects of prolonged time in storage of effusion specimens prior to preservation.

OBJECTIVES: To investigate the pre-analytics of the molecular testing of cytology specimens, we studied the effects of time in refrigerator storage (4°C) of malignant effusions on RNA sequencing (RNAseq) results. METHODS: Ten effusion specimens were stored in a refrigerator (4°C) for different durations (day 0, 1, 4, and 7). All specimens were prepared as cytospins fixed in either Carnoy's solution or 95% ethanol (EtOH) and in an RNA preservative for a fresh frozen (FF) high-quality reference. Whole transcriptome (wt) and targeted (t)RNAseq of two multigene expression signatures were performed. We then compared transcript expression levels (including mutant allele fraction) according to pre-analytical variables using a concordance correlation coefficient (CCC) and a mixed effect model. RESULTS: Sequencing results were mostly stable over increasing time in storage. Cytospins fixed in Carnoy's solution were more concordant with FF samples than cytospins fixed in 95% EtOH at all timepoints. This finding was consistent for both wtRNAseq (averages: day 0 CCC = 0.98 vs 0.91; day 7 CCC = 0.88 vs 0.78) and tRNAseq methods (averages: day 0 CCC = 0.98 vs 0.81; day 7 CCC = 0.98 vs 0.90). Cytospins fixed in Carnoy's solution did not show significant changes in expression over timepoints or between expression signatures, whereas 95% EtOH did. CONCLUSION: RNAseq can be accurately performed on effusion specimens after prolonged refrigerator storage. RNA extracted from scraped cytospin slides fixed in Carnoy's solution was marginally superior to 95% EtOH fixation, but either method had comparable analytic performance to high-quality FF RNA samples.

Targeted RNAseq assay incorporating unique molecular identifiers for improved quantification of gene expression signatures and transcribed mutation fraction in fixed tumor samples.

BACKGROUND: Our objective was to assess whether modifications to a customized targeted RNA sequencing (RNAseq) assay to include unique molecular identifiers (UMIs) that collapse read counts to their source mRNA counts would improve quantification of transcripts from formalin-fixed paraffin-embedded (FFPE) tumor tissue samples. The assay (SET4) includes signatures that measure hormone receptor and PI3-kinase related transcriptional activity (SETER/PR and PI3Kges), and measures expression of selected activating point mutations and key breast cancer genes. METHODS: Modifications included steps to introduce eight nucleotides-long UMIs during reverse transcription (RT) in bulk solution, followed by polymerase chain reaction (PCR) of labeled cDNA in droplets, with optimization of the polymerase enzyme and reaction conditions. We used Lin's concordance correlation coefficient (CCC) to measure concordance, including precision (Rho) and accuracy (Bias), and nonparametric tests (Wilcoxon, Levene's) to compare the modified (NEW) SET4 assay to the original (OLD) SET4 assay and to whole transcriptome RNAseq using RNA from matched fresh frozen (FF) and FFPE samples from 12 primary breast cancers. RESULTS: The modified (NEW) SET4 assay measured single transcripts (p< 0.001) and SETER/PR (p=0.002) more reproducibly in technical replicates from FFPE samples. The modified SET4 assay was more precise for measuring single transcripts (Rho 0.966 vs 0.888, p< 0.01) but not multigene expression signatures SETER/PR (Rho 0.985 vs 0.968) or PI3Kges (Rho 0.985 vs 0.946) in FFPE, compared to FF samples. It was also more precise than wtRNAseq of FFPE for measuring transcripts (Rho 0.986 vs 0.934, p< 0.001) and SETER/PR (Rho 0.993 vs 0.915, p=0.004), but not PI3Kges (Rho 0.988 vs 0.945, p=0.051). Accuracy (Bias) was comparable between protocols. Two samples carried a PIK3CA mutation, and measurements of transcribed mutant allele fraction was similar in FF and FFPE samples and appeared more precise with the modified SET4 assay. Amplification efficiency (reads per UMI) was consistent in FF and FFPE samples, and close to the theoretically expected value, when the library size exceeded 400,000 aligned reads. CONCLUSIONS: Modifications to the targeted RNAseq protocol for SET4 assay significantly increased the precision of UMI-based and reads-based measurements of individual transcripts, multi-gene signatures, and mutant transcript fraction, particularly with FFPE samples.

Analysis of the Applicability of microRNAs in Peripheral Blood Leukocytes as Biomarkers of Sensitivity and Exposure to Fractionated Radiotherapy towards Breast Cancer.

Biomarkers for predicting individual response to radiation and for dose verification are needed to improve radiotherapy. A biomarker should optimally show signal fidelity, meaning that its level is stable and proportional to the absorbed dose. miRNA levels in human blood serum were suggested as promising biomarkers. The aim of the present investigation was to test the miRNA biomarker in leukocytes of breast cancer patients undergoing external beam radiotherapy. Leukocytes were isolated from blood samples collected prior to exposure (control); on the day when a total dose of 2 Gy, 10 Gy, or 20 Gy was reached; and one month after therapy ended (46-50 Gy in total). RNA sequencing was performed and univariate analysis was used to analyse the effect of the radiation dose on the expression of single miRNAs. To check if combinations of miRNAs can predict absorbed dose, a multinomial logistic regression model was built using a training set from eight patients (representing 40 samples) and a validation set with samples from the remaining eight patients (15 samples). Finally, Broadside, an explorative interaction mining tool, was used to extract sets of interacting miRNAs. The most prominently increased miRNA was miR-744-5p, followed by miR-4461, miR-34a-5p, miR-6513-5p, miR-1246, and miR-454-3p. Decreased miRNAs were miR-3065-3p, miR-103a-2-5p, miR-30b-3p, and miR-5690. Generally, most miRNAs showed a relatively strong inter-individual variability and different temporal patterns over the course of radiotherapy. In conclusion, miR-744-5p shows promise as a stable miRNA marker, but most tested miRNAs displayed individual signal variability which, at least in this setting, may exclude them as sensitive biomarkers of radiation response.

BatchI: Batch effect Identification in high-throughput screening data using a dynamic programming algorithm.

MOTIVATION: In contemporary biological experiments, bias, which interferes with the measurements, requires attentive processing. Important sources of bias in high-throughput biological experiments are batch effects and diverse methods towards removal of batch effects have been established. These include various normalization techniques, yet many require knowledge on the number of batches and assignment of samples to batches. Only few can deal with the problem of identification of batch effect of unknown structure. For this reason, an original batch identification algorithm through dynamical programming is introduced for omics data that may be sorted on a timescale. RESULTS: BatchI algorithm is based on partitioning a series of high-throughput experiment samples into sub-series corresponding to estimated batches. The dynamic programming method is used for splitting data with maximal dispersion between batches, while maintaining minimal within batch dispersion. The procedure has been tested on a number of available datasets with and without prior information about batch partitioning. Datasets with a priori identified batches have been split accordingly, measured with weighted average Dice Index. Batch effect correction is justified by higher intra-group correlation. In the blank datasets, identified batch divisions lead to improvement of parameters and quality of biological information, shown by literature study and Information Content. The outcome of the algorithm serves as a starting point for correction methods. It has been demonstrated that omitting the essential step of batch effect control may lead to waste of valuable potential discoveries. AVAILABILITY AND IMPLEMENTATION: The implementation is available within the BatchI R package at http://zaed.aei.polsl.pl/index.php/pl/111-software. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Gene set enrichment for reproducible science: comparison of CERNO and eight other algorithms.

MOTIVATION: Analysis of gene set (GS) enrichment is an essential part of functional omics studies. Here, we complement the established evaluation metrics of GS enrichment algorithms with a novel approach to assess the practical reproducibility of scientific results obtained from GS enrichment tests when applied to related data from different studies. RESULTS: We evaluated eight established and one novel algorithm for reproducibility, sensitivity, prioritization, false positive rate and computational time. In addition to eight established algorithms, we also included Coincident Extreme Ranks in Numerical Observations (CERNO), a flexible and fast algorithm based on modified Fisher P-value integration. Using real-world datasets, we demonstrate that CERNO is robust to ranking metrics, as well as sample and GS size. CERNO had the highest reproducibility while remaining sensitive, specific and fast. In the overall ranking Pathway Analysis with Down-weighting of Overlapping Genes, CERNO and over-representation analysis performed best, while CERNO and GeneSetTest scored high in terms of reproducibility. AVAILABILITY AND IMPLEMENTATION: tmod package implementing the CERNO algorithm is available from CRAN (cran.r-project.org/web/packages/tmod/index.html) and an online implementation can be found at http://tmod.online/. The datasets analyzed in this study are widely available in the KEGGdzPathwaysGEO, KEGGandMetacoreDzPathwaysGEO R package and GEO repository. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Technical Validity of a Customized Assay of Sensitivity to Endocrine Therapy Using Sections from Fixed Breast Cancer Tissue.

BACKGROUND: We translated a multigene expression index to predict sensitivity to endocrine therapy for Stage II-III breast cancer (SET2,3) to hybridization-based expression assays of formalin-fixed paraffin-embedded (FFPE) tissue sections. Here we report the technical validity with FFPE samples, including preanalytical and analytical performance. METHODS: We calibrated SET2,3 from microarrays (Affymetrix U133A) of frozen samples to hybridization-based assays of FFPE tissue, using bead-based QuantiGene Plex (QGP) and slide-based NanoString (NS). The following preanalytical and analytical conditions were tested in controlled studies: replicates within and between frozen and fixed samples, age of paraffin blocks, homogenization of fixed sections versus extracted RNA, core biopsy versus surgically resected tumor, technical replicates, precision over 20 weeks, limiting dilution, linear range, and analytical sensitivity. Lin's concordance correlation coefficient (CCC) was used to measure concordance between measurements. RESULTS: SET2,3 index was calibrated to use with QGP (CCC 0.94) and NS (CCC 0.93) technical platforms, and was validated in two cohorts of older fixed samples using QGP (CCC 0.72, 0.85) and NS (CCC 0.78, 0.78). QGP assay was concordant using direct homogenization of fixed sections versus purified RNA (CCC 0.97) and between core and surgical sample types (CCC 0.90), with 100% accuracy in technical replicates, 1-9% coefficient of variation over 20 weekly tests, linear range 3.0-11.5 (log2 counts), and analytical sensitivity ≥2.0 (log2 counts). CONCLUSIONS: Measurement of the novel SET2,3 assay was technically valid from fixed tumor sections of biopsy or resection samples using simple, inexpensive, hybridization methods, without the need for RNA purification.

The impact of RNA extraction method on accurate RNA sequencing from formalin-fixed paraffin-embedded tissues.

BACKGROUND: Utilization of RNA sequencing methods to measure gene expression from archival formalin-fixed paraffin-embedded (FFPE) tumor samples in translational research and clinical trials requires reliable interpretation of the impact of pre-analytical variables on the data obtained, particularly the methods used to preserve samples and to purify RNA. METHODS: Matched tissue samples from 12 breast cancers were fresh frozen (FF) and preserved in RNAlater or fixed in formalin and processed as FFPE tissue. Total RNA was extracted and purified from FF samples using the Qiagen RNeasy kit, and in duplicate from FFPE tissue sections using three different kits (Norgen, Qiagen and Roche). All RNA samples underwent whole transcriptome RNA sequencing (wtRNAseq) and targeted RNA sequencing for 31 transcripts included in a signature of sensitivity to endocrine therapy. We assessed the effect of RNA extraction kit on the reliability of gene expression levels using linear mixed-effects model analysis, concordance correlation coefficient (CCC) and differential analysis. All protein-coding genes in the wtRNAseq and three gene expression signatures for breast cancer were assessed for concordance. RESULTS: Despite variable quality of the RNA extracted from FFPE samples by different kits, all had similar concordance of overall gene expression from wtRNAseq between matched FF and FFPE samples (median CCC 0.63-0.66) and between technical replicates (median expression difference 0.13-0.22). More than half of genes were differentially expressed between FF and FFPE, but with low fold change (median |LFC| 0.31-0.34). Two out of three breast cancer signatures studied were highly robust in all samples using any kit, whereas the third signature was similarly discordant irrespective of the kit used. The targeted RNAseq assay was concordant between FFPE and FF samples using any of the kits (CCC 0.91-0.96). CONCLUSIONS: The selection of kit to purify RNA from FFPE did not influence the overall quality of results from wtRNAseq, thus variable reproducibility of gene signatures probably relates to the reliability of individual gene selected and possibly to the algorithm. Targeted RNAseq showed promising performance for clinical deployment of quantitative assays in breast cancer from FFPE samples, although numerical scores were not identical to those from wtRNAseq and would require calibration.

Lung cancer survival and comorbidities in lung cancer screening participants of the Gdansk screening cohort.

BACKGROUND: In 2010, the World Health Organisation recommended implementation of screening programmes in four groups of diseases-neoplasms, cardiovascular diseases (CVD), diabetes mellitus (DM) and chronic obstructive pulmonary disease (COPD). It is due to the fact that they share the same, modifiable risk factors. METHODS: Between 2009 and 2011, 8637 heavy smokers (aged 50-75, smoking history >20 pack-years) were screened in the Pomeranian Pilot Lung Cancer Screening Programme (PPP) in Gdansk, Poland. We looked at 5-year follow-up and analysed the medical events and comorbidities of all participants. One health care provider in the Polish health care system provides a unique opportunity to gather most reliable data on all medical events in each person. RESULTS: In 52.0% of lung cancer screening participants CVD (33.5%), DM (26.0%) and COPD (21.0%) were diagnosed. Prevalence of these diseases is higher in lung cancer patients than in the non-cancer screening group (P < 0.0001). One hundred and seven (1.2%) lung cancers were diagnosed during PPP programme performance and another 382 cases (4.4%) in the 5-year follow-up, so the potential mean annual lung cancer detection rate is 0.77%. CONCLUSIONS: Lung cancer screening programme offers a great potential for joint screening of lung cancer, CVD, diabetes and COPD.

Ranking metrics in gene set enrichment analysis: do they matter?

BACKGROUND: There exist many methods for describing the complex relation between changes of gene expression in molecular pathways or gene ontologies under different experimental conditions. Among them, Gene Set Enrichment Analysis seems to be one of the most commonly used (over 10,000 citations). An important parameter, which could affect the final result, is the choice of a metric for the ranking of genes. Applying a default ranking metric may lead to poor results. METHODS AND RESULTS: In this work 28 benchmark data sets were used to evaluate the sensitivity and false positive rate of gene set analysis for 16 different ranking metrics including new proposals. Furthermore, the robustness of the chosen methods to sample size was tested. Using k-means clustering algorithm a group of four metrics with the highest performance in terms of overall sensitivity, overall false positive rate and computational load was established i.e. absolute value of Moderated Welch Test statistic, Minimum Significant Difference, absolute value of Signal-To-Noise ratio and Baumgartner-Weiss-Schindler test statistic. In case of false positive rate estimation, all selected ranking metrics were robust with respect to sample size. In case of sensitivity, the absolute value of Moderated Welch Test statistic and absolute value of Signal-To-Noise ratio gave stable results, while Baumgartner-Weiss-Schindler and Minimum Significant Difference showed better results for larger sample size. Finally, the Gene Set Enrichment Analysis method with all tested ranking metrics was parallelised and implemented in MATLAB, and is available at https://github.com/ZAEDPolSl/MrGSEA . CONCLUSIONS: Choosing a ranking metric in Gene Set Enrichment Analysis has critical impact on results of pathway enrichment analysis. The absolute value of Moderated Welch Test has the best overall sensitivity and Minimum Significant Difference has the best overall specificity of gene set analysis. When the number of non-normally distributed genes is high, using Baumgartner-Weiss-Schindler test statistic gives better outcomes. Also, it finds more enriched pathways than other tested metrics, which may induce new biological discoveries.

Radiation-induced changes in serum lipidome of head and neck cancer patients.

Cancer radiotherapy (RT) induces response of the whole patient's body that could be detected at the blood level. We aimed to identify changes induced in serum lipidome during RT and characterize their association with doses and volumes of irradiated tissue. Sixty-six patients treated with conformal RT because of head and neck cancer were enrolled in the study. Blood samples were collected before, during and about one month after the end of RT. Lipid extracts were analyzed using MALDI-oa-ToF mass spectrometry in positive ionization mode. The major changes were observed when pre-treatment and within-treatment samples were compared. Levels of several identified phosphatidylcholines, including (PC34), (PC36) and (PC38) variants, and lysophosphatidylcholines, including (LPC16) and (LPC18) variants, were first significantly decreased and then increased in post-treatment samples. Intensities of changes were correlated with doses of radiation received by patients. Of note, such correlations were more frequent when low-to-medium doses of radiation delivered during conformal RT to large volumes of normal tissues were analyzed. Additionally, some radiation-induced changes in serum lipidome were associated with toxicity of the treatment. Obtained results indicated the involvement of choline-related signaling and potential biological importance of exposure to clinically low/medium doses of radiation in patient's body response to radiation.

miR-122 and miR-21 are Stable Components of miRNA Signatures of Early Lung Cancer after Validation in Three Independent Cohorts.

Several panels of circulating miRNAs have been reported as potential biomarkers of early lung cancer, yet the overlap of components between different panels is limited, and the universality of proposed biomarkers has been minimal across proposed panels. To assess the stability of the diagnostic potential of plasma miRNA signature of early lung cancer among different cohorts, a panel of 24 miRNAs tested in the frame of one lung cancer screening study (MOLTEST-2013, Poland) was validated with material collected in the frame of two other screening studies (MOLTEST-BIS, Poland; and SMAC, Italy) using the same standardized analytical platform (the miRCURY LNA miRNA PCR assay). On analysis of selected miRNAs, two associated with lung cancer development, miR-122 and miR-21, repetitively differentiated healthy participants from individuals with lung cancer. Additionally, miR-144 differentiated controls from cases specifically in subcohorts with adenocarcinoma. Other tested miRNAs did not overlap in the three cohorts. Classification models based on neither a single miRNA nor multicomponent miRNA panels (24-mer and 7-mer) showed classification performance sufficient for a standalone diagnostic biomarker (AUC, 75%, 71%, and 53% in MOLTEST-2013, SMAC, and MOLTEST-BIS, respectively, in the 7-mer model). The performance of classification in the MOLTEST-BIS cohort with the lowest contribution of adenocarcinomas was increased when only this cancer type was considered (AUC, 60% in 7-mer model).

Adaptive filtering of microarray gene expression data based on Gaussian mixture decomposition.

BACKGROUND: DNA microarrays are used for discovery of genes expressed differentially between various biological conditions. In microarray experiments the number of analyzed samples is often much lower than the number of genes (probe sets) which leads to many false discoveries. Multiple testing correction methods control the number of false discoveries but decrease the sensitivity of discovering differentially expressed genes. Concerning this problem, filtering methods for improving the power of detection of differentially expressed genes were proposed in earlier papers. These techniques are two-step procedures, where in the first step some pool of non-informative genes is removed and in the second step only the pool of the retained genes is used for searching for differentially expressed genes. RESULTS: A very important parameter to choose is the proportion between the sizes of the pools of removed and retained genes. A new method, which we propose, allow to determine close to optimal threshold values for sample means and sample variances for gene filtering. The method is adaptive and based on the decomposition of the histogram of gene expression means or variances into mixture of Gaussian components. CONCLUSIONS: By performing analyses of several publicly available datasets and simulated datasets we demonstrate that our adaptive method increases sensitivity of finding differentially expressed genes compared to previous methods of filtering microarray data based on using fixed threshold values.

Sources of high variance between probe signals in Affymetrix short oligonucleotide microarrays.

High density oligonucleotide microarrays present a big challenge for statistical data processing methods which aim to separate changes induced by experimental factors from those caused by artifacts and measurement inaccuracies. Despite huge advances in the field of microarray probe design methods, the signal variation between probes that target a single transcript is substantially larger than their between-replicate array variability, suggesting a large influence of various probe-specific effects that introduce bias to the data. In this work we present the influence of probe-related design variations on the expression intensities of individual probes, focusing on five potential sources of high probe signal variance: the GC composition of the probe, the distance between individual probe target sites, G-quadruplex formation in the probe sequence, the occurrence of sequence motifs complementary to the oligo(dT) primer, and the specificity of unrecognized alternative splicing probeset assignment. By focusing on two high quality microarray datasets based on two distinct array designs we show the extent of variance between probes that target a specific transcript providing guidelines for the future design of microarrays and data processing methods.

Molecular testing opportunities on cytology effusion specimens: the pre-analytic effects of various body fluid cytology preparation methods on RNA extraction quality and targeted sequencing.

INTRODUCTION: RNA sequencing (RNAseq) analysis is emerging as a clinical research or diagnostic approach for cytologic samples, but there is need for formal comparison of different sample preparation methods in the cytology laboratory to identify which pre-analytic methods could provide alternatives to formalin-fixed paraffin-embedded (FFPE) sections. MATERIALS AND METHODS: We prepared 13 malignant effusions (metastatic estrogen receptor-positive breast cancer) in the cytology laboratory using 6 routine cytologic methods: FFPE cell block, Carnoy's solution, 95% ethanol (EtOH), air-dried and Diff-Quik, ThinPrep, and SurePath preparations. Measurements of RNA quality, expression of 2 multigene expression signatures, molecular subtype, and 4 common activating mutation sites in each preparation were compared with fresh frozen (FF) cell pellet in RNA preservative using distribution of fragment length and concordance correlation coefficient (CCC). RESULTS: The fraction of RNA fragments measuring 200 bases or more (DV200) were 24% higher from cytospins fixed in Carnoy's solution or 95% EtOH than DV200 from FFPE cell blocks. SurePath samples failed RNAseq quality control. There was high concordance of gene expression measurements with FF samples using cytospins fixed in Carnoy's solution, 95% EtOH, Diff-Quik (CCC = 0.829, 0.812, 0.760, respectively), or ThinPrep (CCC = 0.736), but lower using FFPE cell block (CCC = 0.564). The proportion of mutant transcripts was concordant between FF and any cytologic preparation methods. CONCLUSIONS: Cytospin preparations fixed with Carnoy's or 95% ETOH then Papanicolaou stained produced RNAseq results that were equivalent to FF samples and superior to FFPE cell block sections.

Assessment of stained direct cytology smears of breast cancer for whole transcriptome and targeted messenger RNA sequencing.

BACKGROUND: Rather than surgical resection, cytologic specimens are often used as first-line clinical diagnostic procedures due to higher safety, speed, and cost-effectiveness. Archival diagnostic cytology slides containing cancer can be equivalent to tissue biopsies for DNA mutation testing, but the accuracy of transcriptomic profiling by RNA sequencing (RNA-seq) is less understood. METHODS: This study compares the results from whole transcriptome RNA-seq and a targeted RNA-seq assay of stained cytology smears (CS) versus matched tumor tissue samples preserved fresh-frozen (FF) and processed as formalin-fixed paraffin-embedded (FFPE) sections. Cellular cytology scrapes from all 11 breast cancers were fixed and stained using three common protocols: Carnoy's (CS_C) or 95% ethanol (CS_E) fixation and then Papanicolaou stain or air-dried then methanol fixation and DiffQuik stain (CS_DQ). Agreement between samples was assessed using Lin's concordance correlation coefficient. RESULTS: Library yield for CS_DQ was too low, therefore it was not sequenced. The distributions of concordance correlation coefficient of gene expression levels in comparison to FF were comparable between CS_C and CS_E, but expression of genes enriched in stroma was lower in cytosmear samples than in FF or FFPE. Six signatures showed similar concordance to FF for all methods and two were slightly worse in CS_C and CS_E. Genomic signatures were highly concordant using targeted RNA-seq. The allele fraction of selected mutations calculated on cytosmear specimens was highly correlated with FF tissues using both RNA-seq methods. CONCLUSION: RNA can be reliably extracted from cytology smears and is suitable for transcriptome profiling or mutation detection, except for signatures of tumor stroma.

Evaluation of zero counts to better understand the discrepancies between bulk and single-cell RNA-Seq platforms.

Recent advances in sample preparation and sequencing technology have made it possible to profile the transcriptomes of individual cells using single-cell RNA sequencing (scRNA-Seq). Compared to bulk RNA-Seq data, single-cell data often contain a higher percentage of zero reads, mainly due to lower sequencing depth per cell, which affects mostly measurements of low-expression genes. However, discrepancies between platforms are observed regardless of expression level. Using four paired datasets with multiple samples each, we investigated technical and biological factors that can contribute to this expression shift. Using two separate machine learning models we found that, in addition to expression level, RNA integrity, gene or UTR3 length, and the number of transcripts potentially also influence the occurrence of zeros. These findings could enable the development of novel analytical methods for cross-platform expression shift correction. We also identified genes and biological pathways in our diverse datasets that consistently showed differences when assessed at the single cell versus bulk level to assist in interpreting analysis across transcriptomic platforms. At the gene level, 25 genes (0.12%) were found in all datasets as discordant, but at the pathway level, 7 pathways (2.02%) showed shared enrichment in discordant genes.

Short- and long-term effects of radiation exposure at low dose and low dose rate in normal human VH10 fibroblasts.

Introduction: Experimental studies complement epidemiological data on the biological effects of low doses and dose rates of ionizing radiation and help in determining the dose and dose rate effectiveness factor. Methods: Human VH10 skin fibroblasts exposed to 25, 50, and 100 mGy of 137Cs gamma radiation at 1.6, 8, 12 mGy/h, and at a high dose rate of 23.4 Gy/h, were analyzed for radiation-induced short- and long-term effects. Two sample cohorts, i.e., discovery (n = 30) and validation (n = 12), were subjected to RNA sequencing. The pool of the results from those six experiments with shared conditions (1.6 mGy/h; 24 h), together with an earlier time point (0 h), constituted a third cohort (n = 12). Results: The 100 mGy-exposed cells at all abovementioned dose rates, harvested at 0/24 h and 21 days after exposure, showed no strong gene expression changes. DMXL2, involved in the regulation of the NOTCH signaling pathway, presented a consistent upregulation among both the discovery and validation cohorts, and was validated by qPCR. Gene set enrichment analysis revealed that the NOTCH pathway was upregulated in the pooled cohort (p = 0.76, normalized enrichment score (NES) = 0.86). Apart from upregulated apical junction and downregulated DNA repair, few pathways were consistently changed across exposed cohorts. Concurringly, cell viability assays, performed 1, 3, and 6 days post irradiation, and colony forming assay, seeded just after exposure, did not reveal any statistically significant early effects on cell growth or survival patterns. Tendencies of increased viability (day 6) and reduced colony size (day 21) were observed at 12 mGy/h and 23.4 Gy/min. Furthermore, no long-term changes were observed in cell growth curves generated up to 70 days after exposure. Discussion: In conclusion, low doses of gamma radiation given at low dose rates had no strong cytotoxic effects on radioresistant VH10 cells.

Combining Low-Dose Computer-Tomography-Based Radiomics and Serum Metabolomics for Diagnosis of Malignant Nodules in Participants of Lung Cancer Screening Studies.

Radiomics is an emerging approach to support the diagnosis of pulmonary nodules detected via low-dose computed tomography lung cancer screening. Serum metabolome is a promising source of auxiliary biomarkers that could help enhance the precision of lung cancer diagnosis in CT-based screening. Thus, we aimed to verify whether the combination of these two techniques, which provides local/morphological and systemic/molecular features of disease at the same time, increases the performance of lung cancer classification models. The collected cohort consists of 1086 patients with radiomic and 246 patients with serum metabolomic evaluations. Different machine learning techniques, i.e., random forest and logistic regression were applied for each omics. Next, model predictions were combined with various integration methods to create a final model. The best single omics models were characterized by an AUC of 83% in radiomics and 60% in serum metabolomics. The model integration only slightly increased the performance of the combined model (AUC equal to 85%), which was not statistically significant. We concluded that radiomics itself has a good ability to discriminate lung cancer from benign lesions. However, additional research is needed to test whether its combination with other molecular assessments would further improve the diagnosis of screening-detected lung nodules.

POLCOVID: a multicenter multiclass chest X-ray database (Poland, 2020-2021).

The outbreak of the SARS-CoV-2 pandemic has put healthcare systems worldwide to their limits, resulting in increased waiting time for diagnosis and required medical assistance. With chest radiographs (CXR) being one of the most common COVID-19 diagnosis methods, many artificial intelligence tools for image-based COVID-19 detection have been developed, often trained on a small number of images from COVID-19-positive patients. Thus, the need for high-quality and well-annotated CXR image databases increased. This paper introduces POLCOVID dataset, containing chest X-ray (CXR) images of patients with COVID-19 or other-type pneumonia, and healthy individuals gathered from 15 Polish hospitals. The original radiographs are accompanied by the preprocessed images limited to the lung area and the corresponding lung masks obtained with the segmentation model. Moreover, the manually created lung masks are provided for a part of POLCOVID dataset and the other four publicly available CXR image collections. POLCOVID dataset can help in pneumonia or COVID-19 diagnosis, while the set of matched images and lung masks may serve for the development of lung segmentation solutions.