LATENT RANK CHANGE DETECTION FOR ANALYSIS OF SPLICE-JUNCTION MICROARRAYS WITH NONLINEAR EFFECTS

Alternative splicing of gene transcripts greatly expands the functional capacity of the genome, and certain splice isoforms may indicate specific disease states such as cancer. Splice junction microarrays interrogate thousands of splice junctions, but data analysis is difficult and error prone because of the increased complexity compared to differential gene expression analysis. We present Rank Change Detection (RCD) as a method to identify differential splicing events based upon a straightforward probabilistic model comparing the over-or underrepresentation of two or more competing isoforms. RCD has advantages over commonly used methods because it is robust to false positive errors due to nonlinear trends in microarray measurements. Further, RCD does not depend on prior knowledge of splice isoforms, yet it takes advantage of the inherent structure of mutually exclusive junctions, and it is conceptually generalizable to other types of splicing arrays or RNA-Seq. RCD specifically identifies the biologically important cases when a splice junction becomes more or less prevalent compared to other mutually exclusive junctions. The example data is from different cell lines of glioblastoma tumors assayed with Agilent microarrays.

miRNA analysis of prostate cancer by quantitative real time PCR: Comparison between formalin-fixed paraffin embedded and fresh-frozen tissue

Objective: Micro RNA (miRNA) is a class of small noncoding RNA that plays a major role in the regulation of gene expression, which has been related to cancer behavior. The possibility of analyzing miRNA from the archives of pathology laboratories is exciting, as it allows for large retrospective studies. Formalin is the most common fixative used in the surgical pathology routine, and its promotion of nucleic acid degradation is well known. Our aim is to compare miRNA profiles from formalin-fixed paraffin embedded (FFPE) tissues with fresh-frozen prostate cancer tissues. Methods: The expression of 14 miRNAs was determined by quantitative real time polymerase chain reaction (qRT-PCR) in 5 paired fresh-frozen and FFPE tissues, which were representative of prostate carcinoma. Results: There was a very good correlation of the miRNA expression of miR-let7c and miR-32 between the fresh-frozen and FFPE tissues, with Pearson`s correlation coefficients of 0.927 (P = 0.023) and 0.960 (P = 0.010), respectively. For the remaining miRNAs, the correlation was good with Spearman correlation coefficient of 0.638 (P < 0.001). Conclusion: Analysis of miRNAs from routinely processed and stored FFPE prostate tissue is feasible for some miRNAs using qRT-PCR. Further studies should be conducted to confirm the reliability of using stock tissues for miRNA expression determination. (C) 2011 Elsevier Inc. All rights reserved.

Change in expression of miR-let7c, miR-100, and miR-218 from high grade localized prostate cancer to metastasis

Objective: MicroRNAs (miRNAs) are small noncoding regulatory RNAs (19-25 nucleotides) that play a major role in regulation of gene expression. They are responsible for the control of fundamental cellular processes that has been reported to be involved in human tumorigenesis. The characterization of miRNA profiles in human tumors is crucial for the understanding of carcinogenesis processes, finding of new tumor markers, and discovering of specific targets for the development of innovative therapies. The aim of this study is to find miRNAs involved in prostate cancer progression comparing the profile of miRNA expressed by localized high grade carcinoma and bone metastasis. Material and methods: Two groups of tumors where submitted to analyses. The first is characterized by 18 patients who underwent radical prostatectomy for treatment of localized high grade prostate carcinoma (PC) with mean Gleason score 8.6, all staged pT3. The second group is composed of 4 patients with metastatic, androgen-independent prostate carcinoma, and 2 PC cell lines. LNCaP derived from a metastatic PC to a lymph node, and another derived from an obstructive, androgen-independent PC (PcBRA1). Expression analysis of 14 miRNAs was carried out using quantitative RT-PCR. Results: miR-let7c, miR-100, and miR-218 were significantly overexpressed by all localized high GS, pT3 PC in comparison with metastatic carcinoma. (35.065 vs. 0.996 P < 0.001), (55.550 vs. 8.314, P = 0.010), and (33.549 vs. 2.748, P = 0.001), respectively. Conclusion: We hypothesize that miR-let7c, miR-100, and miR-218 may be involved in the process of metastasization of PC, and their role as controllers of the expression of RAS, c-myc, Laminin 5 beta 3, THAP2, SMARCA5, and BAZ2A should be matter of additional studies. (C) 2011 Elsevier Inc. All rights reserved.

Farnesoic acid O-methyl transferase (FAMeT) isoforms: Conserved traits and gene expression patterns related to caste differentiation in the stingless bee, Melipona scutellaris

Farnesoic acid O-methyl transferase (FAMeT) is the enzyme that catalyzes the formation of methyl farnesoate (MF) from farnesoic acid (FA) in the biosynthetic pathway of juvenile hormone (JH). This work reports the cloning, sequencing, and expression of FAMeT gene from the stingless bee Melipona scutellaris (MsFAMeT). The MsFAMeT in silica analysis showed that greatest sequence similarity is found in Apis mellifera and other insects, while relatively less similarity is shown in crustaceans. Evidence of alternative splicing of a 27 nucleotide (nt) microexon explains the presence of the detected isoforms, 1 and 2. The expression analysis of the two isoforms showed a marked difference when castes were compared, suggesting that they could be involved differently in the JH metabolism in M. scutellaris, providing new insights for the comprehension of female plasticity.

Mutation analysis of the CDKN2A promoter in Australian melanoma families

Approximately 50% of all melanoma families worldwide show linkage to 9p21-22, but only about half of these have been shown to contain germ line CDKN2A mutations. It has been hypothesized that a proportion of these families carry mutations in the noncoding regions of CDKN2A. Several Canadian families have been reported to carry a mutation in the 5' UTR, at position -34 relative to the start site, which gives rise to a novel AUG translation initiation codon that markedly decreases translation from the wild-type AUG (Liu et al., 1999). Haplotype sharing in these Canadian families suggested that this mutation is of British origin. We sequenced 1,327 base pairs (bp) of CDKN2A, making up 1,116 bp of the 5' UTR and promoter, all of exon 1, and 61 bp of intron 1, in at least one melanoma case from 110 Australian families with three or more affected members known not to carry mutations within the p16 coding region. In addition, 431 bp upstream of the start codon was sequenced in an additional 253 affected probands from two-case melanoma families for which the CDKN2A mutation status was unknown. Several known polymorphisms at positions -33, -191, -493, and -735 were detected, in addition to four novel variants at positions 120, -252, -347, and -981 relative to the start codon. One of the probands from a two-case family was found to have the previously reported Q50R mutation. No family member was found to carry the mutation at position -34 or any other disease-associated mutation. For further investigation of noncoding CDKN2A mutations that may affect transcription, allele-specific expression analysis was carried out in 31 of the families with at least three affected members who showed either complete or indeterminate 9p haplotype sharing without CDKN2A exonic mutations. Reverse transcription polymerase chain reaction and automated sequencing showed expression of both CDKN2A alleles in all family members tested. The lack of CDKN2A promoter mutations and the absence of transcriptional silencing in the germ line of this cohort of families suggest that mutations in the promoter and 5' UTR play a very limited role in melanoma predisposition. (C) 2001 Wiley-Liss, Inc.

A novel gene expression signature in peripheral blood mononuclear cells for early detection of colorectal cancer.

BACKGROUND: Early detection and treatment of colorectal adenomatous polyps (AP) and colorectal cancer (CRC) is associated with decreased mortality for CRC. However, accurate, non-invasive and compliant tests to screen for AP and early stages of CRC are not yet available. A blood-based screening test is highly attractive due to limited invasiveness and high acceptance rate among patients. AIM: To demonstrate whether gene expression signatures in the peripheral blood mononuclear cells (PBMC) were able to detect the presence of AP and early stages CRC. METHODS: A total of 85 PBMC samples derived from colonoscopy-verified subjects without lesion (controls) (n = 41), with AP (n = 21) or with CRC (n = 23) were used as training sets. A 42-gene panel for CRC and AP discrimination, including genes identified by Digital Gene Expression-tag profiling of PBMC, and genes previously characterised and reported in the literature, was validated on the training set by qPCR. Logistic regression analysis followed by bootstrap validation determined CRC- and AP-specific classifiers, which discriminate patients with CRC and AP from controls. RESULTS: The CRC and AP classifiers were able to detect CRC with a sensitivity of 78% and AP with a sensitivity of 46% respectively. Both classifiers had a specificity of 92% with very low false-positive detection when applied on subjects with inflammatory bowel disease (n = 23) or tumours other than CRC (n = 14). CONCLUSION: This pilot study demonstrates the potential of developing a minimally invasive, accurate test to screen patients at average risk for colorectal cancer, based on gene expression analysis of peripheral blood mononuclear cells obtained from a simple blood sample.

Impact of normalization on miRNA microarray expression profiling.

Profiling miRNA levels in cells with miRNA microarrays is becoming a widely used technique. Although normalization methods for mRNA gene expression arrays are well established, miRNA array normalization has so far not been investigated in detail. In this study we investigate the impact of normalization on data generated with the Agilent miRNA array platform. We have developed a method to select nonchanging miRNAs (invariants) and use them to compute linear regression normalization coefficients or variance stabilizing normalization (VSN) parameters. We compared the invariants normalization to normalization by scaling, quantile, and VSN with default parameters as well as to no normalization using samples with strong differential expression of miRNAs (heart-brain comparison) and samples where only a few miRNAs are affected (by p53 overexpression in squamous carcinoma cells versus control). All normalization methods performed better than no normalization. Normalization procedures based on the set of invariants and quantile were the most robust over all experimental conditions tested. Our method of invariant selection and normalization is not limited to Agilent miRNA arrays and can be applied to other data sets including those from one color miRNA microarray platforms, focused gene expression arrays, and gene expression analysis using quantitative PCR.

Analysis of angiotensin II- and ACTH-driven mineralocorticoid functions and omental adiposity in a non-genetic, hyperadipose female rat phenotype

The hypothalamic damage induced by neonatal treatment with monosodium l-glutamate (MSG) induces several metabolic abnormalities, resulting in a rat hyperleptinemic-hyperadipose phenotype. This study was conducted to explore the impact of the neonatal MSG treatment, in the adult (120 days old) female rat on: (a) the in vivo and in vitro mineralocorticoid responses to ACTH and angiotensin II (AII); (b) the effect of leptin on ACTH- and AII-stimulated mineralocorticoid secretions by isolated corticoadrenal cells; and (c) abdominal adiposity characteristics. Our data indicate that, compared with age-matched controls, MSG rats displayed: (1) enhanced and reduced mineralocorticoid responses to ACTH and AII treatments, respectively, effects observed in both in vivo and in vitro conditions; (2) adrenal refractoriness to the inhibitory effect of exogenous leptin on ACTH-stimulated aldosterone output by isolated adrenocortical cells; and (3) distorted omental adiposity morphology and function. This study supports that the adult hyperleptinemic MSG female rat is characterized by enhanced ACTH-driven mineralocorticoid function, impaired adrenal leptin sensitivity, and disrupted abdominal adiposity function. MSG rats could counteract undesirable effects of glucocorticoid excess, by developing a reduced AII-driven mineralocorticoid function. Thus, chronic hyperleptinemia could play a protective role against ACTH-mediated allostatic loads in the adrenal leptin resistant, MSG female rat phenotype.

Small nucleolar RNA expression profiling identifies potential prognostic markers in peripheral T-cell lymphoma.

Peripheral T-cell lymphoma (PTCL) is a rare, heterogeneous type of non-Hodgkin lymphoma (NHL) that, in general, is associated with a poor clinical outcome. Therefore, a current major challenge is the discovery of new prognostic tools for this disease. In the present study, a cohort of 122 patients with PTCL was collected from a multicentric T-cell lymphoma consortium (TENOMIC). We analyzed the expression of 80 small nucleolar RNAs (snoRNAs) using high-throughput quantitative PCR. We demonstrate that snoRNA expression analysis may be useful in both the diagnosis of some subtypes of PTCL and the prognostication of both PTCL-not otherwise specified (PTCL-NOS; n = 26) and angio-immunoblastic T-cell lymphoma (AITL; n = 46) patients treated with chemotherapy. Like miRNAs, snoRNAs are globally down-regulated in tumor cells compared with their normal counterparts. In the present study, the snoRNA signature was robust enough to differentiate anaplastic large cell lymphoma (n = 32) from other PTCLs. For PTCL-NOS and AITL, we obtained 2 distinct prognostic signatures with a reduced set of 3 genes. Of particular interest was the prognostic value of HBII-239 snoRNA, which was significantly over-expressed in cases of AITL and PTCL-NOS that had favorable outcomes. Our results suggest that snoRNA expression profiles may have a diagnostic and prognostic significance for PTCL, offering new tools for patient care and follow-up.

Correlation between vasoconstrictor roles and mRNA expression of alpha1-adrenoceptor subtypes in blood vessels of genetically engineered mice.

We examined the contribution of each alpha(1)-adrenoceptor (AR) subtype in noradrenaline (NAd)-evoked contraction in the thoracic aortas and mesenteric arteries of mice. Compared with the concentration-response curves (CRCs) for NAd in the thoracic aortas of wild-type (WT) mice, the CRCs of mutant mice showed a significantly lower sensitivity. The pD(2) value in rank order is as follows: WT mice (8.21) > alpha(1B)-adrenoceptor knockout (alpha(1B)-KO) (7.77) > alpha(1D)-AR knockout (alpha(1D)-KO) (6.44) > alpha(1B)- and alpha(1D)-AR double knockout (alpha(1BD)-KO) (5.15). In the mesenteric artery, CRCs for NAd did not differ significantly between either WT (6.52) and alpha(1B)-KO mice (7.12) or alpha(1D)-KO (6.19) and alpha(1BD)-KO (6.29) mice. However, the CRC maximum responses to NAd in alpha(1D)- and alpha(1BD)-KO mice were significantly lower than those in WT and alpha(1B)-KO mice. Except in the thoracic aortas of alpha(1BD)-KO mice, the competitive antagonist prazosin inhibited the contraction response to NAd with high affinity. However, prazosin produced shallow Schild slopes in the vessels of mice lacking the alpha(1D)-AR gene. In the thoracic aorta, pA(2) values in WT mice for KMD-3213 and BMY7378 were 8.25 and 8.46, respectively, and in alpha(1B)-KO mice they were 8.49 and 9.13, respectively. In the mesenteric artery, pA(2) values in WT mice for KMD-3213 and BMY7378 were 8.34 and 7.47, respectively, and in alpha(1B)-KO mice they were 8.11 and 7.82, respectively. These pharmacological findings were in fairly good agreement with findings from comparison of CRCs, with the exception of the mesenteric arteries of WT and alpha(1B)-KO mice, which showed low affinities to BMY7378. We performed a quantitative analysis of the mRNA expression of each alpha(1)-AR subtype in these vessels in order to examine the correlation between mRNA expression level and the predominance of each alpha(1)-AR subtype in mediating vascular contraction. The rank order of each alpha(1)-AR subtype in terms of its vasoconstrictor role was in fairly good agreement with the level of expression of mRNA of each subtype, that is, alpha(1D)-AR > alpha(1B)-AR > alpha(1A)-AR in the thoracic aorta and alpha(1D)-AR > alpha(1A)-AR > alpha(1B)-AR in the mesenteric artery. No dramatic compensatory change of alpha(1)-AR subtype in mutant mice was observed in pharmacological or quantitative mRNA expression analysis.

Semiquantitative RT-PCR measurement of gene expression in rat tissues including a correction for varying cell size and number

Background: Current methodology of gene expression analysis limits the possibilities of comparison between cells/tissues of organs in which cell size and/or number changes as a consequence of the study (e.g. starvation). A method relating the abundance of specific mRNA copies per cell may allow direct comparison or different organs and/or changing physiological conditions. Methods: With a number of selected genes, we analysed the relationship of the number of bases and the fluorescence recorded at a present level using cDNA standards. A lineal relationship was found between the final number of bases and the length of the transcript. The constants of this equation and those of the relationship between fluorescence and number of bases in cDNA were determined and a general equation linking the length of the transcript and the initial number of copies of mRNA was deduced for a given pre-established fluorescence setting. This allowed the calculation of the concentration of the corresponding mRNAs per g of tissue. The inclusion of tissue RNA and the DNA content per cell, allowed the calculation of the mRNA copies per cell. Results: The application of this procedure to six genes: Arbp, cyclophilin, ChREBP, T4 deiodinase 2, acetyl-CoA carboxylase 1 and IRS-1, in liver and retroperitoneal adipose tissue of food-restricted rats allowed precise measures of their changes irrespective of the shrinking of the tissue, the loss of cells or changes in cell size, factors that deeply complicate the comparison between changing tissue conditions. The percentage results obtained with the present methods were essentially the same obtained with the delta-delta procedure and with individual cDNA standard curve quantitative RT-PCR estimation. Conclusion: The method presented allows the comparison (i.e. as copies of mRNA per cell) between different genes and tissues, establishing the degree of abundance of the different molecular species tested.

Expression of O⁶-methylguanine-DNA methyltransferase in childhood medulloblastoma.

Medulloblastomas (MB) are the most common malignant brain tumors in childhood. Alkylator-based drugs are effective agents in the treatment of patients with MB. In several tumors, including malignant glioma, elevated O(6)-methylguanine-DNA methyltransferase (MGMT) expression levels or lack of MGMT promoter methylation have been found to be associated with resistance to alkylating chemotherapeutic agents such as temozolomide (TMZ). In this study, we examined the MGMT status of MB and central nervous system primitive neuroectodermal tumor (PNET) cells and two large sets of primary MB. In seven MB/PNET cell lines investigated, MGMT promoter methylation was detected only in D425 human MB cells as assayed by the qualitative methylation-specific PCR and the more quantitative pyrosequencing assay. In D425 human MB cells, MGMT mRNA and protein expression was clearly lower when compared with the MGMT expression in the other MB/PNET cell lines. In MB/PNET cells, sensitivity towards TMZ and 1-(2-chloroethyl)-3-cyclohexyl-1-nitrosourea (CCNU) correlated with MGMT methylation and MGMT mRNA expression. Pyrosequencing in 67 primary MB samples revealed a mean percentage of MGMT methylation of 3.7-92% (mean: 13.25%, median: 10.67%). Percentage of MGMT methylation and MGMT mRNA expression as determined by quantitative RT-PCR correlated inversely (n = 46; Pearson correlation r (2) = 0.14, P = 0.01). We then analyzed MGMT mRNA expression in a second set of 47 formalin-fixed paraffin-embedded primary MB samples from clinically well-documented patients treated within the prospective randomized multicenter trial HIT'91. No association was found between MGMT mRNA expression and progression-free or overall survival. Therefore, it is not currently recommended to use MGMT mRNA expression analysis to determine who should receive alkylating agents and who should not.

Integrative analysis of the Regulatory Region of the FGFR3 Oncogene

The study of transcriptional regulation often needs the integration of diverse yet independent data. In the present work, sequence conservation, predic-tion of transcription factor binding sites (TFBS) and gene expression analysis have been applied to the detection of putative transcription factor (TF) modules in the regulatory region of the FGFR3 oncogene. Several TFs with conserved binding sites in the FGFR3 regulatory region have shown high positive or negative corre-lation with FGFR3 expression both in urothelial carcinoma and in benign nevi. By means of conserved TF cluster analysis, two different TF modules have been iden-tified in the promoter and first intron of FGFR3 gene. These modules contain acti-vating AP2, E2F, E47 and SP1 binding sites plus motifs for EGR with possible repressor function.

The expression of estrogen receptors as well as GREB1, c-MYC, and cyclin D1, estrogen-regulated genes implicated in proliferation, is increased in peritoneal endometriosis.

OBJECTIVE: To analyze the expression of estrogen receptors α and β as well as their target genes implicated in proliferation, c-myc, cyclin D1, and GREB1, in the endometrium of women with or without endometriosis. DESIGN: Expression analysis in human tissue. SETTING: University hospitals and a clinic. PATIENT(S): Ninety-one premenopausal women (59 patients with endometriosis and 32 controls) undergoing laparoscopic surgery. INTERVENTION(S): Biopsies were obtained at time of surgery, performed during the proliferative phase of the cycle. MAIN OUTCOME MEASURE(S): Estrogen receptors α and β as well as c-myc, cyclin D1, and GREB1 mRNA expression levels were determined by quantitative reverse transcriptase-polymerase chain reaction. Tissue localization of these estrogen-regulated genes was analyzed by immunohistochemistry. RESULT(S): Estrogen receptors α and β as well as c-myc, cyclin D1, and GREB1 mRNA expression levels were increased in ectopic tissue in comparison with both normal and eutopic endometrium. Estrogen receptor mRNA levels also were upregulated in the eutopic peritoneal tissue of patients with endometriosis. Cyclin D1 and GREB1 expression was augmented in eutopic endometrium. c-myc, cyclin D1, and GREB1 proteins exhibited a nuclear localization in ectopic endometrial tissue. CONCLUSION(S): This constitutes the first report of increased expression of GREB1, as well as cyclin D1 and c-myc, in peritoneal endometriotic lesions, implicating these proteins in estrogen-dependent growth in this context.

Metabolic gene expression changes in astrocytes in Multiple Sclerosis cerebral cortex are indicative of immune-mediated signaling.

Emerging as an important correlate of neurological dysfunction in Multiple Sclerosis (MS), extended focal and diffuse gray matter abnormalities have been found and linked to clinical manifestations such as seizures, fatigue and cognitive dysfunction. To investigate possible underlying mechanisms we analyzed the molecular alterations in histopathological normal appearing cortical gray matter (NAGM) in MS. By performing a differential gene expression analysis of NAGM of control and MS cases we identified reduced transcription of astrocyte specific genes involved in the astrocyte-neuron lactate shuttle (ANLS) and the glutamate-glutamine cycle (GGC). Additional quantitative immunohistochemical analysis demonstrating a CX43 loss in MS NAGM confirmed a crucial involvement of astrocytes and emphasizes their importance in MS pathogenesis. Concurrently, a Toll-like/IL-1β signaling expression signature was detected in MS NAGM, indicating that immune-related signaling might be responsible for the downregulation of ANLS and GGC gene expression in MS NAGM. Indeed, challenging astrocytes with immune stimuli such as IL-1β and LPS reduced their ANLS and GGC gene expression in vitro. The detected upregulation of IL1B in MS NAGM suggests inflammasome priming. For this reason, astrocyte cultures were treated with ATP and ATP/LPS as for inflammasome activation. This treatment led to a reduction of ANLS and GGC gene expression in a comparable manner. To investigate potential sources for ANLS and GGC downregulation in MS NAGM, we first performed an adjuvant-driven stimulation of the peripheral immune system in C57Bl/6 mice in vivo. This led to similar gene expression changes in spinal cord demonstrating that peripheral immune signals might be one source for astrocytic gene expression changes in the brain. IL1B upregulation in MS NAGM itself points to a possible endogenous signaling process leading to ANLS and GGC downregulation. This is supported by our findings that, among others, MS NAGM astrocytes express inflammasome components and that astrocytes are capable to release Il-1β in-vitro. Altogether, our data suggests that immune signaling of immune- and/or central nervous system origin drives alterations in astrocytic ANLS and GGC gene regulation in the MS NAGM. Such a mechanism might underlie cortical brain dysfunctions frequently encountered in MS patients.