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

A Detailed Overview About the Single-Cell Analyses of Solid Tumors Focusing on Colorectal Cancer

In recent years, the evolution of the molecular biological technical background led to the widespread application of single-cell sequencing, a versatile tool particularly useful in the investigation of tumor heterogeneity. Even 10 years ago the comprehensive characterization of colorectal cancers by The Cancer Genome Atlas was based on measurements of bulk samples. Nowadays, with single-cell approaches, tumor heterogeneity, the tumor microenvironment, and the interplay between tumor cells and their surroundings can be described in unprecedented detail. In this review article we aimed to emphasize the importance of single-cell analyses by presenting tumor heterogeneity and the limitations of conventional investigational approaches, followed by an overview of the whole single-cell analytic workflow from sample isolation to amplification, sequencing and bioinformatic analysis and a review of recent literature regarding the single-cell analysis of colorectal cancers.

Molecular tools for the characterization of seizure susceptibility in genetic rodent models of epilepsy

Epilepsy is a chronic neurological disorder characterized by abnormal neuronal activity that arises from imbalances between excitatory and inhibitory synapses, which are highly correlated to functional and structural changes in specific brain regions. The difference between the normal and the epileptic brain may harbor genetic alterations, gene expression changes, and/or protein alterations in the epileptogenic nucleus. It is becoming increasingly clear that such differences contribute to the development of distinct epilepsy phenotypes. The current major challenges in epilepsy research include understanding the disease progression and clarifying epilepsy classifications by searching for novel molecular biomarkers. Thus, the application of molecular techniques to carry out comprehensive studies at deoxyribonucleic acid, messenger ribonucleic acid, and protein levels is of utmost importance to elucidate molecular dysregulations in the epileptic brain. The present review focused on the great diversity of technical approaches available and new research methodology, which are already being used to study molecular alterations underlying epilepsy. We have grouped the different techniques according to each step in the flow of information from DNA to RNA to proteins, and illustrated with specific examples in animal models of epilepsy, some of which are our own. Separately and collectively, the genomic and proteomic techniques, each with its own strengths and limitations, provide valuable information on molecular mechanisms underlying seizure susceptibility and regulation of neuronal excitability. Determining the molecular differences between genetic rodent models of epilepsy and their wild-type counterparts might be a key in determining mechanisms of seizure susceptibility and epileptogenesis as well as the discovery and development of novel antiepileptic agents. This article is part of the Special Issue "NEWroscience 2018".

Retrospective analysis of cancer-specific gene expression panel for thyroid fine needle aspiration specimens

BACKGROUND

While molecular testing is a promising strategy for preoperative assessment of cytologically indeterminate thyroid nodules, thyroid fine needle aspiration biopsy (FNA) presents unique challenges for molecular assays, including contaminating peripheral blood mononuclear cells (PBMC) and variable numbers of evaluable epithelial thyroid cells. Moreover, the newly recognized entity, noninvasive follicular thyroid neoplasm with papillary-like nuclear features (NIFTP), has added an additional challenge to the currently available molecular diagnostic platforms. New diagnostic tools are still needed to correctly distinguish benign and malignant thyroid nodules preoperatively.

METHODS

Twenty-two transcript splice variants from 12 genes we previously identified as discriminating benign from malignant thyroid nodules were characterized in 80 frozen thyroid tumors from 8 histological subtypes. Isoforms detectable in PBMC were excluded, and the 5 most discriminating isoforms were further validated by real-time quantitative PCR (qPCR) on intraoperative FNA samples from 59 malignant tumors, 55 benign nodules, and 23 NIFTP samples. The qPCR threshold cycle values for each transcript were normalized to the thyrocyte-specific thyroid peroxidase isoform 1 (TPO1) and z-transformed. Receiver operating characteristic (ROC) analyses of the composite transcript scores were used to evaluate classification of thyroid FNAs by the 5-gene isoform expression panel.

RESULTS

A molecular signature was developed by combining expression levels of specific isoforms of CDH3, FNDC4, HMGA2, KLK7, and PLAG1. FNAs containing at least 12-36 thyrocytes were sufficient for this assay. The 5-gene composite score achieved an area under the ROC curve (AUC) of 0.86 for distinguishing malignant from benign nodules, with a specificity of 91%, sensitivity of 75%, negative predictive value of 91%, and positive predictive value of 74%.

CONCLUSION

Our newly developed 5-gene isoform expression panel distinguishes benign from malignant thyroid tumors and, may help distinguish benign from malignant thyroid nodules in the context of the new NIFTP subtype.

Microfluidics for genome-wide studies involving next generation sequencing

Next-generation sequencing (NGS) has revolutionized how molecular biology studies are conducted. Its decreasing cost and increasing throughput permit profiling of genomic, transcriptomic, and epigenomic features for a wide range of applications. Microfluidics has been proven to be highly complementary to NGS technology with its unique capabilities for handling small volumes of samples and providing platforms for automation, integration, and multiplexing. In this article, we review recent progress on applying microfluidics to facilitate genome-wide studies. We emphasize on several technical aspects of NGS and how they benefit from coupling with microfluidic technology. We also summarize recent efforts on developing microfluidic technology for genomic, transcriptomic, and epigenomic studies, with emphasis on single cell analysis. We envision rapid growth in these directions, driven by the needs for testing scarce primary cell samples from patients in the context of precision medicine.

Transcriptome Profiling of In-Vivo Produced Bovine Pre-implantation Embryos Using Two-color Microarray Platform

Early embryonic loss is a large contributor to infertility in cattle. Moreover, bovine becomes an interesting model to study human preimplantation embryo development due to their similar developmental process. Although genetic factors are known to affect early embryonic development, the discovery of such factors has been a serious challenge. Microarray technology allows quantitative measurement and gene expression profiling of transcript levels on a genome-wide basis. One of the main decisions that have to be made when planning a microarray experiment is whether to use a one- or two-color approach. Two-color design increases technical replication, minimizes variability, improves sensitivity and accuracy as well as allows having loop designs, defining the common reference samples. Although microarray is a powerful biological tool, there are potential pitfalls that can attenuate its power. Hence, in this technical paper we demonstrate an optimized protocol for RNA extraction, amplification, labeling, hybridization of the labeled amplified RNA to the array, array scanning and data analysis using the two-color analysis strategy.

Protein synthetic machinery and mRNA in regenerating tips of spinal cord axons in lamprey

