Multivariate gene expression analysis reveals functional connectivity changes between normal/tumoral prostates

Abstract Background Prostate cancer is a leading cause of death in the male population, therefore, a comprehensive study about the genes and the molecular networks involved in the tumoral prostate process becomes necessary. In order to understand the biological process behind potential biomarkers, we have analyzed a set of 57 cDNA microarrays containing ~25,000 genes. Results Principal Component Analysis (PCA) combined with the Maximum-entropy Linear Discriminant Analysis (MLDA) were applied in order to identify genes with the most discriminative information between normal and tumoral prostatic tissues. Data analysis was carried out using three different approaches, namely: (i) differences in gene expression levels between normal and tumoral conditions from an univariate point of view; (ii) in a multivariate fashion using MLDA; and (iii) with a dependence network approach. Our results show that malignant transformation in the prostatic tissue is more related to functional connectivity changes in their dependence networks than to differential gene expression. The MYLK, KLK2, KLK3, HAN11, LTF, CSRP1 and TGM4 genes presented significant changes in their functional connectivity between normal and tumoral conditions and were also classified as the top seven most informative genes for the prostate cancer genesis process by our discriminant analysis. Moreover, among the identified genes we found classically known biomarkers and genes which are closely related to tumoral prostate, such as KLK3 and KLK2 and several other potential ones. Conclusion We have demonstrated that changes in functional connectivity may be implicit in the biological process which renders some genes more informative to discriminate between normal and tumoral conditions. Using the proposed method, namely, MLDA, in order to analyze the multivariate characteristic of genes, it was possible to capture the changes in dependence networks which are related to cell transformation.

Semantic integration of gene expression analysis tools and data sources using software connectors

Abstract Background The study and analysis of gene expression measurements is the primary focus of functional genomics. Once expression data is available, biologists are faced with the task of extracting (new) knowledge associated to the underlying biological phenomenon. Most often, in order to perform this task, biologists execute a number of analysis activities on the available gene expression dataset rather than a single analysis activity. The integration of heteregeneous tools and data sources to create an integrated analysis environment represents a challenging and error-prone task. Semantic integration enables the assignment of unambiguous meanings to data shared among different applications in an integrated environment, allowing the exchange of data in a semantically consistent and meaningful way. This work aims at developing an ontology-based methodology for the semantic integration of gene expression analysis tools and data sources. The proposed methodology relies on software connectors to support not only the access to heterogeneous data sources but also the definition of transformation rules on exchanged data. Results We have studied the different challenges involved in the integration of computer systems and the role software connectors play in this task. We have also studied a number of gene expression technologies, analysis tools and related ontologies in order to devise basic integration scenarios and propose a reference ontology for the gene expression domain. Then, we have defined a number of activities and associated guidelines to prescribe how the development of connectors should be carried out. Finally, we have applied the proposed methodology in the construction of three different integration scenarios involving the use of different tools for the analysis of different types of gene expression data. Conclusions The proposed methodology facilitates the development of connectors capable of semantically integrating different gene expression analysis tools and data sources. The methodology can be used in the development of connectors supporting both simple and nontrivial processing requirements, thus assuring accurate data exchange and information interpretation from exchanged data.

Gene expression analysis in ‘Candidatus Phytoplasma mali’-resistant and -susceptible Malus genotypes

