Comprehensive Analysis of the Transcriptional and Mutational Landscape of Follicular and Papillary Thyroid Cancers.

Follicular thyroid carcinoma (FTC) and benign follicular adenoma (FA) are indistinguishable by preoperative diagnosis due to their similar histological features. Here we report the first RNA sequencing study of these tumors, with data for 30 minimally invasive FTCs (miFTCs) and 25 FAs. We also compared 77 classical papillary thyroid carcinomas (cPTCs) and 48 follicular variant of PTCs (FVPTCs) to observe the differences in their molecular properties. Mutations in H/K/NRAS, DICER1, EIF1AX, IDH1, PTEN, SOS1, and SPOP were identified in miFTC or FA. We identified a low frequency of fusion genes in miFTC (only one, PAX8-PPARG), but a high frequency of that in PTC (17.60%). The frequencies of BRAFV600E and H/K/NRAS mutations were substantially different in miFTC and cPTC, and those of FVPTC were intermediate between miFTC and cPTC. Gene expression analysis demonstrated three molecular subtypes regardless of their histological features, including Non-BRAF-Non-RAS (NBNR), as well as BRAF-like and RAS-like. The novel molecular subtype, NBNR, was associated with DICER1, EIF1AX, IDH1, PTEN, SOS1, SPOP, and PAX8-PPARG. The transcriptome of miFTC or encapsulated FVPTC was indistinguishable from that of FA, providing a molecular explanation for the similarly indolent behavior of these tumors. We identified upregulation of genes that are related to mitochondrial biogenesis including ESRRA and PPARGC1A in oncocytic follicular thyroid neoplasm. Arm-level copy number variations were correlated to histological and molecular characteristics. These results expanded the current molecular understanding of thyroid cancer and may lead to new diagnostic and therapeutic approaches to the disease.

Development of a common platform for the noninvasive prenatal diagnosis of X-linked diseases.

OBJECTIVE: The aim of this study was to develop a common targeted massively parallel sequencing platform for the noninvasive prenatal diagnosis (NIPD) of multiple X-linked diseases. METHOD: The custom capture probe was designed to target 33 genes and recombination hotspots. We tested the carrier mother and male proband pair of 6 families. Plasma DNA of the pregnant carrier mother was collected at different gestational weeks and sequenced. The fetal genotype of each family was determined by estimating the imbalance between the 2 maternal haplotypes constructed using a common custom-designed platform. RESULTS: The targeted sequencing of the maternal, proband, and fetal genomic DNAs and maternal plasma DNAs resulted in uniform coverage across the target region. Three to 5 recombination points were observed in each sample. However, these recombination points did not affect the haplotype dosage analysis for fetal genotype prediction. Consequently, all fetal genotypes in the 6 families obtained from haplotype dosage analysis of maternal plasma sequencing data were predicted correctly. CONCLUSIONS: Since a single platform that covers multiple diseases may prevent the need for disease-specific probes for the NIPD of individual disorders, this approach may provide a practical advantage for clinically implementing the NIPD of X-linked diseases.

Noninvasive prenatal diagnosis of duchenne muscular dystrophy: comprehensive genetic diagnosis in carrier, proband, and fetus.

BACKGROUND: Noninvasive prenatal diagnosis of monogenic disorders using maternal plasma and targeted massively parallel sequencing is being investigated actively. We previously demonstrated that comprehensive genetic diagnosis of a Duchenne muscular dystrophy (DMD) patient is feasible using a single targeted sequencing platform. Here we demonstrate the applicability of this approach to carrier detection and noninvasive prenatal diagnosis. METHODS: Custom solution-based target enrichment was designed to cover the entire dystrophin (DMD) gene region. Targeted massively parallel sequencing was performed using genomic DNA from 4 mother and proband pairs to test whether carrier status could be detected reliably. Maternal plasma DNA at varying gestational weeks was collected from the same families and sequenced using the same targeted platform to predict the inheritance of the DMD mutation by their fetus. Overrepresentation of an inherited allele was determined by comparing the allele fraction of 2 phased haplotypes after examining and correcting for the recombination event. RESULTS: The carrier status of deletion/duplication and point mutations was detected reliably through using a single targeted massively parallel sequencing platform. Whether the fetus had inherited the DMD mutation was predicted correctly in all 4 families as early as 6 weeks and 5 days of gestation. In one of these, detection of the recombination event and reconstruction of the phased haplotype produced a correct diagnosis. CONCLUSIONS: Noninvasive prenatal diagnosis of DMD is feasible using a single targeted massively parallel sequencing platform with tiling design.

