Large scale real-time PCR validation on gene expression measurements from two commercial long-oligonucleotide microarrays

A transcriptomal analysis of bovine oviductal epithelial cells collected during the follicular phase versus the luteal phase of the estrous cycle

BACKGROUND

Reproductive success depends on a functional oviduct for gamete storage, maturation, fertilization, and early embryonic development. The ovarian-derived steroids estrogen and progesterone are key regulators of oviductal function. The objective of this study was to investigate luteal and follicular phase-specific oviductal epithelial cell function by using microarray-based transcriptional profiling, to increase our understanding of mRNAs regulating epithelial cell processes, and to identify novel genes and biochemical pathways that may be found to affect fertility in the future.

METHODS

Six normally cycling Angus heifers were assigned to either luteal phase (LP, n = 3) or follicular phase (FP, n = 3) treatment groups. Heifers in the LP group were killed between day 11 and 12 after estrus. Heifers in the FP group were treated with 25 mg PGF2α (Lutalyse, Pfizer, NY) at 8 pm on day 6 after estrus and killed 36 h later. Transcriptional profiling by microarray and confirmation of selected mRNAs by real-time RT-PCR analyses was performed using total RNA from epithelial cells isolated from sections of the ampulla and isthmus collected from LP and FP treatment groups. Differentially expressed genes were subjected to gene ontology classification and bioinformatic pathway analyses.

RESULTS

Statistical one-way ANOVA using Benjamini-hochberg multiple testing correction for false discovery rate (FDR) and pairwise comparison of epithelial cells in the ampulla of FP versus LP groups revealed 972 and 597 transcripts up- and down-regulated, respectively (P < 0.05). Within epithelial cells of the isthmus in FP versus LP groups, 946 and 817 transcripts were up- and down-regulated, respectively (P < 0.05). Up-regulated genes from both ampulla and isthmus were found to be largely involved in cholesterol biosynthesis and cell cycle pathways, while down-regulated genes were found in numerous inflammatory response pathways.

CONCLUSIONS

Microarray-based transcriptional profiling revealed phase of the cycle-dependent changes in the expression of mRNA within the epithelium of the oviducts' ampulla and isthmus.

Piwi-interacting RNAs as novel prognostic markers in clear cell renal cell carcinomas

Background

Piwi-interacting RNAs (piRNAs) are small RNAs of 27–30 nucleotides mapping to transposons or clustering in repeat genomic regions. Preliminary studies suggest an important role in cancerogenesis. This study is the first one investigating their prognostic impact in clear cell renal cell cancer (ccRCC) patients.

Methods

Three piRNAs (piR-30924, piR-57125, and piR-38756) selected on the basis of initial piRNA microarray analyses were determined using RT-qPCR in non-metastatic (n = 76) and metastatic (n = 30) ccRCC tissue at the time of nephrectomy in comparison to normal renal tissue (n = 77) and tissue from distant ccRCC metastases (n = 13). Primary clinical end points were recurrence-free and overall survival.

Results

piR-57125 showed lower expression in metastatic than in non-metastatic tumors, whereas the expression of piR-30924 and piR-38756 increased in metastatic tumors. The higher expression of piR-30924 and piR-38756 as well as the lower expression of piR-57125 in metastatic primary tumors were significantly associated with tumor recurrence and overall survival. Multivariate Cox regression analyses revealed both piR-30924 and piR-57125 as independent prognostic predictors. This impact was even more pronounced in non-metastatic patients.

Conclusions

This study demonstrates that the expression levels of these piRNAs in primary non-metastatic and metastatic ccRCC tissue can serve as potential prognostic biomarkers in combination with clinicopathological factors.

Electronic supplementary material

The online version of this article (doi:10.1186/s13046-015-0180-3) contains supplementary material, which is available to authorized users.

Altered expression of Arabidopsis genes in response to a multifunctional geminivirus pathogenicity protein

BACKGROUND

Geminivirus AC2 is a multifunctional protein that acts as a pathogenicity factor. Transcriptional regulation by AC2 appears to be mediated through interaction with a plant specific DNA binding protein, PEAPOD2 (PPD2), that specifically binds to sequences known to mediate activation of the CP promoter of Cabbage leaf curl virus (CaLCuV) and Tomato golden mosaic virus (TGMV). Suppression of both basal and innate immune responses by AC2 in plants is mediated through inactivation of SnRK1.2, an Arabidopsis SNF1 related protein kinase, and adenosine kinase (ADK). An indirect promoter targeting strategy, via AC2-host dsDNA binding protein interactions, and inactivation of SnRK1.2-mediated defense responses could provide the opportunity for geminiviruses to alter host gene expression and in turn, reprogram the host to support virus infection. The goal of this study was to identify changes in the transcriptome of Arabidopsis induced by the transcription activation function of AC2 and the inactivation of SnRK1.2.

RESULTS

Using full-length and truncated AC2 proteins, microarray analyses identified 834 genes differentially expressed in response to the transcriptional regulatory function of the AC2 protein at one and two days post treatment. We also identified 499 genes differentially expressed in response to inactivation of SnRK1.2 by the AC2 protein at one and two days post treatment. Network analysis of these two sets of differentially regulated genes identified several networks consisting of between four and eight highly connected genes. Quantitative real-time PCR analysis validated the microarray expression results for 10 out of 11 genes tested.

CONCLUSIONS

It is becoming increasingly apparent that geminiviruses manipulate the host in several ways to facilitate an environment conducive to infection, predominantly through the use of multifunctional proteins. Our approach of identifying networks of highly connected genes that are potentially co-regulated by geminiviruses during infection will allow us to identify novel pathways of co-regulated genes that are stimulated in response to pathogen infection in general, and virus infection in particular.

Analysis of the arabidopsis REM gene family predicts functions during flower development

BACKGROUND AND AIMS

The REM (Reproductive Meristem) gene family of Arabidopsis thaliana is part of the B3 DNA-binding domain superfamily. Despite the fact that several groups have worked on the REM genes for many years, little is known about the function of this transcription factor family. This study aims to identify a set of REM genes involved in flower development and to characterize their function.

METHODS

In order to provide an overview of the REM gene family, a detailed expression analysis for all REM genes of A. thaliana was performed and combined with a meta-analysis of ChIP-sequencing and microarray experiments.

KEY RESULTS

Two sets of phylogenetically closely related REM genes, namely REM23, REM24 and REM25, and REM34, REM35 and REM36, were identified as possibly being involved in the early stages of flower development. Single- and double-mutant combinations were analysed for these genes, and no phenotypic effects were detected during flower development.

CONCLUSIONS

The data suggest that the REM34, REM35 and REM36 group is the most interesting one, as REM34 is co-expressed with the floral meristem identity (FMI) genes, they are bound by AP1, SVP, AP3 and PI, and they are expressed in the floral meristem and during the earliest stages of flower development. However, it appears that high levels of functional redundancy may conceal the exact function of these transcription factor genes.

Comparative functional genomics of adaptation to muscular disuse in hibernating mammals

Hibernation is an energy-saving adaptation that involves a profound suppression of physical activity that can continue for 6-8 months in highly seasonal environments. While immobility and disuse generate muscle loss in most mammalian species, in contrast, hibernating bears and ground squirrels demonstrate limited muscle atrophy over the prolonged periods of physical inactivity during winter, suggesting that hibernating mammals have adaptive mechanisms to prevent disuse muscle atrophy. To identify common transcriptional programmes that underlie molecular mechanisms preventing muscle loss, we conducted a large-scale gene expression screen in hind limb muscles comparing hibernating and summer-active black bears and arctic ground squirrels using custom 9600 probe cDNA microarrays. A molecular pathway analysis showed an elevated proportion of overexpressed genes involved in all stages of protein biosynthesis and ribosome biogenesis in muscle of both species during torpor of hibernation that suggests induction of translation at different hibernation states. The induction of protein biosynthesis probably contributes to attenuation of disuse muscle atrophy through the prolonged periods of immobility of hibernation. The lack of directional changes in genes of protein catabolic pathways does not support the importance of metabolic suppression for preserving muscle mass during winter. Coordinated reduction in multiple genes involved in oxidation-reduction and glucose metabolism detected in both species is consistent with metabolic suppression and lower energy demand in skeletal muscle during inactivity of hibernation.

Neurotoxin-induced pathway perturbation in human neuroblastoma SH-EP cells

