CAP-RNAseq: an integrated pipeline for functional annotation and prioritization of co-expression clusters.

Cluster analysis is one of the most widely used exploratory methods for visualization and grouping of gene expression patterns across multiple samples or treatment groups. Although several existing online tools can annotate clusters with functional terms, there is no all-in-one webserver to effectively prioritize genes/clusters using gene essentiality as well as congruency of mRNA-protein expression. Hence, we developed CAP-RNAseq that makes possible (1) upload and clustering of bulk RNA-seq data followed by identification, annotation and network visualization of all or selected clusters; and (2) prioritization using DepMap gene essentiality and/or dependency scores as well as the degree of correlation between mRNA and protein levels of genes within an expression cluster. In addition, CAP-RNAseq has an integrated primer design tool for the prioritized genes. Herein, we showed using comparisons with the existing tools and multiple case studies that CAP-RNAseq can uniquely aid in the discovery of co-expression clusters enriched with essential genes and prioritization of novel biomarker genes that exhibit high correlations between their mRNA and protein expression levels. CAP-RNAseq is applicable to RNA-seq data from different contexts including cancer and available at http://konulabapps.bilkent.edu.tr:3838/CAPRNAseq/ and the docker image is downloadable from https://hub.docker.com/r/konulab/caprnaseq.

Long-Term Acetylcholinesterase Depletion Alters the Levels of Key Synaptic Proteins while Maintaining Neuronal Markers in the Aging Zebrafish (Danio rerio) Brain.

INTRODUCTION: Interventions targeting cholinergic neurotransmission like acetylcholinesterase (AChE) inhibition distinguish potential mechanisms to delay age-related impairments and attenuate deficits related to neurodegenerative diseases. However, the chronic effects of these interventions are not well described. METHODS: In the current study, global levels of cholinergic, cellular, synaptic, and inflammation-mediating proteins were assessed within the context of aging and chronic reduction of AChE activity. Long-term depletion of AChE activity was induced by using a mutant zebrafish line, and they were compared with the wildtype group at young and old ages. RESULTS: Results demonstrated that AChE activity was lower in both young and old mutants, and this decrease coincided with a reduction in ACh content. Additionally, an overall age-related reduction in AChE activity and the AChE/ACh ratio was observed, and this decline was more prominent in wildtype groups. The levels of an immature neuronal marker were upregulated in mutants, while a glial marker showed an overall reduction. Mutants had preserved levels of inhibitory and presynaptic elements with aging, whereas glutamate receptor subunit levels declined. CONCLUSION: Long-term AChE activity depletion induces synaptic and cellular alterations. These data provide further insights into molecular targets and adaptive responses following the long-term reduction of AChE activity that was also targeted pharmacologically to treat neurodegenerative diseases in human subjects.

Functional analysis of co-expression networks of zebrafish ace2 reveals enrichment of pathways associated with development and disease.

Human Angiotensin I Converting Enzyme 2 (ACE2) plays an essential role in blood pressure regulation and SARS-CoV-2 entry. ACE2 has a highly conserved, one-to-one ortholog (ace2) in zebrafish, which is an important model for human diseases. However, the zebrafish ace2 expression profile has not yet been studied during early development, between genders, across different genotypes, or in disease. Moreover, a network-based meta-analysis for the extraction of functionally enriched pathways associated with differential ace2 expression is lacking in the literature. Herein, we first identified significant development-, tissue-, genotype-, and gender-specific modulations in ace2 expression via meta-analysis of zebrafish Affymetrix transcriptomics datasets (ndatasets = 107); and the correlation analysis of ace2 meta-differential expression profile revealed distinct positively and negatively correlated local functionally enriched gene networks. Moreover, we demonstrated that ace2 expression was significantly modulated under different physiological and pathological conditions related to development, tissue, gender, diet, infection, and inflammation using additional RNA-seq datasets. Our findings implicate a novel translational role for zebrafish ace2 in organ differentiation and pathologies observed in the intestines and liver.

CoVrimer: A tool for aligning SARS-CoV-2 primer sequences and selection of conserved/degenerate primers.

