Advantages of mRNA amplification for microarray analysis

LINCATRA: Two-cycle method to amplify RNA for transcriptome analysis from formalin-fixed paraffin-embedded tissue

Whole transcriptome analysis (WTA) using RNA extracted from Formalin Fixed Paraffin Embedded (FFPE) tissue is an invaluable tool to understand the molecular pathology of disease. RNA extracted from FFPE tissue is either degraded and/or in very low quantities hampering gene expression analysis. Earlier studies described protocols applied for cellular RNA using poly-A primer-based linear amplification. The current study describes a method, LINCATRA (LINear amplifiCAtion of RNA for whole TRAnscriptome analysis). It employs random nonamer primer based method which can amplify short, fragmented RNA with high fidelity from as low as 5 ng to obtain enough material for WTA. The two-cycle method significantly amplified RNA at ∼1000 folds (p < 0.0001) improving the mean read lengths (p < 0.05) in WTA. Overall, increased mean read length positively correlated with on-target reads (Pearson's r = 0.71, p < 0.0001) in both amplified and unamplified RNA-seq analysis. Gene expression analysis compared between unamplified and amplified group displayed substantial overlap of the differentially expressed genes (DEGs) (log2 fold change cut-off < -2 and >2, p < 0.05) identified between lung cancer and asthma cohorts validating the method developed. This method is applicable in clinical molecular pathology field for both diagnostics and elucidation of disease mechanisms.

Functional genomics in sand fly-derived Leishmania promastigotes

Background

Leishmania development in the sand fly gut leads to highly infective forms called metacyclic promastigotes. This process can be routinely mimicked in culture. Gene expression–profiling studies by transcriptome analysis have been performed with the aim of studying promastigote forms in the sand fly gut, as well as differences between sand fly–and culture-derived promastigotes.

Findings

Transcriptome analysis has revealed the crucial role of the microenvironment in parasite development within the sand fly gut because substantial differences and moderate correlation between the transcriptomes of cultured and sand fly–derived promastigotes have been found. Sand fly–derived metacyclics are more infective than metacyclics in culture. Therefore, some caution should be exercised when using cultured promastigotes, depending on the experimental design. The most remarkable examples are the hydrophilic acidic surface protein/small endoplasmic reticulum protein (HASP/SHERP) cluster, the glycoprotein 63 (gp63), and autophagy genes, which are up-regulated in sand fly–derived promastigotes compared with cultured promastigotes. Because HASP/SHERP genes are up-regulated in nectomonad and metacyclic promastigotes in the sand fly, the encoded proteins are not metacyclic specific. Metacyclic promastigotes are distinguished by morphology and high infectivity. Isolating them from the sand fly gut is not exempt from technical difficulty, because other promastigote forms remain in the gut even 15 days after infection. Leishmania major procyclic promastigotes within the sand fly gut up-regulate genes involved in cell cycle regulation and glucose catabolism, whereas metacyclics increase transcript levels of fatty acid biosynthesis and ATP-coupled proton transport genes. Most parasite's signal transduction pathways remain uncharacterized. Future elucidation may improve understanding of parasite development, particularly signaling molecule-encoding genes in sand fly versus culture and between promastigote forms in the sand fly gut.

Conclusions

Transcriptome analysis has been demonstrated to be technically efficacious to study differential gene expression in sand fly gut promastigote forms. Transcript and protein levels are not well correlated in these organisms (approximately 25% quantitative coincidences), especially under stress situations and at differentiation processes. However, transcript and protein levels behave similarly in approximately 60% of cases from a qualitative point of view (increase, decrease, or no variation). Changes in translational efficiency observed in other trypanosomatids strongly suggest that the differences are due to translational regulation and regulation of the steady-state protein levels. The lack of low-input sample strategies does not allow translatome and proteome analysis of sand fly–derived promastigotes so far.

The successes and future prospects of the linear antisense RNA amplification methodology

This Perspective discusses the development of the linear amplified RNA amplification technique over the last 25 years, and future applications of this important and versatile methodology.

It has been over a quarter of a century since the introduction of the linear RNA amplification methodology known as antisense RNA (aRNA) amplification. Whereas most molecular biology techniques are rapidly replaced owing to the fast-moving nature of development in the field, the aRNA procedure has become a base that can be built upon through varied uses of the technology. The technique was originally developed to assess RNA populations from small amounts of starting material, including single cells, but over time its use has evolved to include the detection of various cellular entities such as proteins, RNA-binding-protein-associated cargoes, and genomic DNA. In this Perspective we detail the linear aRNA amplification procedure and its use in assessing various components of a cell's chemical phenotype. This procedure is particularly useful in efforts to multiplex the simultaneous detection of various cellular processes. These efforts are necessary to identify the quantitative chemical phenotype of cells that underlies cellular function.

Trypanosoma cruzi specific mRNA amplification by in vitro transcription improves parasite transcriptomics in host-parasite RNA mixtures

Background

Trypanosomatids are a group of protozoan parasites that includes the etiologic agents of important human illnesses as Chagas disease, sleeping sickness and leishmaniasis. These parasites have a significant distinction from other eukaryotes concerning mRNA structure, since all mature mRNAs have an identical species-specific sequence of 39 nucleotides at the 5′ extremity, named spliced leader (SL). Considering this peculiar aspect of trypanosomatid mRNA, the aim of the present work was to develop a Trypanosoma cruzi specific in vitro transcription (IVT) linear mRNA amplification method in order to improve parasite transcriptomics analyses.

Methods

We designed an oligonucleotide complementary to the last 21 bases of T. cruzi SL sequence, bearing an upstream T7 promoter (T7SL primer), which was used to direct the synthesis of second-strand cDNA. Original mRNA was then amplified by IVT using T7 RNA polymerase. T7SL-amplified RNA from two distinct T. cruzi stages (epimastigotes and trypomastigotes) were deep sequenced in SOLiD platform. Usual poly(A) + RNA and and T7-oligo(dT) amplified RNA (Eberwine method) were also sequenced. RNA-Seq reads were aligned to our new and improved T. cruzi Dm28c genome assembly (PacBio technology) and resulting transcriptome pattern from these three RNA preparation methods were compared, mainly concerning the conservation of mRNA transcritional levels and DEGs detection between epimastigotes and trypomastigotes.

Results

T7SL IVT method detected more potential differentially expressed genes in comparison to either poly(A) + RNA or T7dT IVT, and was also able to produce reliable quantifications of the parasite transcriptome down to 3 ng of total RNA. Furthermore, amplification of parasite mRNA in HeLa/epimastigote RNA mixtures showed that T7SL IVT generates transcriptome quantification with similar detection of differentially expressed genes when parasite RNA mass was only 0.1% of the total mixture (R = 0.78 when compared to poly(A) + RNA).

Conclusions

The T7SL IVT amplification method presented here allows the detection of more potential parasite differentially expressed genes (in comparison to poly(A) + RNA) in host-parasite mixtures or samples with low amount of RNA. This method is especially useful for trypanosomatid transcriptomics because it produces less bias than PCR-based mRNA amplification. Additionally, by simply changing the complementary region of the T7SL primer, the present method can be applied to any trypanosomatid species.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4163-y) contains supplementary material, which is available to authorized users.

Hypoxic Air Inhalation and Ischemia Interventions Both Elicit Preconditioning Which Attenuate Subsequent Cellular Stress In vivo Following Blood Flow Occlusion and Reperfusion