Apple proliferation (AP) disease is the most important graft-transmissible and vector-borne disease of apple in Europe. ‘Candidatus Phytoplasma mali’ (Ca. P. mali) is the causal agent of AP. Apple (Malus x domestica) and other Malus species are the only known woody hosts. In European apple orchards, the cultivars are mainly grafted on one rootstock, M. x domestica cv. M9. M9 like all other M. x domestica cultivars is susceptible to ‘Ca. P. mali’. Resistance to AP was found in the wild genotype Malus sieboldii (MS) and in MS-derived hybrids but they were characterised by poor agronomic value. The breeding of a new rootstock carrying the resistant and the agronomic traits was the major aim of a project of which this work is a part. The objective was to shed light into the unknown resistance mechanism. The plant-phytoplasma interaction was studied by analysing differences between the ‘Ca. P. mali’-resistant and -susceptible genotypes related to constitutively expressed genes or to induced genes during infection. The cDNA-Amplified Fragment Length Polymorphism (cDNA-AFLP) technique was employed in both approaches. Differences related to constitutively expressed genes were identified between two ‘Ca. P. mali’-resistant hybrid genotypes (4551 and H0909) and the ‘Ca. P. mali’-susceptible M9. 232 cDNA-AFLP bands present in the two resistant genotypes but absent in the susceptible one were isolated but several different products associated to each band were found. Therefore, two different macroarray hybridisation experiments were performed with the cDNA-AFLP fragments yielding 40 sequences encoding for genes of unknown function or a wide array of functions including plant defence. In the second approach, individuation and analysis of the induced genes was carried out exploiting an in vitro system in which healthy and ‘Ca. P. mali’-infected micropropagated plants were maintained under controlled conditions. Infection trials using in vitro grafting of ‘Ca. P. mali’ showed that the resistance phenotype could be reproduced in this system. In addition, ex vitro plants were generated as an independent control of the genes differentially expressed in the in vitro plants. The cDNA-AFLP analysis in in vitro plants yielded 63 bands characterised by over-expression in the infected state of both the H0909 and MS genotypes. The major part (37 %) of the associated sequences showed homology with products of unknown function. The other genes were involved in plant defence, energy transport/oxidative stress response, protein metabolism and cellular growth. Real-time qPCR analysis was employed to validate the differential expression of the genes individuated in the cDNA-AFLP analysis. Since no internal controls were available for the study of the gene expression in Malus, an analysis on housekeeping genes was performed. The most stably expressed genes were the elongation factor-1 α (EF1) and the eukaryotic translation initiation factor 4-A (eIF4A). Twelve out of 20 genes investigated through qPCR were significantly differentially expressed in at least one genotype either in in vitro plants or in ex vitro plants. Overall, about 20% of the genes confirmed their cDNA-AFLP expression pattern in M. sieboldii or H0909. On the contrary, 30 % of the genes showed down-regulation or were not differentially expressed. For the remaining 50 % of the genes a contrasting behaviour was observed. The qPCR data could be interpreted as follows: the phytoplasma infection unbalance photosynthetic activity and photorespiration down-regulating genes involved in photosynthesis and in the electron transfer chain. As result, and in contrast to M. x domestica genotypes, an up-regulation of genes of the general response against pathogens was found in MS. These genes involved the pathway of H2O2 and the production of secondary metabolites leading to the hypothesis that a response based on the accumulation of H2O2 in MS would be at the base of its resistance. This resembles a phenomenon known as “recovery” where the spontaneous remission of the symptoms is observed in old susceptible plants but occurring in a stochastic way while the resistance in MS is an inducible but stable feature. As additional product of this work three cDNA-AFLP-derived markers were developed which showed independent distribution among the seedlings of two breeding progenies and were associated to a genomic region characteristic of MS. These markers will contribute to the development of molecular markers for the resistance as well as to map the resistance on the Malus genome.

Genus Bifidobacterium: taxonomy studies and gene expression analysis on folate pathway

Folates (vitamin B9) are essential water soluble vitamins, whose deficiency in humans may contribute to the onset of several diseases, such as anaemia, cancer, cardiovascular diseases, neurological problems as well as defects in embryonic development. Human and other mammals are unable to synthesize ex novo folate obtaining it from exogenous sources, via intestinal absorption. Recently the gut microbiota has been identified as an important source of folates and the selection and use of folate producing microorganisms represents an innovative strategy to increase human folate levels. The aim of this thesis was to gain a fundamental understanding of folate metabolism in Bifidobacterium adolescentis. The work was subdivided in three main phases, also aimed to solve different problems encountered working with Bifidobacterium strains. First, a new identification method (based on PCR-RFLP of hsp60 gene) was specifically developed to identify Bifidobacterium strains. Secondly, Bifidobacterium adolescentis biodiversity was explored in order to recognize representing strains of this species to be screened for their folate production ability. Results showed that this species is characterized by a wide variability and support the idea that a possible new taxonomic re-organization would be required. Finally B. adolescentis folate metabolism was studied using a double approach. A quantitative analysis of folate content was complemented by the examination of expression levels of genes involved in folate related pathways. For the normalization process, required to increase the robustness of the qRT-PCR analysis, an appropriate set of reference genes was tested using two different algorithms. Results demonstrate that B.adolescentis strains may represent an endogenous source of natural folate and they could be used to fortify fermented dairy products. This bio-fortification strategy presents many advantages for the consumer, providing native folate forms more bio-available, and not implicated in the discussed controversy concerning the safety of high intake of synthetic folic acid.

