Systematic comparison of the fidelity of aRNA, mRNA and T-RNA on gene expression profiling using cDNA microarray

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.

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.

Shewanella baltica Ecotypes Have Wide Transcriptional Variation under the Same Growth Conditions

In bacterial populations, subtle expressional differences may promote ecological specialization through the formation of distinct ecotypes. In a barrier-free habitat, this process most likely precedes population divergence and may predict speciation events. To examine this, we used four sequenced strains of the bacterium Shewanella baltica, OS155, OS185, OS195, and OS223, as models to assess transcriptional variation and ecotype formation within a prokaryotic population. All strains were isolated from different depths throughout a water column of the Baltic Sea, occupying different ecological niches characterized by various abiotic parameters. Although the genome sequences are nearly 100% conserved, when grown in the laboratory under standardized conditions, all strains exhibited different growth rates, suggesting significant expressional variation. Using the Ecotype Simulation algorithm, all strains were considered to be discrete ecotypes when compared to 32 other S. baltica strains isolated from the same water column, suggesting ecological divergence. Next, we employed custom microarray slides containing oligonucleotide probes representing the core genome of OS155, OS185, OS195, and OS223 to detect natural transcriptional variation among strains grown under identical conditions. Significant transcriptional variation was noticed among all four strains. Differentially expressed gene profiles seemed to coincide with the metabolic signatures of the environment at the original isolation depth. Transcriptional pattern variations such as the ones highlighted here may be used as indicators of short-term evolution emerging from the formation of bacterial ecotypes. IMPORTANCE Eukaryotic studies have shown considerable transcriptional variation among individuals from the same population. It has been suggested that natural variation in eukaryotic gene expression may have significant evolutionary consequences and may explain large-scale phenotypic divergence of closely related species, such as humans and chimpanzees (M.-C. King and A. C. Wilson, Science 188:107-116, 1975, http://dx.doi.org/10.1126/science.1090005; M. F. Oleksiak, G. A. Churchill, and D. L. Crawford, Nat Genet 32:261-266, 2002, http://dx.doi.org/10.1038/ng983). However, natural variation in gene expression is much less well understood in prokaryotic organisms. In this study, we used four sequenced strains of the marine bacterium Shewanella baltica to better understand the natural transcriptional divergence of a stratified prokaryotic population. We found substantial low-magnitude expressional variation among the four S. baltica strains cultivated under identical laboratory conditions. Collectively, our results indicate that transcriptional variation is an important factor for ecological speciation.

Transcriptome profiling of human hippocampus dentate gyrus granule cells in mental illness

This study is, to the best of our knowledge, the first application of whole transcriptome sequencing (RNA-seq) to cells isolated from postmortem human brain by laser capture microdissection. We investigated the transcriptome of dentate gyrus (DG) granule cells in postmortem human hippocampus in 79 subjects with mental illness (schizophrenia, bipolar disorder, major depression) and nonpsychiatric controls. We show that the choice of normalization approach for analysis of RNA-seq data had a strong effect on results; under our experimental conditions a nonstandard normalization method gave superior results. We found evidence of disrupted signaling by miR-182 in mental illness. This was confirmed using a novel method of leveraging microRNA genetic variant information to indicate active targeting. In healthy subjects and those with bipolar disorder, carriers of a high- vs those with a low-expressing genotype of miR-182 had different levels of miR-182 target gene expression, indicating an active role of miR-182 in shaping the DG transcriptome for those subject groups. By contrast, comparing the transcriptome between carriers of different genotypes among subjects with major depression and schizophrenia suggested a loss of DG miR-182 signaling in these conditions.

Enlarging cells initiating apomixis in Hieracium praealtum transition to an embryo sac program prior to entering mitosis

Hieracium praealtum forms seeds asexually by apomixis. During ovule development, sexual reproduction initiates with megaspore mother cell entry into meiosis and formation of a tetrad of haploid megaspores. The sexual pathway ceases when a diploid aposporous initial (AI) cell differentiates, enlarges, and undergoes mitosis, forming an aposporous embryo sac that displaces sexual structures. Embryo and endosperm development in aposporous embryo sacs is fertilization independent. Transcriptional data relating to apomixis initiation in Hieracium spp. ovules is scarce and the functional identity of the AI cell relative to other ovule cell types is unclear. Enlarging AI cells with undivided nuclei, early aposporous embryo sacs containing two to four nuclei, and random groups of sporophytic ovule cells not undergoing these events were collected by laser capture microdissection. Isolated amplified messenger RNA samples were sequenced using the 454 pyrosequencing platform and comparatively analyzed to establish indicative roles of the captured cell types. Transcriptome and protein motif analyses showed that approximately one-half of the assembled contigs identified homologous sequences in Arabidopsis (Arabidopsis thaliana), of which the vast majority were expressed during early Arabidopsis ovule development. The sporophytic ovule cells were enriched in signaling functions. Gene expression indicative of meiosis was notably absent in enlarging AI cells, consistent with subsequent aposporous embryo sac formation without meiosis. The AI cell transcriptome was most similar to the early aposporous embryo sac transcriptome when comparing known functional annotations and both shared expressed genes involved in gametophyte development, suggesting that the enlarging AI cell is already transitioning to an embryo sac program prior to mitotic division.

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.

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.

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.

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.

RNA preparation and characterization for gene expression studies

Much information can be obtained from knowledge of the relative expression level of each gene in the transcriptome. With the current advances in technology as little as a single cell is required as starting material for gene expression experiments. The mRNA from a single cell may be linearly amplified to an excess of 10(6)-fold. Reverse transcription and fluorescent labeling of the amplified RNA yields a stable target for subsequent hybridization to DNA microarrays.

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.

Epstein–Barr virus-encoded LMP1 induces a hyperproliferative and inflammatory gene expression programme in cultured keratinocytes

SCC12F cells are a line of keratinocytes that retain the capacity for terminal differentiation in vitro. We showed previously that the Epstein–Barr virus (EBV)-encoded oncogene latent membrane protein 1 (LMP1) altered SCC12F morphology in vitro, downregulated cell–cell-adhesion molecule expression and promoted cell motility. In organotypic raft culture, LMP1-expressing cells failed to stratify and formed poorly organized structures which displayed impaired terminal differentiation. To understand better the mechanism(s) by which LMP1 induces these effects, we generated SCC12F cells in which LMP1 expression is inducible. Following induction, these cells exhibited phenotypic changes similar to those observed previously and allowed us to investigate the effects of LMP1 expression on cellular pathways associated with growth, differentiation and morphology. Using microarrays and a number of confirmatory techniques, we identified sets of differentially expressed genes that are characteristically expressed in inflammatory and hyperproliferative epidermis, including chemokines, cytokines and their receptors, growth factors involved in promoting epithelial cell motility and proliferation and signalling molecules that regulate actin filament reorganization and cell movement. Among the genes whose expression was differentially induced significantly by LMP1, the induction of IL-1β and IL-1α was of particular interest, as many of the LMP1-regulated genes identified are established targets of these cytokines. Our findings suggest that alterations in the IL-1 signalling network may be responsible for many of the changes in host-cell gene expression induced in response to LMP1. Identification of these LMP1-regulated genes helps to define the mechanism(s) by which this oncoprotein influences cellular pathways that regulate terminal differentiation, cell motility and inflammation.

Simultaneously detection of genomic and expression alterations in prostate cancer using cDNA microarray

BACKGROUND

Prostate cancer is a common disease among men but the knowledge of the prostate carcinogenesis is still limited.

METHODS

cDNA microarray-based comparative genomic hybridization (CGH) and expression profiling were performed to screen the genomic and the expression changes in prostate cancer respectively. The two data were integrated to study the influence of genomic aberrations on gene expression and seek for the genes with their expression affected by the genomic aberrations. Real-time PCR was performed to evaluate the array data.

RESULTS

Array-based CGH detected gains at 2q, 3p/q, 5q, 6q, 8q, 9p, 10p/q, 11q, 12p, 14q, and 19p/q and losses at 1p, 2p, 4q, 6p/q, 7p, 11p/q, 12q, 17p/q, 19p/q, and Xp/q in more than 20% prostate tumors and narrowed these aberrations. For example, the gain of 8q was mapped to five minimal regions. Novel aberrations were also identified, such as loss at Xq21.33-q22.2. Expression profiling discovered the significant biological processes involved in the prostate carcinogenesis, such as exogenous antigen presentation via MHC class II and protein ubiquitination. Integration analysis revealed a weak positive correlation between genomic copy number and gene expression level. Fifty-three genes showed their expression directly affected by the genomic aberrations possibly, including more than one member of Ras superfamily and major histocompatibility complex (MHC). These genes are involved in multiple biological processes.

CONCLUSIONS

Integration of the CGH and expression data provided more information than separate analysis. Although the direct influence of genomic aberrations on gene expression seems weak, the influence can be extended by indirect regulation through a few directly affected genes. Because the influence can be persistent, the genes directly affected by the genomic aberrations may play key roles in the prostate carcinogenesis and are worth further analysis.

Technical variables in high-throughput miRNA expression profiling: much work remains to be done

MicroRNA (miRNA) gene expression profiling has provided important insights into plant and animal biology. However, there has not been ample published work about pitfalls associated with technical parameters in miRNA gene expression profiling. One source of pertinent information about technical variables in gene expression profiling is the separate and more well-established literature regarding mRNA expression profiling. However, many aspects of miRNA biochemistry are unique. For example, the cellular processing and compartmentation of miRNAs, the differential stability of specific miRNAs, and aspects of global miRNA expression regulation require specific consideration. Additional possible sources of systematic bias in miRNA expression studies include the differential impact of pre-analytical variables, substrate specificity of nucleic acid processing enzymes used in labeling and amplification, and issues regarding new miRNA discovery and annotation. We conclude that greater focus on technical parameters is required to bolster the validity, reliability, and cultural credibility of miRNA gene expression profiling studies.

Application of laser capture microdissection and differential display technique for screening of pathogenic genes involved in endometrial carcinoma

The objective is to eliminate the interference from other cell types; gene fragments involved in endometrial carcinoma (EC) are screened and cloned. Human normal endometrial glandular epithelia and EC cells were harvested with laser capture microdissection (LCM). The purification and concentration of minimal RNA were used to screen differential expressed gene fragments involved in EC by fluoro differential display polymerase chain reaction (FDD-PCR). The differential gene fragments were cloned, sequenced, and then identified by reverse Northern blot hybridization. Positive fragments were analyzed with basic local alignment search tool (BLAST). Cyclin-dependent kinase 7 (CDK7) expressions in EC and normal endometrial tissue were tested by immunohistochemical staining. Of 38 differential bands, 3 bands were of high expression in normal endometrium and 35 in EC. Six positive differential gene fragments were obtained and BLAST analysis for them suggested that L1.1 was homologous (99% identical) to the CDK7; L1.9 had a 99% homology with protein phosphatase 1 regulatory (inhibitor) subunit 12 A (PPP1R12A); L1.21 and L1.22 showed a 100% homology with cellular repressor of E1A-stimulated genes 1 (CREG); and L1.25 and L1.26 indicated more than 98% homology with solute carrier family 39 (zinc transporter), member 10 (SLC39A10). Immunohistochemistry revealed that CDK7 expression was higher in EC than in normal endometrium. We conclude that pathogenic genes involved in EC are obtained with LCM and FDD-PCR. It has been first found that CDK7, PPP1R12A, CREG, and SLC39A10 are correlative with EC from gene level. CDK7 is strongly associated with EC and can be used as potential molecular marker of EC for further studies

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.

Single molecule transcription profiling with AFM

Established techniques for global gene expression profiling, such as microarrays, face fundamental sensitivity constraints. Due to greatly increasing interest in examining minute samples from micro-dissected tissues, including single cells, unorthodox approaches, including molecular nanotechnologies, are being explored in this application. Here, we examine the use of single molecule, ordered restriction mapping, combined with AFM, to measure gene transcription levels from very low abundance samples. We frame the problem mathematically, using coding theory, and present an analysis of the critical error sources that may serve as a guide to designing future studies. We follow with experiments detailing the construction of high density, single molecule, ordered restriction maps from plasmids and from cDNA molecules, using two different enzymes, a result not previously reported. We discuss these results in the context of our calculations.

Assessing the reliability of amplified RNA used in microarrays: a DUMB table approach

A certain minimal amount of RNA from biological samples is necessary to perform a microarray experiment with suitable replication. In some cases, the amount of RNA available is insufficient, necessitating RNA amplification prior to target synthesis. However, there is some uncertainty about the reliability of targets that have been generated from amplified RNA, because of nonlinearity and preferential amplification. This current work develops a straightforward strategy to assess the reliability of microarray data obtained from amplified RNA. The tabular method we developed, which utilises a Down-Up-Missing-Below (DUMB) classification scheme, shows that microarrays generated with amplified RNA targets are reliable within constraints. There was an increase in false negatives because of the need for increased filtering. Furthermore, this analysis method is generic and can be broadly applied to evaluate all microarray data. A copy of the Microsoft Excel spreadsheet is available upon request from Edward Bearden.

In vitro transcription amplification and labeling methods contribute to the variability of gene expression profiling with DNA microarrays

The effect of different amplification and labeling methods on DNA microarray expression results has not been previously delineated. To analyze the variation associated with widely accepted T7-based RNA amplificationand labeling methods, aliquots of the Stratagene Human Universal Reference RNA were labeled using three eukaryotic target preparation methods followed by uniform replicate array hybridization (Affymetrix U95Av2). Method-dependent variability was observed in the yield and size distribution of labeled products, as well as in the gene expression results. A significant increase in short transcripts, when compared to unamplified mRNA, was observed in methods with long in vitro transcription reactions. Intramethod reproducibility showed correlation coefficients >0.99, whereas intermethod comparisons showed coefficients ranging from 0.94 to 0.98 and a nearly twofold increase in coefficient of variation. Fold amplification for each method positively correlated with the number of genes present. Our experiments uncovered two factors that introduced significant bias in gene expression data: the number of labeled nucleotides, which introduces sequence-dependent bias, and the length of the in vitro transcription reaction, which introduces transcript size-dependent bias. This study provides evidence that variability in expression data may be caused, in part, by differences in amplification and labeling protocols.

Reverse transcription using random pentadecamer primers increases yield and quality of resulting cDNA

Reverse transcription of RNA is an invaluable method for gene expression analysis by real-time PCR or microarray methods. Random primers of varying lengths were compared with respect to their efficiency of priming reverse transcription reactions. The results showed that 15-nucleotide-long random oligonucleotides (pentadecamers) consistently yielded at least 2-fold as much cDNA as did random hexamers using either poly(A) RNA or an amplified version of messenger RNA (aRNA) as a template. The cDNA generated using pentadecamers did not differ in size distribution or the amount of incorporated label compared with cDNA generated with random hexamers. The increased efficiency of priming using random pentadecamers resulted in reverse transcription of >80% of the template aRNA, while random hexamers induced reverse transcription of only 40% of the template aRNA. This suggests a better coverage of the transcriptome when using random pentadecamers over random hexamers. Using the same amount of aRNA as starting material, random pentadecamer-primed reactions resulted in 11-fold more genes being detected in whole transcriptome DNA microarray experiments than random hexamer-primed reactions. The results indicate that random pentadecamers can replace random hexamers in reverse transcription reactions on both poly(A) RNA and amplified RNA, resulting in higher cDNA yields and quality.

Spotting and validation of a genome wide oligonucleotide chip with duplicate measurement of each gene

The quality of DNA microarray based gene expression data relies on the reproducibility of several steps in a microarray experiment. We have developed a spotted genome wide microarray chip with oligonucleotides printed in duplicate in order to minimise undesirable biases, thereby optimising detection of true differential expression. The validation study design consisted of an assessment of the microarray chip performance using the MessageAmp and FairPlay labelling kits. Intraclass correlation coefficient (ICC) was used to demonstrate that MessageAmp was significantly more reproducible than FairPlay. Further examinations with MessageAmp revealed the applicability of the system. The linear range of the chips was three orders of magnitude, the precision was high, as 95% of measurements deviated less than 1.24-fold from the expected value, and the coefficient of variation for relative expression was 13.6%. Relative quantitation was more reproducible than absolute quantitation and substantial reduction of variance was attained with duplicate spotting. An analysis of variance (ANOVA) demonstrated no significant day-to-day variation.

Options available for profiling small samples: a review of sample amplification technology when combined with microarray profiling

The possibility of performing microarray analysis on limited material has been demonstrated in a number of publications. In this review we approach the technical aspects of mRNA amplification and several important implicit consequences, for both linear and exponential procedures. Amplification efficiencies clearly allow profiling of extremely small samples. The conservation of transcript abundance is the most important issue regarding the use of sample amplification in combination with microarray analysis, and this aspect has generally been found to be acceptable, although demonstrated to decrease in highly diluted samples. The fact that variability and discrepancies in microarray profiles increase with minute sample sizes has been clearly documented, but for many studies this does appear to have affected the biological conclusions. We suggest that this is due to the data analysis approach applied, and the consequence is the chance of presenting misleading results. We discuss the issue of amplification sensitivity limits in the light of reports on fidelity, published data from reviewed articles and data analysis approaches. These are important considerations to be reflected in the design of future studies and when evaluating biological conclusions from published microarray studies based on extremely low input RNA quantities.

Laser microdissection of plant tissue: what you see is what you get

Laser microdissection (LM) utilizes a cutting or harvesting laser to isolate specific cells from histological sections; the process is guided by microscopy. This provides a means of removing selected cells from complex tissues, based only on their identification by microscopic appearance, location, or staining properties (e.g., immunohistochemistry, reporter gene expression, etc.). Cells isolated by LM can be a source of cell-specific DNA, RNA, protein or metabolites for subsequent evaluation of DNA modifications, transcript/protein/metabolite profiling, or other cell-specific properties that would be averaged with those of neighboring cell types during analysis of undissected complex tissues. Plants are particularly amenable to the application of LM; the highly regular tissue organization and stable cell walls of plants facilitate the visual identification of most cell types even in unstained tissue sections. Plant cells isolated by LM have been the starting point for a variety of genomic and metabolite studies of specific cell types.

Limitations of mRNA amplification from small-size cell samples

Background

Global mRNA amplification has become a widely used approach to obtain gene expression profiles from limited material. An important concern is the reliable reflection of the starting material in the results obtained. This is especially important with extremely low quantities of input RNA where stochastic effects due to template dilution may be present. This aspect remains under-documented in the literature, as quantitative measures of data reliability are most often lacking. To address this issue, we examined the sensitivity levels of each transcript in 3 different cell sample sizes. ANOVA analysis was used to estimate the overall effects of reduced input RNA in our experimental design. In order to estimate the validity of decreasing sample sizes, we examined the sensitivity levels of each transcript by applying a novel model-based method, TransCount.

Results

From expression data, TransCount provided estimates of absolute transcript concentrations in each examined sample. The results from TransCount were used to calculate the Pearson correlation coefficient between transcript concentrations for different sample sizes. The correlations were clearly transcript copy number dependent. A critical level was observed where stochastic fluctuations became significant. The analysis allowed us to pinpoint the gene specific number of transcript templates that defined the limit of reliability with respect to number of cells from that particular source. In the sample amplifying from 1000 cells, transcripts expressed with at least 121 transcripts/cell were statistically reliable and for 250 cells, the limit was 1806 transcripts/cell. Above these thresholds, correlation between our data sets was at acceptable values for reliable interpretation.

Conclusion

These results imply that the reliability of any amplification experiment must be validated empirically to justify that any gene exists in sufficient quantity in the input material. This finding has important implications for any experiment where only extremely small samples such as single cell analyses or laser captured microdissected cells are available.

Macrodissection versus microdissection of rectal carcinoma: minor influence of stroma cells to tumor cell gene expression profiles

Background

The molecular determinants of carcinogenesis, tumor progression and patient prognosis can be deduced from simultaneous comparison of thousands of genes by microarray analysis. However, the presence of stroma cells in surgically excised carcinoma tissues might obscure the tumor cell-specific gene expression profiles of these samples. To circumvent this complication, laser microdissection can be performed to separate tumor epithelium from the surrounding stroma and healthy tissue. In this report, we compared RNAs isolated from macrodissected, of which only surrounding healthy tissue had been removed, and microdissected rectal carcinoma samples by microarray analysis in order to determine the most reliable approach to detect the expression of tumor cell-derived genes by microarray analysis.

Results

As microdissection yielded low tissue and RNA quantities, extra rounds of mRNA amplification were necessary to obtain sufficient RNA for microarray experiments. These second rounds of amplification influenced the gene expression profiles. Moreover, the presence of stroma cells in macrodissected samples had a minor contribution to the tumor cell gene expression profiles, which can be explained by the observation that more RNA is extracted from tumor epithelial cells than from stroma.

Conclusion

These data demonstrate that the more convenient procedure of macrodissection can be adequately used and yields reliable data regarding the identification of tumor cell-specific gene expression profiles.

Transit tumor retrieval preserves RNA fidelity and obviates snap-freezing

Genetic expression profiling is enabling investigators to discover new diagnostic and possibly therapeutic pathways in sarcoma biology. To draw substantial conclusions from these molecular analyses, adequate tissue samples must be accrued. Beyond cohort size, the most variable and limiting aspect of doing gene expression analyses on fresh human tissue is the preservation of labile ribonucleic acids extracted from clinical specimens. We have developed a novel retrieval protocol that is readily amenable to the clinical constraints placed on surgeons and pathologists that minimizes variables that can corrupt ribonucleic acid fidelity. We evaluate critically genomic message integrity of mesenchymal tumors derived from transcontinental inter-institutional collaboration. Intact total ribonucleic acid was isolated and assessed for quality and quantity. Ribosomal RNA integrity was quantified using a bioanalyzer. Ribonucleic acid from 42 mesenchymal tumors was isolated and quantified, with selected samples amplified. The mean ribosomal ratios for collaborative institutions ranged from 1.0 to 1.18. Samples remained at 4 degrees C before processing from 1 to 17 days. Tumors stabilized using this protocol retained total ribonucleic acid integrity suitable for amplification and genomic expression analysis regardless of the institutional source or preprocessing duration, enabling a potential consortium of investigators to collaborate in the expression profiling of sarcomas.

LEVEL OF EVIDENCE

Diagnostic study, Level III-3 (no consistently applied gold standard). See the Guidelines for Authors for a complete description of levels of evidence.

RNA amplification for successful gene profiling analysis

The study of clinical samples is often limited by the amount of material available to study. While proteins cannot be multiplied in their natural form, DNA and RNA can be amplified from small specimens and used for high-throughput analyses. Therefore, genetic studies offer the best opportunity to screen for novel insights of human pathology when little material is available. Precise estimates of DNA copy numbers in a given specimen are necessary. However, most studies investigate static variables such as the genetic background of patients or mutations within pathological specimens without a need to assess proportionality of expression among different genes throughout the genome. Comparative genomic hybridization of DNA samples represents a crude exception to this rule since genomic amplification or deletion is compared among different specimens directly. For gene expression analysis, however, it is critical to accurately estimate the proportional expression of distinct RNA transcripts since such proportions directly govern cell function by modulating protein expression. Furthermore, comparative estimates of relative RNA expression at different time points portray the response of cells to environmental stimuli, indirectly informing about broader biological events affecting a particular tissue in physiological or pathological conditions. This cognitive reaction of cells is similar to the detection of electroencephalographic patterns which inform about the status of the brain in response to external stimuli. As our need to understand human pathophysiology at the global level increases, the development and refinement of technologies for high fidelity messenger RNA amplification have become the focus of increasing interest during the past decade. The need to increase the abundance of RNA has been met not only for gene specific amplification, but, most importantly for global transcriptome wide, unbiased amplification. Now gene-specific, unbiased transcriptome wide amplification accurately maintains proportionality among all RNA species within a given specimen. This allows the utilization of clinical material obtained with minimally invasive methods such as fine needle aspirates (FNA) or cytological washings for high throughput functional genomics studies. This review provides a comprehensive and updated discussion of the literature in the subject and critically discusses the main approaches, the pitfalls and provides practical suggestions for successful unbiased amplification of the whole transcriptome in clinical samples.

Analysis of global mRNA expression in human skeletal muscle during recovery from endurance exercise

To search for novel transcriptional pathways that are activated in skeletal muscle after endurance exercise, we used cDNA microarrays to measure global mRNA expression after an exhaustive bout of high-intensity cycling (approximately 75 min). Healthy, young, sedentary males performed the cycling bout, and skeletal muscle biopsies were taken from the vastus lateralis before, and at 3 and 48 h after exercise. We examined mRNA expression in individual muscle samples from four subjects using cDNA microarrays, used repeated-measures significance analysis of microarray (SAM) to determine statistically significant expression changes, and confirmed selected results using real-time RT-PCR. In total, the expression of 118 genes significantly increased 3 h postcycling and 8 decreased. At 48 h, the expression of 29 genes significantly increased and 5 decreased. Many of these are potentially important novel genes involved in exercise recovery and adaptation, including several involved in 1) metabolism and mitochondrial biogenesis (FOXO1, PPARdelta, PPARgamma, nuclear receptor binding protein 2, IL-6 receptor, ribosomal protein L2, aminolevulinate delta-synthase 2); 2) the oxidant stress response (metalothioneins 1B, 1F, 1G, 1H, 1L, 2A, 3, interferon regulatory factor 1); and 3) electrolyte transport across membranes [Na+-K+-ATPase (beta3), SERCA3, chloride channel 4]. Others include genes involved in cell stress, proteolysis, apoptosis, growth, differentiation, and transcriptional activation, as well as all three nuclear receptor subfamily 4A family members (Nur77, Nurr1, and Nor1). This study is the first to characterize global mRNA expression during recovery from endurance exercise, and the results provide potential insight into 1) the transcriptional contributions to homeostatic recovery in human skeletal muscle after endurance exercise, and 2) the transcriptional contributions from a single bout of endurance exercise to the adaptive processes that occur after a period of endurance exercise training.

