Genes that code for T cell signaling proteins establish transcriptional regulatory networks during thymus ontogeny

Gene expression profiling by cDNA microarrays during murine thymus ontogeny has contributed to dissecting the large-scale molecular genetics of T cell maturation. Gene profiling, although useful for characterizing the thymus developmental phases and identifying the differentially expressed genes, does not permit the determination of possible interactions between genes. In order to reconstruct genetic interactions, on RNA level, within thymocyte differentiation, a pair of microarrays containing a total of 1,576 cDNA sequences derived from the IMAGE MTB library was applied on samples of developing thymuses (14-17 days of gestation). The data were analyzed using the GeneNetwork program. Genes that were previously identified as differentially expressed during thymus ontogeny showed their relationships with several other genes. The present method provided the detection of gene nodes coding for proteins implicated in the calcium signaling pathway, such as Prrg2 and Stxbp3, and in protein transport toward the cell membrane, such as Gosr2. The results demonstrate the feasibility of reconstructing networks based on cDNA microarray gene expression determinations, contributing to a clearer understanding of the complex interactions between genes involved in thymus/thymocyte development.

Arrow Plot: a new graphical tool for selecting up and down regulated genes and genes differentially expressed on samples subgroups

Background: A common task in analyzing microarray data is to determine which genes are differentially expressed across two (or more) kind of tissue samples or samples submitted under experimental conditions. Several statistical methods have been proposed to accomplish this goal, generally based on measures of distance between classes. It is well known that biological samples are heterogeneous because of factors such as molecular subtypes or genetic background that are often unknown to the experimenter. For instance, in experiments which involve molecular classification of tumors it is important to identify significant subtypes of cancer. Bimodal or multimodal distributions often reflect the presence of subsamples mixtures. Consequently, there can be genes differentially expressed on sample subgroups which are missed if usual statistical approaches are used. In this paper we propose a new graphical tool which not only identifies genes with up and down regulations, but also genes with differential expression in different subclasses, that are usually missed if current statistical methods are used. This tool is based on two measures of distance between samples, namely the overlapping coefficient (OVL) between two densities and the area under the receiver operating characteristic (ROC) curve. The methodology proposed here was implemented in the open-source R software. Results: This method was applied to a publicly available dataset, as well as to a simulated dataset. We compared our results with the ones obtained using some of the standard methods for detecting differentially expressed genes, namely Welch t-statistic, fold change (FC), rank products (RP), average difference (AD), weighted average difference (WAD), moderated t-statistic (modT), intensity-based moderated t-statistic (ibmT), significance analysis of microarrays (samT) and area under the ROC curve (AUC). On both datasets all differentially expressed genes with bimodal or multimodal distributions were not selected by all standard selection procedures. We also compared our results with (i) area between ROC curve and rising area (ABCR) and (ii) the test for not proper ROC curves (TNRC). We found our methodology more comprehensive, because it detects both bimodal and multimodal distributions and different variances can be considered on both samples. Another advantage of our method is that we can analyze graphically the behavior of different kinds of differentially expressed genes. Conclusion: Our results indicate that the arrow plot represents a new flexible and useful tool for the analysis of gene expression profiles from microarrays.

Inheritance of histone H3 variants across mitotic cell divisions

Dissertation presented to obtain the Ph.D degree in Molecular Biology.

Genomic, Epigenomic and Transcriptomic Molecular Characterization of Splenic Marginal Zone Lymphoma Reveals Recurrent Abnormalities in miR-182

Splenic marginal zone lymphoma (SMZL) is a low grade B-cell non-Hodgkin's lymphoma. The molecular pathology of this entity remains poorly understood. To characterise this lymphoma at the molecular level, we performed an integrated analysis of 1) genome wide genetic copy number alterations 2) gene expression profiles and 3) epigenetic DNA methylation profiles.We have previously shown that SMZL is characterised by recurrent alterations of chromosomes 7q, 6q, 3q, 9q and 18; however, gene resolution oligonucleotide array comparative genomic hybridisation did not reveal evidence of cryptic amplification or deletion in these regions. The most frequently lost 7q32 region contains a cluster of miRNAs. qRT-PCR revealed that three of these (miR-182/96/183) show underexpression in SMZL, and miR-182 is somatically mutated in >20% of cases of SMZL, as well as in >20% of cases of follicular lymphoma, and between 5-15% of cases of chronic lymphocytic leukaemia, MALT-lymphoma and hairy cell leukaemia. We conclude that miR-182 is a strong candidate novel tumour suppressor miRNA in lymphoma.The overall gene expression signature of SMZL was found to be strongly distinct fromthose of other lymphomas. Functional analysis of gene expression data revealed SMZL to be characterised by abnormalities in B-cell receptor signalling (especially through the CD19/21-PI3K/AKT pathway) and apoptotic pathways. In addition, genes involved in the response to viral infection appeared upregulated. SMZL shows a unique epigenetic profile, but analysis of differentially methylated genes showed few with methylation related transcriptional deregulation, suggesting that DNA methylation abnormalities are not a critical component of the SMZL malignant phenotype.

