A statistical method for image-mediated association studies discovers genes and pathways associated with four brain disorders.

Brain imaging and genomics are critical tools enabling characterization of the genetic basis of brain disorders. However, imaging large cohorts is expensive and may be unavailable for legacy datasets used for genome-wide association studies (GWASs). Using an integrated feature selection/aggregation model, we developed an image-mediated association study (IMAS), which utilizes borrowed imaging/genomics data to conduct association mapping in legacy GWAS cohorts. By leveraging the UK Biobank image-derived phenotypes (IDPs), the IMAS discovered genetic bases underlying four neuropsychiatric disorders and verified them by analyzing annotations, pathways, and expression quantitative trait loci (eQTLs). A cerebellar-mediated mechanism was identified to be common to the four disorders. Simulations show that, if the goal is identifying genetic risk, our IMAS is more powerful than a hypothetical protocol in which the imaging results were available in the GWAS dataset. This implies the feasibility of reanalyzing legacy GWAS datasets without conducting additional imaging, yielding cost savings for integrated analysis of genetics and imaging.

cLD: Rare-variant linkage disequilibrium between genomic regions identifies novel genomic interactions.

Linkage disequilibrium (LD) is a fundamental concept in genetics; critical for studying genetic associations and molecular evolution. However, LD measurements are only reliable for common genetic variants, leaving low-frequency variants unanalyzed. In this work, we introduce cumulative LD (cLD), a stable statistic that captures the rare-variant LD between genetic regions, which reflects more biological interactions between variants, in addition to lack of recombination. We derived the theoretical variance of cLD using delta methods to demonstrate its higher stability than LD for rare variants. This property is also verified by bootstrapped simulations using real data. In application, we find cLD reveals an increased genetic association between genes in 3D chromatin interactions, a phenomenon recently reported negatively by calculating standard LD between common variants. Additionally, we show that cLD is higher between gene pairs reported in interaction databases, identifies unreported protein-protein interactions, and reveals interacting genes distinguishing case/control samples in association studies.

Interaction-integrated linear mixed model reveals 3D-genetic basis underlying Autism.

Genetic interactions play critical roles in genotype-phenotype associations. We developed a novel interaction-integrated linear mixed model (ILMM) that integrates a priori knowledge into linear mixed models. ILMM enables statistical integration of genetic interactions upfront and overcomes the problems of searching for combinations. To demonstrate its utility, with 3D genomic interactions (assessed by Hi-C experiments) as a priori, we applied ILMM to whole-genome sequencing data for Autism Spectrum Disorders (ASD) and brain transcriptome data, revealing the 3D-genetic basis of ASD and 3D-expression quantitative loci (3D-eQTLs) for brain tissues. Notably, we reported a potential mechanism involving distal regulation between FOXP2 and DNMT3A, conferring the risk of ASD.

Reconstruction of Microbial Haplotypes by Integration of Statistical and Physical Linkage in Scaffolding.

DNA sequencing technologies provide unprecedented opportunities to analyze within-host evolution of microorganism populations. Often, within-host populations are analyzed via pooled sequencing of the population, which contains multiple individuals or "haplotypes." However, current next-generation sequencing instruments, in conjunction with single-molecule barcoded linked-reads, cannot distinguish long haplotypes directly. Computational reconstruction of haplotypes from pooled sequencing has been attempted in virology, bacterial genomics, metagenomics, and human genetics, using algorithms based on either cross-host genetic sharing or within-host genomic reads. Here, we describe PoolHapX, a flexible computational approach that integrates information from both genetic sharing and genomic sequencing. We demonstrated that PoolHapX outperforms state-of-the-art tools tailored to specific organismal systems, and is robust to within-host evolution. Importantly, together with barcoded linked-reads, PoolHapX can infer whole-chromosome-scale haplotypes from 50 pools each containing 12 different haplotypes. By analyzing real data, we uncovered dynamic variations in the evolutionary processes of within-patient HIV populations previously unobserved in single position-based analysis.