Ischemic preconditioning (IPC) is valid technique which elicits reductions in femoral blood flow occlusion mediated reperfusion stress (oxidative stress, Hsp gene transcripts) within the systemic blood circulation and/or skeletal muscle. It is unknown whether systemic hypoxia, evoked by hypoxic preconditioning (HPC) has efficacy in priming the heat shock protein (Hsp) system thus reducing reperfusion stress following blood flow occlusion, in the same manner as IPC. The comparison between IPC and HPC being relevant as a preconditioning strategy prior to orthopedic surgery. In an independent group design, 18 healthy men were exposed to 40 min of (1) passive whole-body HPC (FiO₂ = 0.143; no ischemia. N = 6), (2) IPC (FiO₂ = 0.209; four bouts of 5 min ischemia and 5 min reperfusion. n = 6), or (3) rest (FiO₂ = 0.209; no ischemia. n = 6). The interventions were administered 1 h prior to 30 min of tourniquet derived femoral blood flow occlusion and were followed by 2 h subsequent reperfusion. Systemic blood samples were taken pre- and post-intervention. Systemic blood and gastrocnemius skeletal muscle samples were obtained pre-, 15 min post- (15PoT) and 120 min (120PoT) post-tourniquet deflation. To determine the cellular stress response gastrocnemius and leukocyte Hsp72 mRNA and Hsp32 mRNA gene transcripts were determined by RT-qPCR. The plasma oxidative stress response (protein carbonyl, reduced glutathione/oxidized glutathione ratio) was measured utilizing commercially available kits. In comparison to control, at 15PoT a significant difference in gastrocnemius Hsp72 mRNA was seen in HPC (−1.93-fold; p = 0.007) and IPC (−1.97-fold; p = 0.006). No significant differences were observed in gastrocnemius Hsp32 and Hsp72 mRNA, leukocyte Hsp72 and Hsp32 mRNA, or oxidative stress markers (p > 0.05) between HPC and IPC. HPC provided near identical amelioration of blood flow occlusion mediated gastrocnemius stress response (Hsp72 mRNA), compared to an established IPC protocol. This was seen independent of changes in systemic oxidative stress, which likely explains the absence of change in Hsp32 mRNA transcripts within leukocytes and the gastrocnemius. Both the established IPC and novel HPC interventions facilitate a priming of the skeletal muscle, but not leukocyte, Hsp system prior to femoral blood flow occlusion. This response demonstrates a localized tissue specific adaptation which may ameliorate reperfusion stress.

Working with Human Tissues for Translational Cancer Research

Medical research for human benefit is greatly impeded by the necessity for human tissues and subjects. However, upon obtaining consent for human specimens, precious samples must be handled with the greatest care in order to ensure integrity of organs, tissues, and cells to the highest degree. Unfortunately, tissue processing by definition requires extraction of tissues from the host, a change which can cause great cellular stress and have major repercussions on subsequent analyses. These stresses could result in the specimen being no longer representative of the site from which it was retrieved. Therefore, a strict protocol must be adhered to while processing these specimens to ensure representativeness. The desired assay(s) must also be taken into consideration in order to ensure that an optimal technique is used for sample processing. Outlined here is a protocol for tissue retrieval, processing and various analyses which may be performed on processed tissue in order to maximize downstream production from limited tissue samples.

Proinflammatory Cytokine IL-6 and JAK-STAT Signaling Pathway in Myeloproliferative Neoplasms

The recent JAK1/2 inhibitor trial in myeloproliferative neoplasms (MPNs) showed that reducing inflammation can be more beneficial than targeting gene mutants. We evaluated the proinflammatory IL-6 cytokine and JAK-STAT signaling pathway related genes in circulating CD34⁺ cells of MPNs. Regarding laboratory data, leukocytosis has been observed in polycythemia vera (PV) and JAK2V617F mutation positive versus negative primary myelofibrosis (PMF) patients. Moreover, thrombocytosis was reduced by JAK2V617F allele burden in essential thrombocythemia (ET) and PMF. 261 significantly changed genes have been detected in PV, 82 in ET, and 94 genes in PMF. The following JAK-STAT signaling pathway related genes had augmented expression in CD34⁺ cells of MPNs: CCND3 and IL23A regardless of JAK2V617F allele burden; CSF3R, IL6ST, and STAT1/2 in ET and PV with JAK2V617F mutation; and AKT2, IFNGR2, PIM1, PTPN11, and STAT3 only in PV. STAT5A gene expression was generally reduced in MPNs. IL-6 cytokine levels were increased in plasma, as well as IL-6 protein levels in bone marrow stroma of MPNs, dependent on JAK2V617F mutation presence in ET and PMF patients. Therefore, the JAK2V617F mutant allele burden participated in inflammation biomarkers induction and related signaling pathways activation in MPNs.

Microarray RNA transcriptional profiling: Part I. Platforms, experimental design and standardization

This review summarizes, in a balanced and comprehensive manner, the various components of microarrays and their types, substrate architecture, platforms for microarray probe implementation, standardizations and confounders. The review is intended to familiarize the beginner with the principles of experimental design and the selection of an appropriate microarray platform. This parallel technology has revolutionized transcriptomic approaches to data profiling and has a major role in the identification of expressed genes, classification and diagnosis studies. The technology is still evolving and guidelines for standardization and reporting have been developed and are being improved.

Different responses to mechanical injury in neonatal and adult ovine articular cartilage

Background

Articular cartilage injury remains a major challenge in orthopedic surgery. This study aimed to identify differences in gene expression and molecular responses between neonatal and adult articular cartilage during the healing of an injury.

Methods

An established in vitro model was used to compare the transcriptional response to cartilage injury in neonatal and adult sheep by microarray analysis of gene expression. Total RNA was isolated from tissue samples, linearly amplified, and 15,208 ovine probes were applied to cDNA microarray. Validation for selected genes was obtained by real-time quantitative polymerase chain reaction (RT-qPCR).

Results

We found 1,075 (11.6%) differentially expressed probe sets in adult injured cartilage relative to normal cartilage. A total of 1,016 (11.0%) probe sets were differentially expressed in neonatal injured cartilage relative to normal cartilage. A total of 1,492 (16.1%) probe sets were differentially expressed in adult normal cartilage relative to neonatal normal cartilage. A total of 1,411 (15.3%) probe sets were differentially expressed in adult injured cartilage relative to neonatal injured cartilage. Significant functional clusters included genes associated with wound healing, articular protection, inflammation, and energy metabolism. Selected genes (PPARG, LDH, TOM, HIF1A, SMAD7, and NF-κB) were also found and validated by RT-qPCR.

Conclusions

There are significant differences in gene expression between neonatal and adult ovine articular cartilage following acute injury. They are partly due to intrinsic differences in the process of development, and partly to different biological responses to mechanical trauma between neonatal and adult articular cartilage.

Functional genomics of a generalist parasitic plant: laser microdissection of host-parasite interface reveals host-specific patterns of parasite gene expression

BACKGROUND

Orobanchaceae is the only plant family with members representing the full range of parasitic lifestyles plus a free-living lineage sister to all parasitic lineages, Lindenbergia. A generalist member of this family, and an important parasitic plant model, Triphysaria versicolor regularly feeds upon a wide range of host plants. Here, we compare de novo assembled transcriptomes generated from laser micro-dissected tissues at the host-parasite interface to uncover details of the largely uncharacterized interaction between parasitic plants and their hosts.

RESULTS

The interaction of Triphysaria with the distantly related hosts Zea mays and Medicago truncatula reveals dramatic host-specific gene expression patterns. Relative to above ground tissues, gene families are disproportionally represented at the interface including enrichment for transcription factors and genes of unknown function. Quantitative Real-Time PCR of a T. versicolor β-expansin shows strong differential (120x) upregulation in response to the monocot host Z. mays; a result that is concordant with our read count estimates. Pathogenesis-related proteins, other cell wall modifying enzymes, and orthologs of genes with unknown function (annotated as such in sequenced plant genomes) are among the parasite genes highly expressed by T. versicolor at the parasite-host interface.

