Differential gene expression patterns in vein regions susceptible versus resistant to neointimal hyperplasia.

Arteriovenous hemodialysis graft (AVG) stenosis results in thrombosis and AVG failure, but prevention of stenosis has been unsuccessful due in large part to our limited understanding of the molecular processes involved in neointimal hyperplasia (NH) formation. AVG stenosis develops chiefly as a consequence of highly localized NH formation in the vein-graft anastomosis region. Surprisingly, the vein region just downstream of the vein-graft anastomosis (herein termed proximal vein region) is relatively resistant to NH. We hypothesized that the gene expression profiles of the NH-prone and NH-resistant regions will be different from each other after graft placement, and analysis of their genomic profiles may yield potential therapeutic targets to prevent AVG stenosis. To test this, we evaluated the vein-graft anastomosis (NH-prone) and proximal vein (NH-resistant) regions in a porcine model of AVG stenosis with a porcine microarray. Gene expression changes in these two distinct vein regions, relative to the gene expression in unoperated control veins, were examined at early (5 days) and later (14 days) time points following graft placement. Global genomic changes were much greater in the NH-prone region than in the NH-resistant region at both time points. In the NH-prone region, genes related to regulation of cell proliferation and osteo-/chondrogenic vascular remodeling were most enriched among the significantly upregulated genes, and genes related to smooth muscle phenotype were significantly downregulated. These results provide insights into the spatial and temporal genomic modulation underlying NH formation in AVG and suggest potential therapeutic strategies to prevent and/or limit AVG stenosis.

Temporal Dysynchrony in brain connectivity gene expression following hypoxia.

BACKGROUND: Despite the fundamental biological importance and clinical relevance of characterizing the effects of chronic hypoxia exposure on central nervous system (CNS) development, the changes in gene expression from hypoxia are unknown. It is not known if there are unifying principles, properties, or logic in the response of the developing CNS to hypoxic exposure. Here, we use the small vertebrate zebrafish (Danio rerio) to study the effects of hypoxia on connectivity gene expression across development. We perform transcriptional profiling at high temporal resolution to systematically determine and then experimentally validate the response of CNS connectivity genes to hypoxia exposure. RESULTS: We characterized mRNA changes during development, comparing the effects of chronic hypoxia exposure at different time-points. We focused on changes in expression levels of a subset of 1270 genes selected for their roles in development of CNS connectivity, including axon pathfinding and synapse formation. We found that the majority of CNS connectivity genes were unaffected by hypoxia. However, for a small subset of genes hypoxia significantly affected their gene expression profiles. In particular, hypoxia appeared to affect both the timing and levels of expression, including altering expression of interacting gene pairs in a fashion that would potentially disrupt normal function. CONCLUSIONS: Overall, our study identifies the response of CNS connectivity genes to hypoxia exposure during development. While for most genes hypoxia did not significantly affect expression, for a subset of genes hypoxia changed both levels and timing of expression. Importantly, we identified that some genes with interacting proteins, for example receptor/ligand pairs, had dissimilar responses to hypoxia that would be expected to interfere with their function. The observed dysynchrony of gene expression could impair the development of normal CNS connectivity maps.

Ccdc94 protects cells from ionizing radiation by inhibiting the expression of p53.

DNA double-strand breaks (DSBs) represent one of the most deleterious forms of DNA damage to a cell. In cancer therapy, induction of cell death by DNA DSBs by ionizing radiation (IR) and certain chemotherapies is thought to mediate the successful elimination of cancer cells. However, cancer cells often evolve to evade the cytotoxicity induced by DNA DSBs, thereby forming the basis for treatment resistance. As such, a better understanding of the DSB DNA damage response (DSB-DDR) pathway will facilitate the design of more effective strategies to overcome chemo- and radioresistance. To identify novel mechanisms that protect cells from the cytotoxic effects of DNA DSBs, we performed a forward genetic screen in zebrafish for recessive mutations that enhance the IR-induced apoptotic response. Here, we describe radiosensitizing mutation 7 (rs7), which causes a severe sensitivity of zebrafish embryonic neurons to IR-induced apoptosis and is required for the proper development of the central nervous system. The rs7 mutation disrupts the coding sequence of ccdc94, a highly conserved gene that has no previous links to the DSB-DDR pathway. We demonstrate that Ccdc94 is a functional member of the Prp19 complex and that genetic knockdown of core members of this complex causes increased sensitivity to IR-induced apoptosis. We further show that Ccdc94 and the Prp19 complex protect cells from IR-induced apoptosis by repressing the expression of p53 mRNA. In summary, we have identified a new gene regulating a dosage-sensitive response to DNA DSBs during embryonic development. Future studies in human cancer cells will determine whether pharmacological inactivation of CCDC94 reduces the threshold of the cancer cell apoptotic response.

