Using biplots to interpret gene expression patterns in plants

Sum: Plant biologists in fields of ecology, evolution, genetics and breeding frequently use multivariate methods. This paper illustrates Principal Component Analysis (PCA) and Gabriel's biplot as applied to microarray expression data from plant pathology experiments. Availability: An example program in the publicly distributed statistical language R is available from the web site (www.tpp.uq.edu.au) and by e-mail from the contact. Contact: scott.chapman@csiro.au.

Rhipicephalus (Boophilus) microplus: Clotting time in tick-infested skin varies according to local inflammation and gene expression patterns in tick salivary glands

Ticks deposit saliva at the site of their attachment to a host in order to inhibit haemostasis, inflammation and innate and adaptive immune responses. The anti-haemostatic properties of tick saliva have been described by many studies, but few show that tick infestations or its anti-haemostatic components exert systemic effects in vivo. In the present study, we extended these observations and show that, compared with normal skin, bovine hosts that are genetically susceptible to tick infestations present an increase in the clotting time of blood collected from the immediate vicinity of haemorrhagic feeding pools in skin infested with different developmental stages of Rhipicepahlus microplus; conversely, we determined that clotting time of tick-infested skin from genetically resistant bovines was shorter than that of normal skin. Coagulation and inflammation have many components in common and we determined that in resistant bovines, eosinophils and basophils, which are known to contain tissue factor, are recruited in greater numbers to the inflammatory site of tick bites than in susceptible hosts. Finally, we correlated the observed differences in clotting times with the expression profiles of transcripts for putative anti-haemostatic proteins in different developmental stages of R. microplus fed on genetically susceptible and resistant hosts: we determined that transcripts coding for proteins similar to these molecules are overrepresented in salivary glands from nymphs and males fed on susceptible bovines. Our data indicate that ticks are able to modulate their host`s local haemostatic reactions. In the resistant phenotype, larger amounts of inflammatory cells are recruited and expression of anti-coagulant molecules is decreased tick salivary glands, features that can hamper the tick`s blood meal. (C) 2010 Elsevier Inc. All rights reserved.

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

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

Effects of ploidy and sex-locus genotype on gene expression patterns in the fire ant Solenopsis invicta.

Males in many animal species differ greatly from females in morphology, physiology and behaviour. Ants, bees and wasps have a haplodiploid mechanism of sex determination whereby unfertilized eggs become males while fertilized eggs become females. However, many species also have a low frequency of diploid males, which are thought to develop from diploid eggs when individuals are homozygous at one or more sex determination loci. Diploid males are morphologically similar to haploids, though often larger and typically sterile. To determine how ploidy level and sex-locus genotype affect gene expression during development, we compared expression patterns between diploid males, haploid males and females (queens) at three developmental timepoints in Solenopsis invicta. In pupae, gene expression profiles of diploid males were very different from those of haploid males but nearly identical to those of queens. An unexpected shift in expression patterns emerged soon after adult eclosion, with diploid male patterns diverging from those of queens to resemble those of haploid males, a pattern retained in older adults. The finding that ploidy level effects on early gene expression override sex effects (including genes implicated in sperm production and pheromone production/perception) may explain diploid male sterility and lack of worker discrimination against them during development.

Analysis of gene expression patterns in animals