Image-based deep learning model using DNA methylation data predicts the origin of cancer of unknown primary.

Cancer of unknown primary (CUP) is a rare type of metastatic cancer in which the origin of the tumor is unknown. Since the treatment strategy for patients with metastatic tumors depends on knowing the primary site, accurate identification of the origin site is important. Here, we developed an image-based deep-learning model that utilizes a vision transformer algorithm for predicting the origin of CUP. Using DNA methylation dataset of 8,233 primary tumors from The Cancer Genome Atlas (TCGA), we categorized 29 cancer types into 18 organ classes and extracted 2,312 differentially methylated CpG sites (DMCs) from non-squamous cancer group and 420 DMCs from squamous cell cancer group. Using these DMCs, we created organ-specific DNA methylation images and used them for model training and testing. Model performance was evaluated using 394 metastatic cancer samples from TCGA (TCGA-meta) and 995 samples (693 primary and 302 metastatic cancers) obtained from 20 independent external studies. We identified that the DNA methylation image reveals a distinct pattern based on the origin of cancer. Our model achieved an overall accuracy of 96.95 % in the TCGA-meta dataset. In the external validation datasets, our classifier achieved overall accuracies of 96.39 % and 94.37 % in primary and metastatic tumors, respectively. Especially, the overall accuracies for both primary and metastatic samples of non-squamous cell cancer were exceptionally high, with 96.79 % and 96.85 %, respectively.

LORIS robustly predicts patient outcomes with immune checkpoint blockade therapy using common clinical, pathologic and genomic features.

Despite the revolutionary impact of immune checkpoint blockade (ICB) in cancer treatment, accurately predicting patient responses remains challenging. Here, we analyzed a large dataset of 2,881 ICB-treated and 841 non-ICB-treated patients across 18 solid tumor types, encompassing a wide range of clinical, pathologic and genomic features. We developed a clinical score called LORIS (logistic regression-based immunotherapy-response score) using a six-feature logistic regression model. LORIS outperforms previous signatures in predicting ICB response and identifying responsive patients even with low tumor mutational burden or programmed cell death 1 ligand 1 expression. LORIS consistently predicts patient objective response and short-term and long-term survival across most cancer types. Moreover, LORIS showcases a near-monotonic relationship with ICB response probability and patient survival, enabling precise patient stratification. As an accurate, interpretable method using a few readily measurable features, LORIS may help improve clinical decision-making in precision medicine to maximize patient benefit. LORIS is available as an online tool at https://loris.ccr.cancer.gov/ .

An immunogenetic basis for lung cancer risk.

Cancer risk is influenced by inherited mutations, DNA replication errors, and environmental factors. However, the influence of genetic variation in immunosurveillance on cancer risk is not well understood. Leveraging population-level data from the UK Biobank and FinnGen, we show that heterozygosity at the human leukocyte antigen (HLA)-II loci is associated with reduced lung cancer risk in smokers. Fine-mapping implicated amino acid heterozygosity in the HLA-II peptide binding groove in reduced lung cancer risk, and single-cell analyses showed that smoking drives enrichment of proinflammatory lung macrophages and HLA-II+ epithelial cells. In lung cancer, widespread loss of HLA-II heterozygosity (LOH) favored loss of alleles with larger neopeptide repertoires. Thus, our findings nominate genetic variation in immunosurveillance as a critical risk factor for lung cancer.

Central neurocytoma exhibits radial glial cell signatures with FGFR3 hypomethylation and overexpression.

We explored the genomic events underlying central neurocytoma (CN), a rare neoplasm of the central nervous system, via multiomics approaches, including whole-exome sequencing, bulk and single-nuclei RNA sequencing, and methylation sequencing. We identified FGFR3 hypomethylation leading to FGFR3 overexpression as a major event in the ontogeny of CN that affects crucial downstream events, such as aberrant PI3K-AKT activity and neuronal development pathways. Furthermore, we found similarities between CN and radial glial cells based on analyses of gene markers and CN tumor cells and postulate that CN tumorigenesis is due to dysregulation of radial glial cell differentiation into neurons. Our data demonstrate the potential role of FGFR3 as one of the leading drivers of tumorigenesis in CN.