As mutations in SARS-CoV-2 virus accumulate rapidly, novel primers that amplify this virus sensitively and specifically are in demand. We have developed a webserver named CoVrimer by which users can search for and align existing or newly designed conserved/degenerate primer pair sequences against the viral genome and assess the mutation load of both primers and amplicons. CoVrimer uses mutation data obtained from an online platform established by NGDC-CNCB (12 May 2021) to identify genomic regions, either conserved or with low levels of mutations, from which potential primer pairs are designed and provided to the user for filtering based on generalized and SARS-CoV-2 specific parameters. Alignments of primers and probes can be visualized with respect to the reference genome, indicating variant details and the level of conservation. Consequently, CoVrimer is likely to help researchers with the challenges posed by viral evolution and is freely available at http://konulabapps.bilkent.edu.tr:3838/CoVrimer/.

Increased SGK1 activity potentiates mineralocorticoid/NaCl-induced kidney injury.

Serum and glucocorticoid-regulated kinase 1 (SGK1) stimulates aldosterone-dependent renal Na+ reabsorption and modulates blood pressure. In addition, genetic ablation or pharmacological inhibition of SGK1 limits the development of kidney inflammation and fibrosis in response to excess mineralocorticoid signaling. In this work, we tested the hypothesis that a systemic increase in SGK1 activity would potentiate mineralocorticoid/salt-induced hypertension and kidney injury. To that end, we used a transgenic mouse model with increased SGK1 activity. Mineralocorticoid/salt-induced hypertension and kidney damage was induced by unilateral nephrectomy and treatment with deoxycorticosterone acetate and NaCl in the drinking water for 6 wk. Our results show that although SGK1 activation did not induce significantly higher blood pressure, it produced a mild increase in glomerular filtration rate, increased albuminuria, and exacerbated glomerular hypertrophy and fibrosis. Transcriptomic analysis showed that extracellular matrix- and immune response-related terms were enriched in the downregulated and upregulated genes, respectively, in transgenic mice. In conclusion, we propose that systemically increased SGK1 activity is a risk factor for the development of mineralocorticoid-dependent kidney injury in the context of low renal mass and independently of blood pressure.NEW & NOTEWORTHY Increased activity of the protein kinase serum and glucocorticoid-regulated kinase 1 may be a risk factor for accelerated renal damage. Serum and glucocorticoid-regulated kinase 1 expression could be a marker for the rapid progression toward chronic kidney disease and a potential therapeutic target to slow down the process.

Advanced preclinical models for evaluation of drug-induced liver injury - consensus statement by the European Drug-Induced Liver Injury Network [PRO-EURO-DILI-NET].

Drug-induced liver injury (DILI) is a major cause of acute liver failure (ALF) and one of the leading indications for liver transplantation in Western societies. Given the wide use of both prescribed and over the counter drugs, DILI has become a major health issue for which there is a pressing need to find novel and effective therapies. Although significant progress has been made in understanding the molecular mechanisms underlying DILI, our incomplete knowledge of its pathogenesis and inability to predict DILI is largely due to both discordance between human and animal DILI in preclinical drug development and a lack of models that faithfully recapitulate complex pathophysiological features of human DILI. This is exemplified by the hepatotoxicity of acetaminophen (APAP) overdose, a major cause of ALF because of its extensive worldwide use as an analgesic. Despite intensive efforts utilising current animal and in vitro models, the mechanisms involved in the hepatotoxicity of APAP are still not fully understood. In this expert Consensus Statement, which is endorsed by the European Drug-Induced Liver Injury Network, we aim to facilitate and outline clinically impactful discoveries by detailing the requirements for more realistic human-based systems to assess hepatotoxicity and guide future drug safety testing. We present novel insights and discuss major players in APAP pathophysiology, and describe emerging in vitro and in vivo pre-clinical models, as well as advanced imaging and in silico technologies, which may improve prediction of clinical outcomes of DILI.

Design, synthesis and anticancer/antiestrogenic activities of novel indole-benzimidazoles.