Polyribosomes, mRNA, and other elements of translational machinery have been reported in peripheral nerves and in elongating injured axons of sensory neurons in vitro, primarily in growth cones. Evidence for involvement of local protein synthesis in regenerating central nervous system (CNS) axons is less extensive. We monitored regeneration of back-labeled lamprey spinal axons after spinal cord transection and detected mRNA in axon tips by in situ hybridization and microaspiration of their axoplasm. Poly(A)+mRNA was present in the axon tips, and was more abundant in actively regenerating tips than in static or retracting ones. Target-specific polymerase chain reaction (PCR) and in situ hybridization revealed plentiful mRNA for the low molecular neurofilament subunit and β-tubulin, but very little for β-actin, consistent with the morphology of their tips, which lack filopodia and lamellipodia. Electron microscopy showed ribosomes/polyribosomes in the distal parts of axon tips and in association with vesicle-like membranes, primarily in the tip. In one instance, there were structures with the appearance of rough endoplasmic reticulum. Immunohistochemistry showed patches of ribosomal protein S6 positivity in a similar distribution. The results suggest that local protein synthesis might be involved in the mechanism of axon regeneration in the lamprey spinal cord. J. Comp. Neurol. 524:3614-3640, 2016. © 2016 Wiley Periodicals, Inc.

The Maternal-to-Zygotic Transition in Higher Plants: Available Approaches, Critical Limitations, and Technical Requirements

Fertilization marks the turnover from the gametophyte to sporophyte generation in higher plants. After fertilization, sporophytic development undergoes genetic turnover from maternal to zygotic control: the maternal-to-zygotic transition (MZT). The MZT is thought to be critical for early embryogenesis; however, little is known about the time course or developmental impact of the MZT in higher plants. Here, we discuss what is known in the field and focus on techniques used in relevant studies and their limitations. Some significant questions and technical requirements for further investigations are also discussed.

Global transcriptional profiling using RNA sequencing and DNA methylation patterns in highly enriched mesenchymal cells from young versus elderly women

Age-related bone loss in humans is associated with a decrease in bone formation relative to bone resorption, although the mechanisms for this impairment in bone formation with aging are not well understood. It is known that the precursors for the bone-forming osteoblasts reside in the mesenchymal cell population in bone marrow. Thus, in an effort to identify relevant genetic pathways that are altered with aging, we examined the gene expression and DNA methylation patterns from a highly enriched bone marrow mesenchymal cell population from young (mean age, 28.7 years) versus old (mean age, 73.3 years) women. Bone marrow mononuclear cells from these women were depleted of hematopoietic lineage (lin) and endothelial cells using a combination of magnetic- and fluorescence-activated cell sorting, yielding a previously characterized mesenchymal cell population (lin-/CD34-/CD31- cells) that is capable of osteoblast differentiation. Whole transcriptome RNA sequencing (RNAseq) of freshly isolated cells (without in vitro culture) identified 279 differentially expressed genes (p < 0.05, false discovery rate [q]< 0.10) between the young and old subjects. Pathway analysis revealed statistically significant (all p < 0.05) alterations in protein synthesis and degradation pathways, as well as mTOR, gap junction, calcium, melatonin and NFAT signaling pathways. Further, Reduced Representational Bisulphite sequencing (RRBS DNA methylation sequencing) revealed significant differences in methylation between the young and old subjects surrounding the promoters of 1528 target genes that also exhibited significant differences in gene expression by RNAseq. In summary, these studies provide novel insights into potential pathways affected by aging in a highly enriched human mesenchymal cell population analyzed without the confounding effects of in vitro culture. Specifically, our finding of alterations in several genes and pathways leading to impaired protein synthesis and turnover with aging in bone marrow mesenchymal cells points to the need for further studies examining how these changes, as well as the other alterations with aging that we identified, may contribute to the age-related impairment in osteoblast formation and/or function.

Transcriptional analysis of micro-dissected articular cartilage in post-traumatic murine osteoarthritis

OBJECTIVE

Identify gene changes in articular cartilage of the medial tibial plateau (MTP) at 2, 4 and 8 weeks after destabilisation of the medial meniscus (DMM) in mice. Compare our data with previously published datasets to ascertain dysregulated pathways and genes in osteoarthritis (OA).

DESIGN

RNA was extracted from the ipsilateral and contralateral MTP cartilage, amplified, labelled and hybridized on Illumina WGv2 microarrays. Results were confirmed by real-time polymerase chain reaction (PCR) for selected genes.

RESULTS

Transcriptional analysis and network reconstruction revealed changes in extracellular matrix and cytoskeletal genes induced by DMM. TGFβ signalling pathway and complement and coagulation cascade genes were regulated at 2 weeks. Fibronectin (Fn1) is a hub in a reconstructed network at 2 weeks. Regulated genes decrease over time. By 8 weeks fibromodulin (Fmod) and tenascin N (Tnn) are the only dysregulated genes present in the DMM operated knees. Comparison with human and rodent published gene sets identified genes overlapping between our array and eight other studies.

CONCLUSIONS

Cartilage contributes a minute percentage to the RNA extracted from the whole joint (<0.2%), yet is sensitive to changes in gene expression post-DMM. The post-DMM transcriptional reprogramming wanes over time dissipating by 8 weeks. Common pathways between published gene sets include focal adhesion, regulation of actin cytoskeleton and TGFβ. Common genes include Jagged 1 (Jag1), Tetraspanin 2 (Tspan2), neuroblastoma, suppression of tumourigenicity 1 (Nbl1) and N-myc downstream regulated gene 2 (Ndrg2). The concomitant genes and pathways we identify may warrant further investigation as biomarkers or modulators of OA.

MicroRNA-769-3p down-regulates NDRG1 and enhances apoptosis in MCF-7 cells during reoxygenation

Hypoxia and reoxygenation are common characteristics of solid tumors, which lead to oxidative stress and activation of stress-response genes. Previously, we observed that N-myc downstream-regulated gene 1 (NDRG1) was strongly down-regulated after shifting to reoxygenation, but the regulatory mechanism of NDRG1 remained elusive. Here we focused on the regulation of NDRG1 by microRNAs (miRNAs). Breast cancer MCF-7 cells were cultured under hypoxia for 24 h followed by 24 h of reoxygenation. The miRNA profiles were examined by Nanostring nCounter assays. Forty-three miRNAs had significant changes upon reoxygenation. In silico analysis identified four oxygen-sensitive miRNAs whose seed regions perfectly matched the 3′-UTR of NDRG1. In particular, miR-769-3p was able to inhibit the expression of NDRG1, which caused a significant reduction of NDRG1 protein upon reoxygenation. Furthermore, overexpression of miR-769-3p significantly inhibited cell proliferation and enhanced apoptosis. Our results revealed that miR-769-3p can functionally regulate NDRG1 during changes in oxygen concentration.

Single-cell RNA-seq: advances and future challenges

Phenotypically identical cells can dramatically vary with respect to behavior during their lifespan and this variation is reflected in their molecular composition such as the transcriptomic landscape. Single-cell transcriptomics using next-generation transcript sequencing (RNA-seq) is now emerging as a powerful tool to profile cell-to-cell variability on a genomic scale. Its application has already greatly impacted our conceptual understanding of diverse biological processes with broad implications for both basic and clinical research. Different single-cell RNA-seq protocols have been introduced and are reviewed here—each one with its own strengths and current limitations. We further provide an overview of the biological questions single-cell RNA-seq has been used to address, the major findings obtained from such studies, and current challenges and expected future developments in this booming field.

Isolation and characterization of human osteoblasts from needle biopsies without in vitro culture

SUMMARY