The exact causes of cell death in Parkinson's disease (PD) remain unknown despite extensive studies on PD.The identification of signaling and metabolic pathways involved in PD might provide insight into the molecular mechanisms underlying PD. The neurotoxin 1-methyl-4-phenylpyridinium (MPP(+)) induces cellular changes characteristic of PD, and MPP(+)-based models have been extensively used for PD studies. In this study, pathways that were significantly perturbed in MPP(+)-treated human neuroblastoma SH-EP cells were identified from genome-wide gene expression data for five time points (1.5, 3, 9, 12, and 24 h) after treatment. The mitogen-activated protein kinase (MAPK) signaling pathway and endoplasmic reticulum (ER) protein processing pathway showed significant perturbation at all time points. Perturbation of each of these pathways resulted in the common outcome of upregulation of DNA-damage-inducible transcript 3 (DDIT3). Genes involved in ER protein processing pathway included ubiquitin ligase complex genes and ER-associated degradation (ERAD)-related genes. Additionally, overexpression of DDIT3 might induce oxidative stress via glutathione depletion as a result of overexpression of CHAC1. This study suggests that upregulation of DDIT3 caused by perturbation of the MAPK signaling pathway and ER protein processing pathway might play a key role in MPP(+)-induced neuronal cell death. Moreover, the toxicity signal of MPP(+) resulting from mitochondrial dysfunction through inhibition of complex I of the electron transport chain might feed back to the mitochondria via ER stress. This positive feedback could contribute to amplification of the death signal induced by MPP(+).

Assessing reference genes for accurate transcript normalization using quantitative real-time PCR in pearl millet [Pennisetum glaucum (L.) R. Br]

Pearl millet [Pennisetum glaucum (L.) R.Br.], a close relative of Panicoideae food crops and bioenergy grasses, offers an ideal system to perform functional genomics studies related to C4 photosynthesis and abiotic stress tolerance. Quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) provides a sensitive platform to conduct such gene expression analyses. However, the lack of suitable internal control reference genes for accurate transcript normalization during qRT-PCR analysis in pearl millet is the major limitation. Here, we conducted a comprehensive assessment of 18 reference genes on 234 samples which included an array of different developmental tissues, hormone treatments and abiotic stress conditions from three genotypes to determine appropriate reference genes for accurate normalization of qRT-PCR data. Analyses of Ct values using Stability Index, BestKeeper, ΔCt, Normfinder, geNorm and RefFinder programs ranked PP2A, TIP41, UBC2, UBQ5 and ACT as the most reliable reference genes for accurate transcript normalization under different experimental conditions. Furthermore, we validated the specificity of these genes for precise quantification of relative gene expression and provided evidence that a combination of the best reference genes are required to obtain optimal expression patterns for both endogeneous genes as well as transgenes in pearl millet.

Transcriptional behavior of EUL-related rice lectins toward important abiotic and biotic stresses

The rice genome encodes several genes for putative carbohydrate-binding proteins belonging to the family of Euonymus related lectins (EULs). This lectin family was discovered recently and evidence shows that the expression of these proteins is subject to multiple environmental stresses. In this study, quantitative reverse transcription PCR (qRT-PCR) was conducted on rice seedlings exposed to various abiotic (150mM NaCl, 100mM mannitol, and 100μM abscisic acid (ABA)) and biotic (Xanthomonas oryzae pv. oryzae and Magnaporthe oryzae) stresses to compare the transcriptional behavior of the EULs and a known stress related lectin Orysata belonging to the family of jacalin-related lectins. All EUL transcripts were strongly up-regulated after ABA and NaCl treatments in the roots whereas the overall expression level was generally lower and more variable in the shoots. Moreover, all abiotic stresses induced Orysata in both tissues except for mannitol treatment which failed to show an effect in the roots. Orysata also strongly accumulated after X. oryzae pv. oryzae infection, as were various D-type EUL lectins. In contrast, some of the EUL proteins, including OrysaEULS3, OrysaEULD1A and OrysaEULD2, as well as Orysata were significantly down-regulated upon M. oryzae attack, suggesting fungal manipulation of these genes. Collectively, our results clearly show that rice expresses multiple carbohydrate-binding proteins in response to a wide variety of abiotic and biotic stress conditions. We hypothesize that the Euonymus related proteins fulfill a prominent role in sensing and responding to multiple environmental cues.

Identification of gene-based responses in human blood cells exposed to alpha particle radiation

BACKGROUND

The threat of a terrorist-precipitated nuclear event places humans at danger for radiological exposures. Isotopes which emit alpha (α)-particle radiation pose the highest risk. Currently, gene expression signatures are being developed for radiation biodosimetry and triage with respect to ionizing photon radiation. This study was designed to determine if similar gene expression profiles are obtained after exposures involving α-particles.

METHODS

Peripheral blood mononuclear cells (PBMCs) were used to identify sensitive and robust gene-based biomarkers of α-particle radiation exposure. Cells were isolated from healthy individuals and were irradiated at doses ranging from 0-1.5 Gy. Microarray technology was employed to identify transcripts that were differentially expressed relative to unirradiated cells 24 hours post-exposure. Statistical analysis identified modulated genes at each of the individual doses.

RESULTS

Twenty-nine genes were common to all doses with expression levels ranging from 2-10 fold relative to control treatment group. This subset of genes was further assessed in independent complete white blood cell (WBC) populations exposed to either α-particles or X-rays using quantitative real-time PCR. This 29 gene panel was responsive in the α-particle exposed WBCs and was shown to exhibit differential fold-changes compared to X-irradiated cells, though no α-particle specific transcripts were identified.

CONCLUSION

Current gene panels for photon radiation may also be applicable for use in α-particle radiation biodosimetry.

Overexpression of MAGE-D4 in colorectal cancer is a potentially prognostic biomarker and immunotherapy target

Melanoma-associated antigen D4 (MAGE-D4) is a novel member of MAGE family. This study aimed to examine the expression and immunogenicity of MAGE-D4 in colorectal cancer (CRC) to determine its potential as a prognosis and immunotherapeutic target. The expression of MAGE-D4 mRNA and protein was determined by RT-PCR and immunohistochemistry (IHC) in CRCs with paired adjacent non-tumor tissues, colorectal adenomas and normal colorectal tissues, respectively. Sera from 64 CRC patients were tested for MAGE-D4 antibody by ELISA. MAGE-D4 mRNA was more frequently expressed in CRCs (76.7%, 46/60) than in adjacent non-tumor tissues (15.0%, 9/60). MAGE-D4 protein was detected in all the CRC tissues tested, 70.0% of which showed high expression. There was no MAGE-D4 protein detected in any paired adjacent non-tumor tissue. No MAGE-D4 expression was found in colorectal adenomas and normal colorectal tissues by either RT-PCR or immunohistochemistry. Patients with high MAGE-D4 protein expression had significantly shorter overall survival than those with low MAGE-D4 protein expression (median, 68.6 vs 122.2 months; P=0.030). Furthermore, multivariate analysis exhibited high MAGE-D4 protein expression had a trend toward an independent prognostic factor (hazard ratio: 6.124; P=0.050). Humoral immunity to MAGE-D4 was detected in 12 of 64 (18.8%) CRC patients' sera but not in 77 healthy donors. There was no correlation between MAGE-D4 expression, serum antibody and clinicopathological parameters. These findings suggest MAGE-D4 may serve as a potentially prognostic biomarker and an attractive target of immunotherapy in CRC.

A detailed analysis of next generation sequencing reads of microRNA expression in Barrett's esophagus: absolute versus relative quantification

BACKGROUND

Next generation sequencing (NGS) is a state of the art technology for microRNA (miRNA) analysis. The quantitative interpretation of the primary output of NGS i.e. the read counts for a miRNA sequence that can vary by several orders of magnitude (1 to 107) remains incompletely understood.

FINDINGS

NGS (SOLiD 3 technology) was performed on biopsies from 6 Barrett's esophagus (BE) and 5 Gastroesophageal Reflux Disease (GERD) patients. Read sequences were aligned to miRBase 18.0. Differential expression analysis was adjusted for false discovery rate of 5%. Quantitative real-time polymerase chain reaction (qRT-PCR) was performed for 36 miRNA in a validation cohort of 47 patients (27 BE and 20 GERD). Correlation coefficients, accuracy, precision and recall of NGS compared to qRT-PCR were calculated. Increase in NGS reads was associated with progressively lower Cq values, p < 0.05. Although absolute quantification between NGS reads and Cq values correlated modestly: -0.38, p = 0.01 for BE and -0.32, p = 0.05 for GERD, relative quantification (fold changes) of miRNA expression between BE &GERD by NGS correlated highly with qRT-PCR 0.86, p = 2.45E-11. Fold change correlations were unaffected when different thresholds of NGS read counts were compared (>1000 vs. <1000, >500 vs. <500 and >100 vs. <100). The accuracy, precision and recall of NGS to label a miRNA as differentially expressed were 0.71, 0.88 and 0.74 respectively.

CONCLUSION

