Papillary Thyroid Carcinoma (PTC) is characterized by unique tumor morphology, treatment response, and patient outcomes according to subtype and gender. While previous studies have implicated the intratumor bacterial microbiome in the incidence and progression of PTC, few studies have investigated the potential role of fungal and archaeal species in oncogenesis. In this study, we aimed to characterize the intratumor mycobiome and archaeometry in PTC with respect to its three primary subtypes: Classical (CPTC), Follicular Variant (FVPTC), and Tall Cell (TCPTC), and also with respect to gender. RNA-sequencing data were downloaded from The Cancer Genome Atlas (TCGA), including 453 primary tumor tissue samples and 54 adjacent solid tissue normal samples. The PathoScope 2.0 framework was used to extract fungal and archaeal microbial read counts from raw RNA-sequencing data. Overall, we found that the intratumor mycobiome and archaeometry share significant similarities in CPTC, FVPTC, and TCPTC, although most dysregulated species in CPTC are underabundant compared to normal. Furthermore, differences between the mycobiome and archaeometry were more significant between males and females, with a disproportionate number of fungal species overabundant in female tumor samples. Additionally, the expression of oncogenic PTC pathways was distinct across CPTC, FVPTC, and TCPTC, indicating that these microbes may uniquely contribute to PTC pathogenesis in each subtype. Furthermore, differences in the expression of these pathways were observed between males and females. Finally, we found a specific panel of fungi to be dysregulated in BRAF V600E-positive tumors. This study demonstrates the potential importance of microbial species to PTC incidence and oncogenesis.

Transfer-RNA-derived fragments (tRFs) are a novel class of small non-coding RNAs that have been implicated in oncogenesis. tRFs may act as post-transcriptional regulators by recruiting AGO proteins and binding to highly complementary regions of mRNA at seed regions, resulting in the knockdown of the transcript. Therefore, tRFs may be critical to tumorigenesis and warrant investigation as potential biomarkers. Meanwhile, the incidence of papillary thyroid carcinoma (PTC) has increased in recent decades and current diagnostic technology stands to benefit from new detection methods. Although small non-coding RNAs have been studied for their role in oncogenesis, there is currently no standard for their use as PTC biomarkers, and tRFs are especially underexplored. Accordingly, we aim to identify dysregulated tRFs in PTC that may serve as biomarker candidates. We identified dysregulated tRFs and driver genes between PTC primary tumor samples (n = 511) and adjacent normal tissue samples (n = 59). Expression data were obtained from MINTbase v2.0 and The Cancer Genome Atlas. Dysregulated tRFs and genes were analyzed in tandem to find pairs with anticorrelated expression. Significantly anticorrelated tRF-gene pairs were then tested for potential binding affinity using RNA22—if a heteroduplex can form via complementary binding, this would support the hypothesized RNA silencing mechanism. Four tRFs were significantly dysregulated in PTC tissue (p < 0.05), with only AsnGTT 3′-tRF being upregulated. Binding affinity analysis revealed that tRF-30-RY73W0K5KKOV (AsnGTT 3′-tRF) exhibits sufficient complementarity to potentially bind to and regulate transcripts of SLC26A4, SLC5A8, DIO2, and TPO, which were all found to be downregulated in PTC tissue. In the present study, we identified dysregulated tRFs in PTC and found that AsnGTT 3′-tRF is a potential post-transcriptional regulator and biomarker.