Genetics of Parkinson disease and essential tremor.

Purpose of review: Elucidating the genetic background of Parkinson disease and essential tremor is crucial to understand the pathogenesis and improve diagnostic and therapeutic strategies. Recent findings: A number of approaches have been applied including familial and association studies, and studies of gene expression profiles to identify genes involved in susceptibility to Parkinson disease. These studies have nominated a number of candidate Parkinson disease genes and novel loci including Omi/HtrA2, GIGYF2, FGF20, PDXK, EIF4G1 and PARK16. A recent notable finding has been the confirmation for the role of heterozygous mutations in glucocerebrosidase (GBA) as risk factors for Parkinson disease. Finally, association studies have nominated genetic variation in the leucine-rich repeat and Ig containing 1 gene (LINGO1) as a risk for both Parkinson disease and essential tremor, providing the first genetic evidence of a link between the two conditions. Summary: Although undoubtedly genes remain to be identified, considerable progress has been achieved in the understanding of the genetic basis of Parkinson disease. This same effort is now required for essential tremor. The use of next-generation high-throughput sequencing and genotyping technologies will help pave the way for future insight leading to advances in diagnosis, prevention and cure.

In vivo conditional Pax4 overexpression in mature islet β-cells prevents stress-induced hyperglycemia in mice.

OBJECTIVE To establish the role of the transcription factor Pax4 in pancreatic islet expansion and survival in response to physiological stress and its impact on glucose metabolism, we generated transgenic mice conditionally and selectively overexpressing Pax4 or a diabetes-linked mutant variant (Pax4R129 W) in β-cells. RESEARCH DESIGN AND METHODS Glucose homeostasis and β-cell death and proliferation were assessed in Pax4- or Pax4R129 W-overexpressing transgenic animals challenged with or without streptozotocin. Isolated transgenic islets were also exposed to cytokines, and apoptosis was evaluated by DNA fragmentation or cytochrome C release. The expression profiles of proliferation and apoptotic genes and β-cell markers were studied by immunohistochemistry and quantitative RT-PCR. RESULTS Pax4 but not Pax4R129 W protected animals against streptozotocin-induced hyperglycemia and isolated islets from cytokine-mediated β-cell apoptosis. Cytochrome C release was abrogated in Pax4 islets treated with cytokines. Interleukin-1β transcript levels were suppressed in Pax4 islets, whereas they were increased along with NOS2 in Pax4R129 W islets. Bcl-2, Cdk4, and c-myc expression levels were increased in Pax4 islets while MafA, insulin, and GLUT2 transcript levels were suppressed in both animal models. Long-term Pax4 expression promoted proliferation of a Pdx1-positive cell subpopulation while impeding insulin secretion. Suppression of Pax4 rescued this defect with a concomitant increase in pancreatic insulin content. CONCLUSIONS Pax4 protects adult islets from stress-induced apoptosis by suppressing selective nuclear factor-κB target genes while increasing Bcl-2 levels. Furthermore, it promotes dedifferentiation and proliferation of β-cells through MafA repression, with a concomitant increase in Cdk4 and c-myc expression.

IQRray, a new method for Affymetrix microarray quality control, and the homologous organ conservation score, a new benchmark method for quality control metrics.

MOTIVATION: Microarray results accumulated in public repositories are widely reused in meta-analytical studies and secondary databases. The quality of the data obtained with this technology varies from experiment to experiment, and an efficient method for quality assessment is necessary to ensure their reliability. RESULTS: The lack of a good benchmark has hampered evaluation of existing methods for quality control. In this study, we propose a new independent quality metric that is based on evolutionary conservation of expression profiles. We show, using 11 large organ-specific datasets, that IQRray, a new quality metrics developed by us, exhibits the highest correlation with this reference metric, among 14 metrics tested. IQRray outperforms other methods in identification of poor quality arrays in datasets composed of arrays from many independent experiments. In contrast, the performance of methods designed for detecting outliers in a single experiment like Normalized Unscaled Standard Error and Relative Log Expression was low because of the inability of these methods to detect datasets containing only low-quality arrays and because the scores cannot be directly compared between experiments. AVAILABILITY AND IMPLEMENTATION: The R implementation of IQRray is available at: ftp://lausanne.isb-sib.ch/pub/databases/Bgee/general/IQRray.R. CONTACT: Marta.Rosikiewicz@unil.ch SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Aging of myelinating glial cells predominantly affects lipid metabolism and immune response pathways.