Construction of a novel cuproptosis-related gene signature for predicting prognosis and estimating tumor immune microenvironment status in papillary thyroid carcinoma.

BACKGROUND: Cuproptosis, a new form of programmed cell death, has been recently reported to be closely related to tumor progression. However, the significance of cuproptosis-related genes (CRGs) in papillary thyroid carcinoma (PTC) is still unclear. Therefore, this study aimed to investigate the role of the CRG signature in prognosis prediction and immunotherapeutic effect estimation in patients with PTC. METHODS: RNA-seq data and the corresponding clinical information of patients with PTC were obtained from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Comprehensive analyses, namely, consensus clustering, immune analyses, functional enrichment, least absolute shrinkage and selection operator-multivariate Cox regression, and nomogram analysis, were performed to identify new molecular subgroups, determine the tumor immune microenvironment (TIME) status of the identified subgroups, and construct a clinical model. Independent verification cohort data and quantitative real-time polymerase chain reaction (qPCR) was performed to validate the expression of specific prognosis-related and differentially expressed CRGs (P-DECRGs). RESULTS: In the TCGA database, 476 patients with PTC who had complete clinical and follow-up information were included. Among 135 CRGs, 21 were identified as P-DECRGs. Two molecular subgroups with significantly different disease-free survival and TIME statuses were identified based on these 21 P-DECRGs. The differentially expressed genes between the two subgroups were mainly associated with immune regulation. The risk model and nomogram were constructed based on four specific P-DECRGs and validated as accurate prognostic predictions and TIME status estimation for PTC by TCGA and GEO verification cohorts. Finally, the qPCR results of 20 PTC and paracancerous thyroid tissues validated those in the TCGA database. CONCLUSIONS: Four specific P-DECRGs in PTC were identified, and a clinical model based on them was established, which may be helpful for individualized immunotherapeutic strategies and prognostic prediction in patients with PTC.

Exosomes Derived from Human Pulmonary Artery Endothelial Cells Shift the Balance between Proliferation and Apoptosis of Smooth Muscle Cells.

BACKGROUND: The overproliferation of pulmonary vascular cells is noted in pulmonary hypertension. The role of exosomes from pulmonary artery endothelial cells (PAEC) in the proliferation and apoptosis of pulmonary artery smooth muscle cells (PASMC) remains unclear. METHODS: Exosomes were isolated and purified from the culture medium of PAEC using a commercial kit. Lipopolysaccharide (LPS), hypoxia, and hydrogen peroxide were utilized to induce PAEC injury. Coculture of PAEC and PASMC was conducted using Transwell plates, and GW4869 was applied to inhibit exosome release. The proliferation and apoptosis level of PASMC was assayed by MTT assay, apoptosis staining, and cleaved caspase-3 immunoblotting. Plasma exosomes were isolated by differential ultracentrifugation. RESULTS: LPS or hypoxia enhance exosome release from PAEC. Release of PAEC-derived exosomes positively correlates with LPS concentration. The coculture of LPS-disposed PAEC with PASMC leads to overproliferation and apoptosis resistance in PASMC, and the exosome inhibitor GW4869 can partly cancel out this effect. Exosomes derived from PAEC could be internalized into PASMC, and thus promote proliferation and induce apoptosis resistance in PASMC. Idiopathic pulmonary arterial hypertension patients exhibit a higher circulation level of endothelium-derived exosomes. CONCLUSIONS: Inflammation and hypoxia could induce PAEC to release exosomes. PAEC- derived exosomes are involved in overproliferation and apoptosis resistance in PASMC, by which they may contribute to the pathogenesis of pulmonary hypertension.

rvTWAS: identifying gene-trait association using sequences by utilizing transcriptome-directed feature selection.

