Activated oncogenic pathways and therapeutic targets in extranodal nasal-type NK/T cell lymphoma revealed by gene expression profiling

We performed comprehensive genome-wide gene expression profiling (GEP) of extranodal nasal-type natural killer/T-cell lymphoma (NKTL) using formalin-fixed, paraffin-embedded tissue (n = 9) and NK cell lines (n = 5) in comparison with normal NK cells, with the objective of understanding the oncogenic pathways involved in the pathogenesis of NKTL and to identify potential therapeutic targets. Pathway and network analysis of genes differentially expressed between NKTL and normal NK cells revealed significant enrichment for cell cycle-related genes and pathways, such as PLK1, CDK1, and Aurora-A. Furthermore, our results demonstrated a pro-proliferative and anti-apoptotic phenotype in NKTL characterized by activation of Myc and nuclear factor kappa B (NF-kappa B), and deregulation of p53. In corroboration with GEP findings, a significant percentage of NKTLs (n = 33) overexpressed c-Myc (45.4%), p53 (87.9%), and NF-kappa B p50 (67.7%) on immunohistochemistry using a tissue microarray containing 33 NKTL samples. Notably, overexpression of survivin was observed in 97% of cases. Based on our findings, we propose a model of NKTL pathogenesis where deregulation of p53 together with activation of Myc and NF-kappa B, possibly driven by EBV LMP-1, results in the cumulative up-regulation of survivin. Down-regulation of survivin with Terameprocol (EM-1421, a survivin inhibitor) results in reduced cell viability and increased apoptosis in tumour cells, suggesting that targeting survivin may be a potential novel therapeutic strategy in NKTL. Copyright (C) 2011 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Molecular similarity between myelodysplastic form of chronic myelomonocytic leukemia and refractory anemia with ring sideroblasts

Chronic myelomonocytic leukemia is similar to but a separate entity from both myeloproliferative neoplasms and myelodysplastic syndromes, and shows either myeloproliferative or myelodysplastic features. We ask whether this distinction may have a molecular basis. We established the gene expression profiles of 39 samples of chronic myelomonocytic leukemia (including 12 CD34-positive) and 32 CD34-positive samples of myelodysplastic syndromes by using Affymetrix microarrays, and studied the status of 18 genes by Sanger sequencing and array-comparative genomic hybridization in 53 samples. Analysis of 12 mRNAS from chronic myelomonocytic leukemia established a gene expression signature of 122 probe sets differentially expressed between proliferative and dysplastic cases of chronic myelomonocytic leukemia. As compared to proliferative cases, dysplastic cases over-expressed genes involved in red blood cell biology. When applied to 32 myelodysplastic syndromes, this gene expression signature was able to discriminate refractory anemias with ring sideroblasts from refractory anemias with excess of blasts. By comparing mRNAS from these two forms of myelodysplastic syndromes we derived a second gene expression signature. This signature separated the myelodysplastic and myeloproliferative forms of chronic myelomonocytic leukemias. These results were validated using two independent gene expression data sets. We found that myelodysplastic chronic myelomonocytic leukemias are characterized by mutations in transcription/epigenetic regulators (ASXL1, RUNX1, TET2) and splicing genes (SRSF2) and the absence of mutations in signaling genes. Myelodysplastic chronic myelomonocytic leukemias and refractory anemias with ring sideroblasts share a common expression program suggesting they are part of a continuum, which is not totally explained by their similar but not, however, identical mutation spectrum. © 2013 Ferrata Storti Foundation.

ProP Is Required for the Survival of Desiccated Salmonella enterica Serovar Typimurium Cells on a Stainless Steel Surface