CONCLUSIONS

Laser capture microdissection makes it possible to sample the small region of cells at the epicenter of parasite host interactions. The results of our analysis suggest that T. versicolor's generalist strategy involves a reliance on overlapping but distinct gene sets, depending upon the host plant it is parasitizing. The massive upregulation of a T. versicolor β-expansin is suggestive of a mechanism for parasite success on grass hosts. In this preliminary study of the interface transcriptomes, we have shown that T. versicolor, and the Orobanchaceae in general, provide excellent opportunities for the characterization of plant genes with unknown functions.

Optimization and improvement of oligonucleotide microarray-based detection of tomato viruses and pospiviroids

Tomato (Solanum lycopersicum L.) is a vegetable crop which is affected by many viruses and several viroids, causing significant economic loss. Their detection and identification is of critical importance for plant protection and quarantine and certification programs. The potential was examined of an array based on the Combimatrix platform for the detection of 37 viruses belonging to 13 families, one of which is unassigned, together with six pospiviroid species, genus Pospiviroid, family Pospiviroidae. More than 470 oligonucleotide probes (40-mer) were selected for the microarray diagnostic technique developed in this investigation. Most of the virus probes were highly specific and were able to identify tomato viruses. Most pospiviroid probes, however, were non-specific in terms of species, but were specific at the genus level as they hybridized to members of the genus Pospiviroid. Only one probe of the Tomato apical stem viroid was species specific. The repeatability and specificity of the Combimatrix method showed that it can be considered for routine diagnostic use in suspected tomato germplasm since it detected 37 viruses and one pospiviroid at the species level and 5-6 pospiviroids at the genus level. The estimated cost for testing of a single tomato virus is similar to or less than the cost of using ELISA.

Gene expression profiling of liver cancer stem cells by RNA-sequencing

BACKGROUND

Accumulating evidence supports that tumor growth and cancer relapse are driven by cancer stem cells. Our previous work has demonstrated the existence of CD90(+) liver cancer stem cells (CSCs) in hepatocellular carcinoma (HCC). Nevertheless, the characteristics of these cells are still poorly understood. In this study, we employed a more sensitive RNA-sequencing (RNA-Seq) to compare the gene expression profiling of CD90(+) cells sorted from tumor (CD90(+)CSCs) with parallel non-tumorous liver tissues (CD90(+)NTSCs) and elucidate the roles of putative target genes in hepatocarcinogenesis.

METHODOLOGY/PRINCIPAL FINDINGS

CD90(+) cells were sorted respectively from tumor and adjacent non-tumorous human liver tissues using fluorescence-activated cell sorting. The amplified RNAs of CD90(+) cells from 3 HCC patients were subjected to RNA-Seq analysis. A differential gene expression profile was established between CD90(+)CSCs and CD90(+)NTSCs, and validated by quantitative real-time PCR (qRT-PCR) on the same set of amplified RNAs, and further confirmed in an independent cohort of 12 HCC patients. Five hundred genes were differentially expressed (119 up-regulated and 381 down-regulated genes) between CD90(+)CSCs and CD90(+)NTSCs. Gene ontology analysis indicated that the over-expressed genes in CD90(+)CSCs were associated with inflammation, drug resistance and lipid metabolism. Among the differentially expressed genes, glypican-3 (GPC3), a member of glypican family, was markedly elevated in CD90(+)CSCs compared to CD90(+)NTSCs. Immunohistochemistry demonstrated that GPC3 was highly expressed in forty-two human liver tumor tissues but absent in adjacent non-tumorous liver tissues. Flow cytometry indicated that GPC3 was highly expressed in liver CD90(+)CSCs and mature cancer cells in liver cancer cell lines and human liver tumor tissues. Furthermore, GPC3 expression was positively correlated with the number of CD90(+)CSCs in liver tumor tissues.

CONCLUSIONS/SIGNIFICANCE

The identified genes, such as GPC3 that are distinctly expressed in liver CD90(+)CSCs, may be promising gene candidates for HCC therapy without inducing damages to normal liver stem cells.

Global analysis of cell type-specific gene expression

The tissues and organs of multicellular eukaryotes are frequently observed to comprise complex three-dimensional interspersions of different cell types. It is a reasonable assumption that different global patterns of gene expression are found within these different cell types. This review outlines general experimental strategies designed to characterize these global gene expression patterns, based on a combination of methods of transgenic fluorescent protein (FP) expression and targeting, of flow cytometry and sorting and of high-throughput gene expression analysis.

Transcriptomic analysis of grain amaranth (Amaranthus hypochondriacus) using 454 pyrosequencing: comparison with A. tuberculatus, expression profiling in stems and in response to biotic and abiotic stress

BACKGROUND

Amaranthus hypochondriacus, a grain amaranth, is a C4 plant noted by its ability to tolerate stressful conditions and produce highly nutritious seeds. These possess an optimal amino acid balance and constitute a rich source of health-promoting peptides. Although several recent studies, mostly involving subtractive hybridization strategies, have contributed to increase the relatively low number of grain amaranth expressed sequence tags (ESTs), transcriptomic information of this species remains limited, particularly regarding tissue-specific and biotic stress-related genes. Thus, a large scale transcriptome analysis was performed to generate stem- and (a)biotic stress-responsive gene expression profiles in grain amaranth.

RESULTS

A total of 2,700,168 raw reads were obtained from six 454 pyrosequencing runs, which were assembled into 21,207 high quality sequences (20,408 isotigs + 799 contigs). The average sequence length was 1,064 bp and 930 bp for isotigs and contigs, respectively. Only 5,113 singletons were recovered after quality control. Contigs/isotigs were further incorporated into 15,667 isogroups. All unique sequences were queried against the nr, TAIR, UniRef100, UniRef50 and Amaranthaceae EST databases for annotation. Functional GO annotation was performed with all contigs/isotigs that produced significant hits with the TAIR database. Only 8,260 sequences were found to be homologous when the transcriptomes of A. tuberculatus and A. hypochondriacus were compared, most of which were associated with basic house-keeping processes. Digital expression analysis identified 1,971 differentially expressed genes in response to at least one of four stress treatments tested. These included several multiple-stress-inducible genes that could represent potential candidates for use in the engineering of stress-resistant plants. The transcriptomic data generated from pigmented stems shared similarity with findings reported in developing stems of Arabidopsis and black cottonwood (Populus trichocarpa).

CONCLUSIONS

This study represents the first large-scale transcriptomic analysis of A. hypochondriacus, considered to be a highly nutritious and stress-tolerant crop. Numerous genes were found to be induced in response to (a)biotic stress, many of which could further the understanding of the mechanisms that contribute to multiple stress-resistance in plants, a trait that has potential biotechnological applications in agriculture.

Brain transcriptomic response of threespine sticklebacks to cues of a predator

Predation pressure represents a strong selective force that influences the development and evolution of living organisms. An increasing number of studies have shown that both environmental and social factors, including exposure to predators, substantially shape the structure and function of the brain. However, our knowledge about the molecular mechanisms underlying the response of the brain to environmental stimuli is limited. In this study, we used whole-genome comparative oligonucleotide microarrays to investigate the brain transcriptomic response to cues of a predator in the threespine stickleback, Gasterosteus aculeatus. We found that repeated exposure to olfactory, visual and tactile cues of a predator (rainbow trout, Oncorrhynchus mykiss) for 6 days resulted in subtle but significant transcriptomic changes in the brain of sticklebacks. Gene functional analysis and gene ontology enrichment revealed that the majority of the transcripts differentially expressed between the fish exposed to cues of a predator and the control group were related to antigen processing and presentation involving the major histocompatibility complex, transmission of synaptic signals, brain metabolic processes, gene regulation and visual perception. The top four identified pathways were synaptic long-term depression, RAN signaling, relaxin signaling and phototransduction. Our study demonstrates that exposure of sticklebacks to cues of a predator results in the activation of a wide range of biological and molecular processes and lays the foundation for future investigations on the molecular factors that modulate the function and evolution of the brain in response to stressors.