Zinc-dependent regulation of the Adh1 antisense transcript in fission yeast.

In yeast, Adh1 (alcohol dehydrogenase 1) is an abundant zinc-binding protein that is required for the conversion of acetaldehyde to ethanol. Through transcriptome profiling of the Schizosaccharomyces pombe genome, we identified a natural antisense transcript at the adh1 locus that is induced in response to zinc limitation. This antisense transcript (adh1AS) shows a reciprocal expression pattern to that of the adh1 mRNA partner. In this study, we show that increased expression of the adh1AS transcript in zinc-limited cells is necessary for the repression of adh1 gene expression and that the increased level of the adh1AS transcript in zinc-limited cells is a result of two mechanisms. At the transcriptional level, the adh1AS transcript is expressed at a high level in zinc-limited cells. In addition to this transcriptional control, adh1AS transcripts preferentially accumulate in zinc-limited cells when the adh1AS transcript is expressed from a constitutive promoter. This secondary mechanism requires the simultaneous expression of adh1. Our studies reveal how multiple mechanisms can synergistically control the ratio of sense to antisense transcripts and highlight a novel mechanism by which adh1 gene expression can be controlled by cellular zinc availability.

Chemical genetic screen reveals a role for desmosomal adhesion in mammary branching morphogenesis.

During the process of branching morphogenesis, the mammary gland undergoes distinct phases of remodeling to form an elaborate ductal network that ultimately produces and delivers milk to newborn animals. These developmental events rely on tight regulation of critical cellular pathways, many of which are probably disrupted during initiation and progression of breast cancer. Transgenic mouse and in vitro organoid models previously identified growth factor signaling as a key regulator of mammary branching, but the functional downstream targets of these pathways remain unclear. Here, we used purified primary mammary epithelial cells stimulated with fibroblast growth factor-2 (FGF2) to model mammary branching morphogenesis in vitro. We employed a forward chemical genetic approach to identify modulators of this process and describe a potent compound, 1023, that blocks FGF2-induced branching. In primary mammary epithelial cells, we used lentivirus-mediated knockdown of the aryl hydrocarbon receptor (AHR) to demonstrate that 1023 acts through AHR to block branching. Using 1023 as a tool, we identified desmosomal adhesion as a novel target of AHR signaling and show that desmosomes are critical for AHR agonists to block branching. Our findings support a functional role for desmosomes during mammary morphogenesis and also in blocking FGF-induced invasion.

Systematic generation of high-resolution deletion coverage of the Drosophila melanogaster genome.

In fruit fly research, chromosomal deletions are indispensable tools for mapping mutations, characterizing alleles and identifying interacting loci. Most widely used deletions were generated by irradiation or chemical mutagenesis. These methods are labor-intensive, generate random breakpoints and result in unwanted secondary mutations that can confound phenotypic analyses. Most of the existing deletions are large, have molecularly undefined endpoints and are maintained in genetically complex stocks. Furthermore, the existence of haplolethal or haplosterile loci makes the recovery of deletions of certain regions exceedingly difficult by traditional methods, resulting in gaps in coverage. Here we describe two methods that address these problems by providing for the systematic isolation of targeted deletions in the D. melanogaster genome. The first strategy used a P element-based technique to generate deletions that closely flank haploinsufficient genes and minimize undeleted regions. This deletion set has increased overall genomic coverage by 5-7%. The second strategy used FLP recombinase and the large array of FRT-bearing insertions described in the accompanying paper to generate 519 isogenic deletions with molecularly defined endpoints. This second deletion collection provides 56% genome coverage so far. The latter methodology enables the generation of small custom deletions with predictable endpoints throughout the genome and should make their isolation a simple and routine task.

