Network Pharmacological Analysis and Experimental Validation of the Effects of Silybin on Proliferation, Migration, and Immunotherapy of Papillary Thyroid Cancer

AIM

The study aimed to use network pharmacology research and in vitro experiments to investigate the material basis and molecular mechanisms of silybin in the treatment of papillary thyroid carcinoma.

BACKGROUND

Papillary thyroid cancer (PTC) has a decent prognosis; however, recurrence and metastasis are the leading causes of death in patients with PTC. A key research focus in thyroid cancer treatment is the inhibition of PTC proliferation, invasion, and migration. Silybin, the major active element in the traditional Chinese herb silymarin, has been used to treat a range of diseases, including cancer, but no study has been undertaken to determine whether it can help prevent PTC.

OBJECTIVE

In this study, we attempted to determine through network pharmacology and in vitro experiments if silybin inhibits the development of papillary thyroid cancer by inhibiting cell cycle and invasive migration.

METHODS

To predict the probable targets and underlying mechanisms of silybin against PTC, a network pharmacology research was performed. In vitro experiments were conducted to further evaluate silybin's anti-cancer properties and priority targets against PTC.

RESULTS

The datasets revealed a total of 489 silybin targets acting on PTC, and functional enrichment analysis suggested that the target genes were enriched in functions and pathways related to PTC development, invasion, migration, and immunotherapy. By constructing these target PPI networks, the seven hub genes, fibronectin 1 (FN1), tissue inhibitor of metalloproteinases 1 (TIMP1), N-cadherin (CDH2), collagen type III alpha 1 chain (COL3A1), cyclin D1 (CCND1), AP-1 transcription factor subunit (JUN), and hepatocyte growth factor receptor (MET) were found. These hub genes were determined to be highly linked to a worse clinicopathological form, a higher risk of metastatic recurrence, and a worse prognosis of PTC. The common immunological checkpoint gene expression levels were positively correlated with the expression levels of the hub genes. Silybin decreased the proliferative and metastatic capacity of PTC cells, according to in vitro investigations. When PTC was treated with silybin, the FN1/AKT signaling pathway was blocked, CCND1 expression was reduced, and CDH2, Vimentin, Snail, Slug and PD-L1 expressions were dramatically reduced, while E-cadherin expression was significantly elevated.

CONCLUSION

These findings provide preliminary evidence that silybin inhibits PTC cell proliferation, metastasis, and invasion by altering the FN1/AKT signaling pathway and inhibiting the EMT process. Silybin can reverse immunosuppression in papillary thyroid cancer by affecting immunological checkpoint gene expression levels. These studies provide a theoretical and experimental scientific basis for the potential anticancer effects of silybin on PTC.

Identification of prognostic biomarkers of breast cancer based on the immune-related gene module

Breast cancer (BC) is highly malignant and its mortality rate remains high. The development of immunotherapy has gradually improved the prognosis and survival rate of patients. Therefore, identifying molecular markers concerned with BC immunity is of great importance for the treatment of this disease. The Cancer Genome Atlas-breast invasive carcinoma (TCGA-BRCA) was utilized as the training set while the BC expression dataset from the gene expression omnibus database was taken as the validation set here. Weighted gene co-expression network analysis combined with Pearson analysis and Tumor immune estimation resource (TIMER) was used to obtain immune cell-related hub gene module. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed on this module. Then, receiver operating characteristic curves combining Kaplan-Meier was used to evaluate the effectiveness of the model. Feature genes were screened and the independence of risk score was evaluated by univariate and multivariate Cox analyses. Differences in immune characteristics were analyzed via single-sample gene set enrichment analysis and CIBERSORT, and differences in gene mutation frequency were assessed via GenVisR analysis. Finally, the expression levels of prognostic feature genes in BC cells were validated by quantitative reverse transcription polymerase chain reaction (qRT-PCR). In this study, cell immune-related gene modules in TCGA-BRCA were successfully excavated, and a five-gene (TNFRSF14, NFKBIA, DLG3, IRF2, and CYP27A1) prognostic model was established. The prognostic model could effectively forecast the prognosis and survival rate of BC patients. The result showed that human leukocyte antigen-related proteins and macrophage M2 scores were remarkably highly expressed in the high-risk group, whereas CD8+ T cells, natural killer cells, M1, and other anti-tumor cells were lowly expressed. The model could be used as an independent prognostic factor to predict the prognosis of BC patients. The results of qRT-PCR validation were consistent with the results in the database, that is, except DLG3, the other four feature genes were lowly expressed in BC. The five-gene model established in this study can predict the prognostic and immune mode of BC patients effectively, which is anticipated to become a feasible molecular target for BC therapy.

