Next Article in Journal
The Role of NK Cells in EBV Infection and Related Diseases: Current Understanding and Hints for Novel Therapies
Previous Article in Journal
Treatment Costs and Social Burden of Pancreatic Cancer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Areas of Crush Nuclear Streaming Should Be Included as Tumor Content in the Era of Molecular Diagnostics

1
Department of Pathology, Kansai Medical University, 2-5-1 Shin-machi, Hirakata 573-1010, Osaka, Japan
2
Department of Pathology and Laboratory Medicine, Kansai Medical University Hospital, 2-3-1 Shin-machi, Hirakata 573-1191, Osaka, Japan
3
Department of Genome Analysis, Institute of Biomedical Science, Kansai Medical University, 2-5-1 Shin-machi, Hirakata 573-1010, Osaka, Japan
*
Author to whom correspondence should be addressed.
Cancers 2023, 15(6), 1910; https://doi.org/10.3390/cancers15061910
Submission received: 29 January 2023 / Revised: 14 March 2023 / Accepted: 21 March 2023 / Published: 22 March 2023
(This article belongs to the Section Methods and Technologies Development)

Abstract

:

Simple Summary

To expand the pool of samples available for genetic analysis, the quality and utility of DNA and RNA extracted from degenerated tumor tissues were examined. The DNA obtained from the nuclear streaming samples was preserved, enabling the identification of reliable variants using next-generation sequencing (NGS) analysis. The NGS metrics and variant allele frequency of the nuclear streaming samples did not differ from those of the control. This suggests that non-tumor cells exhibiting nuclear streaming could be used for genetical analysis and therefore, should be considered as tumor content to avoid misinterpreting the variant allele frequency.

Abstract

Degenerated tissues are frequently observed in malignant tumors, but are not analyzed. We investigated whether nuclear streaming and necrosis samples could be used for genetic analysis to expand the sample pool. A total of 81 samples were extracted from small cell carcinoma and lymphoma FFPE tissue blocks and classified into three histological cohorts: 33 materials with well-preserved tumor morphology, 31 nuclear streaming samples, and 17 necrosis samples. DNA and RNA integrity numbers, percentage of RNA fragments with >200 nucleotides, and next-generation sequencing quality metrics were compared among the cohorts. DNA quality did not significantly differ between nuclear streaming materials and materials with well-preserved morphology, whereas that of the necrosis samples was inferior. RNA quality decreased in the following order: materials with well-preserved morphology > nuclear streaming > necrosis. The sequencing metrics did not differ significantly between the nuclear streaming samples and materials with well-preserved morphology, and reliable variants were detected. The necrosis samples extracted from resections exhibited sequencing failure and showed significantly fewer on-target aligned reads and variants. However, variant allele frequency did not differ among the cohorts. We revelated that DNA in nuclear streaming samples, especially within biopsies, could be used for genetic analysis. Moreover, degenerated non-tumor cells should be counted when evaluating tumor content to avoid misinterpreting the variant allele frequency.

1. Introduction

Next-generation sequencing (NGS) to detect molecular alterations has emerged as a precise approach for assessing malignant tumors [1,2,3,4,5,6,7,8]. Although surgical resection is the most effective therapy for early-stage solid malignancies, some advanced cases are inoperable. In such cases, histological and molecular diagnoses are conducted using small amounts of materials, such as biopsies and/or cytological specimens [7]. Although liquid components of blood are available for NGS analysis, tissue and cell materials are still required, as they allow for larger numbers of panels and provide more accurate analytical results than liquid materials [7].
The amount of DNA required for NGS is 50–300 ng [4,7], depending on the platform and target enrichment method. Larger gene panels require a greater amount of materials, and low DNA input negatively affects the sensitivity of variant detection [1,2,3,4,5,6,7]. However, DNA yield depends on the size and cellularity of the selected tumor [7,8,9,10]. The expected average DNA yield from a single nucleated cell is 6 pg. Therefore, to obtain 10 ng DNA, approximately 2000 whole cells are required [4]. Sample quantity can be increased when tumor areas with degradation, such as areas of necrosis or crushed nuclear streaming, which have previously been considered difficult to analyze via NGS [7,8,9,10], are considered suitable for molecular diagnosis.
The recommended histological selection criteria for target sequencing (TS) include at least 20% tumor content without degenerative features; accordingly, areas with tumor necrosis and areas abundant in wild-type non-tumor cells (e.g., inflammatory cells) are excluded [7,8,9,10]. Increasing tumor content and DNA yield would increase the NGS success rate [11]; however, the quality of DNA and RNA in tumor tissues exhibiting degenerative changes has not been examined. In fact, necrotic and crushed samples encountered in routine pathology work are excluded from genetic analysis, despite tumor samples derived from patients with malignant tumors, who specifically require genetic treatment, frequently exhibiting such histological changes. Determining whether necrosis and crushed nuclear streaming samples exhibit genomic integrity could facilitate the improvement of tumor content estimation.
In this study, we examined the quality of nucleic acids extracted from histologically degenerated tissues, such as samples with necrosis and nuclear streaming, and performed targeted NGS and evaluated sequencing metrics to confirm whether these samples could be used in genetic therapy. We evaluated their DNA integrity number (DIN), RNA integrity number (RIN), and percentage of RNA fragments with >200 nucleotides (DV200), which are commonly used indices to assess nucleic acid quality. For these evaluations, we used materials from small cell lung carcinoma (SC) and diffuse large B-cell lymphoma (DLBCL), in which morphological nuclear degeneration is frequently observed. Our findings can potentially help in increasing the number of samples available for genetic analysis, thereby improving diagnosis and treatment outcomes.

2. Materials and Methods

2.1. Patient Eligibility and Sample Preparation