Absolute NGS reads correlated modestly with qRT-PCR but fold changes correlated highly. NGS is robust at relative but not absolute quantification of miRNA levels and accurate for high-throughput identification of differentially expressed miRNA.

mRNA-Seq and microarray development for the Grooved Carpet shell clam, Ruditapes decussatus: a functional approach to unravel host-parasite interaction

Background

The Grooved Carpet shell clam Ruditapes decussatus is the autochthonous European clam and the most appreciated from a gastronomic and economic point of view. The production is in decline due to several factors such as Perkinsiosis and habitat invasion and competition by the introduced exotic species, the manila clam Ruditapes philippinarum. After we sequenced R. decussatus transcriptome we have designed an oligo microarray capable of contributing to provide some clues on molecular response of the clam to Perkinsiosis.

Results

A database consisting of 41,119 unique transcripts was constructed, of which 12,479 (30.3%) were annotated by similarity. An oligo-DNA microarray platform was then designed and applied to profile gene expression in R. decussatus heavily infected by Perkinsus olseni. Functional annotation of differentially expressed genes between those two conditionswas performed by gene set enrichment analysis. As expected, microarrays unveil genes related with stress/infectious agents such as hydrolases, proteases and others. The extensive role of innate immune system was also analyzed and effect of parasitosis upon expression of important molecules such as lectins reviewed.

Conclusions

This study represents a first attempt to characterize Ruditapes decussatus transcriptome, an important marine resource for the European aquaculture. The trancriptome sequencing and consequent annotation will increase the available tools and resources for this specie, introducing the possibility of high throughput experiments such as microarrays analysis. In this specific case microarray approach was used to unveil some important aspects of host-parasite interaction between the Carpet shell clam and Perkinsus, two non-model species, highlighting some genes associated with this interaction. Ample information was obtained to identify biological processes significantly enriched among differentially expressed genes in Perkinsus infected versus non-infected gills. An overview on the genes related with the immune system on R. decussatus transcriptome is also reported.

Comparative transcriptional profiling of 3 murine models of SLE nephritis reveals both unique and shared regulatory networks

Objective

To define shared and unique features of SLE nephritis in mouse models of proliferative and glomerulosclerotic renal disease.

Methods

Perfused kidneys from NZB/W F1, NZW/BXSB and NZM2410 mice were harvested before and after nephritis onset. Affymetrix based gene expression profiles of kidney RNA were analyzed using Genomatix Pathway Systems and Ingenuity Pathway Analysis software. Gene expression patterns were confirmed using real-time PCR.

Results

955, 1168 and 755 genes were regulated in the kidneys of nephritic NZB/W F1, NZM2410 and NZW/BXSB mice respectively. 263 genes were regulated concordantly in all three strains reflecting immune cell infiltration, endothelial cell activation, complement activation, cytokine signaling, tissue remodeling and hypoxia. STAT3 was the top associated transcription factor, having a binding site in the gene promoter of 60/263 regulated genes. The two strains with proliferative nephritis shared a macrophage/DC infiltration and activation signature. NZB/W and NZM2410 mice shared a mitochondrial dysfunction signature. Dominant T cell and plasma cell signatures in NZB/W mice reflected lymphoid aggregates; this was the only strain with regulatory T cell infiltrates. NZW/BXSB mice manifested tubular regeneration and NZM2410 mice had the most metabolic stress and manifested loss of nephrin, indicating podocyte loss.

Conclusions

These findings identify shared inflammatory mechanisms of SLE nephritis that can be therapeutically targeted. Nevertheless, the heterogeneity of effector mechanisms suggests that individualized therapy might need to be based on biopsy findings. Some common mechanisms are shared with non-immune–mediated renal diseases, suggesting that strategies to prevent tissue hypoxia and remodeling may be useful in SLE nephritis.

Capturing drug responses by quantitative promoter activity profiling

Quantitative analysis of cellular responses to drugs is of major interest in pharmaceutical research. Microarray technologies have been widely used for monitoring genome-wide expression changes. However, this approach has several limitations in terms of coverage of targeted RNAs, sensitivity, and quantitativeness, which are crucial for accurate monitoring of cellular responses. In this article, we report an application of genome-wide and quantitative profiling of cellular responses to drugs. We monitored promoter activities in MCF-7 cells by Cap Analysis of Gene Expression using a single-molecule sequencer. We identified a distinct set of promoters affected even by subtle inhibition of the Ras-ERK and phosphatidylinositol-3-kinase-Akt signal-transduction pathways. Furthermore, we succeeded in explaining the majority of promoter responses to inhibition of the upstream epidermal growth factor receptor kinase quantitatively based on the promoter profiles upon inhibition of the two individual downstream signaling pathways. Our results demonstrate unexplored utility of highly quantitative promoter activity profiling in drug research.

Tpl2 knockout keratinocytes have increased biomarkers for invasion and metastasis

Skin cancer is the most common form of cancer in the USA, with an estimated two million cases diagnosed annually. Tumor progression locus 2 (Tpl2), also known as MAP3K8, is a serine/threonine protein kinase in the mitogen-activated protein kinase signal transduction cascade. Tpl2 was identified by our laboratory as having a tumor suppressor function in skin carcinogenesis, with the absence of this gene contributing to heightened inflammation and increased skin carcinogenesis. In this study, we used gene expression profiling to compare expression levels between Tpl2 (+/+) and Tpl2 (-) (/-) keratinocytes. We identified over 2000 genes as being differentially expressed between genotypes. Functional annotation analysis identified cancer, cell growth/proliferation, cell death, cell development, cell movement and cell signaling as the top biological processes to be differentially regulated between genotypes. Further microarray analysis identified several candidate genes, including Mmp1b, Mmp2, Mmp9 and Mmp13, involved in migration and invasion to be upregulated in Tpl2 (-) (/-) keratinocytes. Moreover, Tpl2 (-/-) keratinocytes had a significant downregulation in the matrix metalloproteinase (MMP) inhibitor Timp3. Real-time PCR validated the upregulation of the MMPs in Tpl2 (-/-) keratinocytes and zymography confirmed that MMP2 and MMP9 activity was higher in conditioned media from Tpl2 (-/-) keratinocytes. Immunohistochemistry confirmed higher MMP9 staining in 12-O-tetradecanoylphorbol-13-acetate-treated skin from Tpl2 (-/-) mice and grafted tumors formed from v-ras(Ha) retrovirus-infected Tpl2 (-/-) keratinocytes. Additionally, Tpl2 (-/-) keratinocytes had significantly higher invasion, malignant conversion rates and increased endothelial cell tube formation when compared with Tpl2 (+/+) keratinocytes. In summary, our studies reveal that keratinocytes from Tpl2 (-/-) mice demonstrate a higher potential to be invasive and metastatic.

Physiological and Molecular Timing of the Glucose to Acetate Transition in Escherichia coli

The glucose-acetate transition in Escherichia coli is a classical model of metabolic adaptation. Here, we describe the dynamics of the molecular processes involved in this metabolic transition, with a particular focus on glucose exhaustion. Although changes in the metabolome were observed before glucose exhaustion, our results point to a massive reshuffling at both the transcriptome and metabolome levels in the very first min following glucose exhaustion. A new transcriptional pattern, involving a change in genome expression in one-sixth of the E. coli genome, was established within 10 min and remained stable until the acetate was completely consumed. Changes in the metabolome took longer and stabilized 40 min after glucose exhaustion. Integration of multi-omics data revealed different modifications and timescales between the transcriptome and metabolome, but both point to a rapid adaptation of less than an hour. This work provides detailed information on the order, timing and extent of the molecular and physiological events that occur during the glucose-acetate transition and that are of particular interest for the development of dynamic models of metabolism.

Comparative analysis of gene transcripts for cell signaling receptors in bone marrow-derived hematopoietic stem/progenitor cell and mesenchymal stromal cell populations

Introduction

Knowing the repertoire of cell signaling receptors would provide pivotal insight into the developmental and regenerative capabilities of bone marrow cell (BMC)-derived hematopoietic stem/progenitor cells (HSPCs) and bone marrow mesenchymal stromal cells (BMMSCs).

Methods

Murine HSPCs were enriched from fluorescence-activated cell sorting (FACS)-sorted Lin^–c-Kit⁺Sca-1⁺ BMCs isolated from the tibia and femoral marrow compartments. Purified BMMSCs (CD73⁺, CD90⁺, CD105⁺, and CD45^–, CD34^–, CD31^–, c-Kit^–) with extensive self-renewal potential and multilineage differentiation capacity (into different mesodermal cell lineages including osteocytes, chrondrocytes, adipocytes) were derived from adherent BMC cultures after CD45⁺ cell depletion. Adherent colony-forming cells were passaged two to three times and FACS analysis was used to assess cell purity and validate cell-specific surface marker phenotype prior to experimentation. Gene transcripts for a number of cell signaling molecules were assessed using a custom quantitative real-time RT-PCR low-density microarray (94 genes; TaqMan^® technology).