Indole-benzimidazoles have recently gained attention due to their antiproliferative and antiestrogenic effects. However, their structural similarities and molecular mechanisms shared with selective estrogen receptor modulators (SERMs) have not yet been investigated. In this study, we synthesized novel ethylsulfonyl indole-benzimidazole derivatives by substituting the first (R1) and fifth (R2) positions of benzimidazole and indole groups, respectively. Subsequently, we performed 1H NMR, 13C NMR, and Mass spectral and in silico docking analyses, and anticancer activity screening studies of these novel indole-benzimidazoles. The antiproliferative effects of indole-benzimidazoles were found to be more similar between the estrogen (E2) responsive cell lines MCF-7 and HEPG2 in comparison to the Estrogen Receptor negative (ER-) cell line MDA-MB-231. R1:p-fluorobenzyl group members were selected as lead compounds for their potent anticancer effects and moderate structural affinity to ER. Microarray expression profiling and gene enrichment analyses (GSEA) of the selected compounds (R1:p-fluorobenzyl: 48, 49, 50, 51; R1:3,4-difluorobenzyl: 53) helped determine the similarly modulated cellular signaling pathways among derivatives. Moreover, we identified known compounds that have significantly similar gene signatures to that of 51 via queries performed in LINCS database; and further transcriptomics comparisons were made using public GEO datasets (GSE35428, GSE7765, GSE62673). Our results strongly demonstrate that these novel indole-benzimidazoles can modulate ER target gene expression as well as dioxin-mediated aryl hydrocarbon receptor and amino acid deprivation-mediated integrated stress response signaling in a dose-dependent manner.

Lower connectivity of tumor coexpression networks is not specific to cancer.

Global level network analysis of molecular links is necessary for systems level view of complex diseases like cancer. Using genome-wide expression datasets, we constructed and compared gene co-expression based specific networks of pre-cancerous tumors (adenoma) and cancerous tumors (carcinoma) with paired normal networks to assess for any possible changes in network connectivity. Previously, loss of connectivity was reported as a characteristics of cancer samples. Here, we observed that pre-cancerous conditions also had significantly less connections than paired normal samples. We observed a loss of connectivity trend for colorectal adenoma, aldosterone producing adenoma and uterine leiomyoma. We also showed that the loss of connectivity trend is not specific to positive or negative correlation based networks. Differential hub genes, which were the most highly differentially less connected genes in tumor, were mostly different between different datasets. No common gene list could be defined which underlies the lower connectivity of tumor specific networks. Connectivity of colorectal cancer methylation targets was different from other genes. Extracellular space related terms were enriched in negative correlation based differential hubs and common methylation targets of colorectal carcinoma. Our results indicate a systems level change of lower connectivity as cells transform to not only cancer but also pre-cancerous conditions. This systems level behavior could not be attributed to a group of genes.

PKNOX2 expression and regulation in the bone marrow mesenchymal stem cells of Fanconi anemia patients and healthy donors.

HOX and TALE transcription factors are important regulators of development and homeostasis in determining cellular identity. Deregulation of this process may drive cancer progression. The aim of this study was to investigate the expression of these transcription factors in the bone marrow derived mesenchymal stem cells (BM-MSCs) of Fanconi anemia (FA) patients, which is a cancer-predisposing disease. Expression levels of HOX and TALE genes in BM-MSCs were obtained from FA patients and healthy donors by RT-qPCR and highly conserved expression levels were observed between patient and donor cells, except PKNOX2, which is a member of TALE class. PKNOX2 was significantly downregulated in FA cells compared to donors (P < 0.05). PKNOX2 expression levels did not change with diepoxybutane (DEB), a DNA crosslinking agent, in either donor or FA cells except one patient's with a truncation mutation of FANCA. A difference of PKNOX2 protein level was not obtained between FA patient and donor BM-MSCs by western blot analysis. When human TGF-ß1 (rTGF-ß1) recombinant protein was provided to the cultures, PKNOX2 as well as TGF-ß1 expression increased both in FA and donor BM-MSCs in a dose dependent manner. 5 ng/mL rTGF-ß stimulation had more dominant effect on the gene expression of donor BM-MSCs compared to FA cells. Decreased PKNOX2 expression in FA BM-MSCs may provide new insights into the molecular pathophysiology of the disease and TGF-ß1 levels of the microenvironment may be the cause of PKNOX2 downregulation.

Enhancer cooperativity as a novel mechanism underlying the transcriptional regulation of E-cadherin during mesenchymal to epithelial transition.