Consumers trust commercial food production to be safe, and it is important to strive to improve food safety at every level. Several outbreaks of food-borne disease have been caused by Salmonella strains associated with dried food. Currently we do not know the mechanisms used by Salmonella enterica serovar Typhimurium to survive in desiccated environments. The aim of this study was to discover the responses of S. Typhimurium ST4/74 at the transcriptional level to desiccation on a stainless steel surface and to subsequent rehydration. Bacterial cells were dried onto the same steel surfaces used during the production of dry foods, and RNA was recovered for transcriptomic analysis. Subsequently, dried cells were rehydrated and were again used for transcriptomic analysis. A total of 266 genes were differentially expressed under desiccation stress compared with a static broth culture. The osmoprotectant transporters proP, proU, and osmU (STM1491 to STM1494) were highly upregulated by drying. Deletion of any one of these transport systems resulted in a reduction in the long-term viability of S. Typhimurium on a stainless steel food contact surface. The proP gene was critical for survival; proP deletion mutants could not survive desiccation for long periods and were undetectable after 4 weeks. Following rehydration, 138 genes were differentially expressed, with upregulation observed for genes such as proP, proU, and the phosphate transport genes (pstACS). In time, this knowledge should prove valuable for understanding the underlying mechanisms involved in pathogen survival and should lead to improved methods for control to ensure the safety of intermediate-and low-moisture foods. © 2013, American Society for Microbiology.

Connectivity Mapping for Candidate Therapeutics Identification Using Next Generation Sequencing RNA-Seq Data

The advent of next generation sequencing technologies (NGS) has expanded the area of genomic research, offering high coverage and increased sensitivity over older microarray platforms. Although the current cost of next generation sequencing is still exceeding that of microarray approaches, the rapid advances in NGS will likely make it the platform of choice for future research in differential gene expression. Connectivity mapping is a procedure for examining the connections among diseases, genes and drugs by differential gene expression initially based on microarray technology, with which a large collection of compound-induced reference gene expression profiles have been accumulated. In this work, we aim to test the feasibility of incorporating NGS RNA-Seq data into the current connectivity mapping framework by utilizing the microarray based reference profiles and the construction of a differentially expressed gene signature from a NGS dataset. This would allow for the establishment of connections between the NGS gene signature and those microarray reference profiles, alleviating the associated incurring cost of re-creating drug profiles with NGS technology. We examined the connectivity mapping approach on a publicly available NGS dataset with androgen stimulation of LNCaP cells in order to extract candidate compounds that could inhibit the proliferative phenotype of LNCaP cells and to elucidate their potential in a laboratory setting. In addition, we also analyzed an independent microarray dataset of similar experimental settings. We found a high level of concordance between the top compounds identified using the gene signatures from the two datasets. The nicotine derivative cotinine was returned as the top candidate among the overlapping compounds with potential to suppress this proliferative phenotype. Subsequent lab experiments validated this connectivity mapping hit, showing that cotinine inhibits cell proliferation in an androgen dependent manner. Thus the results in this study suggest a promising prospect of integrating NGS data with connectivity mapping. © 2013 McArt et al.

Trait directed de novo population transcriptome dissects genetic regulation of a balanced polymorphism in phosphorus nutrition/arsenate tolerance in a wild grass Holcus lanatus

The aim of this study was to characterize the transcriptome of a balanced polymorphism, under the regulation of a single gene, for phosphate fertilizer responsiveness/arsenate toler- ance in wild grass Holcus lanatus genotypes screened from the same habitat.

De novo transcriptome sequencing, RNAseq (RNA sequencing) and single nucleotide poly- morphism (SNP) calling were conducted on RNA extracted from H.lanatus. Roche 454 sequencing data were assembled into c. 22 000 isotigs, and paired-end Illumina reads for phosphorus-starved (P) and phosphorus-treated (P+) genovars of tolerant (T) and nontoler- ant (N) phenotypes were mapped to this reference transcriptome.

Heatmaps of the gene expression data showed strong clustering of each P+/P treated genovar, as well as clustering by N/T phenotype. Statistical analysis identified 87 isotigs to be significantly differentially expressed between N and T phenotypes and 258 between P+ and P treated plants. SNPs and transcript expression that systematically differed between N and T phenotypes had regulatory function, namely proteases, kinases and ribonuclear RNA- binding protein and transposable elements.

A single gene for arsenate tolerance led to distinct phenotype transcriptomes and SNP pro- files, with large differences in upstream post-translational and post-transcriptional regulatory genes rather than in genes directly involved in P nutrition transport and metabolism per se.

Genetic and genomic analyses of musculoskeletal differences between BEH and BEL strains

Berlin high (BEH) and Berlin low (BEL) strains selected for divergent growth differ 3-fold in body weight. We aimed at examining muscle mass, which is a major contributor to body weight, by exploring anatomical characteristics of the soleus muscle, its fiber numbers and their cross sectional area (CSA), by analysing transcriptome of the gastrocnemius and by initiating quantitative trait locus (QTL) mapping. BEH muscles were 4-to-8 times larger compared to BEL strain. In sub-strain BEH+/+, mutant myostatin was replaced with a wild type allele, however, BEH+/+muscles still were 2-to-4 times larger compared to the BEL strain. BEH soleus contained 2-times more (P<0.0001) and 2-times larger in CSA (P<0.0001) fibers compared to BEL strain. In addition, soleus femoral attachment anomaly (SFAA) was observed in all BEL mice. One significant (chromosome 1) and four suggestive (chromosomes 3, 4, 6 and 9) muscle weight QTLs were mapped in 21-day old F2 intercross (n=296) between BEH and BEL strains. The frequency of SFAA incidence in the F2 and in the backcross to BEL strain (BCL) suggested the presence of more than one causative gene. Two suggestive SFAA QTLs were mapped in BCL, however, their peak markers were not associated with the phenotype in F2. RNA-Seq analysis revealed 2,148 differentially expressed (P<0.1) genes and 45,673 SNPs and >2,000 indels between BEH+/+ and BEL males. In conclusion, contrasting muscle traits, genomic and gene expression differences between BEH and BEL strains provide a promising model for the search of genes involved in muscle growth and musculoskeletal morphogenesis.

Re-analysis of microarray data reveals insights into altered transcriptional activity of TH17 and TReg signalling in psoriasis

Identifying differential expression of genes in psoriatic and healthy skin by microarray data analysis is a key approach to understand the pathogenesis of psoriasis. Analysis of more than one dataset to identify genes commonly upregulated reduces the likelihood of false positives and narrows down the possible signature genes. Genes controlling the critical balance between T helper 17 and regulatory T cells are of special interest in psoriasis. Our objectives were to identify genes that are consistently upregulated in lesional skin from three published microarray datasets. We carried out a reanalysis of gene expression data extracted from three experiments on samples from psoriatic and nonlesional skin using the same stringency threshold and software and further compared the expression levels of 92 genes related to the T helper 17 and regulatory T cell signaling pathways. We found 73 probe sets representing 57 genes commonly upregulated in lesional skin from all datasets. These included 26 probe sets representing 20 genes that have no previous link to the etiopathogenesis of psoriasis. These genes may represent novel therapeutic targets and surely need more rigorous experimental testing to be validated. Our analysis also identified 12 of 92 genes known to be related to the T helper 17 and regulatory T cell signaling pathways, and these were found to be differentially expressed in the lesional skin samples.

Gene Sets Net Correlations Analysis (GSNCA): a multivariate differential coexpression test for gene sets

Motivation: To date, Gene Set Analysis (GSA) approaches primarily focus on identifying differentially expressed gene sets (pathways). Methods for identifying differentially coexpressed pathways also exist but are mostly based on aggregated pairwise correlations, or other pairwise measures of coexpression. Instead, we propose Gene Sets Net Correlations Analysis (GSNCA), a multivariate differential coexpression test that accounts for the complete correlation structure between genes.

Results: In GSNCA, weight factors are assigned to genes in proportion to the genes' cross-correlations (intergene correlations). The problem of finding the weight vectors is formulated as an eigenvector problem with a unique solution. GSNCA tests the null hypothesis that for a gene set there is no difference in the weight vectors of the genes between two conditions. In simulation studies and the analyses of experimental data, we demonstrate that GSNCA, indeed, captures changes in the structure of genes' cross-correlations rather than differences in the averaged pairwise correlations. Thus, GSNCA infers differences in coexpression networks, however, bypassing method-dependent steps of network inference. As an additional result from GSNCA, we define hub genes as genes with the largest weights and show that these genes correspond frequently to major and specific pathway regulators, as well as to genes that are most affected by the biological difference between two conditions. In summary, GSNCA is a new approach for the analysis of differentially coexpressed pathways that also evaluates the importance of the genes in the pathways, thus providing unique information that may result in the generation of novel biological hypotheses.

The prognostic significance of PD-L1 in bladder cancer

Immunotherapy is a promising strategy for the treatment of various types of cancer. An antibody that targets programmed death ligand-1 (PD-L1) pathway has been shown to be active towards various types of cancer, including melanoma and lung cancer. MPDL3280A, an anti‑PD-L1 antibody, has shown clear clinical activity in PD-L1-overexpressing bladder cancer with an objective response rate of 40-50%, resulting in a breakthrough therapy designation granted by FDA. These events pronounce the importance of targeting the PD-L1 pathway in the treatment of bladder cancer. In the present study, we investigated the prognostic significance of the expression of three genes in the PD-L1 pathway, including PD-L1, B7.1 and PD-1, in three independent bladder cancer datasets in the Gene Expression Omnibus database. PD-L1, B7.1 and PD-1 were significantly associated with clinicopathological parameters indicative of a more aggressive phenotype of bladder cancer, such as a more advanced stage and a higher tumor grade. In addition, a high level expression of PD-L1 was associated with reduced patient survival. Of note, the combination of PD-L1 and B7.1 expression, but not other combinations of the three genes, were also able to predict patient survival. Our findings support the development of anti-PD-L1, which blocks PD-L1-PD-1 and B7.1-PD-L1 interactions, in treatment of bladder cancer. The observations were consistent in the three independent bladder cancer datasets consisting of a total of 695 human bladder specimens. The datasets were then assessed and it was found that the expression levels of the chemokine CC-motif ligand (CCL), CCL3, CCL8 and CCL18, were correlated with the PD-L1 expression level, while ADAMTS13 was differentially expressed in patients with a different survival status (alive or deceased). Additional investigations are required to elucidate the role of these genes in the PD-L1-mediated immune system suppression and bladder cancer progression. In conclusion, findings of this study suggested that PD-L1 is an important prognostic marker and a therapeutic target for bladder cancer.

Pro-neural transcription factors as cancer markers

BACKGROUND: The aberrant transcription in cancer of genes normally associated with embryonic tissue differentiation at various organ sites may be a hallmark of tumour progression. For example, neuroendocrine differentiation is found more commonly in cancers destined to progress, including prostate and lung. We sought to identify proteins which are involved in neuroendocrine differentiation and differentially expressed in aggressive/metastatic tumours.

RESULTS: Expression arrays were used to identify up-regulated transcripts in a neuroendocrine (NE) transgenic mouse model of prostate cancer. Amongst these were several genes normally expressed in neural tissues, including the pro-neural transcription factors Ascl1 and Hes6. Using quantitative RT-PCR and immuno-histochemistry we showed that these same genes were highly expressed in castrate resistant, metastatic LNCaP cell-lines. Finally we performed a meta-analysis on expression array datasets from human clinical material. The expression of these pro-neural transcripts effectively segregates metastatic from localised prostate cancer and benign tissue as well as sub-clustering a variety of other human cancers.

CONCLUSION: By focussing on transcription factors known to drive normal tissue development and comparing expression signatures for normal and malignant mouse tissues we have identified two transcription factors, Ascl1 and Hes6, which appear effective markers for an aggressive phenotype in all prostate models and tissues examined. We suggest that the aberrant initiation of differentiation programs may confer a selective advantage on cells in all contexts and this approach to identify biomarkers therefore has the potential to uncover proteins equally applicable to pre-clinical and clinical cancer biology.

Genome Expression Profiling-Based Identification and Administration Efficacy of Host-Directed Antimicrobial Drugs against Respiratory Infection by Nontypeable Haemophilus influenzae

Therapies that are safe, effective, and not vulnerable to developing resistance are highly desirable to counteract bacterial infections. Host-directed therapeutics is an antimicrobial approach alternative to conventional antibiotics based on perturbing host pathways subverted by pathogens during their life cycle by using host-directed drugs. In this study, we identified and evaluated the efficacy of a panel of host-directed drugs against respiratory infection by nontypeable Haemophilus influenzae (NTHi). NTHi is an opportunistic pathogen that is an important cause of exacerbation of chronic obstructive pulmonary disease (COPD). We screened for host genes differentially expressed upon infection by the clinical isolate NTHi375 by analyzing cell whole-genome expression profiling and identified a repertoire of host target candidates that were pharmacologically modulated. Based on the proposed relationship between NTHi intracellular location and persistence, we hypothesized that drugs perturbing host pathways used by NTHi to enter epithelial cells could have antimicrobial potential against NTHi infection. Interfering drugs were tested for their effects on bacterial and cellular viability, on NTHi-epithelial cell interplay, and on mouse pulmonary infection. Glucocorticoids and statins lacked in vitro and/or in vivo efficacy. Conversely, the sirtuin-1 activator resveratrol showed a bactericidal effect against NTHi, and the PDE4 inhibitor rolipram showed therapeutic efficacy by lowering NTHi375 counts intracellularly and in the lungs of infected mice. PDE4 inhibition is currently prescribed in COPD, and resveratrol is an attractive geroprotector for COPD treatment. Together, these results expand our knowledge of NTHi-triggered host subversion and frame the antimicrobial potential of rolipram and resveratrol against NTHi respiratory infection.

Natural Killer Cell-like signature observed in locally advanced rectal cancer after neoadjuvant chemoradiotherapy in a tumour regression grade dependent manner

Background: Around 10-15% of patients with locally advanced rectal cancer (LARC) undergo a pathologically complete response (TRG4) to neoadjuvant chemoradiotherapy; the rest of patients exhibit a spectrum of tumour regression (TRG1-3). Understanding therapy-related genomic alterations may help us to identify underlying biology or novel targets associated with response that could increase the efficacy of therapy in patients that do not benefit from the current standard of care.
Methods: 48 FFPE rectal cancer biopsies and matched resections were analysed using the WG-DASL HumanHT-12_v4 Beadchip array on the illumina iScan. Bioinformatic analysis was conducted in Partek genomics suite and R studio. Limma and glmnet packages were used to identify genes differentially expressed between tumour regression grades. Validation of microarray results will be carried out using IHC, RNAscope and RT-PCR.
Results: Immune response genes were observed from supervised analysis of the biopsies which may have predictive value. Differential gene expression from the resections as well as pre and post therapy analysis revealed induction of genes in a tumour regression dependent manner. Pathway mapping and Gene Ontology analysis of these genes suggested antigen processing and natural killer mediated cytotoxicity respectively. The natural killer-like gene signature was switched off in non-responders and on in the responders. IHC has confirmed the presence of Natural killer cells through CD56+ staining.
Conclusion: Identification of NK cell genes and CD56+ cells in patients responding to neoadjuvant chemoradiotherapy warrants further investigation into their association with tumour regression grade in LARC. NK cells are known to lyse malignant cells and determining whether their presence is a cause or consequence of response is crucial. Interrogation of the cytokines upregulated in our NK-like signature will help guide future in vitro models.

An integrated and comparative approach to genetic selection of the aquaculture species Dicentrarchus labrax

The European sea bass, Dicentrarchus labrax, is one of the most important marine species cultivated in Southern Europe and has not benefited from selective breeding. One of the major goals in the sea bass (D. labrax) aquaculture industry is to understand and control the complexity of growth associated traits. The aim of the methodology developed for the studies reported in the thesis was not only to establish genetic and genomic resources for sea bass, but to also develop a conceptual strategy to efficiently create knowledge in a research environment that can easily be transferred to the aquaculture industry. The strategy involved; i) establishing an annotated sea bass transcriptome and then using it to, ii) identify new genetic markers for target QTL regions so that, iii) new QTL analysis could be performed and marker based resolution of the DNA regions of interest increased, and then iv) to merge the linkage map and the physical map in order to map the QTL confidence intervals to the sea bass genome and identify genes underlying the targeted traits. Finally to test if genes in the QTL regions that are candidates for divergent growth phenotypes have modified patterns of transcription that reflects the modified whole organism physiology SuperSAGE-SOLiD4 gene expression was used with sea bass with high growth heterogeneity. The SuperSAGE contributed to significantly increase the transcriptome information for sea bass muscle, brain and liver and also led to the identification of putative candidate genes lying in the genomic region of growth related QTL. Lastly all differentially expressed transcripts in brain, liver and muscle of the European sea bass with divergent specific growth rates were mapped to gene pathways and networks and the regulatory pathways most affected identified and established the tissue specific changes underlying the divergent SGR. Owing to the importance of European sea bass to Mediterranean aquaculture and the developed genomics resources from the present thesis and from other studies it should be possible to implement genetic selection programs using marker assisted selection.

Molecular structure and functional analysis of runx3 in zebrafish

Tese de doutoramento, Ciências Biomédicas, Universidade do Algarve, Departamento de Ciências Biomédicas e Medicina, 2014

Identification of novel molecular determinants of tissue mineralization in fish

The identification of genes involved in signaling and regulatory pathways, and matrix formation is paramount to the better understanding of the complex mechanisms of bone formation and mineralization, and critical to the successful development of therapies for human skeletal disorders. To achieve this objective, in vitro cell systems derived from skeletal tissues and able to mineralize their extracellular matrix have been used to identify genes differentially expressed during mineralization and possibly new markers of bone and cartilage homeostasis. Using cell systems of fish origin and techniques such as suppression subtractive hybridization and microarray hybridization, three genes never associated with mechanisms of calcification were identified: the calcium binding protein S100-like, the short-chain dehydrogenase/reductase sdr-like and the betaine homocysteine S-methyltransferase bhmt3. Analysis of the spatial-temporal expression of these 3 genes by qPCR and in situ hybridization revealed: (1) the up-regulation of sdr-like transcript during in vitro mineralization of gilthead seabream cell lines and its specificity for calcified tissues and differentiating osteoblasts; (2) the up-regulation of S100-like and the down-regulation of bhmt3 during in vitro mineralization and the central role of both genes in cartilaginous tissues undergoing endo/perichondral mineralization in juvenile fish. While expression of S100-like and bhmt3 was restricted to calcified tissues, sdr-like transcript was also detected in soft tissues, in particular in tissues of the gastrointestinal tract. Functional analysis of gene promoters revealed the transcriptional regulation of the 3 genes by known regulators of osteoblast and chondrocyte differentiation/mineralization: RUNX2 and RAR (sdr-like), ETS1 (s100-like; bhmt3), SP1 and MEF2c (bhmt3). The evolutionary relationship of the different orthologs and paralogs identified within the scope of this work was also inferred from taxonomic and phylogenetic analyses and revealed novel protein subfamilies (S100-like and Sdr-like) and the explosive diversity of Bhmt family in particular fish groups (Neoteleostei). Altogether our results contribute with new data on SDR, S100 and BHMT proteins, evidencing for the first time the role for these three proteins in mechanisms of mineralization in fish and emphasized their potential as markers of mineralizing cartilage and bone in developing fish.