Next-generation tissue microarray (ngTMA) increases the quality of biomarker studies: an example using CD3, CD8, and CD45RO in the tumor microenvironment of six different solid tumor types

Background Tissue microarray (TMA) technology revolutionized the investigation of potential biomarkers from paraffin-embedded tissues. However, conventional TMA construction is laborious, time-consuming and imprecise. Next-generation tissue microarrays (ngTMA) combine histological expertise with digital pathology and automated tissue microarraying. The aim of this study was to test the feasibility of ngTMA for the investigation of biomarkers within the tumor microenvironment (tumor center and invasion front) of six tumor types, using CD3, CD8 and CD45RO as an example. Methods Ten cases each of malignant melanoma, lung, breast, gastric, prostate and colorectal cancers were reviewed. The most representative H&E slide was scanned and uploaded onto a digital slide management platform. Slides were viewed and seven TMA annotations of 1 mm in diameter were placed directly onto the digital slide. Different colors were used to identify the exact regions in normal tissue (n = 1), tumor center (n = 2), tumor front (n = 2), and tumor microenvironment at invasion front (n = 2) for subsequent punching. Donor blocks were loaded into an automated tissue microarrayer. Images of the donor block were superimposed with annotated digital slides. Exact annotated regions were punched out of each donor block and transferred into a TMA block. 420 tissue cores created two ngTMA blocks. H&E staining and immunohistochemistry for CD3, CD8 and CD45RO were performed. Results All 60 slides were scanned automatically (total time < 10 hours), uploaded and viewed. Annotation time was 1 hour. The 60 donor blocks were loaded into the tissue microarrayer, simultaneously. Alignment of donor block images and digital slides was possible in less than 2 minutes/case. Automated punching of tissue cores and transfer took 12 seconds/core. Total ngTMA construction time was 1.4 hours. Stains for H&E and CD3, CD8 and CD45RO highlighted the precision with which ngTMA could capture regions of tumor-stroma interaction of each cancer and the T-lymphocytic immune reaction within the tumor microenvironment. Conclusion Based on a manual selection criteria, ngTMA is able to precisely capture histological zones or cell types of interest in a precise and accurate way, aiding the pathological study of the tumor microenvironment. This approach would be advantageous for visualizing proteins, DNA, mRNA and microRNAs in specific cell types using in situ hybridization techniques.

SURVIVAL PREDICTION FOR BRAIN TUMOR PATIENTS USING GENE EXPRESSION DATA

Brain tumor is one of the most aggressive types of cancer in humans, with an estimated median survival time of 12 months and only 4% of the patients surviving more than 5 years after disease diagnosis. Until recently, brain tumor prognosis has been based only on clinical information such as tumor grade and patient age, but there are reports indicating that molecular profiling of gliomas can reveal subgroups of patients with distinct survival rates. We hypothesize that coupling molecular profiling of brain tumors with clinical information might improve predictions of patient survival time and, consequently, better guide future treatment decisions. In order to evaluate this hypothesis, the general goal of this research is to build models for survival prediction of glioma patients using DNA molecular profiles (U133 Affymetrix gene expression microarrays) along with clinical information. First, a predictive Random Forest model is built for binary outcomes (i.e. short vs. long-term survival) and a small subset of genes whose expression values can be used to predict survival time is selected. Following, a new statistical methodology is developed for predicting time-to-death outcomes using Bayesian ensemble trees. Due to a large heterogeneity observed within prognostic classes obtained by the Random Forest model, prediction can be improved by relating time-to-death with gene expression profile directly. We propose a Bayesian ensemble model for survival prediction which is appropriate for high-dimensional data such as gene expression data. Our approach is based on the ensemble "sum-of-trees" model which is flexible to incorporate additive and interaction effects between genes. We specify a fully Bayesian hierarchical approach and illustrate our methodology for the CPH, Weibull, and AFT survival models. We overcome the lack of conjugacy using a latent variable formulation to model the covariate effects which decreases computation time for model fitting. Also, our proposed models provides a model-free way to select important predictive prognostic markers based on controlling false discovery rates. We compare the performance of our methods with baseline reference survival methods and apply our methodology to an unpublished data set of brain tumor survival times and gene expression data, selecting genes potentially related to the development of the disease under study. A closing discussion compares results obtained by Random Forest and Bayesian ensemble methods under the biological/clinical perspectives and highlights the statistical advantages and disadvantages of the new methodology in the context of DNA microarray data analysis.