Results

We identified 16 mRNA transcripts that were specifically expressed in BMC-derived HSPC (including Ptprc, c-Kit, Csf3r, Csf2rb2, Ccr4, Cxcr3 and Tie-1), and 14 transcripts specifically expressed in BMMSCs (including Pdgfra, Ddr2, Ngfr, Mst1r, Fgfr2, Epha3, and Ephb3). We also identified 27 transcripts that were specifically upregulated (≥2-fold expression) in BMMSCs relative to HSPCs (Axl, Bmpr1a, Met, Pdgfrb, Fgfr1, Mertk, Cmkor1, Egfr, Epha7, and Ephb4), and 19 transcripts that were specifically upregulated in HSPCs relative to BMMSCs (Ccr1, Csf1r, Csf2ra, Epor, IL6ra, and IL7r). Eleven transcripts were equally expressed (<2-fold upregulation) in HSPCs and BMMSCs (Flt1, Insr, Kdr, Jak1, Agtrl1, Ccr3, Ednrb, Il3ra, Hoxb4, Tnfrsf1a, and Abcb1b), whilst another seven transcripts (Epha6, Epha8, Musk, Ntrk2, Ros1, Srms, and Tnk1) were not expressed in either cell population.

Conclusions

We demonstrate that besides their unique immunophenotype and functional differences, BMC-derived HSPCs and BMMSCs have different molecular receptor signaling transcript profiles linked to cell survival, growth, cell differentiation status, growth factor/cytokine production and genes involved in cell migration/trafficking/adhesion that may be critical to maintain their pluripotency, plasticity, and stem cell function.

The use of transcriptomics to unveil the role of nutrients in Mammalian liver

Liver is the organ primarily responding to diet, and it is crucial in determining plasma carbohydrate, protein, and lipid levels. In addition, it is mainly responsible for transformation of xenobiotics. For these reasons, it has been a target of transcriptomic analyses. In this review, we have covered the works dealing with the response of mammalian liver to different nutritional stimuli such as fasting/feeding, caloric restriction, dietary carbohydrate, cholesterol, fat, protein, bile acid, salt, vitamin, and oligoelement contents. Quality of fats or proteins has been equally addressed, and has the influence of minor dietary components. Other compounds, not purely nutritional as those represented by alcohol and food additives, have been included due to their relevance in processed food. The influence has been studied not only on mRNA but also on miRNA. The wide scope of the technology clearly reflects that any simple intervention has profound changes in many metabolic parameters and that there is a synergy in response when more compounds are included in the intervention. Standardized arrays to systematically test the same genes in all studies and analyzing data to establish patterns of response are required, particularly for RNA sequencing. Moreover, RNA is a valuable, easy-screening ally but always requires further confirmation.

Accurate and precise DNA quantification in the presence of different amplification efficiencies using an improved Cy0 method

Quantitative real-time PCR represents a highly sensitive and powerful technology for the quantification of DNA. Although real-time PCR is well accepted as the gold standard in nucleic acid quantification, there is a largely unexplored area of experimental conditions that limit the application of the Ct method. As an alternative, our research team has recently proposed the Cy0 method, which can compensate for small amplification variations among the samples being compared. However, when there is a marked decrease in amplification efficiency, the Cy0 is impaired, hence determining reaction efficiency is essential to achieve a reliable quantification. The proposed improvement in Cy0 is based on the use of the kinetic parameters calculated in the curve inflection point to compensate for efficiency variations. Three experimental models were used: inhibition of primer extension, non-optimal primer annealing and a very small biological sample. In all these models, the improved Cy0 method increased quantification accuracy up to about 500% without affecting precision. Furthermore, the stability of this procedure was enhanced integrating it with the SOD method. In short, the improved Cy0 method represents a simple yet powerful approach for reliable DNA quantification even in the presence of marked efficiency variations.

Sex-linked transcriptional divergence in the hermaphrodite fungus Neurospora tetrasperma

In the filamentous ascomycete Neurospora tetrasperma, a large (approx. 7 Mbp) region of suppressed recombination surrounds the mating-type (mat) locus. While the remainder of the genome is largely homoallelic, this region of recombinational suppression, extending over 1500 genes, is associated with sequence divergence. Here, we used microarrays to examine how the molecular phenotype of gene expression level is linked to this divergent region, and thus to the mating type. Culturing N. tetrasperma on agar media that induce sexual/female or vegetative/male tissue, we found 196 genes significantly differentially expressed between mat A and mat a mating types. Our data show that the genes exhibiting mat-linked expression are enriched in the region genetically linked to mating type, and sequence and expression divergence are positively correlated. Our results indicate that the phenotype of mat A strains is optimized for traits promoting sexual/female development and the phenotype of mat a strains for vegetative/male development. This discovery of differentially expressed genes associated with mating type provides a link between genotypic and phenotypic divergence in this taxon and illustrates a fungal analogue to sexual dimorphism found among animals and plants.

miRNAs mediate SnRK1-dependent energy signaling in Arabidopsis

The SnRK1 protein kinase, the plant ortholog of mammalian AMPK and yeast Snf1, is activated by the energy depletion caused by adverse environmental conditions. Upon activation, SnRK1 triggers extensive transcriptional changes to restore homeostasis and promote stress tolerance and survival partly through the inhibition of anabolism and the activation of catabolism. Despite the identification of a few bZIP transcription factors as downstream effectors, the mechanisms underlying gene regulation, and in particular gene repression by SnRK1, remain mostly unknown. microRNAs (miRNAs) are 20-24 nt RNAs that regulate gene expression post-transcriptionally by driving the cleavage and/or translation attenuation of complementary mRNA targets. In addition to their role in plant development, mounting evidence implicates miRNAs in the response to environmental stress. Given the involvement of miRNAs in stress responses and the fact that some of the SnRK1-regulated genes are miRNA targets, we postulated that miRNAs drive part of the transcriptional reprogramming triggered by SnRK1. By comparing the transcriptional response to energy deprivation between WT and dcl1-9, a mutant deficient in miRNA biogenesis, we identified 831 starvation genes misregulated in the dcl1-9 mutant, out of which 155 are validated or predicted miRNA targets. Functional clustering analysis revealed that the main cellular processes potentially co-regulated by SnRK1 and miRNAs are translation and organelle function and uncover TCP transcription factors as one of the most highly enriched functional clusters. TCP repression during energy deprivation was impaired in miR319 knockdown (MIM319) plants, demonstrating the involvement of miR319 in the stress-dependent regulation of TCPs. Altogether, our data indicates that miRNAs are components of the SnRK1 signaling cascade contributing to the regulation of specific mRNA targets and possibly tuning down particular cellular processes during the stress response.

Where statistics and molecular microarray experiments biology meet

This review chapter presents a statistical point of view to microarray experiments with the purpose of understanding the apparent contradictions that often appear in relation to their results. We give a brief introduction of molecular biology for nonspecialists. We describe microarray experiments from their construction and the biological principles the experiments rely on, to data acquisition and analysis. The role of epidemiological approaches and sample size considerations are also discussed.

Global effect of inauhzin on human p53-responsive transcriptome

Background

Previously, we reported that Inauhzin (INZ) induces p53 activity and suppresses tumor growth by inhibiting Sirt1. However, it remains unknown whether INZ may globally affect p53-dependent gene expression or not. Herein, we have conducted microarray and real-time PCR analyses of gene expression to determine the global effect of INZ on human p53-responsive transcriptome.

Methodology/Principal Findings

In this study, we conducted microarray analysis followed by PCR validation of general gene expression in HCT116^p53+/+ and HCT116^p53−/− cells treated with or without INZ. Microarray data showed that 324 genes were up-regulated by ≥2.3-fold and 266 genes were down-regulated by ≥2-fold in response to INZ treatment in a p53-dependent manner. GO analysis for these genes further revealed that INZ affects several biological processes, including apoptosis (GO:0006915), cell cycle (GO:0007049), immune system process (GO:0002376), and cell adhesion (GO:0007155), which are in line with p53 functions in cells. Also, pathway and STRING analyses of these genes indicated that the p53-signaling pathway is the most significant pathway responsive to INZ treatment as predicted, since a number of these p53 target genes have been previously reported and some of them were validated by RT-qPCR. Finally, among the 9 tested and highly expressed genes, ACBD4, APOBEC3C, and FLJ14327 could be novel p53 target genes, for they were up-regulated by INZ in HCT116^p53+/+ cells, but not in HCT116^p53−/− cells.

Conclusions/Significance