Validation of a Machine Learning Model to Predict Immunotherapy Response in Head and Neck Squamous Cell Carcinoma.

Head and neck squamous-cell carcinoma (HNSCC) is a disease with a generally poor prognosis; half of treated patients eventually develop recurrent and/or metastatic (R/M) disease. Patients with R/M HNSCC generally have incurable disease with a median survival of 10 to 15 months. Although immune-checkpoint blockade (ICB) has improved outcomes in patients with R/M HNSCC, identifying patients who are likely to benefit from ICB remains a challenge. Biomarkers in current clinical use include tumor mutational burden and immunohistochemistry for programmed death-ligand 1, both of which have only modest predictive power. Machine learning (ML) has the potential to aid in clinical decision-making as an approach to estimate a tumor's likelihood of response or a patient's likelihood of experiencing clinical benefit from therapies such as ICB. Previously, we described a random forest ML model that had value in predicting ICB response using 11 or 16 clinical, laboratory, and genomic features in a pan-cancer development cohort. However, its applicability to certain cancer types, such as HNSCC, has been unknown, due to a lack of cancer-type-specific validation. Here, we present the first validation of a random forest ML tool to predict the likelihood of ICB response in patients with R/M HNSCC. The tool had adequate predictive power for tumor response (area under the receiver operating characteristic curve = 0.65) and was able to stratify patients by overall (HR = 0.53 [95% CI 0.29-0.99], p = 0.045) and progression-free (HR = 0.49 [95% CI 0.27-0.87], p = 0.016) survival. The overall accuracy was 0.72. Our study validates an ML predictor in HNSCC, demonstrating promising performance in a novel cohort of patients. Further studies are needed to validate the generalizability of this algorithm in larger patient samples from additional multi-institutional contexts.

A whole-genome reference panel of 14,393 individuals for East Asian populations accelerates discovery of rare functional variants.

Underrepresentation of non-European (EUR) populations hinders growth of global precision medicine. Resources such as imputation reference panels that match the study population are necessary to find low-frequency variants with substantial effects. We created a reference panel consisting of 14,393 whole-genome sequences including more than 11,000 Asian individuals. Genome-wide association studies were conducted using the reference panel and a population-specific genotype array of 72,298 subjects for eight phenotypes. This panel yields improved imputation accuracy of rare and low-frequency variants within East Asian populations compared with the largest reference panel. Thirty-nine previously unidentified associations were found, and more than half of the variants were East Asian specific. We discovered genes with rare protein-altering variants, including LTBP1 for height and GPR75 for body mass index, as well as putative regulatory mechanisms for rare noncoding variants with cell type-specific effects. We suggest that this dataset will add to the potential value of Asian precision medicine.

Pre-treatment serum albumin and mutational burden as biomarkers of response to immune checkpoint blockade.

The effects of cytokine and protein stabilizing carriers, such as serum albumin, on tumor response to immune checkpoint blockade (ICB) is not well understood. By examining 1714 patients across 16 cancer types, we found that high pretreatment serum albumin level predicts favorable tumor radiographic response following ICB treatment in a dose-dependent fashion. Serum albumin is a candidate biomarker that can be combined with tumor mutational burden (TMB) for additional predictive capacity, and the tumor response rate to ICB was ~49% in the albumin-high/TMB-high group.

Analysis of surface characteristics of (Y, Nb)-TZP after finishing and polishing.

