Open systems: panoramic views of gene expression

Identification of Five Genes as a Potential Biomarker for Predicting Progress and Prognosis in Adrenocortical Carcinoma

Background: Adrenocortical carcinoma (ACC) is a limited endocrine fatality with a minor diagnosis and rare remedial options. The progressive and predictive meaning of message RNA (mRNA) expression oddity in ACC has been studied extensively in recent years. However, differences in measurement platforms and lab protocols as well as small sample sizes can render gene expression levels incomparable.

Methods: An extensive study of GEO datasets was conducted to define potential mRNA biomarkers for ACC. The study compared the mRNA expression profiles of ACC tissues and neighboring noncancerous adrenal tissues in the pair. The study covered a sum of 165 tumors and 36 benign control samples. Hub genes were identified through a protein-protein interaction (PPI) network and Robust Rank Aggregation method. Then the Cancer Genome Atlas (TCGA) and Oncomine database were used to perform the validation of hub genes. 4 ACC tissues and 4 normal tissues were collected and then Polymerase Chain Reaction (PCR), Western-blot and immunofluorescence were conducted to validate the expression of five hub genes.

Results: We identified five statistically significant genes (TOP2A, NDC80, CEP55, CDKN3, CDK1) corrected with clinical features. The expression of five hub genes in TCGA and Oncomine database were significantly overexpressed in ACC compared with normal ones. Among all the TCGA ACC cases, the strong expression of TOP2A (logrank p=1.4e-04, HR=4.7), NDC80 (logrank p=8.8e-05, HR=4.9), CEP55 (logrank p=5.2e-07, HR=8.6), CDKN3 (log rank p=2.3e-06, HR=7.6) and CDK1 (logrank p=7e-08, HR=11) were correlated with low comprehensive survival, disease free survival (logrank p < 0.001), pathology stage and pathology T stage (FDR < 0.001). PCR results showed that the transcriptional levels of these five genes were significantly higher in ACC tissues than in normal tissues. The western blotting results also showed that the translational level of TOP2A was significantly higher in tumor tissues than in normal tissues. The results of immunofluorescence showed that TOP2A was abundantly observed in the adrenal cortical cell membrane and nucleus and its expression in ACC tissues was significantly higher than that in normal tissues.

Conclusions: The distinguished five genes may be utilized to form a board of progressive and predictive biomarkers for ACC for clinical purpose.

Down but Not Out: The Role of MicroRNAs in Hibernating Bats

MicroRNAs (miRNAs) regulate many physiological processes through post-transcriptional control of gene expression and are a major part of the small noncoding RNAs (snRNA). As hibernators can survive at low body temperatures (Tb) for many months without suffering tissue damage, understanding the mechanisms that enable them to do so are of medical interest. Because the brain integrates peripheral physiology and white adipose tissue (WAT) is the primary energy source during hibernation, we hypothesized that both of these organs play a crucial role in hibernation, and thus, their activity would be relatively increased during hibernation. We carried out the first genomic analysis of small RNAs, specifically miRNAs, in the brain and WAT of a hibernating bat (Myotis ricketti) by comparing deeply torpid with euthermic individual bats using high-throughput sequencing (Solexa) and qPCR validation of expression levels. A total of 196 miRNAs (including 77 novel bat-specific miRNAs) were identified, and of these, 49 miRNAs showed significant differences in expression during hibernation, including 33 in the brain and 25 in WAT (P≤0.01 &│logFC│≥1). Stem-loop qPCR confirmed the miRNA expression patterns identified by Solexa sequencing. Moreover, 31 miRNAs showed tissue- or state-specific expression, and six miRNAs with counts >100 were specifically expressed in the brain. Putative target gene prediction combined with KEGG pathway and GO annotation showed that many essential processes of both organs are significantly correlated with differentially expressed miRNAs during bat hibernation. This is especially evident with down-regulated miRNAs, indicating that many physiological pathways are altered during hibernation. Thus, our novel findings of miRNAs and Interspersed Elements in a hibernating bat suggest that brain and WAT are active with respect to the miRNA expression activity during hibernation.

Transcriptome profile of human neuroblastoma cells in the hypomagnetic field

Research has shown that the hypomagnetic field (HMF) can affect embryo development, cell proliferation, learning and memory, and in vitro tubulin assembly. In the present study, we aimed to elucidate the molecular mechanism by which the HMF exerts its effect, by comparing the transcriptome profiles of human neuroblastoma cells exposed to either the HMF or the geomagnetic field. A total of 2464 differentially expressed genes (DEGs) were identified, 216 of which were up-regulated and 2248 of which were down-regulated after exposure to the HMF. These DEGs were found to be significantly clustered into several key processes, namely macromolecule localization, protein transport, RNA processing, and brain function. Seventeen DEGs were verified by real-time quantitative PCR, and the expression levels of nine of these DEGs were measured every 6 h. Most notably, MAPK1 and CRY2, showed significant up- and down-regulation, respectively, during the first 6 h of HMF exposure, which suggests involvement of the MAPK pathway and cryptochrome in the early bio-HMF response. Our results provide insights into the molecular mechanisms underlying the observed biological effects of the HMF.

Transcriptome analysis of HeLa cells response to Brucella melitensis infection: a molecular approach to understand the role of the mucosal epithelium in the onset of the Brucella pathogenesis

Brucella spp. infect hosts primarily by adhering and penetrating mucosal surfaces, however the initial molecular phenomena of this host:pathogen interaction remain poorly understood. We hypothesized that characterizing the epithelial-like human HeLa cell line molecular response to wild type Brucella melitensis infection would help to understand the role of the mucosal epithelium at the onset of the Brucella pathogenesis. RNA samples from B. melitensis-infected HeLa cells were taken at 4 and 12 h of infection and hybridized in a cDNA microarray. The analysis using a dynamic Bayesian network modeling approach (DBN) identified several pathways, biological processes, cellular components and molecular functions altered due to infection at 4 h p.i., but almost none at 12 h p.i. The in silico modeling results were experimentally tested by knocking down the expression of MAPK1 by siRNA technology. MAPK1-siRNA transfected cell cultures decreased the internalization and impaired the intracellular replication of the pathogen in HeLa cells after 4 h p.i. DBN analysis provides important insights into the role of the epithelial cells response to Brucella infection and guide research to novel mechanisms identification.

