MicroRNAs miR-16 and miR-519 control meningioma cell proliferation via overlapping transcriptomic programs shared with the RNA-binding protein HuR

Introduction

Meningiomas are the most common type of primary central nervous system tumors. In about 80% cases, these tumors are benign and grow very slowly, but the remainder 20% can unlock higher proliferation rates and become malignant. In this study we examined two miRs, miR-16 and miR-519, and evaluated their role in tumorigenesis and cell growth in human meningioma.

Methods

A cohort of 60 intracranial grade 1 and grade 2 human meningioma plus 20 healthy meningeal tissues was used to quantify miR-16 and miR-519 expressions. Cell growth and dose-response assays were performed in two human meningioma cell lines, Ben-Men-1 (benign) and IOMM-Lee (aggressive). Transcriptomes of IOMM-lee cells were measured after both miR-mimics transfection, followed by integrative bioinformatics to expand on available data.

Results

In tumoral tissues, we detected decreased levels of miR-16 and miR-519 when compared with arachnoid cells of healthy patients (miR-16: P=8.7e-04; miR-519: P=3.5e-07). When individually overexpressing these miRs in Ben-Men-1 and IOMM-Lee, we observed that each showed reduced growth (P<0.001). In IOMM-Lee cell transcriptomes, downregulated genes, among which ELAVL1/HuR (miR-16: P=6.1e-06; miR-519:P=9.38e-03), were linked to biological processes such as mitotic cell cycle regulation, pre-replicative complex, and brain development (FDR<1e-05). Additionally, we uncovered a specific transcriptomic signature of miR-16/miR-519-dysregulated genes which was highly enriched in HuR targets (>6-fold; 79.6% of target genes).

Discussion

These results were confirmed on several public transcriptomic and microRNA datasets of human meningiomas, hinting that the putative tumor suppressor effect of these miRs is mediated, at least in part, via HuR direct or indirect inhibition.

Integrative genomics analysis of nasal intestinal-type adenocarcinomas demonstrates the major role of CACNA1C and paves the way for a simple diagnostic tool in male woodworkers

Background

Nasal intestinal-type adenocarcinomas (ITAC) are strongly related to chronic wood dust exposure: The intestinal phenotype relies on CDX2 overexpression but underlying molecular mechanisms remain unknown. Our objectives were to investigate transcriptomic and methylation differences between healthy non-exposed and tumor olfactory cleft mucosae and to compare transcriptomic profiles between non-exposed, wood dust-exposed and ITAC mucosa cells.

Methods

We conducted a prospective monocentric study (NCT0281823) including 16 woodworkers with ITAC, 16 healthy exposed woodworkers and 13 healthy, non-exposed, controls. We compared tumor samples with healthy non-exposed samples, both in transcriptome and in methylome analyses. We also investigated wood dust-induced transcriptome modifications of exposed (without tumor) male woodworkers’ samples and of contralateral sides of woodworkers with tumors. We conducted in parallel transcriptome and methylome analysis, and then, the transcriptome analysis was focused on the genes highlighted in methylome analysis. We replicated our results on dataset GSE17433.

Results

Several clusters of genes enabled the distinction between healthy and ITAC samples. Transcriptomic and IHC analysis confirmed a constant overexpression of CDX2 in ITAC samples, without any specific DNA methylation profile regarding the CDX2 locus. ITAC woodworkers also exhibited a specific transcriptomic profile in their contralateral (non-tumor) olfactory cleft, different from that of other exposed woodworkers, suggesting that they had a different exposure or a different susceptibility. Two top-loci (CACNA1C/CACNA1C-AS1 and SLC26A10) were identified with a hemimethylated profile, but only CACNA1C appeared to be overexpressed both in transcriptomic analysis and in immunohistochemistry.

Conclusions

Several clusters of genes enable the distinction between healthy mucosa and ITAC samples even in contralateral nasal fossa thus paving the way for a simple diagnostic tool for ITAC in male woodworkers. CACNA1C might be considered as a master gene of ITAC and should be further investigated.

Trial registration: NIH ClinicalTrials, NCT0281823, registered May 23^d 2016, https://www.clinicaltrials.gov/NCT0281823.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-021-01122-5.

Characteristics of Lymphoma in Patients with Inflammatory Bowel Disease: A Systematic Review

BACKGROUND

Lymphoma is a dreaded complication of inflammatory bowel diseases [IBD]. Knowledge about lymphoma in patients with IBD is limited to epidemiological data and the description of risk factors. We performed a systematic review to describe the clinical characteristics and prognosis of lymphoma in patients with IBD.

METHODS

Electronic databases were searched up to June 1, 2020. All published clinical characteristics of lymphoma occurring in patients with IBD were collected.

RESULTS

Eleven studies were included. A total of 589 lymphomas were described in patients with IBD. As seen in de novo lymphoma, non-Hodgkin's lymphoma [NHL] was the most common histological subtype [83.9%]. Diffuse large B-cell lymphoma [DLBCL] and follicular lymphoma were the most well-represented NHL in patients with IBD [30% and 13% respectively]. Two main differences were observed in comparison with de novo lymphoma: primary intestinal lymphoma [PIL] represented a large proportion of lymphoma in patients with IBD [22-75%] whereas mucosa-associated lymphoid tissue [MALT] lymphoma was under-represented. Epstein-Barr virus [EBV]-positive status was observed in a large proportion of tumours [44-75%]. Survival data of lymphoma in patients with IBD were similar to those of de novo lymphoma.

DISCUSSION

This systematic review first highlights that PIL [especially DLBCL subtype] is significantly more frequent in patients with IBD and represents the most common entity. Conversely, MALT lymphoma is extremely rare in the IBD population. However, the overall quality of the evidence is low. Further studies are required to better define lymphoma characteristics in patients with IBD.

Genomic and molecular alterations in human inflammatory bowel disease-associated colorectal cancer

Patients with inflammatory bowel disease are at increased risk of colorectal cancer, which has worse prognosis than sporadic colorectal cancer. Until recently, understanding of pathogenesis in inflammatory bowel disease-associated colorectal cancer was restricted to the demonstration of chromosomic/microsatellite instabilities and aneuploidy. The advance of high-throughput sequencing technologies has highlighted the complexity of the pathobiology and revealed recurrently mutated genes involved in the RTK/RAS, PI3K, WNT, and TGFβ pathways, leading to potentially new targetable mutations. Moreover, alterations of mitochondrial DNA and the dysregulation of non-coding sequences have also been described, as well as several epigenetic modifications. Although recent studies have brought new insights into pathobiology and raised the prospect of innovative therapeutic approaches, the understanding of colorectal carcinogenesis in inflammatory bowel disease and how it differs from sporadic colorectal cancer remains not fully elucidated. Further studies are required to better understand the pathogenesis and molecular alterations leading to human inflammatory bowel disease-associated colorectal cancer.

Ki-67 and MCM6 labeling indices are correlated with overall survival in anaplastic oligodendroglioma, IDH1-mutant and 1p/19q-codeleted: a multicenter study from the French POLA network

Anaplastic oligodendroglioma (AO), IDH-mutant and 1p/19q codeleted (IDHmut+/1p19qcodel), is a high-grade glioma with only limited prognostic markers. The primary objective of this study was to evaluate, by immunohistochemistry, the prognostic value of two proliferation markers, MCM6 and Ki-67, in a large series of IDHmut+/1p19qcodel AO included in the POLA ("Prise en charge des Oligodendrogliomes Anaplasiques") French national multicenter network. We additionally examined the transcriptome obtained from this series to understand the functional pathways dysregulated with the mRNA overexpression of these two markers. The labeling indices (LI) of MCM6 and Ki-67 were obtained via computer-assisted color image analyses on immunostained AO tissues of the cohort (n = 220). Furthermore, a subgroup of AO (n = 68/220) was used to perform transcriptomic analyses. A high LI of either MCM6 (≥50%) or Ki-67 (≥15%) correlated with shorter overall survival, both in univariate (P = 0.013 and P = 0.004, respectively) and multivariate analyses (P = 0.027; multivariate Cox model including age, mitotic index, MCM6 and Ki-67). MCM6 and Ki-67 LI also correlated with overall survival in an additional retrospective cohort of 30 grade II IDHmut+/1p19qcodel oligodendrogliomas. The prognostic value of MCM6 mRNA level was confirmed in The Cancer Genome Atlas (TCGA) IDHmut+/1p19qcodel gliomas. The transcriptomic approach revealed that high transcriptional expressions of MCM6 and MKI67 were both linked positively with cell cycle progression, DNA replication, mitosis, pro-neural phenotype as well as neurogenesis, and negatively with microglial cell activation, immune response, positive regulation of myelination, oligodendrocyte development, beta-amyloid binding and postsynaptic specialization. In conclusion, the overexpression of MCM6 and/or Ki-67 is independently associated to shorter overall survival in IDHmut+/1p19qcodel AO. These two easy-to-use and cost-effective markers could thus be used concurrently in routine pathology practice. Additionally, the transcriptomic analyses showed that AO with high proliferation index have down-regulated immune response and lower microglial cells activation, and bears pro-neural phenotype.