From our whole genome microarray analysis followed by validation with RT-qPCR, we found that INZ can indeed induce the expression of p53 target genes at a larger scale or globally. Our findings not only verify that INZ indeed activates the p53 signaling pathway, but also provide useful information for identifying novel INZ and/or p53 targets. The global effect of INZ on human p53-responsive transcriptome could also be instrumental to the future design of INZ clinical trials.

A comparison of the accuracy of Illumina HumanHT-12 v3 Expression BeadChip and TaqMan qRT-PCR gene expression results in patient samples from the Tampere Vascular Study

AIMS

This study was set up to compare the accuracy, sensitivity and specificity of gene expression results between the Illumina HumanHT-12 v3 Expression BeadChip (GWE) and the TaqMan qRT-PCR low-density array (LDA) in atherosclerotic plaques.

METHODS

Gene expression levels of 196 genes were determined from 22 atherosclerotic samples and 6 controls with both the GWE and the LDA.

RESULTS

The accuracies of GWE in comparison to that of qRT-PCR for absolute fold changes (FC) 1.2, 1.6, 2.0 and 3.0 between cases and controls were 73.5%, 79.1%, 72.4% and 60.7%, respectively. The correlation of expression measurements between these methods was good (r = 0.87, y = 0.151 + 0.586x), and the Bland-Altman plot showed that for highly up- or down-regulated transcripts, GWE yields lower absolute FC values than LDA.

CONCLUSION

Gene expression results obtained with the GWE are replicated with high accuracy in LDA, even for technically demanding atherosclerotic samples.

DAMTC regulates cytoskeletal reorganization and cell motility in human lung adenocarcinoma cell line: an integrated proteomics and transcriptomics approach

DAMTC (7,8-diacetoxy-4-methylcoumarin) is a thioderivative of 4-methyl coumarin, and previously we have shown that DAMTC is a potent inhibitor of cell growth and an inducer of apoptosis in non-small cell lung cancer (A549) cells. It induces apoptosis through mitochondrial pathway by modulating NF-κB, mitogen-activated protein kinase (MAPK) and p53 pathways. Herein, we explored the genome-wide effects of DAMTC in A549 cells using the concerted approach of transcriptomics and proteomics. In addition to apoptotic pathways, which have been validated earlier, the bioinformatic analysis of microarray data identified small GTPase-mediated signal transduction among the significantly altered biological processes. Interestingly, we observed significant downregulation of some members of the Rho family GTPases in the proteomics data too. Downregulation of Rho GTPases (RhoGDIα (Rho GDP dissociation inhibitor-α, also known as ARHGDIA), Ras homolog family member A, Ras-related C3 botulinum toxin substrate 1 and cell division cycle 42) was validated by western blotting. The Rho protein family is implicated in maintaining the actin filament assembly and cell motility, and we also observed that DAMTC treatment causes actin cytoskeletal reorganization, promotes filopodia formation and inhibits cell motility in A549 cells. The effect of DAMTC treatment on cytoskeleton was reversed after the overexpression of RhoGDIα. In addition, DAMTC augmented the apoptotic effect of etoposide, a proapoptotic chemotherapeutic drug. This elucidation of the mechanism behind DAMTC-induced apoptosis and inhibition of cell motility in A549 cells may make it a potential therapeutic for lung cancer.

Genomic and metabolic responses to methionine-restricted and methionine-restricted, cysteine-supplemented diets in Fischer 344 rat inguinal adipose tissue, liver and quadriceps muscle

BACKGROUND/AIMS

Methionine restriction (MR) is a dietary intervention that increases lifespan, reduces adiposity and improves insulin sensitivity. These effects are reversed by supplementation of the MR diet with cysteine (MRC). Genomic and metabolomic studies were conducted to identify potential mechanisms by which MR induces favorable metabolic effects, and that are reversed by cysteine supplementation.

METHODS

Gene expression was examined by microarray analysis and TaqMan quantitative PCR. Levels of selected proteins were measured by Western blot and metabolic intermediates were analyzed by mass spectrometry.

RESULTS

MR increased lipid metabolism in inguinal adipose tissue and quadriceps muscle while it decreased lipid synthesis in liver. In inguinal adipose tissue, MR not only caused the transcriptional upregulation of genes associated with fatty acid synthesis but also of Lpin1, Pc, Pck1 and Pdk1, genes that are associated with glyceroneogenesis. MR also upregulated lipolysis-associated genes in inguinal fat and led to increased oxidation in this tissue, as suggested by higher levels of methionine sulfoxide and 13-HODE + 9-HODE compared to control-fed (CF) rats. Moreover, MR caused a trend toward the downregulation of inflammation-associated genes in inguinal adipose tissue. MRC reversed most gene and metabolite changes induced by MR in inguinal adipose tissue, but drove the expression of Elovl6, Lpin1, Pc, and Pdk1 below CF levels. In liver, MR decreased levels of a number of long-chain fatty acids, glycerol and glycerol-3-phosphate corresponding with the gene expression data. Although MR increased the expression of genes associated with carbohydrate metabolism, levels of glycolytic intermediates were below CF levels. MR, however, stimulated gluconeogenesis and ketogenesis in liver tissue. As previously reported, sulfur amino acids derived from methionine were decreased in liver by MR, but homocysteine levels were elevated. Increased liver homocysteine levels by MR were associated with decreased cystathionine β-synthase (CBS) protein levels and lowered vitamin B6 and 5-methyltetrahydrofolate (5MeTHF) content. Finally, MR upregulated fibroblast growth factor 21 (FGF21) gene and protein levels in both liver and adipose tissues. MRC reversed some of MR's effects in liver and upregulated the transcription of genes associated with inflammation and carcinogenesis such as Cxcl16, Cdh17, Mmp12, Mybl1, and Cav1 among others. In quadriceps muscle, MR upregulated lipid metabolism-associated genes and increased 3-hydroxybutyrate levels suggesting increased fatty acid oxidation as well as stimulation of gluconeogenesis and glycogenolysis in this tissue.

CONCLUSION

Increased lipid metabolism in inguinal adipose tissue and quadriceps muscle, decreased triglyceride synthesis in liver and the downregulation of inflammation-associated genes are among the factors that could favor the lean phenotype and increased insulin sensitivity observed in MR rats.

Factors regulated by interferon gamma and hypoxia-inducible factor 1A contribute to responses that protect mice from Coccidioides immitis infection

Background

Coccidioidomycosis results from airborne infections caused by either Coccidioides immitis or C. posadasii. Both are pathogenic fungi that live in desert soil in the New World and can infect normal hosts, but most infections are self-limited. Disseminated infections occur in approximately 5% of cases and may prove fatal. Mouse models of the disease have identified strains that are resistant (e.g. DBA/2) or susceptible (e.g. C57BL/6) to these pathogens. However, the genetic and immunological basis for this difference has not been fully characterized.

Results

Microarray technology was used to identify genes that were differentially expressed in lung tissue between resistant DBA/2 and sensitive C57BL/6 mice after infection with C. immitis. Differentially expressed genes were mapped onto biological pathways, gene ontologies, and protein interaction networks, which revealed that innate immune responses mediated by Type II interferon (i.e., IFNG) and the signal transducer and activator of transcription 1 (STAT1) contribute to the resistant phenotype. In addition, upregulation of hypoxia inducible factor 1A (HIF1A), possibly as part of a larger inflammatory response mediated by tumor necrosis factor alpha (TNFA), may also contribute to resistance. Microarray gene expression was confirmed by real-time quantitative PCR for a subset of 12 genes, which revealed that IFNG HIF1A and TNFA, among others, were significantly differentially expressed between the two strains at day 14 post-infection.

Conclusion

These results confirm the finding that DBA/2 mice express more Type II interferon and interferon stimulated genes than genetically susceptible strains and suggest that differential expression of HIF1A may also play a role in protection.

Critical selection of internal control genes for quantitative real-time RT-PCR studies in lipopolysaccharide-stimulated human THP-1 and K562 cells

The choice of internal control genes is important since it may affect the study outcome in RT-qPCR. Indeed, it is well-known that expression levels of traditional internal control genes can vary across tissue types and across experimental settings within one specific tissue type. The aim of this study is an evaluation of a set of housekeeping genes (HKGs) to be used in the normalization of gene expression in vitro different cultured cells, THP-1 and K562. The transcriptional stability of eleven potential internal control genes (RPL37A, ACTB, GAPDH, B(2)M, PPIB, PGK1, PPIA, SDHA, TBP, HPRT1 and RPL13A) were evaluated using RT-qPCR and were compared in different treatment, that was un-stimulated or LPS-stimulated cells. The raw Ct values were determined for each candidate gene at different time points following LPS-stimulated or unstimulated cells. Furthermore, all data were analyzed by the geNorm, BestKeeper, and NormFinder validation programs. Results indicated that PPIB and PGK1 were the most stable internal control genes in this study. RPL13A was found to be the least stable. This study provides the comprehensive reported assessment of internal control genes for use in expression studies in vitro cultured cells. These findings further emphasize the need to accurately validate candidate internal control genes in the study before use in gene expression studies using RT-qPCR.