During my PhD, my aim was to provide new tools to increase our capacity to analyse gene expression patterns, and to study on a large-scale basis the evolution of gene expression in animals. Gene expression patterns (when and where a gene is expressed) are a key feature in understanding gene function, notably in development. It appears clear now that the evolution of developmental processes and of phenotypes is shaped both by evolution at the coding sequence level, and at the gene expression level.Studying gene expression evolution in animals, with complex expression patterns over tissues and developmental time, is still challenging. No tools are available to routinely compare expression patterns between different species, with precision, and on a large-scale basis. Studies on gene expression evolution are therefore performed only on small genes datasets, or using imprecise descriptions of expression patterns.The aim of my PhD was thus to develop and use novel bioinformatics resources, to study the evolution of gene expression. To this end, I developed the database Bgee (Base for Gene Expression Evolution). The approach of Bgee is to transform heterogeneous expression data (ESTs, microarrays, and in-situ hybridizations) into present/absent calls, and to annotate them to standard representations of anatomy and development of different species (anatomical ontologies). An extensive mapping between anatomies of species is then developed based on hypothesis of homology. These precise annotations to anatomies, and this extensive mapping between species, are the major assets of Bgee, and have required the involvement of many co-workers over the years. My main personal contribution is the development and the management of both the Bgee database and the web-application.Bgee is now on its ninth release, and includes an important gene expression dataset for 5 species (human, mouse, drosophila, zebrafish, Xenopus), with the most data from mouse, human and zebrafish. Using these three species, I have conducted an analysis of gene expression evolution after duplication in vertebrates.Gene duplication is thought to be a major source of novelty in evolution, and to participate to speciation. It has been suggested that the evolution of gene expression patterns might participate in the retention of duplicate genes. I performed a large-scale comparison of expression patterns of hundreds of duplicated genes to their singleton ortholog in an outgroup, including both small and large-scale duplicates, in three vertebrate species (human, mouse and zebrafish), and using highly accurate descriptions of expression patterns. My results showed unexpectedly high rates of de novo acquisition of expression domains after duplication (neofunctionalization), at least as high or higher than rates of partitioning of expression domains (subfunctionalization). I found differences in the evolution of expression of small- and large-scale duplicates, with small-scale duplicates more prone to neofunctionalization. Duplicates with neofunctionalization seemed to evolve under more relaxed selective pressure on the coding sequence. Finally, even with abundant and precise expression data, the majority fate I recovered was neither neo- nor subfunctionalization of expression domains, suggesting a major role for other mechanisms in duplicate gene retention.

Gene expression patterns unveil a new level of molecular heterogeneity in colorectal cancer.