Toward the identification of genetic basis of complex traits, transcriptome-wide association study (TWAS) is successful in integrating transcriptome data. However, TWAS is only applicable for common variants, excluding rare variants in exome or whole-genome sequences. This is partly because of the inherent limitation of TWAS protocols that rely on predicting gene expressions. Our previous research has revealed the insight into TWAS: the 2 steps in TWAS, building and applying the expression prediction models, are essentially genetic feature selection and aggregations that do not have to involve predictions. Based on this insight disentangling TWAS, rare variants' inability of predicting expression traits is no longer an obstacle. Herein, we developed "rare variant TWAS," or rvTWAS, that first uses a Bayesian model to conduct expression-directed feature selection and then uses a kernel machine to carry out feature aggregation, forming a model leveraging expressions for association mapping including rare variants. We demonstrated the performance of rvTWAS by thorough simulations and real data analysis in 3 psychiatric disorders, namely schizophrenia, bipolar disorder, and autism spectrum disorder. We confirmed that rvTWAS outperforms existing TWAS protocols and revealed additional genes underlying psychiatric disorders. Particularly, we formed a hypothetical mechanism in which zinc finger genes impact all 3 disorders through transcriptional regulations. rvTWAS will open a door for sequence-based association mappings integrating gene expressions.

Large-scale alternative polyadenylation-wide association studies to identify putative cancer susceptibility genes.

Alternative polyadenylation (APA) modulates mRNA processing in the 3' untranslated regions (3' UTR), affecting mRNA stability and translation efficiency. Research into genetically regulated APA has the potential to provide insights into cancer risk. Herein, we conducted large alternative polyadenylation-wide association studies (APA-WAS) to investigate associations of APA levels with cancer risk. Genetic models were built to predict APA levels in multiple tissues using genotype and RNA-sequencing data from 1,337 samples from the Genotype-Tissue Expression Project. Associations of genetically predicted APA levels with cancer risk were assessed by applying the prediction models to data from large genome-wide association studies of six common cancers among European-ancestry populations, including breast, ovary, prostate, colorectum, lung, and pancreas. A total of 58 risk genes (corresponding to 76 APA sites) were associated with at least one type of cancer, including 25 genes previously not linked to cancer susceptibility. Of the identified risk APAs, 97.4% and 26.3% were supported by 3' UTR APA quantitative trait loci and co-localization analyses, respectively. Luciferase reporter assays for four selected putative regulatory 3' UTR variants demonstrated that the risk alleles of 3' UTR variants, rs324015 (STAT6), rs2280503 (DIP2B), rs1128450 (FBXO38), and rs145220637 (LDHA), significantly increased the post-transcriptional activities of their target genes compared to reference alleles. Furthermore, knockdown of the target genes confirmed their ability to promote proliferation and migration. Overall, this study provides insights into the role of APA in the genetic susceptibility to common cancers.

Enhancing Disease Risk Gene Discovery by Integrating Transcription Factor-Linked Trans-located Variants into Transcriptome-Wide Association Analyses.

Transcriptome-wide association studies (TWAS) have been successful in identifying putative disease susceptibility genes by integrating gene expression predictions with genome-wide association studies (GWAS) data. However, current TWAS models only consider cis-located variants to predict gene expression. Here, we introduce transTF-TWAS, which includes transcription factor (TF)-linked trans-located variants for model building. Using data from the Genotype-Tissue Expression project, we predict alternative splicing and gene expression and applied these models to large GWAS datasets for breast, prostate, and lung cancers. Our analysis revealed 887 putative cancer susceptibility genes, including 465 in regions not yet reported by previous GWAS and 137 in known GWAS loci but not yet reported previously, at Bonferroni-corrected P < 0.05. We demonstrate that transTF-TWAS surpasses other approaches in both building gene prediction models and identifying disease-associated genes. These results have shed new light on several genetically driven key regulators and their associated regulatory networks underlying disease susceptibility.

Large-scale alternative polyadenylation (APA)-wide association studies to identify putative susceptibility genes in human common cancers.