Epithelial-mesenchymal transition (EMT) and mesenchymal-epithelial transition (MET) highlight crucial steps during embryogenesis and tumorigenesis. Induction of dramatic changes in gene expression and cell features is reflected by modulation of Cdh1 (E-cadherin) expression. We show that Cdh1 activity during MET is governed by two enhancers at +7.8 kb and at +11.5 kb within intron 2 that are activated by binding of Grhl3 and Hnf4α, respectively. Recruitment of Grhl3 and Hnf4α to the enhancers is crucial for activating Cdh1 and accomplishing MET in non-tumorigenic mouse mammary gland cells (NMuMG). Moreover, the two enhancers cooperate via Grhl3 and Hnf4α binding, induction of DNA-looping and clustering at the promoter to orchestrate E-cadherin re-expression. Our results provide novel insights into the cellular mechanisms whereby cells respond to MET signals and re-establish an epithelial phenotype by enhancer cooperativity. A general importance of our findings including MET-mediated colonization of metastasizing tumor cells is suggested.

Characterization of a novel zebrafish (Danio rerio) gene, wdr81, associated with cerebellar ataxia, mental retardation and dysequilibrium syndrome (CAMRQ).

BACKGROUND: WDR81 (WD repeat-containing protein 81) is associated with cerebellar ataxia, mental retardation and disequilibrium syndrome (CAMRQ2, [MIM 610185]). Human and mouse studies suggest that it might be a gene of importance during neurodevelopment. This study aimed at fully characterizing the structure of the wdr81 transcript, detecting the possible transcript variants and revealing its expression profile in zebrafish, a powerful model organism for studying development and disease. RESULTS: As expected in human and mouse orthologous proteins, zebrafish wdr81 is predicted to possess a BEACH (Beige and Chediak-Higashi) domain, a major facilitator superfamily domain and WD40-repeats, which indicates a conserved function in these species. We observed that zebrafish wdr81 encodes one open reading frame while the transcript has one 5' untranslated region (UTR) and the prediction of the 3' UTR was mainly confirmed along with a detected insertion site in the embryo and adult brain. This insertion site was also found in testis, heart, liver, eye, tail and muscle, however, there was no amplicon in kidney, intestine and gills, which might be the result of possible alternative polyadenylation processes among tissues. The 5 and 18 hpf were critical timepoints of development regarding wdr81 expression. Furthermore, the signal of the RNA probe was stronger in the eye and brain at 18 and 48 hpf, then decreased at 72 hpf. Finally, expression of wdr81 was detected in the adult brain and eye tissues, including but not restricted to photoreceptors of the retina, presumptive Purkinje cells and some neurogenic brains regions. CONCLUSIONS: Taken together these data emphasize the importance of this gene during neurodevelopment and a possible role for neuronal proliferation. Our data provide a basis for further studies to fully understand the function of wdr81.

Red emitting, cucurbituril-capped, pH-responsive conjugated oligomer-based nanoparticles for drug delivery and cellular imaging.

Here we report the synthesis of nanoparticles based on a conjugated oligomer which is synthesized through Heck-coupling of divinylfluorene and dibromobenzothiodiazole monomers. These water dispersible nanoparticles emit in the region of red tailing to the near-infrared region of the spectrum with high fluorescent quantum yield and brightness. The nanoparticles were found to be stable in water for a prolonged time without forming any aggregates and could carry camptothecin, an anticancer drug with high loading efficiency. MTT cell viability studies performed with breast cancer cell lines showed that half-maximal inhibitory concentration (IC50) values of nanoparticles for MCF7 and MDA-MB-231 were 44.7 µM and 24.8 µM, respectively. In order to further decrease the cytotoxicity and increase the stability of nanoparticles, amine groups were disguised by capping with cucurbit[7]uril (CB7). Drug release studies showed that drugs were released at low pH (at 5.0) faster than physiological pH (7.4) confirming the pH-responsive nature of the nanoparticles. On the other hand, CB7-capped drug-loaded nanoparticles regulated the release rate by providing slower release at pH 7.4 than the nanoparticles in the absence of CB7s. IC50 values for camptothecin in the presence of nanoparticles with or without CB7 were significantly reduced in MCF7 and MDA-MB-231 cells.

