Design, construction, and analysis of cell line arrays and tissue microarrays for gene expression analysis

Cell line array (CMA) and tissue microarray (TMA) technologies are high-throughput methods for analysing both the abundance and distribution of gene expression in a panel of cell lines or multiple tissue specimens in an efficient and cost-effective manner. The process is based on Kononen's method of extracting a cylindrical core of paraffin-embedded donor tissue and inserting it into a recipient paraffin block. Donor tissue from surgically resected paraffin-embedded tissue blocks, frozen needle biopsies or cell line pellets can all be arrayed in the recipient block. The representative area of interest is identified and circled on a haematoxylin and eosin (H&E)-stained section of the donor block. Using a predesigned map showing a precise spacing pattern, a high density array of up to 1,000 cores of cell pellets and/or donor tissue can be embedded into the recipient block using a tissue arrayer from Beecher Instruments. Depending on the depth of the cell line/tissue removed from the donor block 100-300 consecutive sections can be cut from each CMA/TMA block. Sections can be stained for in situ detection of protein, DNA or RNA targets using immunohistochemistry (IHC), fluorescent in situ hybridisation (FISH) or mRNA in situ hybridisation (RNA-ISH), respectively. This chapter provides detailed methods for CMA/TMA design, construction and analysis with in-depth notes on all technical aspects including tips to deal with common pitfalls the user may encounter. © Springer Science+Business Media, LLC 2011.

Gene expression analysis of diagnostic biopsies predicts pathological response to neoadjuvant chemoradiotherapy of esophageal cancer

OBJECTIVE: This study explored gene expression differences in predicting response to chemoradiotherapy in esophageal cancer. PURPOSE:: A major pathological response to neoadjuvant chemoradiation is observed in about 40% of esophageal cancer patients and is associated with favorable outcomes. However, patients with tumors of similar histology, differentiation, and stage can have vastly different responses to the same neoadjuvant therapy. This dichotomy may be due to differences in the molecular genetic environment of the tumor cells. BACKGROUND DATA: Diagnostic biopsies were obtained from a training cohort of esophageal cancer patients (13), and extracted RNA was hybridized to genome expression microarrays. The resulting gene expression data was verified by qRT-PCR. In a larger, independent validation cohort (27), we examined differential gene expression by qRT-PCR. The ability of differentially-regulated genes to predict response to therapy was assessed in a multivariate leave-one-out cross-validation model. RESULTS: Although 411 genes were differentially expressed between normal and tumor tissue, only 103 genes were altered between responder and non-responder tumor; and 67 genes differentially expressed >2-fold. These included genes previously reported in esophageal cancer and a number of novel genes. In the validation cohort, 8 of 12 selected genes were significantly different between the response groups. In the predictive model, 5 of 8 genes could predict response to therapy with 95% accuracy in a subset (74%) of patients. CONCLUSIONS: This study has identified a gene microarray pattern and a set of genes associated with response to neoadjuvant chemoradiation in esophageal cancer. The potential of these genes as biomarkers of response to treatment warrants further investigation. Copyright © 2009 by Lippincott Williams & Wilkins.

Genome analysis reveals insights into physiology and longevity of the Brandt’s bat Myotis brandtii

Bats account for one-fifth of mammalian species, are the only mammals with powered flight, and are among the few animals that echolocate. The insect-eating Brandt’s bat (Myotis brandtii) is the longest-lived bat species known to date (lifespan exceeds 40 years) and, at 4–8 g adult body weight, is the most extreme mammal with regard to disparity between body mass and longevity. Here we report sequencing and analysis of the Brandt’s bat genome and transcriptome, which suggest adaptations consistent with echolocation and hibernation, as well as altered metabolism, reproduction and visual function. Unique sequence changes in growth hormone and insulin-like growth factor 1 receptors are also observed. The data suggest that an altered growth hormone/insulin-like growth factor 1 axis, which may be common to other long-lived bat species, together with adaptations such as hibernation and low reproductive rate, contribute to the exceptional lifespan of the Brandt’s bat.

A multi-genome analysis approach enables tracking of the invasion of a single Russian wheat aphid (Diuraphis noxia) clone throughout the New World