Sperm mRNA transcripts are indicators of sub-chronic low dose testicular injury in the Fischer 344 rat

Current human reproductive risk assessment methods rely on semen and serum hormone analyses, which are not easily comparable to the histopathological endpoints and mating studies used in animal testing. Because of these limitations, there is a need to develop universal evaluations that reliably reflect male reproductive function. We hypothesized that toxicant-induced testicular injury can be detected in sperm using mRNA transcripts as indicators of insult. To test this, we exposed adult male Fischer 344 rats to low doses of model testicular toxicants and classically characterized the testicular injury while simultaneously evaluating sperm mRNA transcripts from the same animals. Overall, this study aimed to: 1) identify sperm transcripts altered after exposure to the model testicular toxicant, 2,5-hexanedione (HD) using microarrays; 2) expand on the HD-induced transcript changes in a comprehensive time course experiment using qRT-PCR arrays; and 3) test these injury indicators after exposure to another model testicular toxicant, carbendazim (CBZ). Microarray analysis of HD-treated adult Fischer 344 rats identified 128 altered sperm mRNA transcripts when compared to control using linear models of microarray analysis (q<0.05). All transcript alterations disappeared after 3 months of post-exposure recovery. In the time course experiment, time-dependent alterations were observed for 12 candidate transcripts selected from the microarray data based upon fold change and biological relevance, and 8 of these transcripts remained significantly altered after the 3-month recovery period (p<0.05). In the last experiment, 8 candidate transcripts changed after exposure to CBZ (p<0.05). The two testicular toxicants produced distinct molecular signatures with only 4 overlapping transcripts between them, each occurring in opposite directions. Overall, these results suggest that sperm mRNA transcripts are indicators of low dose toxicant-induced testicular injury in the rat.

RanBPM expression regulates transcriptional pathways involved in development and tumorigenesis

RanBPM is a ubiquitous protein that has been reported to regulate several cellular processes through interactions with various proteins. However, it is not known whether RanBPM may regulate gene expression patterns. As it has been shown that RanBPM interacts with a number of transcription factors, we hypothesized that it may have wide ranging effects on gene expression that may explain its function. To test this hypothesis, we generated stable RanBPM shRNA cell lines to analyze the effect of RanBPM on global gene expression. Microarray analyses were conducted comparing the gene expression profile of Hela and HCT116 RanBPM shRNA cells versus control shRNA cells. We identified 167 annotated genes significantly up- or down-regulated in the two cell lines. Analysis of the gene set revealed that down-regulation of RanBPM led to gene expression changes that affect regulation of cell, tissue, and organ development and morphology, as well as biological processes implicated in tumorigenesis. Analysis of Transcription Factor Binding Sites (TFBS) present in the gene set identified several significantly over-represented transcription factors of the Forkhead, HMG, and Homeodomain families of transcription factors, which have previously been demonstrated as having important roles in development and tumorigenesis. In addition, the combined results of these analyses suggested that several signaling pathways were affected by RanBPM down-regulation, including ERK1/2, Wnt, Notch, and PI3K/Akt pathways. Lastly, analysis of selected target genes by quantitative RT-qPCR confirmed the changes revealed by microarray. Several of the genes up-regulated in RanBPM shRNA cells encode proteins with known oncogenic functions, such as the RON tyrosine kinase, the adhesion molecule L1CAM, and transcription factor ELF3/ESE-1, suggesting that RanBPM functions as a tumor suppressor to prevent deregulated expression of these genes. Altogether, these results suggest that RanBPM does indeed function to regulate many genomic events that regulate embryonic, tissue, and cellular development as well as those involved in cancer development and progression.

Reliable and sample saving gene expression analysis approach for diagnostic tool development

This work answers the question of whether it is necessary to hybridize individual instead of pooled RNA samples on microarrays for screening gene targets suitable as diagnostic tools for radiation exposure scenarios, while at the same time meeting comparable microarray quality criteria. For developing new clinical diagnostic tools, a two-stage study design was employed in five projects. At first, pooled and not individual RNA samples were hybridized on microarrays for screening purposes. Potential gene candidates were selected based on their fold-change only. This was followed by a validation/quantification step using individual RNA samples and quantitative RT-PCR. Quality criteria from the screening approach with pooled RNA samples were compared with published data from the MicroArray Quality Control (MAQC) consortium that hybridized each reference RNA sample separately and established quality criteria for microarrays. When comparing both approaches, only insignificant differences for quality criteria such as false positives, sensitivity, specificity, and overall agreement were found. However, material, costs, and time were drastically reduced when hybridizing pooled RNA and gene targets applicable for clinical diagnostic purposes could be successfully selected. In search of new diagnostic tools for radiation exposure scenarios, the two stage study design using either pooled or individual RNA samples on microarrays shows comparable quality criteria, but the RNA pooling approach saves unique material, costs, and efforts and successfully selects gene targets that can be used for the desired diagnostic purposes.

Effects of chronic stress on prefrontal cortex transcriptome in mice displaying different genetic backgrounds

There is increasing evidence that depression derives from the impact of environmental pressure on genetically susceptible individuals. We analyzed the effects of chronic mild stress (CMS) on prefrontal cortex transcriptome of two strains of mice bred for high (HA)and low (LA) swim stress-induced analgesia that differ in basal transcriptomic profiles and depression-like behaviors. We found that CMS affected 96 and 92 genes in HA and LA mice, respectively. Among genes with the same expression pattern in both strains after CMS, we observed robust upregulation of Ttr gene coding transthyretin involved in amyloidosis, seizures, stroke-like episodes, or dementia. Strain-specific HA transcriptome affected by CMS was associated with deregulation of genes involved in insulin secretion (Acvr1c, Nnat, and Pfkm), neuropeptide hormone activity (Nts and Trh), and dopamine receptor mediated signaling pathway (Clic6, Drd1a, and Ppp1r1b). LA transcriptome affected by CMS was associated with genes involved in behavioral response to stimulus (Fcer1g, Rasd2, S100a8, S100a9, Crhr1, Grm5, and Prkcc), immune effector processes (Fcer1g, Mpo, and Igh-VJ558), diacylglycerol binding (Rasgrp1, Dgke, Dgkg, and Prkcc), and long-term depression (Crhr1, Grm5, and Prkcc) and/or coding elements of dendrites (Crmp1, Cntnap4, and Prkcc) and myelin proteins (Gpm6a, Mal, and Mog). The results indicate significant contribution of genetic background to differences in stress response gene expression in the mouse prefrontal cortex.

Sex-specific gene expression during asexual development of Neurospora crassa

The impact of loci that determine sexual identity upon the asexual, dominant stage of fungal life history has been well studied. To investigate their impact, expression differences between strains of different mating type during asexual development were assayed, with RNA sampled from otherwise largely isogenic mat A and mat a strains of Neurospora crassa at early, middle, and late clonal stages of development. We observed significant differences in overall gene expression between mating types across clonal development, especially at late development stages. The expression levels of mating-type genes and pheromone genes were assayed by reverse transcription and quantitative PCR, revealing expression of pheromone and receptor genes in strains of both mating types in all development stages, and revealing that mating type (mat) genes were increasingly expressed over the course of asexual development. Interestingly, among differentially expressed genes, the mat A genotype more frequently exhibited a higher expression level than mat a, and demonstrated greater transcriptional regulatory dynamism. Significant up-regulation of expression was observed for many late light-responsive genes at late asexual development stages. Further investigation of the impact of light and the roles of light response genes in asexual development of both mating types are warranted.

Transcriptional homeostasis of a mangrove species, Ceriops tagal, in saline environments, as revealed by microarray analysis

BACKGROUND

Differential responses to the environmental stresses at the level of transcription play a critical role in adaptation. Mangrove species compose a dominant community in intertidal zones and form dense forests at the sea-land interface, and although the anatomical and physiological features associated with their salt-tolerant lifestyles have been well characterized, little is known about the impact of transcriptional phenotypes on their adaptation to these saline environments.

METHODOLOGY AND PRINCIPAL FINDINGS

We report the time-course transcript profiles in the roots of a true mangrove species, Ceriops tagal, as revealed by a series of microarray experiments. The expression of a total of 432 transcripts changed significantly in the roots of C. tagal under salt shock, of which 83 had a more than 2-fold change and were further assembled into 59 unigenes. Global transcription was stable at the early stage of salt stress and then was gradually dysregulated with the increased duration of the stress. Importantly, a pair-wise comparison of predicted homologous gene pairs revealed that the transcriptional regulations of most of the differentially expressed genes were highly divergent in C. tagal from that in salt-sensitive species, Arabidopsis thaliana.