Synthesis and Structure of Novel Phenothiazine Derivatives, and Compound Prioritization via In Silico Target Search and Screening for Cytotoxic and Cholinesterase Modulatory Activities in Liver Cancer Cells and In Vivo in Zebrafish.

Phenothiazines (PTZ) are antipsychotics known to modulate a variety of neurotransmitter activities that include dopaminergic and cholinergic signaling and have been identified as potential anticancer agents in vitro. However, it is important to also test whether a highly cytotoxic, repurposed, or novel PTZ has low toxicity and neuromodulatory activity in vivo using vertebrate model organisms, such as zebrafish. In this study, we synthesized novel phenothiazines and screened them in vitro in liver cancer and in vivo in zebrafish embryos/larvae. The syntheses of several intermediate PTZ 10-yl acyl chlorides were followed by elemental analysis and determination of 1H NMR and 13C NMR mass (ESI+) spectra of a large number of novel PTZ 10-carboxamides. Cytotoxicities of 28 PTZ derivatives (1-28) screened against Hep3B and SkHep1 liver cancer cell lines revealed five intermediate and five novel leads along with trifluoperazine (TFP), prochlorperazine (PCP), and perphenazine, which are relatively more cytotoxic than the basic PTZ core. Overall, the derivatives were more cytotoxic to Hep3B than SkHep1 cells. Moreover, in silico target screening identified cholinesterases as some of the commonest targets of the screened phenothiazines. Interestingly, molecular docking studies with acetylcholinesterase (AChE) and butyrylcholinesterase proteins showed that the most cytotoxic compounds 1, 3, PCP, and TFP behaved similar to Huprin W in their amino acid interactions with the AChE protein. The highly cytotoxic intermediate PTZ derivative 1 exhibited a relatively lower toxicity profile than those of 2 and 3 during the zebrafish development. It also modulated in vivo the cholinesterase activity in a dose-dependent manner while significantly increasing the total cholinesterase activity and/or ACHE mRNA levels, independent of the liver cancer cell type. Our screen also identified novel phenothiazines, i.e., 8 and 10, with significant cytotoxic and cholinesterase modulatory effects in liver cancer cells; yet both compounds had low levels of toxicity in zebrafish. Moreover, they modulated the cholinesterase activity or expression of ACHE in a cancer cell line-specific manner, and compound 10 significantly inhibited the cholinesterase activity in zebrafish. Accordingly, using a successful combination of in silico, in vitro, and in vivo approaches, we identified several lead anticancer and cholinesterase modulatory PTZ derivatives for future research.

Cancer-testis gene expression is associated with the methylenetetrahydrofolate reductase 677 C>T polymorphism in non-small cell lung carcinoma.

BACKGROUND: Tumor-specific, coordinate expression of cancer-testis (CT) genes, mapping to the X chromosome, is observed in more than 60% of non-small cell lung cancer (NSCLC) patients. Although CT gene expression has been unequivocally related to DNA demethylation of promoter regions, the underlying mechanism leading to loss of promoter methylation remains elusive. Polymorphisms of enzymes within the 1-carbon pathway have been shown to affect S-adenosyl methionine (SAM) production, which is the sole methyl donor in the cell. Allelic variants of several enzymes within this pathway have been associated with altered SAM levels either directly, or indirectly as reflected by altered levels of SAH and Homocysteine levels, and altered levels of DNA methylation. We, therefore, asked whether the five most commonly occurring polymorphisms in four of the enzymes in the 1-carbon pathway associated with CT gene expression status in patients with NSCLC. METHODS: Fifty patients among a cohort of 763 with NSCLC were selected based on CT gene expression status and typed for five polymorphisms in four genes known to affect SAM generation by allele specific q-PCR and RFLP. RESULTS: We identified a significant association between CT gene expression and the MTHFR 677 CC genotype, as well as the C allele of the SNP, in this cohort of patients. Multivariate analysis revealed that the genotype and allele strongly associate with CT gene expression, independent of potential confounders. CONCLUSIONS: Although CT gene expression is associated with DNA demethylation, in NSCLC, our data suggests this is unlikely to be the result of decreased MTHFR function.

mESAdb: microRNA expression and sequence analysis database.