PURPOSE: This in vitro study aimed to evaluate the surface characteristics of a full veneer crown fabricated chairside (CS) from a (Y, Nb)-TZP zirconia block in response to conventional zirconia grinding and polishing. MATERIALS AND METHODS: Zirconia crowns (n = 40) were first prepared and divided into two groups of materials: Labside (LS) and CS, after which each specimen went through a five-step grinding and polishing procedure. Following each surface treatment, surface characteristics were analyzed using confocal laser microscopy (CLSM), average surface roughness (Ra) values were processed from the profile data through Gaussian filtering, and X-ray diffraction pattern analysis was performed to evaluate the monoclinic (M) phase content. Then, a representative specimen was selected for field-emission scanning electron microscopy (FE-SEM), followed by a final analysis of the roughness and X-ray diffraction of the specimens using the independent t-test and repeated measures analysis of variance (RM-ANOVA). RESULTS: In every group, polishing significantly reduced the Ra values (P < .001). There was no significant difference in Ra between the polished state CS and LS. Furthermore, CLSM and FE-SEM investigations revealed that even though grain exposure was visible in CS specimens throughout the as-delivered and ground states, the exposure was reduced after polishing. Moreover, while no phase transformation was visible in the LS, phase transformation was visible in CS after every surface treatment, with the M phase content of the CS group showing a significant reduction after polishing (P < .001). CONCLUSION: Within the limits of this study, clinically acceptable level of surface finishing of (Y, Nb)-TZP can be achieved after conventional zirconia polishing sequence.

Improved prediction of immune checkpoint blockade efficacy across multiple cancer types.

Only a fraction of patients with cancer respond to immune checkpoint blockade (ICB) treatment, but current decision-making procedures have limited accuracy. In this study, we developed a machine learning model to predict ICB response by integrating genomic, molecular, demographic and clinical data from a comprehensively curated cohort (MSK-IMPACT) with 1,479 patients treated with ICB across 16 different cancer types. In a retrospective analysis, the model achieved high sensitivity and specificity in predicting clinical response to immunotherapy and predicted both overall survival and progression-free survival in the test data across different cancer types. Our model significantly outperformed predictions based on tumor mutational burden, which was recently approved by the U.S. Food and Drug Administration for this purpose1. Additionally, the model provides quantitative assessments of the model features that are most salient for the predictions. We anticipate that this approach will substantially improve clinical decision-making in immunotherapy and inform future interventions.

Indigenous Ancestry and Admixture in the Uruguayan Population.

The Amerindian group known as the Charrúas inhabited Uruguay at the timing of European colonial contact. Even though they were extinguished as an ethnic group as a result of a genocide, Charrúan heritage is part of the Uruguayan identity both culturally and genetically. While mitochondrial DNA studies have shown evidence of Amerindian ancestry in living Uruguayans, here we undertake whole-genome sequencing of 10 Uruguayan individuals with self-declared Charruan heritage. We detect chromosomal segments of Amerindian ancestry supporting the presence of indigenous genetic ancestry in living descendants. Specific haplotypes were found to be enriched in "Charrúas" and rare in the rest of the Amerindian groups studied. Some of these we interpret as the result of positive selection, as we identified selection signatures and they were located mostly within genes related to the infectivity of specific viruses. Historical records describe contacts of the Charrúas with other Amerindians, such as Guaraní, and patterns of genomic similarity observed here concur with genomic similarity between these groups. Less expected, we found a high genomic similarity of the Charrúas to Diaguita from Argentinian and Chile, which could be explained by geographically proximity. Finally, by fitting admixture models of Amerindian and European ancestry for the Uruguayan population, we were able to estimate the timing of the first pulse of admixture between European and Uruguayan indigenous peoples in approximately 1658 and the second migration pulse in 1683. Both dates roughly concurring with the Franciscan missions in 1662 and the foundation of the city of Colonia in 1680 by the Spanish.

Recent Improvements in Genomic and Transcriptomic Understanding of Anaplastic and Poorly Differentiated Thyroid Cancers.

Anaplastic thyroid cancer (ATC) is a lethal human cancer with a 5-year survival rate of less than 10%. Recently, its genomic and transcriptomic characteristics have been extensively elucidated over 5 years owing to advance in high throughput sequencing. These efforts have extended molecular understandings into the progression mechanisms and therapeutic vulnerabilities of aggressive thyroid cancers. In this review, we provide an overview of genomic and transcriptomic alterations in ATC and poorly-differentiated thyroid cancer, which are distinguished from differentiated thyroid cancers. Clinically relevant genomic alterations and deregulated signaling pathways will be able to shed light on more effective prevention and stratified therapeutic interventions for affected patients.

Genomic and Transcriptomic Characteristics According to Size of Papillary Thyroid Microcarcinoma.

