Construction of a global map of human gene expression : the process, tools and analysis

This thesis studies human gene expression space using high throughput gene expression data from DNA microarrays. In molecular biology, high throughput techniques allow numerical measurements of expression of tens of thousands of genes simultaneously. In a single study, this data is traditionally obtained from a limited number of sample types with a small number of replicates. For organism-wide analysis, this data has been largely unavailable and the global structure of human transcriptome has remained unknown. This thesis introduces a human transcriptome map of different biological entities and analysis of its general structure. The map is constructed from gene expression data from the two largest public microarray data repositories, GEO and ArrayExpress. The creation of this map contributed to the development of ArrayExpress by identifying and retrofitting the previously unusable and missing data and by improving the access to its data. It also contributed to creation of several new tools for microarray data manipulation and establishment of data exchange between GEO and ArrayExpress. The data integration for the global map required creation of a new large ontology of human cell types, disease states, organism parts and cell lines. The ontology was used in a new text mining and decision tree based method for automatic conversion of human readable free text microarray data annotations into categorised format. The data comparability and minimisation of the systematic measurement errors that are characteristic to each lab- oratory in this large cross-laboratories integrated dataset, was ensured by computation of a range of microarray data quality metrics and exclusion of incomparable data. The structure of a global map of human gene expression was then explored by principal component analysis and hierarchical clustering using heuristics and help from another purpose built sample ontology. A preface and motivation to the construction and analysis of a global map of human gene expression is given by analysis of two microarray datasets of human malignant melanoma. The analysis of these sets incorporate indirect comparison of statistical methods for finding differentially expressed genes and point to the need to study gene expression on a global level.

Distress call-induced gene expression in the brain of the Indian short-nosed fruit bat, Cynopterus sphinx

Individuals in distress emit audible vocalizations to either warn or inform conspecifics. The Indian short-nosed fruit bat, Cynopterus sphinx, emits distress calls soon after becoming entangled in mist nets, which appear to attract conspecifics. Phase I of these distress calls is longer and louder, and includes a secondary peak, compared to phase II. Activity-dependent expression of egr-1 was examined in free-ranging C. sphinx following the emissions and responses to a distress call. We found that the level of expression of egr-1 was higher in bats that emitted a distress call, in adults that responded, and in pups than in silent bats. Up-regulated cDNA was amplified to identify the target gene (TOE1) of the protein Egr-1. The observed expression pattern Toe1 was similar to that of egr-1. These findings suggest that the neuronal activity related to recognition of a distress call and an auditory feedback mechanism induces the expression of Egr-1. Co-expression of egr-1 with Toe1 may play a role in initial triggering of the genetic mechanism that could be involved in the consolidation or stabilization of distress call memories.

Deregulation of miR-126 expression in colorectal cancer pathogenesis and its clinical significance

In this study, we investigated the expression profiles and clinicopathological significance of miR-126 in large cohort of patients with colorectal cancers as well the cellular repercussions of miR-126 in colon cancer cells along with its targets in-vitro. Down regulation of miR-126 expression was associated with histological subtypes, peri-neural tumour infiltration, microsatellite instability and pathological staging of colorectal cancers (p<0.05). Low miR-126 expression was also associated with poorer survival in patients with colorectal cancer. Analysis of matched tissues from the same patient revealed that approximately 70% of the tested patients had similar levels of expression of miR-126 in primary cancer and cancer metastases in both lymph node and distant metastases. In addition, induced overexpression of miR-126 showed reduced cell proliferation, increased apoptosis and decreased accumulation of cells in the G0-G1 phase of the colon cancer cells. Furthermore, SW480(+miR-126) cells showed reduced BCL-2 and increased P53 protein expression. To conclude, deregulation of miR-126 in colorectal cancer at the tissue and cellular levels as well as its correlation with various clinicopathological parameters confirm the cancer suppressive role of miR-126 in colorectal cancer.

Expression of the receptor guanylyl cyclase C and its ligands in reproductive tissues of the rat: A potential role for a novel signaling pathway in the epididymis