Bicarbonate enhances expression of the endocarditis and biofilm associated pilus locus, ebpR-ebpABC, in Enterococcus faecalis.

BACKGROUND: We previously identified ebpR, encoding a potential member of the AtxA/Mga transcriptional regulator family, and showed that it is important for transcriptional activation of the Enterococcus faecalis endocarditis and biofilm associated pilus operon, ebpABC. Although ebpR is not absolutely essential for ebpABC expression (100-fold reduction), its deletion led to phenotypes similar to those of an ebpABC mutant such as absence of pili at the cell surface and, consequently, reduced biofilm formation. A non-piliated ebpABC mutant has been shown to be attenuated in a rat model of endocarditis and in a murine urinary tract infection model, indicating an important participation of the ebpR-ebpABC locus in virulence. However, there is no report relating to the environmental conditions that affect expression of the ebpR-ebpABC locus. RESULTS: In this study, we examined the effect of CO2/HCO3(-), pH, and the Fsr system on the ebpR-ebpABC locus expression. The presence of 5% CO2/0.1 M HCO3(-) increased ebpR-ebpABC expression, while the Fsr system was confirmed to be a weak repressor of this locus. The mechanism by which the Fsr system repressed the ebpR-ebpABC locus expression appears independent of the effects of CO2(-) bicarbonate. Furthermore, by using an ebpA::lacZ fusion as a reporter, we showed that addition of 0.1 M sodium bicarbonate to TSBG (buffered at pH 7.5), but not the presence of 5% CO2, induced ebpA expression in TSBG broth. In addition, using microarray analysis, we found 73 genes affected by the presence of sodium bicarbonate (abs(fold) > 2, P < 0.05), the majority of which belong to the PTS system and ABC transporter families. Finally, pilus production correlated with ebpA mRNA levels under the conditions tested. CONCLUSIONS: This study reports that the ebp locus expression is enhanced by the presence of bicarbonate with a consequential increase in the number of cells producing pili. Although the molecular basis of the bicarbonate effect remains unclear, the pathway is independent of the Fsr system. In conclusion, E. faecalis joins the growing family of pathogens that regulates virulence gene expression in response to bicarbonate and/or CO2.

Genetic networks controlling retinal injury.

PURPOSE: The present study defines genomic loci underlying coordinate changes in gene expression following retinal injury. METHODS: A group of acute phase genes expressed in diverse nervous system tissues was defined by combining microarray results from injury studies from rat retina, brain, and spinal cord. Genomic loci regulating the brain expression of acute phase genes were identified using a panel of BXD recombinant inbred (RI) mouse strains. Candidate upstream regulators within a locus were defined using single nucleotide polymorphism databases and promoter motif databases. RESULTS: The acute phase response of rat retina, brain, and spinal cord was dominated by transcription factors. Three genomic loci control transcript expression of acute phase genes in brains of BXD RI mouse strains. One locus was identified on chromosome 12 and was highly correlated with the expression of classic acute phase genes. Within the locus we identified the inhibitor of DNA binding 2 (Id2) as a candidate upstream regulator. Id2 was upregulated as an acute phase transcript in injury models of rat retina, brain, and spinal cord. CONCLUSIONS: We defined a group of transcriptional changes associated with the retinal acute injury response. Using genetic linkage analysis of natural transcript variation, we identified regulatory loci and candidate regulators that control transcript levels of acute phase genes.

LDOC1, A NOVEL BIOMARKER OF PROGNOSIS IN CHRONIC LYMPHOCYTIC LEUKEMIA