Alternative polyadenylation (APA) modulates mRNA processing in the 3' untranslated regions (3'UTR), which affect mRNA stability and translation efficiency. Here, we build genetic models to predict APA levels in multiple tissues using sequencing data of 1,337 samples from the Genotype-Tissue Expression, and apply these models to assess associations between genetically predicted APA levels and cancer risk with data from large genome-wide association studies of six common cancers, including breast, ovary, prostate, colorectum, lung, and pancreas among European-ancestry populations. At a Bonferroni-corrected P □<□0.05, we identify 58 risk genes, including seven in newly identified loci. Using luciferase reporter assays, we demonstrate that risk alleles of 3'UTR variants, rs324015 ( STAT6 ), rs2280503 ( DIP2B ), rs1128450 ( FBXO38 ) and rs145220637 ( LDAH ), could significantly increase post-transcriptional activities of their target genes compared to reference alleles. Further gene knockdown experiments confirm their oncogenic roles. Our study provides additional insight into the genetic susceptibility of these common cancers.

Identification and Validation of Prognostic Related Hallmark ATP-Binding Cassette Transporters Associated With Immune Cell Infiltration Patterns in Thyroid Carcinoma.

ATP-binding cassette (ABC) transporters are a large superfamily of membrane proteins that facilitate the translocation of heterogeneous substrates. Studies indicate that ABC transporters may play important roles in various carcinomas. However, the correlation between ABC transporters and immunomodulation in thyroid carcinoma (TC), as well as the prognoses for this disease, is poorly understood.TC data from The Cancer Genome Atlas (TCGA) database were used to identify prognostic hallmark ABC transporters associated with immune cell infiltration patterns via multiple bioinformatic analyses. Thereafter, quantitative real-time polymerase chain reaction (qRT-PCR) was performed to validate the expression of these selected hallmark ABC transporters in both TC and para-cancerous thyroid tissues. Of a total of 49 ABC transporters, five (ABCA8, ABCA12, ABCB6, ABCB8, and ABCC10) were identified as hallmark ABC transporters. All five were differentially expressed in TC and associated with the relapse-free survival rates of patients with TC. Immunoregulation by these five hallmark ABC transporters involved the modulation of various aspects of immune cell infiltration, such as hot or cold tumor subsets and the abundances of infiltrating immune cells, as well as specific immunomodulators and chemokines. Besides the diverse significantly correlated factors, the five hallmark ABC transporters and correlated genes were most highly enriched in plasma membrane, transporter activity, and transmembrane transport of small molecules. In addition, many chemicals, namely bisphenol A and vincristine, affected the expression of these five transporters. The qRT-PCR results of collected TC and para-cancerous thyroid tissues were consistent with those of TCGA. The findings in this study may reveal the role played by these five hallmark ABC transporters in regulating immune cell infiltration patterns in TC as well as the molecular mechanisms underlying their functions, leading to a better understanding of their potential prognostic and immunotherapeutic values.

Tripartite motif containing 63, regulated by E26 transformation specific variant 4, facilitates the thyroid carcinoma progression and the AKT, p38, and ERK signaling pathways.

The current study aimed to investigate the role of tripartite motif containing 63 (TRIM63) in the progression of thyroid carcinomas. The results showed that TRIM63 was highly expressed in thyroid carcinomas tissues. TRIM63 knockdown inhibited the proliferation and induced the apoptosis of thyroid carcinoma cells, and overexpression of TRIM63 promoted the cell proliferation ability. These results were further confirmed by the in vivo growth of xenograft tumors. Subsequently, the underlying mechanism was explored. TRIM63 silencing repressed the AKT, p38, and ERK signaling pathways in thyroid carcinoma cells, and the contrary results were observed in TRIM63-upregulated thyroid carcinoma cells. Furthermore, we found that E26 transformation specific variant 4 (ETV4) regulated the transcription of TRIM63. The loss of TRIM63 reversed the ETV4 overexpression-induced promotion of proliferation in thyroid carcinomas cells. In conclusion, TRIM63, regulated by ETV4, activates the AKT, p38, and ERK signaling pathways and facilitates the thyroid carcinoma development.