Deep small RNA-Seq reveals microRNAs expression profiles in lactating mammary gland of 2 sheep breeds with different milk performance

Micro ribonucleic acid (miRNA) is a type of noncoding RNA, and it has been revealed to play important roles in the activity of the mammary gland (MG) in some species. However, the function of miRNAs in MG of sheep is poorly understood. In the study, Gansu Alpine Merino (GAM; n = 9) and Small-tailed Han sheep (STH; n = 9) with different milk production traits were investigated. Microstructures and the expression profile of miRNAs of MG tissues at peak lactation were studied. Mature alveolar lumens of MG in appearance were larger in STH than GAM. The expression levels of CSN2 and the content of rough endoplasmic reticulum were also higher in STH ewes than GAM ewes. A total of 124 mature miRNAs were expressed, and 18 of these were differentially expressed between the 2 breeds. The KEGG analysis results showed that the targeted genes of differentially expressed miRNAs were mainly involved in some metabolic pathways and signaling pathways related to MG development, milk protein, and fat synthesis. The findings in the study can improve our understanding of the roles of miRNAs in the development and lactation of MG in sheep.

Tailoring spontaneous infrared emission of HgTe quantum dots with laser-printed plasmonic arrays

Chemically synthesized near-infrared to mid-infrared (IR) colloidal quantum dots (QDs) offer a promising platform for the realization of devices including emitters, detectors, security, and sensor systems. However, at longer wavelengths, the quantum yield of such QDs decreases as the radiative emission rate drops following Fermi's golden rule, while non-radiative recombination channels compete with light emission. Control over the radiative and non-radiative channels of the IR-emitting QDs is crucially important to improve the performance of IR-range devices. Here, we demonstrate strong enhancement of the spontaneous emission rate of near- to mid-IR HgTe QDs coupled to periodically arranged plasmonic nanoantennas, in the form of nanobumps, produced on the surface of glass-supported Au films via ablation-free direct femtosecond laser printing. The enhancement is achieved by simultaneous radiative coupling of the emission that spectrally matches the first-order lattice resonance of the arrays, as well as more efficient photoluminescence excitation provided by coupling of the pump radiation to the local surface plasmon resonances of the isolated nanoantennas. Moreover, coupling of the HgTe QDs to the lattice plasmons reduces the influence of non-radiative decay losses mediated by the formation of polarons formed between QD surface-trapped carriers and the IR absorption bands of dodecanethiol used as a ligand on the QDs, allowing us to improve the shape of the emission spectrum through a reduction in the spectral dip related to this ligand coupling. Considering the ease of the chemical synthesis and processing of the HgTe QDs combined with the scalability of the direct laser fabrication of nanoantennas with tailored plasmonic responses, our results provide an important step towards the design of IR-range devices for various applications.

Methods for Structural and Functional Analyses of Intramembrane Prenyltransferases in the UbiA Superfamily

The UbiA superfamily is a group of intramembrane prenyltransferases that generate lipophilic compounds essential in biological membranes. These compounds, which include various quinones, hemes, chlorophylls, and vitamin E, participate in electron transport and function as antioxidants, as well as acting as structural lipids of microbial cell walls and membranes. Prenyltransferases producing these compounds are involved in important physiological processes and human diseases. These UbiA superfamily members differ significantly in their enzymatic activities and substrate selectivities. This chapter describes examples of methods that can be used to group these intramembrane enzymes, analyze their activity, and screen and crystallize homolog proteins for structure determination. Recent structures of two archaeal homologs are compared with structures of soluble prenyltransferases to show distinct mechanisms used by the UbiA superfamily to control enzymatic activity in membranes.

Detection and quantitation of radiolabeled proteins in gels and blots

Permanent records of the results of electrophoretic separations of radiolabeled proteins and membrane-bound proteins (and RNA and DNA) can be made using autoradiography, fluorography, and phosphor imaging. These images can subsequently be quantified using densitometry to obtain a relative measure of the amount of radioactivity in a sample. This unit also includes descriptions of image-enhancement techniques and guidelines for selection of appropriate recording media.