We classified 17 SC and 28 DLBCL samples collected at the Kansai Medical University Hospital (Osaka, Japan) between 2017 and 2020 into three histological areas: material with a well-preserved tumor morphological component (solid tumor component with a clear nuclear shape and recognizable outline; Figure 1A), material with nuclear streaming components (component showing linear or lumpy basophilic nuclei, with no clear cell to cell border, no clear outline of individual nuclei, no cell membrane, and no eosinophilic cytoplasm, usually considered to be histological “artifacts”; Figure 1B), and material containing necrosis components (solid tumor component with an eosinophilic morphology, lacking basophilic nuclei, a clear outline of the cell membrane, and cell adhesion; Figure 1C). The three histological components were extracted from the 45 cases that were included in the study. Histological representative areas with >10% wild-type non-tumor cells and <90% tumor cells were excluded from this study. If all three components were unavailable in the same case, or the amount of the components present in that case was low, only the available components were included in the analysis. The final cohort comprised 33 samples without degeneration (8 SC and 25 DLBCL), including 31 with nuclear streaming (13 SC and 18 DLBCL) and 17 with necrosis (3 SC and 14 DLBCL).
Moreover, we examine whether the sample volume under the formalin fixed process used in the present study influences the quality of nucleic acid. Considering that the fixation time in the present study ranged from 12–24 h, the penetration rate of 10% neutral-buffered formalin was 1 mm/h [7], and the SC and diffuse large B-cell lymphoma (DLBCL) samples were divided into three groups, according to the sample size, according to diameter and height: small (≤12 × 12 mm), large (≥24 × 24 mm), and medium (with dimensions other than those of the small and large groups).
This study was approved by the Institutional Review Board of the Kansai Medical University Hospital (approval no. 2020271). Informed consent was obtained from all patients, and the opt-out approach was applied owing to the retrospective design of the study, with no new risk to the participants. Information regarding this study, such as the inclusion criteria and the opportunity to opt-out, was provided on the hospital’s website. The inclusion and exclusion criteria, sample selection, and sample preparation are described in Supplemental Information S1–S4.

2.2. Assessment of DNA and RNA Quality and Variant Detection Using NGS

DNA and RNA were isolated from each sample, and their purities were measured; their quality was based on RIN, DIN, and DV200. DNA was isolated using a QIAamp DNA formalin-fixed, paraffin-embedded (FFPE) Tissue Kit (Qiagen, Hilden, Germany). RNA was isolated using an miRNeasy FFPE Tissue Kit (Qiagen). The purity of the nucleic acid samples was assessed based on their spectrophotometric absorbance at 260–280 nm (A260/A280) using a NanoDrop 3300 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). RIN, DIN, and DV200 indices and the concentrations of the extracted DNA and RNA samples were quantified using TapeStation 4150 (2015, Agilent Technologies, Santa Clara, CA, USA) based on the genomic DNA (gDNA) and RNA ScreenTape assays. Owing to instrument performance, RIN, DV200, and DIN can provide quality measurements in samples with RNA concentrations >2.0 ng/μL and gDNA concentrations >3.0 ng/μL. Statistical analysis was performed on cases with evaluable RIN, DV200, and DIN. Subsequently, to asses DNA quality and detect variants, NGS was conducted using an MiSeq sequencing platform. The AmpliSeq for Illumina Cancer HotSpot Panel v2 (50 target genes, 207 amplicons; Illumina Inc., San Diego, CA, USA) protocol was applied (Supplemental Information S5).

2.3. NGS Data Analysis and Validation of Observed Variants

NGS data analysis is summarized in Supplemental Information S6. Recommended total read depth corresponding to variant allele frequency (VAF) was calculated using https://github.com/mvasinek/olgen-coverage-limit, accessed on 28 January 2022 [12]. When considering clinical reliability, samples with an amplicon mean coverage of <500× were considered non-informative for detecting somatic mutations with a frequency of 5% [8,13].

2.4. Statistical Analysis

Because DIN, RIN, and DV200 did not differ between DLBCL and SC samples (p > 0.05, Supplemental Information S7), these two tumors were analyzed together. Differences in DIN, RIN, and DV200 were assessed using Welch’s t-test, followed by the F test. The percentage of variants that passed the filter with a quality of 100 was calculated as:
n u m b e r o f v a r i a n t s t h a t p a s s e d t h e f i l t e r w i t h a q u a l i t y o f 100 n u m b e r o f a l l d e t e c t e d v a r i a n t s
The percentage of variants that did not pass the filter or did not show a quality of 100 was calculated as:
n u m b e r o f v a r i a n t s t h a t d i d n o t p a s s t h e f i l t e r + n u m b e r o f r e a d s w h i c h d i d n o t h a v e a q u a l i t y o f 100 n u m b e r o f a l l d e t e c t e d v a r i a n t s
Chi-square (χ2) analysis and Tukey–Kramer or Gomes–Howell tests were used to assess the relationships between materials without degeneration, those with nuclear streaming, and those with necrosis, as well as the differences based on the sample volume (small, medium, and large). Spearman’s correlation was used to assess the relationship between the percentage of on-target aligned reads and DIN. Statistical analysis was performed using SPSS (v20.0; IBM Corp., Armonk, NY, USA), with significance set at p < 0.05.

3. Results

3.1. DIN, RIN, and DV200 of Materials without Degeneration, with Nuclear Streaming, and with Necrosis

3.1.1. DIN

DIN was available from 71 samples (range = 1.6–5.5, mean = 2.8 ± 0.12, Supplemental Information S2, S3, and S8). DIN was significantly lower in necrosis samples than in nuclear streaming samples (p < 0.05, Table 1, and Figure 2A). Although no significant difference in DIN was observed between necrosis samples and samples with well-preserved morphology, necrosis samples showed the lowest average DIN among cohorts (average DIN: materials without degeneration 2.78 ± 0.69; nuclear streaming 3.16 ± 1.12; necrosis 2.29 ± 0.53).

3.1.2. RIN

RIN was available from 72 samples (range = 1.1–3.7, mean = 1.87 ± 0.06, Supplemental Information S2 and S3). Nuclear streaming and necrosis samples showed a lower RIN than materials with well-preserved morphology (vs. nuclear streaming, p = 0.05; vs. necrosis, p = 0.026; Table 1 and Figure 2B). The average RIN decreased in the following order: materials with well-preserved morphology (2.05 ± 0.60) > nuclear streaming samples (1.76 ± 0.32) > necrosis samples (1.68 ± 0.28).

3.1.3. DV200

