Evaluation of reference-based two-color methods for measurement of gene expression ratios using spotted cDNA microarrays

Circulating Tumor Cell Transcriptomics as Biopsy Surrogates in Metastatic Breast Cancer

BACKGROUND

Metastatic breast cancer (MBC) and the circulating tumor cells (CTCs) leading to macrometastases are inherently different than primary breast cancer. We evaluated whether whole transcriptome RNA-Seq of CTCs isolated via an epitope-independent approach may serve as a surrogate for biopsies of macrometastases.

METHODS

We performed RNA-Seq on fresh metastatic tumor biopsies, CTCs, and peripheral blood (PB) from 19 newly diagnosed MBC patients. CTCs were harvested using the ANGLE Parsortix microfluidics system to isolate cells based on size and deformability, independent of a priori knowledge of cell surface marker expression.

RESULTS

Gene expression separated CTCs, metastatic biopsies, and PB into distinct groups despite heterogeneity between patients and sample types. CTCs showed higher expression of immune oncology targets compared with corresponding metastases and PB. Predictive biomarker (n = 64) expression was highly concordant for CTCs and metastases. Repeat observation data post-treatment demonstrated changes in the activation of different biological pathways. Somatic single nucleotide variant analysis showed increasing mutational complexity over time.

CONCLUSION

We demonstrate that RNA-Seq of CTCs could serve as a surrogate biomarker for breast cancer macrometastasis and yield clinically relevant insights into disease biology and clinically actionable targets.

FGF23 modulates the effects of erythropoietin on gene expression in renal epithelial cells

Background

FGF23 plays an important role in calcium–phosphorus metabolism. Other roles of FGF23 have recently been reported, such as commitment to myocardium enlargement and immunological roles in the spleen. In this study, we aimed to identify the roles of FGF23 in the kidneys other than calcium–phosphorus metabolism.

Methods

DNA microarrays and bioinformatics tools were used to analyze gene expression in mIMCD3 mouse renal tubule cells following treatment with FGF23, erythropoietin and/or an inhibitor of ERK.

Results

Three protein-coding genes were upregulated and 12 were downregulated in response to FGF23. Following bioinformatics analysis of these genes, PPARγ and STAT3 were identified as candidate transcript factors for mediating their upregulation, and STAT1 as a candidate for mediating their downregulation. Because STAT1 and STAT3 also mediate erythropoietin signaling, we investigated whether FGF23 and erythropoietin might show interactive effects in these cells. Of the 15 genes regulated by FGF23, 11 were upregulated by erythropoietin; 10 of these were downregulated following cotreatment with FGF23. Inhibition of ERK, an intracellular mediator of FGF23, reversed the effects of FGF23. However, FGF23 did not influence STAT1 phosphorylation, suggesting that it impinges on erythropoietin signaling through other mechanisms.

Conclusion

Our results suggest cross talk between erythropoietin and FGF23 signaling in the regulation of renal epithelial cells.

Comparison of performance of one-color and two-color gene-expression analyses in predicting clinical endpoints of neuroblastoma patients

Microarray-based prediction of clinical endpoints may be performed using either a one-color approach reflecting mRNA abundance in absolute intensity values or a two-color approach yielding ratios of fluorescent intensities. In this study, as part of the MAQC-II project, we systematically compared the classification performance resulting from one- and two-color gene-expression profiles of 478 neuroblastoma samples. In total, 196 classification models were applied to these measurements to predict four clinical endpoints, and classification performances were compared in terms of accuracy, area under the curve, Matthews correlation coefficient and root mean-squared error. Whereas prediction performance varied with distinct clinical endpoints and classification models, equivalent performance metrics were observed for one- and two-color measurements in both internal and external validation. Furthermore, overlap of selected signature genes correlated inversely with endpoint prediction difficulty. In summary, our data strongly substantiate that the choice of platform is not a primary factor for successful gene expression based-prediction of clinical endpoints.

Development of a common oligonucleotide reference standard for microarray data normalization and comparison across different microbial communities