NR4A1, 2, 3--an orphan nuclear hormone receptor family involved in cell apoptosis and carcinogenesis

Three members of the NR4A1/Nur77/ NGFIB orphan nuclear hormone receptor subfamily (NR4A1, NR4A2, and NR4A3) belong to the steroid nuclear hormone receptor superfamily. They share similar structural features and have no known natural ligand. They constitute immediate early genes that are induced by serum, growth factors and receptor engagement and are thus implicated in cell mitogenic responses. These nuclear receptors are transcription factors that exert their functions through activation and subsequent induction of the downstream pathways. They have been shown to play a role in complex pathways of cell survival and apoptosis. Although the expression of these genes have been shown to be pro-survival, it has also been reported that NR4A1 expression can cause apoptosis. These two opposite effects apparently result from distinct mechanisms: either transcriptional activation of genes responsible for cell survival or cell apoptosis, or translocation into the cytoplasm where they target the mitochondria and cause cell apoptosis via Bcl-2 binding. The latter mechanism constitutes a new paradigm of cellular apoptosis. In vitro functional studies using over-expression (gain of function) or gene inactivation (loss of function) type assays, combined with transgenic or knockout animal data in vivo, have revealed these effects and their physiological roles, including thymocyte development for NR4A1/3 and pro-survival in CNS for NR4A2. Recent studies have also suggested an important role of these receptors in cell transformation and tumorigenicity via both their anti-apoptotic and pro-apoptotic functions. In particular, the recent identification of a functional ligand for NR4A1 suggests that these members could potentially serve as drug targets for disease indications such as cancer. While many aspects of these receptors have been previously reviewed, this article focuses on new experimentation and discovery of their apoptotic and carcinogenic roles, and discusses their potential roles as therapeutic targets.

Detection and quantitation of radiolabeled proteins in gels and blots

This unit provides procedures for detecting an quantifying radiolabeled proteins in SDS-polyacrylamide gels or blots. Autoradiography requires fixed, dried gels before exposure to X-ray film. The results may be enhanced by using intensifying screens and/or preflashing the film. Alternatively fluorography with organic scintillant-impregnated gels may be used to enhance signal. The autoradiogram may be quantified by densitometry. Phosphorimaging is a useful alternative to film exposure.

Bcl-B, a novel Bcl-2 family member that differentially binds and regulates Bax and Bak

A novel human member of the Bcl-2 family was identified, Bcl-B, which is closest in amino acid sequence homology to the Boo (Diva) protein. The Bcl-B protein contains four Bcl-2 homology (BH) domains (BH1, BH2, BH3, BH4) and a predicted carboxyl-terminal transmembrane (TM) domain. The BCL-B mRNA is widely expressed in adult human tissues. The Bcl-B protein binds Bcl-2, Bcl-X(L), and Bax but not Bak. In transient transfection assays, Bcl-B suppresses apoptosis induced by Bax but not Bak. Deletion of the TM domain of Bcl-B impairs its association with intracellular organelles and diminishes its anti-apoptotic function. Bcl-B thus displays a unique pattern of selectivity for binding and regulating the function of other members of the Bcl-2 family.

Drosophila pro-apoptotic Bcl-2/Bax homologue reveals evolutionary conservation of cell death mechanisms

Genetic analysis of programmed cell death in Drosophila reveals many similarities with mammals. Heretofore, a missing link in the fly has been the absence of any Bcl-2/Bax family members, proteins that function in mammals as regulators of mitochondrial cytochrome c release. A Drosophila homologue of the human killer protein Bok (DBok) was identified. The predicted structure of DBok is similar to pore-forming Bcl-2/Bax family members. DBok induces apoptosis in insect and human cells, which is suppressible by anti-apoptotic human Bcl-2 family proteins. A caspase inhibitor suppressed DBok-induced apoptosis but did not prevent DBok-induced cell death. Moreover, DBok targets mitochondria and triggers cytochrome c release through a caspase-independent mechanism. These characteristics of DBok reveal evolutionary conservation of cell death mechanisms in flies and humans.

Epicubenol synthase. Origin of the oxygen atom of a bacterial sesquiterpene alcohol