Use of mutagenesis, genetic mapping and next generation transcriptomics to investigate insecticide resistance mechanisms

Insecticide resistance is a worldwide problem with major impact on agriculture and human health. Understanding the underlying molecular mechanisms is crucial for the management of the phenomenon; however, this information often comes late with respect to the implementation of efficient counter-measures, particularly in the case of metabolism-based resistance mechanisms. We employed a genome-wide insertional mutagenesis screen to Drosophila melanogaster, using a Minos-based construct, and retrieved a line (MiT[w⁻]3R2) resistant to the neonicotinoid insecticide Imidacloprid. Biochemical and bioassay data indicated that resistance was due to increased P450 detoxification. Deep sequencing transcriptomic analysis revealed substantial over- and under-representation of 357 transcripts in the resistant line, including statistically significant changes in mixed function oxidases, peptidases and cuticular proteins. Three P450 genes (Cyp4p2, Cyp6a2 and Cyp6g1) located on the 2R chromosome, are highly up-regulated in mutant flies compared to susceptible Drosophila. One of them (Cyp6g1) has been already described as a major factor for Imidacloprid resistance, which validated the approach. Elevated expression of the Cyp4p2 was not previously documented in Drosophila lines resistant to neonicotinoids. In silico analysis using the Drosophila reference genome failed to detect transcription binding factors or microRNAs associated with the over-expressed Cyp genes. The resistant line did not contain a Minos insertion in its chromosomes, suggesting a hit-and-run event, i.e. an insertion of the transposable element, followed by an excision which caused the mutation. Genetic mapping placed the resistance locus to the right arm of the second chromosome, within a ∼1 Mb region, where the highly up-regulated Cyp6g1 gene is located. The nature of the unknown mutation that causes resistance is discussed on the basis of these results.

Identification of stress-responsive genes in plants using suppression subtraction hybridization: ozone stress as an example

Among the open-ended techniques for identifying differentially expressed genes in response to stress, the PCR-based suppression subtraction hybridization (SSH) is widely used. The popularity of this technique stems from the ease of conducting this procedure in any laboratory set up for basic molecular biology research. Further, the availability of a comprehensive kit for conducting suppression subtractions from BD Biosciences has made this technique easy to adapt and adopt to any biological system. In this chapter we describe in detail the SSH procedure and explain the subtle changes that have been incorporated to make this technique adaptable for identifying stress-responsive genes in plants.

Transcriptome response to pollutants and insecticides in the dengue vector Aedes aegypti using next-generation sequencing technology

BACKGROUND

The control of mosquitoes transmitting infectious diseases relies mainly on the use of chemical insecticides. However, mosquito control programs are now threatened by the emergence of insecticide resistance. Hitherto, most research efforts have been focused on elucidating the molecular basis of inherited resistance. Less attention has been paid to the short-term response of mosquitoes to insecticides and pollutants which could have a significant impact on insecticide efficacy. Here, a combination of LongSAGE and Solexa sequencing was used to perform a deep transcriptome analysis of larvae of the dengue vector Aedes aegypti exposed for 48 h to sub-lethal doses of three chemical insecticides and three anthropogenic pollutants.

RESULTS

Thirty millions 20 bp cDNA tags were sequenced, mapped to the mosquito genome and clustered, representing 6850 known genes and 4868 additional clusters not located within predicted genes. Mosquitoes exposed to insecticides or anthropogenic pollutants showed considerable modifications of their transcriptome. Genes encoding cuticular proteins, transporters, and enzymes involved in the mitochondrial respiratory chain and detoxification processes were particularly affected. Genes and molecular mechanisms potentially involved in xenobiotic response and insecticide tolerance were identified.

CONCLUSIONS

The method used in the present study appears as a powerful approach for investigating fine transcriptome variations in genome-sequenced organisms and can provide useful informations for the detection of novel transcripts. At the biological level, despite low concentrations and no apparent phenotypic effects, the significant impact of these xenobiotics on mosquito transcriptomes raise important questions about the 'hidden impact' of anthropogenic pollutants on ecosystems and consequences on vector control.

Bioprospecting metagenomes: glycosyl hydrolases for converting biomass

Throughout immeasurable time, microorganisms evolved and accumulated remarkable physiological and functional heterogeneity, and now constitute the major reserve for genetic diversity on earth. Using metagenomics, namely genetic material recovered directly from environmental samples, this biogenetic diversification can be accessed without the need to cultivate cells. Accordingly, microbial communities and their metagenomes, isolated from biotopes with high turnover rates of recalcitrant biomass, such as lignocellulosic plant cell walls, have become a major resource for bioprospecting; furthermore, this material is a major asset in the search for new biocatalytics (enzymes) for various industrial processes, including the production of biofuels from plant feedstocks. However, despite the contributions from metagenomics technologies consequent upon the discovery of novel enzymes, this relatively new enterprise requires major improvements. In this review, we compare function-based metagenome screening and sequence-based metagenome data mining, discussing the advantages and limitations of both methods. We also describe the unusual enzymes discovered via metagenomics approaches, and discuss the future prospects for metagenome technologies.

A combination of LongSAGE with Solexa sequencing is well suited to explore the depth and the complexity of transcriptome

Background

"Open" transcriptome analysis methods allow to study gene expression without a priori knowledge of the transcript sequences. As of now, SAGE (Serial Analysis of Gene Expression), LongSAGE and MPSS (Massively Parallel Signature Sequencing) are the mostly used methods for "open" transcriptome analysis. Both LongSAGE and MPSS rely on the isolation of 21 pb tag sequences from each transcript. In contrast to LongSAGE, the high throughput sequencing method used in MPSS enables the rapid sequencing of very large libraries containing several millions of tags, allowing deep transcriptome analysis. However, a bias in the complexity of the transcriptome representation obtained by MPSS was recently uncovered.

Results

In order to make a deep analysis of mouse hypothalamus transcriptome avoiding the limitation introduced by MPSS, we combined LongSAGE with the Solexa sequencing technology and obtained a library of more than 11 millions of tags. We then compared it to a LongSAGE library of mouse hypothalamus sequenced with the Sanger method.