Purity for clarity: the need for purification of tumor cells in DNA microarray studies

It is now well established that gene expression profiling using DNA microarrays can provide novel information about various types of hematological malignancies, which may lead to identification of novel diagnostic markers. However, to successfully use microarrays for this purpose, the quality and reproducibility of the procedure need to be guaranteed. The quality of microarray analyses may be severely reduced, if variable frequencies of nontarget cells are present in the starting material. To systematically investigate the influence of different types of impurity, we determined gene expression profiles of leukemic samples containing different percentages of nonleukemic leukocytes. Furthermore, we used computer simulations to study the effect of different kinds of impurity as an alternative to conducting hundreds of microarray experiments on samples with various levels of purity. As expected, the percentage of erroneously identified genes rose with the increase of contaminating nontarget cells in the samples. The simulations demonstrated that a tumor load of less than 75% can lead to up to 25% erroneously identified genes. A tumor load of at least 90% leads to identification of at most 5% false-positive genes. We therefore propose that in order to draw well-founded conclusions, the percentage of target cells in microarray experiment samples should be at least 90%.

Comparison of two probe preparation methods using long oligonucleotide microarrays

The use of oligonucleotides as a capture platform for microarray-based experiments is gaining popularity. Oligonucleotide-based microarrays involving various probe preparations have been compared by a number of researchers. Limited data are available, however, regarding the concordances and efficacies of various probe preparations on long oligonucleotide-based microarrays. Accordingly, the current investigation assesses two labeling methods, namely Atlas™ PowerScript™ fluorescent cDNA (random priming) and T7 in vitro transcription cRNA [poly(T) priming] labeling kits. Our data revealed that a high degree of reproducibility among the examined genes for each assay used with correlation coefficients of 0.93 and 0.94 for random priming and poly(T) priming, respectively. It is worthy of note, however, that when the two assaying methods were compared, the data showed a poor correlation coefficient. A confirmatory step involving real-time reverse transcription PCR (RT-PCR) of 18 selected genes favors the superiority of the cDNA fluorescent labeling over the T7 labeling method. Overall, the microarray results generated by the poly(T) priming methodology should be viewed cautiously even when high reproducibility is evident.

Evaluation of sense-strand mRNA amplification by comparative quantitative PCR

Background

RNA amplification is required for incorporating laser-capture microdissection techniques into microarray assays. However, standard oligonucleotide microarrays contain sense-strand probes, so traditional T7 amplification schemes producing anti-sense RNA are not appropriate for hybridization when combined with conventional reverse transcription labeling methods. We wished to assess the accuracy of a new sense-strand RNA amplification method by comparing ratios between two samples using quantitative real-time PCR (qPCR), mimicking a two-color microarray assay.

Results

We performed our validation using qPCR. Three samples of rat brain RNA and three samples of rat liver RNA were amplified using several kits (Ambion messageAmp, NuGen Ovation, and several versions of Genisphere SenseAmp). Results were assessed by comparing the liver/brain ratio for 192 mRNAs before and after amplification. In general, all kits produced strong correlations with unamplified RNAs. The SenseAmp kit produced the highest correlation, and was also able to amplify a partially degraded sample accurately.

Conclusion

We have validated an optimized sense-strand RNA amplification method for use in comparative studies such as two-color microarrays.