Guanylyl cyclase C (GC-C) is a membrane-associated form of guanylyl cyclase and serves as the receptor for the heat-stable enterotoxin (ST) peptide and endogenous ligands guanylin, uroguanylin, and lymphoguanylin. The major site of expression of GC-C is the intestinal epithelial cell, although GC-C is also expressed in extraintestinal tissue such as the kidney, airway epithelium, perinatal liver, stomach, brain, and adrenal glands. Binding of ligands to GC-C leads to accumulation of intracellular cGMP, the activation of protein kinases G and A, and phosphorylation of the cystic fibrosis transmembrane conductance regulator (CFTR), a chloride channel that regulates salt and water secretion. We examined the expression of GC-C and its ligands in various tissues of the reproductive tract of the rat. Using reverse transcriptase and the polymerase chain reaction, we demonstrated the presence of GC-C, uroguanylin, and guanylin mRNA in both male and female reproductive organs. Western blot analysis using a monoclonal antibody to GC-C revealed the presence of differentially glycosylated forms of GC-C in the caput and cauda epididymis. Exogenous addition of uroguanylin to minced epididymal tissue resulted in cGMP accumulation, suggesting an autocrine or endocrine activation of GC-C in this tissue. Immunohistochemical analyses demonstrated expression of GC-C in the tubular epithelial cells of both the caput epididymis and cauda epididymis. Our results suggest that the GC-C signaling pathway could converge on CFTR in the epididymis and perhaps control fluid and ion balance for optimal sperm maturation and storage in this tissue.

Toll-Like Receptors 7, 8, and 9 Expression and Function in Primary Human Cervical Cancer Langerhans Cells Evidence of Anergy

Objective: Human papillomavirus oncoproteins E6 and E7 down modulate Toll-like receptor (TLR) 9 expression in infected keratinocytes. We explored the status of expression and function of TLR7, TLR8, and TLR9 in primary human Langerhans cells (LCs) isolated from cervical tumors. Methodology: Single-cell suspensions were made from fresh tissues of squamous cell carcinoma (International Federation of Gynecology and Obstetrics stage IB2); myeloid dendritic cells were purified using CD1c magnetic activated cell separation kits. Langerhans cells were further flow sorted into CD1a(+)CD207(+) cells. Acute monocytic leukemia cell line THP-1-derived LCs (moLCs) formed the controls. mRNA from flow-sorted LCs was reverse transcribed to cDNA and TLR7, TLR8, and TLR9 amplified. Monocyte-derived Langerhans cells and cervical tumor LCs were stimulated with TLR7, TLR8, and TLR9 ligands. Culture supernatants were assayed for interleukin (IL) 1 beta, IL-6, IL-10, IL-12p70, interferon (IFN) alpha, interferon gamma, and tumor necrosis factor (TNF) alpha by Luminex multiplex bead array. Human papillomavirus was genotyped. Results: We have for the first time demonstrated that the acute monocytic leukemia cell line THP-1 can be differentiated into LCs in vitro. Although these moLCs. expressed all the 3 TLRs, tumor LCs expressed TLR7 and TLR8, but uniformly lacked TLR9. Also, moLCs secreted IL-6, IL-1 beta, and tumor necrosis factor alpha to TLR8 ligand and interferon alpha in response to TLR9 ligand; in contrast, tumor LCs did not express any cytokine to any of the 3 TLR ligands. Human papillomavirus type 16 was one of the common human papillomavirus types in all cases. Conclusions: Cervical tumor LCs lacked TLR9 expression and were functionally anergic to all the 3: TLR7, TLR8, and TLR9 ligands, which may play a crucial role in immune tolerance. The exact location of block(s) in TLR7 and TLR8 signaling needs to be investigated, which would have important immunotherapeutic implications.

An in vitro transposon system for highly regulated gene expression: construction of Escherichia coli strains with arabinose-dependent growth at low temperatures.

Placing a gene of interest under the control of an inducible promoter greatly aids the purification, localization and functional analysis of proteins but usually requires the sub-cloning of the gene of interest into an appropriate expression vector. Here, we describe an alternative approach employing in vitro transposition of Tn Omega P(BAD) to place the highly regulable, arabinose inducible P(BAD) promoter upstream of the gene to be expressed. The method is rapid, simple and facilitates the optimization of expression by producing constructs with variable distances between the P(BAD) promoter and the gene. To illustrate the use of this approach, we describe the construction of a strain of Escherichia coli in which growth at low temperatures on solid media is dependent on threshold levels of arabinose. Other uses of the transposable promoter are also discussed.