Conclusion

We found that Solexa sequencing technology combined with LongSAGE is perfectly suited for deep transcriptome analysis. In contrast to MPSS, it gives a complex representation of transcriptome as reliable as a LongSAGE library sequenced by the Sanger method.

Metagenomics in animal gastrointestinal ecosystem: a microbiological and biotechnological perspective

Metagenomics- the application of the genomics technologies to nonculturable microbial communities, is coming of age. These approaches can be used for the screening and selection of nonculturable rumen microbiota for assessing their role in gastrointestinal (GI) nutrition, plant material fermentation and the health of the host. The technologies designed to access this wealth of genetic information through environmental nucleic acid extraction have provided a means of overcoming the limitations of culture-dependent microbial genetic exploitation. The molecular procedures and techniques will result in reliable insights into the GI microbial structure and activity of the livestock gut microbes in relation to functional interactions, temporal and spatial relationships among different microbial consortia and dietary ingredients. Future developments and applications of these methods promise to provide the first opportunity to link distribution and identity of rumen microbes in their natural habitats with their genetic potential and in situ activities.

Isolation of mycoparasitic-related transcripts by SSH during interaction of the mycoparasite Stachybotrys elegans with its host Rhizoctonia solani

Mycoparasitism by antagonistic fungi involves changes in the biochemistry and physiology of both partners. Analysis of genes that are expressed during mycoparasite-host interaction represents a powerful strategy to obtain insight into the molecular events underlying these changes. The aim of this study is to identify genes whose expression is upregulated when the mycoparasite Stachybotrys elegans is in direct confrontation with its host Rhizoctonia solani. Suppression subtractive hybridization (SSH) was used to create a subtracted cDNA library, and differential screening was applied to identify the over-expressed transcripts. We report the analysis of 2,166 clones, among which 47% were upregulated during mycoparasitism. Two hundred and sixty-one clones were sequenced that corresponded to 94 unique genes. Forty-four of these were identified as novel genes, while the remainder showed similarity to a broad diversity of genes with putative functions related to toxin production, pathogenicity, and metabolism. As a result of mycoparasitism, 15 genes belonged to R. solani among which 9 genes were assigned putative functions. Quantitative RT-PCR was used to examine the upregulation of 12 genes during the course of mycoparasitism. Seven genes showed significant upregulation at least at one-time point during interaction of the mycoparasite with its host. This study describes a first step toward knowledge of S. elegans genome. The results present the useful application of EST analysis on S. elegans and provide preliminary indication of gene expression putatively involved in mycoparasitism.

Deep analysis of cellular transcriptomes - LongSAGE versus classic MPSS

BACKGROUND

Deep transcriptome analysis will underpin a large fraction of post-genomic biology. 'Closed' technologies, such as microarray analysis, only detect the set of transcripts chosen for analysis, whereas 'open' e.g. tag-based technologies are capable of identifying all possible transcripts, including those that were previously uncharacterized. Although new technologies are now emerging, at present the major resources for open-type analysis are the many publicly available SAGE (serial analysis of gene expression) and MPSS (massively parallel signature sequencing) libraries. These technologies have never been compared for their utility in the context of deep transcriptome mining.

RESULTS

We used a single LongSAGE library of 503,431 tags and a "classic" MPSS library of 1,744,173 tags, both prepared from the same T cell-derived RNA sample, to compare the ability of each method to probe, at considerable depth, a human cellular transcriptome. We show that even though LongSAGE is more error-prone than MPSS, our LongSAGE library nevertheless generated 6.3-fold more genome-matching (and therefore likely error-free) tags than the MPSS library. An analysis of a set of 8,132 known genes detectable by both methods, and for which there is no ambiguity about tag matching, shows that MPSS detects only half (54%) the number of transcripts identified by SAGE (3,617 versus 1,955). Analysis of two additional MPSS libraries shows that each library samples a different subset of transcripts, and that in combination the three MPSS libraries (4,274,992 tags in total) still only detect 73% of the genes identified in our test set using SAGE. The fraction of transcripts detected by MPSS is likely to be even lower for uncharacterized transcripts, which tend to be more weakly expressed. The source of the loss of complexity in MPSS libraries compared to SAGE is unclear, but its effects become more severe with each sequencing cycle (i.e. as MPSS tag length increases).

CONCLUSION

We show that MPSS libraries are significantly less complex than much smaller SAGE libraries, revealing a serious bias in the generation of MPSS data unlikely to have been circumvented by later technological improvements. Our results emphasize the need for the rigorous testing of new expression profiling technologies.

Oxidative stress response of European flounder (Platichthys flesus) to cadmium determined by a custom cDNA microarray

The monitoring of the impact of chemical pollutants upon marine ecosystems commonly employs a multi-biomarker approach. Functional genomics, using cDNA microarrays, allows for a comprehensive view of how an organism is responding to an exposure, with respect to changes in gene expression. Differentially expressed mRNAs were first isolated from livers of European flounder by means of suppressive, subtractive hybridisation. A clone set containing a total of 284 different potentially differentially expressed mRNAs was produced, of which 84 were tentatively identified. These were combined with previously cloned known stress genes isolated by degenerate PCR to produce a custom 500-clone microarray platform with each clone arrayed to four spots. Subsequent array experiments using cadmium-treated flounder detected up-regulation of 27 transcripts, including Cu/Zn superoxide dismutase, thioredoxin, a peroxiredoxin and a glutathione-S-transferase, reflecting oxidative stress in exposed flounder, while CYP1A expression was down-regulated. These changes were confirmed by real-time PCR. The array experiment highlighted a number of candidate genes for further analysis as potential novel biomarkers of cadmium exposure and demonstrated the applicability of the custom microarray approach in the study of the effects of toxicants.

Identification of Differentially Expressed Genes by cDNA-Amplified Fragment Length Polymorphism in the Biocontrol Agent Pichia anomala (Strain Kh5)