We isolate and characterize osteoblasts from humans without in vitro culture. These techniques should be broadly applicable to studying the pathogenesis of osteoporosis and other bone disorders.

INTRODUCTION

There is currently no data regarding the expression of specific genes or pathways in human osteoblasts that have not been subjected to extensive in vitro culture. Thus, we developed methods to rapidly isolate progressively enriched osteoblast populations from humans and characterized these cells.

METHODS

Needle bone biopsies of the posterior iliac crest were subjected to sequential collagenase digests. The cells from the second digest were stained with an alkaline phosphatase (AP) antibody, and the AP+ cells were isolated using magnetic cell sorting.

RESULTS

Relative to AP- cells, the AP+ cells contained virtually all of the mineralizing cells and were enriched for key osteoblast marker genes. The AP+ cells were further purified by depletion of cells expressing CD45, CD34, or CD31 (AP+/CD45/34/31- cells), which represented a highly enriched human osteoblast population devoid of hematopoietic/endothelial cells. These cells expressed osteoblast marker genes but very low to undetectable levels of SOST. We next used high-throughput RNA sequencing to compare the transcriptome of the AP+/CD45/34/31- cells to human fibroblasts and identified genes and pathways expressed only in human osteoblasts in vivo, but not in fibroblasts, including 448 genes unique to human osteoblasts.

CONCLUSIONS

We provide a detailed characterization of highly enriched human osteoblast populations without in vitro culture. These techniques should be broadly applicable to studying the pathogenesis of osteoporosis and other bone disorders.

RNA amplification for pseudogene detection using RNA-seq

Transcriptomic research using microarrays and RNA-Sequencing (RNA-seq) is now possible starting from minute biological samples, such as clinical specimens or embryos, due to the development of highly sensitive and reproducible cDNA synthesis methods. Here, we describe a quick method of RNA amplification and double-stranded cDNA synthesis starting with 10 ng of high-quality total RNA extracted from porcine embryos. The final product (double-stranded DNA) is adequate for the detection by RNA-seq of protein-coding transcripts, as well as of all the other classes of noncoding RNAs, including pseudogenes.

Genome-wide comparison of two RNA-stabilizing reagents for transcriptional profiling of peripheral blood

Peripheral whole blood is relatively easily obtained for monitoring gene expression for biomarker discovery using transcriptomic platforms such as genome-wide microarrays. However, whole blood provides challenges caused by sensitivity for ex vivo incubation and overrepresentation of globin mRNAs. We compared the performance of 2 commercial whole blood preservation methods, TEMPUS (Applied Biosystems, Foster City, CA) and PAXgene (PreAnalytiX, Qiagen BD, Valencia, CA), using 2 RNA amplification protocols and high-density microarrays. Performance of commercial globin mRNA reduction protocol also was studied. Human peripheral blood samples collected with TEMPUS and PAXgene Blood RNA tubes were amplified with the RiboAmp OA 1 Round RNA Amplification Kit (Arcturus; Applied Biosystems) and the Affymetrix (Santa Clara, CA) small sample protocol. Affymetrix globin reduction protocol was applied for total RNA samples. Samples amplified with RiboAmp were hybridized on Illumina Sentrix HumanRef-8 Expression BeadChips (Illumina Inc, San Diego, CA) and subjected to statistical analyses. RiboAmp mRNA amplification did not notably amplify globin mRNA that is overrepresented in RNA isolated by both TEMPUS and PAXgene preservation. Enzymatic depletion of globin transcript reduced the quality of total RNA and is thus not recommendable. Microarray analysis showed acceptable correlation within and between the RNA preservation methods, but altogether 443 transcripts were differentially expressed between RNA samples preserved in TEMPUS and PAXgene tubes. We demonstrated that the 2 tested blood RNA-preservation methods combined with RiboAmp mRNA amplification may be used for microarray experiments without the need for a prior globin RNA reduction. However, because genes involved in immune cell functions and gene regulatory pathways were differentially expressed as a result of the technical bias between the preservation methods, they should not be used in the same analytic setting.

Adjustment method for microarray data generated using two-cycle RNA labeling protocol

Background

Microarray technology is widely utilized for monitoring the expression changes of thousands of genes simultaneously. However, the requirement of relatively large amount of RNA for labeling and hybridization makes it difficult to perform microarray experiments with limited biological materials, thus leads to the development of many methods for preparing and amplifying mRNA. It is addressed that amplification methods usually bring bias, which may strongly hamper the following interpretation of the results. A big challenge is how to correct for the bias before further analysis.

Results

In this article, we observed the bias in rice gene expression microarray data generated with the Affymetrix one-cycle, two-cycle RNA labeling protocols, followed by validation with Real Time PCR. Based on these data, we proposed a statistical framework to model the processes of mRNA two-cycle linear amplification, and established a linear model for probe level correction. Maximum Likelihood Estimation (MLE) was applied to perform robust estimation of the Retaining Rate for each probe. After bias correction, some known pre-processing methods, such as PDNN, could be combined to finish preprocessing. Then, we evaluated our model and the results suggest that our model can effectively increase the quality of the microarray raw data: (i) Decrease the Coefficient of Variation for PM intensities of probe sets; (ii) Distinguish the microarray samples of five stages for rice stamen development more clearly; (iii) Improve the correlation coefficients among stamen microarray samples. We also discussed the necessity of model adjustment by comparing with another simple adjustment method.

Conclusion

We conclude that the adjustment model is necessary and could effectively increase the quality of estimation for gene expression from the microarray raw data.

Extraction of bacterial RNA from soil: challenges and solutions

Detection of bacterial gene expression in soil emerged in the early 1990s and provided information on bacterial responses in their original soil environments. As a key procedure in the detection, extraction of bacterial RNA from soil has attracted much interest, and many methods of soil RNA extraction have been reported in the past 20 years. In addition to various RT-PCR-based technologies, new technologies for gene expression analysis, such as microarrays and high-throughput sequencing technologies, have recently been applied to examine bacterial gene expression in soil. These technologies are driving improvements in RNA extraction protocols. In this mini-review, progress in the extraction of bacterial RNA from soil is summarized with emphasis on the major difficulties in the development of methodologies and corresponding strategies to overcome them.

Characterization of in vitro transcription amplification linearity and variability in the low copy number regime using External RNA Control Consortium (ERCC) spike-ins

Using spike-in controls designed to mimic mammalian mRNA species, we used the quantitative reverse transcription polymerase chain reaction (RT-qPCR) to assess the performance of in vitro transcription (IVT) amplification process of small samples. We focused especially on the confidence of the transcript level measurement, which is essential for differential gene expression analyses. IVT reproduced gene expression profiles down to approximately 100 absolute input copies. However, a RT-qPCR analysis of the antisense RNA showed a systematic bias against low copy number transcripts, regardless of sequence. Experiments also showed that noise increases with decreasing copy number. First-round IVT preserved the gene expression information within a sample down to the 100 copy level, regardless of total input sample amount. However, the amplification was nonlinear under low total RNA input/long IVT conditions. Variability of the amplification increased predictably with decreasing input copy number. For the small enrichments of interest in typical differential gene expression studies (e.g., twofold changes), the bias from IVT reactions is unlikely to affect the results. In limited cases, some transcript-specific differential gene expression values will need adjustment to reflect this bias. Proper experimental design with reasonable detection limits will yield differential gene expression capability even between low copy number transcripts.