A ratiometric-based measure of gene co-expression

Current measures of global gene expression analyses, such as correlation and mutual information-based approaches, largely depend on the degree of association between mRNA levels and to a lesser extent on variability. I develop and implement a new approach, called the Ratiometric method, which is based on the coefficient of variation of the expression ratio of two genes, relying more on variation than previous methods. The advantage of such modus operandi is the ability to detect possible gene pair interactions regardless of the degree of expression dispersion across the sample group. Gene pairs with low expression dispersion, i.e., their absolute expressions remain constant across the sample group, are systematically missed by correlation and mutual information analyses. The superiority of the Ratiometric method in finding these gene pair interactions is demonstrated in a data set of RNA-seq B-cell samples from the 1000 Genomes Project Consortium. The Ratiometric method renders a more comprehensive recovery of KEGG pathways and GO-terms.

Cloning, characterization, and gene expression analysis of a novel cadmium metallothionein gene in Tetrahymena pigmentosa

A novel cadmium-inducible metallothionein (MT) gene (Tpig-MT1) was cloned and sequenced from the ciliate Tetrahymena pigmentosa. The number of deduced amino acids is 118. The polypeptide possesses CCC and CC clusters characteristic of typical Tetrahymena Cd-inducible MTs. The structure of Tpig-MT1 is different from the reported Cd-MT in T. pyriformis, T. thermophila and T. pigmentosa. Tpig-MT1 contains two intragenic tandem repeats with 72.9% identity described as Tpig-MT1 (repeat A1) and Tpig-MT1 (repeat A2). The transcriptional response of Tpig-MT1 gene to different heavy metals (Cd, Cu, Zn, Hg, Pb) and oxidative stress (H2O2) was measured using real-time quantitative PCR. The results showed that the gene was quickly induced (1 h) by the five heavy metals and the order of expression level was Hg>Pb>Cd>Cu>Zn. The induction effect of H2O2 was 5-fold after about 15 min, but soon decreased to a non-significant level (30 min). The genetic diversity of Tetrahymena MT genes is discussed in relation to the unique structure of the Tpig-MT1 gene and other reported Cd-MT isoforms. (C) 2008 Elsevier B.V. All rights reserved.

Differential expression of two Carassius auratus Mx genes in cultured CAB cells induced by grass carp hemorrhage virus and interferon

UV-inactivated GCHV (grass carp hemorrhage virus) is able to induce an antiviral state in cultured CAB cells (crucian carp Carassius auratus blastulae embryonic cells) via the production of interferon (IFN). In the current work, the full-length cDNAs of two Mx genes, termed CaMx1 and CaMx2, have been cloned and sequenced from UV-inactivated GCHV-infected and still IFN-producing CAB cells by suppression subtractive hybridization. Their putative proteins show the characteristically structural features of mammalian IFN-induced Mx proteins, including GTP-binding motif, dynamin family signature and leucine zipper motif. CaMx1 exhibits 85% sequence identity to zebrafish MxA and 72-74% to three Atlantic salmon Mx proteins. CaMx2 is most similar to zebrafish MxE, with 80% identity, and then rainbow trout Mx3, with 52%. Constitutive expression was detected by RT-PCR for CaMx1, but not for CaMx2, in normal CAB cells, but their up-regulations could be induced after treatment with active GCHV, UV-inactivated GCHV and CAB IFN. Distinct kinetics of expression was observed for either CaMx1 or CaMx2 corresponding to the three stimuli, and even between CaMx1 and CaMx2, corresponding to the same stimulus. Upon virus infection, the transcriptional induction was strongly blocked for CaMx2 by cycloheximide (CHX), whereas almost nothing was observed for CaMx1. By contrast, following treatment with CAB IFN, CHX did not inhibit either gene transcription. Collectively, these results suggest that there are very distinct mechanisms for modulating the expression of both CaMx1 and CaMx2 in normal and GCHV-infected CAB cells.

Identification and expression analysis of two IFN-inducible genes in crucian carp (Carassius auratus L.)