Functional comparison of innate immune signaling pathways in primates

Humans respond differently than other primates to a large number of infections. Differences in susceptibility to infectious agents between humans and other primates are probably due to inter-species differences in immune response to infection. Consistent with that notion, genes involved in immunity-related processes are strongly enriched among recent targets of positive selection in primates, suggesting that immune responses evolve rapidly, yet providing only indirect evidence for possible inter-species functional differences. To directly compare immune responses among primates, we stimulated primary monocytes from humans, chimpanzees, and rhesus macaques with lipopolysaccharide (LPS) and studied the ensuing time-course regulatory responses. We find that, while the universal Toll-like receptor response is mostly conserved across primates, the regulatory response associated with viral infections is often lineage-specific, probably reflecting rapid host–virus mutual adaptation cycles. Additionally, human-specific immune responses are enriched for genes involved in apoptosis, as well as for genes associated with cancer and with susceptibility to infectious diseases or immune-related disorders. Finally, we find that chimpanzee-specific immune signaling pathways are enriched for HIV–interacting genes. Put together, our observations lend strong support to the notion that lineage-specific immune responses may help explain known inter-species differences in susceptibility to infectious diseases.

We know of a large number of diseases or medical conditions that affect humans more severely than non-human primates, such as AIDS, malaria, hepatitis B, and cancer. These differences likely arise from different immune responses to infection among species. However, due to the lack of comparative functional data across species, it remains unclear how the immune system of humans and other primates differ. In this work, we present the first genome-wide characterization of functional differences in innate immune responses between humans and our closest evolutionary relatives. Our results indicate that “core” immune responses, those that are critical to fight any invading pathogen, are the most conserved across primates and that much of the divergence in immune responses is observed in genes that are involved in response to specific microbial and viral agents. In addition, we show that human-specific immune responses are enriched for genes involved in apoptosis and cancer biology, as well as with genes previously associated with susceptibility to infectious diseases or immune-related disorders. Finally, we find that chimpanzee-specific immune signaling pathways are enriched for HIV–interacting genes. Our observations may therefore help explain known inter-species differences in susceptibility to infectious diseases.

Location-specific responses to thermal stress in larvae of the reef-building coral Montastraea faveolata

BACKGROUND

The potential to adapt to a changing climate depends in part upon the standing genetic variation present in wild populations. In corals, the dispersive larval phase is particularly vulnerable to the effects of environmental stress. Larval survival and response to stress during dispersal and settlement will play a key role in the persistence of coral populations.

METHODOLOGY/PRINCIPAL FINDINGS

To test the hypothesis that larval transcription profiles reflect location-specific responses to thermal stress, symbiont-free gametes from three to four colonies of the scleractinian coral Montastraea faveolata were collected from Florida and Mexico, fertilized, and raised under mean and elevated (up 1 to 2 degrees C above summer mean) temperatures. These locations have been shown to exchange larvae frequently enough to prevent significant differentiation of neutral loci. Differences among 1,310 unigenes were simultaneously characterized using custom cDNA microarrays, allowing investigation of gene expression patterns among larvae generated from wild populations under stress. Results show both conserved and location-specific variation in key processes including apoptosis, cell structuring, adhesion and development, energy and protein metabolism, and response to stress, in embryos of a reef-building coral.

CONCLUSIONS/SIGNIFICANCE

These results provide first insights into location-specific variation in gene expression in the face of gene flow, and support the hypothesis that coral host genomes may house adaptive potential needed to deal with changing environmental conditions.

Transcriptional comparisons between equine articular repair tissue, neonatal cartilage, cultured chondrocytes and mesenchymal stromal cells

Human and equine cell transplant strategies for cartilage lesions usually result in scar tissue that is similar to what is produced naturally during the repair process. In this study, culture-expanded de-differentiated chondrocytes and primary bone marrow stromal cells at a pre-transplantation time-point were compared along with neonatal cartilage to repair tissue. Transcriptional profiling using a 9413-probeset equine-specific cDNA microarray and targeted real-time quantitative polymerase chain reaction validation were used to characterize relationships between these cell types and repair tissue both broadly and for individual cartilage biomarkers. The greatest divergence in expression was detected for transcripts encoding matrix proteins that typically define the differentiation status of normal articular cartilage and fibrocartilage repair tissue. Expression patterns and gene ontology analyses indicated that while the repair cells were more chondrogenic than bone marrow stromal cells and de-differentiated cultured chondrocytes, steady-state levels of transcripts encoding cartilage biomarkers were substantially lower than the amounts found in neonatal articular cartilage. By characterizing gene expression differences amongst these tissues, we present important targets to monitor when developing improvements to cartilage engineering therapies.

Reproducible RNA preparation from sugarcane and citrus for functional genomic applications

High-throughput functional genomic procedures depend on the quality of the RNA used. Copurifying molecules can negatively impact the functionality of some plant RNA preparations employed in these procedures. We present a simplified, rapid, and scalable SDS/phenol-based method that provides the high-quantity and -quality RNA required by the newly emerging biotechnology applications. The method is applied to isolating RNA from tissues of two biotechnologically important crop plants, sugarcane and citrus, which provide a challenge due to the presence of fiber, polysaccharides, or secondary metabolites. The RNA isolated by this method is suitable for several downstream applications including northern blot hybridization, microarray analysis, and quantitative RT-PCR. This method has been used in a diverse range of projects ranging from screening plant lines overexpressing mammalian genes to analyzing plant responses to viral infection and defense signaling molecules.

Manual microdissection combined with antisense RNA-longSAGE for the analysis of limited cell numbers

Establishing a gene expression profile of defined subtypes of cells within an organ is still challenging because it frequently requires microdissection and subsequent amplification of the limited amount of messenger RNA (mRNA) isolated from the microdissected tissue in order to be able to proceed with comprehensive gene expression analyses via microarray or serial analysis of gene expression (SAGE) technology. Here we describe a manual microdissection strategy for the isolation of high-quality RNA. Furthermore, a strategy for combining linear amplification of RNA with longSAGE is described that allows the use of antisense RNA (aRNA) generated via the well-established linear amplification of RNA procedure together with the conventional SAGE or longSAGE technology.

Analysis of the degradome with the CLIP-CHIP microarray

The degradome microarray - CLIP-CHIP - is a dedicated and focused array that allows the analysis of all proteases, non-proteolytic homologs, and protease inhibitor gene transcripts in the human and murine genomes at the mRNA transcript level.Based on unique 70-mer oligonucleotides, designed to match parts of the sequence of known or predicted protease and inhibitor mRNAs in both species and printed on a glass-matrix surface, the CLIP-CHIP microarray can be used to analyze differentially expressed protease and inhibitor gene products and give expression profiles for any analyzed sample.

Gene expression using microarrays in transplant recipients at risk of EBV lymphoproliferation after organ transplantation: preliminary proof-of-concept