CONCLUSIONS/SIGNIFICANCE

This work suggests that transcriptional homeostasis and specific transcriptional regulation are major events in the roots of C. tagal when subjected to salt shock, which could contribute to the establishment of adaptation to saline environments and, thus, facilitate the salt-tolerant lifestyle of this mangrove species. Furthermore, the candidate genes underlying the adaptation were identified through comparative analyses. This study provides a foundation for dissecting the genetic basis of the adaptation of mangroves to intertidal environments.

Accurate measurement of the relative abundance of different DNA species in complex DNA mixtures

A molecular tool that can compare the abundances of different DNA sequences is necessary for comparing intergenic or interspecific gene expression. We devised and verified such a tool using a quantitative competitive polymerase chain reaction approach. For this approach, we adapted a competitor array, an artificially made plasmid DNA in which all the competitor templates for the target DNAs are arranged with a defined ratio, and melting analysis for allele quantitation for accurate quantitation of the fractional ratios of competitively amplified DNAs. Assays on two sets of DNA mixtures with explicitly known compositional structures of the test sequences were performed. The resultant average relative errors of 0.059 and 0.021 emphasize the highly accurate nature of this method. Furthermore, the method's capability of obtaining biological data is demonstrated by the fact that it can illustrate the tissue-specific quantitative expression signatures of the three housekeeping genes G6pdx, Ubc, and Rps27 by using the forms of the relative abundances of their transcripts, and the differential preferences of Igf2 enhancers for each of the multiple Igf2 promoters for the transcription.

Microarray expression analysis in delayed cardioprotection: the effect of exercise, AICAR, or metformin and the possible role of AMP-activated protein kinase (AMPK)

AMP-activated protein kinase (AMPK) is an enzyme which may be involved in cardioprotective mechanisms in the ischemic heart. Exercise, AICAR, and metformin, all known activators of AMPK, induce delayed cardioprotection which protects the heart against ischemia-reperfusion injury. The objective was to determine the effect of exercise, AICAR, and metformin on gene expression profile and to demonstrate possible interactions in different genes and functions. Rats were divided into either an exercise, AICAR, metformin, or control group. 3, 12, and 24 h after either a single bout of exercise training, a single injection of AICAR or a single dose of metformin, hearts were removed and gene expression profiles were analyzed in tissue from the left ventricle using Affymetrix gene chip probe arrays. Ingenuity Pathway Analysis (IPA) tool was used to analyze the regulated genes for relevant functions and diseases. Each gene chip identified up to 30,000 different probesets of which Ingenuity identified approximately up to 12,000 genes. A total of 147, 304, and 114 different genes in the left ventricle whose expressions were altered >2.0-fold were identified in the exercise, AICAR, and metformin group, respectively. Seventy eight different genes were overlapping the exercise and AICAR group at 24 h. Ingenuity identified six overlapping genes between the exercise, AICAR, and metformin groups including NR4A3, TNFRSF12A, HBB, PENK, PAP, and MAP4K4. IPA software revealed an overabundance of focus molecules in all three intervention groups involving functions related to cell death, cellular growth and proliferation, gene expression and cancer. Exercise, AICAR, and metformin regulate several genes in the rat myocardium with the majority of overlapping genes observed in the exercise and AICAR group. Changes in gene programming mainly involved inflammatory and opioid systems recognized as cardioprotective pathways. Some of these genes may represent possible candidate genes involved in the molecular mechanisms of AMPK-induced delayed PC.

Endothelial-like cells in chronic thromboembolic pulmonary hypertension: crosstalk with myofibroblast-like cells

Background

Chronic thromboembolic pulmonary hypertension (CTEPH) is characterized by intravascular thrombus formation in the pulmonary arteries.

Recently, it has been shown that a myofibroblast cell phenotype was predominant within endarterectomized tissues from CTEPH patients. Indeed, our recent study demonstrated the existence of not only myofibroblast-like cells (MFLCs), but also endothelial-like cells (ELCs). Under in vitro conditions, a few transitional cells (co-expressing both endothelial- and SM-cell markers) were observed in the ELC population. We hypothesized that MFLCs in the microenvironment created by the unresolved clot may promote the endothelial-mesenchymal transition and/or induce endothelial cell (EC) dysfunction.

Methods

We isolated cells from these tissues and identified them as MFLCs and ELCs. In order to test whether the MFLCs provide the microenvironment which causes EC alterations, ECs were incubated in serum-free medium conditioned by MFLCs, or were grown in co-culture with the MFLCs.

Results

Our experiments demonstrated that MFLCs promoted the commercially available ECs to transit to other mesenchymal phenotypes and/or induced EC dysfunction through inactivation of autophagy, disruption of the mitochondrial reticulum, alteration of the SOD-2 localization, and decreased ROS production. Indeed, ELCs included a few transitional cells, lost the ability to form autophagosomes, and had defective mitochondrial structure/function. Moreover, rapamycin reversed the phenotypic alterations and the gene expression changes in ECs co-cultured with MFLCs, thus suggesting that this agent had beneficial therapeutic effects on ECs in CTEPH tissues.

Conclusions

It is possible that the microenvironment created by the stabilized clot stimulates MFLCs to induce EC alterations.

Multiple platform assessment of the EGF dependent transcriptome by microarray and deep tag sequencing analysis

Background

Epidermal Growth Factor (EGF) is a key regulatory growth factor activating many processes relevant to normal development and disease, affecting cell proliferation and survival. Here we use a combined approach to study the EGF dependent transcriptome of HeLa cells by using multiple long oligonucleotide based microarray platforms (from Agilent, Operon, and Illumina) in combination with digital gene expression profiling (DGE) with the Illumina Genome Analyzer.

Results

By applying a procedure for cross-platform data meta-analysis based on RankProd and GlobalAncova tests, we establish a well validated gene set with transcript levels altered after EGF treatment. We use this robust gene list to build higher order networks of gene interaction by interconnecting associated networks, supporting and extending the important role of the EGF signaling pathway in cancer. In addition, we find an entirely new set of genes previously unrelated to the currently accepted EGF associated cellular functions.

Conclusions

We propose that the use of global genomic cross-validation derived from high content technologies (microarrays or deep sequencing) can be used to generate more reliable datasets. This approach should help to improve the confidence of downstream in silico functional inference analyses based on high content data.

Overexpression of Pax6 in mouse cornea directly alters corneal epithelial cells: changes in immune function, vascularization, and differentiation

PURPOSE

To assess whether Pax6 functions directly in the cornea, a corneal-preferred promoter was used to overexpress Pax6 specifically in the cornea.

METHODS

Transgenic mice harboring a construct containing mouse Pax6 coding sequences fused downstream of the aldehyde dehydrogenase 3a1 (Aldh3a1) promoter were generated (Pax6 Tg). Pax6 expression was analyzed by Western blot and immunohistochemistry. Eye sections were stained with hematoxylin and eosin, Schiff reagent, and fluorescein, to assess morphologic changes, the presence of goblet cells, and barrier integrity, respectively. Gene expression changes in mildly affected Pax6 Tg corneas were compared to age-matched, wild-type (WT) corneas by microarray analysis and quantitative PCR. Promoter regulation of several differentially expressed genes was examined by monitoring luciferase activity of reporter constructs after cotransfection with Pax6 in COS7 cells.

RESULTS

Corneal overexpression of Pax6 produces an abnormal cornea with altered epithelial cell morphology, neovascularization, immune cell invasion, and a compromised barrier; the lens appeared normal. Major changes in expression of genes involved in immune function, vascularization, and epithelial differentiation occurred in corneas from Pax6 Tg versus WT mice. The keratin (K) profile was dramatically altered in the Pax6 Tg corneas, as were several components of the Wnt signaling pathway. In severely affected Pax6 Tg corneas, K12 was reduced, and Pax6 was redistributed into the cytoplasm. Promoters from the chitinase 3-like 3, Wnt inhibitory factor 1, and fms-related tyrosine kinase 1/soluble VEGF receptor genes were upregulated five-, seven-, and threefold, respectively, by Pax6 in transfected COS7 cells.

CONCLUSIONS

Pax6 functions directly to maintain normal, corneal epithelial cells.

Global analysis of gene expression in mineralizing fish vertebra-derived cell lines: new insights into anti-mineralogenic effect of vanadate

BACKGROUND