Quantitative single-cell ion-channel gene expression profiling through an improved qRT-PCR technique combined with whole cell patch clamp

Cellular excitability originates from a concerted action of different ion channels. The genomic diversity of ion channels (over 100 different genes) underlies the functional diversity of neurons in the central nervous system (CNS) and even within a specific type of neurons large differences in channel expression have been observed. Patch-clamp is a powerful technique to study the electrophysiology of excitability at the single cell level, allowing exploration of cell-to-cell variability. Only a few attempts have been made to link electrophysiological profiling to mRNA transcript levels and most suffered from experimental noise precluding conclusive quantitative correlations. Here we describe a refinement to the technique that combines patch-clamp analysis with quantitative real-time (qRT) PCR at the single cell level. Hereto the expression of a housekeeping gene was used to normalize for cell-to-cell variability in mRNA isolation and the subsequent processing steps for performing qRT-PCR. However, the mRNA yield from a single cell was insufficient for performing a valid qRT-PCR assay; this was resolved by including a RNA amplification step. The technique was validated on a stable Ltk(-) cell line expressing the Kv2.1 channel and on embryonic dorsal root ganglion (DRG) cells probing for the expression of Kv2.1. Current density and transcript quantity displayed a clear correlation when the qRT-PCR assay was done in twofold and the data normalized to the transcript level of the housekeeping gene GAPD. Without this normalization no significant correlation was obtained. This improved technique should prove very valuable for studying the molecular background of diversity in cellular excitability.

Options available for labelling nucleic acid samples in DNA microarray-based detection methods

DNA microarrays are considered by many researchers to be the platform of choice for the high-throughput analysis of nucleic acids. Since the past two decades, they have been used constantly as powerful tools in differential gene expression, SNP genotyping, DNA sequencing, gene discovery, disease diagnostic and pathways reconstruction. Several methods have been developed to enable samples of limited amounts of RNA to be quantified. Here we evaluate classical and up-to-date assays made available for labelling those samples. This review also sheds light on the recently developed strategies that ensure high sensitivity such as sample and signal amplification, quantum dot, surface plasmom resonance, nanoparticles and cationinc polythiophenes.

Systems biology of fungal infection

Elucidation of pathogenicity mechanisms of the most important human-pathogenic fungi, Aspergillus fumigatus and Candida albicans, has gained great interest in the light of the steadily increasing number of cases of invasive fungal infections. A key feature of these infections is the interaction of the different fungal morphotypes with epithelial and immune effector cells in the human host. Because of the high level of complexity, it is necessary to describe and understand invasive fungal infection by taking a systems biological approach, i.e., by a comprehensive quantitative analysis of the non-linear and selective interactions of a large number of functionally diverse, and frequently multifunctional, sets of elements, e.g., genes, proteins, metabolites, which produce coherent and emergent behaviors in time and space. The recent advances in systems biology will now make it possible to uncover the structure and dynamics of molecular and cellular cause-effect relationships within these pathogenic interactions. We review current efforts to integrate omics and image-based data of host-pathogen interactions into network and spatio-temporal models. The modeling will help to elucidate pathogenicity mechanisms and to identify diagnostic biomarkers and potential drug targets for therapy and could thus pave the way for novel intervention strategies based on novel antifungal drugs and cell therapy.

Characterization of mesenchymal progenitor cells isolated from human bone marrow by negative selection

Studies on the pathogenesis of osteoporosis and other metabolic bone diseases would be greatly facilitated by the development of approaches to assess changes in gene expression in osteoblast/osteoprogenitor populations in vivo without the potentially confounding effects of in vitro culture and expansion of the cells. While positive selection to identify a progenitor population in human marrow can be used to select for cells capable of osteoblast differentiation, each of the markers that have been used to identify marrow mesenchymal populations (alkaline phosphatase [AP], Stro-1, CD29, CD49a, CD73, CD90, CD105, CD166, CD44, CD146 and CD271) may be expressed on distinct subsets of marrow mesenchymal cells. Thus, positive selection with one or more of these markers could exclude a possibly relevant cell population that may undergo important changes in various clinical conditions. In the present report, we describe the isolation and characterization of human osteoprogenitor cells obtained by depletion of bone marrow cells of all hematopoietic lineage/hematopoietic stem cells and endothelial/endothelial precursor cells (lin-/CD34/CD31-). The yield of lin-/CD34/CD31- cells from ~10 mL of bone marrow (~80 million mononuclear cells) was ~80,000 cells (0.1% of mononuclear cells). While not selected on the basis of expression for the mesenchymal marker, Stro-1, 68% of these cells were Stro-1+. Using linear whole transcriptome amplification followed by quantitative polymerase chain reaction (QPCR) analysis, we also demonstrated that, compared to lin- cells (which are already depleted of hematopoietic cells), lin-/CD34/31- cells expressed markedly lower mRNA levels for the endothelial/hematopoietic markers, CD34, CD31, CD45, and CD133. Lin-/CD34/31- cells were also enriched for the expression of mesenchymal/osteoblastic markers, with a further increase in runx2, osterix, and AP mRNA expression following in vitro culture under osteogenic conditions. Importantly, lin-/CD34/31- cells contained virtually all of the mineralizing cells in human marrow: while these cells displayed robust calcium deposition in vitro, lin-/CD34/31+ cells demonstrated little or no mineralization when cultured under identical osteogenic conditions. Lin-/CD34/31- cells thus represent a human bone marrow population highly enriched for mesenchymal/osteoblast progenitor cells that can be analyzed without in vitro culture in various metabolic bone disorders, including osteoporosis and aging.

Use of non-amplified RNA samples for microarray analysis of gene expression

Demand for high quality gene expression data has driven the development of revolutionary microarray technologies. The quality of the data is affected by the performance of the microarray platform as well as how the nucleic acid targets are prepared. The most common method for target nucleic acid preparation includes in vitro transcription amplification of the sample RNA. Although this method requires a small amount of starting material and is reported to have high reproducibility, there are also technical disadvantages such as amplification bias and the long, laborious protocol. Using RNA derived from human brain, breast and colon, we demonstrate that a non-amplification method, which was previously shown to be inferior, could be transformed to a highly quantitative method with a dynamic range of five orders of magnitude. Furthermore, the correlation coefficient calculated by comparing microarray assays using non-amplified samples with qRT-PCR assays was approximately 0.9, a value much higher than when samples were prepared using amplification methods. Our results were also compared with data from various microarray platforms studied in the MicroArray Quality Control (MAQC) project. In combination with micro-columnar 3D-Gene™ microarray, this non-amplification method is applicable to a variety of genetic analyses, including biomarker screening and diagnostic tests for cancer.