High-density functional gene arrays have become a powerful tool for environmental microbial detection and characterization. However, microarray data normalization and comparison for this type of microarray remain a challenge in environmental microbiology studies because some commonly used normalization methods (e.g., genomic DNA) for the study of pure cultures are not applicable. In this study, we developed a common oligonucleotide reference standard (CORS) method to address this problem. A unique 50-mer reference oligonucleotide probe was selected to co-spot with gene probes for each array feature. The complementary sequence was synthesized and labeled for use as the reference target, which was then spiked and cohybridized with each sample. The signal intensity of this reference target was used for microarray data normalization and comparison. The optimal amount or concentration were determined to be ca. 0.5 to 2.5% of a gene probe for the reference probe and ca. 0.25 to 1.25 fmol/microl for the reference target based on our evaluation with a pilot array. The CORS method was then compared to dye swap and genomic DNA normalization methods using the Desulfovibrio vulgaris whole-genome microarray, and significant linear correlations were observed. This method was then applied to a functional gene array to analyze soil microbial communities, and the results demonstrated that the variation of signal intensities among replicates based on the CORS method was significantly lower than the total intensity normalization method. The developed CORS provides a useful approach for microarray data normalization and comparison for studies of complex microbial communities.

Oligonucleotide microarrays for bacteriophage expression studies

Gene expression microarrays offer the ability to monitor the expression of all phage genes over an infection cycle. However, there are relatively few reports to date of microarrays being used to investigate phage biology. This chapter aims to provide an overview of how to design and implement a microarray experiment to investigate phage biology. Given the nature of microarrays being specific to an organism, each will provide a number of unique issues. In this chapter, we outline the basic theory behind microarrays and provide details on how to implement a microarray experiment from the design of oligonucleotide probes through to the hybridisation of microarrays. The matter of designing oligonucleotide probes will be discussed with regards to how probe length, secondary structure, free energy, probe orientation and amplification all have to be taken into account. As means of an example, the conditions used for the hybridisation of an array designed to be specific to the cyanophage S-PM2 is detailed.

Optimization of experimental design parameters for high-throughput chromatin immunoprecipitation studies

High-throughput, microarray-based chromatin immunoprecipitation (ChIP-chip) technology allows in vivo elucidation of transcriptional networks. However this complex is not yet readily accessible, in part because its many parameters have not been systematically evaluated and optimized. We address this gap by systematically assessing experimental-design parameters including antibody purity, dye-bias, array-batch, inter-day hybridization bias, amplification method and choice of hybridization control. The combined performance of these optimized parameters shows a 90% validation rate in ChIP-chip analysis of Myc genomic binding in HL60 cells using two different microarray platforms. Increased sensitivity and decreased noise in ChIP-chip assays will enable wider use of this methodology to accurately and affordably elucidate transcriptional networks.

Splice variants of TLE family genes and up-regulation of a TLE3 isoform in prostate tumors

The TLE genes constitute a family of important transcriptional co-repressors involved in many cellular processes. We found evidence of alternatively spliced mRNAs for human TLE1-4 containing premature stop codons, thus encoding putative shortened proteins. Microarray experiments and Real-time RT-PCR assays showed that alternatively spliced isoforms of TLE1, TLE2 an d TLE3 were preferentially expressed in prostate in comparison to liver and kidney tissues. We identified by orientation-specific R T-PCR an antisense partially intronic non-coding RNA that overlaps a novel exon of the TLE3 gene, raising the possibility of regulation of alternative splicing by this non-coding transcript. The alternatively spliced isoform of TLE3 was up-regulated (6- to 17-fo ld) in prostate tumors in comparison to matched non-tumor adjacent tissue from 7 out of 11 (64%) patients and in four prostate tumor cell lines in comparison to a normal prostate cell line. These results demonstrate that different isoforms of TLE genes are commonly transcribed in human tissues and suggest that TLE3 could be involved in prostate cancer development.

Third-generation human mitochondria-focused cDNA microarray and its bioinformatic tools for analysis of gene expression

To facilitate profiling mitochondrial transcriptomes, we developed a third-generation human mitochondria-focused cDNA microarray (hMitChip3) and its bioinformatic tools. hMitChip3 consists of the 37 mitochondrial DNA-encoded genes, 1098 nuclear DNA-encoded and mitochondria-related genes, and 225 controls, each in triplicate. The bioinformatic tools included data analysis procedures and customized database for interpretation of results. The database associated 645 molecular functions with 946 hMitChip3 genes, 612 biological processes with 930 genes, 172 cellular components with 869 genes, 107 biological chemistry pathways with 476 genes, 23 reactome events with 227 genes, 320 genetic disorders with 237 genes, and 87 drugs targets with 55 genes. To test these tools, hMitChip3 was used to compare expression profiles between human melanoma cell lines UACC903 (rapidly dividing) and UACC903(+6) (slowly dividing). Our results demonstrated internal gene-set consistency (correlation R > or = 0.980 +/- 0.005) and interexperimental reproducibility (R > or = 0.931 +/- 0.013). Expression patterns of 16 genes, involved in DNA, RNA, or protein biosyntheses in mitochondrial and other organelles, were consistent with the proliferation rates of both cell lines, and the patterns of 6 tested genes were verified by quantitative reverse transcription PCR (RT-PCR). Thus, hMitChip3 and its bioinformatics software provide an integrated tool for profiling mitochondria-focused gene expression.