It is controversial as to whether papillary thyroid microcarcinoma (PTMC) has some genomic and transcriptomic characteristics that differentiate between an early-stage lesion that would eventually evolve into the larger papillary thyroid cancer (PTC), and an occult indolent cancer in itself. To investigate this, we comprehensively elucidated the genomic and transcriptomic landscapes of PTMCs of different sizes, using a large-scaled database. This study included 3435 PTCs, 1985 of which were PTMCs. We performed targeted next-generation sequencing for 221 PTCs and integrated these data with the data including The Cancer Genome Atlas (TCGA) project. The frequency of v-raf murine sarcoma viral oncogene homolog B (BRAF)V600E mutation was higher in PTMCs >0.5 cm than that in very small PTMCs (≤0.5 cm) and decreased again in PTCs >2 cm. Among PTMCs, the prevalence of mutations in rat sarcoma (RAS) and telomerase reverse transcriptase (TERT) promoter was not significantly different according to their size, but lower than in large PTCs. There was no change in the tumor mutational burden, the number of driver mutations, and transcriptomic profiles with tumor size, among PTMCs and all PTCs. Although a few genes with differential expression and TERT promoter mutations were found in a few PTMCs, our findings showed that there were no useful genomic or transcriptomic characteristics for the prediction of the future progression of PTMC.

BRAFV600E Transduction of an SV40-Immortalized Normal Human Thyroid Cell Line Induces Dedifferentiated Thyroid Carcinogenesis in a Mouse Xenograft Model.

Background: Despite active studies of the clinical importance of BRAFV600E, suitable research models to investigate the role of this mutation in the etiopathogenesis of human thyroid cancers are limited. Thus, we generated cell lines by transducing the simian virus (SV)-40 immortalized human thyroid cell line Nthy-ori 3-1 (Nthy) with lentiviral vectors expressing either BRAFWT (Nthy/WT) or BRAFV600E. Nthy/WT and Nthy/V600E cells were then xenografted into mice to evaluate the carcinogenic role of BRAFV600E. Methods: Each cell line was subcutaneously injected into NOD.Cg-Prkdcscid Il2rgtm1Wjl/SzJ mice, and a pathological analysis was performed. The effects of the mutation were further verified by using a BRAFV600E-selective inhibitor (PLX-4032, vemurafenib). The transcriptome was analyzed by RNA sequencing and compared with data from The Cancer Cell Line Encyclopedia and Gene Expression Omnibus. Results: While Nthy/WT was not tumorigenic in vivo, Nthy/V600E formed tumors reaching 2784.343 mm3 in 4 weeks, on average. A pathological analysis indicated that Nthy/V600E tumors were dedifferentiated thyroid cancer. We found metastases in the lung, liver, and relevant lymph nodes. A transcriptomic analysis revealed 5512 differentially expressed genes (DEGs) between the mutant and wild-type cell lines, and more DEGs were shared with anaplastic thyroid cancer than with papillary thyroid cancer. BRAFV600E activated the cell cycle mainly by regulating G1/S phases. PLX-4032 treatment significantly inhibited tumor growth and metastasis. Conclusions: Our data show that BRAFV600E plays a pivotal role in the carcinogenic transformation of an SV40-transfected immortalized normal human thyroid cell line. This xenograft model is expected to contribute to studies of the etiopathogenesis and treatment of highly malignant thyroid cancers.

CXCL16 positively correlated with M2-macrophage infiltration, enhanced angiogenesis, and poor prognosis in thyroid cancer.

Although various chemokines have pro-tumorigenic actions in cancers, the effects of CXCL16 remain controversial. The aim of this study was to investigate the molecular characteristics of CXCL16-expressing papillary thyroid cancers (PTCs). CXCL16 expressions were significantly higher in PTCs than benign or normal thyroid tissues. In the TCGA dataset for PTCs, a higher CXCL16 expression was associated with M2 macrophage- and angiogenesis-related genes and poor prognostic factors including a higher TNM staging and the BRAFV600E mutation. PTCs with a higher expression of 3-gene panel including CXCL16, AHNAK2, and THBS2 showed poor recurrence-free survivals than that of the lower expression group. Next, shCXCL16 was introduced into BHP10-3SCp cells to deplete the endogenous CXCL16, and then, the cells were subcutaneously injected to athymic mice. Tumors from the BHP10-3SCpshCXCL16 exhibited a delayed tumor growth with decreased numbers of ERG+ endothelial cells and F4/80+ macrophages than those from the BHP10-3SCpcontrol. CXCL16-related genes including AHNAK2 and THBS2 were downregulated in the tumors from the BHP10-3SCpshCXCL16 compared with that from the BHP10-3SCpcontrol. In conclusion, a higher CXCL16 expression was associated with macrophage- and angiogenesis-related genes and aggressive phenotypes in PTC. Targeting CXCL16 may be a good therapeutic strategy for advanced thyroid cancer.