Highly parallel genome-wide expression analysis of single mammalian cells

BACKGROUND

We have developed a high-throughput amplification method for generating robust gene expression profiles using single cell or low RNA inputs.

METHODOLOGY/PRINCIPAL FINDINGS

The method uses tagged priming and template-switching, resulting in the incorporation of universal PCR priming sites at both ends of the synthesized cDNA for global PCR amplification. Coupled with a whole-genome gene expression microarray platform, we routinely obtain expression correlation values of R(2)~0.76-0.80 between individual cells and R(2)~0.69 between 50 pg total RNA replicates. Expression profiles generated from single cells or 50 pg total RNA correlate well with that generated with higher input (1 ng total RNA) (R(2)~0.80). Also, the assay is sufficiently sensitive to detect, in a single cell, approximately 63% of the number of genes detected with 1 ng input, with approximately 97% of the genes detected in the single-cell input also detected in the higher input.

CONCLUSIONS/SIGNIFICANCE

In summary, our method facilitates whole-genome gene expression profiling in contexts where starting material is extremely limiting, particularly in areas such as the study of progenitor cells in early development and tumor stem cell biology.

Single-cell and regional gene expression analysis in Alzheimer's disease

The clinical manifestations of Alzheimer's disease (AD) are secondary to the substantial loss of cortical neurons. To be effective, neuroprotective strategies will need to target the primary pathogenic mechanisms of AD prior to cell loss. The differences between neurons are largely determined by their specific repertoire of mRNAs. Thus, transcriptomic analyses that do not assume a priori etiological hypotheses are potentially powerful tools that can be used to understand the pathogenesis of complex diseases, including AD. The human brain comprises thousands of different cell types of both neuronal and non-neuronal origins. Information about individual cell-type-specific gene expression patterns will allow for a better understanding of the mechanisms that govern the progression of AD, which may lead to new therapeutic targets for prevention and treatment of the disease. This review provides an overview of the current technologies in use and the developments for single-cell extraction and transcriptome analysis. Recent transcriptome profiling studies on individual AD-afflicted brain cells are also discussed.

Laser capture microdissection of Candida albicans from host tissue

Laser microdissection is a technique in which specific populations of cells are acquired from sections of complex tissue under direct microscopic visualization. The technique can be used to selectively harvest or ablate host and/or fungal cells from a variety of biological specimens, including human, animal, or plant tissue sections. When coupled with downstream applications such as proteomic and molecular analyses, laser microdissection can address a variety of important biological questions specifically related to the in vivo host-fungus interaction. In this chapter, we describe how laser microdissection enables researchers to selectively isolate Candida albicans cells from host-infected tissue. Detailed protocols are provided for tissue handling and processing, slide preparation, and laser capture microdissection (LCM). Using these methods, we highlight the use of LCM to examine infection-related C. albicans gene expression.

Methyl-binding domain protein-based DNA isolation from human blood serum combines DNA analyses and serum-autoantibody testing

Background

Circulating cell free DNA in serum as well as serum-autoantibodies and the serum proteome have great potential to contribute to early cancer diagnostics via non invasive blood tests. However, most DNA preparation protocols destroy the protein fraction and therefore do not allow subsequent protein analyses. In this study a novel approach based on methyl binding domain protein (MBD) is described to overcome the technical difficulties of combining DNA and protein analysis out of one single serum sample.

Methods

Serum or plasma samples from 98 control individuals and 54 breast cancer patients were evaluated upon silica membrane- or MBD affinity-based DNA isolation via qPCR targeting potential DNA methylation markers as well as by protein-microarrays for tumor-autoantibody testing.

Results

In control individuals, an average DNA level of 22.8 ± 25.7 ng/ml was detected applying the silica membrane based protocol and 8.5 ± 7.5 ng/ml using the MBD-approach, both values strongly dependent on the serum sample preparation methods used. In contrast to malignant and benign tumor serum samples, cell free DNA concentrations were significantly elevated in sera of metastasizing breast cancer patients. Technical evaluation revealed that serum upon MBD-based DNA isolation is suitable for protein-array analyses when data are consistent to untreated serum samples.

Conclusion

MBD affinity purification allows DNA isolations under native conditions retaining the protein function, thus for example enabling combined analyses of DNA methylation and autoantigene-profiles from the same serum sample and thereby improving minimal invasive diagnostics.

Effects of intermittent parathyroid hormone treatment on osteoprogenitor cells in postmenopausal women

Intermittent parathyroid hormone (PTH) 1-34 treatment stimulates bone formation, but the molecular mechanisms mediating this effect have not been previously studied in humans. Thus, we used magnetic activated cell sorting to isolate hematopoietic lineage negative (lin-)/alkaline phosphatase positive (AP+) osteoprogenitor cells from bone marrow of 20 postmenopausal women treated with PTH (1-34) for 14 days and 19 control subjects. Serum PINP and CTX increased in PTH-treated subjects (by 97% and 30%, respectively, P<0.001). Bone marrow lin-/AP+ cells from PTH-treated subjects showed an increase in the RANKL/OPG mRNA ratio (by 7.5-fold, P=0.011) and in the mRNAs for c-fos (a known PTH-responsive gene, by 42%, P=0.035) and VEGF-C (by 57%, P=0.046). Gene Set Enrichment Analysis (GSEA, testing for changes in pre-specified pathways) demonstrated that PTH had no effect on osteoblast proliferation, apoptosis, or differentiation markers. However, PTH treatment resulted in a significant decrease (GSEA P-value, 0.005) in a panel of BMP target genes in the lin-/AP+ cells. Our findings thus identify several future directions for studying mechanisms of PTH action in humans. First, given the increasing evidence that PTH induces angiogenesis, the role of increased VEGF-C production by bone marrow osteoprogenitor cells in mediating this effect and the anabolic response to PTH warrants further study. Second, while the observed inhibition of BMP target gene expression by PTH is not consistent with the anabolic effects of PTH on bone and requires further validation, these data do generate the hypothesis that an inhibition of BMP signaling by PTH may, over time, limit the availability of mature osteoblasts on bone surfaces and thereby contribute to the observed waning of the anabolic response to PTH.

Effects of estrogen on osteoprogenitor cells and cytokines/bone-regulatory factors in postmenopausal women