Primaquine-induced differential gene expression analysis in mice liver using DNA microarrays

Primaquine (PQ). a clinically important derivative of 8-aminoquinoline used against the hepatic stages (hypnozoites) of Plasmodium vivax and Plasmodium ova Ie. was studied to evaluate and compare between mRNA expression. and biochemical and histological parameters of hepatic stress in adult Swiss mice (Mus musculus). Following single oral dose of PQ (40 mglkg. bw). alanine aminotransferase (ALT) and aspartate aminotransferase (AST) along with hematoxylin and eosin stained liver sections did not show any signs of hepatic stress at 6. 12 and 24 h except for ALT activity at 6 h. However. analysis at RNA transcript level revealed consistent and significant deregulation (p<0.01 and twofold) of 16 probes corresponding to important cellular processes such as protein transportation. transcription regulation. intracellular signaling. protein synthesis, hematopoiesis, cell adhesion and cell proliferation. Pathway analysis identified large number of affected genes corresponding to 40 Gene Ontology terms having a z score greaibr than 2. These results indicate that PQ at high doses may affect gene expression in liver and may produce undesirable outcomes if consumed for longer durations.

Early gene expression analysis in 9L orthotopic tumor-bearing rats identifies immune modulation in molecular response to synchrotron microbeam radiation therapy

Synchrotron Microbeam Radiation Therapy (MRT) relies on the spatial fractionation of the synchrotron photon beam into parallel micro-beams applying several hundred of grays in their paths. Several works have reported the therapeutic interest of the radiotherapy modality at preclinical level, but biological mechanisms responsible for the described efficacy are not fully understood to date. The aim of this study was to identify the early transcriptomic responses of normal brain and glioma tissue in rats after MRT irradiation (400Gy). The transcriptomic analysis of similarly irradiated normal brain and tumor tissues was performed 6 hours after irradiation of 9 L orthotopically tumor-bearing rats. Pangenomic analysis revealed 1012 overexpressed and 497 repressed genes in the irradiated contralateral normal tissue and 344 induced and 210 repressed genes in tumor tissue. These genes were grouped in a total of 135 canonical pathways. More than half were common to both tissues with a predominance for immunity or inflammation (64 and 67% of genes for normal and tumor tissues, respectively). Several pathways involving HMGB1, toll-like receptors, C-type lectins and CD36 may serve as a link between biochemical changes triggered by irradiation and inflammation and immunological challenge. Most immune cell populations were involved: macrophages, dendritic cells, natural killer, T and B lymphocytes. Among them, our results highlighted the involvement of Th17 cell population, recently described in tumor. The immune response was regulated by a large network of mediators comprising growth factors, cytokines, lymphokines. In conclusion, early response to MRT is mainly based on inflammation and immunity which appear therefore as major contributors to MRT efficacy.

GLUCAGON GENE EXPRESSION: ANALYSIS OF PREPROGLUCAGON