Inherited disorders of cobalamin metabolism disrupt nucleocytoplasmic transport of mRNA through impaired methylation/phosphorylation of ELAVL1/HuR

The molecular mechanisms that underlie the neurological manifestations of patients with inherited diseases of vitamin B12 (cobalamin) metabolism remain to date obscure. We observed transcriptomic changes of genes involved in RNA metabolism and endoplasmic reticulum stress in a neuronal cell model with impaired cobalamin metabolism. These changes were related to the subcellular mislocalization of several RNA binding proteins, including the ELAVL1/HuR protein implicated in neuronal stress, in this cell model and in patient fibroblasts with inborn errors of cobalamin metabolism and Cd320 knockout mice. The decreased interaction of ELAVL1/HuR with the CRM1/exportin protein of the nuclear pore complex and its subsequent mislocalization resulted from hypomethylation at R-217 produced by decreased S-adenosylmethionine and protein methyl transferase CARM1 and dephosphorylation at S221 by increased protein phosphatase PP2A. The mislocalization of ELAVL1/HuR triggered the decreased expression of SIRT1 deacetylase and genes involved in brain development, neuroplasticity, myelin formation, and brain aging. The mislocalization was reversible upon treatment with siPpp2ca, cobalamin, S-adenosylmethionine, or PP2A inhibitor okadaic acid. In conclusion, our data highlight the key role of the disruption of ELAVL1/HuR nuclear export, with genomic changes consistent with the effects of inborn errors of Cbl metabolisms on brain development, neuroplasticity and myelin formation.

Transcriptional orchestration of mitochondrial homeostasis in a cellular model of PGC-1-related coactivator-dependent thyroid tumor

The PGC-1 (Peroxisome proliferator-activated receptor Gamma Coactivator-1) family of coactivators (PGC-1α, PGC-1β, and PRC) plays a central role in the transcriptional control of mitochondrial biogenesis and oxidative phosphorylation (OXPHOS) processes. These coactivators integrate mitochondrial energy production into cell metabolism using complementary pathways. The XTC.UC1 cell line is a mitochondria-rich model of thyroid tumors whose biogenesis is almost exclusively dependent on PRC. Here we aim to propose an integrative view of the cellular pathways regulated by PRC through integration of cDNA and miRNA microarray data and chromatin immunoprecipitation results obtained from XTC.UC1 cells invalidated for PRC. This study showes that PRC induces a complex network of cellular functions interacting with at least one to five of the studied transcription factors (Estrogen Related Receptor alpha, ERR1; Nuclear-Respiratory Factors, NRF1 and NRF2; cAMP Response Element Binding, CREB; and Ying Yang, YY1). Our data confirm that ERR1 is a key partner of PRC in the regulation of mitochondrial functions and suggest a potential role of this complex in RNA processing. PRC is also involved in transcriptional regulatory complexes targeting 12 miRNAs, five of which are involved in the control of the OXPHOS process. Our findings demonstrate that the PRC coactivator can act in complex with several transcription factors and regulate miRNA expression to control the fine regulation of main metabolic functions in the cell. Therefore, in PGC-1α/β-associated pathologies, PRC, as a metabolic sensor, may ensure mitochondrial homeostasis.

Identification of genomic differences among peripheral arterial beds in atherosclerotic and healthy arteries

Calcification is independently associated with cardiovascular events and morbidity. The calcification burden in atherosclerotic lesions quantitatively and qualitatively differs between arterial beds. Cardiovascular risk factors (CVRF) differentially affect plaque development between arterial beds. The aim of this study was to evaluate the impact of CVRF on atherosclerotic plaque calcification and to further study the molecular arterial heterogeneity that could account for these differences. Histological analysis was performed on atherosclerotic plaques from 153 carotid, 97 femoral and 28 infrapopliteal arteries. CVRF showed minor associations with plaque calcification: age and hypertension affected only the overall presence of calcification but not the type of the calcification, which significantly differed between arterial beds. Transcriptome analysis revealed distinct gene expression profiles associated with each territory in atherosclerotic and healthy arteries. Canonical pathway analysis showed the preferential involvement of immune system-related processes in both atherosclerotic and healthy carotid arteries. Bone development-related genes were among those mostly enriched in atherosclerotic and healthy femoral arteries, which are more prone to developing endochondral calcification. This study highlights the heterogeneous nature of arteries from different peripheral vascular beds and contributes to a better understanding of atherosclerosis formation and evolution.

Relapsed diffuse large B-cell lymphoma present different genomic profiles between early and late relapses

Despite major advances in first-line treatment, a significant proportion of patients with diffuse large B-cell lymphoma (DLBCL) will experience treatment failure. Prognosis is particularly poor for relapses occurring less than one year after the end of first-line treatment (early relapses/ER) compared to those occurring more than one year after (late relapses/LR). To better understand genomic alterations underlying the delay of relapse, we identified copy number variations (CNVs) on 39 tumor samples from a homogeneous series of patients included in the Collaborative Trial in Relapsed Aggressive Lymphoma (CORAL) prospective study. To identify CNVs associated with ER or LR, we devised an original method based on Significance Analysis of Microarrays, a permutation-based method which allows control of false positives due to multiple testing. Deletions of CDKN2A/B (28%) and IBTK (23%) were frequent events in relapsed DLBCLs. We identified 56 protein-coding genes and 25 long non-coding RNAs with significantly differential CNVs distribution between ER and LR DLBCLs, with a false discovery rate < 0.05. In ER DLBCLs, CNVs were related to transcription regulation, cell cycle and apoptosis, with duplications of histone H1T (31%), deletions of DIABLO (26%), PTMS (21%) and CK2B (15%). In LR DLBCLs, CNVs were related to immune response, with deletions of B2M (20%) and CD58 (10%), cell proliferation regulation, with duplications of HES1 (25%) and DVL3 (20%), and transcription regulation, with MTERF4 deletions (20%). This study provides new insights into the genetic aberrations in relapsed DLBCLs and suggest pathway-targeted therapies in ER and LR DLBCLs.

Implication of molecular vascular smooth muscle cell heterogeneity among arterial beds in arterial calcification

Vascular calcification is a strong and independent predictive factor for cardiovascular complications and mortality. Our previous work identified important discrepancies in plaque composition and calcification types between carotid and femoral arteries. The objective of this study is to further characterize and understand the heterogeneity in vascular calcification among vascular beds, and to identify molecular mechanisms underlying this process. We established ECLAGEN biocollection that encompasses human atherosclerotic lesions and healthy arteries from different locations (abdominal, thoracic aorta, carotid, femoral, and infrapopliteal arteries) for histological, cell isolation, and transcriptomic analysis. Our results show that lesion composition differs between these locations. Femoral arteries are the most calcified arteries overall. They develop denser calcifications (sheet-like, nodule), and are highly susceptible to osteoid metaplasia. These discrepancies may derive from intrinsic differences between SMCs originating from these locations, as microarray analysis showed specific transcriptomic profiles between primary SMCs isolated from each arterial bed. These molecular differences translated into functional disparities. SMC from femoral arteries showed the highest propensity to mineralize due to an increase in basal TGFβ signaling. Our results suggest that biological heterogeneity of resident vascular cells between arterial beds, reflected by our transcriptomic analysis, is critical in understanding plaque biology and calcification, and may have strong implications in vascular therapeutic approaches.

VEGF121, is predictor for survival in activated B-cell-like diffuse large B-cell lymphoma and is related to an immune response gene signature conserved in cancers