microRNA expression and sequence analysis database (http://konulab.fen.bilkent.edu.tr/mirna/) (mESAdb) is a regularly updated database for the multivariate analysis of sequences and expression of microRNAs from multiple taxa. mESAdb is modular and has a user interface implemented in PHP and JavaScript and coupled with statistical analysis and visualization packages written for the R language. The database primarily comprises mature microRNA sequences and their target data, along with selected human, mouse and zebrafish expression data sets. mESAdb analysis modules allow (i) mining of microRNA expression data sets for subsets of microRNAs selected manually or by motif; (ii) pair-wise multivariate analysis of expression data sets within and between taxa; and (iii) association of microRNA subsets with annotation databases, HUGE Navigator, KEGG and GO. The use of existing and customized R packages facilitates future addition of data sets and analysis tools. Furthermore, the ability to upload and analyze user-specified data sets makes mESAdb an interactive and expandable analysis tool for microRNA sequence and expression data.

Global miRNA expression of bone marrow mesenchymal stem/stromal cells derived from Fanconi anemia patients.

Fanconi anemia (FA) is a rare genetic disorder characterized by genomic instability, developmental defects, and bone marrow (BM) failure. Hematopoietic stem cells (HSCs) in BM interact with the mesenchymal stem/stromal cells (MSCs); and this partly sustains the tissue homeostasis. MicroRNAs (miRNAs) can play a critical role during these interactions possibly via paracrine mechanisms. This is the first study addressing the miRNA profile of FA BM-MSCs obtained before and after BM transplantation (preBMT and postBMT, respectively). Non-coding RNA expression profiling and quality control analyses were performed in Donors (n = 13), FA preBMT (n = 11), and FA postBMT (n = 6) BM-MSCs using GeneChip miRNA 2.0 Array. Six Donor-FA preBMT pairs were used to identify a differentially expressed miRNA expression signature containing 50 miRNAs, which exhibited a strong correlation with the signature obtained from unpaired samples. Five miRNAs (hsa-miR-146a-5p, hsa-miR-148b-3p, hsa-miR-187-3p, hsa-miR-196b-5p, and hsa-miR-25-3p) significantly downregulated in both the paired and unpaired analyses were used to generate the BM-MSCs' miRNA-BM mononuclear mRNA networks upon integration of a public dataset (GSE16334; studying Donor versus FA samples). Functionally enriched KEGG pathways included cellular senescence, miRNAs, and pathways in cancer. Here, we showed that hsa-miR-146a-5p and hsa-miR-874-3p were rescued upon BMT (n = 3 triplets). The decrease in miR-146a-5p was also validated using RT-qPCR and emerged as a strong candidate as a modulator of BM mRNAs in FA patients.

Doxorubicin induces prolonged DNA damage signal in cells overexpressing DEK isoform-2.

DEK has a short isoform (DEK isoform-2; DEK2) that lacks amino acid residues between 49-82. The full-length DEK (DEK isoform-1; DEK1) is ubiquitously expressed and plays a role in different cellular processes but whether DEK2 is involved in these processes remains elusive. We stably overexpressed DEK2 in human bone marrow stromal cell line HS-27A, in which endogenous DEKs were intact or suppressed via short hairpin RNA (sh-RNA). We have found that contrary to ectopic DEK1, DEK2 locates in the nucleus and nucleolus, causes persistent γH2AX signal upon doxorubicin treatment, and couldn't functionally compensate for the loss of DEK1. In addition, DEK2 overexpressing cells were more sensitive to doxorubicin than DEK1-cells. Expressions of DEK1 and DEK2 in cell lines and primary tumors exhibit tissue specificity. DEK1 is upregulated in cancers of the colon, liver, and lung compared to normal tissues while both DEK1 and DEK2 are downregulated in subsets of kidney, prostate, and thyroid carcinomas. Interestingly, only DEK2 was downregulated in a subset of breast tumors suggesting that DEK2 can be modulated differently than DEK1 in specific cancers. In summary, our findings show distinct expression patterns and subcellular location and suggest non-overlapping functions between the two DEK isoforms.

Dose- and time-dependent expression patterns of zebrafish orthologs of selected E2F target genes in response to serum starvation/replenishment.

Targets of E2F transcription factors effectively regulate the cell cycle from worms to humans. Furthermore, the dysregulation of E2F transcription modules plays a highly conserved role in cancers of human and zebrafish. Studying E2F target expression under a given cellular state, such as quiescence, might lead to a better understanding of the conserved patterns of expression in different taxa. In the present study, we used literature searches and phylogeny to identify several targets of E2F transcription factors that are known to be serum-responsive; namely, PCNA, MYBL2, MCM7, TYMS, and CTGF. The transcriptional serum response of zebrafish orthologs of these genes were quantified under different doses (i.e., 0, 0.1, 1, 3, and 10% FBS) and time points (i.e., 6, 24 and 48 hours, h) using quantitative RT-PCR (qRT-PCR) in the zebrafish fibroblast cells (ZF4). Our results indicated that mRNA expression of zebrafish pcna, mybl2, mcm7 and tyms drastically decreased while that of ctgf increased with decreasing serum levels as observed in mammals. These genes responded to serum starvation at 24 and 48 h and to the mitogenic stimuli as early as 6 h except for ctgf whose expression was significantly altered at 24 h. The zebrafish Mcm7 protein levels also were modulated by serum starvation/replenishment. The present study provides a foundation for the comparative analysis of quantitative expression patterns for genes involved in regulation of cell cycle using a zebrafish serum response model.

Zebrafish Xenotransplantation Models for Studying Gene Function and Drug Treatment in Hepatocellular Carcinoma.

INTRODUCTION: Zebrafish is a promising model organism for human disease including hepatocellular cancer (HCC). Recently, zebrafish has emerged also as a host for xenograft studies of liver cancer cell lines and patient derived tumors of HCC. Zebrafish embryos enable drug screening and gene function studies of xenografted cells via ease of microinjection and visualization of tumor growth and metastasis. OBJECTIVES: In this review, we aimed to overview zebrafish HCC and liver cancer xenotransplantation studies focusing on 'gene functional analysis' and 'drug/chemical screening'. METHODS: Herein, a comprehensive literature search was performed for liver and HCC xenografts in zebrafish on PubMed using different key words and filters for molecular modifications or drug exposure. RESULTS: Our literature search revealed around 250 studies which were filtered and summarized in a table (Table 1) revealing comprehensive collection of experimental and technical details on microinjection, injected cell lines, molecular modifications of injected cells, types and doses of drug treatments as well as biological assessments. CONCLUSION: This review provides a platform for HCC and liver xenografts and highlights studies performed to understand gene functionality and drug efficacy in vivo in zebrafish.

Discovery of Cancer-Specific and Independent Prognostic Gene Subsets of the Slit-Robo Family Using TCGA-PANCAN Datasets.

The Slit-Robo family of axon guidance molecules works in concert, playing important roles in organ development and cancer. Expressions of individual Slit-Robo genes have been used in calculating univariable hazard ratios (HRuni) for predicting cancer prognosis in the literature. However, Slit-Robo members do not act independently; hence, hazard ratios from multivariable Cox regression (HRmulti) on the whole gene set can further lead to identification of cancer-specific, novel, and independent prognostic gene pairs or modules. Herein, we obtained mRNA expressions of the Slit-Robo family consisting of four Robos (ROBO1/2/3/4) and three Slits (SLIT1/2/3), along with four types of survival outcome across cancers found in the Cancer Genome Atlas (TCGA). We used cluster heat maps to visualize closely associated pairs/modules of prognostic genes across 33 different cancers. We found a smaller number of significant genes in HRmulti than in HRuni, suggesting that the former analysis was less redundant. High ROBO4 expression emerged as relatively protective within the family, in both types of HR analyses. Multivariable Cox regression, on the other hand, revealed significantly more HR signatures containing Slit-Robo pairs acting in opposing directions than those containing Slit-Slit or Robo-Robo pairs for disease-specific survival. Furthermore, we discovered, through the online app SmulTCan's lasso regression, Slit-Robo gene subsets that significantly differentiated between high- versus low-risk prognosis patient groups, particularly for renal cancers and low-grade glioma. The statistical pipeline reported herein can help test independent and significant pairs/modules within a codependent gene family for cancer prognostication, and thus should also prove useful in personalized/precision medicine research.