This study investigated the population genetics, demographic history and pathway of invasion of the Russian wheat aphid (RWA) from its native range in Central Asia, the Middle East and Europe to South Africa and the Americas. We screened microsatellite markers, mitochondrial DNA and endosymbiont genes in 504 RWA clones from nineteen populations worldwide. Following pathway analyses of microsatellite and endosymbiont data, we postulate that Turkey and Syria were the most likely sources of invasion to Kenya and South Africa, respectively. Furthermore, we found that one clone transferred between South Africa and the Americas was most likely responsible for the New World invasion. Finally, endosymbiont DNA was found to be a high resolution population genetic marker, extremely useful for studies of invasion over a relatively short evolutionary history time frame. This study has provided valuable insights into the factors that may have facilitated the recent global invasion by this damaging pest.

Determining epithelial contribution to in vivo mesenchymal tumour expression signature using species-specific microarray profiling analysis of xenografts

Gene expression profiling using microarrays and xenograft transplants of human cancer cell lines are both popular tools to investigate human cancer. However, the undefined degree of cross hybridization between the mouse and human genomes hinders the use of microarrays to characterize gene expression of both the host and the cancer cell within the xenograft. Since an increasingly recognized aspect of cancer is the host response (or cancer-stroma interaction), we describe here a bioinformatic manipulation of the Affymetrix profiling that allows interrogation of the gene expression of both the mouse host and the human tumour. Evidence of microenvironmental regulation of epithelial mesenchymal transition of the tumour component in vivo is resolved against a background of mesenchymal gene expression. This tool could allow deeper insight to the mechanism of action of anti-cancer drugs, as typically novel drug efficacy is being tested in xenograft systems.

Molecular and cellular analysis of basement membrane invasion by human breast cancer cells in Matrigel-based in vitro assays

In vitro analyses of basement membrane invasiveness employing Matrigel (a murine tumor extract rich in basement membrane components) have been performed on human breast cancer model systems. Constitutive invasiveness of different human breast cancer (HBC) cell lines has been examined as well as regulation by steroid hormones, growth factors, and oncogenes. Carcinoma cells exhibiting a mesenchymal-like phenotype (vimentin expression, lack of cell border associated uvomorulin) show dramatically increased motility, invasiveness, and metastatic potential in nude mice. These findings support the hypothesis that epithelial to mesenchymal transition (EMT)-like events may be instrumental in the metastatic progression of human breast cancer. The MCF-7 subline MCF-7ADR appears to have undergone such a transition. The importance of such a transition may be reflected in the emergence of vimentin expression as an indicator of poor prognosis in HBC. Matrix degradation and laminin recognition are highlighted as potential targets for antimetastatic therapy, and analyses of laminin attachment and the matrix metalloproteinase (MMP) family in HBC cell lines are summarized. Matrigel-based assays have proved useful in the study of the molecular mechanisms of basement membrane invasiveness, their regulation in HBC cells, and their potential as targets for antimetastatic therapy.

Integrated analysis of epigenomic and genomic changes by DNA methylation dependent mechanisms provides potential novel biomarkers for prostate cancer

Epigenetic silencing mediated by CpG methylation is a common feature of many cancers. Characterizing aberrant DNA methylation changes associated with tumor progression may identify potential prognostic markers for prostate cancer (PCa). We treated two PCa cell lines, 22Rv1 and DU-145 with the demethylating agent 5-Aza 2’–deoxycitidine (DAC) and global methylation status was analyzed by performing methylation-sensitive restriction enzyme based differential methylation hybridization strategy followed by genome-wide CpG methylation array profiling. In addition, we examined gene expression changes using a custom microarray. Gene Set Enrichment Analysis (GSEA) identified the most significantly dysregulated pathways. In addition, we assessed methylation status of candidate genes that showed reduced CpG methylation and increased gene expression after DAC treatment, in Gleason score (GS) 8 vs. GS6 patients using three independent cohorts of patients; the publically available The Cancer Genome Atlas (TCGA) dataset, and two separate patient cohorts. Our analysis, by integrating methylation and gene expression in PCa cell lines, combined with patient tumor data, identified novel potential biomarkers for PCa patients. These markers may help elucidate the pathogenesis of PCa and represent potential prognostic markers for PCa patients.

A comparative analysis of Mitochondrial genomes in Coleoptera (Arthropoda: Insecta) and genome descriptions of six new beetles