Tumor microenvironment including endothelial and immune cells plays a crucial role in tumor progression and has been shown to dramatically influence cancer survival. In this study, we investigated the clinical relevance of the gene expression of key mediators of angiogenesis, VEGF isoforms 121, 165, and 189, and their receptors (VEGFR-1 and R-2) in a cohort of patients (n = 37) with relapsed/refractory diffuse large B-cell lymphoma (DLBCL) from the Collaborative Trial in Relapsed Aggressive Lymphoma (CORAL). In patients with ABC-like DLBCL, but not in patients with GCB-like DLBCL, low VEGF₁₂₁ expression was associated with a significantly better survival than in those with high VEGF₁₂₁ level: 4-year overall survival at 100% vs 36% (p = .011), respectively. A specific gene signature including 57 genes was correlated to VEGF₁₂₁ expression level and was analyzed using a discovery process in 1,842 GSE datasets of public microarray studies. This gene signature was significantly expressed in other cancer datasets and was associated with immune response. In conclusion, low VEGF₁₂₁ expression level was significantly associated with a good prognosis in relapsed/refractory ABC-like DLBCL, and with a well-conserved gene-expression profiling signature related to immune response. These findings pave the way for rationalization of drugs targeting immune response in refractory/relapsed ABC-like DLBCL.

Cytoplasmic overexpression of RNA-binding protein HuR is a marker of poor prognosis in meningioma, and HuR knockdown decreases meningioma cell growth and resistance to hypoxia

HuR regulates cytoplasmic mRNA stability and translatability, and the HuR expression level has been shown to correlate with poor disease outcome in several cancer types; however, the prognostic value and potential pro-oncogenic properties of HuR in meningioma remain unclear. Thus, in the present study, we analysed 85 meningioma tissue samples to establish the relationship between HuR expression, tumour cell proliferation, and/or patient survival. In addition, we examined the anti-proliferative effects of HuR knockdown in two meningioma cell lines (IOMM-Lee and Ben-Men-1) and conducted transcriptome-wide analyses (IOMM-Lee cells) to elucidate the molecular consequences of HuR knockdown. The results of the present study showed HuR cytoplasmic expression to correlate positively with tumour grade (p = 1.2 × 10^-8 ) and negatively with progression-free and overall survival (p = 0.01) time in human meningioma tissues. In vitro, siHuR-induced HuR knockdown was shown to reduce the growth of both Ben-Men-1 (p = 2 × 10^-8 ) and IOMM-Lee (p = 4 × 10^-9 ) cells. Transcriptome analyses revealed HuR knockdown in IOMM-Lee cells to deregulate the HIF1A signalling pathway (p = 1.5 × 10^-6 ) and to up-regulate the expression of genes essential for the assembly of the cytoplasmic mRNA processing body, global genome nucleotide-excision repair, poly(A)-specific ribonuclease activity, the positive regulation of apoptosis and of cell cycle arrest, and the negative regulation of RNA splicing [p(FDR) < 0.001]. Interestingly, HuR knockdown under hypoxic culture conditions further potentiated the effects of HuR knockdown on cell growth, apoptosis, and HIF1A expression. We thus conclude that cytoplasmic HuR expression is a marker of poor prognosis in meningioma and that HuR is a promising potential therapeutic target for use in tumours refractory to standard therapies. Copyright © 2017 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

A composite score associated with spontaneous operational tolerance in kidney transplant recipients

New challenges in renal transplantation include using biological information to devise a useful clinical test for discerning high- and low-risk patients for individual therapy and ascertaining the best combination and appropriate dosages of drugs. Based on a 20-gene signature from a microarray meta-analysis performed on 46 operationally tolerant patients and 266 renal transplant recipients with stable function, we applied the sparse Bolasso methodology to identify a minimal and robust combination of six genes and two demographic parameters associated with operational tolerance. This composite score of operational tolerance discriminated operationally tolerant patients with an area under the curve of 0.97 (95% confidence interval 0.94-1.00). The score was not influenced by immunosuppressive treatment, center of origin, donor type, or post-transplant lymphoproliferative disorder history of the patients. This composite score of operational tolerance was significantly associated with both de novo anti-HLA antibodies and tolerance loss. It was validated by quantitative polymerase chain reaction using independent samples and demonstrated specificity toward a model of tolerance induction. Thus, our score would allow clinicians to improve follow-up of patients, paving the way for individual therapy.

Neuromuscular junction immaturity and muscle atrophy are hallmarks of the ColQ-deficient mouse, a model of congenital myasthenic syndrome with acetylcholinesterase deficiency

The collagen ColQ anchors acetylcholinesterase (AChE) in the synaptic cleft of the neuromuscular junction (NMJ). It also binds MuSK and perlecan/dystroglycan, 2 signaling platforms of the postsynaptic domain. Mutations in ColQ cause a congenital myasthenic syndrome (CMS) with AChE deficiency. Because the absence of AChE does not fully explain the complexity of the syndrome and there is no curative treatment for the disease, we explored additional potential targets of ColQ by conducting a large genetic screening of ColQ-deficient mice, a model for CMS with AChE deficiency, and analyzed their NMJ and muscle phenotypes. We demonstrated that ColQ controls the development and the maturation of the postsynaptic domain by regulating synaptic gene expression. Notably, ColQ deficiency leads to an up-regulation of the 5 subunits of the nicotinic acetylcholine receptor (AChR), leading to mixed mature and immature AChRs at the NMJ of adult mice. ColQ also regulates the expression of extracellular matrix (ECM) components. However, whereas the ECM mRNAs were down-regulated in vitro, compensation seemed to occur in vivo to maintain normal levels of these mRNAs. Finally, ColQ deficiency leads to a general atrophic phenotype and hypoplasia that affect fast muscles. This study points to new specific hallmarks for this CMS.-Sigoillot, S. M., Bourgeois, F., Karmouch, J., Molgó, J., Dobbertin, A., Chevalier, C., Houlgatte, R., Léger, J., Legay, C. Neuromuscular junction immaturity and muscle atrophy are hallmarks of the ColQ-deficient mouse, a model of congenital myasthenic syndrome with acetylcholinesterase deficiency.

Identification of master genes involved in liver key functions through transcriptomics and epigenomics of methyl donor deficiency in rat: relevance to nonalcoholic liver disease

SCOPE

Our study aims to investigate molecular events associated to methyl donor deficiency (MDD) by analyzing the transcriptome and the methylome of MDD rats in liver.

METHODS AND RESULTS

Twenty-one-day-old rats born to mothers fed either with a standard diet or a MDD diet during gestation and lactation were compared. From a total of 44 000 probes for 26 456 genes, we found two gene clusters in MDD rats whose expression levels had significant differences compared with controls: 3269 overexpressed (p < 0.0009) and 2841 underexpressed (p < 0.0004) genes. Modifications of DNA methylation were found in the promoter regions of 1032 genes out of 14 981 genes. Ontological analyses revealed that these genes are mainly involved in glucose and lipid metabolism, nervous system, coagulation, ER stress, and mitochondrial function.

CONCLUSION

Putative master genes exhibiting changes in both gene expression and DNA methylation are limited to 266 genes and are mainly involved in the renin-angiotensin system (n = 3), mitochondrion metabolism (n = 18), and phospholipid homeostasis (n = 3). Most of these master genes participate in nonalcoholic fatty liver disease. The adverse effects of MDD on the metabolic process indicate the beneficial impact of folate and vitamin B12, especially during the perinatal period.

Profiling of androgen response in rainbow trout pubertal testis: relevance to male gonad development and spermatogenesis

The capacity of testicular somatic cells to promote and sustain germ cell differentiation is largely regulated by sexual steroids and notably androgens. In fish species the importance of androgens is emphasized by their ability to induce sex reversal of the developing fries and to trigger spermatogenesis. Here we studied the influence of androgens on testicular gene expression in trout testis using microarrays. Following treatment of immature males with physiological doses of testosterone or 11-ketotestosterone, 418 genes that exhibit changes in expression were identified. Interestingly, the activity of testosterone appeared stronger than that of 11-ketotestosterone. Expression profiles of responsive genes throughout testis development and in isolated germ cells confirmed androgens to mainly affect gene expression in somatic cells. Furthermore, specific clusters of genes that exhibit regulation coincidently with changes in the natural circulating levels of androgens during the reproductive cycle were highlighted, reinforcing the physiological significance of these data. Among somatic genes, a phylogenetic footprinting study identified putative androgen response elements within the proximal promoter regions of 42 potential direct androgen target genes. Finally, androgens were also found to alter the germ line towards meiotic expression profiles, supporting the hypothesis of a role for the somatic responsive genes in driving germ cell fate. This study significantly increases our understanding of molecular pathways regulated by androgens in vertebrates. The highly cyclic testicular development in trout together with functions associated with regulated genes reveal potential mechanisms for androgen actions in tubule formation, steroid production, germ cell development and sperm secretion.