In chronic lymphocytic leukemia (CLL), one of the best predictors of outcome is the somatic mutation status of the immunoglobulin heavy chain variable region (IGHV) genes. Patients whose CLL cells have unmutated IGHV genes have a median survival of 8 years; those with mutated IGHV genes have a median survival of 25 years. To identify new prognostic biomarkers and molecular targets for therapy in untreated CLL patients, we reanalyzed the raw data from four published gene expression profiling microarray studies. Of 88 candidate biomarkers associated with IGHV somatic mutation status, we identified LDOC1 (Leucine Zipper, Down-regulated in Cancer 1), as one of the most significantly differentially expressed genes that distinguished mutated from unmutated CLL cases. LDOC1 is a putative transcription factor of unknown function in B-cell development and CLL pathophysiology. Using a highly sensitive quantitative RT-PCR (QRT-PCR) assay, we confirmed that LDOC1 mRNA was dramatically down-regulated in mutated compared to unmutated CLL cases. Expression of LDOC1 mRNA was also vii strongly associated with other markers of poor prognosis, including ZAP70 protein and cytogenetic abnormalities of poor prognosis (deletions of chromosomes 6q21, 11q23, and 17p13.1, and trisomy 12). CLL cases positive for LDOC1 mRNA had significantly shorter overall survival than negative cases. Moreover, in a multivariate model, LDOC1 mRNA expression predicted overall survival better than IGHV mutation status or ZAP70 protein, among the best markers of prognosis in CLL. We also discovered LDOC1S, a new LDOC1 splice variant. Using isoform-specific QRT-PCR assays that we developed, we found that both isoforms were expressed in normal B cells (naïve > memory), unmutated CLL cells, and in B-cell non-Hodgkin lymphomas with unmutated IGHV genes. To investigate pathways in which LDOC1 is involved, we knocked down LDOC1 in HeLa cells and performed global gene expression profiling. GFI1 (Growth Factor-Independent 1) emerged as a significantly up-regulated gene in both HeLa cells and CLL cells that expressed high levels of LDOC1. GFI1 oncoprotein is implicated in hematopoietic stem cell maintenance, lymphocyte development, and lymphomagenesis. Our findings indicate that LDOC1 mRNA is an excellent biomarker of overall survival in CLL, and may contribute to B-cell differentiation and malignant transformation.

Variability in Accumulation Rates from GPR Profiling on the West Antarctic Plateau

lsochronal layers in firn detected with ground-penetrating radar (GPR) and dated using results from ice-core analyses are used to calculate accumulation rates along a 100 km across-flow profile in West Antarctica. Accumulation rates are shown to be highly variable over short distances. Elevation measurements from global positioning system surveys show that accumulation rates derived from shallow horizons correlate well with surface undulations, which implies that wind redistribution of snow is the leading cause of this variability. Temporal changes in accumulation rate over 25-185 year intervals are smoothed to along-track length scales comparable to surface undulations in order to identify trends in accumulation that are likely related to changes in climate. Results show that accumulation rates along this profile have decreased in recent decades, which is consistent with core-derived time series of annual accumulation rates measured at the two ends of the radar profile. These results suggest that temporal variability observed in accumulation-rate records from ice cores and GPR profiles can be obscured by spatial influences, although it is possible to resolve temporal signals if the effects of local topography and ice flow are quantified and removed.

Identification of cryptic Anopheles mosquito species by molecular protein profiling.

Vector control is the mainstay of malaria control programmes. Successful vector control profoundly relies on accurate information on the target mosquito populations in order to choose the most appropriate intervention for a given mosquito species and to monitor its impact. An impediment to identify mosquito species is the existence of morphologically identical sibling species that play different roles in the transmission of pathogens and parasites. Currently PCR diagnostics are used to distinguish between sibling species. PCR based methods are, however, expensive, time-consuming and their development requires a priori DNA sequence information. Here, we evaluated an inexpensive molecular proteomics approach for Anopheles species: matrix assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS). MALDI-TOF MS is a well developed protein profiling tool for the identification of microorganisms but so far has received little attention as a diagnostic tool in entomology. We measured MS spectra from specimens of 32 laboratory colonies and 2 field populations representing 12 Anopheles species including the A. gambiae species complex. An important step in the study was the advancement and implementation of a bioinformatics approach improving the resolution over previously applied cluster analysis. Borrowing tools for linear discriminant analysis from genomics, MALDI-TOF MS accurately identified taxonomically closely related mosquito species, including the separation between the M and S molecular forms of A. gambiae sensu stricto. The approach also classifies specimens from different laboratory colonies; hence proving also very promising for its use in colony authentication as part of quality assurance in laboratory studies. While being exceptionally accurate and robust, MALDI-TOF MS has several advantages over other typing methods, including simple sample preparation and short processing time. As the method does not require DNA sequence information, data can also be reviewed at any later stage for diagnostic or functional patterns without the need for re-designing and re-processing biological material.

Observing Biogeochemical Cycles at Global Scales With Profiling Floats and Gliders Prospects for a Global Array