ABSTRACT cDNA-amplified fragment length polymorphism (cDNA-AFLP) analysis was used to identify genes potentially involved in biological control, by strain Kh5 (Pichia anomala), of Botrytis cinerea, an important post-harvest pathogen on apples. Strain Kh5 was grown in yeast nitrogen base (YNB) plus glucose (G medium) or YNB plus cell walls of B. cinerea (B medium). Thirty-five primer pairs were used in AFLP amplifications, resulting in a total of more than 2,450 bands derived from the mRNA of strain Kh5 grown in B medium. Eighty-six bands (3.5%) corresponded to genes upregulated in B medium compared with G medium. Of these 86 bands, 28 were selected, cloned, sequenced, and subjected to real-time reverse transcription-polymerase chain reaction (RT-PCR) to confirm their differential expression. An appropriate housekeeping gene, G2, was selected and used to normalize the results of RT-PCR. Eleven genes presented an increased gene expression in the presence of B. cinerea cell walls (expression >1). Statistical analysis showed a significant increase for 5 of these 11 genes. The overexpressed genes show homologies to yeast genes with various functions, including beta-glucosidase, transmembrane transport, citrate synthase, and external amino acid sensing and transport. Some of these functions could be related to cell wall metabolism and potentially involved in mycoparasitic properties.

Metagenomic gene discovery: past, present and future

It is now widely accepted that the application of standard microbiological methods for the recovery of microorganisms from the environment has had limited success in providing access to the true extent of microbial biodiversity. It follows that much of the extant microbial genetic diversity (collectively termed the metagenome) remains unexploited, an issue of considerable relevance to a wider understanding of microbial communities and of considerable importance to the biotechnology industry. The recent development of technologies designed to access this wealth of genetic information through environmental nucleic acid extraction has provided a means of avoiding the limitations of culture-dependent genetic exploitation.

HIV-1-derived lentiviral vectors and fetal route of administration on transgene biodistribution and expression in rhesus monkeys

The gene transfer efficiency of lentiviral vectors pseudotyped with vesicular stomatitis virus-glycoprotein (VSV-G) driven by the MND or CMV promoters and expressing the enhanced green fluorescent protein (EGFP) was investigated in fetal rhesus monkeys (Macaca mulatta) (N=21). Fetuses (50+/-10 days gestation; term 165+/-10 days) were injected under ultrasound guidance using an intraperitoneal (i.p.) or intrahepatic (i.h.) approach with a range of 1 x 10(7)-2 x 10(8) infectious particles/fetus. Analysis of transgene biodistribution and expression was performed in multiple tissues at 3-7 months postgene delivery using quantitative techniques. Overall, results indicated the following: (1) i.p. gene transfer at 40 days gestation resulted in a more diffuse distribution of the vector compared to administration at 60 days gestation; (2) vector biodistribution was similar after administration by the i.p. or i.h. routes; and (3) gene expression analysis in transduced tissues showed the presence of mRNA transcripts that correlated with the level of gene transfer. These studies suggest that fetal gene transfer using the i.p. and i.h. routes results in prolonged transduction and expression of the transgene in multiple tissues.

Identification of two novel human genes, DIPLA1 and DIPAS, expressed in placenta tissue

Here we report the identification and expression analysis of two novel human genes--DIPLA1 (Differentially expressed in placenta 1) and DIPAS (DIPLA1 Antisense). These genes are located at chromosomal region 9q33.1, in opposite orientations, and are flanked by the pregnancy-associated plasma protein-A (PAPP-A) and astrotactin 2 (ASTN2) genes. The mRNA sequences of both genes contain several upstream AUGs (uAUG) and various potential open reading frames (ORFs). DIPLA1 mRNA is 1.8 kb long and contains a 285 nt ORF coding for a polypeptide designated as replicative senescence up-regulated (RSU) protein. Antisense DIPAS mRNA is 2.7 kb long and contains a 309 nt ORF coding for a protein with partial similitude to the gamma isoform variant of the human Ca(2+)/calmodulin (CaM)-dependent protein kinase II. Both genes are conserved in placental-species and are presumably transcribed from initiator (Inr) promoter elements located at opposite strands. In 20 human normal tissues tested, DIPLA1 mRNA expression was placenta-specific, whereas DIPAS mRNA expression was higher in placenta, brain, kidney and testis. In addition, DIPAS mRNA hybridizes with the 3'UTR region from PAPP-A mRNA, which spans over 4 kb more than previously reported, forming a potential sense-antisense double stranded RNA (dsRNA) duplex. Our results are of interest for placenta gene expression regulation and for the identification of novel genes in the human genome.

Isolation of candidate genes for apomixis in Poa pratensis L

The essential feature of apomixis is that an embryo is formed autonomously by parthenogenesis from an unreduced egg of an embryo sac generated through apomeiosis. The genetic constitution of the offspring is, therefore, usually identical to the maternal parent, a trait of great interest to plant breeders. If apomixis were well understood and harnessed, it could be exploited to indefinitely propagate superior hybrids or specific genotypes bearing complex gene sets. A fundamental contribution to the understanding of the genetic control of the apomictic pathway could be provided by a deep knowledge of molecular mechanisms that regulate the reproductive events. In Poa pratensis the cDNA-AFLP method of mRNA profiling allowed us to visualize a total of 2248 transcript-derived fragments and to isolate 179 sequences that differed qualitatively or quantitatively between apomictic and sexual genotypes at the time of flowering when the primary stages of apomixis occur. Three ESTs were chosen for further molecular characterization because of their cDNA-AFLP expression pattern and BLAST information retrieval. The full-lengths of the newly isolated genes were recovered by RACE and their temporal expression patterns were assessed by RT-PCR. Their putative role in cell signaling transduction cascades and trafficking events required during sporogenesis, gametogenesis and embryogenesis in plants is reported and discussed.

Transcriptomic analysis in the leech Theromyzon tessulatum: involvement of cystatin B in innate immunity