Immune response and mitochondrial metabolism are commonly deregulated in DMD and aging skeletal muscle

Duchenne Muscular Dystrophy (DMD) is a complex process involving multiple pathways downstream of the primary genetic insult leading to fatal muscle degeneration. Aging muscle is a multifactorial neuromuscular process characterized by impaired muscle regeneration leading to progressive atrophy. We hypothesized that these chronic atrophying situations may share specific myogenic adaptative responses at transcriptional level according to tissue remodeling. Muscle biopsies from four young DMD and four AGED subjects were referred to a group of seven muscle biopsies from young subjects without any neuromuscular disorder and explored through a dedicated expression microarray. We identified 528 differentially expressed genes (out of 2,745 analyzed), of which 328 could be validated by an exhaustive meta-analysis of public microarray datasets referring to DMD and Aging in skeletal muscle. Among the 328 validated co-expressed genes, 50% had the same expression profile in both groups and corresponded to immune/fibrosis responses and mitochondrial metabolism. Generalizing these observed meta-signatures with large compendia of public datasets reinforced our results as they could be also identified in other pathological processes and in diverse physiological conditions. Focusing on the common gene signatures in these two atrophying conditions, we observed enrichment in motifs for candidate transcription factors that may coordinate either the immune/fibrosis responses (ETS1, IRF1, NF1) or the mitochondrial metabolism (ESRRA). Deregulation in their expression could be responsible, at least in part, for the same transcriptome changes initiating the chronic muscle atrophy. This study suggests that distinct pathophysiological processes may share common gene responses and pathways related to specific transcription factors.

Nitric oxide and calcium participate in the fine regulation of mitochondrial biogenesis in follicular thyroid carcinoma cells

Members of the peroxisome proliferator-activated receptor γ coactivator-1 family (i.e. PGC-1α, PGC-1β, and the PGC-1-related coactivator (PRC)) are key regulators of mitochondrial biogenesis and function. These regulators serve as mediators between environmental or endogenous signals and the transcriptional machinery governing mitochondrial biogenesis. The FTC-133 and RO82 W-1 follicular thyroid carcinoma cell lines, which present significantly different numbers of mitochondria, metabolic mechanisms, and expression levels of PRC and PGC-1α, may employ retrograde signaling in response to respiratory dysfunction. Nitric oxide (NO) and calcium have been hypothesized to participate in this activity. We investigated the effects of the S-nitroso-N-acetyl-DL-penicillamine-NO donor, on the expression of genes involved in mitochondrial biogenesis and cellular metabolic functions in FTC-133 and RO82 W-1 cells by measuring lactate dehydrogenase and cytochrome c oxidase (COX) activities. We studied the action of ionomycin and 1,2-bis(2-aminophenoxy)ethane-N,N,N',N'-tetraacetic acid-acetoxymethyl ester (BAPTA/AM) (i.e. a calcium ionophore and a cytosolic calcium chelator) on whole genome expression and mitochondrial biogenesis in RO82 W-1 cells. COX activity and the dynamics of endoplasmic reticulum and mitochondrial networks were analyzed in regard to calcium-modulating treatments. In the FTC-133 and RO82 W-1 cells, the mitochondrial biogenesis induced by NO was mainly related to PRC expression as a retrograde mitochondrial signaling. Ionomycin diminished COX activity and negatively regulated PRC-mediated mitochondrial biogenesis in RO82 W-1 cells, whereas BAPTA/AM produced the opposite effects with a reorganization of the mitochondrial network. This is the first demonstration that NO and calcium regulate mitochondrial biogenesis through the PRC pathway in thyroid cell lines.

Meta-analysis of muscle transcriptome data using the MADMuscle database reveals biologically relevant gene patterns

Background

DNA microarray technology has had a great impact on muscle research and microarray gene expression data has been widely used to identify gene signatures characteristic of the studied conditions. With the rapid accumulation of muscle microarray data, it is of great interest to understand how to compare and combine data across multiple studies. Meta-analysis of transcriptome data is a valuable method to achieve it. It enables to highlight conserved gene signatures between multiple independent studies. However, using it is made difficult by the diversity of the available data: different microarray platforms, different gene nomenclature, different species studied, etc.

Description

We have developed a system tool dedicated to muscle transcriptome data. This system comprises a collection of microarray data as well as a query tool. This latter allows the user to extract similar clusters of co-expressed genes from the database, using an input gene list. Common and relevant gene signatures can thus be searched more easily. The dedicated database consists in a large compendium of public data (more than 500 data sets) related to muscle (skeletal and heart). These studies included seven different animal species from invertebrates (Drosophila melanogaster, Caenorhabditis elegans) and vertebrates (Homo sapiens, Mus musculus, Rattus norvegicus, Canis familiaris, Gallus gallus). After a renormalization step, clusters of co-expressed genes were identified in each dataset. The lists of co-expressed genes were annotated using a unified re-annotation procedure. These gene lists were compared to find significant overlaps between studies.

Conclusions

Applied to this large compendium of data sets, meta-analyses demonstrated that conserved patterns between species could be identified. Focusing on a specific pathology (Duchenne Muscular Dystrophy) we validated results across independent studies and revealed robust biomarkers and new pathways of interest. The meta-analyses performed with MADMuscle show the usefulness of this approach. Our method can be applied to all public transcriptome data.

Identification of proteomic signatures of mantle cell lymphoma, small lymphocytic lymphoma, and marginal zone lymphoma biopsies by surface enhanced laser desorption/ionization-time of flight mass spectrometry

Mantle cell lymphoma (MCL), small lymphocytic lymphoma (SLL), and marginal zone lymphoma (MZL) are small B-cell non-Hodgkin lymphomas (NHLs) that may be difficult to distinguish. In order to identify specific proteomic biomarkers, differential proteomic analysis of these three NHLs was performed using surface enhanced laser desorption/ionization-time of flight mass spectrometry (SELDI-TOF-MS). Whole cell lysates obtained from 18 MCL, 20 SLL, and 20 MZL biopsies were applied on two different ProteinChips (cationic and anionic). Hierarchical clustering and discriminating scores combined with an innovative bio-informatics microdissection strategy allowed us to distinguish specific lymphoma proteomic signatures based on the expression of 37 protein peaks. SELDI-assisted protein purification combined with nano-liquid chromatography (LC) quadrupole-time of flight tandem mass spectrometry (Q-TOF MS/MS) was used to identify proteins overexpressed in both MCL and SLL tumors. Among them two histones, H2B and H4, were identified in MCL tumor biopsies and the signal recognition particle 9 kDa protein, SRP9, in SLL tumor biopsies.

MADGene: retrieval and processing of gene identifier lists for the analysis of heterogeneous microarray datasets

Summary: MADGene is a software environment comprising a web-based database and a java application. This platform aims at unifying gene identifiers (ids) and performing gene set analysis. MADGene allows the user to perform inter-conversion of clone and gene ids over a large range of nomenclatures relative to 17 species. We propose a set of 23 functions to facilitate the analysis of gene sets and we give two microarray applications to show how MADGene can be used to conduct meta-analyses.

Availability: The MADGene resources are freely available online from http://www.madtools.org, a website dedicated to the analysis and annotation of DNA microarray data.

Contact:gerard.ramstein@univ-nantes.fr

Supplementary information:Supplementary data are available at Bioinformatics online.

A parallel scheme for comparing transcription factor binding sites matrices