Integrating transcription factor occupancy with transcriptome-wide association analysis identifies susceptibility genes in human cancers.

Transcriptome-wide association studies (TWAS) have successfully discovered many putative disease susceptibility genes. However, TWAS may suffer from inaccuracy of gene expression predictions due to inclusion of non-regulatory variants. By integrating prior knowledge of susceptible transcription factor occupied elements, we develop sTF-TWAS and demonstrate that it outperforms existing TWAS approaches in both simulation and real data analyses. Under the sTF-TWAS framework, we build genetic models to predict alternative splicing and gene expression in normal breast, prostate and lung tissues from the Genotype-Tissue Expression project and apply these models to data from large genome-wide association studies (GWAS) conducted among European-ancestry populations. At Bonferroni-corrected P < 0.05, we identify 354 putative susceptibility genes for these cancers, including 189 previously unreported in GWAS loci and 45 in loci unreported by GWAS. These findings provide additional insight into the genetic susceptibility of human cancers. Additionally, we show the generalizability of the sTF-TWAS on non-cancer diseases.

Disentangling genetic feature selection and aggregation in transcriptome-wide association studies.

The success of transcriptome-wide association studies (TWAS) has led to substantial research toward improving the predictive accuracy of its core component of genetically regulated expression (GReX). GReX links expression information with genotype and phenotype by playing two roles simultaneously: it acts as both the outcome of the genotype-based predictive models (for predicting expressions) and the linear combination of genotypes (as the predicted expressions) for association tests. From the perspective of machine learning (considering SNPs as features), these are actually two separable steps-feature selection and feature aggregation-which can be independently conducted. In this study, we show that the single approach of GReX limits the adaptability of TWAS methodology and practice. By conducting simulations and real data analysis, we demonstrate that disentangled protocols adapting straightforward approaches for feature selection (e.g., simple marker test) and aggregation (e.g., kernel machines) outperform the standard TWAS protocols that rely on GReX. Our development provides more powerful novel tools for conducting TWAS. More importantly, our characterization of the exact nature of TWAS suggests that, instead of questionably binding two distinct steps into the same statistical form (GReX), methodological research focusing on optimal combinations of feature selection and aggregation approaches will bring higher power to TWAS protocols.

Tumor vasculature-targeting nanomedicines.

Uncontrolled tumor growth and subsequent distant metastasis are highly dependent on an adequate nutrient supply from tumor blood vessels, which have relatively different pathophysiological characteristics from those of normal vasculature. Obviously, strategies targeting tumor vasculature, such as anti-angiogenic drugs and vascular disrupting agents, are attractive methods for cancer therapy. However, the off-target effects and high dose administration of these drug regimens critically restrict their clinical applications. In recent years, nanomedicines focused on tumor vasculature have been shown to be superior to traditional therapeutic methods and do not induce side effects. This review will first highlight the recent development of tumor vasculature-targeting nanomedicines from the following four aspects: 1) angiogenesis-inhibiting nanomedicines (AINs); 2) vasculature-disrupting nanomedicines (VDNs); 3) vasculature infarction nanomedicines (VINs); and 4) vasculature-regulating nanomedicines (VRNs). Furthermore, the design principles, limitations, and future directions are also discussed. STATEMENT OF SIGNIFICANCE: Based on the essential roles of tumor blood vessels, the therapeutic strategies targeting tumor vasculature have exhibited good clinical therapeutic outcomes. However, poor patient adherence to free drug administration limits their clinical usage. Nanomedicines have great potential to overcome the abovementioned obstacle. This review summarizes the tumor-vasculature targeting nanomedicines from four aspects: 1) angiogenesis-inhibiting nanomedicines (AINs); 2) vasculature-disrupting nanomedicines (VDNs); 3) vasculature infarction nanomedicines (VINs); and 4) vasculature regulating nanomedicines (VRNs). In addition, this review provides perspectives on this research field.

Neurophysiologic Profiling of At-Risk Low and Very Low Birth-Weight Infants Using Magnetic Resonance Imaging.