At the present time, there is little information on mechanisms of innate immunity in invertebrate groups other than insects, especially annelids. In the present study, we have performed a transcriptomic study of the immune response in the leech Theromyzon tessulatum after bacterial challenge, by a combination of differential display RT (reverse transcriptase)-PCR and cDNA microarrays. The results show relevant modulations concerning several known and unknown genes. Indeed, threonine deaminase, malate dehydrogenase, cystatin B, polyadenylate-binding protein and alpha-tubulin-like genes are up-regulated after immunostimulation. We focused on cystatin B (stefin B), which is an inhibitor of cysteine proteinases involved in the vertebrate immune response. We have cloned the full-length cDNA and named the T. tessulatum gene as Tt-cysb. Main structural features of cystatins were identified in the derived amino acid sequence of Tt-cysb cDNA; namely, a glycine residue in the N-terminus and a consensus sequence of Gln-Xaa-Val-Xaa-Gly (QXVXG) corresponding to the catalytic site. Moreover, Tt-cysb is the first cystatin B gene characterized in invertebrates. We have determined by in situ hybridization and immunocytochemistry that Tt-cysb is only expressed in large coelomic cells. In addition, this analysis confirmed that Tt-cysb is up-regulated after bacterial challenge, and that increased expression occurs only in coelomic cells. These data demonstrate that the innate immune response in the leech involves a cysteine proteinase inhibitor that is not found in ecdysozoan models, such as Drosophila melanogaster or Caenorhabditis elegans, and so underlines the great need for information about innate immunity mechanisms in different invertebrate groups.

Molecular biomarkers in drug development

Arguably, the most immediately promising reverberation of the genomics era has been the application of biomarkers to drug development. The promise of applying biomarkers to early drug development is that they might aid in preclinical and early clinical decisions such as dose ranging, definition of treatment regimen, or even a preview of efficacy. Later in the clinic, biomarkers could be used to facilitate patient stratification, selection and the description of surrogate endpoints. Information derived from biomarkers should result in a better understanding of preclinical and clinical data, which ultimately benefits patients and drug developers. If the promise of biomarkers is realized, they will become a routine component of drug development and companions to newly discovered therapies.

From differential gene expression to differential gene function and back

Following the decoding of several plant and animal genomes, the identification of all corresponding transcripts and proteins and understanding how their expression corresponds to physiological and pathological states is the obvious next step. Nucleic acid quantification methods have become increasingly high-throughput and relatively low-cost, and moving ahead, combinations of technologies monitoring differential gene expression and those defining differential cellular function will yield maximum benefit in furthering biology and for drug target identification and validation.:

Differential display RT-PCR analysis of human choriocarcinoma cell lines and normal term trophoblast cells: identification of new genes expressed in placenta

In this study, we performed the differential display technique to identify genes specifically expressed in human choriocarcinoma cell lines (JEG-3, JAR and BeWo) and normal placental term cells. Few differences were found among the expression profiles of the three choriocarcinoma cell lines and most of the differentially expressed genes were detected in normal term placenta. A total of 36 cDNA fragments were isolated and analysed. Of these, 19 sequences corresponded to regions in the human genome coding for potential novel genes. We confirmed by RT-PCR, the placental mRNA expression of three selected new human genes, on chromosomes 16q12, 9q32 and 6q22. The other 17 sequences showed high similarity to known human genes (like PSG3, FN1, PAI-2). Interestingly, the functions of five known proteins (from genes IK, TRA-1, HERPUD1, UBA-2, and TRAP240) have not yet been well characterized in placenta tissue. In addition, new alternative spliced mRNAs were detected for IK, TRAP240 and PLAC3 genes. The differential expression of the PAI-2 gene among the choriocarcinoma cell lines was also confirmed. The genes identified in this analysis will be of interest for future studies regarding both a better understanding of the biology of the trophoblast cell and the formation of placental tumors.

Modified serial analysis of gene expression method for construction of gene expression profiles of microbial eukaryotic species

Serial analysis of gene expression (SAGE) is a powerful approach for the identification of differentially expressed genes, providing comprehensive and quantitative gene expression profiles in the form of short tag sequences. Each tag represents a unique transcript, and the relative frequencies of tags in the SAGE library are equal to the relative proportions of the transcripts they represent. One of the major obstacles in the preparation of SAGE libraries from microorganisms is the requirement for large amounts of starting material (i.e., mRNA). Here, we present a novel approach for the construction of SAGE libraries from small quantities of total RNA by using Y linkers to selectively amplify 3' cDNA fragments. To validate this method, we constructed comprehensive gene expression profiles of the toxic dinoflagellate Pfiesteria shumwayae. SAGE libraries were constructed from an actively toxic fish-fed culture of P. shumwayae and from a recently toxic alga-fed culture. P. shumwayae-specific gene transcripts were identified by comparison of tag sequences in the two libraries. Representative tags with frequencies ranging from 0.026 to 3.3% of the total number of tags in the libraries were chosen for further analysis. Expression of each transcript was confirmed in separate control cultures of toxic P. shumwayae. The modified SAGE method described here produces gene expression profiles that appear to be both comprehensive and quantitative, and it is directly applicable to the study of gene expression in other environmentally relevant microbial species.

Gene expression analysis of the critical period in the visual cortex

The development of the primary visual cortex in animals possessing binocular vision is a classical paradigm for the study of activity-dependent neuronal plasticity. To elucidate the genetic determinants of this period of substantial plasticity, we conducted an unbiased and comprehensive transcript profiling analysis with differential display and DNA array techniques. We characterized the transcripts that change significantly between the critical and postcritical periods in the rat binocular visual cortex. We determined if these changes are specific for the visual cortex by simultaneously profiling the hippocampus and examined the impact of sensory experience on the accumulation of the identified transcripts. Our results uncover visual cortex-specific and unspecific transcription programs. Transcripts for protein kinases and phosphatases are particularly regulated. The identified transcripts support the notion that the critical period provides a permissive state for plasticity.

Quantitative comparison of cDNA-AFLP, microarrays, and GeneChip expression data in Saccharomyces cerevisiae

cDNA-AFLP is a genome-wide expression analysis technology that does not require any prior knowledge of gene sequences. This PCR-based technique combines a high sensitivity with a high specificity, allowing detection of rarely expressed genes and distinguishing between homologous genes. In this report, we validated quantitative expression data of 110 cDNA-AFLP fragments in yeast with DNA microarrays and GeneChip data. The best correlation was found between cDNA-AFLP and GeneChip data. The cDNA-AFLP data revealed a low number of inconsistent profiles that could be explained by gel artifact, overexposure, or mismatch amplification. In addition, 18 cDNA-AFLP fragments displayed homology to genomic yeast DNA, but could not be linked unambiguously to any known ORF. These fragments were most probably derived from 5' or 3' noncoding sequences or might represent previously unidentified ORFs. Genes liable to cross hybridization showed identical results in cDNA-AFLP and GeneChip analysis. Three genes, which were readily detected with cDNA-AFLP, showed no significant expression in GeneChip experiments. We show that cDNA-AFLP is a very good alternative to microarrays and since no preexisting biological or sequence information is required, it is applicable to any species.