A complementary transposon tool kit for Drosophila melanogaster using P and piggyBac.

With the availability of complete genome sequence for Drosophila melanogaster, one of the next strategic goals for fly researchers is a complete gene knockout collection. The P-element transposon, the workhorse of D. melanogaster molecular genetics, has a pronounced nonrandom insertion spectrum. It has been estimated that 87% saturation of the approximately 13,500-gene complement of D. melanogaster might require generating and analyzing up to 150,000 insertions. We describe specific improvements to the lepidopteran transposon piggyBac and the P element that enabled us to tag and disrupt genes in D. melanogaster more efficiently. We generated over 29,000 inserts resulting in 53% gene saturation and a more diverse collection of phenotypically stronger insertional alleles. We found that piggyBac has distinct global and local gene-tagging behavior from that of P elements. Notably, piggyBac excisions from the germ line are nearly always precise, piggyBac does not share chromosomal hotspots associated with P and piggyBac is more effective at gene disruption because it lacks the P bias for insertion in 5' regulatory sequences.

MicroRNAs and colon and rectal cancer: differential expression by tumor location and subtype.

MicroRNAs are thought to have an impact on cell proliferation, apoptosis, stress responses, maintenance of stem cell potency, and metabolism and are, therefore, important in the carcinogenic process. In this study, we examined 40 colon tumors, 30 rectal tumors, and 30 normal tissue samples (10 proximal colon, 10 distal colon, and 10 rectal paired with cancer cases) to examine miRNA expression profiles in colon and rectal tumors. MiRNA expression levels were adjusted for multiple comparisons; tumor tissue was compared with noncancerous tissue from the same site. A comparison of normal tissue showed 287 unique miRNAs that were significantly differentially expressed at the 1.5-fold level and 73 with over a two-fold difference in expression between colon and rectal tissue. Examination of miRNAs that were significantly differentially expressed at the 1.5-fold level by tumor phenotype showed 143 unique miRNAs differentially expression for microsatellite instability positive (MSI+) colon tumors; 129 unique miRNAs differentially expressed for CpG Island Methylator Phenotype positive (CIMP+) colon tumors; 135 miRNAs were differentially expressed for KRAS2-mutated colon tumors, and 139 miRNAs were differentially expressed for TP53-mutated colon tumors. Similar numbers of differentially expressed miRNAs were observed for rectal tumors, although the miRNAs differentially expressed differed. There were 129 unique miRNAs for CIMP+, 143 unique miRNAs for KRAS2-mutated, and 136 unique miRNAs for TP53-mutated rectal tumors. These results suggest the importance of miRNAs in colorectal cancer and the need for studies that can confirm these results and provide insight into the diet, lifestyle, and genetic factors that influence miRNA expression.

Forkhead proteins control the outcome of transcription factor binding by antiactivation.

Transcription factors with identical DNA-binding specificity often activate different genes in vivo. Yeast Ace2 and Swi5 are such activators, with targets we classify as Swi5-only, Ace2-only, or both. We define two unique regulatory modes. Ace2 and Swi5 both bind in vitro to Swi5-only genes such as HO, but only Swi5 binds and activates in vivo. In contrast, Ace2 and Swi5 both bind in vivo to Ace2-only genes, such as CTS1, but promoter-bound Swi5 fails to activate. We show that activation by Swi5 is prevented by the binding of the Forkhead factors Fkh1 and Fkh2, which recruit the Rpd3(Large) histone deacetylase complex to the CTS1 promoter. Global analysis shows that all Ace2-only genes are bound by both Ace2 and Swi5, and also by Fkh1/2. Genes normally activated by either Ace2 or Swi5 can be converted to Ace2-only genes by the insertion of Fkh-binding sites. Thus Fkh proteins, which function initially to activate SWI5 and ACE2, subsequently function as Swi5-specific antiactivators.

Next generation tools for genomic data generation, distribution, and visualization.