Gene regulation implies many mechanisms. Their identification is a crucial task to construct regulatory networks, and is necessary to understand the pathology in many cases. This requires the identification of transcription factors that play a role in regulation. Numerous motif discovery tools are now available. Combining efficiently their results appears useful for comparing and clustering these motifs in order to reduce redundancies and to identify the corresponding transcription factor. We develop a method that produces, compares and clusters a set of motifs and identifies some close motifs in databases like JASPAR and the public version of Transfac. Unlike previous comparison methods, where each matrix column is compared independently, we have developed a global method to compare motifs that also helps to reduce the number of false positives. We also propose an original graph motif model that generalizes the classical position specific pattern matrices. Finally, we present an application of our method to study ChIP-chip data sets in the context of an eukaryotic organism.

Nutritional programming in the rat is linked to long-lasting changes in nutrient sensing and energy homeostasis in the hypothalamus

Background

Nutrient deficiency during perinatal development is associated with an increased risk to develop obesity, diabetes and hypertension in the adulthood. However, the molecular mechanisms underlying the developmental programming of the metabolic syndrome remain largely unknown.

Methodology/Principal Findings

Given the essential role of the hypothalamus in the integration of nutritional, endocrine and neuronal cues, here we have analyzed the profile of the hypothalamus transcriptome in 180 days-old rats born to dams fed either a control (200 g/kg) or a low-protein (80 g/kg) diet through pregnancy and lactation. From a total of 26 209 examined genes, 688 were up-regulated and 309 down-regulated (P<0.003) by early protein restriction. Further bioinformatic analysis of the data revealed that perinatal protein restriction permanently alters the expression of two gene clusters regulating common cellular processes. The first one includes several gate keeper genes regulating insulin signaling and nutrient sensing. The second cluster encompasses a functional network of nuclear receptors and co-regulators of transcription involved in the detection and use of lipid nutrients as fuel which, in addition, link temporal and nutritional cues to metabolism through their tight interaction with the circadian clock.

Conclusions/Significance

Collectively, these results indicate that the programming of the hypothalamic circuits regulating energy homeostasis is a key step in the development of obesity associated with malnutrition in early life and provide a valuable resource for further investigating the role of the hypothalamus in the programming of the metabolic syndrome.

Inhibition of GST-pi nuclear transfer increases mantle cell lymphoma sensitivity to cisplatin, cytarabine, gemcitabine, bortezomib and doxorubicin

PURPOSE

Mantle cell lymphoma (MCL) is a chemoresistant lymphoma overexpressing the class pi glutathione-S-transferase (GST-pi). The nuclear localisation of GST-pi is induced by chemotherapy and is correlated to cell resistance. In this study, the effect of the Agaricus bisporus lectin (ABL), a GST-pi nuclear transfer inhibitor, on the chemosensitivity of MCL cells was investigated.

METHODS

The proliferation of three MCL cell lines was evaluated in the presence of doxorubicin (DOX), cisplatin (CDDP), cytarabine (Ara-C), gemcitabine (GEM) or bortezomib with or without ABL pre-treatment.

RESULTS

The cytotoxic activities of CDDP, Ara-C, GEM and bortezomib were increased in all cell lines. The DOX cytotoxic activity was enhanced in two of three cell lines. The inhibition of GST-pi nuclear transfer led to the potentialisation of all drug combinations.

CONCLUSION

The inhibition of the nuclear transfer of GST-pi increases the MCL sensitivity to DOX, CDDP, Ara-C, GEM and bortezomib, alone or in combination.

Expression profiling of prospero in the Drosophila larval chemosensory organ: Between growth and outgrowth

Background

The antenno-maxilary complex (AMC) forms the chemosensory system of the Drosophila larva and is involved in gustatory and olfactory perception. We have previously shown that a mutant allele of the homeodomain transcription factor Prospero (prosVoila1, V1), presents several developmental defects including abnormal growth and altered taste responses. In addition, many neural tracts connecting the AMC to the central nervous system (CNS) were affected. Our earlier reports on larval AMC did not argue in favour of a role of pros in cell fate decision, but strongly suggested that pros could be involved in the control of other aspect of neuronal development. In order to identify these functions, we used microarray analysis of larval AMC and CNS tissue isolated from the wild type, and three other previously characterised prospero alleles, including the V1 mutant, considered as a null allele for the AMC.

Results

A total of 17 samples were first analysed with hierarchical clustering. To determine those genes affected by loss of pros function, we calculated a discriminating score reflecting the differential expression between V1 mutant and other pros alleles. We identified a total of 64 genes in the AMC. Additional manual annotation using all the computed information on the attributed role of these genes in the Drosophila larvae nervous system, enabled us to identify one functional category of potential Prospero target genes known to be involved in neurite outgrowth, synaptic transmission and more specifically in neuronal connectivity remodelling. The second category of genes found to be differentially expressed between the null mutant AMC and the other alleles concerned the development of the sensory organs and more particularly the larval olfactory system. Surprisingly, a third category emerged from our analyses and suggests an association of pros with the genes that regulate autophagy, growth and insulin pathways. Interestingly, EGFR and Notch pathways were represented in all of these three functional categories. We now propose that Pros could perform all of these different functions through the modulation of these two antagonistic and synergic pathways.

Conclusions

The current data contribute to the clarification of the prospero function in the larval AMC and show that pros regulates different function in larvae as compared to those controlled by this gene in embryos. In the future, the possible mechanism by which Pros could achieve its function in the AMC will be explored in detail.

PGC-1-related coactivator modulates mitochondrial-nuclear crosstalk through endogenous nitric oxide in a cellular model of oncocytic thyroid tumours

Background

The PGC-1 related coactivator (PRC), which shares structural and functional features with PGC-1α, is believed to regulate several metabolic pathways as well as mitochondrial biogenesis. Its involvement in the early programming of cell proliferation suggests the existence of finely regulated crosstalk between mitochondrial functions and the cell cycle status.

Methodology/Principal Findings

PRC-regulated pathways were explored in a cell-line model derived from mitochondrial-rich tumours with an essentially oxidative metabolism and specifically high PRC expression. The functional status of mitochondria was compared to the results of microarray analysis under conditions of temporal PRC inhibition. To specify the fine PRC regulation, the expression levels of the genes and proteins involved in the oxidative phosphorylation process were studied by real time quantitative PCR and western blotting. As in earlier studies on PGC-1α, we investigated the role of nitric oxide in PRC-regulated mitochondrial biogenesis and determined its action in the control of the phosphorylation status of the mitogen-activated protein kinase pathway.

Conclusion/Significance

We found that nitric oxide rapidly influences PRC expression at the transcriptional level. Focusing on mitochondrial energetic metabolism, we observed that PRC differentially controls respiratory chain complexes and coupling efficiency in a time-dependent manner to maintain mitochondrial homeostasis. Our results highlight the key role of PRC in the rapid modulation of metabolic functions in response to the status of the cell cycle.

Expression profiling of rainbow trout testis development identifies evolutionary conserved genes involved in spermatogenesis

BACKGROUND

Spermatogenesis is a late developmental process that involves a coordinated expression program in germ cells and a permanent communication between the testicular somatic cells and the germ-line. Current knowledge regarding molecular factors driving male germ cell proliferation and differentiation in vertebrates is still limited and mainly based on existing data from rodents and human. Fish with a marked reproductive cycle and a germ cell development in synchronous cysts have proven to be choice models to study precise stages of the spermatogenetic development and the germ cell-somatic cell communication network. In this study we used 9K cDNA microarrays to investigate the expression profiles underlying testis maturation during the male reproductive cycle of the trout, Oncorhynchus mykiss.

RESULTS

Using total testis samples at various developmental stages and isolated spermatogonia, spermatocytes and spermatids, 3379 differentially expressed trout cDNAs were identified and their gene activation or repression patterns throughout the reproductive cycle were reported. We also performed a tissue-profiling analysis and highlighted many genes for which expression signals were restricted to the testes or gonads from both sexes. The search for orthologous genes in genome-sequenced fish species and the use of their mammalian orthologs allowed us to provide accurate annotations for trout cDNAs. The analysis of the GeneOntology terms therefore validated and broadened our interpretation of expression clusters by highlighting enriched functions that are consistent with known sequential events during male gametogenesis. Furthermore, we compared expression profiles of trout and mouse orthologs and identified a complement of genes for which expression during spermatogenesis was maintained throughout evolution.

CONCLUSION