Chemical and biological sensor technologies have advanced rapidly in the past five years. Sensors that require low power and operate for multiple years are now available for oxygen, nitrate, and a variety of bio-optical properties that serve as proxies for important components of the carbon cycle (e.g., particulate organic carbon). These sensors have all been deployed successfully for long periods, in some cases more than three years, on platforms such as profiling floats or gliders. Technologies for pH, pCO(2), and particulate inorganic carbon are maturing rapidly as well. These sensors could serve as the enabling technology for a global biogeochemical observing system that might operate on a scale comparable to the current Argo array. Here, we review the scientific motivation and the prospects for a global observing system for ocean biogeochemistry.

Protein profiling of organic stone matrix and urine from dogs with urolithiasis

Two-thirds of the organic matrix in urinary stones consists of proteins. Their relationship to calculogenesis remains controversial with regard to their effect as inhibitors or promoters during stone formation. The purpose of the present study was to determine the differences in peptide and protein pattern between the urine of stone formers (n = 23) and control dogs (n = 12), as well as between organic matrix of different urinary stones (struvite n = 11, calcium oxalate n = 8, uric acid n = 4) using surface-enhanced laser desorption/ionization time-of-flight mass spectrometry. Specific differences in protein and peptide profiles were found in the organic matrix of different mineral compositions. Characteristic differences were also found in urinary peptide and protein pattern especially in molecular masses below 20 kDa between affected and healthy dogs. Based on the obtained molecular masses they were in some cases tentatively identified as proteins that are known to be involved in stone formation in humans. The study shows that in dogs, specific-urinary peptides and proteins might be associated with urolithiasis. It indicates the importance to further characterize those proteins for possible diagnostic purposes in prognosis and therapy

Protein profiling of urine from dogs with renal disease using ProteinChip analysis

Measurement of total urinary proteins in individuals that tested positive by urinary dipstick is a typical method for assessing the presence of potentially serious renal disorders. In the absence of such overt proteinuria, however, measurement of specific urinary proteins may be useful in the diagnosis of nephropathies and may provide greater insight into the pathogenesis. The urine of 28 dogs (16 with renal disease and 12 healthy) was evaluated to determine whether specific low-molecular-weight proteins or the pattern of protein excretion could also be used as a marker of tubular dysfunction in dogs. Specific proteins were assessed by immunological methods, whereas protein profiles were determined by surface-enhanced laser desorption/ionization time-of-flight mass spectrometry (MS). In particular, changes in the excretion of retinol-binding protein (RBP) and Tamm-Horsfall protein (THP) appear to be of clinical relevance in the diagnosis of canine kidney diseases. The pattern of urinary protein and peptides revealed specific changes in abundance in dogs with renal disease at molecular masses (kD) of 11.58, 12.41, 12.60, 14.58, 20.95 (RBP), 27.85, and 65.69 (albumin). In conclusion, comparable proteins as in humans might be used as urinary markers for proximal (RBP) and distal (THP) tubular dysfunction in dogs. Surface-enhanced laser desorption/ionization time-of-flight MS is a promising tool for the study of kidney physiology and pathophysiology and might aid in the discovery of new biomarkers of renal disease.

A next-generation tissue microarray (ngTMA) protocol for biomarker studies

Biomarker research relies on tissue microarrays (TMA). TMAs are produced by repeated transfer of small tissue cores from a 'donor' block into a 'recipient' block and then used for a variety of biomarker applications. The construction of conventional TMAs is labor intensive, imprecise, and time-consuming. Here, a protocol using next-generation Tissue Microarrays (ngTMA) is outlined. ngTMA is based on TMA planning and design, digital pathology, and automated tissue microarraying. The protocol is illustrated using an example of 134 metastatic colorectal cancer patients. Histological, statistical and logistical aspects are considered, such as the tissue type, specific histological regions, and cell types for inclusion in the TMA, the number of tissue spots, sample size, statistical analysis, and number of TMA copies. Histological slides for each patient are scanned and uploaded onto a web-based digital platform. There, they are viewed and annotated (marked) using a 0.6-2.0 mm diameter tool, multiple times using various colors to distinguish tissue areas. Donor blocks and 12 'recipient' blocks are loaded into the instrument. Digital slides are retrieved and matched to donor block images. Repeated arraying of annotated regions is automatically performed resulting in an ngTMA. In this example, six ngTMAs are planned containing six different tissue types/histological zones. Two copies of the ngTMAs are desired. Three to four slides for each patient are scanned; 3 scan runs are necessary and performed overnight. All slides are annotated; different colors are used to represent the different tissues/zones, namely tumor center, invasion front, tumor/stroma, lymph node metastases, liver metastases, and normal tissue. 17 annotations/case are made; time for annotation is 2-3 min/case. 12 ngTMAs are produced containing 4,556 spots. Arraying time is 15-20 hr. Due to its precision, flexibility and speed, ngTMA is a powerful tool to further improve the quality of TMAs used in clinical and translational research.