Decreases in estrogen levels contribute not only to early postmenopausal bone loss but also to bone loss with aging. While estrogen is critical for the maintenance of bone formation, the mechanism(s) of this effect remain unclear. Thus, we assessed the effects of 4months of transdermal estradiol treatment (0.05mg/day) of postmenopausal women as compared to no treatment (n=16 per group) on the expression of genes in pre-specified pathways in freshly isolated bone marrow osteoprogenitor cells (hematopoietic lineage [lin]-/Stro1+). We also evaluated whether estrogen treatment modulated peripheral blood or bone marrow plasma levels of the Wnt antagonists, sclerostin and DKK1, as well as serotonin, OPG, RANKL, adiponectin, oxytocin, and inflammatory cytokines (TNFα, IL-1β, and IL-6), as each of these molecules have recently been shown to play an important role in regulating osteoblast function and/or being responsive to estrogen. We observed a significant decrease in the expression of several proliferation markers (cyclin B1, cyclin E1, E2F1) and increase in adhesion molecules (N-cadherin) in bone marrow lin-/Stro1+ cells from estrogen-treated compared to control women. None of the peripheral blood or bone marrow plasma marker levels differed between the two groups, with the exception of sclerostin levels, which were significantly lower in the estrogen-treated as compared to the control women in peripheral serum (by 32%, P=0.009) and in bone marrow plasma (by 34%, P=0.017). There were significant differences in bone marrow versus peripheral plasma levels of several factors: sclerostin and OPG levels were higher in bone marrow as compared to peripheral plasma, whereas serotonin and adiponectin levels were higher in peripheral as compared to bone marrow plasma. In summary, our data directly assessing possible regulation by estrogen of osteoprogenitor cells in humans indicate that, consistent with previous studies in mice, estrogen suppresses the proliferation of human bone marrow lin-/Stro1+ cells, which likely represent early osteoprogenitor cells. Further animal and human studies are needed to define the role of the changes we observed in mRNAs for adhesion molecules in these cells and in local sclerostin production in bone in mediating the effects of estrogen on bone metabolism in humans.

Oligonucleotide microarray expression profiling of contrasting invasive phenotypes in colorectal cancer

This chapter refers to the application of laser-capture microdissection with oligonucleotide microarray analysis. The protocol described has been successfully used to identify differential transcript expression between contrasting colorectal cancer invasive phenotypes. Tissue processing, RNA extraction, quality control, amplification, fluorescent labelling, purification, hybridisation, and elements of data analysis are covered.

Functional genomics: new insights into the 'function' of low levels of gene expression in stem cells

Understanding the global gene expression profile of stem cells and their multilineage differentiation will be essential for their ultimate therapeutic application. Efforts to characterize stem cells have relied on analyzing the genome-wide expression profiles that are biased towards the identification of genes that display the most pronounced differential expression. Rather than being viewed as a “blank” state, recent studies suggest that stem cells express low levels of multiple lineage specific genes prior to differentiation, a phenomenon known as “lineage priming.” It is not likely that low levels of lineage-specific genes produce sufficient amounts of differentiation factors, but rather to provide rapid transcription to a wide range of lineage programs prior to differentiation. Thus, stem cell differentiation may involve the elimination of other potential pathways and the activation of a specific lineage program.

Microgenomic analysis in skeletal muscle: expression signatures of individual fast and slow myofibers

Background

Skeletal muscle is a complex, versatile tissue composed of a variety of functionally diverse fiber types. Although the biochemical, structural and functional properties of myofibers have been the subject of intense investigation for the last decades, understanding molecular processes regulating fiber type diversity is still complicated by the heterogeneity of cell types present in the whole muscle organ.

Methodology/Principal Findings

We have produced a first catalogue of genes expressed in mouse slow-oxidative (type 1) and fast-glycolytic (type 2B) fibers through transcriptome analysis at the single fiber level (microgenomics). Individual fibers were obtained from murine soleus and EDL muscles and initially classified by myosin heavy chain isoform content. Gene expression profiling on high density DNA oligonucleotide microarrays showed that both qualitative and quantitative improvements were achieved, compared to results with standard muscle homogenate. First, myofiber profiles were virtually free from non-muscle transcriptional activity. Second, thousands of muscle-specific genes were identified, leading to a better definition of gene signatures in the two fiber types as well as the detection of metabolic and signaling pathways that are differentially activated in specific fiber types. Several regulatory proteins showed preferential expression in slow myofibers. Discriminant analysis revealed novel genes that could be useful for fiber type functional classification.

Conclusions/Significance

As gene expression analyses at the single fiber level significantly increased the resolution power, this innovative approach would allow a better understanding of the adaptive transcriptomic transitions occurring in myofibers under physiological and pathological conditions.

Early in vitro transcription termination in human H5 influenza viral RNA synthesis

Rapid diagnostic identification of the human H5 influenza virus is a strategic cornerstone for outbreak prevention. We recently reported a method for direct detection of viral RNA from a highly pathogenic human H5 influenza strain (A/Hanoi/30408/2005(H5N1)), which necessarily was transcribed in vitro from non-viral sources. This article provides an in-depth analysis of the reaction conditions for in vitro transcription (IVT) of full-length influenza H5 RNA, which is needed for diagnostic RNA production, for the T7 and SP6 phage promoter systems. Gel analysis of RNA transcribed from plasmids containing the H5 sequence between a 5' SP6 promoter and 3' restriction site (BsmBI) showed that three sequence-verified bands at 1,776, 784, and 591 bases were consistently produced, whereas only one 1,776-base band was expected. These fragments were not observed in H1 or H3 influenza RNA transcribed under similar conditions. A reverse complement of the sequence produced only a single band at 1,776 bases, which suggested either self-cleavage or early termination. Aliquots of the IVT reaction were quenched with EDTA to track the generation of the bands over time, which maintained a constant concentration ratio. The H5 sequence was cloned with T7 and SP6 RNA polymerase promoters to allow transcription in either direction with either polymerase. The T7 transcription product from purified, restricted plasmids in the vRNA direction only produced the 1,776-base full-length sequence and the 784-base fragment, instead of the three bands generated by the SP6 system, suggesting an early termination mechanism. Additionally, the T7 system produced a higher fraction of full-length vRNA transcripts than the SP6 system did under similar reaction conditions. By sequencing we identified a type II RNA hairpin loop terminator, which forms in a transcription direction-dependent fashion. Variation of the magnesium concentration produced the greatest impact on termination profiles, where some reaction mixtures were unable to produce full-length transcripts. Optimized conditions are presented for the T7 and SP6 phage polymerase systems to minimize these early termination events during in vitro transcription of H5 influenza vRNA.

RNA sequencing: advances, challenges and opportunities

In the few years since its initial application, massively parallel cDNA sequencing, or RNA-seq, has allowed many advances in the characterization and quantification of transcriptomes. Recently, several developments in RNA-seq methods have provided an even more complete characterization of RNA transcripts. These developments include improvements in transcription start site mapping, strand-specific measurements, gene fusion detection, small RNA characterization and detection of alternative splicing events. Ongoing developments promise further advances in the application of RNA-seq, particularly direct RNA sequencing and approaches that allow RNA quantification from very small amounts of cellular materials.

Accurate expression profiling of very small cell populations

BACKGROUND

Expression profiling, the measurement of all transcripts of a cell or tissue type, is currently the most comprehensive method to describe their physiological states. Given that accurate profiling methods currently available require RNA amounts found in thousands to millions of cells, many fields of biology working with specialized cell types cannot use these techniques because available cell numbers are limited. Currently available alternative methods for expression profiling from nanograms of RNA or from very small cell populations lack a broad validation of results to provide accurate information about the measured transcripts.

METHODS AND FINDINGS