Low birth-weight (LBW) and very low birth-weight (VLBW) newborns have increased risks of brain injuries, growth failure, motor difficulties, developmental coordination disorders or delay, and adult-onset vascular diseases. However, relatively little is known of the neurobiologic underpinnings. To clarify the pathophysiologic vulnerabilities of such neonates, we applied several advanced techniques for assessing brain physiology, namely T2-relaxation-under-spin-tagging (TRUST) magnetic resonance imaging (MRI) and phase-contrast (PC) MRI. This enabled quantification of oxygen extraction fraction (OEF), global cerebral blood flow (CBF), and cerebral metabolic rate of oxygen (CMRO2). A total of 50 neonates (LBW-VLBW, 41; term controls, 9) participated in this study. LBW-VLBW neonates were further stratified as those with (LBW-VLBW-a, 24) and without (LBW-VLBW-n, 17) structural MRI (sMRI) abnormalities. TRUST and PC MRI studies were undertaken to determine OEF, CBF, and CMRO2. Ultimately, CMRO2 proved significantly lower (p = 0.01) in LBW-VLBW (vs term) neonates, both LBW-VLBW-a and LBW-VLBW-n subsets showing significantly greater physiologic deficits than term controls (p = 0.03 and p = 0.04, respectively). CMRO2 and CBF in LBW-VLBW-a and LBW-VLBW-n subsets did not differ significantly (p > 0.05), although OEF showed a tendency to diverge (p = 0.15). However, OEF values in the LBW-VLBW-n subset differed significantly from those of term controls (p = 0.02). Compared with brain volume or body weight, these physiologic parameters yield higher area-under-the-curve (AUC) values for distinguishing neonates of the LBW-VLBW-a subset. The latter displayed distinct cerebral metabolic and hemodynamic, whereas changes were marginal in the LBW-VLBW-n subset (i.e., higher OEF and lower CBF and CMRO2) by comparison. Physiologic imaging may therefore be useful in identifying LBW-VLBW newborns at high risk of irreversible brain damage.

Functions and Molecular Mechanisms of Deltex Family Ubiquitin E3 Ligases in Development and Disease.

Ubiquitination is a posttranslational modification of proteins that significantly affects protein stability and function. The specificity of substrate recognition is determined by ubiquitin E3 ligase during ubiquitination. Human Deltex (DTX) protein family, which functions as ubiquitin E3 ligases, comprises five members, namely, DTX1, DTX2, DTX3, DTX3L, and DTX4. The characteristics and functional diversity of the DTX family proteins have attracted significant attention over the last decade. DTX proteins have several physiological and pathological roles and are closely associated with cell signal transduction, growth, differentiation, and apoptosis, as well as the occurrence and development of various tumors. Although they have been extensively studied in various species, data on structural features, biological functions, and potential mechanisms of action of the DTX family proteins remain limited. In this review, recent research progress on each member of the DTX family is summarized, providing insights into future research directions and potential strategies in disease diagnosis and therapy.

Fruquintinib beneficial in elderly patient with neoplastic pericardial effusion from rectal cancer: A case report.

BACKGROUND: Neoplastic pericardial effusion (NPE) is a rare consequence of rectal cancer and carries a poor prognosis. Optimal management has yet to be determined. Fruquintinib is an oral anti-vascular endothelial growth factor receptor tyrosine kinase inhibitor approved by the China Food and Drug Administration in September 2018 as third-line treatment of metastatic colorectal cancer. CASE SUMMARY: Herein, we report an elderly patient with NPE from rectal cancer who responded to the use of fruquintinib. In March 2015, a 65-year-old Chinese woman diagnosed with KRAS-mutated adenocarcinoma of the rectum was subjected to proctectomy, adjuvant concurrent chemoradiotherapy, and adjuvant chemotherapy. By October 2018, a mediastinal mass was detected via computed tomography. The growth had invaded parietal pericardium and left hilum, displaying features of rectal adenocarcinoma in a bronchial biopsy. FOLFIRI and FOLFOX chemotherapeutic regimens were administered as first- and second-line treatments. After two cycles of second-line agents, a sizeable pericardial effusion resulting in tamponade was drained by pericardial puncture. Fluid cytology showed cells consistent with rectal adenocarcinoma. Single-agent fruquintinib was initiated on January 3, 2019, as a third-line therapeutic. Ten cycles were delivered before the NPE recurred and other lesions progressed. The recurrence-free interval for NPE was 9.2 mo, attesting to the efficacy of fruquintinib. Ultimately, the patient entered a palliative care unit for best supportive care. CONCLUSION: Fruquintinib may confer good survival benefit in elderly patients with NPEs due to rectal cancer.