Androgen responsive intronic non-coding RNAs

BACKGROUND

Transcription of large numbers of non-coding RNAs originating from intronic regions of human genes has been recently reported, but mechanisms governing their biosynthesis and biological functions are largely unknown. In this work, we evaluated the existence of a common mechanism of transcription regulation shared by protein-coding mRNAs and intronic RNAs by measuring the effect of androgen on the transcriptional profile of a prostate cancer cell line.

RESULTS

Using a custom-built cDNA microarray enriched in intronic transcribed sequences, we found 39 intronic non-coding RNAs for which levels were significantly regulated by androgen exposure. Orientation-specific reverse transcription-PCR indicated that 10 of the 13 were transcribed in the antisense direction. These transcripts are long (0.5-5 kb), unspliced and apparently do not code for proteins. Interestingly, we found that the relative levels of androgen-regulated intronic transcripts could be correlated with the levels of the corresponding protein-coding gene (asGAS6 and asDNAJC3) or with the alternative usage of exons (asKDELR2 and asITGA6) in the corresponding protein-coding transcripts. Binding of the androgen receptor to a putative regulatory region upstream from asMYO5A, an androgen-regulated antisense intronic transcript, was confirmed by chromatin immunoprecipitation.

CONCLUSION

Altogether, these results indicate that at least a fraction of naturally transcribed intronic non-coding RNAs may be regulated by common physiological signals such as hormones, and further corroborate the notion that the intronic complement of the transcriptome play functional roles in the human gene-expression program.

Genome mapping and expression analyses of human intronic noncoding RNAs reveal tissue-specific patterns and enrichment in genes related to regulation of transcription

BACKGROUND

RNAs transcribed from intronic regions of genes are involved in a number of processes related to post-transcriptional control of gene expression. However, the complement of human genes in which introns are transcribed, and the number of intronic transcriptional units and their tissue expression patterns are not known.

RESULTS

A survey of mRNA and EST public databases revealed more than 55,000 totally intronic noncoding (TIN) RNAs transcribed from the introns of 74% of all unique RefSeq genes. Guided by this information, we designed an oligoarray platform containing sense and antisense probes for each of 7,135 randomly selected TIN transcripts plus the corresponding protein-coding genes. We identified exonic and intronic tissue-specific expression signatures for human liver, prostate and kidney. The most highly expressed antisense TIN RNAs were transcribed from introns of protein-coding genes significantly enriched (p = 0.002 to 0.022) in the 'Regulation of transcription' Gene Ontology category. RNA polymerase II inhibition resulted in increased expression of a fraction of intronic RNAs in cell cultures, suggesting that other RNA polymerases may be involved in their biosynthesis. Members of a subset of intronic and protein-coding signatures transcribed from the same genomic loci have correlated expression patterns, suggesting that intronic RNAs regulate the abundance or the pattern of exon usage in protein-coding messages.

CONCLUSION

We have identified diverse intronic RNA expression patterns, pointing to distinct regulatory roles. This gene-oriented approach, using a combined intron-exon oligoarray, should permit further comparative analysis of intronic transcription under various physiological and pathological conditions, thus advancing current knowledge about the biological functions of these noncoding RNAs.

Gender biased differential alternative splicing patterns of the transcriptional cofactor CA150 gene in Schistosoma mansoni

The complex molecular systems involved in the process of sex-differentiation and fertility in Schistosoma mansoni have not yet been completely described. Using a 4608-element cDNA microarray, we have now determined 90 and 139 genes with significantly (q-value</=0.06) higher expression levels in adult males and females, respectively. Eight out of eleven (73%) selected transcripts had their differential expression levels validated by real-time RT-PCR. One of these transcripts was extended by RT-PCR and was shown to span the intronic region between exons 9 and 11 of the S. mansoni CA150 gene, a transcriptional cofactor known in humans to interact with both RNA polymerase II and the spliceosome complex. The longer transcript probably represents a novel isoform of S. mansoni CA150. Additionally, we obtained full-length sequences for three other isoforms of the SmCA150 gene, coding for proteins of different lengths and domain compositions. Semi-quantitative RT-PCR showed different expression ratios among these isoforms between male and female. Due to the role of CA150 in RNA transcription and processing, we hypothesize that these differential expression events may be important in the generation and maintenance of the different phenotypes between male and female.