Interferon (IFN) exerts its antiviral effects mainly through activation of a subset of IFN-stimulated genes (ISG), but relatively few of fish ISGs have been isolated and characterized so far. Here, we report two fish ISGs, termed CaIF158 and CaIF156, cloned from a subtractive cDNA library constructed with mRNAs obtained from crucian carp (Carassius auratus L.) blastulae embryonic (CAB) cells infected by UV-inactivated GCHV and mock-infected cells. Database search revealed that both ISGs had a high-level homology with all members of a well conserved gene family with multiple tetratricopeptide repeat (TPR) motifs, including human IF160, IF158, IF156, IFI54 and their homologues in some other mammalian species. The transcripts of CaIF158 and CaIF156 were undetectable in CAB cells but could be induced by active GCHV, UV-inactivated GCHV or CAB IFN. Analysis of expression difference between them and IFN signal factors, CaSTAT1 and CaIRF7, indicated that their transcriptions were mediated possibly through JAK-STAT signal pathway, which was further supported by the induction analysis in UV-inactivated GCHV infected, IFN-treated and untreated cells in the presence or absence of cycloheximide (CHX), a potent inhibitor of protein synthesis. In addition, a pufferfish (Fugu rubrides) DNA sequence representing putative FrIFI56 was also revealed when CalF158 and CalF156 were used to search the pufferfish genome database. Phylogenetic analysis showed that these fish ISGs form a unique clad independent of mammalian homologues, reflecting a distant evolutionary relationship from mammals. These studies identified the first teleost IFI56 and IFI58 orthologues. (C) 2003 Elsevier B.V All rights reserved.

Expression pattern of dmrt4 from olive flounder (Paralichthys olivaceus) in adult gonads and during embryogenesis

The dmrt (doublesex and mab-3 related transcription factor) gene family comprises several transcription factors that share a conserved DM domain. Dmrt1 is considered to be involved in sexual development, but the precise function of other family members is unclear. In this study, we isolated genomic DNA and cDNA sequences of dmrt4, a member of the dmrt gene family, from olive flounder, Paralichthys olivaceus, through genome walking and real-time reverse transcriptase (RT)-PCR. Sequence analysis indicated that its genomic DNA contains two exons and one intron. A transcriptional factor binding sites prediction program identified a sexual development-related protein, Sox9 (Sry-like HMG box containing 9) in its 5' promoter. Protein alignment and phylogenetic analysis suggested that flounder Dmrt4 is closely related to tilapia Dmo (DM domain gene in ovary). The expression of dmrt4 in adult flounder was sexually dimorphic, as shown by real-time RT-PCR analysis, with strong expression in the testis but very weak expression in the ovary. Its expression was also strong in the brain and gill, but there was only weak or no expression at all in some of the other tissues tested of both sexes. During embryogenesis, its expression was detected in most developmental stages, although the level of expression was distinctive of the various stages. Whole mount in situ hybridization revealed that the dmrt4 was expressed in the otic placodes, forebrain, telencephalon and olfactory placodes of embryos at different developmental stages. These results will improve our understanding of the possible role of flounder dmrt4 in the development of the gonads, nervous system and sense organs.

Expression profiles of penaeidin from Fenneropenaeus chinensis in response to WSSV and vibrio infection by real-time PCR

Penaeidin from Chinese shrimp (Fenneropenaeus chinensis) has proved to be one of the most important antimicrobial peptides in the bodies of animals. The relative quantitative real-time PCR method is developed to study through time, the mRNA expression profile of penaeidin in the muscle and haemocyte tissue of Chinese shrimp infected with vibrio (Vibrio anguillarum) and WSSV (white spot syndrome virus). Research results showed that the same pathogens infection experiments produced similar gene expression profile in different tissues while different expression profiles appeared in the same tissues infected by different exterior pathogens. In vibrio infection experiments, a "U" Re expression profile resulted. Expression levels of penaeidin increased and surpassed the non-stimulated level, indicating that penaeidin from Chinese shrimp has noticeable antimicrobial activities. In WSSV infection experiments, the expression profile appeared as an inverse "U" with the expression of penaeidin gradually decreasing to below baseline level after 24 h. The expression of antimicrobial peptides gene in mRNA level in response to virus infection in shrimp showed that international mechanisms of virus to haemocytes and microbial to haemocytes are completely different. Decline of penaeidins expression levels may be due to haemocytes being destroyed by WSSV or that the virus can inhibit the expression of penaeidins by yet undiscovered modes. The expression profiles of penaeidin in response to exterior pathogen and the difference of expression profiles between vibrio and WSSV infection provided some clues to further understanding the complex innate immune mechanism in shrimp.