Circular RNA hsa_circ_0001666 sponges miR­330­5p, miR­193a­5p and miR­326, and promotes papillary thyroid carcinoma progression via upregulation of ETV4.

Circular RNAs (circRNAs) are a group of regulators that affect the aggressive behaviors of several types of cancer. Hsa_circ_0001666 (also referred to as hsa_circ_000742) is a newly discovered circRNA that is upregulated in human papillary thyroid carcinoma (PTC) based on microarray analysis. However, the role of hsa_circ_0001666 in PTC progression remains unknown. Thus, the aim of the present study was to determine the potential function and underlying mechanism of hsa_circ_0001666 in PTC. The results demonstrated that hsa_circ_0001666 was upregulated in both PTC clinical samples and cell lines. Its expression was associated with lymph node metastasis of patients with PTC. Knocking down hsa_circ_0001666 expression inhibited cell proliferation, as evidenced by decreased cell viability, arrest of cell cycle progression at the G1 phase and an increase in cell cycle­associated proteins. Apoptosis rates and expression levels of pro­apoptotic proteins were also increased by silencing hsa_circ_0001666. In xenograft experiments, the oncogenic effect of hsa_circ_0001666 on tumor growth was verified. Additionally, luciferase reporter assays showed that hsa_circ_0001666 and ETS variant transcription factor 4 (ETV4) shared common binding sites with three microRNAs [(miRNA/miR)­330­5p, miR­193a­5p and miR­326]. Knockdown of these miRNAs separately reversed the inhibitory effect of hsa_circ_0001666 small interfering RNAs on PTC tumor aggressiveness, and ETV4 overexpression also induced a similar effect to that of miRNA inhibitors. Thus, hsa_circ_0001666 may function as an oncogene, promoting PTC tumorigenesis via the miR­330­5p/miR­193a­5p/miR­326/ETV4 pathway. This provides a basis for identifying potential novel therapeutic targets for PTC.

The Knockdown of ETV4 Inhibits the Papillary Thyroid Cancer Development by Promoting Ferroptosis Upon SLC7A11 Downregulation.

Papillary thyroid cancer (PTC) is one of the most prevalent endocrine malignancies. Herein, we aimed to provide a new viewpoint for the PTC progression and explore a new target for the effective therapy for PTC. We found that E26 transformation specific (ETS) variant 4 (ETV4, an ETS family transcription factor) was upregulated in PTC tissues and cells. In vitro experiments exhibited that silencing ETV4 suppressed PTC cell proliferation and cell cycle progression, while the overexpression of ETV4 gained the opposite results. Dual-luciferase reporter assay highlighted that ETV4 could upregulate the solute carrier family 7 member 11 (SLC7A11, a key role for cysteine uptake in ferroptosis) transcription by binding to its promoter region directly. Moreover, the viability inhibition of PTC cells induced by the knockdown of ETV4 was at least partly through the promotion of ferroptosis upon the downregulation of SLC7A11. In in vivo experiment, the results showed that the downregulation of ETV4 repressed the tumor development through the low expression of SLC7A11, and the ETV4 overexpression obtained the contrary effects. Overall, the data suggested that the knockdown of ETV4 suppressed the PTC progression by promoting ferroptosis upon SLC7A11 downregulation.