Both the central and the peripheral nervous systems are prone to multiple age-dependent neurological deficits, often attributed to still unknown alterations in the function of myelinating glia. To uncover the biological processes affected in glial cells by aging, we analyzed gene expression of the Schwann cell-rich mouse sciatic nerve at 17 time points throughout life, from day of birth until senescence. By combining these data with the gene expression data of myelin mouse mutants carrying deletions of either Pmp22, SCAP, or Lpin1, we found that the majority of age-related transcripts were also affected in myelin mutants (54.4%) and were regulated during PNS development (59.5%), indicating a high level of overlap in implicated molecular pathways. The expression profiles in aging copied the direction of transcriptional changes observed in neuropathy models; however, they had the opposite direction when compared with PNS development. The most significantly altered biological processes in aging involved the inflammatory/immune response and lipid metabolism. Interestingly, both these pathways were comparably changed in the aging optic nerve, suggesting that similar biological processes are affected in aging of glia-rich parts of the central and peripheral nervous systems. Our comprehensive comparison of gene expression in three distinct biological conditions including development, aging, and myelin disease thus revealed a previously unanticipated relationship among themselves and identified lipid metabolism and inflammatory/immune response pathways as potential therapeutical targets to prevent or delay so far incurable age-related and inherited forms of neuropathies.

Bone anabolic potential of soy isoflavones and plant extracts and genomic changes during differentiation of hPOB-tert osteoblast-like cells

Summary For the nutritional management of bone health and the prevention of osteoporosis it is important to identify nutrients that positively influence the bone remodeling process at the cellular level. Soy isoflavones show promising osteoprotective effects in animals and humans but their mechanism of action in bone cells is yet poorly understood. Firstly, soy tissue cultures were characterized as a new and optimized source of isoflavones. A large variability in the isoflavone content was observed and high-producing strains (46.3 mg/g dry wt isoflavones) were identified. In the Ishikawa cells bioassay, the estrogenicity of isoflavones was confirmed to be 1000 to 10000 less than 17Mestradiol and that of the malonyl forms was shown for the first time (EC50 of 350 nM and 1880 nM for malonylgenistin and malonyldaidzin, respectively). The estrogenic activity of soya tissue culture extracts correlated to their isoflavone content. Secondly, the effects of phytonutrients on BMP-2 gene expression and on the mevalonate synthesis pathway, as key mediators of bone formation, were investigated. Dietary achievable concentrations of genistein and daidzein (10vM), and statins (4xM) but not 17M estradiol (10nM), induced BMP-2 gene expression (by up to 3-fold) and inhibited the cholesterol biosynthetic pathway (by up to 50%) in the human osteoblastic cell line hP0B¬tert. In addition, several plant extracts (Cyperus rotundus, Lindera benzoin and Cnidium monnieri) induced BMP-2 gene expression but this induction was not restricted to the inhibition of the cholesterol synthesis pathway neither to the estrogenicity. Finally, the gene expression profiles during hP0B-tert differentiation induced by vitamin D and dexamethasone were analyzed with the Affymetrix human GeneChip. 1665 different genes and 98 ESTs were significantly regulated. The expression profiles of bone-related genes was largely in agreement with previously documented patterns, supporting the physiological relevance of the genomic results and the hP0B-tert cell line as a valid model of human osteoblast differentiation. The expression of alternative differentiation markers during the osteogenic treatment of hP0B-tert cells indicated that the adipocyte and myoblast differentiation pathways were repressed, confirming that these culture conditions allowed only osteoblast differentiation. The gene ontology analysis identified further sub-groups of genes that may be involved in the bone formation process. Aims of the thesis In order to define new strategies for the nutritional management of bone health and for the prevention of osteoporosis the major goal of the present work was to investigate the potential of phytonutrients to positively modulate the bone formation process at the cellular level and, in particular: 1.To select and optimise alternative plant sources containing high levels of isoflavones with estrogenic activity (Chapter 3). 2.To compare the effects of statins and phytonutrients on BMP-2 gene expression and on the mevalonate synthesis pathway and to select new plant extracts with a bone anabolic potential (Chapter 4). 3.To further characterize the new human periosteal cell line, hP0B-tert, as a bone- formation model, by elucidating its gene expression profile during differentiation induced by vitamin D and dexamethasone (Chapter 5).