Isolation of differentially expressed genes from contaminant exposed European flounder by suppressive, subtractive hybridisation

Commercially farmed European flounder (Platichthys flesus) were treated with benzo[a]pyrene or cadmium under controlled laboratory conditions. In addition, animals were collected from the polluted Tyne and relatively clean Alde estuaries in the UK. Potential differentially expressed genes were isolated from liver tissue by means of suppressive, subtractive hybridisation. A total of 284 different mRNAs were identified as having the potential for being differentially expressed, of which 84 could be tentatively identified by homology with database sequences. These clones will form part of a flounder cDNA array for validation and future applications, with full identification being pursued for those expressed sequence tags for which differential expression can be confirmed by array technology.

The future of genetic research on appetitive behavior

The risk for developing obesity has a significant genetic component. Several quantitative trait loci and candidate genes have been identified using current methodological approaches however the information gained thus far is insufficient to adequately explain the genetics underlying human obesity. The completion of a draft of the human genome sequence, the potential benefit of single nucleotide polymorphisms association studies for identifying risk conferring alleles, and developing functional genomics technologies promise to accelerate obesity gene discovery. These advances, used with current evaluative tools (murine molecular genetic techniques), may increase our understanding of human obesity, and ultimately provide better approaches to diagnosis and treatment.

Common genomic response in different mouse models of beta-adrenergic-induced cardiomyopathy

BACKGROUND

Although beta-adrenergic receptor (AR) blockade therapy is beneficial in the treatment of heart failure, little is known regarding the transcriptional mechanisms underlying this salutary action.

METHODS AND RESULTS

In the present study, we screened mice overexpressing Gsalpha, beta1AR, beta2AR, or protein kinase A to test if a common genomic pathway exists in different models with enhanced beta-adrenergic signaling. In mice overexpressing Gsalpha, differentially expressed genes were identified by mRNA profiling. In addition to well-known markers of cardiac hypertrophy (atrial natriuretic factor, CARP, and beta-myosin heavy chain), uncoupling protein 2 (UCP2), a protein involved in the control of mitochondrial membrane potential, and four-and-a-half LIM domain protein-1 (FHL1), a member of the LIM protein family, were predicted to be upregulated. Upregulation of these genes was confirmed by quantitative reverse transcriptase-polymerase chain reaction at all time points tested during the development of cardiomyopathy in mice overexpressing Gsalpha. In mice overexpressing beta1AR, beta2AR, or protein kinase A, increased UCP2 and FHL1 expression was also observed at the onset of cardiomyopathy. BetaAR blockade treatment reversed the cardiomyopathy and suppressed the increased expression of UCP2 and FHL1 in mice overexpressing Gsalpha.

CONCLUSIONS

UCP2 and FHL1 are important candidate genes that correlate with the development of betaAR-induced cardiomyopathy in different mouse models with enhanced betaAR signaling. In addition to preserving cardiac function, betaAR blockade treatment also prevents the genomic regulation that correlates with the onset of heart failure.

A sensitive transcriptome analysis method that can detect unknown transcripts

We have developed an AFLP-based gene expression profiling method called 'high coverage expression profiling' (HiCEP) analysis. By making improvements to the selective PCR technique we have reduced the rate of false positive peaks to approximately 4% and consequently the number of peaks, including overlapping peaks, has been markedly decreased. As a result we can determine the relationship between peaks and original transcripts unequivocally. This will make it practical to prepare a database of all peaks, allowing gene assignment without having to isolate individual peaks. This precise selection also enables us to easily clone peaks of interest and predict the corresponding gene for each peak in some species. The procedure is highly reproducible and sensitive enough to detect even a 1.2-fold difference in gene expression. Most importantly, the low false positive rate enables us to analyze gene expression with wide coverage by means of four instead of six nucleotide recognition site restriction enzymes for fingerprinting mRNAs. Therefore, the method detects 70-80% of all transcripts, including non-coding transcripts, unknown and known genes. Moreover, the method requires no sequence information and so is applicable even to eukaryotes for which there is no genome information available.

The T cell surface--how well do we know it?

The overall degree of complexity of the T cell surface has been unclear, constraining our understanding of its biology. Using global gene expression analysis, we show that 111 of 374 genes encoding well-characterized leukocyte surface antigens are expressed by a resting cytotoxic T cell. Unexpectedly, of 97 stringently defined, T cell-specific transcripts with unknown functions that we identify, none encode proteins with the modular architecture characteristic of 80% of leukocyte surface antigens. Only two encode proteins with membrane topologies found exclusively in cell surface molecules. Our analysis indicates that the cell type-specific composition of the resting CD8+ T cell surface is now largely defined, providing an insight into the overall compositional complexity of the mammalian cell surface and a framework for formulating systematic models of T cell surface-dependent processes, such as T cell receptor triggering.

Customized molecular phenotyping by quantitative gene expression and pattern recognition analysis

Description of the molecular phenotypes of pathobiological processes in vivo is a pressing need in genomic biology. We have implemented a high-throughput real-time PCR strategy to establish quantitative expression profiles of a customized set of target genes. It enables rapid, reproducible data acquisition from limited quantities of RNA, permitting serial sampling of mouse blood during disease progression. We developed an easy to use statistical algorithm--Global Pattern Recognition--to readily identify genes whose expression has changed significantly from healthy baseline profiles. This approach provides unique molecular signatures for rheumatoid arthritis, systemic lupus erythematosus, and graft versus host disease, and can also be applied to defining the molecular phenotype of a variety of other normal and pathological processes.

Assessing gene expression variation in normal human tissues using GeneTag, a novel, global, sensitive profiling method