A comprehensive study of gene expression and associated functions during testis maturation and germ cell differentiation in the rainbow trout is presented. The study identifies new pathways involved during spermatogonia self-renewal or rapid proliferation, meiosis and gamete differentiation, in fish and potentially in all vertebrates. It also provides the necessary basis to further investigate the hormonal and molecular networks that trigger puberty and annual testicular recrudescence in seasonally breeding species.

Regulation of intestinal epithelial cells transcriptome by enteric glial cells: impact on intestinal epithelial barrier functions

BACKGROUND

Emerging evidences suggest that enteric glial cells (EGC), a major constituent of the enteric nervous system (ENS), are key regulators of intestinal epithelial barrier (IEB) functions. Indeed EGC inhibit intestinal epithelial cells (IEC) proliferation and increase IEB paracellular permeability. However, the role of EGC on other important barrier functions and the signalling pathways involved in their effects are currently unknown. To achieve this goal, we aimed at identifying the impact of EGC upon IEC transcriptome by performing microarray studies.

RESULTS

EGC induced significant changes in gene expression profiling of proliferating IEC after 24 hours of co-culture. 116 genes were identified as differentially expressed (70 up-regulated and 46 down-regulated) in IEC cultured with EGC compared to IEC cultured alone. By performing functional analysis of the 116 identified genes using Ingenuity Pathway Analysis, we showed that EGC induced a significant regulation of genes favoring both cell-to-cell and cell-to-matrix adhesion as well as cell differentiation. Consistently, functional studies showed that EGC induced a significant increase in cell adhesion. EGC also regulated genes involved in cell motility towards an enhancement of cell motility. In addition, EGC profoundly modulated expression of genes involved in cell proliferation and cell survival, although no clear functional trend could be identified. Finally, important genes involved in lipid and protein metabolism of epithelial cells were shown to be differentially regulated by EGC.

CONCLUSION

This study reinforces the emerging concept that EGC have major protective effects upon the IEB. EGC have a profound impact upon IEC transcriptome and induce a shift in IEC phenotype towards increased cell adhesion and cell differentiation. This concept needs to be further validated under both physiological and pathophysiological conditions.

Increasing the number of thyroid lesions classes in microarray analysis improves the relevance of diagnostic markers

Background

Genetic markers for thyroid cancers identified by microarray analysis have offered limited predictive accuracy so far because of the few classes of thyroid lesions usually taken into account. To improve diagnostic relevance, we have simultaneously analyzed microarray data from six public datasets covering a total of 347 thyroid tissue samples representing 12 histological classes of follicular lesions and normal thyroid tissue. Our own dataset, containing about half the thyroid tissue samples, included all categories of thyroid lesions.

Methodology/Principal Findings

Classifier predictions were strongly affected by similarities between classes and by the number of classes in the training sets. In each dataset, sample prediction was improved by separating the samples into three groups according to class similarities. The cross-validation of differential genes revealed four clusters with functional enrichments. The analysis of six of these genes (APOD, APOE, CLGN, CRABP1, SDHA and TIMP1) in 49 new samples showed consistent gene and protein profiles with the class similarities observed. Focusing on four subclasses of follicular tumor, we explored the diagnostic potential of 12 selected markers (CASP10, CDH16, CLGN, CRABP1, HMGB2, ALPL2, ADAMTS2, CABIN1, ALDH1A3, USP13, NR2F2, KRTHB5) by real-time quantitative RT-PCR on 32 other new samples. The gene expression profiles of follicular tumors were examined with reference to the mutational status of the Pax8-PPARγ, TSHR, GNAS and NRAS genes.

Conclusion/Significance

We show that diagnostic tools defined on the basis of microarray data are more relevant when a large number of samples and tissue classes are used. Taking into account the relationships between the thyroid tumor pathologies, together with the main biological functions and pathways involved, improved the diagnostic accuracy of the samples. Our approach was particularly relevant for the classification of microfollicular adenomas.

Molecular risk stratification in advanced heart failure patients

Risk stratification in advanced heart failure (HF) is crucial for the individualization of therapeutic strategy, in particular for heart transplantation and ventricular assist device implantation. We tested the hypothesis that cardiac gene expression profiling can distinguish between HF patients with different disease severity. We obtained tissue samples from both left (LV) and right (RV) ventricle of explanted hearts of 44 patients undergoing cardiac transplantation or ventricular assist device placement. Gene expression profiles were obtained using an in-house microarray containing 4217 muscular organ-relevant genes. Based on their clinical status, patients were classified into three HF-severity groups: deteriorating (n= 12), intermediate (n= 19) and stable (n= 13). Two-class statistical analysis of gene expression profiles of deteriorating and stable patients identified a 170-gene and a 129-gene predictor for LV and RV samples, respectively. The LV molecular predictor identified patients with stable and deteriorating status with a sensitivity of 88% and 92%, and a specificity of 100% and 96%, respectively. The RV molecular predictor identified patients with stable and deteriorating status with a sensitivity of 100% and 96%, and a specificity of 100% and 100%, respectively. The molecular prediction was reproducible across biological replicates in LV and RV samples. Gene expression profiling has the potential to reproducibly detect HF patients with highest HF severity with high sensitivity and specificity. In addition, not only LV but also RV samples could be used for molecular risk stratification with similar predictive power.

Genetically determined differences in sodium current characteristics modulate conduction disease severity in mice with cardiac sodium channelopathy

Conduction slowing of the electric impulse that drives the heartbeat may evoke lethal cardiac arrhythmias. Mutations in SCN5A, which encodes the pore-forming cardiac sodium channel alpha subunit, are associated with familial arrhythmia syndromes based on conduction slowing. However, disease severity among mutation carriers is highly variable. We hypothesized that genetic modifiers underlie the variability in conduction slowing and disease severity. With the aim of identifying such modifiers, we studied the Scn5a(1798insD/+) mutation in 2 distinct mouse strains, FVB/N and 129P2. In 129P2 mice, the mutation resulted in more severe conduction slowing particularly in the right ventricle (RV) compared to FVB/N. Pan-genomic mRNA expression profiling in the 2 mouse strains uncovered a drastic reduction in mRNA encoding the sodium channel auxiliary subunit beta4 (Scn4b) in 129P2 mice compared to FVB/N. This corresponded to low to undetectable beta4 protein levels in 129P2 ventricular tissue, whereas abundant beta4 protein was detected in FVB/N. Sodium current measurements in isolated myocytes from the 2 mouse strains indicated that sodium channel activation in myocytes from 129P2 mice occurred at more positive potentials compared to FVB/N. Using computer simulations, this difference in activation kinetics was predicted to explain the observed differences in conduction disease severity between the 2 strains. In conclusion, genetically determined differences in sodium current characteristics on the myocyte level modulate disease severity in cardiac sodium channelopathies. In particular, the sodium channel subunit beta4 (SCN4B) may constitute a potential genetic modifier of conduction and cardiac sodium channel disease.

Molecular characterisation and phylogenetic analysis of Chronic bee paralysis virus, a honey bee virus

The complete sequences of the two major RNAs of Chronic bee paralysis virus (CBPV) have been determined. RNA 1 (3674nt long) and RNA 2 (2305nt long) are positive single-stranded RNAs that are capped but not polyadenylated. The 3' ends of both RNAs are unreactive to polymerisation or ligation even in denaturing conditions, a feature already observed in alphanodavirus RNAs. The three previously described smaller RNAs [Overton, H.A., Buck, K.W., Bailey, L., et al., 1982. Relationships between the RNA components of Chronic bee-paralysis virus and those of chronic bee-paralysis virus associate. J. Gen. Virol. 63, 171-179], were not detected in this study, supporting the hypothesis that they would correspond to the three RNAs of the Chronic bee paralysis satellite virus (CBPSV). RNA 1 and RNA 2 encoded three and four overlapping open reading frames (ORFs), respectively. The amino acid sequences deduced from the ORF 3 on RNA 1 shared the conserved motifs of the RNA-dependent RNA polymerase (RdRp) sequence and presented similarities with members of the Nodaviridae and Tombusviridae families. However, no similarities were found between the other CBPV deduced amino acid sequences and sequences in the NCBI databases, suggesting that CBPV is the prototype of a new family of positive single-stranded RNA viruses.

Expression profiling of candidate genes during ovary-to-testis trans-differentiation in rainbow trout masculinized by androgens