Cellular origin of T-cell lymphomas.

In this issue of Blood, Iqbal et al, having compiled gene expression profiles from >300 peripheral T-cell lymphomas, expand previous findings on the diagnostic value of molecular signatures that correlate with different histological types of T-cell lymphomas. They report the discovery of 2 molecular subgroups of peripheral T-cell lymphomas, not otherwise specified (PTCL, NOS), characterized by high expression of either GATA-binding protein 3 (GATA-3) or t-box 21 (TBX21) transcription factors and corresponding target genes, with the GATA3 subgroup being associated with distinctly worse prognosis. In an independent study, Wang et al(2) also show that GATA3 expression in a subset of PTCL, NOS identifies a subgroup of patients with inferior survival.

Cloning and characterization of a V-ATPase subunit C from the American visceral leishmaniasis vector Lutzomyia longipalpis modulated during development and blood ingestion

Visceral leishmaniasis (VL) is a serious tropical disease that affects approximately 500 thousand people worldwide every year. In the Americas, VL is caused by the parasite Leishmania (Leishmania) infantum chagasi mainly transmitted by the bite of the sand fly vector Lutzomyia longipalpis. Despite recent advances in the study of interaction between Leishmania and sand flies, very little is known about sand fly protein expression profiles. Understanding how the expression of proteins may be affected by blood feeding and/or presence of parasite in the vector's midgut might allow us to devise new strategies for controlling the spread of leishmaniasis. In this work, we report the characterization of a vacuolar ATPase subunit C from L. longipalpis by screening of a midgut cDNA library with a 220 bp fragment identified by means of differential display reverse transcriptase-polymerase chain reaction analysis. The expression of the gene varies along insect development and is upregulated in males and bloodfed L. longipalpis, compared to unfed flies.

Changes in the transcriptional profile of transporters in the intestine along the anterior-posterior and crypt-villus axes.

BACKGROUND: The purpose of this work was to characterize the expression of drug and nutrient carriers along the anterior-posterior and crypt-villus axes of the intestinal epithelium and to study the validity of utilizing whole gut tissue rather than purified epithelial cells to examine regional variations in gene expression. RESULTS: We have characterized the mRNA expression profiles of 76 % of all currently known transporters along the anterior-posterior axis of the gut. This is the first study to describe the expression profiles of the majority of all known transporters in the intestine. The expression profiles of transporters, as defined according to the Gene Ontology consortium, were measured in whole tissue of the murine duodenum, jejunum, ileum and colon using high-density microarrays. For nine transporters (Abca1, Abcc1, Abcc3, Abcg8, Slc10a2, Slc28a2, Slc2a1, Slc34a2 and Slc5a8), the mRNA profiles were further measured by RT-PCR in laser micro-dissected crypt and villus epithelial cells corresponding to the aforementioned intestinal regions. With respect to differentially regulated transporters, the colon had a distinct expression profile from small intestinal segments. The majority (59 % for p cutoff < or = 0.05) of transporter mRNA levels were constant across the intestinal sections studied. For the transporter subclass "carrier activity", which contains the majority of known carriers for biologically active compounds, a significant change (p < or = 0.05) along the anterior-posterior axis was observed. CONCLUSION: All nine transporters examined in laser-dissected material demonstrated good replication of the region-specific profiles revealed by microarray. Furthermore, we suggest that the distribution characteristics of Slc5a8 along the intestinal tract render it a suitable candidate carrier for monocarboxylate drugs in the posterior portion of the intestine. Our findings also predict that there is a significant difference in the absorption of carrier-mediated compounds in the different intestinal segments. The most pronounced differences can be expected between the adjoining segments ileum and colon, but the differences between the other adjoining segments are not negligible. Finally, for the examined genes, profiles measured in whole intestinal tissue extracts are representative of epithelial cell-only gene expression.

Species differences in the response of liver drug-metabolizing enzymes to (S)-4-O-tolylsulfanyl-2-(4-trifluormethyl-phenoxy)-butyric acid (EMD 392949) in vivo and in vitro.