Incubation of epicubenol synthase with farnesyl pyrophosphate in the presence of 11.1 atom% H2(18)O gave epicubenol (2) in which the hydroxyl oxygen atom was shown to be derived exclusively from water, as established by GC-selected ion monitoring MS of the derived TMS-epicubenol derivative (15).

The yeast retrotransposon Ty5 uses the anticodon stem-loop of the initiator methionine tRNA as a primer for reverse transcription

Retrotransposons and retroviruses replicate by reverse transcription of an mRNA intermediate. Most retroelements initiate reverse transcription from a host-encoded tRNA primer. DNA synthesis typically extends from the 3'-OH of the acceptor stem, which is complementary to sequences on the retroelement mRNA (the primer binding site, PBS). However, for some retrotransposons, including the yeast Ty5 elements, sequences in the anticodon stem-loop of the initiator methionine tRNA (IMT) are complementary to the PBS. We took advantage of the genetic tractability of the yeast system to investigate the mechanism of Ty5 priming. We found that transposition frequencies decreased at least 800-fold for mutations in the Ty5 PBS that disrupt complementarity with the IMT. Similarly, transposition was reduced at least 200-fold for IMT mutations in the anticodon stem-loop. Base pairing between the Ty5 PBS and IMT is essential for transposition, as compensatory changes that restored base pairing between the two mutant RNAs restored transposition significantly. An analysis of 12 imt mutants with base changes outside of the region of complementarity failed to identify other tRNA residues important for transposition. In addition, assays carried out with heterologous IMTs from Schizosaccharomyces pombe and Arabidopsis thaliana indicated that residues outside of the anticodon stem-loop have at most a fivefold effect on transposition. Our genetic system should make it possible to further define the components required for priming and to understand the mechanism by which Ty5's novel primer is generated.

cDNA of the yeast retrotransposon Ty5 preferentially recombines with substrates in silent chromatin

The yeast retrotransposon Ty5 preferentially integrates into regions of silent chromatin. Ty5 cDNA also recombines with homologous sequences, generating tandem elements or elements that have exchanged markers between cDNA and substrate. In this study, we demonstrate that Ty5 integration depends upon the conserved DD(35)E domain of integrase and cis-acting sequences at the end of the long terminal repeat (LTR) implicated in integrase binding. cDNA recombination requires Rad52p, which is responsible for homologous recombination. Interestingly, Ty5 cDNA recombines at least three times more frequently with substrates in silent chromatin than with a control substrate at an internal chromosomal locus. This preference depends upon the Ty5 targeting domain that is responsible for integration specificity, suggesting that localization of cDNA to silent chromatin results in the enhanced recombination. Recombination with a telomeric substrate occasionally generates highly reiterated Ty5 arrays, and mechanisms for tandem element formation were explored by using a plasmid-based recombination assay. Point mutations were introduced into plasmid targets, and recombination products were characterized to determine recombination initiation sites. Despite our previous observation of the importance of the LTR in forming tandem elements, recombination cannot simply be explained by crossover events between the LTRs of substrate and cDNA. We propose an alternative model based on single-strand annealing, where single-stranded cDNA initiates tandem element formation and the LTR is required for strand displacement to form a looped intermediate. Retrotransposons are increasingly found associated with chromosome ends, and amplification of Ty5 by both integration and recombination exemplifies how retroelements can contribute to telomere dynamics.

The pheromone response pathway activates transcription of Ty5 retrotransposons located within silent chromatin of Saccharomyces cerevisiae

The Saccharomyces retrotransposon Ty5 integrates preferentially into transcriptionally inactive regions (silent chromatin) at the HM loci and telomeres. We found that silent chromatin represses basal Ty5 transcription, indicating that these elements are encompassed by silent chromatin in their native genomic context. Because transcription is a requirement for transposition, integration into silent chromatin would appear to prevent subsequent rounds of replication. Using plasmid-borne Ty5-lacZ constructs, we found that Ty5 expression is haploid specific and is repressed 10-fold in diploid strains. Ty5 transcription is also regulated by the pheromone response pathway and is induced approximately 20-fold upon pheromone treatment. Deletion analysis of the Ty5 LTR promoter revealed that a 33 bp region with three perfect matches to the pheromone response element is responsible for both mating pheromone and cell-type regulation. Transcriptional repression of Ty5 by silent chromatin can be reversed by pheromone treatment, which leads to transcription and transposition. Ty5 replication, therefore, is normally repressed by silent chromatin and appears to be induced during mating. This is the first example of transcriptional activation of a gene that naturally resides within silent chromatin.