Fish gonadal phenotype is very sensitive to sex steroid and functional masculinizations can be obtained in most species using androgen treatments. To gain insight into the molecular effects of androgen-induced masculinization we characterized, in the rainbow trout, the gonadal expression profiles of 103 candidate genes involved in sex differentiation and early gametogenesis. The androgen treatment (11beta-hydroxyandrostenedione, 10 mg/kg of food for 3 months) was administered in a genetic all-female population. Gonads were sampled at different time points in genetic all-male and all-female control populations and in the androgen-treated group. Gene expression profiles were recorded by real-time RT-PCR and biological samples and gene expressions were compared using a global clustering analysis. This analysis revealed that masculinization with androgens acts firstly by repressing granulosa cell related genes, including genes involved in ovarian differentiation (foxl2a, fst, cyp19a1a), and subsequently by repressing genes important for early oogenesis (gdf9, bcl2lb, fancl, gcl, fshb, lhb, sox23, sox24, nup62 and vtgr). However, this masculinizing treatment did not induce a testicular differentiation similar to what was observed in the control male population. This was especially noticeable for many Leydig cell genes encoding proteins involved in steroidogenesis or its control (hsd3b1, star, cyp17a1, cyp11b2.1 and nr5a1b) that were down-regulated in the androgen-treated group. Concomitantly some Sertoli cells marker genes were up-regulated by the androgen treatment (sox9a.1, nr0b1, cldn11, dmrt1) whereas others were down-regulated (amh, sox9a.2), suggesting a partial differentiation of the Sertoli cell lineage. Overall, this suggests that the crucial step of this masculinization process is the de-differentiation of the granulosa cells.

Androgen-induced masculinization in rainbow trout results in a marked dysregulation of early gonadal gene expression profiles

Background

Fish gonadal sex differentiation is affected by sex steroids treatments providing an efficient strategy to control the sexual phenotype of fish for aquaculture purposes. However, the biological effects of such treatments are poorly understood. The aim of this study was to identify the main effects of an androgen masculinizing treatment (11β-hydroxyandrostenedione, 11βOHΔ4, 10 mg/kg of food for 3 months) on gonadal gene expression profiles of an all-female genetic population of trout. To characterize the most important molecular features of this process, we used a large scale gene expression profiling approach using rainbow trout DNA microarrays combined with a detailed gene ontology (GO) analysis.

Results

2,474 genes were characterized as up-regulated or down-regulated in trout female gonads masculinized by androgen in comparison with control male or female gonads from untreated all-male and all-female genetic populations. These genes were classified in 13 k-means clusters of temporally correlated expression profiles. Gene ontology (GO) data mining revealed that androgen treatment triggers a marked down-regulation of genes potentially involved in early oogenesis processes (GO 'mitotic cell cycle', 'nucleolus'), an up-regulation of the translation machinery (GO 'ribosome') along with a down-regulation of proteolysis (GO 'proteolysis', 'peptidase' and 'metallopeptidase activity'). Genes considered as muscle fibres markers (GO 'muscle contraction') and genes annotated as structural constituents of the extracellular matrix (GO 'extracellular matrix') or related to meiosis (GO 'chromosome' and 'meiosis') were found significantly enriched in the two clusters of genes specifically up-regulated in androgen-treated female gonads. GO annotations 'Sex differentiation' and 'steroid biosynthesis' were enriched in a cluster of genes with high expression levels only in control males. Interestingly none of these genes were stimulated by the masculinizing androgen treatment.

Conclusion

This study provides evidence that androgen masculinization results in a marked dysregulation of early gene expression profiles when compared to natural testicular or ovarian differentiation. Based on these results we suggest that, in our experimental conditions, androgen masculinization proceeds mainly through an early inhibition of female development.

Gene expression profiling of human melanoma cell lines with distinct metastatic potential identifies new progression markers

BACKGROUND

Identification of markers associated with melanoma progression is crucial to identify new prognostic and/or therapeutic targets.

MATERIALS AND METHODS

By using DNA microarrays, two human melanoma cell lines, M4Be and Tw12, derived from the same tumor, but with different metastatic potential, were profiled. Western blot of cell lines, immunohistochemistry on melanoma biopsies and in silico analyses validated and extended our results.

RESULTS

Thirty-six clones were differentially-expressed between the two cell lines, representing 33 named genes and 2 expressed sequence tags. The most up-regulated gene in the strongly metastatic clone Tw12 was CD10. Protein analysis with anti-CD1O antibody confirmed this finding in cell lines and clinical samples with expression being more frequent in metastatic compared to primary tumors. Many up-regulated genes were involved in angiogenesis, invasion, growth and apoptosis. Down-regulated genes included tumor suppressor genes and those were involved in differentiation.

CONCLUSION

We identified several genes the expression of which is associated with metastatic progression in human melanoma cells. Although further analyses are warranted to clarify their exact role in tumor progression, they might lead to new prognostic markers and/or molecular therapeutic targets in metastatic melanoma.

Gene expression profiling identifies molecular subgroups among nodal peripheral T-cell lymphomas

The classification of peripheral T-cell lymphomas (PTCL) is still a matter of debate. To establish a molecular classification of PTCL, we analysed 59 primary nodal T-cell lymphomas using cDNA microarrays, including 56 PTCL and three T-lymphoblastic lymphoma (T-LBL). The expression profiles could discriminate angioimmunoblastic lymphoma, anaplastic large-cell lymphoma and T-LBL. In contrast, cases belonging to the broad category of 'PTCL, unspecified' (PTCL-U) did not share a single molecular profile. Using a multiclass predictor, we could separate PTCL-U into three molecular subgroups called U1, U2 and U3. The U1 gene expression signature included genes known to be associated with poor outcome in other tumors, such as CCND2. The U2 subgroup was associated with overexpression of genes involved in T-cell activation and apoptosis, including NFKB1 and BCL-2. The U3 subgroup was mainly defined by overexpression of genes involved in the IFN/JAK/STAT pathway. It comprised a majority of histiocyte-rich PTCL samples. Gene Ontology annotations revealed different functional profile for each subgroup. These results suggest the existence of distinct subtypes of PTCL-U with specific molecular profiles, and thus provide a basis to improve their classification and to develop new therapeutic targets.

Gene profiling reveals specific oncogenic mechanisms and signaling pathways in oncocytic and papillary thyroid carcinoma

The oncogenic pathways in mitochondrial-rich thyroid carcinomas are not clearly understood. To investigate the possible implication of mitochondrial abundance in the genesis of thyroid tumors, we have explored the gene expression profile of six oncocytic carcinomas and six mitochondrial-rich papillary carcinomas using cDNA-microarray technology. A supervised approach allowed us to identify 83 genes differentially expressed in the two types of carcinoma. These genes were classified according to their ontologic profiles. Three genes, NOS3, alpha-actinin-2 and alpha-catenin, suspected of playing a role in tumor genesis, were explored by quantitative RT-PCR analysis and immunohistochemistry. Of the 59 genes overexpressed in papillary carcinomas, 51% were involved in cell communication. Of the 24 genes overexpressed in oncocytic carcinomas, 84% were involved in mitochondrial and cellular metabolism. Our results suggest that mitochondrial respiratory chain complexes III and IV play a significant role in the regulation of reactive oxygen species production by oncocytic tumors.

Gene Expression Profiling Identifies Molecular Subtypes of Inflammatory Breast Cancer

Breast cancer is a heterogeneous disease. Comprehensive gene expression profiles obtained using DNA microarrays have revealed previously indistinguishable subtypes of noninflammatory breast cancer (NIBC) related to different features of mammary epithelial biology and significantly associated with survival. Inflammatory breast cancer (IBC) is a rare, particular, and aggressive form of disease. Here we have investigated whether the five molecular subtypes described for NIBC (luminal A and B, basal, ERBB2 overexpressing, and normal breast-like) were also present in IBC. We monitored the RNA expression of approximately 8,000 genes in 83 breast tissue samples including 37 IBC, 44 NIBC, and 2 normal breast samples. Hierarchical clustering identified the five subtypes of breast cancer in both NIBC and IBC samples. These subtypes were highly similar to those defined in previous studies and associated with similar histoclinical features. The robustness of this classification was confirmed by the use of both alternative gene set and analysis method, and the results were corroborated at the protein level. Furthermore, we show that the differences in gene expression between NIBC and IBC and between IBC with and without pathologic complete response that we have recently reported persist in each subtype. Our results show that the expression signatures defining molecular subtypes of NIBC are also present in IBC. Obtained using different patient series and different microarray platforms, they reinforce confidence in the expression-based molecular taxonomy but also give evidence for its universality in breast cancer, independently of a specific clinical form.