BACKGROUND: With the rapidly falling cost and availability of high throughput sequencing and microarray technologies, the bottleneck for effectively using genomic analysis in the laboratory and clinic is shifting to one of effectively managing, analyzing, and sharing genomic data. RESULTS: Here we present three open-source, platform independent, software tools for generating, analyzing, distributing, and visualizing genomic data. These include a next generation sequencing/microarray LIMS and analysis project center (GNomEx); an application for annotating and programmatically distributing genomic data using the community vetted DAS/2 data exchange protocol (GenoPub); and a standalone Java Swing application (GWrap) that makes cutting edge command line analysis tools available to those who prefer graphical user interfaces. Both GNomEx and GenoPub use the rich client Flex/Flash web browser interface to interact with Java classes and a relational database on a remote server. Both employ a public-private user-group security model enabling controlled distribution of patient and unpublished data alongside public resources. As such, they function as genomic data repositories that can be accessed manually or programmatically through DAS/2-enabled client applications such as the Integrated Genome Browser. CONCLUSIONS: These tools have gained wide use in our core facilities, research laboratories and clinics and are freely available for non-profit use. See http://sourceforge.net/projects/gnomex/, http://sourceforge.net/projects/genoviz/, and http://sourceforge.net/projects/useq.

Retrospective clinical trial experimentally validates glioblastoma genome-wide pattern of DNA copy-number alterations predictor of survival.

Modeling of genomic profiles from the Cancer Genome Atlas (TCGA) by using recently developed mathematical frameworks has associated a genome-wide pattern of DNA copy-number alterations with a shorter, roughly one-year, median survival time in glioblastoma (GBM) patients. Here, to experimentally test this relationship, we whole-genome sequenced DNA from tumor samples of patients. We show that the patients represent the U.S. adult GBM population in terms of most normal and disease phenotypes. Intratumor heterogeneity affects ≈ 11 % and profiling technology and reference human genome specifics affect <1% of the classifications of the tumors by the pattern, where experimental batch effects normally reduce the reproducibility, i.e., precision, of classifications based upon between one to a few hundred genomic loci by >30%. With a 2.25-year Kaplan-Meier median survival difference, a 3.5 univariate Cox hazard ratio, and a 0.78 concordance index, i.e., accuracy, the pattern predicts survival better than and independent of age at diagnosis, which has been the best indicator since 1950. The prognostic classification by the pattern may, therefore, help to manage GBM pseudoprogression. The diagnostic classification may help drugs progress to regulatory approval. The therapeutic predictions, of previously unrecognized targets that are correlated with survival, may lead to new drugs. Other methods missed this relationship in the roughly 3B-nucleotide genomes of the small, order of magnitude of 100, patient cohorts, e.g., from TCGA. Previous attempts to associate GBM genotypes with patient phenotypes were unsuccessful. This is a proof of principle that the frameworks are uniquely suitable for discovering clinically actionable genotype-phenotype relationships.

Global expression profiling identifies a novel biosignature for protein aggregation R120GCryAB cardiomyopathy in mice.

Protein aggregation cardiomyopathy is a life-threatening manifestation of a multisystem disorder caused by the exchange mutation in the gene encoding the human small heat shock protein alphaB-crystallin (hR120GCryAB). Genetic studies in mice have established cardiac hR120GCryAB expression causes increased activity of glucose 6-phosphate dehydrogenase (G6PD) and "reductive stress" (Rajasekaran et al., Cell 130: 427-439, 2007). However, the initiating molecular events in the pathogenesis of this novel toxic gain-of-function mechanism remain poorly defined. In an integrated systems approach using gene expression profiling, we identified a "biosignature," whose features can be validated to predict the onset, rate of progression, and clinical outcome of R120GCryAB cardiomyopathy. At the 3 mo disease-related but compensated stage, we demonstrate that transcripts were only upregulated in three distinct pathways: stress response (e.g., Hsp70, Hsp90), glutathione metabolism (Gpx1, Gpx3, glutathione S-transferase), and complement and coagulation cascades in hR120GCryAB transgenic mouse hearts compared with either hCryAB WT transgenic mice or nontransgenic controls. In 6 mo old myopathic hearts, ribosomal synthesis and cellular remodeling associated with increased cardiac hypertrophy were additional upregulated pathways. In contrast, the predominant downregulated pathways were for oxidative phosphorylation, fatty acid metabolism, intermediate metabolism, and energetic balance, supporting their primary pathogenic roles by which G6PD-dependent reductive stress causes cardiac decompensation and overt heart failure in hR120GCryAB cardiomyopathy. This study extends and confirms our previous findings that reductive stress is a causal mechanism for hR120G CryAB cardiomyopathy and demonstrates that alteration in glutathione pathway gene expression is an early biosignature with utility for presymptomatic detection.