DV200 was available from 79 samples (range = 34.89–90.77, mean = 62.7 ± 1.64, Supplemental Information S2 and S3). The average DV200 value of the DLBCL and SC samples decreased in the following order: materials with well-preserved morphology (68.3 ± 16.4) > nuclear streaming samples (59.6 ± 13.7) > necrosis samples (57.5 ± 8.16; Table 1 and Figure 2C). Nuclear streaming and necrosis samples showed a significantly lower DV200 than materials with well-preserved morphology (both p < 0.05).

3.2. Results of NGS via AmpliSeq for Illumina Cancer HotSpot Panel v2

The average A260/A280 ratio in the three histological cohorts was within the reference range (1.8–2.0), with no significant difference between the cohorts (all p > 0.05, Supplemental Information S9). RIN values indicated that the samples were too degraded for NGS analysis; RIN > 8 is regarded as usable [7]. Therefore, TS was performed with AmpliSeq, which can be achieved using low-quality DNA.
TS was performed by selecting cases with high and low DIN values. Considering that DNA degradation during storage can affect the TS results [7], TS was also performed using two SC samples obtained in 2015. The NGS results, metrics, and data of the detected variants are summarized in Table 2 and Table 3, and the detailed data on the FASTQ file, and the detected variants are presented in Supplemental Information S10 and S11.
A total of 25 samples (10 materials with well-preserved morphology, 9 with nuclear streaming, and 6 with necrosis) were obtained from 10 cases and evaluated using NGS. Both samples obtained in 2015 showed an amplicon coverage >500×, suggesting that the storage time did not greatly affect nucleic acid quality. However, amplicons with a coverage <500× were detected in two necrosis samples (SC2018-2, amplicon mean coverage 29.2×; DL2018-19, 82.6×; Figure 3A,B).
Statistical analyses of the 25 samples revealed significantly lower percentages of on-target aligned reads for necrosis samples than for materials with well-preserved morphology and nuclear streaming samples (materials with well-preserved morphology, 95.7 ± 3.02; with nuclear streaming, 94.6 ± 6.44; and with necrosis, 54.1 ± 40.6; all p < 0.05, Table 4 and Figure 3C). No significant differences in other NGS metrics were observed among materials with well-preserved morphology, those with nuclear streaming, or those showing necrosis (all p > 0.05).
The number of variants detected is shown in Table 3. All 25 samples exhibited the recommended read depth corresponding to VAF. Informative variants were detected in 23 cases, including nuclear streaming and necrosis samples; these included variants in KDR, TP53, KIT, PTEN, KRAS, ERBB4, MET, STK11 in SCs, and TP53, PDGFR, and KIT in DLBCLs. The percentage of variants that passed the filter with a quality of 100, percentage of variants that did not achieve a quality of 100 (thus not passing the filter), and VAF did not differ significantly among the three histological groups (all p > 0.05, Table 4). However, the percentage of variants that passed the filter with a quality of 100 detected using TS decreased in the following order: materials with well-preserved morphology (81.7 ± 15.6) > nuclear streaming samples (80.0 ± 18.9) > necrosis samples (60.8 ± 30.8) (all p > 0.05). Furthermore, the percentage of variants that did not pass the filter or achieve a quality of 100 increased in the following order: materials with well-preserved morphology (18.3 ± 15.6) < nuclear streaming samples (20.2 ± 17.9) < necrosis samples (39.2 ± 30.8) (all p > 0.05). TP53 mutations were detected in all 10 cases (materials with well-preserved morphology: 100% (9/9), nuclear streaming samples: 100% (9/9), and necrosis samples: 50% (2/4). The VAFs of TP53 mutations in all histological cohorts did not differ significantly, being approximately 80% for all cohorts, and occurred in the following order: materials with well-preserved morphology (84.0 ± 22.7) > nuclear streaming samples (83.8 ± 28.7) > necrosis samples (85.8 ± 28.5) (all p > 0.05). The high VAFs suggest that the NGS metrics reflect the quality of DNA obtained within the tumor and are not derived from the non-tumor component.

3.3. Relationship between Performance in NGS Analysis Using MiSeq and Histological Changes

In the 25 samples, DIN and the percentage of on-target aligned reads were positively correlated (r = 0.671, p = 0.003, Figure 3D). All samples except four—one nuclear streaming (DL2018-14) and three necrosis (SC2018-2, DL2018-14, and DL2018-19) samples—showed >80% on-target aligned reads, indicating good alignment with the target region. Although the four samples with <80% on-target aligned reads were all resections, all materials with well-preserved morphology in the same case showed >80% aligned reads, indicating that the results were influenced by the degenerative status in each of the cohorts.

3.4. Differences in DIN, RIN, and DV200 Based on Sample Volume

As the samples with <80% on-target aligned reads were all resections, we investigated whether the formalin fixed tissue volumes influenced the nucleic acid quality. DIN was significantly higher in the group with small tissue volumes compared to that in the groups with medium and large tissue volumes (both p < 0.05, Figure 4A, Supplemental Information S12). No significant differences in RIN were detected among the three groups, but the mean RIN in the groups with large, medium, and small sample volumes was 2.05, 2.03, and 1.78, respectively (Figure 4B). DV200 was significantly higher in the group with small tissue volumes compared to that in the groups with medium and large tissue volumes (both p < 0.05, Figure 4C). These results indicate that the DNA extracted from smaller tissue samples and RNA extracted from larger tissue samples were better for genetical analysis during 12–24 h of formalin fixation time.

4. Discussion

The DNA quality of nuclear streaming samples did not significantly differ from that of materials with well-preserved morphology, whereas necrosis samples were inferior to the other two histological groups. The quality of both DNA and RNA was poor in necrosis samples. The NGS metrics of the nuclear streaming samples did not differ from those of materials with well-preserved morphology, and informative variants were detected. The percentage of on-target aligned reads and variants passing through the filter with high quality was lower in necrosis samples than in the other two sample types. Although some informative variants were detected in the necrosis samples, all failed samples were from the necrosis group, and the quality of DNA was better in smaller samples, such as biopsies.
DIN is a scale used to assess the quality of DNA, and it ranges from 1 (completely degraded) to 10 (completely intact), based on gDNA degradation [7,10,14]. Samples with DIN > 3 are recommended for NGS analysis [7]. The DINs of the samples without degeneration and the nuclear streaming samples were close to 3, similar to those of samples evaluated in the last 4 years, according to the Japanese guidelines (DIN 3.77 ± 1.77) [7,8,10]. The DINs of the necrosis samples were significantly lower. The difference in DIN between nuclear streaming and necrosis samples may be associated with their degenerative processes; necrosis is a biological response that leads to nucleic acid fragmentation [4], whereas nuclear streaming is a histological change caused by sampling stimulation [15,16,17,18]. To enable genetic analysis, particularly for small samples, the nuclear streaming component should be included in the material, without trimming during the pre-analysis evaluation process.
The percentage of on-target aligned reads and variants passing through the filter with high quality was lower in the necrosis samples than in the materials with well-preserved morphology and nuclear streaming samples. Some informative variants were detected in the necrosis samples; however, most samples with an amplicon mean coverage <500× and on-target aligned reads <80% that were considered unsuitable for analysis were necrotic [8]. Because nucleic acid quality most strongly influences NGS quality metrics [1,3,5,6,7], degraded DNA limits the depth of coverage and number of supporting reads, leading to limited variant detection [1,19]. These features are consistent with those of necrosis samples, but not nuclear streaming samples. Only one study has evaluated the usefulness of necrotic tumor samples using TS, reporting that DNA samples obtained from the area of necrosis in colorectal adenocarcinomas are comparable to those of viable tumors [4,20]. However, that was a preliminary finding, and no detailed description of the research results is available [20]. Our results indicate that variants might be identified from necrotic DNA, but these results should be carefully evaluated for false negatives. However, nuclear streaming samples should be considered, as they may contain preserved target DNA and provide genetic information to facilitate therapy.
Furthermore, DIN and DV200 were significantly different regarding formalin fixed tissue sample volumes. Formalin promotes the RNA fragmentation, but inhibits the DNA degradation [8,10]. However, the quality and yield of DNA in large surgical specimens tend to be poorer and less than those of biopsy specimens due to the poor fixation status [10]. Therefore, when using DNA in nuclear streaming samples for genetical analysis, it is preferable to extract it from small tissue sample, such as biopsies.
Although the DNA quality of the degenerative component was low, the VAF and total passing filter (PF) reads of this component did not differ from those of samples with well-preserved morphology. VAF is expressed as the percentage of sequence reads observed by matching a specific mutation derived from the overall coverage of a given locus [21,22]. When used in conjunction with the percentage of cancer cells, VAF can be beneficial in interpreting the mutation type and distribution [7,21,22]. In histological evaluations, non-tumor cells with degenerated features are not counted as tumor content because their DNA is considered incapable of amplification, and it does not enhance the total PF read number. However, our results indicate that not counting degenerated tumor cells will falsely elevate the VAF. Therefore, degenerated non-tumor cells should preferably be assessed when determining the tumor percentage to avoid underestimating the VAF.
Consistent with a previous study, we observed a positive correlation between the percentages of on-target aligned reads and DIN [14]. DIN is a useful indicator of the ability to capture target areas. Particularly, determining DIN is cost-effective and less invasive than re-sampling, as it can be performed during the pre-analysis process. Our results also suggest that evaluating DIN for limited samples with histological degeneration may be necessary to provide more optimal treatment.
Comparative analysis based on RIN and DV200 showed that RNA quality in nuclear streaming and necrosis samples was significantly inferior to materials with well-preserved morphology. In NGS analysis, samples with a RIN > 8.0 are of acceptable quality, whereas those with a DV200 < 30% are considered unsuitable for analysis [7,8,17]. Herein, the RIN and DV200 of the three histological groups did not meet the criteria for NGS, consistent with a previous study conducted in Japan [8] (present RIN, 2.05 ± 0.60; previous RIN, 2.23 ± 0.45). As RNA is less stable than DNA [23], extraction from specimens frozen in liquid nitrogen is recommended [7]; however, this method is not routinely used in Japan. Indicative of poor RNA quality, significantly lower RIN and DV200 values were observed for nuclear streaming and necrosis samples; therefore, the use of degenerated samples for RNA-based analyses was challenging. Therefore, DNA-based analysis should be prioritized for samples with histological degeneration, or nucleic acid quality evaluation should be performed before NGS.
This study had some limitations. First, the DNA quality in the samples was low. Furthermore, we were unable to extract all three tissue materials from same patient in all cases, and some samples exceeded the recommended ischemic time of 3 h (Supplemental Information S4) [7]. However, a comparison of DIN revealed significant differences among cohorts, and this result is consistent with the NGS metrics of high VAF; the DIN was also positively correlated with the percentages of on-target aligned reads. Therefore, the results can be considered reliable. Second, we did not examine whether RNA extracted from degenerated materials can be used for NGS. Third, we did not make any deep references to germline or somatic loss-of-function mutations in the gene of interest. Future studies should aim to use samples with better nucleic acid quality; obtain all three tissue materials from the same patients in all cases; and employ different library preparations, various genetic analysis modalities, and several variant types.
Based on the results of this study, we recommend pre-analytically evaluating the DIN of nuclear streaming samples and using their DNA for NGS when the sample volume is insufficient. Moreover, if samples with degenerated non-tumor cells, such as lymphocytes, are not completely dissected, wild-type gene reads derived from non-tumor cells could increase the total PF reads and may be misinterpreted as the tumor variant having a lower VAF.

5. Conclusions

This is the first study to investigate the use of histologically degenerated samples to determine if these can be used to increase the number of samples available, thus improving diagnosis and treatment. We suggest that DNA extracted from the well-fixed nuclear streaming component can provide useful sequencing data for some target regions, and that degenerated non-tumor cells showing nuclear streaming should be counted when evaluating tumor content to avoid misinterpreting the variant allele frequency.

Supplementary Materials

The following supporting information can be downloaded at: https://drive.google.com/drive/folders/10tJ1voQbhMpYJs3d1h8l0Jnud5lCznii?usp=sharing, Supplemental Information S1. Case and sample selection algorithm for analyzing nucleic acid quality. Supplemental Information S2. Status of the small cell lung carcinoma samples. Supplemental Information S3. Status of the diffuse large B-cell lymphoma samples. Supplemental Information S4. Sample preparation status. Supplemental Information S5. Assessment of DNA quality and variant detection using NGS. Supplemental information S6. NGS data analysis and validation of observed variants. Supplemental Information S7. Differences in DNA integrity number, RNA integrity number, and DV200 based on tumor histology. Supplemental Information S8. DNA integrity data analyzed using the Tape Station 4150 system. Supplementary Information S9. Concentration comparison of purity based on the A260/A280 ratio (DNA) among the three histological groups. Supplemental Information S10. FASTq file details of next-generation sequencing (Excel). Supplemental Information S11. Detected variants of next-generation sequencing. Supplemental Information S12. Analysis of differences in DNA integrity number, RNA integrity number, and DV200 based on sample volume.

Author Contributions

All authors contributed to the conception and design of the study. Y.N. developed the theory and performed the computations. R.Y., N.A. and K.H. supported the preparation of tissue specimens and NGS analysis procedures. K.T. conceived the idea and supervised the study. The first draft of the manuscript was written by Y.N., and R.Y., N.A., K.H., K.T. and Y.N. commented on the previous version of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

This study was approved by the Institutional Review Board of the Kansai Medical University Hospital (approval no. 2020271).

Informed Consent Statement

Informed consent was obtained from all patients, and the opt-out approach was applied, owing to the retrospective design of the study, with no new risk to the participants. Information regarding this study, such as the inclusion criteria and the opportunity to opt-out, was provided on the hospital’s website. The inclusion and exclusion criteria, sample selection, and sample preparation are described in Supplemental Information S1–S4.

Data Availability Statement

All data were included in manuscript and Supplemental Information.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bettoni, F.; Koyama, F.C.; Carpinetti, P.D.A.; Galante, P.A.F.; Camargo, A.A.; Asprino, P.F. A straightforward assay to evaluate DNA integrity and optimize next-generation sequencing for clinical diagnosis in oncology. Exp. Mol. Pathol. 2017, 103, 294–299. [Google Scholar] [CrossRef] [PubMed]
  2. Beigh, M. Next-Generation Sequencing: The Translational Medicine Approach from “Bench to Bedside to Population”. Medicines 2016, 3, 14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Endrullat, C.; Glökler, J.; Franke, P.; Frohme, M. Standardization and quality management in next-generation sequencing. Appl. Transl. Genom. 2016, 10, 2–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Chen, H.; Luthra, R.; Goswami, R.; Singh, R.; Roy-Chowdhuri, S. Analysis of Pre-Analytic Factors Affecting the Success of Clinical Next-Generation Sequencing of Solid Organ Malignancies. Cancers 2015, 7, 1699–1715. [Google Scholar] [CrossRef] [PubMed]
  5. A Cree, I.; Deans, Z.; Ligtenberg, M.J.L.; Normanno, N.; Edsjö, A.; Rouleau, E.; Solé, F.; Thunnissen, E.; Timens, W.; Schuuring, E.; et al. Guidance for laboratories performing molecular pathology for cancer patients. J. Clin. Pathol. 2014, 67, 923–931. [Google Scholar] [CrossRef] [Green Version]
  6. Nagahashi, M.; Shimada, Y.; Ichikawa, H.; Nakagawa, S.; Sato, N.; Kaneko, K.; Homma, K.; Kawasaki, T.; Kodama, K.; Lyle, S.; et al. Formalin-fixed paraffin-embedded sample conditions for deep next generation sequencing. J. Surg. Res. 2017, 220, 125–132. [Google Scholar] [CrossRef]
  7. Kanai, Y.; Nishihara, H.; Miyagi, Y.; Tsuruyama, T.; Taguchi, K.; Katoh, H.; Takeuchi, T.; Gotoh, M.; Kuramoto, J.; Arai, E.; et al. The Japanese Society of Pathology Guidelines on the handling of pathological tissue samples for genomic research: Standard operating procedures based on empirical analyses. Pathol. Int. 2018, 68, 63–90. [Google Scholar] [CrossRef]
  8. Fujii, T.; Uchiyama, T.; Matsuoka, M.; Myojin, T.; Sugimoto, S.; Nitta, Y.; Okabe, F.; Sugimoto, A.; Sekita-Hatakeyama, Y.; Morita, K.; et al. Evaluation of DNA and RNA quality from archival formalin-fixed paraffin-embedded tissue for next-generation sequencing—Retrospective study in Japanese single institution. Pathol. Int. 2020, 70, 602–611. [Google Scholar] [CrossRef]
  9. Kuwata, T.; Wakabayashi, M.; Hatanaka, Y.; Morii, E.; Oda, Y.; Taguchi, K.; Noguchi, M.; Ishikawa, Y.; Nakajima, T.; Sekine, S.; et al. Impact of DNA integrity on the success rate of tissue-based next-generation sequencing: Lessons from nationwide cancer genome screening project SCRUM-Japan GI-SCREEN. Pathol. Int. 2020, 70, 932–942. [Google Scholar] [CrossRef]
  10. Hatanaka, Y.; Kuwata, T.; Morii, E.; Kanai, Y.; Ichikawa, H.; Kubo, T.; Hatanaka, K.C.; Sakai, K.; Nishio, K.; Fujii, S.; et al. The Japanese Society of Pathology Practical Guidelines on the handling of pathological tissue samples for cancer genomic medicine. Pathol. Int. 2021, 71, 725–740. [Google Scholar] [CrossRef]
  11. Goswami, R.; Luthra, R.; Singh, R.R.; Patel, K.P.; Routbort, M.J.; Aldape, K.D.; Yao, H.; Dang, H.D.; Barkoh, B.A.; Manekia, J.; et al. Identification of Factors Affecting the Success of Next-Generation Sequencing Testing in Solid Tumors. Am. J. Clin. Pathol. 2016, 145, 222–237. [Google Scholar] [CrossRef]
  12. Petrackova, A.; Vasinek, M.; Sedlarikova, L.; Dyskova, T.; Schneiderova, P.; Novosad, T.; Papajik, T.; Kriegova, E. Standardization of Sequencing Coverage Depth in NGS: Recommendation for Detection of Clonal and Subclonal Mutations in Cancer Diagnostics. Front. Oncol. 2019, 9, 851. [Google Scholar] [CrossRef] [PubMed]
  13. D’Haene, N.; Le Mercier, M.; De Nève, N.; Blanchard, O.; Delaunoy, M.; El Housni, H.; Dessars, B.; Heimann, P.; Remmelink, M.; Demetter, P.; et al. Clinical Validation of Targeted Next Generation Sequencing for Colon and Lung Cancers. PLoS ONE 2015, 10, e0138245. [Google Scholar] [CrossRef]
  14. Morlote, D.; Janowski, K.M.; Siniard, R.C.; Guo, R.J.; Winokur, T.; DeFrank, G.; Harada, S. Effects of Improved DNA Integrity by Punch From Tissue Blocks as Compared to Pinpoint Extraction From Unstained Slides on Next-Generation Sequencing Quality Metrics. Am. J. Clin. Pathol. 2019, 152, 27–35. [Google Scholar] [CrossRef] [PubMed]
  15. Travis, W.D. Update on small cell carcinoma and its differentiation from squamous cell carcinoma and other non-small cell carcinomas. Mod. Pathol. 2012, 25 (Suppl. S1), S18–S30. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Davenport, R.D. Diagnostic value of crush artifact in cytologic specimens. Occurrence in small cell carcinoma of the lung. Acta Cytol. 1990, 34, 502–504. [Google Scholar]
  17. Wehmas, L.C.; E Wood, C.; Chorley, B.N.; Yauk, C.; Nelson, G.M.; Hester, S.D. Enhanced Quality Metrics for Assessing RNA Derived From Archival Formalin-Fixed Paraffin-Embedded Tissue Samples. Toxicol. Sci. 2019, 170, 357–373. [Google Scholar] [CrossRef]
  18. Flamminio, F.; Tosi, A.L.; Fellegara, G. Crushing Artifacts Resulting in Small Blue Cellular Clusters That Simulate Small Cell Carcinoma. Int. J. Surg. Pathol. 2011, 19, 487–491. [Google Scholar] [CrossRef]
  19. Dehghani, M.; Rosenblatt, K.P.; Li, L.; Rakhade, M.; Amato, R.J. Validation and Clinical Applications of a Comprehensive Next Generation Sequencing System for Molecular Characterization of Solid Cancer Tissues. Front. Mol. Biosci. 2019, 6, 82. [Google Scholar] [CrossRef]
  20. Jabbar, K.; Routbort, M.; Singh, C.; Rashid, A.; Broaddus, R. Impact of tumor necrosis on success of clinical next generation sequencing. In Laboratory Investigation; Nature Publishing Group: New York, NY, USA, 2015; p. 501A. [Google Scholar]
  21. Li, Y.; Xie, X. Deconvolving tumor purity and ploidy by integrating copy number alterations and loss of heterozygosity. Bioinformatics 2014, 30, 2121–2129. [Google Scholar] [CrossRef] [Green Version]
  22. Strom, S.P. Current practices and guidelines for clinical next-generation sequencing oncology testing. Cancer Biol. Med. 2016, 13, 3–11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Maiväli, Ü.; Paier, A.; Tenson, T. When stable RNA becomes unstable: The degradation of ribosomes in bacteria and beyond. Biol. Chem. 2013, 394, 845–855. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Histological features of the materials with a well-preserved tumor morphological component (A), nuclear streaming samples (B) (the arrows show the shape-preserved cells), and necrotic samples (C). (AC) Hematoxylin and eosin staining; original magnification, ×200; bar, 100 μm.
Figure 1. Histological features of the materials with a well-preserved tumor morphological component (A), nuclear streaming samples (B) (the arrows show the shape-preserved cells), and necrotic samples (C). (AC) Hematoxylin and eosin staining; original magnification, ×200; bar, 100 μm.
Cancers 15 01910 g001
Figure 2. Differences in DIN (A), RIN (B), and DV200 (C) among the materials with well-preserved morphology, nuclear streaming, and necrosis groups in SC and DLBCL samples. DIN, DNA integrity number; RIN, RNA integrity number; SC, small cell lung carcinoma; DLBCL, diffuse large-B-cell lymphoma. (*, p < 0.05).
Figure 2. Differences in DIN (A), RIN (B), and DV200 (C) among the materials with well-preserved morphology, nuclear streaming, and necrosis groups in SC and DLBCL samples. DIN, DNA integrity number; RIN, RNA integrity number; SC, small cell lung carcinoma; DLBCL, diffuse large-B-cell lymphoma. (*, p < 0.05).
Cancers 15 01910 g002
Figure 3. Amplicons with a coverage <×500 (A) (necrosis in SC 2018-2). All samples that failed NGS were necrotic (B). Differences in percentage of on-target aligned reads among three histological groups (C). Positive correlation between the percentage of on-target aligned reads and DIN (D) (green, materials with well-preserved tumor morphology; red, nuclear streaming; blue, necrosis). DIN, DNA integrity number; SC, small cell lung carcinoma; DLBCL, diffuse large-B-cell lymphoma. (**, p < 0.01).
Figure 3. Amplicons with a coverage <×500 (A) (necrosis in SC 2018-2). All samples that failed NGS were necrotic (B). Differences in percentage of on-target aligned reads among three histological groups (C). Positive correlation between the percentage of on-target aligned reads and DIN (D) (green, materials with well-preserved tumor morphology; red, nuclear streaming; blue, necrosis). DIN, DNA integrity number; SC, small cell lung carcinoma; DLBCL, diffuse large-B-cell lymphoma. (**, p < 0.01).
Cancers 15 01910 g003
Figure 4. The relationships of DIN (A), RIN (B), and DV200 (C) relative to the sample volume: small (≤12 × 12 mm), large (≥24 × 24 mm), and medium (with dimensions other than those of the small and large groups) in SC and DLBCL samples. DIN, DNA integrity number; RIN, RNA integrity number; SC, small cell lung carcinoma; DLBCL, diffuse large-B-cell lymphoma (*, p < 0.05; **, p < 0.01).
Figure 4. The relationships of DIN (A), RIN (B), and DV200 (C) relative to the sample volume: small (≤12 × 12 mm), large (≥24 × 24 mm), and medium (with dimensions other than those of the small and large groups) in SC and DLBCL samples. DIN, DNA integrity number; RIN, RNA integrity number; SC, small cell lung carcinoma; DLBCL, diffuse large-B-cell lymphoma (*, p < 0.05; **, p < 0.01).
Cancers 15 01910 g004
Table 1. Comparison of DNA integrity number, RNA integrity number, and DV200 between the three histological groups.
Table 1. Comparison of DNA integrity number, RNA integrity number, and DV200 between the three histological groups.
nMeanSDSE95% CIvs. Without
Degeneration
p-Value
vs. Nuclear
Streaming
p-Value
vs. Necrosis
p-Value
DIN (n = 71)
Well-preserved morphology302.780.687020.125432.523–3.037-0.2220.212
Nuclear streaming283.160711.119630.211592.727–3.5930.222-0.011
Necrosis132.292300.528300.146521.973–2.6120.2130.011
RIN (n = 72)
Well-preserved morphology312.051610.596020.107041.833–2.270-0.050.026
Nuclear streaming301.756660.316970.057871.638–1.8750.05-0.747
Necrosis111.681810.278630.084011.495–1.8690.0260.747-
DV200 (n = 79)
Well-preserved morphology3268.3409316.404202.8998862.4266–74.2553-0.0410.031
Nuclear streaming3059.5603313.705922.3953454.4425–64.68720.041-0.878
Necrosis1757.494118.157001.9783653.3002–61.68810.0310.878-
n, number of samples; σ 2, non-diversification; SD, standard deviation; SE, standard error; CI, confidence interval; DIN, DNA integrity number; RIN, RNA integrity number; DV200, percentage of RNA fragments >200 nt; well-preserved morphology, materials with well-preserved tumor morphology; nuclear streaming, materials with nuclear streaming; necrosis, materials with necrosis; bold, statistically significant values.
Table 2. Overall sequencing results for the 25 samples obtained from 10 cases.
Table 2. Overall sequencing results for the 25 samples obtained from 10 cases.
Year-No.StatusOrganStatusPercent
Q30 Bases
Total PF ReadsPercentage On-target
Aligned Reads
Uniformity of Coverage
[Pct > 0.2 × Mean]
Amplicon Mean CoverageDIN
(≥2.3)
SC
2015-1
Ruteruswell-preserved morphology90.3690,94295.8193.242803.92.0
nuclear streaming90.37396,17692.3390.8215471.2
necrosis89.64543,10071.4860.871494.51.4
SC
2015-2
Blungwell-preserved morphology79.69695,09888.7967.151253.52.3
nuclear streaming80.92654,33897.3582.611358.22.5
SC
2018-2
Rlungwell-preserved morphology80.63553,47497.1588.891173.5-
nuclear streaming79.97482,59096.9785.991009.72.4
necrosis86.36270,7543.1978.7429.21.7
SC
2018-3
Rlungwell-preserved morphology90.33677,40096.4575.362737.32.8
necrosis90.84668,77096.5271.012673.82.5
SC
2018-4
Rbrainwell-preserved morphology79.74665,55895.385.991347.4-
nuclear streaming79.88501,74878.1482.61780.72.1
SC
2019-8
Rbrainwell-preserved morphology80.94662,30097.0982.131402.62.2
nuclear streaming92.26396,97696.1866.671626.72.1
SC
2019-9
Buteruswell-preserved morphology91.89949,21497.3988.413901.84.1
nuclear streaming81.02945,08297.8793.241947.3-
DL
2018-14
Rintestinewell-preserved morphology80.36641,64298.4783.571375.62.6
nuclear streaming80.84648,07298.1479.231387.92.7
necrosis76.66737,33253.6482.13756.42.5
DL
2018-19
Rbrainwell-preserved morphology80.46421,50692.2379.23845.92
nuclear streaming80.13593,10295.9382.131236.5-
necrosis75.26517,6068.0576.3382.61.9
DL
2020-2
Rcolonwell-preserved morphology80.83464,44698.3480.681003.72.4
nuclear streaming80.25645,40498.4687.441393.12.5
necrosis79.49555,93491.6482.611079.71.9
SC, small cell lung carcinoma; DL, diffuse large B-cell lymphoma; R, resection; B, biopsy; PF, passing filter reads; DIN, DNA integrity number; well-preserved morphology, materials with well-preserved tumor morphology; nuclear streaming, materials with nuclear streaming; necrosis, materials with necrosis; -, unable to measure because the concentration was below the measurement range of the instrument; bold, sample with amplicon mean coverage less than ×500; Alt variant freq, variant allele frequency.
Table 3. Mutations in 25 samples obtained from 10 cases.
Table 3. Mutations in 25 samples obtained from 10 cases.
Year-No.StatusDetected
Mutation (n)
Mutations Passed
Filter with
Quality 100 * (n)
Mutations *
Having PP (n)
Gene (Variant)PPVAFTotal RDRecommended RD
-SC
2015-1
preserved morphology31233KDR(T > T/A)B66.814739
TP53(A > A/G)D40.412229
TP53(G > G/C)B24.457342
nuclear streaming44152KDR(T > A/A)B98.11602
TP53(A > A/G)D26.930937
necrosis28132KDR(T > T/A)B84.8334
TP53(A > A/G)D4346717
SC
2015-2
preserved morphology13121TP53(G > C/C)B100481
nuclear streaming14132TP53(A > C/C)D92.73553
TP53(G > C/C)B1001381
SC
2018-2
preserved morphology15154KIT(G > G/C)D47.8146815
PTEN(A > A/G)D52.524013
KRAS(G >G/T)D27.488037
TP53(G > C/C)D90.13333
nuclear streaming16144KIT(G > G/C)D47.8189515
PTEN(A > A/G)D50.231714
PIK3CA(A >A/G)D23.1161744
KRAS(G > G/T)D30.5118932
necrosis1521PTEN(A > A/G)D75125
SC
2018-3
preserved morphology30225ERBB4(C > C/A)D35.534627
KIT(A > A/C)B34.5204128
KDR(T > T/A)B63.629410
MET(A > A/G)B41.433618
TP53(G > C/C)D90.11542
necrosis26205ERBB4(C > C/A)D76.61845
KIT(A > A/C)B91699195
KDR(T > T/A)B90.12033
MET(A > A/G)B16.464475
TP53(G > C/C)D1001141
SC
2018-4
preserved morphology15142TP53(T > T/C)D687479
TP53(G > C/C)D1002171
nuclear streaming15142TP53(T > T/C)D91.51773
TP53(G > C/C)B96.9642
SC
2019-8
preserved morphology1050-----
nuclear streaming21173KIT(G > G/C)B48.5133315
TP53(T > T/C)D62.75910
TP53(G > C/C)B100121
SC
2019-9
preserved morphology27182TP53(G > C/C)B99.62332
STK11(C > C/G)B47.435215
nuclear streaming16142TP53(G > C/C)B99.53602
STK11(C > C/G)B1005061
DL
2018-14
preserved morphology16152PDGFR(C > C/G)B38.1120325
TP53(G > C/C)B99.88322
nuclear streaming17142PDGFR(C > C/G)B34.194728
TP53(G > C/C)B98.82572
necrosis16162PDGFR(C > C/G)B32.434530
TP53(G > C/C)B100551
DL
2018-19
preserved morphology1191TP53(G > G/C)B58.27912
nuclear streaming1191TP53(G > G/C)B40.817219
necrosis1050-----
DL
2020-2
preserved morphology12112KIT(A > A/C)B45.1136716
TP53(G > C/C)B99.66342
nuclear streaming14112KIT(A > A/C)B45.1220916
TP53(G > C/C)B99.86042
necrosis14112KIT(A > A/C)B46.7153716
TP53(G > C/C)B99.85662
* Mutations passed filter with quality 100, Recommended read depth was calculated using https://github.com/mvasinek/olgen-coverage-limit, accessed on 28 January 2022, where the variant has sufficient corresponding read depth; n, number; PP, PolyPhem; VAF, alt variant number; B, benign; D, probable damage, RD; read depth; -, not detected; preserved morphology, materials with well-preserved tumor morphology; nuclear streaming, materials with nuclear streaming; necrosis, materials with necrosis.
Table 4. Comparison of features between materials with well-preserved morphology, with nuclear streaming, and with necrosis.
Table 4. Comparison of features between materials with well-preserved morphology, with nuclear streaming, and with necrosis.
MetricnMeanSDSE95% CIvs.
With Well-Preserved
Morphology
p-Value
vs.
Nuclear Streaming
p-Value
vs.
Necrosis
p-Value
Percentage of Q30 bases
Well-preserved morphology1083.5175.0871.60879.8776–87.1564-0.9540.988
Nuclear streaming982.8484.8411.61379.1277–86.57000.954-0.998
Necrosis683.0416.7712.76475.9355–90.14780.9880.998-
Total PF reads
Well-preserved morphology10642.158145.218.28345.922.05353,8275.10–7406,040.9-0.7160.498
Nuclear streaming9584.832169.773.856,591.263454,332.31–715,331.690.716-0.910
Necrosis6548.916160.189.30265,397.009380,807.64–717,024.360.4980.910-
Percentage of on-target aligned reads
Well-preserved morphology1095.7023.0190.95493.5420–97.8620-0.9920.001
Nuclear streaming994.5966.4382.14689.6475–99.54590.992-0.02
Necrosis654.08640.55416.55611.5271–96.64620.0010.02-
Uniformity of coverage
Well-preserved morphology1082.4657.4912.36977.1056–87.8244-0.9600.237
Nuclear streaming983.4157.7052.56877.4922–89.33890.960-0.181
Necrosis675.2818.2413.36466.6328–83.93050.2370.181-
Variants with filter pass and quality 100 (%)
Well-preserved morphology1081.7015.5714.92470.56–92.84-0.9740.332
Nuclear streaming980.018.0006.00066.16–93.840.974-0.401
Necrosis660.8330.82512.68428.48–93.180.3320.401-
Variants without filter pass or with quality < 100 (%)
Well-preserved morphology1018.3015.5714.9247.16–29.44-0.9670.332
Nuclear streaming920.2217.8945.9656.47–33.980.967-0.408
Necrosis639.1730.82512.5846.82–71.520.3320.408-
Variant allele frequency (TP53)
Well-preserved morphology984.0122.6957.565366.566–101.457-0.6500.960
Nuclear streaming983.8028.6579.552461.772–105.8280.650-0.552.
Necrosis485.7728.46614.233440.403–130.9970.9600.552-
n, number of samples; PF, passing filter reads; SD, standard deviation; SE, standard error; CI, confidence interval; preserved morphology, materials with well-preserved tumor morphology; nuclear streaming, materials with nuclear streaming; necrosis, materials with necrosis; bold, statistically significant values.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Noda, Y.; Yamaka, R.; Atsumi, N.; Higasa, K.; Tsuta, K. Areas of Crush Nuclear Streaming Should Be Included as Tumor Content in the Era of Molecular Diagnostics. Cancers 2023, 15, 1910. https://doi.org/10.3390/cancers15061910

AMA Style

Noda Y, Yamaka R, Atsumi N, Higasa K, Tsuta K. Areas of Crush Nuclear Streaming Should Be Included as Tumor Content in the Era of Molecular Diagnostics. Cancers. 2023; 15(6):1910. https://doi.org/10.3390/cancers15061910

Chicago/Turabian Style

Noda, Yuri, Ryosuke Yamaka, Naho Atsumi, Koichiro Higasa, and Koji Tsuta. 2023. "Areas of Crush Nuclear Streaming Should Be Included as Tumor Content in the Era of Molecular Diagnostics" Cancers 15, no. 6: 1910. https://doi.org/10.3390/cancers15061910

APA Style

Noda, Y., Yamaka, R., Atsumi, N., Higasa, K., & Tsuta, K. (2023). Areas of Crush Nuclear Streaming Should Be Included as Tumor Content in the Era of Molecular Diagnostics. Cancers, 15(6), 1910. https://doi.org/10.3390/cancers15061910

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop