Regulation and consequences of differential gene expression in diabetic kidney disease

DIN (diabetic nephropathy) is the leading cause of end-stage renal disease worldwide and develops in 25-40% of patients with Type 1 or Type 2 diabetes mellitus. Elevated blood glucose over long periods together with glomerular hypertension leads to progressive glomerulosclerosis and tubulointerstitial fibrosis in susceptible individuals. Central to the pathology of DIN are cytokines and growth factors such as TGF-beta (transforming growth factor beta) superfamily members, including BMPs (bone morphogenetic protein) and TGF-beta 1, which play key roles in fibrogenic responses of the kidney, including podocyte loss, mesangial cell hypertrophy, matrix accumulation and tubulointerstitial fibrosis. Many of these responses can be mimicked in in vitro models of cells cultured in high glucose. We have applied differential gene expression technologies to identify novel genes expressed in in vitro and in vivo models of DN and, importantly, in human renal tissue. By mining these datasets and probing the regulation of expression and actions of specific molecules, we have identified novel roles for molecules such as Gremlin, IHG-1 (induced in high glucose-1) and CTGF (connective tissue growth factor) in DIN and potential regulators of their bioactions.

A statistical framework for the design of microarray experiments and effective detection of differential gene expression

Motivation: Microarray experiments generate a high data volume. However, often due to financial or experimental considerations, e.g. lack of sample, there is little or no replication of the experiments or hybridizations. These factors combined with the intrinsic variability associated with the measurement of gene expression can result in an unsatisfactory detection rate of differential gene expression (DGE). Our motivation was to provide an easy to use measure of the success rate of DGE detection that could find routine use in the design of microarray experiments or in post-experiment assessment.

Novel Venom Proteins Produced by Differential Domain-Expression Strategies in Beaded Lizards and Gila Monsters (genus Heloderma)

The origin and evolution of venom proteins in helodermatid lizards were investigated by multidisciplinary techniques. Our analyses elucidated novel toxin types resultant from three unique domain-expression processes: 1) The first full-length sequences of lethal toxin isoforms (helofensins) revealed this toxin type to be constructed by an ancestral monodomain, monoproduct gene (beta-defensin) that underwent three tandem domain duplications to encode a tetradomain, monoproduct with a possible novel protein fold; 2) an ancestral monodomain gene (encoding a natriuretic peptide) was medially extended to become a pentadomain, pentaproduct through the additional encoding of four tandemly repeated proline-rich peptides (helokinestatins), with the five discrete peptides liberated from each other by posttranslational proteolysis; and 3) an ancestral multidomain, multiproduct gene belonging to the vasoactive intestinal peptide (VIP)/glucagon family being mutated to encode for a monodomain, monoproduct (exendins) followed by duplication and diversification into two variant classes (exendins 1 and 2 and exendins 3 and 4). Bioactivity characterization of exendin and helokinestatin elucidated variable cardioactivity between isoforms within each class. These results highlight the importance of utilizing evolutionary-based search strategies for biodiscovery and the virtually unexplored potential of lizard venoms in drug design and discovery.

Differential Hox Expression in Murine Embryonic Stem Cell Models of Normal and Malignant Hematopoiesis.

The Hox family are master transcriptional regulators of developmental processes, including hematopoiesis. The Hox regulators, caudal homeobox factors (Cdx1-4), and Meis1, along with several individual Hox proteins, are implicated in stem cell expansion during embryonic development, with gene dosage playing a significant role in the overall function of the integrated Hox network. To investigate the role of this network in normal and aberrant, early hematopoiesis, we employed an in vitro embryonic stem cell differentiation system, which recapitulates mouse developmental hematopoiesis. Expression profiles of Hox, Pbx1, and Meis1 genes were quantified at distinct stages during the hematopoietic differentiation process and compared with the effects of expressing the leukemic oncogene Tel/PDGFRß. During normal differentiation the Hoxa cluster, Pbx1 and Meis1 predominated, with a marked reduction in the majority of Hox genes (27/39) and Meis1 occurring during hematopoietic commitment. Only the posterior Hoxa cluster genes (a9, a10, a11, and a13) maintained or increased expression at the hematopoietic colony stage. Cdx4, Meis1, and a subset of Hox genes, including a7 and a9, were differentially expressed after short-term oncogenic (Tel/PDGFRß) induction. Whereas Hoxa4-10, b1, b2, b4, and b9 were upregulated during oncogenic driven myelomonocytic differentiation. Heterodimers between Hoxa7/Hoxa9, Meis1, and Pbx have previously been implicated in regulating target genes involved in hematopoietic stem cell (HSC) expansion and leukemic progression. These results provide direct evidence that transcriptional flux through the Hox network occurs at very early stages during hematopoietic differentiation and validates embryonic stem cell models for gaining insights into the genetic regulation of normal and malignant hematopoiesis.

Parsing parallel evolution:Ecological divergence and differential gene expression in the adaptive radiations of thick-lipped Midas cichlid fishes from Nicaragua

The study of parallel evolution facilitates the discovery of common rules of diversification. Here, we examine the repeated evolution of thick lips in Midas cichlid fishes (the Amphilophus citrinellus species complex) - from two Great Lakes and two crater lakes in Nicaragua - to assess whether similar changes in ecology, phenotypic trophic traits and gene expression accompany parallel trait evolution. Using next-generation sequencing technology, we characterize transcriptome-wide differential gene expression in the lips of wild-caught sympatric thick- and thin-lipped cichlids from all four instances of repeated thick-lip evolution. Six genes (apolipoprotein D, myelin-associated glycoprotein precursor, four-and-a-half LIM domain protein 2, calpain-9, GTPase IMAP family member 8-like and one hypothetical protein) are significantly underexpressed in the thick-lipped morph across all four lakes. However, other aspects of lips' gene expression in sympatric morphs differ in a lake-specific pattern, including the magnitude of differentially expressed genes (97-510). Generally, fewer genes are differentially expressed among morphs in the younger crater lakes than in those from the older Great Lakes. Body shape, lower pharyngeal jaw size and shape, and stable isotopes (dC and dN) differ between all sympatric morphs, with the greatest differentiation in the Great Lake Nicaragua. Some ecological traits evolve in parallel (those related to foraging ecology; e.g. lip size, body and head shape) but others, somewhat surprisingly, do not (those related to diet and food processing; e.g. jaw size and shape, stable isotopes). Taken together, this case of parallelism among thick- and thin-lipped cichlids shows a mosaic pattern of parallel and nonparallel evolution. © 2012 Blackwell Publishing Ltd.

Association of gene expression with sequential proliferation, differentiation and tumor formation in murine skin

Differential gene expression in two established initiation and promotion skin carcinogenesis models during promotion and tumor formation was determined by microarray technology with the purpose of distinguishing the genes more associated with neoplastic transformation from those linked with proliferation and differentiation. The first model utilized dimethylbenz[a]anthracene initiation and 12-O-tetradecanoylphorbol 13-acetate (TPA) promotion in the FVB/N mouse, and the second TPA promotion of the Tg.Ac mouse, which is endogenously initiated by virtue of an activated Ha-ras transgene. Comparison of gene expression profiles across the two models identified genes whose altered expression was associated with papilloma formation rather than TPA-induced proliferation and differentiation. DMBA suppressed TPA-induced differentiation which allowed identification of those genes associated more specifically with differentiation rather than proliferation. EASE (Expression Analysis Systemic Explorer) indicated a correlation between muscle-associated genes and skin differentiation, whereas genes involved with protein biosynthesis were strongly correlated with proliferation. For verification the altered expression of selected genes were confirmed by RT-PCR; Carbonic anhydrase 2, Thioredoxin 1 and Glutathione S-transferase omega 1 associated with papilloma formation and Enolase 3, Cystatin 6 and Filaggrin associated with TPA-induced proliferation and differentiation. In situ analysis located the papillomas Glutathione S-transferase omega 1 expression to the proliferating areas of the papillomas. Thus we have identified profiles of differential gene expression associated with the tumorigenesis and promotion stages for skin carcinogenesis in the mouse.

cudaMap: a GPU accelerated program for gene expression connectivity mapping

Background: Modern cancer research often involves large datasets and the use of sophisticated statistical techniques. Together these add a heavy computational load to the analysis, which is often coupled with issues surrounding data accessibility. Connectivity mapping is an advanced bioinformatic and computational technique dedicated to therapeutics discovery and drug re-purposing around differential gene expression analysis. On a normal desktop PC, it is common for the connectivity mapping task with a single gene signature to take >2h to complete using sscMap, a popular Java application that runs on standard CPUs (Central Processing Units). Here, we describe new software, cudaMap, which has been implemented using CUDA C/C++ to harness the computational power of NVIDIA GPUs (Graphics Processing Units) to greatly reduce processing times for connectivity mapping.

Results: cudaMap can identify candidate therapeutics from the same signature in just over thirty seconds when using an NVIDIA Tesla C2050 GPU. Results from the analysis of multiple gene signatures, which would previously have taken several days, can now be obtained in as little as 10 minutes, greatly facilitating candidate therapeutics discovery with high throughput. We are able to demonstrate dramatic speed differentials between GPU assisted performance and CPU executions as the computational load increases for high accuracy evaluation of statistical significance.

Conclusion: Emerging 'omics' technologies are constantly increasing the volume of data and information to be processed in all areas of biomedical research. Embracing the multicore functionality of GPUs represents a major avenue of local accelerated computing. cudaMap will make a strong contribution in the discovery of candidate therapeutics by enabling speedy execution of heavy duty connectivity mapping tasks, which are increasingly required in modern cancer research. cudaMap is open source and can be freely downloaded from http://purl.oclc.org/NET/cudaMap.

CAMP factor homologues in Propionibacterium acnes: a new protein family differentially expressed by types I and II

Analysis of the draft genome sequence of the opportunistic pathogen Propionibacterium acnes type strain NCTC 737 (=ATCC 6919) revealed five genes with sequence identity to the co-haemolytic Christie-Atkins-Munch-Peterson (CAMP) factor of Streptococcus agalactiae. The predicted molecular masses for the expressed proteins ranged from 28 to 30 kDa. The genes were present in each of the three recently identified recA-based phylogenetic groupings of P. acnes (IA, IB and 11), as assessed by PCR amplification. Conserved differences in CAMP factor gene sequences between these three groups were also consistent with their previous phylogenetic designations. All type IA, IB and 11 isolates were positive for the co-haemolytic; reaction on sheep blood agar. Immunoblotting and silver staining of SIDS-PAGE gels, however, revealed differential protein expression of CAMP factors amongst the different groups. Type IB and 11 isolates produced an abundance of CAMP factor 1, detectable by specific antibody labelling and silver staining of SDS-PAGE gels. In contrast, abundant CAMP factor production was lacking in type A isolates, although larger amounts of CAMP factor 2 were detectable by immunoblotting compared with type 11 isolates. While the potential role of the abundant CAMP factor 1 in host colonization or virulence remains to be determined, it should be noted that the type strain of P. acnes used in much of the published literature is a type A isolate and is, therefore, lacking in this attribute.

Hypoxia selects for androgen independent LNCaP cells with a more malignant geno- and phenotype.

Hypoxia confers resistance to common cancer therapies, however, it has also has been shown to result in genetic alterations which may allow a survival advantage and increase the tumorigenic properties of cancer cells. Additionally, it may exert a selection pressure, allowing expansion of tumor cells with a more aggressive phenotype. To further assess the role of hypoxia in malignant progression in prostate cancer we exposed human androgen dependent prostate cancer cells (LNCaP) to cycles of chronic hypoxia and isolated a subline, LNCaP-H1. This article describes the partial characterization of this cell line. The LNCaP-H1 subline showed altered growth characteristics and exhibited androgen independent growth both in vitro and in vivo. Furthermore, these cells were resistant to mitochondrial-mediated apoptosis, probably since the endogenous levels of Bax was lower and Bcl-2 higher than in the parental LNCaP cells. Microarray analysis revealed that a complex array of pathways had differential gene expression between the 2 cell lines, with LNCaP-H1 cells exhibiting a genetic profile which suggests that they may be more likely metastasize to distant organs, especially bone. This was supported by an in vitro invasion assay, and an in vivo metastasis study. This study shows that hypoxia can select for androgen independent prostate cancer cells which have a survival advantage and are more likely to invade and metastasize.

Androgen deprivation results in time-dependent hypoxia in LNCaP prostate tumours; informed scheduling of the bioreductive drug AQ4N improves treatment response

Androgen withdrawal induces hypoxia in androgen-sensitive tissue; this is important as in the tumour microenvironment hypoxia is known to drive malignant progression. This study examined the time-dependent effect of androgen deprivation therapy (ADT) on tumour oxygenation and investigated the role of ADT-induced hypoxia on malignant progression in prostate tumours. LNCaP xenografted tumours were treated with anti-androgens and tumour oxygenation measured. Dorsal skin fold chambers (DSF) were used to image tumour vasculature in vivo. Quantitative PCR (QPCR) identified differential gene expression following treatment with bicalutamide. Bicalutamide and vehicle-only treated tumours were re-established in vitro and invasion and sensitivity to docetaxel were measured. Tumour growth delay was calculated following treatment with bicalutamide combined with the bioreductive drug AQ4N. Tumour oxygenation measurements showed a precipitate decrease following initiation of ADT. A clinically relevant dose of bicalutamide (2mg/kg/day) decreased tumour oxygenation by 45% within 24h, reaching a nadir of 0.09% oxygen (0.67±0.06 mmHg) by day 7; this persisted until day 14 when it increased up to day 28. Using DSF chambers, LNCaP tumours treated with bicalutamide showed loss of small vessels at days 7 and 14 with revascularization occurring by day 21. QPCR showed changes in gene expression consistent with the vascular changes and malignant progression. Cells from bicalutamide-treated tumours were more malignant than vehicle-treated controls. Combining bicalutamide with AQ4N (50mg/kg; single dose) caused greater tumour growth delay than bicalutamide alone. This study shows that bicalutamide-induced hypoxia selects for cells that show malignant progression; targeting hypoxic cells may provide greater clinical benefit.

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.

Novel Hypoxia-Associated Markers of Chemoresistance in Ovarian Cancer

Background:
Ovarian cancer is the fifth leading cause of cancer in women and has poor
long-term survival, in part, due to chemoresistance. Tumour hypoxia is associated with
chemoresistance in ovarian cancer. However, relatively little is known about the genes
activated in ovarian cancer which cause chemoresistance due to hypoxia. This study
aimed to firstly identify genes whose expression is associated with hypoxia-induced
chemoresistance, and secondly select hypoxia-associated biomarkers and evaluate their
expression in ovarian tumours.
Design:
Cisplatin-sensitive (A2780) and cisplatin-resistant (A2780cis) ovarian cancer
cell lines were exposed to combinations of hypoxia and/or cisplatin as part of a matrix
designed to reflect clinically relevant scenarios. RNA was extracted and interrogated
on Affymetrix Human Gene arrays. Differential gene expression was analysed for cells
exposed to hypoxia and/or treated with cisplatin. Potential markers of chemoresistance
were selected for evaluation in a cohort of ovarian tumour samples by R
T-PCR.
Results:
A wide range of genes associated with chemoresistance were differentially
expressed in cells exposed to hypoxia and/or cisplatin. Selected genes [ANGPTL4,
HER3 and HIF-1
α
] were chosen for further validation in a cohort of ovarian tumour
samples. High expression of ANGPTL4 trended towards reduced progression-free and
overall survival. High expression of HER3 trended to increased progression-free but
reduced overall survival, while high expression of HIF-1
α
trended towards reduced
progression-free and increased overall survival.
Conclusions:
In conclusion, this study has further characterized the relationship between
hypoxia and chemoresistance in an ovarian cancer model. We have also identified many
potential biomarkers of hypoxia and platinum resistance and provided initial validation
of a subset of these markers in ovarian cancer tissues.