Identification of new classes among acute myelogenous leukaemias with normal karyotype using gene expression profiling

Conventional cytogenetic analysis currently stratifies acute myelogenous leukaemia (AML) into prognostically relevant groups. However, approximately 50% of adult AMLs have normal cytogenetics (NC-AMLs), and represent a heterogeneous and poorly understood group. We analysed gene expression in 55 AML samples including 53 cases from adult patients with NC-AML (n = 36), trisomy 8, t(15;17), t(8;21), t(11;19), 7q deletion, and two cell lines using 9000-gene DNA microarrays. Global hierarchical clustering showed that NC-AMLs are a heterogeneous group. Supervised analysis distinguished two subgroups of NC-AML: one subgroup constituted a homogeneous NC cluster ('pure NC-AML'), and the other NC-AMLs were close to the AML cases with translocations ('translocation like'). Gene expression signatures were also derived for patients with trisomy 8, as well as FLT3 and MLL gene duplications. Importantly, samples from 24 NC-AML patients who could be evaluated for clinical outcome were analysed. In all, 43 genes that discriminated two classes of patients with significantly different prognosis were identified. The poor prognosis class contained a majority of 'pure NC-AMLs', whereas the 'translocation-like' AMLs were in the good prognosis class. Discriminator genes included genes involved in drug resistance (TOP2B), protein transport (MTX2, SLC35A2), and cell signalling (MAPK1, PRKAB2). Our results demonstrate the transcriptional heterogeneity of NC-AMLs, and suggest the existence of 'translocation-like' NC-AMLs and of a gene expression signature that may predict response to chemotherapy.

Gene expression profiling of human adrenocortical tumors using complementary deoxyribonucleic Acid microarrays identifies several candidate genes as markers of malignancy

The aim of this study was to identify predictor sets of genes whose over- or underexpression in human sporadic adrenocortical tumors would help to identify malignant vs. benign tumors and to predict postsurgical metastatic recurrence. For this, we analyzed the expression of 230 candidate genes using cDNA microarrays in a series of 57 well-characterized human sporadic adrenocortical tumors (33 adenomas and 24 carcinomas). We identified two clusters of genes (the IGF-II cluster containing eight genes, including IGF-II, and the steroidogenesis cluster containing six genes encoding steroidogenic enzymes plus eight other genes) whose combined levels of expression appeared to be good predictors of malignancy. This predictive value was as strong as that of the pathological score of Weiss. The analysis of the population of carcinomas (13 tumors) for genes whose expression would be strongly different between recurring and nonrecurring tumors allowed identification of 14 genes meeting these criteria. Among these genes, there are probably new markers of tumor evolution that will deserve additional validation on a larger scale. Taken together, these results show that the parallel analysis of the expression levels of a selected group of genes on microgram quantities of tumor RNA (a quantity that can be obtained from fine needle aspirations) appears as a complementary method to histopathology for the diagnosis and prognosis of evolution of adrenocortical carcinomas.

Gene expression profiling for molecular characterization of inflammatory breast cancer and prediction of response to chemotherapy

Inflammatory breast cancer (IBC) is a rare but aggressive form of breast cancer with a 5-year survival limited to approximately 40%. Diagnosis, based on clinical and/or pathological criteria, may be difficult. Optimal systemic neoadjuvant therapy and accurate predictors of pathological response have yet to be defined for increasing response rate and survival. Using DNA microarrrays containing approximately 8,000 genes, we profiled breast cancer samples from 81 patients, including 37 with IBC and 44 with noninflammatory breast cancer (NIBC). Global unsupervised hierarchical clustering was able to some extent to distinguish IBC and NIBC cases and revealed subclasses of IBC. Supervised analysis identified a 109-gene set the expression of which discriminated IBC from NIBC samples. This molecular signature was validated in an independent series of 26 samples, with an overall performance accuracy of 85%. Discriminator genes were associated with various cellular processes possibly related to the aggressiveness of IBC, including signal transduction, cell motility, adhesion, and angiogenesis. A similar approach, with leave-one-out cross-validation, identified an 85-gene set that divided IBC patients with significantly different pathological complete response rate (70% in one group and 0% in the other group). These results show the potential of gene expression profiling to contribute to a better understanding of IBC, and to provide new diagnostic and predictive factors for IBC, as well as for potential therapeutic targets.

Feature extraction and signal processing for nylon DNA microarrays

Background

High-density DNA microarrays require automatic feature extraction methodologies and softwares. These can be a potential source of non-reproducibility of gene expression measurements. Variation in feature location or in signal integration methodology may be a significant contribution to the observed variance in gene expression levels.

Results

We explore sources of variability in feature extraction from DNA microarrays on Nylon membrane with radioactive detection. We introduce a mathematical model of the signal emission and derive methods for correcting biases such as overshining, saturation or variation in probe amount. We also provide a quality metric which can be used qualitatively to flag weak or untrusted signals or quantitatively to modulate the weight of each experiment or gene in higher level analyses (clustering or discriminant analysis).

Conclusions

Our novel feature extraction methodology, based on a mathematical model of the radioactive emission, reduces variability due to saturation, neighbourhood effects and variable probe amount. Furthermore, we provide a fully automatic feature extraction software, BZScan, which implements the algorithms described in this paper.

Small lymphocytic lymphoma, marginal zone B-cell lymphoma, and mantle cell lymphoma exhibit distinct gene-expression profiles allowing molecular diagnosis

Non-germinal center small B-cell lymphomas represent a heterogeneous group of non-Hodgkin lymphomas, the most frequent histologic subtypes being small lymphocytic lymphoma (SLL), splenic marginal zone B-cell lymphoma (MZL), and mantle cell lymphoma (MCL). In order to identify genomic signatures specific for each disease, we analyzed 128 primary tumors using high-density microarrays. Several clusters of genes significantly discriminated the 3 histologic subtypes. Genes associated with cell adhesion, angiogenesis, and inhibition of apoptosis were up-regulated in SLL. Genes associated with intracellular signaling via the AKT1 pathway were up-regulated in splenic MZL. Genes associated with cell cycle control and multidrug resistance were up-regulated in MCL. Using 44 genes selected within the gene clusters discriminant for the 3 lymphoma subtypes, we generated a class prediction score that allowed us to classify the 3 entities in 96% of the cases, including borderline cases. Whereas specific transcriptional profiles easily distinguished all MZL samples, SLL samples, and most of the MCL samples into separate groups, few MCL cases exhibited MZL-type transcriptional profiles. This study demonstrates that SLL, splenic MZL, and MCL possess specific transcriptional profiles that may be relevant to the pathogenesis and the diagnosis of these histologic subtypes.

Gene expression profiling of colon cancer by DNA microarrays and correlation with histoclinical parameters

Different diagnostic and prognostic groups of colorectal carcinoma (CRC) have been defined. However, accurate diagnosis and prediction of survival are sometimes difficult. Gene expression profiling might improve these classifications and bring new insights into underlying molecular mechanisms. We profiled 50 cancerous and noncancerous colon tissues using DNA microarrrays consisting of approximately 8000 spotted human cDNA. Global hierarchical clustering was to some extent able to distinguish clinically relevant subgroups, normal versus cancer tissues and metastatic versus nonmetastatic tumours. Supervised analyses improved these segregations by identifying sets of genes that discriminated between normal and tumour tissues, tumours associated or not with lymph node invasion or genetic instability, and tumours from the right or left colon. A similar approach identified a gene set that divided patients with significantly different 5-year survival (100% in one group and 40% in the other group; P=0.005). Discriminator genes were associated with various cellular processes. An immunohistochemical study on 382 tumour and normal samples deposited onto a tissue microarray subsequently validated the upregulation of NM23 in CRC and a downregulation in poor prognosis tumours. These results suggest that microarrays may provide means to improve the classification of CRC, provide new potential targets against carcinogenesis and new diagnostic and/or prognostic markers and therapeutic targets.