A novel universal DNA labeling and amplification system for rapid microarray-based detection of 117 antibiotic resistance genes in Gram-positive bacteria.

A rapid and simple DNA labeling system has been developed for disposable microarrays and has been validated for the detection of 117 antibiotic resistance genes abundant in Gram-positive bacteria. The DNA was fragmented and amplified using phi-29 polymerase and random primers with linkers. Labeling and further amplification were then performed by classic PCR amplification using biotinylated primers specific for the linkers. The microarray developed by Perreten et al. (Perreten, V., Vorlet-Fawer, L., Slickers, P., Ehricht, R., Kuhnert, P., Frey, J., 2005. Microarray-based detection of 90 antibiotic resistance genes of gram-positive bacteria. J.Clin.Microbiol. 43, 2291-2302.) was improved by additional oligonucleotides. A total of 244 oligonucleotides (26 to 37 nucleotide length and with similar melting temperatures) were spotted on the microarray, including genes conferring resistance to clinically important antibiotic classes like β-lactams, macrolides, aminoglycosides, glycopeptides and tetracyclines. Each antibiotic resistance gene is represented by at least 2 oligonucleotides designed from consensus sequences of gene families. The specificity of the oligonucleotides and the quality of the amplification and labeling were verified by analysis of a collection of 65 strains belonging to 24 species. Association between genotype and phenotype was verified for 6 antibiotics using 77 Staphylococcus strains belonging to different species and revealed 95% test specificity and a 93% predictive value of a positive test. The DNA labeling and amplification is independent of the species and of the target genes and could be used for different types of microarrays. This system has also the advantage to detect several genes within one bacterium at once, like in Staphylococcus aureus strain BM3318, in which up to 15 genes were detected. This new microarray-based detection system offers a large potential for applications in clinical diagnostic, basic research, food safety and surveillance programs for antimicrobial resistance.

Circulating microRNA profiling in patients with advanced non-squamous NSCLC receiving bevacizumab/erlotinib followed by platinum-based chemotherapy at progression (SAKK 19/05).

OBJECTIVES Molecular subclassification of non small-cell lung cancer (NSCLC) is essential to improve clinical outcome. This study assessed the prognostic and predictive value of circulating micro-RNA (miRNA) in patients with non-squamous NSCLC enrolled in the phase II SAKK (Swiss Group for Clinical Cancer Research) trial 19/05, receiving uniform treatment with first-line bevacizumab and erlotinib followed by platinum-based chemotherapy at progression. MATERIALS AND METHODS Fifty patients with baseline and 24 h blood samples were included from SAKK 19/05. The primary study endpoint was to identify prognostic (overall survival, OS) miRNA's. Patient samples were analyzed with Agilent human miRNA 8x60K microarrays, each glass slide formatted with eight high-definition 60K arrays. Each array contained 40 probes targeting each of the 1347 miRNA. Data preprocessing included quantile normalization using robust multi-array average (RMA) algorithm. Prognostic and predictive miRNA expression profiles were identified by Spearman's rank correlation test (percentage tumor shrinkage) or log-rank testing (for time-to-event endpoints). RESULTS Data preprocessing kept 49 patients and 424 miRNA for further analysis. Ten miRNA's were significantly associated with OS, with hsa-miR-29a being the strongest prognostic marker (HR=6.44, 95%-CI 2.39-17.33). Patients with high has-miR-29a expression had a significantly lower survival at 10 months compared to patients with a low expression (54% versus 83%). Six out of the 10 miRNA's (hsa-miRN-29a, hsa-miR-542-5p, hsa-miR-502-3p, hsa-miR-376a, hsa-miR-500a, hsa-miR-424) were insensitive to perturbations according to jackknife cross-validation on their HR for OS. The respective principal component analysis (PCA) defined a meta-miRNA signature including the same 6 miRNA's, resulting in a HR of 0.66 (95%-CI 0.53-0.82). CONCLUSION Cell-free circulating miRNA-profiling successfully identified a highly prognostic 6-gene signature in patients with advanced non-squamous NSCLC. Circulating miRNA profiling should further be validated in external cohorts for the selection and monitoring of systemic treatment in patients with advanced NSCLC.