Identification and expression of the amphioxus Cks1 gene

A full-length Cks1 homologue gene, AmphiCks1, was identified in amphioxus, Branchiostoma belcheri tsingtauense. Sequence characteristics, phylogeny and patterns of expression during embryonic and larval development were established. The protein predicted from AmphiCks1 showed high sequence identity with vertebrate and invertebrate homologues. Protein structural studies and phylogenetic analysis suggested that Cks homologues are evolutionarily conserved. The AmphiCks1 transcript was detected in most early developmental stages by northern blotting and whole-mount in situ hybridization, suggesting a role for the gene in cell division. (c) 2005 International Federation for Cell Biology. Published by Elsevier Ltd. All rights reserved.

A Dicer-1 gene from white shrimp Litopenaeus vannamei: Expression pattern in the processes of immune response and larval development

Dicer is a member of the RNAase III family which catalyzes the cleavage of double-stranded RNA to small interfering RNAs and micro RNAs, and then directs sequence-specific gene silencing. In this paper, the full-length cDNA of Dicer-1 was cloned from white shrimp Litopenaeus vannamei (designated as LvDcr1). It was of 7636 bp, including a poly A tail, a 5' UTR of 136 bp, a 3' UTR of 78 bp, and an open reading frame (ORF) of 7422 bp encoding a putative protein of 2473 amino acids. The predicted amino acid sequence comprised all recognized functional domains found in other Dicer-1 homologues and showed the highest (97.7%) similarity to the Dicer-1 from tiger shrimp Penaeus mondon. Quantitative real-time PCR was employed to investigate the tissue distribution of LvDcr1 mRNA, and its expression in shrimps under virus challenge and larvae at different developmental stages. The LvDcr1 mRNA could be detected in all examined tissues with the highest expression level in hemocyte, and was up-regulated in hemocytes and gills after virus injection. These results indicated that LvDcr1 was involved in antiviral defense in adult shrimp. During the developmental stages from fertilized egg to postlarva VII, LvDcr1 was constitutively expressed at all examined development stages, but the expression level varied significantly. The highest expression level was observed in fertilized eggs and followed a decrease from fertilized egg to nauplius I stage. Then, the higher levels of expression were detected at nauplius V and postlarva stages. LvDcr1 expression regularly increased at the upper phase of nauplius, zoea and mysis stages than their prophase. The different expression of LvDcr1 in the larval stages could provide clues for understanding the early innate immunity in the process of shrimp larval development. (C) 2010 Elsevier Ltd. All rights reserved.

Molecular structure and expression patterns of flounder (Paralichthys olivaceus) Myf-5, a myogenic regulatory factor

Myf-5, a member of the myogenic regulatory factors (MRF), has been shown to be expressed in muscle precursors in early stage zebrafish embryos. The MRFs, including MyoD, Myf-5, Myogenin and MR-F4, belong to the basic Helix-Loop-Helix transcription factors that contain a conserved basic Helix-Loop-Helix (bHLH) domain. To better understand the role of Myf-5 in the development of fish muscles, we have isolated the Myf-5 genomic sequence and cDNA from Flounder (Paralichthys olivaceus), and analyzed its structures and patterns of expression. Promoter analysis identified several putative transcription factor binding sites such as an E-box, NF-Y sites that might confer muscle-specific expression. Myf-5 transcripts were first detected in the paraxial mesoderm that gives rise to slow muscles. During somitogenesis, Myf-5 expression was found in developing somites. Myf-5 expression decreased gradually in somites in the anterior region, but remained strong in the newly formed somites. In the hatching stage, the expression was also detected in other muscle cells such as head muscle and fin muscle. In the growing fish, RT-PCR results showed that Myf-5 was expressed in the skeletal muscle and intestine. (c) 2006 Elsevier Inc. All rights reserved.