Coleoptera is the most diverse group of insects with over 360,000 described species divided into four suborders: Adephaga, Archostemata, Myxophaga, and Polyphaga. In this study, we present six new complete mitochondrial genome (mtgenome) descriptions, including a representative of each suborder, and analyze the evolution of mtgenomes from a comparative framework using all available coleopteran mtgenomes. We propose a modification of atypical cox1 start codons based on sequence alignment to better reflect the conservation observed across species as well as findings of TTG start codons in other genes. We also analyze tRNA-Ser(AGN) anticodons, usually GCU in arthropods, and report a conserved UCU anticodon as a possible synapomorphy across Polyphaga. We further analyze the secondary structure of tRNA-Ser(AGN) and present a consensus structure and an updated covariance model that allows tRNAscan-SE (via the COVE software package) to locate and fold these atypical tRNAs with much greater consistency. We also report secondary structure predictions for both rRNA genes based on conserved stems. All six species of beetle have the same gene order as the ancestral insect. We report noncoding DNA regions, including a small gap region of about 20 bp between tRNA-Ser(UCN) and nad1 that is present in all six genomes, and present results of a base composition analysis.

Rapid analysis of GSTM1, GSTT1 and GSTP1 polymorphisms using real-time polymerase chain reaction

Polymorphisms of glutathione transferases (GST) are important genetic determinants of susceptibility to environmental carcinogens (Rebbeck, 1997). The GSTs are a multigene family of dimeric enzymes involved in detoxification, and, in a few cases, the bioactivation of a variety of xenobiotics (Hayes et al., 1995). The cytosolic GST enzyme family consists of four major classes of enzymes, referred to as alpha, mu, pi and theta. Several members of this family (for example, GSTM1, GSTT1 and GSTP1) are polymorphic in human populations (Wormhoudt et al., 1999). Molecular epidemiology studies have examined the role of GST polymorphisms as susceptibility factors for environmentally and/or occupationally induced cancers (Wormhoudt et al., 1999). In particular, case-control studies showed a relationship between the GSTM1 null genotype and the development of cancer in association with smoking habits, which has been shown for cancers of the respiratory and gastrointestinal tracts as well as other cancer types (Miller et al., 1997). Only a few molecular epidemiological studies addressed the role of GSTT1 and GSTP1 polymorphisms in cancer susceptibility. Since GSTP1 is a key player in biotransformation/bioactivation of benzo(a)pyrene, GSTP1 may be even more important than GSTM1 in the prevention of tobacco-induced cancers (Harries et al., 1997; Harris et al., 1998). To date, this relationship has not been sufficiently addressed in humans. Comprehensive molecular epidemiological studies may add to the current knowledge of the role of GST polymorphisms in cancer susceptibility and extent of the knowledge gained from approaches that used phenotyping, such as GSTM1 activity as it relates to trans-stilbene oxide, or polymerase chain reaction (PCR) based genotyping of polymorphic isoenzymes (Bell et al., 1993; Pemble et al., 1994; Harries et al., 1997).

Conservation genetics of the water mouse, Xeromys myoides Thomas, 1889

The water mouse, Xeromys myoides, is currently recognised as a vulnerable species in Australia, inhabiting a small number of distinct and isolated coastal regions of Queensland and the Northern Territory. An examination of the evolutionary history and contemporary influences shaping the genetic structure of this species is required to make informed conservation management decisions. Here, we report the first analysis undertaken on the phylogeography and population genetics of the water mouse across its mainland Australian distribution. Genetic diversity was assessed at two mitochondrial DNA (Cytochrome b, 1000 bp; D-loop, 400 bp) and eight microsatellite DNA loci. Very low genetic diversity was found, indicating that water mice underwent a recent expansion throughout their Australian range and constitute a single evolutionarily significant unit. Microsatellite analyses revealed that the highest genetic diversity was found in the Mackay region of central Queensland; population substructure was also identified, suggesting that local populations may be isolated in this region. Conversely, genetic diversity in the Coomera region of south-east Queensland was very low and the population in this region has experienced a significant genetic bottleneck. These results have significant implications for future management, particularly in terms of augmenting populations through translocations or reintroducing water mice in areas where they have gone extinct.

Associations between ghrelin and ghrelin receptor polymorphisms and cancer in Caucasian populations: A meta-analysis