Radiosensitization of cervical cancer cells via double-strand DNA break repair inhibition.

PURPOSE: LY294002, a phosphatidylinositol 3-kinase (PI3K) inhibitor, has been found to radiosensitize various human cancer cells. However, its potential to act as an effective therapeutic agent is diminished by its toxicity levels. The purposes of this study were to determine the mechanism by which LY294002 radiosensitizes. MATERIALS AND METHODS: Cell growth curves and clonogenic assays were performed with increasing LY294002 exposure times proximate to the radiation dose. Protein levels of downstream PI3K effectors were analyzed. Detection of phosphorylated histone H2AX (gammaH2AX) was used to identify DNA double-strand breaks at various time points post-radiation. RESULTS: LY294002 significantly radiosensitized HeLa cervical cancer cells when administered for just 12 h following radiation. Cell growth curves also decreased with brief LY294002 application. DNA double-strand breaks are typically repaired within 2-6 h following radiation. Interestingly, at 48, 72, and 96 h post-irradiation, gammaH2AX was still significantly elevated in cells radiated in combination with LY294002. Protein expressions of ATM and ATR downstream effectors showed no differences among the treated groups, however, DNA-PK activity was significantly inhibited by LY294002. CONCLUSIONS: These results lead us to conclude that the central mechanism by which LY294002 radiosensitizes is via DNA-PK inhibition which induces DNA double-strand break repair inhibition. We are currently investigating radiosensitization induced by DNA-PK-specific inhibition in efforts to find a less toxic, yet equally effective, chemotherapeutic agent than LY294002.

Gene expression signatures in the newt irises during lens regeneration.

Lens regeneration in adult newts is possible by transdifferentiation of the pigment epithelial cells (PECs) of the dorsal iris. The same cells in the ventral iris are not capable of such a process. To understand this difference in regenerative competency, we examined gene expression of 373 genes in the intact dorsal and ventral irises as well as in irises during the process of lens regeneration. We found similar signatures of gene expression in dorsal and ventral with several cases of even higher levels in the ventral iris. Such transcriptional activity in the regeneration-incompetent ventral iris was unexpected and calls for a revision of our views about mechanisms of lens regeneration induction.

The kinetics and redox state of nitric oxide determine the biological consequences in lung adenocarcinoma.

Few studies have explored the mechanistic basis for the apparent paradoxical effects of nitric oxide and its interrelated redox species (NO(X)) in cancer biology. Our aim was to determine the differential effects of the redox state and kinetics of nitrosative species on the key cancer processes of apoptosis. Therefore, a murine lung adenocarcinoma cell line was exposed to various NO(X) donor compounds differing in redox state and delivery kinetics. DNA strand breaks (DSBs) were measured by the alkaline single-cell gel electrophoresis assay (the COMET assay) and correlated with cell viability by the MTT and soft agar colony assays, while caspase enzymatic activity was measured using an in vitro fluorogenic caspase assay. Finally, cDNA microarrays defined apoptosis-related gene expression alterations resultant from these NO(X) donors. Exogenous NO(X) differentially influences DSBs, and apoptosis-related cell death and expression based on the redox state and kinetics of NO(X) delivery. In our murine lung adenocarcinoma model we have demonstrated differential effects of NO(X) based on the mode of delivery and redox state. These data suggest that the development of NO(X)-based cancer chemotherapy must consider the redox state and kinetics of delivery into their logical design.

Spontaneous shaker rat mutant - a new model for X-linked tremor/ataxia.