Androgen receptor driven transcription in molecular apocrine breast cancer is mediated by FoxA1

Breast cancer is a heterogeneous disease and several distinct subtypes exist based on differential gene expression patterns. Molecular apocrine tumours were recently identified as an additional subgroup, characterised as oestrogen receptor negative and androgen receptor positive (ER- AR+), but with an expression profile resembling ER+ luminal breast cancer. One possible explanation for the apparent incongruity is that ER gene expression programmes could be recapitulated by AR. Using a cell line model of ER- AR+ molecular apocrine tumours (termed MDA-MB-453 cells), we map global AR binding events and find a binding profile that is similar to ER binding in breast cancer cells. We find that AR binding is a near-perfect subset of FoxA1 binding regions, a level of concordance never previously seen with a nuclear receptor. AR functionality is dependent on FoxA1, since silencing of FoxA1 inhibits AR binding, expression of the majority of the molecular apocrine gene signature and growth cell growth. These findings show that AR binds and regulates ER cis-regulatory elements in molecular apocrine tumours, resulting in a transcriptional programme reminiscent of ER-mediated transcription in luminal breast cancers.

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.

Proteomic analysis of coronary sinus serum reveals leucine-rich α2-glycoprotein as a novel biomarker of ventricular dysfunction and heart failure

BACKGROUND: Heart failure (HF) prevention strategies require biomarkers that identify disease manifestation. Increases in B-type natriuretic peptide (BNP) correlate with increased risk of cardiovascular events and HF development. We hypothesize that coronary sinus serum from a high BNP hypertensive population reflects an active pathological process and can be used for biomarker exploration. Our aim was to discover differentially expressed disease-associated proteins that identify patients with ventricular dysfunction and HF.

METHODS AND RESULTS: Coronary sinus serum from 11 asymptomatic, hypertensive patients underwent quantitative differential protein expression analysis by 2-dimensional difference gel electrophoresis. Proteins were identified using mass spectrometry and then studied by enzyme-linked immunosorbent assay in sera from 40 asymptomatic, hypertensive patients and 105 patients across the spectrum of ventricular dysfunction (32 asymptomatic left ventricular diastolic dysfunction, 26 diastolic HF, and 47 systolic HF patients). Leucine-rich α2-glycoprotein (LRG) was consistently overexpressed in high BNP serum. LRG levels correlate significantly with BNP in hypertensive, asymptomatic left ventricular diastolic dysfunction, diastolic HF, and systolic HF patient groups (P≤0.05). LRG levels were able to identify HF independent of BNP. LRG correlates with coronary sinus serum levels of tumor necrosis factor-α (P=0.009) and interleukin-6 (P=0.021). LRG is expressed in myocardial tissue and correlates with transforming growth factor-βR1 (P<0.001) and α-smooth muscle actin (P=0.025) expression.

CONCLUSIONS: LRG was identified as a serum biomarker that accurately identifies patients with HF. Multivariable modeling confirmed that LRG is a stronger identifier of HF than BNP and this is independent of age, sex, creatinine, ischemia, β-blocker therapy, and BNP.