NARD: whole-genome reference panel of 1779 Northeast Asians improves imputation accuracy of rare and low-frequency variants.

Here, we present the Northeast Asian Reference Database (NARD), including whole-genome sequencing data of 1779 individuals from Korea, Mongolia, Japan, China, and Hong Kong. NARD provides the genetic diversity of Korean (n = 850) and Mongolian (n = 384) ancestries that were not present in the 1000 Genomes Project Phase 3 (1KGP3). We combined and re-phased the genotypes from NARD and 1KGP3 to construct a union set of haplotypes. This approach established a robust imputation reference panel for Northeast Asians, which yields the greatest imputation accuracy of rare and low-frequency variants compared with the existing panels. NARD imputation panel is available at https://nard.macrogen.com/ .

Interaction of BRAF-induced ETS factors with mutant TERT promoter in papillary thyroid cancer.

Synergistic effects of BRAFV600E and TERT promoter mutations on the poor clinical outcomes in papillary thyroid cancer (PTC) have been demonstrated. The potential mechanism of this phenomenon has been proposed: MAPK pathway activation by the BRAFV600E mutation may upregulate E-twenty six (ETS) transcription factors, increasing TERT expression by binding to the ETS-binding site generated by the TERT promoter mutation; however, it has not yet been fully proven. This article provides transcriptomic insights into the interaction between BRAFV600E and TERT promoter mutations mediated by ETS factors in PTC. RNA sequencing data on 266 PTCs from The Cancer Genome Atlas and 65 PTCs from our institute were analyzed for gene expression changes and related molecular pathways, and the results of transcriptomic analyses were validated by in vitro experiments. TERT mRNA expression was increased by the coexistence of BRAFV600E and TERT promoter mutations (fold change, 16.17; q-value = 7.35 × 10-12 vs no mutation). In the ETS family of transcription factors, ETV1, ETV4 and ETV5 were upregulated by the BRAFV600E/MAPK pathway activation. These BRAFV600E-induced ETS factors selectively bound to the mutant TERT promoter. The molecular pathways activated by BRAFV600E were further augmented by adding the TERT promoter mutation, and the pathways related to immune responses or adhesion molecules were upregulated by TERT expression. The mechanism of the synergistic effect between BRAFV600E and TERT promoter mutations on cancer invasiveness and progression in PTC may be explained by increased TERT expression, which may result from the BRAF-induced upregulation of several ETS transcription factors.

Integrative analysis of genomic and transcriptomic characteristics associated with progression of aggressive thyroid cancer.

Anaplastic thyroid cancer (ATC) and advanced differentiated thyroid cancers (DTCs) show fatal outcomes, unlike DTCs. Here, we demonstrate mutational landscape of 27 ATCs and 86 advanced DTCs by massively-parallel DNA sequencing, and transcriptome of 13 ATCs and 12 advanced DTCs were profiled by RNA sequencing. TERT, AKT1, PIK3CA, and EIF1AX were frequently co-mutated with driver genes (BRAFV600E and RAS) in advanced DTCs as well as ATC, but tumor suppressors (e.g., TP53 and CDKN2A) were predominantly altered in ATC. CDKN2A loss was significantly associated with poor disease-specific survival in patients with ATC or advanced DTCs, and up-regulation of CD274 (PD-L1) and PDCD1LG2 (PD-L2). Transcriptome analysis revealed a fourth molecular subtype of thyroid cancer (TC), ATC-like, which hardly reflects the molecular signatures in DTC. Furthermore, the activation of JAK-STAT signaling pathway could be a potential druggable target in RAS-positive ATC. Our findings provide insights for precision medicine in patients with advanced TCs.