The shaker rat is an X-linked recessive spontaneous model of progressive Purkinje cell (PC) degeneration exhibiting a shaking ataxia and wide stance. Generation of Wistar Furth (WF)/Brown Norwegian (BN) F1 hybrids and genetic mapping of F2 sib-sib offspring using polymorphic markers narrowed the candidate gene region to 26 Mbp denoted by the last recombinant genetic marker DXRat21 at 133 Mbp to qter (the end of the long arm). In the WF background, the shaker mutation has complete penetrance, results in a stereotypic phenotype and there is a narrow window for age of disease onset; by contrast, the F2 hybrid phenotype was more varied, with a later age of onset and likely non-penetrance of the mutation. By deep RNA-sequencing, five variants were found in the candidate region; four were novel without known annotation. One of the variants caused an arginine (R) to cysteine (C) change at codon 35 of the ATPase, Ca(2+) transporting, plasma membrane 3 (Atp2b3) gene encoding PMCA3 that has high expression in the cerebellum. The variant was well supported by hundreds of overlapping reads, and was found in 100% of all affected replicas and 0% of the wild-type (WT) replicas. The mutation segregated with disease in all affected animals and the amino acid change was found in an evolutionarily conserved region of PMCA3. Despite strong genetic evidence for pathogenicity, in vitro analyses of PMCA3(R35C) function did not show any differences to WT PMCA3. Because Atp2b3 mutation leads to congenital ataxia in humans, the identified Atp2b3 missense change in the shaker rat presents a good candidate for the shaker rat phenotype based on genetic criteria, but cannot yet be considered a definite pathogenic variant owing to lack of functional changes.

Multiple abnormalities detected by dye reversal genomic microarrays in prostate cancer: a much greater sensitivity than conventional cytogenetics.

Prostate cancer remains the most common male malignancy in Western countries, yet limited information exists regarding genetic changes and clinical correlations. The advent of comparative genomic hybridization microarray (GM) technology has recently allowed for precise screening of DNAs for genetic copy number changes; this offers an advantage over previous techniques, including conventional cytogenetics. A problem with cytogenetic prostate cancer analysis has been the study of the appropriate cell types because this is a highly heterogeneous tumor. We have performed GM using the Spectral Genomics Inc. dye reversal platform on 20 primary prostate tumors. These tumor samples were from frozen tissue collected over the last 10 years and multiple clinical parameters, including follow-up were collected on these patients; cytogenetic analysis was previously attempted on all patients. Eighty percent (16/20) of specimens showed copy number changes, 65% of which were losses and 35% were gains of genetic material. The most common changes observed were loss of an interstitial region of 2q (8 cases, 40%), followed by loss of interstitial 6q (6 cases, 30%), loss at 8p and 13q (5 cases each, 25%), gain at 3p and loss at 5q, 16q, and Xq (4 cases each, 20%), and gain at 8p (3 cases, 15%). There was evidence of correlation of loss at 5q with a positive node status. Cytogenetic studies on these same patients only detected clonal changes in 40% (8/20) specimens and did not detect the majority of abnormalities seen by the GM technique. We propose this technology for the evaluation of prostate and other heterogeneous cancers as a rapid and efficient way to detect genetic copy number changes.

Novel venom peptides from the cone snail Conus pulicarius discovered through next-generation sequencing of its venom duct transcriptome.

The venom peptides (i.e., conotoxins or conopeptides) that species in the genus Conus collectively produce are remarkably diverse, estimated to be around 50,000 to 140,000, but the pace of discovery and characterization of these peptides have been rather slow. To date, only a minor fraction have been identified and studied. However, the advent of next-generation DNA sequencing technologies has opened up opportunities for expediting the exploration of this diversity. The whole transcriptome of a venom duct from the vermivorous marine snail C. pulicarius was sequenced using the 454 sequencing platform. Analysis of the data set resulted in the identification of over eighty unique putative conopeptide sequences, the highest number discovered so far from a Conus venom duct transcriptome. More importantly, majority of the sequences were potentially novel, many with unexpected structural features, hinting at the vastness of the diversity of Conus venom peptides that remains to be explored. The sequences represented at least 14 major superfamilies/types (disulfide- and non-disulfide-rich), indicating the structural and functional diversity of conotoxins in the venom of C. pulicarius. In addition, the contryphans were surprisingly more diverse than what is currently known. Comparative analysis of the O-superfamily sequences also revealed insights into the complexity of the processes that drive the evolution and diversification of conotoxins.