We provide evidence that currently available methods for expression profiling of very small cell populations are prone to technical noise and therefore cannot be used efficiently as discovery tools. Furthermore, we present Pico Profiling, a new expression profiling method from as few as ten cells, and we show that this approach is as informative as standard techniques from thousands to millions of cells. The central component of Pico Profiling is Whole Transcriptome Amplification (WTA), which generates expression profiles that are highly comparable to those produced by others, at different times, by standard protocols or by Real-time PCR. We provide a complete workflow from RNA isolation to analysis of expression profiles.

CONCLUSIONS

Pico Profiling, as presented here, allows generating an accurate expression profile from cell populations as small as ten cells.

Diversity of phytoplankton nitrate transporter sequences from isolated single cells and mixed samples from the East China Sea and mRNA quantification

The transcript abundances of nitrate transporter genes (Nrt2) were proposed as potential markers for nitrogen deficiency in marine diatoms. To correctly quantify diatom Nrt2 mRNA in the East China Sea (ECS), we utilized both mixed-species sequencing and single-cell PCR to expand the sequence database for this region. Using the single-cell method of PCR, 9 new diatom Nrt2 sequences belonging to 5 genera, the Nrt2 sequences of which have never been reported before, were obtained. On the other hand, 291 sequences homologous to Nrt2 were retrieved from mixed-species sequencing using degenerate primers, and these sequences were clustered into 12 major groups according to a phylogenetic analysis. Based on sequence alignments, 11 pairs of group-specific PCR primers were designed to detect Nrt2 mRNA levels in the ECS, and 3 of these primer pairs showed high specificity to target species. In ECS phytoplankton samples, environmental RNA was amplified via antisense RNA amplification followed by cDNA production. Subsequently, Nrt2 transcript levels were readily detected using quantitative PCR. Our results indicated that investigating sequence diversity followed by careful primer design and evaluation is a good strategy to quantify the expression of genes of ecologically important phytoplankton.

Effects of chronic estrogen treatment on modulating age-related bone loss in female mice

While female mice do not have the equivalent of a menopause, they do undergo reproductive senescence. Thus, to dissociate the effects of aging versus estrogen deficiency on age-related bone loss, we sham-operated, ovariectomized, or ovariectomized and estrogen-replaced female C57/BL6 mice at 6 months of age and followed them to age 18 to 22 months. Lumbar spines and femurs were excised for analysis, and bone marrow hematopoietic lineage negative (lin-) cells (enriched for osteoprogenitor cells) were isolated for gene expression studies. Six-month-old intact control mice were euthanized to define baseline parameters. Compared with young mice, aged/sham-operated mice had a 42% reduction in lumbar spine bone volume/total volume (BV/TV), and maintaining constant estrogen levels over life in ovariectomized/estrogen-treated mice did not prevent age-related trabecular bone loss at this site. By contrast, lifelong estrogen treatment of ovariectomized mice completely prevented the age-related reduction in cortical volumetric bone mineral density (vBMD) and thickness at the tibial diaphysis present in the aged/sham-operated mice. As compared with cells from young mice, lin- cells from aged/sham-operated mice expressed significantly higher mRNA levels for osteoblast differentiation and proliferation marker genes. These data thus demonstrate that, in mice, age-related loss of cortical bone in the appendicular skeleton, but not loss of trabecular bone in the spine, can be prevented by maintaining constant estrogen levels over life. The observed increase in osteoblastic differentiation and proliferation marker gene expression in progenitor bone marrow cells from aged versus young mice may represent a compensatory mechanism in response to ongoing bone loss.

Circulating osteogenic cells: characterization and relationship to rates of bone loss in postmenopausal women

There is increasing evidence that osteogenic cells are present not only in bone marrow (BM) but also in peripheral blood (PB). Since staining for alkaline phosphatase (AP) identifies osteoprogenitor cells in BM, we sought to further characterize BM versus PB hematopoietic lineage negative (lin-)/AP+ cells and to compare gene expression in PB lin-/AP+ cells from postmenopausal women undergoing rapid versus slow bone loss. PB lin-/AP+ cells were smaller than their BM counterparts, and both were negative for the pan-hematopoietic marker, CD45. BM and PB lin-/AP+ cells were capable of mineralization in vitro. Using whole genome linear amplification followed by quantitative polymerase chain reaction (QPCR) analysis, we found that relative to the BM cells, PB lin-/AP+ cells expressed similar levels of a number of key osteoblast marker genes (runx2, osterix, osteopontin, OPG, periostin), consistent with the PB cells being in the osteoblastic lineage. Importantly, however, compared to the BM cells, PB lin-/AP+ cells expressed lower levels of mRNAs for AP, type I collagen, and for a panel of proliferation markers, but higher levels of osteocalcin, osteonectin, and PTHR1 mRNAs, as well as those for RANKL and ICAM-1, both of which are important in supporting osteoclastogenesis. Using microarray followed by QPCR analysis, we further demonstrated that, compared to postmenopausal women undergoing slow bone loss, PB lin-/AP+ cells from women undergoing rapid bone loss expressed lower levels of mRNAs for hydroxyprostaglandin dehydrogenase, interferon regulator factor 3, Wnt1-induced secreted protein 1, and TGFbeta2, but higher levels of the Smad3 interacting protein, zinc finger DHHC-type containing 4 and col1alpha2. These data thus demonstrate that while PB lin-/AP+ cells express a number of osteoblastic genes and are capable of mineralization, they are a relatively quiescent cell population, both in terms of cell proliferation and matrix synthesis. However, their higher expression of RANKL and ICAM-1 mRNAs as compared to BM lin-/AP+ cells suggests a role for the PB lin-/AP+ cells in regulating osteoclastogenesis that warrants further investigation. Our study also provides "proof-of-concept" for the use of PB lin-/AP+ cells in clinical-investigative studies, and identifies several pathways that could potentially regulate rates of bone loss in postmenopausal women.

An optimized RNA amplification method for prokaryotic expression profiling analysis

DNA microarray technology has been extensively used for gene expression analysis of both eukaryotic and prokaryotic organisms. For eukaryotic gene expression profiling, the poly(A)-based reverse transcription of messenger RNA (mRNA) followed by T7 RNA polymerase-based in vitro transcription is generally required to produce enough RNA targets for hybridization with the microarray chips. However, the same method cannot be directly applied to prokaryotic mRNAs due to the lack of poly(A) sequences at the 3' ends. Conventional methods usually require large amounts of starting RNAs and lead to high background noise. Recently developed amplification methods enable smaller amounts of prokaryotic RNA to be used from samples with species-specific primers, oligo(dT) primers, or random primers. In this study, three target preparation methods, including the direct labeling, polyadenylation-involved oligo-dT priming, and random priming amplification (respectively referred to as DL, PAOD, and RPA hereafter) were evaluated through expression profiling of a heat shock model of Escherichia coli. The PAOD method was found to be more sensitive and more specific in differential gene expression measurements than either DL and RPA, even when the E. coli RNA was only a small proportion of the simulated eukaryotic host RNA. The results suggest that PAOD is the preferred target preparation method for prokaryotic transcriptome.