Fish has been deemed suitable to study the complex mechanisms of vertebrate skeletogenesis and gilthead seabream (Sparus aurata), a marine teleost with acellular bone, has been successfully used in recent years to study the function and regulation of bone and cartilage related genes during development and in adult animals. Tools recently developed for gilthead seabream, e.g. mineralogenic cell lines and a 4 × 44K Agilent oligo-array, were used to identify molecular determinants of in vitro mineralization and genes involved in anti-mineralogenic action of vanadate.

RESULTS

Global analysis of gene expression identified 4,223 and 4,147 genes differentially expressed (fold change - FC > 1.5) during in vitro mineralization of VSa13 (pre-chondrocyte) and VSa16 (pre-osteoblast) cells, respectively. Comparative analysis indicated that nearly 45% of these genes are common to both cell lines and gene ontology (GO) classification is also similar for both cell types. Up-regulated genes (FC > 10) were mainly associated with transport, matrix/membrane, metabolism and signaling, while down-regulated genes were mainly associated with metabolism, calcium binding, transport and signaling. Analysis of gene expression in proliferative and mineralizing cells exposed to vanadate revealed 1,779 and 1,136 differentially expressed genes, respectively. Of these genes, 67 exhibited reverse patterns of expression upon vanadate treatment during proliferation or mineralization.

CONCLUSIONS

Comparative analysis of expression data from fish and data available in the literature for mammalian cell systems (bone-derived cells undergoing differentiation) indicate that the same type of genes, and in some cases the same orthologs, are involved in mechanisms of in vitro mineralization, suggesting their conservation throughout vertebrate evolution and across cell types. Array technology also allowed identification of genes differentially expressed upon exposure of fish cell lines to vanadate and likely involved in its anti-mineralogenic activity. Many were found to be unknown or they were never associated to bone homeostasis previously, thus providing a set of potential candidates whose study will likely bring insights into the complex mechanisms of tissue mineralization and bone formation.

Genome-wide expression profiling of human lymphoblastoid cell lines identifies CHL1 as a putative SSRI antidepressant response biomarker

AIMS

Selective serotonin reuptake inhibitors (SSRIs) are the most commonly used class of antidepressants for treating major depression. However, approximately 30% of patients do not respond sufficiently to first-line antidepressant drug treatment and require alternative therapeutics. Genome-wide studies searching for SSRI response DNA biomarkers or studies of candidate serotonin-related genes so far have given inconclusive or contradictory results. Here, we present an alternative transcriptome-based genome-wide approach for searching antidepressant drug-response biomarkers by using drug-effect phenotypes in human lymphoblastoid cell lines (LCLs).

MATERIALS & METHODS

We screened 80 LCLs from healthy adult female individuals for growth inhibition by paroxetine. A total of 14 LCLs with reproducible high and low sensitivities to paroxetine (seven from each phenotypic group) were chosen for genome-wide expression profiling with commercial microarrays.

RESULTS

The most notable genome-wide transcriptome difference between LCLs displaying high versus low paroxetine sensitivities was a 6.3-fold lower (p = 0.0000256) basal expression of CHL1, a gene coding for a neuronal cell adhesion protein implicated in correct thalamocortical circuitry, schizophrenia and autism. The microarray findings were confirmed by real-time PCR (36-fold lower CHL1 expression levels in the high paroxetine sensitivity group). Several additional genes implicated in synaptogenesis or in psychiatric disorders, including ARRB1, CCL5, DDX60, DDX60L, ENDOD1, ENPP2, FLT1, GABRA4, GAP43, MCTP2 and SPRY2, also differed by more than 1.5-fold and a p-value of less than 0.005 between the two paroxetine sensitivity groups, as confirmed by real-time PCR experiments.

CONCLUSION

Genome-wide transcriptional profiling of in vitro phenotyped LCLs identified CHL1 and additional genes implicated in synaptogenesis and brain circuitry as putative SSRI response biomarkers. This method might be used as a preliminary tool for searching for potential depression treatment biomarkers.

Micro-RNA expression in cisplatin resistant germ cell tumor cell lines

Background

We compared microRNA expression patterns in three cisplatin resistant sublines derived from paternal cisplatin sensitive germ cell tumor cell lines in order to improve our understanding of the mechanisms of cisplatin resistance.

Methods

Three cisplatin resistant sublines (NTERA-2-R, NCCIT-R, 2102EP-R) showing 2.7-11.3-fold increase in drug resistance after intermittent exposure to increasing doses of cisplatin were compared to their parental counterparts, three well established relatively cisplatin sensitive germ cell tumor cell lines (NTERA-2, NCCIT, 2102EP). Cells were cultured and total RNA was isolated from all 6 cell lines in three independent experiments. RNA was converted into cDNA and quantitative RT-PCR was run using 384 well low density arrays covering almost all (738) known microRNA species of human origin.

Results

Altogether 72 of 738 (9.8%) microRNAs appeared differentially expressed between sensitive and resistant cell line pairs (NTERA-2R/NTERA-2 = 43, NCCIT-R/NCCIT = 53, 2102EP-R/2102EP = 15) of which 46.7-95.3% were up-regulated (NTERA-2R/NTERA-2 = 95.3%, NCCIT-R/NCCIT = 62.3%, 2102EP-R/2102EP = 46.7%). The number of genes showing differential expression in more than one of the cell line pairs was 34 between NTERA-2R/NTERA-2 (79%) and NCCIT-R/NCCIT (64%), and 3 and 4, respectively, between these two cell lines and 2102EP-R/2102EP (about 27%). Only the has-miR-10b involved in breast cancer invasion and metastasis and has-miR-512-3p appeared to be up-regulated (2-3-fold) in all three cell lines. The hsa-miR-371-373 cluster (counteracting cellular senescence and linked with differentiation potency), as well as hsa-miR-520c/-520h (inhibiting the tumor suppressor p21) were 3.9-16.3 fold up-regulated in two of the three cisplatin resistant cell lines. Several new micro-RNA species missing an annotation towards cisplatin resistance could be identified. These were hsa-miR-512-3p/-515/-517/-518/-525 (up to 8.1-fold up-regulated) and hsa-miR-99a/-100/-145 (up to 10-fold down-regulated).

Conclusion

Examining almost all known human micro-RNA species confirmed the miR-371-373 cluster as a promising target for explaining cisplatin resistance, potentially by counteracting wild-type P53 induced senescence or linking it with the potency to differentiate. Moreover, we describe for the first time an association of the up-regulation of micro-RNA species such as hsa-miR-512-3p/-515/-517/-518/-525 and down-regulation of hsa-miR-99a/-100/-145 with a cisplatin resistant phenotype in human germ cell tumors. Further functional analyses are warranted to gain insight into their role in drug resistance.

Development of commercial gene-expression-based signatures: review of the scientific strategies

Many scientific articles have been published that use gene-expression-based technologies to discriminate a trait of interest, typically a disease subgroup, within a patient population. However, few gene-expression-based signatures have at present reached the market and become a financially and clinically successful product. The technological, scientific and medical challenges, the regulatory environment and the financial considerations are all essential parts of the development process. Here we discuss the scientific aspects of successfully developing a gene-expression-based signature and review the global strategy of six products that made it to the market. We also present a point-to-point guide that should help researchers to successfully develop genomic signatures, thus paving the way towards personalized medicine.

Transcriptome profiling of embryonic development rate in rainbow trout advanced backcross introgression lines

In rainbow trout (Oncorhynchus mykiss) and other fishes, embryonic development rate is an ecologically and evolutionarily important trait that is closely associated with survival and physiological performance later in life. To identify genes differentially regulated in fast and slow-developing embryos of rainbow trout, we examined gene expression across developmental time points in rainbow trout embryos possessing alleles linked to a major quantitative trait loci (QTL) for fast versus slow embryonic development rate. Whole genome expression microarray analyses were conducted using embryos from a fourth generation backcross family, whereby each backcross generation involved the introgression of the fast-developing alleles for a major development rate QTL into a slow-developing clonal line of rainbow trout. Embryos were collected at 15, 19, and 28 days post-fertilization; sex and QTL genotype were determined using molecular markers, and cDNA from 48 embryos were used for microarray analysis. A total of 183 features were identified with significant differences between embryonic development rate genotypes. Genes associated with cell cycle growth, muscle contraction and protein synthesis were expressed significantly higher in embryos with the fast-developing allele (Clearwater) than those with the slow-developing allele (Oregon State University), which may associate with fast growth and early body mass construction in embryo development. Across time points, individuals with the fast-developing QTL allele appeared to have earlier onset of these developmental processes when compared to individuals with the slow development alleles, even as early as 15 days post-fertilization. Differentially expressed candidate genes chosen for linkage mapping were localized primarily to regions outside of the major embryonic development rate QTL, with the exception of a single gene (very low-density lipoprotein receptor precursor).