Background There is growing evidence that the ghrelin axis, including ghrelin (GHRL) and its receptor, the growth hormone secretagogue receptor (GHSR), play a role in cancer progression. Ghrelin gene and ghrelin receptor gene polymorphisms have been reported to have a range of effects in cancer, from increased risk, to protection from cancer, or having no association. In this study we aimed to clarify the role of ghrelin and ghrelin receptor polymorphisms in cancer by performing a meta-analysis of published case–control studies. We conducted searches of the literature published up to January 2013 in MEDLINE using the PubMed search engine. Individual data on 8,430 cases and 14,008 controls from six case–control studies of an all Caucasian population were evaluated for three ghrelin gene (GHRL; rs696217, rs4684677, rs2075356) and one ghrelin receptor (GHSR; rs572169) polymorphism in breast cancer, esophageal cancer, colorectal cancer and non-Hodgkins lymphoma. Results In the overall analysis, homozygous and recessive associations indicated that the minor alleles of rs696217 and rs2075356 GHRL polymorphisms conferred reduced cancer risk (odds ratio [OR] 0.61-0.78). The risk was unchanged for breast cancer patients when analysed separately (OR 0.73-0.83). In contrast, the rs4684677 GHRL and the rs572169 GHSR polymorphisms conferred increased breast cancer risk (OR 1.97-1.98, p = 0.08 and OR 1.42-1.43, p = 0.08, respectively). All dominant and co-dominant effects showed null effects (OR 0.96-1.05), except for the rs572169 co-dominant effect, with borderline increased risk (OR 1.08, p = 0.05). Conclusions This study suggests that the rs696217 and rs2075356 ghrelin gene (GHRL) polymorphisms may protect carriers against breast cancer, and the rs4684677 GHRL and rs572169 GHSR polymorphisms may increase the risk among carriers. In addition, larger studies are required to confirm these findings.

Genetic characterisation and structural analysis of the O-specific polysaccharide of Yersinia pseudotuberculosis serotype O:1c

Many, but not all, of the current 21 serotypes of Yersinia pseudotuberculosis have been investigated with regard to the chemical structures of their O-specific polysaccharide (OPS) and the genetic basis of their biosynthesis. Completion of the genetics and structures of the remaining serotypes will enhance our understanding of the emerging relationship between genetics and structures within this species. Here, we present a structural and genetic analysis of the Y. pseudotuberculosis serotype O:1c OPS. Our results showed that this OPS has the same backbone as Y. pseudotuberculosis O:2b, but with a 3,6-dideoxy-D-ribo-hexofuranose (paratofuranose, Parf) side-branch instead of a 3,6-dideoxy-D-xylo-hexopyranose (abequopyranose, Abep). The 3'-end of the gene cluster is the same as for O:2b and has the genes for synthesis of the backbone and for processing the completed repeat unit. The 5'-end of the cluster consists of the same genes as O:1b for synthesis of Parf and a related gene for its transfer to the repeating unit backbone.

Genetic risk score analysis indicates migraine with and without comorbid depression are genetically different disorders

Migraine and major depressive disorder (MDD) are comorbid, moderately heritable and to some extent influenced by the same genes. In a previous paper, we suggested the possibility of causality (one trait causing the other) underlying this comorbidity. We present a new application of polygenic (genetic risk) score analysis to investigate the mechanisms underlying the genetic overlap of migraine and MDD. Genetic risk scores were constructed based on data from two discovery samples in which genome-wide association analyses (GWA) were performed for migraine and MDD, respectively. The Australian Twin Migraine GWA study (N = 6,350) included 2,825 migraine cases and 3,525 controls, 805 of whom met the diagnostic criteria for MDD. The RADIANT GWA study (N = 3,230) included 1,636 MDD cases and 1,594 controls. Genetic risk scores for migraine and for MDD were used to predict pure and comorbid forms of migraine and MDD in an independent Dutch target sample (NTR-NESDA, N = 2,966), which included 1,476 MDD cases and 1,058 migraine cases (723 of these individuals had both disorders concurrently). The observed patterns of prediction suggest that the 'pure' forms of migraine and MDD are genetically distinct disorders. The subgroup of individuals with comorbid MDD and migraine were genetically most similar to MDD patients. These results indicate that in at least a subset of migraine patients with MDD, migraine may be a symptom or consequence of MDD. © 2013 Springer-Verlag Berlin Heidelberg.