Washout DNA copy number analysis by low-coverage whole genome sequencing for assessment of thyroid FNAs

Abstract

Background: Papillary thyroid microcarcinoma (PTMC) is defined as a papillary carcinoma measuring ≤ 10 mm. The current management of PTMC has become more conservative; however, there are high-risk tumor features that can be revealed only postoperatively. For thyroid cancer, BRAF mutations and somatic copy number variation (CNV) are the most common genetic events. Molecular testing may contribute to clinical decision-making by molecular risk stratification, for example predicting lymph node (LN) metastasis. Here, we build a risk stratification model based on molecular profiling of thyroid fine needle aspiration (FNA) washout DNA (wDNA) for the differential diagnosis of thyroid nodules.

Methods: Fifty-eight patients were recruited, FNA wDNA samples were analyzed using CNV profiling through low-coverage whole genome sequencing (LC-WGS) and BRAF mutation was analyzed using quantitative PCR. FNA pathology was reported as a Bethesda System for Reporting Thyroid Cytopathology (BSRTC) score. Ultrasound examination produced a Thyroid Imaging Reporting and Data System (TIRADS) score.

Results: In total, 37 (63.8%) patients with a TIRADS score of 4A, 13 (22.4%) patients with a TIRADS score of 4B, and 8 (13.8%) patients with a TIRADS score of 4C were recruited after ultrasound examination. All patients underwent FNA with wDNA profiling. CNVs were identified in 17 (29.3%) patients. CNVs were frequent in patients with a BSRTC score of V or VI, including eight (47.1%) patients with a score of VI and five (29.4%) with a score of V, but not in patients with a score of III, II, or I (0%). BRAF mutation was not significantly correlated with BSRTC score. LN metastasis was found more frequently in CNV-positive (CNV+) than in CNV-negative (CNV-) patients (85.7% vs. 34.6%, odds ratio = 11.33, p = 0.002). In total, three molecular subtypes of thyroid nodules were identified in this study: 1) CNV+, 2) CNV- and BRAF positive (BRAF+), and 3) CNV- and BRAF negative (BRAF-). For the CNV+ subtype, 10 (83.3%) lesions with LN metastasis were found, including four (100%) small lesions (i.e. ≤ 5 mm). For the CNV- and BRAF+ nodules, LN metastases were detected in only seven (60.0%) larger tumors (i.e. > 5 mm). For CNV- and BRAF- tumors, LN metastasis was also frequently found in larger tumors only.

Conclusions: It is feasible to identify high-risk LN metastasis thyroid cancer from FNA washout samples preoperatively using wDNA CNV profiling using LC-WGS.

Keywords: LC-WGS; chromosome instability; genome sequencing; thyroid cancer; washout-DNA.