We hypothesized that aspects of the virus-host interaction could be measured to help predict progression to EBV-PTLD. We examined the spectrum of host genes differentially expressed and any relevant clustering in children at risk of EBV lymphoproliferation after organ transplantation. We compared the genes expressed among patients with different levels of viral loads. Gene expression was measured by microarray analysis of RNA from CD19+ B lymphocytes using the Human Genome U133 Plus 2.0 GeneChip. Among 27 samples from 26 transplant recipients, the viral load categories were: low or undetectable loads (LU), n = 14; high or intermediate loads (HI), n = 13. There were seven healthy EBV-seropositive (P) and -seronegative controls (N). Median time of post-transplantation was 0.5 yr (range 0.1-3.8). We identified 24-54 differentially expressed genes in each of four comparisons of HI vs. P, LU vs. P, HI vs. LU, and P vs. N. We identified patterns of 563 gene expressions, creating five clusters aligned with levels of viral load. PTLD occurred in four of five clusters. In summary, we demonstrated varying degrees of alignment between levels of VL and gene clusters. Analyses for differential expression of genes showed genes that could be implicated in the pathogenesis of EBV-PTLD.

Transcriptional profiling differences for articular cartilage and repair tissue in equine joint surface lesions

Background

Full-thickness articular cartilage lesions that reach to the subchondral bone yet are restricted to the chondral compartment usually fill with a fibrocartilage-like repair tissue which is structurally and biomechanically compromised relative to normal articular cartilage. The objective of this study was to evaluate transcriptional differences between chondrocytes of normal articular cartilage and repair tissue cells four months post-microfracture.

Methods

Bilateral one-cm² full-thickness defects were made in the articular surface of both distal femurs of four adult horses followed by subchondral microfracture. Four months postoperatively, repair tissue from the lesion site and grossly normal articular cartilage from within the same femorotibial joint were collected. Total RNA was isolated from the tissue samples, linearly amplified, and applied to a 9,413-probe set equine-specific cDNA microarray. Eight paired comparisons matched by limb and horse were made with a dye-swap experimental design with validation by histological analyses and quantitative real-time polymerase chain reaction (RT-qPCR).

Results

Statistical analyses revealed 3,327 (35.3%) differentially expressed probe sets. Expression of biomarkers typically associated with normal articular cartilage and fibrocartilage repair tissue corroborate earlier studies. Other changes in gene expression previously unassociated with cartilage repair were also revealed and validated by RT-qPCR.

Conclusion

The magnitude of divergence in transcriptional profiles between normal chondrocytes and the cells that populate repair tissue reveal substantial functional differences between these two cell populations. At the four-month postoperative time point, the relative deficiency within repair tissue of gene transcripts which typically define articular cartilage indicate that while cells occupying the lesion might be of mesenchymal origin, they have not recapitulated differentiation to the chondrogenic phenotype of normal articular chondrocytes.

Reproducibility of quantitative RT-PCR array in miRNA expression profiling and comparison with microarray analysis

BACKGROUND

MicroRNAs (miRNAs) have critical functions in various biological processes. MiRNA profiling is an important tool for the identification of differentially expressed miRNAs in normal cellular and disease processes. A technical challenge remains for high-throughput miRNA expression analysis as the number of miRNAs continues to increase with in silico prediction and experimental verification. Our study critically evaluated the performance of a novel miRNA expression profiling approach, quantitative RT-PCR array (qPCR-array), compared to miRNA detection with oligonucleotide microchip (microarray).

RESULTS

High reproducibility with qPCR-array was demonstrated by comparing replicate results from the same RNA sample. Pre-amplification of the miRNA cDNA improved sensitivity of the qPCR-array and increased the number of detectable miRNAs. Furthermore, the relative expression levels of miRNAs were maintained after pre-amplification. When the performance of qPCR-array and microarrays were compared using different aliquots of the same RNA, a low correlation between the two methods (r=-0.443) indicated considerable variability between the two assay platforms. Higher variation between replicates was observed in miRNAs with low expression in both assays. Finally, a higher false positive rate of differential miRNA expression was observed using the microarray compared to the qPCR-array.

CONCLUSION

Our studies demonstrated high reproducibility of TaqMan qPCR-array. Comparison between different reverse transcription reactions and qPCR-arrays performed on different days indicated that reverse transcription reactions did not introduce significant variation in the results. The use of cDNA pre-amplification increased the sensitivity of miRNA detection. Although there was variability associated with pre-amplification in low abundance miRNAs, the latter did not involve any systemic bias in the estimation of miRNA expression. Comparison between microarray and qPCR-array indicated superior sensitivity and specificity of qPCR-array.

Systemic treatment of xenografts with vaccinia virus GLV-1h68 reveals the immunologic facet of oncolytic therapy

BACKGROUND

GLV-1h68 is an attenuated recombinant vaccinia virus (VACV) that selectively colonizes established human xenografts inducing their complete regression.

RESULTS

Here, we explored xenograft/VACV/host interactions in vivo adopting organism-specific expression arrays and tumor cell/VACV in vitro comparing VACV replication patterns. There were no clear-cut differences in vitro among responding and non-responding tumors, however, tumor rejection was associated in vivo with activation of interferon-stimulated genes (ISGs) and innate immune host's effector functions (IEFs) correlating with VACV colonization of the xenografts. These signatures precisely reproduce those observed in humans during immune-mediated tissue-specific destruction (TSD) that causes tumor or allograft rejection, autoimmunity or clearance of pathogens. We recently defined these common pathways in the "immunologic constant of rejection" hypothesis (ICR).

CONCLUSION

This study provides the first prospective validation of a universal mechanism associated with TSD. Thus, xenograft infection by oncolytic VACV, beyond offering a promising therapy of established cancers, may represent a reliable pre-clinical model to test therapeutic strategies aimed at modulating the central pathways leading to TSD; this information may lead to the identification of principles that could refine the treatment of cancer and chronic infection by immune stimulation or autoimmunity and allograft rejection through immune tolerance.

Assessment of the swine protein-annotated oligonucleotide microarray