GeneTag is a novel expression profiling method that allows the visualization, quantification and identification of expressed genes-whether known or novel-in any species, tissue or cell type, independent of knowledge of the underlying sequence. Here we describe the application of this method to determine variation of gene expression in individual human liver samples and the identification of tissue-specific genes by comparing expression patterns across several human organs. Expression data are stored in a database for future reference and data analysis relies on proprietary software, which allows complex comparisons to be performed. Differentially expressed genes are quickly identified through a link to a sequence database. The results from our study underscore the importance of knowledge of individual variation of gene expression for the design and interpretation of transcript profiling experiments in the context of any biological question.

Gene discovery using a human vestibular schwannoma cDNA library constructed from a patient with neurofibromatosis type 2 (NF2)

BACKGROUND

Despite a strong association of schwannomin/merlin gene mutations with vestibular schwannoma formation, the regulatory mechanisms and biologic pathways involved are still largely unknown. The hypothesis of this study is that the genesis and growth characteristics of neurofibromatosis type 2 (NF2)-associated vestibular schwannomas are determined by genetic alterations that vary in gene transcript expression; this transcript expression includes oncogenic gene products that may be identified by construction and sequencing of a cDNA library from NF2-associated vestibular schwannoma.

METHODS

Approximately 3 mL of fresh tumor was obtained during resection of a 4-cm vestibular schwannoma from a patient with NF2. Poly(A)(+) mRNA was isolated, synthesized into double-stranded cDNA, and unidirectionally inserted into Uni-Zap XR (Stratagene, La Jolla, CA) bacteriophage vectors. Bacteriophage vectors containing cDNA inserts were processed into phagemids according to Uni-Zap XR protocol, and inserted vectors were sequenced and analyzed using BLAST software (National Institutes of Health, Bethesda, MD) with GenBank, EMBL, DDBJ, and PBD databases.

RESULTS

The cDNA library contained 2.4 million primary plaques. Inserts averaged 1.8 kilobases (kb) in length, with a range of 0.8 to 3.0 kb. BLAST multidatabase comparison of the sequence data obtained from 50 randomly selected clones yielded identification of 13 sequences representing known human genes and 17 sequences representing cloned sequences with unknown function. Three clones represented sequences not previously described in vestibular schwannomas but strongly implicated in oncogenesis within other tissues.

CONCLUSIONS

These data have implications for understanding the molecular mechanisms of vestibular schwannoma tumor biology. Identified genes may provide future diagnostic/prognostic markers and therapeutic targets.

Application of high-resolution capillary array electrophoresis with automated fraction collection for GeneCalling trade mark analysis of the yeast genomic DNA

Capillary array instrument was applied to transcript profiling of the yeast genomic DNA using GeneCalling trade mark chemistry. The instrument integrated a 12-capillary array for DNA separation with a replaceable sieving matrix, laser-induced fluorescence detection and an automated microfraction collector. The DNA fractions, exiting the separation capillaries, were continuously deposited in a 1536-well collection plate made of agarose gel. DNA fragments recovered from selected fractions were cloned and then sequenced. Over 80% of theoretically predicted fragments could be recovered in the collected fractions, cloned and sequenced with an average redundancy of threefold. Excellent correlation of the experimentally obtained sequences with the theoretically predicted gene fragments demonstrated the suitability of capillary array electrophoresis for micropreparative recovery of DNA fragments. This approach, useful especially for rapid DNA expression profiling of unknown genes for nonsequenced organisms, demonstrates the practical capability of the prototype multicapillary fraction collector.

Distinct transcriptional profiles of adrenocortical tumors uncovered by DNA microarray analysis

Comprehensive expression profiling of tumors using DNA microarrays has been used recently for molecular classification and biomarker discovery, as well as a tool to identify and investigate genes involved in tumorigenesis. Application of this approach to a cohort of benign and malignant adrenocortical tissues would be potentially informative in all of these aspects. In this study, we generated transcriptional profiles of 11 adrenocortical carcinomas (ACCs), 4 adrenocortical adenomas (ACAs), 3 normal adrenal cortices (NCs), and 1 macronodular hyperplasia (MNH) using Affymetrix HG_U95Av2 oligonucleotide arrays representing approximately 10,500 unique genes. The expression data set was used for unsupervised hierarchical cluster analysis as well as principal component analysis to visually represent the expression data. An analysis of variance on the three classes (NC, ACA plus MNH, and ACC) revealed 91 genes that displayed at least threefold differential expression between the ACC cohort and both the NC and ACA cohorts at a significance level of P < 0.01. Included in these 91 genes were those known to be up-regulated in adrenocortical tumors, such as insulin-like growth factor (IGF2), as well as novel differentially expressed genes such as osteopontin (SPP) and serine threonine kinase 15 (STK15). Increased expression of IGF2 was identified in 10 of 11 ACCs (90.9%) and was verified by quantitative reverse transcriptase-polymerase chain reaction. Select proliferation-related genes (TOP2A and Ki-67) were validated at the protein level using immunohistochemistry and adrenocortical tissue microarrays. Our results demonstrated significant and consistent gene expression changes in ACCs compared to benign adrenocortical lesions. Moreover, we identified several genes that represent potential diagnostic markers and may play a role in the pathogenesis of ACC.

Gene expression profile of human endothelial cells exposed to sustained fluid shear stress

Biomechanical forces can modulate endothelial phenotype through changes in gene expression. We hypothesized that physiological laminar shear stresses (LSS) act as differentiative stimuli on endothelial cells (EC) to alter gene expression, creating an antioxidant, anti-apoptotic and anti-proliferative environment. The transcriptional profile of cultured human umbilical vein endothelial cells (HUVEC) exposed to LSS was evaluated by GeneCalling; 107 genes demonstrated at least a twofold change in expression at 24 h (LSS vs. static). These flow-responsive genes represent a limited number of functional clusters that include transcription factors, antioxidants, signaling molecules, cell cycle regulators, and genes involved in cellular differentiation. Immunohistochemistry and in situ hybridization confirmed that many of these flow-responsive genes, including the novel basic helix-loop-helix transcription factor Hath6, are expressed in EC in vivo. Thus these data identify a limited set of flow-responsive genes expressed in the endothelium that may be responsible for the establishment and maintenance of the flow-adapted endothelial phenotype in vivo.