The recognition that colorectal cancer (CRC) is a heterogeneous disease in terms of clinical behaviour and response to therapy translates into an urgent need for robust molecular disease subclassifiers that can explain this heterogeneity beyond current parameters (MSI, KRAS, BRAF). Attempts to fill this gap are emerging. The Cancer Genome Atlas (TGCA) reported two main CRC groups, based on the incidence and spectrum of mutated genes, and another paper reported an EMT expression signature defined subgroup. We performed a prior free analysis of CRC heterogeneity on 1113 CRC gene expression profiles and confronted our findings to established molecular determinants and clinical, histopathological and survival data. Unsupervised clustering based on gene modules allowed us to distinguish at least five different gene expression CRC subtypes, which we call surface crypt-like, lower crypt-like, CIMP-H-like, mesenchymal and mixed. A gene set enrichment analysis combined with literature search of gene module members identified distinct biological motifs in different subtypes. The subtypes, which were not derived based on outcome, nonetheless showed differences in prognosis. Known gene copy number variations and mutations in key cancer-associated genes differed between subtypes, but the subtypes provided molecular information beyond that contained in these variables. Morphological features significantly differed between subtypes. The objective existence of the subtypes and their clinical and molecular characteristics were validated in an independent set of 720 CRC expression profiles. Our subtypes provide a novel perspective on the heterogeneity of CRC. The proposed subtypes should be further explored retrospectively on existing clinical trial datasets and, when sufficiently robust, be prospectively assessed for clinical relevance in terms of prognosis and treatment response predictive capacity. Original microarray data were uploaded to the ArrayExpress database (http://www.ebi.ac.uk/arrayexpress/) under Accession Nos E-MTAB-990 and E-MTAB-1026. © 2013 Swiss Institute of Bioinformatics. Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Identification of molecular subtypes and gene expression patterns of breast cancer analysing RNA-seq data

Breast cancer is the most common diagnosed cancer and the leading cause of cancer death among females worldwide. It is considered a highly heterogeneous disease and it must be classified into more homogeneous groups. Hence, the purpose of this study was to classify breast tumors based on variations in gene expression patterns derived from RNA sequencing by using different class discovery methods. 42 breast tumors paired-samples were sequenced by Illumine Genome Analyzer and the data was analyzed and prepared by TopHat2 and htseq-count. As reported previously, breast cancer could be grouped into five main groups known as basal epithelial-like group, HER2 group, normal breast-like group and two Luminal groups with a distinctive expression profile. Classifying breast tumor samples by using PAM50 method, the most common subtype was Luminal B and was significantly associated with ESR1 and ERBB2 high expression. Luminal A subtype had ESR1 and SLC39A6 significant high expression, whereas HER2 subtype had a high expression of ERBB2 and CNNE1 genes and low luminal epithelial gene expression. Basal-like and normal-like subtypes were associated with low expression of ESR1, PgR and HER2, and had significant high expression of cytokeratins 5 and 17. Our results were similar compared with TGCA breast cancer data results and with known studies related with breast cancer classification. Classifying breast tumors could add significant prognostic and predictive information to standard parameters, and moreover, identify marker genes for each subtype to find a better therapy for patients with breast cancer.

Use of global gene expression patterns in mechanistic studies of oestrogen action in MCF7 human breast cancer cells

Over the years, the MCF7 human breast cancer cell line has provided a model system for the study of cellular and molecular mechanisms in oestrogen regulation of cell proliferation and in progression to oestrogen and antioestrogen independent growth. Global gene expression profiling has shown that oestrogen action in MCF7 cells involves the coordinated regulation of hundreds of genes across a wide range of functional groupings and that more genes are down regulated than upregulated. Adaptation to long-term oestrogen deprivation, which results in loss of oestrogen-responsive growth, involves alterations to gene patterns not only at early time points (0-4 weeks) but continuing through to later times (20-55 weeks), and even involves alterations to patterns of oestrogen-regulated gene expression. Only 48% of the genes which were regulated >= 2-fold by oestradiol in oestrogen-responsive cells retained this responsiveness after long-term oestrogen deprivation but other genes developed de novo oestrogen regulation. Long-term exposure to fulvestrant, which resulted in loss of growth inhibition by the antioestrogen, resulted in some very large fold changes in gene expression up to 10,000-fold. Comparison of gene profiles produced by environmental chemicals with oestrogenic properties showed that each ligand gave its own unique expression profile which suggests that environmental oestrogens entering the human breast may give rise to a more complex web of interference in cell function than simply mimicking oestrogen action at inappropriate times. (C) 2009 Elsevier Ltd. All rights reserved.

Relationships Among Murine CD11chigh Dendritic Cell Subsets as Revealed by Baseline Gene Expression Patterns

The functional relationships and properties of different subtypes of dendritic cells (DC) remain largely undefined. To better characterize these cells, we used global gene analysis to determine gene expression patterns among murine CD11c(high) DC subsets. CD4(+), CD8alpha(+), and CD8alpha(-) CD4(-) (double negative (DN)) DC were purified from spleens of normal C57/BL6 mice and analyzed using Affymetrix microarrays. The CD4(+) and CD8alpha(+) DC subsets showed distinct basal expression profiles differing by >200 individual genes. These included known DC subset markers as well as previously unrecognized, differentially expressed CD Ags such as CD1d, CD5, CD22, and CD72. Flow cytometric analysis confirmed differential expression in nine of nine cases, thereby validating the microarray analysis. Interestingly, the microarray expression profiles for DN cells strongly resembled those of CD4(+) DC, differing from them by <25 genes. This suggests that CD4(+) and DN DC are closely related phylogenetically, whereas CD8alpha(+) DC represent a more distant lineage, supporting the historical distinction between CD8alpha(+) and CD8alpha(-) DC. However, staining patterns revealed that in contrast to CD4(+) DC, the DN subset is heterogeneous and comprises at least two subpopulations. Gene Ontology and literature mining analyses of genes expressed differentially among DC subsets indicated strong associations with immune response parameters as well as cell differentiation and signaling. Such associations offer clues to possible unique functions of the CD11c(high) DC subsets that to date have been difficult to define as rigid distinctions.

Genomic DNA microarray comparison of gene expression patterns in Paracoccidioides brasiliensis mycelia and yeasts in vitro

Paracoccidioides brasiliensis is a thermally dimorphic fungus, and causes the most prevalent systemic mycosis in Latin America. Infection is initiated by inhalation of conidia or mycelial fragments by the host, followed by further differentiation into the yeast form. Information regarding gene expression by either form has rarely been addressed with respect to multiple time points of growth in culture. Here, we report on the construction of a genomic DNA microarray, covering approximately 25% of the genome of the organism, and its utilization in identifying genes and gene expression patterns during growth in vitro. Cloned, amplified inserts from randomly sheared genomic DNA (gDNA) and known control genes were printed onto glass slides to generate a microarray of over 12 000 elements. To examine gene expression, mRNA was extracted and amplified from mycelial or yeast cultures grown in semi-defined medium for 5, 8 and 14 days. Principal components analysis and hierarchical clustering indicated that yeast gene expression profiles differed greatly from those of mycelia, especially at earlier time points, and that mycelial gene expression changed less than gene expression in yeasts over time. Genes upregulated in yeasts were found to encode proteins shown to be involved in methionine/cysteine metabolism, respiratory and metabolic processes (of sugars, amino acids, proteins and lipids), transporters (small peptides, sugars, ions and toxins), regulatory proteins and transcription factors. Mycelial genes involved in processes such as cell division, protein catabolism, nucleotide biosynthesis and toxin and sugar transport showed differential expression. Sequenced clones were compared with Histoplasma capsulatum and Coccidioides posadasii genome sequences to assess potentially common pathways across species, such as sulfur and lipid metabolism, amino acid transporters, transcription factors and genes possibly related to virulence. We also analysed gene expression with time in culture and found that while transposable elements and components of respiratory pathways tended to increase in expression with time, genes encoding ribosomal structural proteins and protein catabolism tended to sharply decrease in expression over time, particularly in yeast. These findings expand our knowledge of the different morphological forms of P. brasiliensis during growth in culture.

Cryotolerance and global gene-expression patterns of Bos taurus indicus and Bos taurus taurus in vitro- and in vivo-produced blastocysts

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Pheochromocytoma in rats with multiple endocrine neoplasia (MENX) shares gene expression patterns with human pheochromocytoma

Pheochromocytomas are rare neoplasias of neural crest origin arising from chromaffin cells of the adrenal medulla and sympathetic ganglia (extra-adrenal pheochromocytoma). Pheochromocytoma that develop in rats homozygous for a loss-of-function mutation in p27Kip1 (MENX syndrome) show a clear progression from hyperplasia to tumor, offering the possibility to gain insight into tumor pathobiology. We compared the gene-expression signatures of both adrenomedullary hyperplasia and pheochromocytoma with normal rat adrenal medulla. Hyperplasia and tumor show very similar transcriptome profiles, indicating early determination of the tumorigenic signature. Overrepresentation of developmentally regulated neural genes was a feature of the rat lesions. Quantitative RT-PCR validated the up-regulation of 11 genes, including some involved in neural development: Cdkn2a, Cdkn2c, Neurod1, Gal, Bmp7, and Phox2a. Overexpression of these genes precedes histological changes in affected adrenal glands. Their presence at early stages of tumorigenesis indicates they are not acquired during progression and may be a result of the lack of functional p27Kip1. Adrenal and extra-adrenal pheochromocytoma development clearly follows diverged molecular pathways in MENX rats. To correlate these findings to human pheochromocytoma, we studied nine genes overexpressed in the rat lesions in 46 sporadic and familial human pheochromocytomas. The expression of GAL, DGKH, BMP7, PHOX2A, L1CAM, TCTE1, EBF3, SOX4, and HASH1 was up-regulated, although with different frequencies. Immunohistochemical staining detected high L1CAM expression selectively in 27 human pheochromocytomas but not in 140 nonchromaffin neuroendocrine tumors. These studies reveal clues to the molecular pathways involved in rat and human pheochromocytoma and identify previously unexplored biomarkers for clinical use.