The specificity and utility of the swine protein-annotated oligonucleotide microarray, or Pigoligoarray (http://www.pigoligoarray.org), has been evaluated by profiling the expression of transcripts from four porcine tissues. Tools for comparative analyses of expression on the Pigoligoarray were developed including HGNC identities and comparative mapping alignments with human orthologs. Hybridization results based on the Pigoligoarray's sets of control, perfect match (PM) and deliberate mismatch (MM) probes provide an important means of assessing non-specific hybridization. Simple descriptive diagnostic analyses of PM/MM probe sets are introduced in this paper as useful tools for detecting non-specific hybridization. Samples of RNA from liver, brain stem, longissimus dorsi muscle and uterine endothelium from four pigs were prepared and hybridized to the arrays. Of the total 20,400 oligonucleotides on the Pigoligoarray, 12,429 transcripts were putatively differentially expressed (DE). Analyses for tissue-specific expression [over-expressed in one tissue with respect to all the remaining three tissues (q < 0.01)] identified 958 DE transcripts in liver, 726 in muscle, 286 in uterine endothelium and 1027 in brain stem. These hybridization results were confirmed by quantitative PCR (QPCR) expression patterns for a subset of genes after affirming that cDNA and amplified antisense RNA (aRNA) exhibited similar QPCR results. Comparison to human ortholog expression confirmed the value of this array for experiments of both agricultural importance and for tests using pigs as a biomedical model for human disease.

Novel serum markers of fibrosis progression for the follow-up of hepatitis C virus-infected patients

Liver biopsy is considered the gold-standard method for the assessment of liver fibrosis during follow-up of hepatitis C virus-infected patients, but this invasive procedure is not devoid of complications. The aim of the present study was to identify novel non-invasive markers of fibrosis progression. By microarray analysis, we compared transcript levels in two extreme stages of fibrosis from 16 patients. Informative transcripts were validated by real-time PCR and used for the assessment of fibrosis in 23 additional patients. Sixteen transcripts were found to be dysregulated during the fibrogenesis process. Among them, some were of great interest because their corresponding proteins could be serologically measured. Thus, the protein levels of inter-alpha inhibitor H1, serpin peptidase inhibitor clade F member 2, and transthyretin were all significantly different according to the four Metavir stages of fibrosis. In conclusion, we report here that dysregulation, at both the transcriptional and protein levels, exists during the fibrogenesis process. Our description of three novel serum markers and their potential use as serological tests for the non-invasive diagnosis of liver fibrosis open new opportunities for better follow-up of hepatitis C virus-infected patients.

From mouse to man: the 5-HT3 receptor modulates physical dependence on opioid narcotics

OBJECTIVES

Addiction to opioid narcotics represents a major public health challenge. Animal models of one component of addiction, physical dependence, show this trait to be highly heritable. The analysis of opioid dependence using contemporary in-silico techniques offers an approach to discover novel treatments for dependence and addiction.

METHODS

In these experiments, opioid withdrawal behavior in 18 inbred strains of mice was assessed. Mice were treated for 4 days with escalating doses of morphine before the administration of naloxone allowing the quantification of opioid dependence. After haplotypic analysis, experiments were designed to evaluate the top gene candidate as a modulator of physical dependence. Behavioral studies as well as measurements of gene expression on the mRNA and protein levels were completed. Finally, a human model of opioid dependence was used to quantify the effects of the 5-HT3 antagonist ondansetron on signs and symptoms of withdrawal.

RESULTS

The Htr3a gene corresponding to the 5-HT3 receptor emerged as the leading candidate. Pharmacological studies using the selective 5-HT3 antagonist ondansetron supported the link in mice. Morphine strongly regulated the expression of the Htr3a gene in various central nervous system regions including the amygdala, dorsal raphe, and periaqueductal gray nuclei, which have been linked to opioid dependence in previous studies. Using an acute morphine administration model, the role of 5-HT3 in controlling the objective signs of withdrawal in humans was confirmed.

CONCLUSION

These studies show the power of in-silico genetic mapping, and reveal a novel target for treating an important component of opioid addiction.

Synthesis of biotin-labeled RNA for gene expression measurements using oligonucleotide arrays

Using gene arrays, it is currently possible to simultaneously measure mRNA levels of many genes in any tissue of interest. Undoubtedly, comprehensive measurements of gene expression as part of carefully designed experiments will continue to further our understanding of audition and have the potential to open up new avenues of research. This chapter describes a reliable protocol to prepare high-quality biotin-labeled RNA target, specifically for oligonucleotide array experiments. The procedure includes isolation of high-quality total RNA, synthesis of double-stranded cDNA engineered for in vitro transcription with T7 RNA polymerase, subsequent in vitro transcription in the presence of biotin-labeled ribonucleotides, and fractionation of the RNA to approximately 500 bp fragments, suitable for oligonucleotide array experiments. Because the membranous labyrinth is composed of functionally interdependent cellular structures, which themselves contain numerous, highly differentiated cell types, comprehensive analysis of gene expression in the cochlea is best complemented by immunohistochemical studies or, if no suitable antibodies are available, by in situ hybridization studies. Either one of these techniques will identify the specific cell types that express the genes of interests.

Differential gene expression associated with postnatal equine articular cartilage maturation

Background

Articular cartilage undergoes an important maturation process from neonate to adult that is reflected by alterations in matrix protein organization and increased heterogeneity of chondrocyte morphology. In the horse, these changes are influenced by exercise during the first five months of postnatal life. Transcriptional profiling was used to evaluate changes in articular chondrocyte gene expression during postnatal growth and development.

Methods

Total RNA was isolated from the articular cartilage of neonatal (0–10 days) and adult (4–5 years) horses, subjected to one round of linear RNA amplification, and then applied to a 9,367-element equine-specific cDNA microarray. Comparisons were made with a dye-swap experimental design. Microarray results for selected genes (COL2A1, COMP, P4HA1, TGFB1, TGFBR3, TNC) were validated by quantitative polymerase chain reaction (qPCR).

Results

Fifty-six probe sets, which represent 45 gene products, were up-regulated (p < 0.01) in chondrocytes of neonatal articular cartilage relative to chondrocytes of adult articular cartilage. Conversely, 586 probe sets, which represent 499 gene products, were up-regulated (p < 0.01) in chondrocytes of adult articular cartilage relative to chondrocytes of neonatal articular cartilage. Collagens, matrix-modifying enzymes, and provisional matrix non-collagenous proteins were expressed at higher levels in the articular cartilage of newborn foals. Those genes with increased mRNA abundance in adult chondrocytes included leucine-rich small proteoglycans, matrix assembly, and cartilage maintenance proteins.

Conclusion

Differential expression of genes encoding matrix proteins and matrix-modifying enzymes between neonates and adults reflect a cellular maturation process in articular chondrocytes. Up-regulated transcripts in neonatal cartilage are consistent with growth and expansion of the articular surface. Expression patterns in mature articular cartilage indicate a transition from growth to homeostasis, and tissue function related to withstanding shear and weight-bearing stresses.

Gene expression profiling of single epidermal, basal and trichome cells of Arabidopsis thaliana

Samples of single epidermal, basal and trichome cells were collected by glass microcapillaries from 7-week-old Arabidopsis thaliana leaves. Transcript amplification of these single-cell samples was performed by RT PCR. For gene expression profiling, we hybridized the amplified transcriptome of each individual cell type to nylon membranes spotted with 16,000 Arabidopsis expressed sequence tags (ESTs). Initial analysis of the array filter data enabled us to functionally categorize transcripts that were present in each individual cell type. In order to confirm the filter array data, we used RT PCR. Results of this RT PCR approach confirmed the presence of 12 selected candidate genes in agreement with array filter hybridization data. Further, transcripts involved in detoxification and sulfur metabolism could be identified in epidermal cell extracts. Together, the results of our study provide the localization of approximately 1000 expressed genes to either pavement, basal or trichome cells. To cluster transcripts with similar expression levels, we developed a novel mathematical algorithm. Based on the mean and standard deviation, ratios of expression levels of a transcript were defined for pairs of the three cell types. This numerical analysis enabled subdivision into 67 categories of genes differentially expressed in epidermal, basal and trichome cells. Transcripts in each category displayed similar ratios of expression levels in the three cell types. Examples of these clusters are presented and discussed in Appendix A.

Improving the prediction accuracy in classification using the combined data sets by ranks of gene expressions

Background

The information from different data sets experimented under different conditions may be inconsistent even though they are performed with the same research objectives. More than that, even when the data sets were generated from the same platform, the data agreement may be affected by the technical variation among the laboratories. In this case, it is necessary to use the combined data set after adjusting the differences between such data sets, for detecting the more reliable information.

Results

The proposed method combines data sets posterior to the discretization of data sets based on the ranks of the gene expression ratios, and the statistical method is applied to the combined data set for predictive gene selection. The efficiency of the proposed method was evaluated using five colon cancer related data sets, which were experimented using cDNA microarrays with different RNA sources, and one experiment utilized oligonucleotide arrays. NCI-60 cell lines data sets were used, which were performed with two different platforms of cDNA microarrays and Affymetrix HU6800 oligonucleotide arrays. The combined data set by the proposed method predicted the test data sets more accurately than the separated data sets did. The biological significant genes were detected from the combined data set, which were missed on the separated data sets.

Conclusion

By transforming gene expressions using ranks, the proposed method is not influenced by systematic bias among chips and normalization method. The method may be especially more useful to find predictive genes from data sets which have different scale in gene expressions.

Comparison of RNA amplification methods and chip platforms for microarray analysis of samples processed by laser capture microdissection

Laser capture microdissection (LCM) permits isolation of pure cell populations from which RNA can be extracted, amplified, and subjected to microarray analysis, allowing information to be obtained on the gene expression profile of defined cell types. To avoid amplification artifacts and detect genes expressed at different levels, it is important to optimize the choice of both RNA amplification step and microarray platform. We captured by LCM the same colon cancer biopsy and conducted a cross comparison of distinct RNA amplification methods and different chip platforms. We tested two RNA amplification methods with different chemistry: the one-cycle Ovation system (NuGEN) and the two-cycle Ribo OA method (Arcturus). We also compared two different whole genome platforms, based on Affymetrix technology: the U133 plus 2.0 and the X3P array, with probe sets closer to the 3' end of transcripts. After RNA amplification, microarray analysis, and data normalization, we investigated reproducibility and correlation of different methods and arrays. Our results indicate that the Arcturus Ribo OA method is superior for both array choices, especially in combination with X3P arrays, showing the lowest variance and Spearman correlation of 0.986. The quicker NuGEN procedure, when coupled with X3P arrays, also yielded excellent results (correlation of 0.951). These observations will be useful for planning large-scale analyses of LCM-dissected clinical samples.

Microbial community gene expression in ocean surface waters

Metagenomics is expanding our knowledge of the gene content, functional significance, and genetic variability in natural microbial communities. Still, there exists limited information concerning the regulation and dynamics of genes in the environment. We report here global analysis of expressed genes in a naturally occurring microbial community. We first adapted RNA amplification technologies to produce large amounts of cDNA from small quantities of total microbial community RNA. The fidelity of the RNA amplification procedure was validated with Prochlorococcus cultures and then applied to a microbial assemblage collected in the oligotrophic Pacific Ocean. Microbial community cDNAs were analyzed by pyrosequencing and compared with microbial community genomic DNA sequences determined from the same sample. Pyrosequencing-based estimates of microbial community gene expression compared favorably to independent assessments of individual gene expression using quantitative PCR. Genes associated with key metabolic pathways in open ocean microbial species-including genes involved in photosynthesis, carbon fixation, and nitrogen acquisition-and a number of genes encoding hypothetical proteins were highly represented in the cDNA pool. Genes present in the variable regions of Prochlorococcus genomes were among the most highly expressed, suggesting these encode proteins central to cellular processes in specific genotypes. Although many transcripts detected were highly similar to genes previously detected in ocean metagenomic surveys, a significant fraction ( approximately 50%) were unique. Thus, microbial community transcriptomic analyses revealed not only indigenous gene- and taxon-specific expression patterns but also gene categories undetected in previous DNA-based metagenomic surveys.

aRNA-LongSAGE: SAGE with antisense RNA

In order to generate serial analysis of gene expression (SAGE) libraries from very small samples such as microdissected cells, the starting material must first be amplified via PCR or linear amplification of RNA. In microarray experiments, it has been shown that linear amplification of RNA can be used to generate reliable gene expression profiles and leads to the detection of expression differences that are not seen with nonamplified starting material. As the product of the amplification is amplified antisense RNA (aRNA), linear amplification of RNA cannot be used in combination with the conventional SAGE protocol. The aRNA-LongSAGE protocol described herein is an adaptation of the MicroSAGE protocol to the use of aRNA as starting material.

Conservation of genetic alterations in recurrent melanoma supports the melanoma stem cell hypothesis

It is generally accepted that human cancers derive from a mutated single cell. However, the genetic steps characterizing various stages of progression remain unclear. Studying a unique case of metastatic melanoma, we observed that cell lines derived from metachronous metastases arising over a decade retained a central core of genetic stability in spite of divergent phenotypes. In the present study, we expanded our previous observations comparing these autologous cell lines of clonal derivation with allogeneic ones and correlated array comparative genomic hybridization (aCGH) with gene expression profiling to determine their relative contribution to the dynamics of disease progression. aCGH and gene expression profiling were performed on autologous cell lines and allogeneic melanoma cell lines originating from other patients. A striking correlation existed between total extent of genetic imbalances, global transcriptional patterns, and cellular phenotypes. They did not follow a strict temporal progression but stemmed independently at various time points from a central core of genetic stability best explained according to the cancer stem cell hypothesis. Although their contribution was intertwined, genomic imbalances detectable by aCGH contributed only 25% of the transcriptional traits determining autologous tumor distinctiveness. Our study provides important insights about the dynamics of cancer progression and supports the development of targeted anticancer therapies aimed against stable genetic factors that are maintained throughout the end stage of disease.

Microarray analysis reveals distinctive signaling between the bed nucleus of the stria terminalis, nucleus accumbens, and dorsal striatum

To identify distinct transcriptional patterns between the major subcortical dopamine targets commonly studied in addiction we studied differences in gene expression between the bed nucleus of the stria terminalis (BNST), nucleus accumbens (NAc), and dorsal striatum (dStr) using microarray analysis. We first tested for differences in expression of genes encoding transcripts for common neurotransmitter systems as well as calcium binding proteins routinely used in neuroanatomical delineation of brain regions. This a priori method revealed differential expression of corticotropin releasing hormone ( Crh), the GABA transporter ( Slc6a1), and prodynorphin ( Pdyn) mRNAs as well as several others. Using a gene ontology tool, functional scoring analysis, and Ingenuity Pathway Analysis, we further identified several physiological pathways that were distinct among these brain regions. These two different analyses both identified calcium signaling, G-coupled protein receptor signaling, and adenylate cyclase-related signaling as significantly different among the BNST, NAc, and dStr. These types of signaling pathways play important roles in, amongst other things, synaptic plasticity. Investigation of differential gene expression revealed several instances that may provide insight into reported differences in synaptic plasticity between these brain regions. The results support other studies suggesting that crucial pathways involved in neurotransmission are distinct among the BNST, NAc, and dStr and provide insight into the potential use of pharmacological agents that may target region-specific signaling pathways. Furthermore, these studies provide a framework for future mouse-mouse comparisons of transcriptional profiles after behavioral/pharmacological manipulation.

Expression of a tumor-related gene network increases in the mammalian hypothalamus at the time of female puberty

Much has been learned in recent years about the central mechanisms controlling the initiation of mammalian puberty. It is now clear that this process requires the interactive participation of several genes. Using a combination of high throughput, molecular, and bioinformatics strategies, in combination with a system biology approach, we singled out from the hypothalamus of nonhuman primates and rats a group of related genes whose expression increases at the time of female puberty. Although these genes [henceforth termed tumor-related genes (TRGs)] have diverse cellular functions, they share the common feature of having been earlier identified as involved in tumor suppression/tumor formation. A prominent member of this group is KiSS1, a gene recently shown to be essential for the occurrence of puberty. Cis-regulatory analysis revealed the presence of a hierarchically arranged gene set containing five major hubs (CDP/CUTL1, MAF, p53, YY1, and USF2) controlling the network at the transcriptional level. In turn, these hubs are heavily connected to non-TRGs involved in the transcriptional regulation of the pubertal process. TRGs may be expressed in the mammalian hypothalamus as components of a regulatory gene network that facilitates and integrates cellular and cell-cell communication programs required for the acquisition of female reproductive competence.

Gene expression assays

The study of the profile of gene expression in a cell or tissue at a particular moment gives an insight into the plans of the cell for protein synthesis. Recent technological advances make it possible to analyze the expression of the entire genome in a single experiment. These "gene expression assays" complement or replace previous assays which measured the gene expression of only one gene, or a select group of genes. Within this chapter we outline the development of the gene expression assay and provide examples of the wide range of disciplines in which it is used. An overview of the current technologies is given, and includes an introduction to laser capture microdissection and linear amplification of RNA, both of which have extended the application of gene expression assays. Illustrative examples in the field of cancer and neuroscience highlight the scientific achievements. This technology has made in understanding the pathogenesis of diseases, including breast cancer, Huntington's disease, and schizophrenia. With recent advances including exon arrays to investigate alternative splicing, tiling arrays to investigate novel transcription start sites, and on-chip chromatin immunoprecipitation to investigate DNA-protein interactions, the future of gene expression assays is set to further our understanding of the complexities of gene expression.

Complementary techniques: RNA amplification for gene profiling analysis

The study of clinical samples is often limited by the amount of material available. DNA and RNA can be amplified from small specimens and, therefore, used for high-throughput analyses. While precise estimates of the level of DNA concentration in a given specimen is rarely studied (with the exception of relatively crude analyses of gene amplification or loss in cancer specimens), it is critical to know the proportional expression of various RNA transcripts since this proportion governs cell function by modulating the expression of various proteins. In addition, accurate estimates of relative RNA expression in biological conditions portray the reaction of cells to environmental stimuli shedding light on the characteristics of the microenvironment associated with particular physiologic or pathologic conditions. For this reason, the development of technologies for high fidelity messenger RNA amplification have been focused of extreme interest in the past decade with specific aim not only of increasing the abundance of RNA available to study but to accurately maintain the proportionality of expression of various RNA species among each other within a given specimen. This chapter will discuss various approaches to proportional RNA amplification focusing on amplification of the whole transcriptome (all transcripts in a given samples) rather than individual genes. These methods are suitable for high-throughput transcriptional profiling studies.

Detection of HHV-6B in post-mortem central nervous system tissue of a post-bone marrow transplant recipient: a multi-virus array analysis

Background

HHV-6 has been implicated in a number of neurological disorders. Recent evidence has suggested high incidence of HHV-6 infection in patients (46%) undergoing allogeneic bone marrow transplant (BMT).

Objective

To investigate whether HHV-6 plays a role in the development of fatal encephalopathy in an allogeneic post-BMT patient using an unbiased approach.

Results

Detection of HHV-6 viral DNA sequence and RNA expression were demonstrated in fresh frozen post-mortem autopsy material derived from the insular cortex using a multi-virus array platform. In addition, PCR analysis by real-time quantitative TaqMan demonstrated high viral burden in multiple brain regions tested. Sequencing analysis of PCR product confirmed the virus to be HHV-6 variant B.

Conclusions

Active infection as demonstrated by expression of viral RNA and high viral load in the CNS suggest a possible pathogenic role of HHV-6 in development neurologic complications post-BMT.

In vivo transcript profiling of Candida albicans identifies a gene essential for interepithelial dissemination

Candida albicans is the most common oral fungal pathogen of humans, but the mechanisms by which C. albicans invades and persists within mucosal epithelium are not clear. To understand oral pathogenesis, we characterized the cellular and molecular mechanisms of epithelial-fungus interactions using reconstituted human oral epithelium (RHE). We observed that hyphal formation facilitates epithelial invasion via both active (physical penetration) and passive (induced endocytosis) processes. Genome wide transcript profiling of C. albicans experimental RHE infection was compared with that from 11 patient samples with pseudomembranous candidiasis to identify genes associated with disease development in vivo. Expression profiles reflected the morphological switch and an adaptive response to neutral pH, non-glucose carbon sources and nitrosative stress. We identified several novel infection-associated genes with unknown function. One gene, upregulated in both RHE infection and patients, named EED1, was essential for maintenance of hyphal elongation. Mutants lacking EED1 showed transient cell elongation on epithelial tissue, which enabled only superficial invasion of epithelial cells. Once inside an epithelial cell, Deltaeed1 cells could proliferate as yeasts or pseudohyphae but remained trapped intracellularly. Our results suggest that the adaptive response and morphology of C. albicans play specific roles for host-fungal interactions during mucosal infections.

An attempt for combining microarray data sets by adjusting gene expressions

PURPOSE

The diverse experimental environments in microarray technology, such as the different platforms or different RNA sources, can cause biases in the analysis of multiple microarrays. These systematic effects present a substantial obstacle for the analysis of microarray data, and the resulting information may be inconsistent and unreliable. Therefore, we introduced a simple integration method for combining microarray data sets that are derived from different experimental conditions, and we expected that more reliable information can be detected from the combined data set rather than from the separated data sets.

MATERIALS AND METHODS

This method is based on the distributions of the gene expression ratios among the different microarray data sets and it transforms, gene by gene, the gene expression ratios into the form of the reference data set. The efficiency of the proposed integration method was evaluated using two microarray data sets, which were derived from different RNA sources, and a newly defined measure, the mixture score.

RESULTS

The proposed integration method intermixed the two data sets that were obtained from different RNA sources, which in turn reduced the experimental bias between the two data sets, and the mixture score increased by 24.2%. A data set combined by the proposed method preserved the inter-group relationship of the separated data sets.

CONCLUSION

The proposed method worked well in adjusting systematic biases, including the source effect. The ability to use an effectively integrated microarray data set yields more reliable results due to the larger sample size and this also decreases the chance of false negatives.

In vivo and ex vivo comparative transcriptional profiling of invasive and non-invasive Candida albicans isolates identifies genes associated with tissue invasion

The human pathogenic fungus Candida albicans can cause a wide range of infections and invade multiple organs. To identify C. albicans genes that are expressed during invasion of the liver, we used genome-wide transcriptional profiling in vivo and ex vivo. By analysing the different phases of intraperitoneal infection from attachment to tissue penetration in a time-course experiment and by comparing the profiles of an invasive with those of a non-invasive strain, we identified genes and transcriptional pattern which are associated with the invasion process. This includes genes involved in metabolism, stress, and nutrient uptake, as well as transcriptional programmes regulating morphology and environmental sensing. One of the genes identified as associated with liver invasion was DFG16, a gene crucial for pH-dependent hyphal formation, correct pH sensing, invasion at physiological pH and systemic infection.

Novel and simple transformation algorithm for combining microarray data sets

Background

With microarray technology, variability in experimental environments such as RNA sources, microarray production, or the use of different platforms, can cause bias. Such systematic differences present a substantial obstacle to the analysis of microarray data, resulting in inconsistent and unreliable information. Therefore, one of the most pressing challenges in the field of microarray technology is how to integrate results from different microarray experiments or combine data sets prior to the specific analysis.

Results

Two microarray data sets based on a 17k cDNA microarray system were used, consisting of 82 normal colon mucosa and 72 colorectal cancer tissues. Each data set was prepared from either total RNA or amplified mRNA, and the difference of RNA source between these two data sets was detected by ANOVA (Analysis of variance) model. A simple integration method was introduced which was based on the distributions of gene expression ratios among different microarray data sets. The method transformed gene expression ratios into the form of a reference data set on a gene by gene basis. Hierarchical clustering analysis, density and box plots, and mixture scores with correlation coefficients revealed that the two data sets were well intermingled, indicating that the proposed method minimized the experimental bias. In addition, any RNA source effect was not detected by the proposed transformation method. In the mixed data set, two previously identified subgroups of normal and tumor were well separated, and the efficiency of integration was more prominent in tumor groups than normal groups. The transformation method was slightly more effective when a data set with strong homogeneity in the same experimental group was used as a reference data set.

Conclusion

Proposed method is simple but useful to combine several data sets from different experimental conditions. With this method, biologically useful information can be detectable by applying various analytic methods to the combined data set with increased sample size.

Technological platforms for microarray gene expression profiling

By using gene microarray technology, scientists can determine in a single experiment, the expression levels of thousands of genes within a given sample. DNA microarray technology is evolving rapidly and there are now numerous high-density platforms available which differ in terms of probe content, design, deposition technology, labeling and hybridization protocols. However, two major platforms for high-density microarray manufacture are in common use. The first utilizes robotic deposition or "spotting" of DNA molecules, while the second uses short oligonucleotides synthesized in situ.