Generation of kidney transcriptomes using serial analysis of gene expression

Chronic renal disease initiation and progression remain incompletely understood. Genomewide expression monitoring should clarify the mechanisms which cause progressive renal disease by determining how clusters of genes coordinately change their activity. Serial analysis of gene expression (SAGE) is a technique of expression profiling which permits simultaneous and quantitative analysis of 9- to 13-bp sequence tags that correspond to unique mRNAs. Key principles of the technique are use of PCR in a manner to minimize distortion and serial concatenation of tags which facilitates sequencing and permits identification of many expressed genes in a single cDNA molecule. Tags are extracted from many concatenated sequences, counted using software, and identified by comparison with existing gene databases. In aggregate, gene expression profiles generated from a tag library comprise a transcriptome which represents a comprehensive and quantitative profile of genes expressed at the time of analysis. These global snapshots of gene expression patterns can better define basic cell biology and provide insights into disease pathogenesis by simultaneously determining the net consequences of gene-gene and gene-environment interactions on expression of thousands of genes. Rather than applying a priori assumptions (i.e., hypothesis testing), transcriptome analysis is hypothesis generating and requires no prior knowledge of gene expression. SAGE kidney transcriptomes, from normal animals and animals with progressive kidney disease, are being produced and can be analyzed for novel pathogenetic mechanisms. The use of SAGE and other genomic and proteomic tools should result in a better understanding of kidney disease pathogenesis and in identification of new therapeutic targets.

Identification of differentially expressed genes in Dreissena polymorpha exposed to contaminants

Development of transcriptome analysis methods such as differential display PCR and construction of subtractive libraries now makes it possible to profile gene expression in response to xenobiotic exposure. As an example of application of these methods, zebra mussels (Dreissena polymorpha) were treated with various contaminants such as Aroclor 1254, 3-methylcholanthrene, chrysene and atrazine. A total of 242 mRNAs were identified as differentially expressed. Analysis of these mRNAs should provide valuable information regarding detoxification mechanisms in this bivalve species. In addition, the use of cDNA array technology applied to these gene products may constitute a multi-marker approach to monitor the effect of contamination on this aquatic species.

Staking out novelty on the genomic frontier

Recent advances in genomics include global assessment and classification of genome content, high-throughput biological pathway construction, systematic identification of previously unpredicted genes and the in vitro creation of novel motifs with biological function not found in nature (extra-genomic gene discovery). The ability to make global surveys of transcriptomes has given rise to fields such as pharmacogenomics and toxicogenomics. These applications of genomics technologies, with conventional drug assessment methodologies, will lead to more tolerable drugs and a better understanding of clinical populations. Integration of pathway mapping, using proteomics married to expression, will also significantly affect how new therapeutics are discovered as cross-biological cross-pathway interactions lead to novel drug targets and better predictions of drug tolerance.

Application of toxicogenomics to drug development

The application of advanced modern biomedical and chemical research technologies in the pharmaceutical industry has led to a significant increase in the number of potential drug targets and lead candidates. Whereas the drug discovery process is enhanced significantly, the failure rate of new compounds due to toxicity remains very high. The pharmaceutical industry is setting high hopes on the new discipline of toxicogenomics to revolutionize the process of drug toxicity assessment by reducing the bottleneck of new drug candidates and minimizing late-stage developmental failures. Toxicogenomics is expected to facilitate the efficient screening of new compounds at an early stage, resulting in significant savings of time and cost associated with new drug development. In this review, a general description of the new discipline of toxicogenomics and its potential impact on the safety assessment of new drugs in the pharmaceutical industry is provided. An overview of the key issues and questions that are confronting investigators in this field today is also given as well as a prospective view of the future of this new discipline.

Heterogeneity of gene expression in human atheroma unmasked using cDNA representational difference analysis

The rupture of an atherosclerotic plaque can have profound consequences, such as myocardial or cerebrovascular infarction. The complex interactions of vascular smooth muscle cells (VSMCs) with inflammatory and immune cells are thought to contribute to both plaque genesis and stability. Key to our understanding of these processes is the identification of genes expressed in human atheromatous lesions. We have employed cDNA representational difference analysis (RDA) to investigate the differences in gene expression between normal and atherosclerotic human vessels. Thirty-one cDNA clones representing sequences expressed in atheroma were isolated, many of which encoded components of inflammatory and immune pathways. The reciprocal experiment, to identify genes expressed in the healthy vasculature, identified two genes associated with the contractile functions of VSMCs. Semiquantitative RT-PCR analysis of expression of these genes in forty samples, derived from healthy and atheromatous vessels, demonstrated marked heterogeneity of gene expression between lesions, although several of the genes were preferentially expressed in atherosclerotic lesions. In situ hybridization identified subsets of macrophages at sites of neovascularization within the lesion and intimal VSMCs as expressing the disease-associated genes. In conclusion, cDNA RDA is a useful, fast, and efficient technique for studying differential gene expression particularly when clinical material is limiting.

Differential display in the time of microarrays

A brief overview of major methods used for genome-wide expression profiling is presented. Special attention is devoted to ordered differential display, subtractive hybridization and DNA microarrays. Future prospects of comparative gene expression studies using combinations of differential display methods and microarray technology are outlined.

The use of animal models in expression pharmacogenomic analyses

Expression pharmacogenomics applies genome/proteome scale differential expression technologies to both in vivo and in vitro models of drug response to identify candidate markers correlative with and predictive of drug toxicity and efficacy. It is anticipated to streamline drug development by triaging towards lead compounds and clinical candidates that maximize efficacy while minimizing safety risks. As the majority of expression pharmacogenomics will be performed on preclinical therapeutic candidates, compatibility with favored preclinical animal model systems will be essential. This review will address expression pharmacogenomics in the context of those animal model systems commonly used for pharmacokinetic, pharmacodynamic and toxicologic analyses. Specific discussions will cover: (A) relative robustness of genomic and proteomic technology platforms used to generate drug response data in critical model systems; (B) animal handling, treatment and other experimental design optimizations; (C) data analysis strategies for extracting and validating candidate pharmacogenomic markers; and (D) overarching limitations in applying expression pharmacogenomics to animal model systems.