Induction of drug-metabolizing enzymes (DMEs) is highly species-specific and can lead to drug-drug interaction and toxicities. In this series of studies we tested the species specificity of the antidiabetic drug development candidate and mixed peroxisome proliferator-activated receptor (PPAR) alpha/gamma agonist (S)-4-O-tolylsulfanyl-2-(4-trifluormethyl-phenoxy)-butyric acid (EMD 392949, EMD) with regard to the induction of gene expression and activities of DMEs, their regulators, and typical PPAR target genes. EMD clearly induced PPARalpha target genes in rats in vivo and in rat hepatocytes but lacked significant induction of DMEs, except for cytochrome P450 (P450) 4A. CYP2C and CYP3A were consistently induced in livers of EMD-treated monkeys. Interestingly, classic rodent peroxisomal proliferation markers were induced in monkeys after 17 weeks but not after a 4-week treatment, a fact also observed in human hepatocytes after 72 h but not 24 h of EMD treatment. In human hepatocyte cultures, EMD showed similar gene expression profiles and induction of P450 activities as in monkeys, indicating that the monkey is predictive for human P450 induction by EMD. In addition, EMD induced a similar gene expression pattern as the PPARalpha agonist fenofibrate in primary rat and human hepatocyte cultures. In conclusion, these data showed an excellent correlation of in vivo data on DME gene expression and activity levels with results generated in hepatocyte monolayer cultures, enabling a solid estimation of human P450 induction. This study also clearly highlighted major differences between primates and rodents in the regulation of major inducible P450s, with evidence of CYP3A and CYP2C inducibility by PPARalpha agonists in monkeys and humans.

Insights into the transcriptome of oenocytes from Aedes aegypti pupae

Oenocytes are ectodermic cells present in the fat body of several insect species and these cells are considered to be analogous to the mammalian liver, based on their role in lipid storage, metabolism and secretion. Although oenocytes were identified over a century ago, little is known about their messenger RNA expression profiles. In this study, we investigated the transcriptome of Aedes aegypti oenocytes. We constructed a cDNA library from Ae. aegypti MOYO-R strain oenocytes collected from pupae and randomly sequenced 687 clones. After sequences editing and assembly, 326 high-quality contigs were generated. The most abundant transcripts identified corresponded to the cytochrome P450 superfamily, whose members have roles primarily related to detoxification and lipid metabolism. In addition, we identified 18 other transcripts with putative functions associated with lipid metabolism. One such transcript, a fatty acid synthase, is highly represented in the cDNA library of oenocytes. Moreover, oenocytes expressed several immunity-related genes and the majority of these genes were lysozymes. The transcriptional profile suggests that oenocytes play diverse roles, such as detoxification and lipid metabolism, and increase our understanding of the importance of oenocytes in Ae. aegypti homeostasis and immune competence.

Mapping the proteome of Leishmania Viannia parasites using two-dimensional polyacrylamide gel electrophoresis and associated technologies.

In this study we have demonstrated the potential of two-dimensional electrophoresis (2DE)-based technologies as tools for characterization of the Leishmania proteome (the expressed protein complement of the genome). Standardized neutral range (pH 5-7) proteome maps of Leishmania (Viannia) guyanensis and Leishmania (Viannia) panamensis promastigotes were reproducibly generated by 2DE of soluble parasite extracts, which were prepared using lysis buffer containing urea and nonidet P-40 detergent. The Coomassie blue and silver nitrate staining systems both yielded good resolution and representation of protein spots, enabling the detection of approximately 800 and 1,500 distinct proteins, respectively. Several reference protein spots common to the proteomes of all parasite species/strains studied were isolated and identified by peptide mass spectrometry (LC-ES-MS/MS), and bioinformatics approaches as members of the heat shock protein family, ribosomal protein S12, kinetoplast membrane protein 11 and a hypothetical Leishmania-specific 13 kDa protein of unknown function. Immunoblotting of Leishmania protein maps using a monoclonal antibody resulted in the specific detection of the 81.4 kDa and 77.5 kDa subunits of paraflagellar rod proteins 1 and 2, respectively. Moreover, differences in protein expression profiles between distinct parasite clones were reproducibly detected through comparative proteome analyses of paired maps using image analysis software. These data illustrate the resolving power of 2DE-based proteome analysis. The production and basic characterization of good quality Leishmania proteome maps provides an essential first step towards comparative protein expression studies aimed at identifying the molecular determinants of parasite drug resistance and virulence, as well as discovering new drug and vaccine targets.