Metagenomics for the discovery of novel human viruses

Modern laboratory techniques for the detection of novel human viruses are greatly needed as physicians and epidemiologists increasingly deal with infectious diseases caused by new or previously unrecognized pathogens. There are many clinical syndromes in which viruses are suspected to play a role, but for which traditional microbiology techniques routinely fail in uncovering the etiologic agent. In addition, new viruses continue to challenge the human population owing to the encroachment of human settlements into animal and livestock habitats, globalization, climate change, growing numbers of immunocompromised people and bioterrorism. Metagenomics-based tools, such as microarrays and high-throughput sequencing are ideal for responding to these challenges. Pan-viral microarrays, containing representative sequences from all known viruses, have been used to detect novel and distantly-related variants of known viruses. Sequencing-based methods have also been successfully employed to detect novel viruses and have the potential to detect the full spectrum of viruses, including those present in low numbers.

Standards affecting the consistency of gene expression arrays in clinical applications

The use of microarray technology to measure gene expression has created optimism for the feasibility of using molecular assessments of tumors routinely in the clinical management of cancer. Gene expression arrays have been pioneers in the development of standards; both for research use and now for clinical application. Some of the existing standards have been driven by the early perception that microarray technology was inconsistent and perhaps unreliable. More recent experimentation has shown that reproducible data can be achieved and clinical standards are beginning to emerge. For the transcriptional assessment of tumors, this means a system that correctly samples a tumor, isolates RNA and processes this for microarray analysis, evaluates the data, and communicates findings in a consistent and timely fashion. The most important standard is to show that a clinically important assessment can be made with microarray data. The standards emerging from work on various parts of the entire process could guide the development of a workable system. However, the final standard for each component of the process depends on the accuracy required when the assay becomes part of the clinical routine: a routine that now includes the molecular evaluation of tumors.

Statistical evaluation of transcriptomic data generated using the Affymetrix one-cycle, two-cycle and IVT-Express RNA labelling protocols with the Arabidopsis ATH1 microarray

BACKGROUND

Microarrays are a powerful tool used for the determination of global RNA expression. There is an increasing requirement to focus on profiling gene expression in tissues where it is difficult to obtain large quantities of material, for example individual tissues within organs such as the root, or individual isolated cells. From such samples, it is difficult to produce the amount of RNA required for labelling and hybridisation in microarray experiments, thus a process of amplification is usually adopted. Despite the increasing use of two-cycle amplification for transcriptomic analyses on the Affymetrix ATH1 array, there has been no report investigating any potential bias in gene representation that may occur as a result.

RESULTS

Here we compare transcriptomic data generated using Affymetrix one-cycle (standard labelling protocol), two-cycle (small-sample protocol) and IVT-Express protocols with the Affymetrix ATH1 array using Arabidopsis root samples. Results obtained with each protocol are broadly similar. However, we show that there are 35 probe sets (of a total of 22810) that are misrepresented in the two-cycle data sets. Of these, 33 probe sets were classed as mis-amplified when comparisons of two independent publicly available data sets were undertaken.

CONCLUSIONS

Given the unreliable nature of the highlighted probes, we caution against using data associated with the corresponding genes in analyses involving transcriptomic data generated with two-cycle amplification protocols. We have shown that the Affymetrix IVT-E labelling protocol produces data with less associated bias than the two-cycle protocol, and as such, would recommend this kit for new experiments that involve small samples.

Comparative gene expression profiling between human cultured myotubes and skeletal muscle tissue

Background

A high-sensitivity DNA microarray platform requiring nanograms of RNA input facilitates the application of transcriptome analysis to individual skeletal muscle (SM) tissue samples. Culturing myotubes from SM-biopsies enables investigating transcriptional defects and assaying therapeutic strategies. This study compares the transcriptome of aneurally cultured human SM cells versus that of tissue biopsies.

Results

We used the Illumina expression BeadChips to determine the transcriptomic differences between tissue and cultured SM samples from five individuals. Changes in the expression of several genes were confirmed by QuantiGene Plex assay or reverse transcription real-time PCR. In cultured myotubes compared to the tissue, 1216 genes were regulated: 583 down and 633 up. Gene ontology analysis showed that downregulated genes were mainly associated with cytoplasm, particularly mitochondria, and involved in metabolism and the muscle-system/contraction process. Upregulated genes were predominantly related to cytoplasm, endoplasmic reticulum, and extracellular matrix. The most significantly regulated pathway was mitochondrial dysfunction. Apoptosis genes were also modulated. Among the most downregulated genes detected in this study were genes encoding metabolic proteins AMPD1, PYGM, CPT1B and UCP3, muscle-system proteins TMOD4, MYBPC1, MYOZ1 and XIRP2, the proteolytic CAPN3 and the myogenic regulator MYF6. Coordinated reduced expression of five members of the GIMAP gene family, which form a cluster on chromosome 7, was shown, and the GIMAP4-reduction was validated. Within the most upregulated group were genes encoding senescence/apoptosis-related proteins CDKN1A and KIAA1199 and potential regulatory factors HIF1A, TOP2A and CCDC80.

Conclusions

Cultured muscle cells display reductive metabolic and muscle-system transcriptome adaptations as observed in muscle atrophy and they activate tissue-remodeling and senescence/apoptosis processes.

Challenges in defining predictive markers for response to endocrine therapy in breast cancer

Endocrine therapy is a major treatment modality for hormone-dependent breast cancer. It has a relatively low morbidity, and there is evidence that antihormonal treatments have had a significant effect in reducing mortality for breast cancer. Despite this, resistance to endocrine therapy, either primary or acquired during treatment, occurs in the majority of patients, and is a major obstacle to optimal clinical management. There is therefore an urgent need to identify, on an individual basis, those tumors that are most likely to respond to endocrine therapy (so sparing patients with resistant tumors the needless side effects of ineffective therapy), and the mechanisms of resistance in tumors that are nonresponsive to treatment (so these can be bypassed). These needs are the focus of this review, which discusses the particular issues encountered when investigating the potential of multigene expression signatures as predictive factors for response to aromatase inhibitors, which have recently become front-line endocrine therapies for postmenopausal patients with breast cancer.

Global expression profiling in epileptogenesis: does it add to the confusion?

Since the inception of global gene expression profiling platforms in the mid-1990s, there has been a significant increase in publications of differentially expressed genes in the process of epileptogenesis. In particular for mesial temporal lobe epilepsy, the presence of a latency period between the first manifestation of seizures to chronic epilepsy provides the opportunity for therapeutic interventions at the molecular biology level. Using global expression profiling techniques, approximately 2000 genes have been published demonstrating differential expression in mesial temporal epilepsy. The majority of these changes, however, are specific to laboratory or experimental conditions with only 53 genes demonstrating changes in more than two publications. To this end, we review the current status of gene expression profiling in epileptogenesis and suggest standard guidelines to be followed for greater accuracy and reproducibility of results.