Glucagon is a 29 amino acid polypeptide hormone produced in the (alpha) cells of the pancreatic islets. The purpose of this research was to understand better the role of glucagon in the regulation of metabolic processes. As with other polypeptide hormones, the synthesis of glucagon is thought to involve a larger precursor, which is then enzymatically cleaved to the functional form. The specific research objectives were to obtain cloned copies of the messenger RNA (mRNA) for pancreatic glucagon, to determine their primary sequences, and from this coding information to deduce the amino acid sequence of the initial glucagon precursor. From this suggested preproglucagon sequence and prior information on possible proglucagon intermediate processing products, the overall objective of this research is to propose a possible pathway for the biosynthesis of pancreatic glucagon.^ Synthetic oligodeoxynucleotide probes of 14-nucleotides (14-mer) and 17-nucleotides (a 17-mer) complementary to codons specifying a unique sequence of mature glucagon were synthesized. The ('32)P-labeled-14-mer was hybridized with size-fractionated fetal bovine pancreatic poly(A('+))RNA bound to nitrocellulose. RNA fractions of (TURN)14S were found to hybridize specifically, resulting in an (TURN)10-fold enrichment for these sequences. These poly(A('+))RNAs were translated in a cell-free system and the products analyzed by gel electrophoresis. The translation products were found to be enriched for a protein of the putative size of mammalian preproglucagon ((TURN)21 kd). These enriched RNA fractions were used to construct a complementary DNA (cDNA) library is plasmid pBR322.^ Screening of duplicate colony filters with the ('32)P-labeled-17-mer and a ('32)P-labeled-17-mer-primed cDNA probe indicated 25 possible glucagon clones from 3100 colonies screened. Restriction mapping of 6 of these clones suggested that they represented a single mRNA species. Primary sequence analysis of one clone containing a 1200 base pair DNA insert revealed that it contained essentially a full-length copy of glucagon cDNA.^ Analaysis of the cDNA suggested that it encoded an initial translation product of 180 amino acids with an M(,r) = 21 kd. The first initiation codon (ATG, methionine) followed by the longest open reading frame of 540 nucleotides was preceded by a 5'-untranslated region of 90 nucleotides, and was followed by a longer 3'-untranslated region of 471 nucleotides, resulting in a total of 1101 nucleotides. . . . (Author's abstract exceeds stipulated maximum length. Discontinued here with permission of author.) UMI ^

Yeast microarrays for genome wide parallel genetic and gene expression analysis

We have developed high-density DNA microarrays of yeast ORFs. These microarrays can monitor hybridization to ORFs for applications such as quantitative differential gene expression analysis and screening for sequence polymorphisms. Automated scripts retrieved sequence information from public databases to locate predicted ORFs and select appropriate primers for amplification. The primers were used to amplify yeast ORFs in 96-well plates, and the resulting products were arrayed using an automated micro arraying device. Arrays containing up to 2,479 yeast ORFs were printed on a single slide. The hybridization of fluorescently labeled samples to the array were detected and quantitated with a laser confocal scanning microscope. Applications of the microarrays are shown for genetic and gene expression analysis at the whole genome level.

Using the R package Shiny to create web applications that facilitate quantitative gene expression analysis

Thesis (Master's)--University of Washington, 2016-06

Comparative gene expression analysis of NKT cell subpopulations

Natural killer T (NKT) cells are a lymphocyte lineage, which has diverse immune regulatory activities in many disease settings. Most previous studies have investigated the functions of this family of cells as a single entity, but more recent evidence highlights the distinct functional and phenotypic properties of NKT cell subpopulations. It is likely that the diverse functions of NKT cells are regulated and coordinated by these different NKT subsets. Little is known about how NKT subsets differ in their interactions with the host. We have undertaken the first microarray analysis comparing the gene expression profiles of activated human NKT cell subpopulations, including CD8(+) NKT cells, which have often been overlooked. We describe the significant gene expression differences among NKT cell subpopulations and some of the molecules likely to confer their distinct functional roles. Several genes not associated previously with NKT cells were shown to be expressed differentially in specific NKT cell subpopulations. Our findings provide new insights into the NKT cell family, which may direct further research toward better manipulation of NKT cells for therapeutic applications.

Genetic variation of temperature-regulated curd induction in cauliflower: elucidation of floral transition by genome-wide association mapping and gene expression analysis

Cauliflower (Brassica oleracea var. botrytis) is a vernalization-responsive crop. High ambient temperatures delay harvest time. The elucidation of the genetic regulation of floral transition is highly interesting for a precise harvest scheduling and to ensure stable market supply. This study aims at genetic dissection of temperature-dependent curd induction in cauliflower by genome-wide association studies and gene expression analysis. To assess temperature dependent curd induction, two greenhouse trials under distinct temperature regimes were conducted on a diversity panel consisting of 111 cauliflower commercial parent lines, genotyped with 14,385 SNPs. Broad phenotypic variation and high heritability (0.93) were observed for temperature-related curd induction within the cauliflower population. GWA mapping identified a total of 18 QTL localized on chromosomes O1, O2, O3, O4, O6, O8, and O9 for curding time under two distinct temperature regimes. Among those, several QTL are localized within regions of promising candidate flowering genes. Inferring population structure and genetic relatedness among the diversity set assigned three main genetic clusters. Linkage disequilibrium (LD) patterns estimated global LD extent of r(2) = 0.06 and a maximum physical distance of 400 kb for genetic linkage. Transcriptional profiling of flowering genes FLOWERING LOCUS C (BoFLC) and VERNALIZATION 2 (BoVRN2) was performed, showing increased expression levels of BoVRN2 in genotypes with faster curding. However, functional relevance of BoVRN2 and BoFLC2 could not consistently be supported, which probably suggests to act facultative and/or might evidence for BoVRN2/BoFLC-independent mechanisms in temperature regulated floral transition in cauliflower. Genetic insights in temperature-regulated curd induction can underpin genetically informed phenology models and benefit molecular breeding strategies toward the development of thermo-tolerant cultivars.

A comparison of collection techniques for gene expression analysis of human oral taste tissue

Variability in human taste perception is associated with both genetic and environmental factors. The influence of taste receptor expression on this variability is unknown, in part, due to the difficulty in obtaining human oral tissue that enables quantitative expression measures of taste genes. In a comparison of six current techniques (Oragene RNeasy Kit, Isohelix swab, Livibrush cytobrush, tongue saliva, cheek saliva collection, and fungiform papillae biopsy), we identify the fungiform papillae biopsy is the optimal sampling technique to analyse human taste gene expression. The fungiform papillae biopsy resulted in the highest RNA integrity, enabling amplification of all the assessed taste receptor genes (TAS1R1, TAS1R2, TAS1R3, SCNN1A and CD36) and taste tissue marker genes (NCAM1, GNAT3 and PLCβ2). Furthermore, quantitative expression was observed in a subset of taste genes assessed from the saliva collection techniques (cheek saliva, tongue saliva and Oragene RNA kit). These saliva collection techniques may be useful as a non-invasive alternative sampling technique to the fungiform papillae biopsy. Identification of the fungiform papillae biopsy as the optimal collection method will facilitate further research into understanding the effect of gene expression on variability in human taste perception.

Analysis of differential gene expression based on Bayesian estimation of variance

Gene expression is arguably the most important indicator of biological function. Thus identifying differentially expressed genes is one of the main aims of high throughout studies that use microarray and RNAseq platforms to study deregulated cellular pathways. There are many tools for analysing differentia gene expression from transciptomic datasets. The major challenge of this topic is to estimate gene expression variance due to the high amount of ‘background noise’ that is generated from biological equipment and the lack of biological replicates. Bayesian inference has been widely used in the bioinformatics field. In this work, we reveal that the prior knowledge employed in the Bayesian framework also helps to improve the accuracy of differential gene expression analysis when using a small number of replicates. We have developed a differential analysis tool that uses Bayesian estimation of the variance of gene expression for use with small numbers of biological replicates. Our method is more consistent when compared to the widely used cyber-t tool that successfully introduced the Bayesian framework to differential analysis. We also provide a user-friendly web based Graphic User Interface for biologists to use with microarray and RNAseq data. Bayesian inference can compensate for the instability of variance caused when using a small number of biological replicates by using pseudo replicates as prior knowledge. We also show that our new strategy to select pseudo replicates will improve the performance of the analysis. - See more at: http://www.eurekaselect.com/node/138761/article#sthash.VeK9xl5k.dpuf

GEDI: a user-friendly toolbox for analysis of large-scale gene expression data

Abstract Background Several mathematical and statistical methods have been proposed in the last few years to analyze microarray data. Most of those methods involve complicated formulas, and software implementations that require advanced computer programming skills. Researchers from other areas may experience difficulties when they attempting to use those methods in their research. Here we present an user-friendly toolbox which allows large-scale gene expression analysis to be carried out by biomedical researchers with limited programming skills. Results Here, we introduce an user-friendly toolbox called GEDI (Gene Expression Data Interpreter), an extensible, open-source, and freely-available tool that we believe will be useful to a wide range of laboratories, and to researchers with no background in Mathematics and Computer Science, allowing them to analyze their own data by applying both classical and advanced approaches developed and recently published by Fujita et al. Conclusion GEDI is an integrated user-friendly viewer that combines the state of the art SVR, DVAR and SVAR algorithms, previously developed by us. It facilitates the application of SVR, DVAR and SVAR, further than the mathematical formulas present in the corresponding publications, and allows one to better understand the results by means of available visualizations. Both running the statistical methods and visualizing the results are carried out within the graphical user interface, rendering these algorithms accessible to the broad community of researchers in Molecular Biology.