High frequency cDNA recombination of the saccharomyces retrotransposon Ty5: The LTR mediates formation of tandem elements

Retroelement cDNA can integrate into the genome using the element-encoded integrase or it can recombine with preexisting elements using the recombination system of the host. Recombination is a particularly important pathway for the yeast retrotransposon Ty5 and accounts for approximately 30% of the putative transposition events when a homologous substrate is carried on a plasmid and approximately 7% when the substrate is located at the chromosomal URA3 locus. Characterization of recombinants revealed that they are either simple replacements of the marker gene tandem elements. Using an assay system in which the donor element and recombination substrates are separated, we found that the long terminal repeats (LTRs) are critical for tandem element formation. LTR-containing substrates generate tandem elements at frequencies more than 10-fold higher than similarly sized internal Ty5 sequences. Internal sequences, however, facilitate tandem element formation when associated with an LTR, and there is a linear relationship between frequencies of tandem element formation and the length of LTR-containing substrates. We propose that recombination is initiated between the LTRs of the cDNA and substrate and that internal sequences promote tandem element formation by facilitating sequence alignment. Because of its location in subtelomeric regions, recombinational amplification of Ty5 may contribute to the organizations of chromosome ends.

The Saccharomyces retrotransposon Ty5 integrates preferentially into regions of silent chromatin at the telomeres and mating loci

The nonrandom integration of retrotransposons and retroviruses suggests that chromatin influences target choice. Targeted integration, in turn, likely affects genome organization. In Saccharomyces, native Ty5 retrotransposons are located near telomeres and the silent mating locus HMR. To determine whether this distribution is a consequence of targeted integration, we isolated a transposition-competent Ty5 element from S. paradoxus, a species closely related to S. cerevisiae. This Ty5 element was used to develop a transposition assay in S. cerevisiae to investigate target preference of de novo transposition events. Of 87 independent Ty5 insertions, approximately 30% were located on chromosome III, indicating this small chromosome (approximately 1/40 of the yeast genome) is a highly preferred target. Mapping of the exact location of 19 chromosome III insertions showed that 18 were within or adjacent to transcriptional silencers flanking HML and HMR or the type X subtelomeric repeat. We predict Ty5 target preference is attributable to interactions between transposition intermediates and constituents of silent chromatin assembled at these sites. Ty5 target preference extends the link between telomere structure and reverse transcription as carried out by telomerase and Drosophila retrotransposons.

Multiple non-LTR retrotransposons in the genome of Arabidopsis thaliana

DNA sequence analysis near the Arabidopsis thaliana ABI3 gene revealed the presence of a non-LTR retrotransposon insertion that we have designated Ta11-1. This insertion is 6.2 kb in length and encodes two overlapping reading frames with similarity to non-LTR retrotransposon proteins, including reverse transcriptase. A polymerase chain reaction assay was developed based on conserved amino acid sequences shared between the Ta11-1 reverse transcriptase and those of non-LTR retrotransposons from other species. Seventeen additional A. thaliana reverse transcriptases were identified that range in nucleotide similarity from 48-88% (Ta12-Ta28). Phylogenetic analyses indicated that the A. thaliana sequences are more closely related to each other than to elements from other organisms, consistent with the vertical evolution of these sequences over most of their evolutionary history. One sequence, Ta17, is located in the mitochondrial genome. The remaining are nuclear and of low copy number among 17 diverse A. thaliana ecotypes tested, suggesting that they are not highly active in transposition. The paucity of retrotransposons and the small genome size of A. thaliana support the hypothesis that most repetitive sequences have been lost from the genome and that mechanisms may exist to prevent amplification of extant element families.

Destabilization of rbcS sense transcripts by antisense RNA

Steady-state rbcS mRNA levels are drastically reduced in transgenic tobacco plants that express rbcS antisense RNAs. We have found that these reductions are not due to an effect of the antisense RNA at the level of rbcS transcription; rather, the sense mRNAs are more actively degraded in the mutant than wild-type plants. We have examined the kinetics of this turnover process by inhibiting transcription with cordycepin, and have found that rbcS sense mRNA decay is accelerated about five-fold in the antisense plants. This provides direct evidence that antisense RNAs can serve to destabilize sense transcripts in plants.