query: De novo assembly of sequence data has fewer specific contigs than unassembled sequence data.

title: A culture-independent sequence-based metagenomics approach to the investigation of an outbreak of Shiga-toxigenic Escherichia coli O104:H4.

text:
IMPORTANCE Identification of the bacterium responsible for an outbreak can aid in disease management. However, traditional culture-based diagnosis can be difficult, particularly if no specific diagnostic test is available for an outbreak strain.
 OBJECTIVE To explore the potential of metagenomics, which is the direct sequencing of DNA extracted from microbiologically complex samples, as an open-ended clinical discovery platform capable of identifying and characterizing bacterial strains from an outbreak without laboratory culture.
 DESIGN, SETTING, AND PATIENTS In a retrospective investigation, 45 samples were selected from fecal specimens obtained from patients with diarrhea during the 2011 outbreak of Shiga-toxigenic Escherichia coli (STEC) O104:H4 in Germany. Samples were subjected to high-throughput sequencing (August-September 2012), followed by a 3-phase analysis (November 2012-February 2013). In phase 1, a de novo assembly approach was developed to obtain a draft genome of the outbreak strain. In phase 2, the depth of coverage of the outbreak strain genome was determined in each sample. In phase 3, sequences from each sample were compared with sequences from known bacteria to identify pathogens other than the outbreak strain.
 MAIN OUTCOMES AND MEASURES The recovery of genome sequence data for the purposes of identification and characterization of the outbreak strain and other pathogens from fecal samples.
 RESULTS During phase 1, a draft genome of the STEC outbreak strain was obtained. During phase 2, the outbreak strain genome was recovered from 10 samples at greater than 10-fold coverage and from 26 samples at greater than 1-fold coverage. Sequences from the Shiga-toxin genes were detected in 27 of 40 STEC-positive samples (67%). In phase 3, sequences from Clostridium difficile, Campylobacter jejuni, Campylobacter concisus, and Salmonella enterica were recovered.
 CONCLUSIONS AND RELEVANCE These results suggest the potential of metagenomics as a culture-independent approach for the identification of bacterial pathogens during an outbreak of diarrheal disease. Challenges include improving diagnostic sensitivity, speeding up and simplifying workflows, and reducing costs.

title: Effects of GC Bias in Next-Generation-Sequencing Data on De Novo Genome Assembly

text: Next-generation-sequencing (NGS) has revolutionized the field of genome assembly because of its much higher data throughput and much lower cost compared with traditional Sanger sequencing. However, NGS poses new computational challenges to de novo genome assembly. Among the challenges, GC bias in NGS data is known to aggravate genome assembly. However, it is not clear to what extent GC bias affects genome assembly in general. In this work, we conduct a systematic analysis on the effects of GC bias on genome assembly. Our analyses reveal that GC bias only lowers assembly completeness when the degree of GC bias is above a threshold. At a strong GC bias, the assembly fragmentation due to GC bias can be explained by the low coverage of reads in the GC-poor or GC-rich regions of a genome. This effect is observed for all the assemblers under study. Increasing the total amount of NGS data thus rescues the assembly fragmentation because of GC bias. However, the amount of data needed for a full rescue depends on the distribution of GC contents. Both low and high coverage depths due to GC bias lower the accuracy of assembly. These pieces of information provide guidance toward a better de novo genome assembly in the presence of GC bias.

title: The miniature genome of a carnivorous plant Genlisea aurea contains a low number of genes and short non-coding sequences

text:
BACKGROUND Genlisea aurea (Lentibulariaceae) is a carnivorous plant with unusually small genome size - 63.6 Mb - one of the smallest known among higher plants. Data on the genome sizes and the phylogeny of Genlisea suggest that this is a derived state within the genus. Thus, G. aurea is an excellent model organism for studying evolutionary mechanisms of genome contraction.
 RESULTS Here we report sequencing and de novo draft assembly of G. aurea genome. The assembly consists of 10,687 contigs of the total length of 43.4 Mb and includes 17,755 complete and partial protein-coding genes. Its comparison with the genome of Mimulus guttatus, another representative of higher core Lamiales clade, reveals striking differences in gene content and length of non-coding regions.
 CONCLUSIONS Genome contraction was a complex process, which involved gene loss and reduction of lengths of introns and intergenic regions, but not intron loss. The gene loss is more frequent for the genes that belong to multigenic families indicating that genetic redundancy is an important prerequisite for genome size reduction.

title: RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome

text:
BACKGROUND RNA-Seq is revolutionizing the way transcript abundances are measured. A key challenge in transcript quantification from RNA-Seq data is the handling of reads that map to multiple genes or isoforms. This issue is particularly important for quantification with de novo transcriptome assemblies in the absence of sequenced genomes, as it is difficult to determine which transcripts are isoforms of the same gene. A second significant issue is the design of RNA-Seq experiments, in terms of the number of reads, read length, and whether reads come from one or both ends of cDNA fragments.
 RESULTS We present RSEM, an user-friendly software package for quantifying gene and isoform abundances from single-end or paired-end RNA-Seq data. RSEM outputs abundance estimates, 95% credibility intervals, and visualization files and can also simulate RNA-Seq data. In contrast to other existing tools, the software does not require a reference genome. Thus, in combination with a de novo transcriptome assembler, RSEM enables accurate transcript quantification for species without sequenced genomes. On simulated and real data sets, RSEM has superior or comparable performance to quantification methods that rely on a reference genome. Taking advantage of RSEM's ability to effectively use ambiguously-mapping reads, we show that accurate gene-level abundance estimates are best obtained with large numbers of short single-end reads. On the other hand, estimates of the relative frequencies of isoforms within single genes may be improved through the use of paired-end reads, depending on the number of possible splice forms for each gene.
 CONCLUSIONS RSEM is an accurate and user-friendly software tool for quantifying transcript abundances from RNA-Seq data. As it does not rely on the existence of a reference genome, it is particularly useful for quantification with de novo transcriptome assemblies. In addition, RSEM has enabled valuable guidance for cost-efficient design of quantification experiments with RNA-Seq, which is currently relatively expensive.

title: Human artificial chromosomes with alpha satellite-based de novo centromeres show increased frequency of nondisjunction and anaphase lag.

text: Human artificial chromosomes have been used to model requirements for human chromosome segregation and to explore the nature of sequences competent for centromere function. Normal human centromeres require specialized chromatin that consists of alpha satellite DNA complexed with epigenetically modified histones and centromere-specific proteins. While several types of alpha satellite DNA have been used to assemble de novo centromeres in artificial chromosome assays, the extent to which they fully recapitulate normal centromere function has not been explored. Here, we have used two kinds of alpha satellite DNA, DXZ1 (from the X chromosome) and D17Z1 (from chromosome 17), to generate human artificial chromosomes. Although artificial chromosomes are mitotically stable over many months in culture, when we examined their segregation in individual cell divisions using an anaphase assay, artificial chromosomes exhibited more segregation errors than natural human chromosomes (P < 0.001). Naturally occurring, but abnormal small ring chromosomes derived from chromosome 17 and the X chromosome also missegregate more than normal chromosomes, implicating overall chromosome size and/or structure in the fidelity of chromosome segregation. As different artificial chromosomes missegregate over a fivefold range, the data suggest that variable centromeric DNA content and/or epigenetic assembly can influence the mitotic behavior of artificial chromosomes.

title: Full-length transcriptome assembly from RNA-Seq data without a reference genome.

text: Massively parallel sequencing of cDNA has enabled deep and efficient probing of transcriptomes. Current approaches for transcript reconstruction from such data often rely on aligning reads to a reference genome, and are thus unsuitable for samples with a partial or missing reference genome. Here we present the Trinity method for de novo assembly of full-length transcripts and evaluate it on samples from fission yeast, mouse and whitefly, whose reference genome is not yet available. By efficiently constructing and analyzing sets of de Bruijn graphs, Trinity fully reconstructs a large fraction of transcripts, including alternatively spliced isoforms and transcripts from recently duplicated genes. Compared with other de novo transcriptome assemblers, Trinity recovers more full-length transcripts across a broad range of expression levels, with a sensitivity similar to methods that rely on genome alignments. Our approach provides a unified solution for transcriptome reconstruction in any sample, especially in the absence of a reference genome.

title: Whole-genome sequencing of Oryza brachyantha reveals mechanisms underlying Oryza genome evolution

text: The wild species of the genus Oryza contain a largely untapped reservoir of agronomically important genes for rice improvement. Here we report the 261-Mb de novo assembled genome sequence of Oryza brachyantha. Low activity of long-terminal repeat retrotransposons and massive internal deletions of ancient long-terminal repeat elements lead to the compact genome of Oryza brachyantha. We model 32,038 protein-coding genes in the Oryza brachyantha genome, of which only 70% are located in collinear positions in comparison with the rice genome. Analysing breakpoints of non-collinear genes suggests that double-strand break repair through non-homologous end joining has an important role in gene movement and erosion of collinearity in the Oryza genomes. Transition of euchromatin to heterochromatin in the rice genome is accompanied by segmental and tandem duplications, further expanded by transposable element insertions. The high-quality reference genome sequence of Oryza brachyantha provides an important resource for functional and evolutionary studies in the genus Oryza.

title: BLAST+: architecture and applications

text:
BACKGROUND Sequence similarity searching is a very important bioinformatics task. While Basic Local Alignment Search Tool (BLAST) outperforms exact methods through its use of heuristics, the speed of the current BLAST software is suboptimal for very long queries or database sequences. There are also some shortcomings in the user-interface of the current command-line applications.
 RESULTS We describe features and improvements of rewritten BLAST software and introduce new command-line applications. Long query sequences are broken into chunks for processing, in some cases leading to dramatically shorter run times. For long database sequences, it is possible to retrieve only the relevant parts of the sequence, reducing CPU time and memory usage for searches of short queries against databases of contigs or chromosomes. The program can now retrieve masking information for database sequences from the BLAST databases. A new modular software library can now access subject sequence data from arbitrary data sources. We introduce several new features, including strategy files that allow a user to save and reuse their favorite set of options. The strategy files can be uploaded to and downloaded from the NCBI BLAST web site.
 CONCLUSION The new BLAST command-line applications, compared to the current BLAST tools, demonstrate substantial speed improvements for long queries as well as chromosome length database sequences. We have also improved the user interface of the command-line applications.

title: CENP-B Controls Centromere Formation Depending on the Chromatin Context

text: The centromere is a chromatin region that serves as the spindle attachment point and directs accurate inheritance of eukaryotic chromosomes during cell divisions. However, the mechanism by which the centromere assembles and stabilizes at a specific genomic region is not clear. The de novo formation of a human/mammalian artificial chromosome (HAC/MAC) with a functional centromere assembly requires the presence of alpha-satellite DNA containing binding motifs for the centromeric CENP-B protein. We demonstrate here that de novo centromere assembly on HAC/MAC is dependent on CENP-B. In contrast, centromere formation is suppressed in cells expressing CENP-B when alpha-satellite DNA was integrated into a chromosomal site. Remarkably, on those integration sites CENP-B enhances histone H3-K9 trimethylation and DNA methylation, thereby stimulating heterochromatin formation. Thus, we propose that CENP-B plays a dual role in centromere formation, ensuring de novo formation on DNA lacking a functional centromere but preventing the formation of excess centromeres on chromosomes.

title: Metagenomic analysis of microbial consortium from natural crude oil that seeps into the marine ecosystem offshore Southern California

text: Crude oils can be major contaminants of the marine ecosystem and microorganisms play a significant role in the degradation of its main constituents. To increase our understanding of the microbial hydrocarbon degradation process in the marine ecosystem, we collected crude oil from an active seep area located in the Santa Barbara Channel (SBC) and generated a total of about 52 Gb of raw metagenomic sequence data. The assembled data comprised ~500 Mb, representing ~1.1 million genes derived primarily from chemolithoautotrophic bacteria. Members of Oceanospirillales, a bacterial order belonging to the Deltaproteobacteria, recruited less than 2% of the assembled genes within the SBC metagenome. In contrast, the microbial community associated with the oil plume that developed in the aftermath of the Deepwater Horizon (DWH) blowout in 2010, was dominated by Oceanospirillales, which comprised more than 60% of the metagenomic data generated from the DWH oil plume. This suggests that Oceanospirillales might play a less significant role in the microbially mediated hydrocarbon conversion within the SBC seep oil compared to the DWH plume oil. We hypothesize that this difference results from the SBC oil seep being mostly anaerobic, while the DWH oil plume is aerobic. Within the Archaea, the phylum Euryarchaeota, recruited more than 95% of the assembled archaeal sequences from the SBC oil seep metagenome, with more than 50% of the sequences assigned to members of the orders Methanomicrobiales and Methanosarcinales. These orders contain organisms capable of anaerobic methanogenesis and methane oxidation (AOM) and we hypothesize that these orders - and their metabolic capabilities - may be fundamental to the ecology of the SBC oil seep.

title: Junctional actin assembly is mediated by Formin-like 2 downstream of Rac1

text: Epithelial integrity is vitally important, and its deregulation causes early stage cancer. De novo formation of an adherens junction (AJ) between single epithelial cells requires coordinated, spatial actin dynamics, but the mechanisms steering nascent actin polymerization for cell-cell adhesion initiation are not well understood. Here we investigated real-time actin assembly during daughter cell-cell adhesion formation in human breast epithelial cells in 3D environments. We identify formin-like 2 (FMNL2) as being specifically required for actin assembly and turnover at newly formed cell-cell contacts as well as for human epithelial lumen formation. FMNL2 associates with components of the AJ complex involving Rac1 activity and the FMNL2 C terminus. Optogenetic control of Rac1 in living cells rapidly drove FMNL2 to epithelial cell-cell contact zones. Furthermore, Rac1-induced actin assembly and subsequent AJ formation critically depends on FMNL2. These data uncover FMNL2 as a driver for human epithelial AJ formation downstream of Rac1.

title: STAR: ultrafast universal RNA-seq aligner.

text:
MOTIVATION Accurate alignment of high-throughput RNA-seq data is a challenging and yet unsolved problem because of the non-contiguous transcript structure, relatively short read lengths and constantly increasing throughput of the sequencing technologies. Currently available RNA-seq aligners suffer from high mapping error rates, low mapping speed, read length limitation and mapping biases.
 RESULTS To align our large (>80 billon reads) ENCODE Transcriptome RNA-seq dataset, we developed the Spliced Transcripts Alignment to a Reference (STAR) software based on a previously undescribed RNA-seq alignment algorithm that uses sequential maximum mappable seed search in uncompressed suffix arrays followed by seed clustering and stitching procedure. STAR outperforms other aligners by a factor of >50 in mapping speed, aligning to the human genome 550 million 2 × 76 bp paired-end reads per hour on a modest 12-core server, while at the same time improving alignment sensitivity and precision. In addition to unbiased de novo detection of canonical junctions, STAR can discover non-canonical splices and chimeric (fusion) transcripts, and is also capable of mapping full-length RNA sequences. Using Roche 454 sequencing of reverse transcription polymerase chain reaction amplicons, we experimentally validated 1960 novel intergenic splice junctions with an 80-90% success rate, corroborating the high precision of the STAR mapping strategy.    AVAILABILITY AND IMPLEMENTATION STAR is implemented as a standalone C++ code. STAR is free open source software distributed under GPLv3 license and can be downloaded from http://code.google.com/p/rna-star/.

title: Direct-coupling analysis of residue co-evolution captures native contacts across many protein families

text: The similarity in the three-dimensional structures of homologous proteins imposes strong constraints on their sequence variability. It has long been suggested that the resulting correlations among amino acid compositions at different sequence positions can be exploited to infer spatial contacts within the tertiary protein structure. Crucial to this inference is the ability to disentangle direct and indirect correlations, as accomplished by the recently introduced Direct Coupling Analysis (DCA) (Weigt et al. (2009) Proc Natl Acad Sci 106:67). Here we develop a computationally efficient implementation of DCA, which allows us to evaluate the accuracy of contact prediction by DCA for a large number of protein domains, based purely on sequence information. DCA is shown to yield a large number of correctly predicted contacts, recapitulating the global structure of the contact map for the majority of the protein domains examined. Furthermore, our analysis captures clear signals beyond intra- domain residue contacts, arising, e.g., from alternative protein conformations, ligand- mediated residue couplings, and inter-domain interactions in protein oligomers. Our findings suggest that contacts predicted by DCA can be used as a reliable guide to facilitate computational predictions of alternative protein conformations, protein complex formation, and even the de novo prediction of protein domain structures, provided the existence of a large number of homologous sequences which are being rapidly made available due to advances in genome sequencing.

title: Nucleation of nuclear bodies by RNA

text: The biogenesis of the many functional compartments contained in the mammalian cell nucleus is poorly understood. More specifically, little is known regarding the initial nucleation step required for nuclear body formation. Here we show that RNA can function as a structural element and a nucleator of nuclear bodies. We find that several types of coding and noncoding RNAs are sufficient to de novo assemble, and are physiologically enriched in, histone locus bodies (with associated Cajal bodies), nuclear speckles, paraspeckles and nuclear stress bodies. Formation of nuclear bodies occurs through recruitment and accumulation of proteins resident in the nuclear bodies by nucleating RNA. These results demonstrate that transcription is a driving force in nuclear body formation and RNA transcripts can function as a scaffold in the formation of major nuclear bodies. Together, these data suggest that RNA-primed biogenesis of nuclear bodies is a general principle of nuclear organization.

title: Detecting selective sweeps: a new approach based on hidden markov models.

text: Detecting and localizing selective sweeps on the basis of SNP data has recently received considerable attention. Here we introduce the use of hidden Markov models (HMMs) for the detection of selective sweeps in DNA sequences. Like previously published methods, our HMMs use the site frequency spectrum, and the spatial pattern of diversity along the sequence, to identify selection. In contrast to earlier approaches, our HMMs explicitly model the correlation structure between linked sites. The detection power of our methods, and their accuracy for estimating the selected site location, is similar to that of competing methods for constant size populations. In the case of population bottlenecks, however, our methods frequently showed fewer false positives.

title: Effects of an opal termination codon preceding the nsP4 gene sequence in the O'Nyong-Nyong virus genome on Anopheles gambiae infectivity.

text: The genomic RNA of an alphavirus encodes four different nonstructural proteins, nsP1, nsP2, nsP3, and nsP4. The polyprotein P123 is produced when translation terminates at an opal termination codon between nsP3 and nsP4. The polyprotein P1234 is produced when translational readthrough occurs or when the opal termination codon has been replaced by a sense codon in the alphavirus genome. Evolutionary pressures appear to have maintained genomic sequences encoding both a stop codon (opal) and an open reading frame (arginine) as a general feature of the O'nyong-nyong virus (ONNV) genome, indicating that both are required at some point. Alternate replication of ONNVs in both vertebrate and invertebrate hosts may determine predominance of a particular codon at this locus in the viral quasispecies. However, no systematic study has previously tested this hypothesis in whole animals. We report here the results of the first study to investigate in a natural mosquito host the functional significance of the opal stop codon in an alphavirus genome. We used a full-length cDNA clone of ONNV to construct a series of mutants in which the arginine between nsP3 and nsP4 was replaced with an opal, ochre, or amber stop codon. The presence of an opal stop codon upstream of nsP4 nearly doubled (75.5%) the infectivity of ONNV over that of virus possessing a codon for the amino acid arginine at the corresponding position (39.8%). Although the frequency with which the opal virus disseminated from the mosquito midgut did not differ significantly from that of the arginine virus on days 8 and 10, dissemination did began earlier in mosquitoes infected with the opal virus. Although a clear fitness advantage is provided to ONNV by the presence of an opal codon between nsP3 and nsP4 in Anopheles gambiae, sequence analysis of ONNV RNA extracted from mosquito bodies and heads indicated codon usage at this position corresponded with that of the virus administered in the blood meal. These results suggest that while selection of ONNV variants is occurring, de novo mutation at the position between nsP3 and nsP4 does not readily occur in the mosquito. Taken together, these results suggest that the primary fitness advantage provided to ONNV by the presence of an opal codon between nsP3 and nsP4 is related to mosquito infectivity.

title: Circulating tumor DNA as a non-invasive substitute to metastasis biopsy for tumor genotyping and personalized medicine in a prospective trial across all tumor types.

text: Cell-free tumor DNA (ctDNA) has the potential to enable non-invasive diagnostic tests for personalized medicine in providing similar molecular information as that derived from invasive tumor biopsies. The histology-independent phase II SHIVA trial matches patients with targeted therapeutics based on previous screening of multiple somatic mutations using metastatic biopsies. To evaluate the utility of ctDNA in this trial, as an ancillary study we performed de novo detection of somatic mutations using plasma DNA compared to metastasis biopsies in 34 patients covering 18 different tumor types, scanning 46 genes and more than 6800 COSMIC mutations with a multiplexed next-generation sequencing panel. In 27 patients, 28 of 29 mutations identified in metastasis biopsies (97%) were detected in matched ctDNA. Among these 27 patients, one additional mutation was found in ctDNA only. In the seven other patients, mutation detection from metastasis biopsy failed due to inadequate biopsy material, but was successful in all plasma DNA samples providing three more potential actionable mutations. These results suggest that ctDNA analysis is a potential alternative and/or replacement to analyses using costly, harmful and lengthy tissue biopsies of metastasis, irrespective of cancer type and metastatic site, for multiplexed mutation detection in selecting personalized therapies based on the patient's tumor genetic content.

title: Production of de novo cardiomyocytes: human pluripotent stem cell differentiation and direct reprogramming.

text: Cardiovascular disease is a leading cause of death worldwide. The limited capability of heart tissue to regenerate has prompted methodological developments for creating de novo cardiomyocytes, both in vitro and in vivo. Beyond uses in cell replacement therapy, patient-specific cardiomyocytes may find applications in drug testing, drug discovery, and disease modeling. Recently, approaches for generating cardiomyocytes have expanded to encompass three major sources of starting cells: human pluripotent stem cells (hPSCs), adult heart-derived cardiac progenitor cells (CPCs), and reprogrammed fibroblasts. We discuss state-of-the-art methods for generating de novo cardiomyocytes from hPSCs and reprogrammed fibroblasts, highlighting potential applications and future challenges.

title: Human artificial chromosomes containing chromosome 17 alphoid DNA maintain an active centromere in murine cells but are not stable.

text: Human artificial chromosomes (HACs) are autonomous molecules that can function and segregate as normal chromosomes in human cells. De novo HACs have successfully been used as gene expression vectors to complement genetic deficiencies in human cultured cells. HACs now offer the possibility of studying the regulation and expression of large genes in a variety of cell types from different tissues and correcting gene deficiencies caused by human inherited diseases. Complementary gene expression studies in mice, especially in mouse models of human genetic diseases, are also important in determining if large human transgenes can be expressed appropriately from artificial chromosomes. Toward this aim we are establishing artificial chromosomes in murine cells as novel gene expression vectors. Initially we transferred HAC vectors into murine cells, but were unable to generate de novo HACs at a reasonable frequency. We then transferred HACs previously established in human HT1080 cells to three different murine cell types by microcell fusion, followed by positive selection. We observed that the HACs in murine cells bound centromere protein C (CENP-C), a marker of active centromeres, and were detected under selection but rapidly lost when selection was removed. These results suggest that the HACs maintain at least a partially functional centromere complex in murine cells, but other factors are required for stability and segregation. Artificial chromosomes containing mouse centromeric sequences may be required for better stability and maintenance in murine cells.

title: Expression profiles during dedifferentiation in newt lens regeneration revealed by expressed sequence tags

text:
PURPOSE The adult newt can regenerate lens from pigmented epithelial cells (PECs) of the dorsal iris via dedifferentiation. The purpose of this research is to obtain sequence resources for a newt lens regeneration study and to obtain insights of dedifferentiation at the molecular level.
 METHODS mRNA was purified from iris during dedifferentiation and its cDNA library was constructed. From the cDNA library 10,449 clones were sequenced and analyzed.
 RESULTS From 10,449 reads, 780 contigs and 1,666 singlets were annotated. The presence of several cancer- and apoptosis-related genes during newt dedifferentiation was revealed. Moreover, several candidate genes, which might participate in reprogramming during dedifferentiation, were also found.
 CONCLUSIONS The expression of cancer- and apoptosis-related genes could be hallmarks during dedifferentiation. The expression sequence tag (EST) resource is useful for the future study of newt dedifferentiation, and the sequence information is available in GenBank (accession numbers; FS290155-FS300559).


query: De novo assembly of sequence data has shorter contigs than unassembled sequence data.

title: A culture-independent sequence-based metagenomics approach to the investigation of an outbreak of Shiga-toxigenic Escherichia coli O104:H4.

text:
IMPORTANCE Identification of the bacterium responsible for an outbreak can aid in disease management. However, traditional culture-based diagnosis can be difficult, particularly if no specific diagnostic test is available for an outbreak strain.
 OBJECTIVE To explore the potential of metagenomics, which is the direct sequencing of DNA extracted from microbiologically complex samples, as an open-ended clinical discovery platform capable of identifying and characterizing bacterial strains from an outbreak without laboratory culture.
 DESIGN, SETTING, AND PATIENTS In a retrospective investigation, 45 samples were selected from fecal specimens obtained from patients with diarrhea during the 2011 outbreak of Shiga-toxigenic Escherichia coli (STEC) O104:H4 in Germany. Samples were subjected to high-throughput sequencing (August-September 2012), followed by a 3-phase analysis (November 2012-February 2013). In phase 1, a de novo assembly approach was developed to obtain a draft genome of the outbreak strain. In phase 2, the depth of coverage of the outbreak strain genome was determined in each sample. In phase 3, sequences from each sample were compared with sequences from known bacteria to identify pathogens other than the outbreak strain.
 MAIN OUTCOMES AND MEASURES The recovery of genome sequence data for the purposes of identification and characterization of the outbreak strain and other pathogens from fecal samples.
 RESULTS During phase 1, a draft genome of the STEC outbreak strain was obtained. During phase 2, the outbreak strain genome was recovered from 10 samples at greater than 10-fold coverage and from 26 samples at greater than 1-fold coverage. Sequences from the Shiga-toxin genes were detected in 27 of 40 STEC-positive samples (67%). In phase 3, sequences from Clostridium difficile, Campylobacter jejuni, Campylobacter concisus, and Salmonella enterica were recovered.
 CONCLUSIONS AND RELEVANCE These results suggest the potential of metagenomics as a culture-independent approach for the identification of bacterial pathogens during an outbreak of diarrheal disease. Challenges include improving diagnostic sensitivity, speeding up and simplifying workflows, and reducing costs.

title: Effects of GC Bias in Next-Generation-Sequencing Data on De Novo Genome Assembly

text: Next-generation-sequencing (NGS) has revolutionized the field of genome assembly because of its much higher data throughput and much lower cost compared with traditional Sanger sequencing. However, NGS poses new computational challenges to de novo genome assembly. Among the challenges, GC bias in NGS data is known to aggravate genome assembly. However, it is not clear to what extent GC bias affects genome assembly in general. In this work, we conduct a systematic analysis on the effects of GC bias on genome assembly. Our analyses reveal that GC bias only lowers assembly completeness when the degree of GC bias is above a threshold. At a strong GC bias, the assembly fragmentation due to GC bias can be explained by the low coverage of reads in the GC-poor or GC-rich regions of a genome. This effect is observed for all the assemblers under study. Increasing the total amount of NGS data thus rescues the assembly fragmentation because of GC bias. However, the amount of data needed for a full rescue depends on the distribution of GC contents. Both low and high coverage depths due to GC bias lower the accuracy of assembly. These pieces of information provide guidance toward a better de novo genome assembly in the presence of GC bias.

text:
BACKGROUND Genlisea aurea (Lentibulariaceae) is a carnivorous plant with unusually small genome size - 63.6 Mb - one of the smallest known among higher plants. Data on the genome sizes and the phylogeny of Genlisea suggest that this is a derived state within the genus. Thus, G. aurea is an excellent model organism for studying evolutionary mechanisms of genome contraction.
 RESULTS Here we report sequencing and de novo draft assembly of G. aurea genome. The assembly consists of 10,687 contigs of the total length of 43.4 Mb and includes 17,755 complete and partial protein-coding genes. Its comparison with the genome of Mimulus guttatus, another representative of higher core Lamiales clade, reveals striking differences in gene content and length of non-coding regions.
 CONCLUSIONS Genome contraction was a complex process, which involved gene loss and reduction of lengths of introns and intergenic regions, but not intron loss. The gene loss is more frequent for the genes that belong to multigenic families indicating that genetic redundancy is an important prerequisite for genome size reduction.

text:
BACKGROUND RNA-Seq is revolutionizing the way transcript abundances are measured. A key challenge in transcript quantification from RNA-Seq data is the handling of reads that map to multiple genes or isoforms. This issue is particularly important for quantification with de novo transcriptome assemblies in the absence of sequenced genomes, as it is difficult to determine which transcripts are isoforms of the same gene. A second significant issue is the design of RNA-Seq experiments, in terms of the number of reads, read length, and whether reads come from one or both ends of cDNA fragments.
 RESULTS We present RSEM, an user-friendly software package for quantifying gene and isoform abundances from single-end or paired-end RNA-Seq data. RSEM outputs abundance estimates, 95% credibility intervals, and visualization files and can also simulate RNA-Seq data. In contrast to other existing tools, the software does not require a reference genome. Thus, in combination with a de novo transcriptome assembler, RSEM enables accurate transcript quantification for species without sequenced genomes. On simulated and real data sets, RSEM has superior or comparable performance to quantification methods that rely on a reference genome. Taking advantage of RSEM's ability to effectively use ambiguously-mapping reads, we show that accurate gene-level abundance estimates are best obtained with large numbers of short single-end reads. On the other hand, estimates of the relative frequencies of isoforms within single genes may be improved through the use of paired-end reads, depending on the number of possible splice forms for each gene.
 CONCLUSIONS RSEM is an accurate and user-friendly software tool for quantifying transcript abundances from RNA-Seq data. As it does not rely on the existence of a reference genome, it is particularly useful for quantification with de novo transcriptome assemblies. In addition, RSEM has enabled valuable guidance for cost-efficient design of quantification experiments with RNA-Seq, which is currently relatively expensive.

title: BLAST+: architecture and applications

title: Spatial genome organization: contrasting views from chromosome conformation capture and fluorescence in situ hybridization.

text: Although important for gene regulation, most studies of genome organization use either fluorescence in situ hybridization (FISH) or chromosome conformation capture (3C) methods. FISH directly visualizes the spatial relationship of sequences but is usually applied to a few loci at a time. The frequency at which sequences are ligated together by formaldehyde cross-linking can be measured genome-wide by 3C methods, with higher frequencies thought to reflect shorter distances. FISH and 3C should therefore give the same views of genome organization, but this has not been tested extensively. We investigated the murine HoxD locus with 3C carbon copy (5C) and FISH in different developmental and activity states and in the presence or absence of epigenetic regulators. We identified situations in which the two data sets are concordant but found other conditions under which chromatin topographies extrapolated from 5C or FISH data are not compatible. We suggest that products captured by 3C do not always reflect spatial proximity, with ligation occurring between sequences located hundreds of nanometers apart, influenced by nuclear environment and chromatin composition. We conclude that results obtained at high resolution with either 3C methods or FISH alone must be interpreted with caution and that views about genome organization should be validated by independent methods.

title: Nucleation of nuclear bodies by RNA

title: Cutadapt removes adapter sequences from high-throughput sequencing reads

text: When small RNA is sequenced on current sequencing machines, the resulting reads are usually longer than the RNA and therefore contain parts of the 3' adapter. That adapter must be found and removed error-tolerantly from each read before read mapping. Previous solutions are either hard to use or do not offer required features, in particular support for color space data. As an easy to use alternative, we developed the command-line tool cutadapt, which supports 454, Illumina and SOLiD (color space) data, offers two adapter trimming algorithms, and has other useful features. Cutadapt, including its MIT-licensed source code, is available for download at http://code.google.com/p/cutadapt/


query: Decrease of p62 in prostate tumor stroma results in defective autophagy.

title: p62/SQSTM1 by Binding to Vitamin D Receptor Inhibits Hepatic Stellate Cell Activity, Fibrosis, and Liver Cancer.

text: Hepatic stellate cells (HSCs) play critical roles in liver fibrosis and hepatocellular carcinoma (HCC). Vitamin D receptor (VDR) activation in HSCs inhibits liver inflammation and fibrosis. We found that p62/SQSTM1, a protein upregulated in liver parenchymal cells but downregulated in HCC-associated HSCs, negatively controls HSC activation. Total body or HSC-specific p62 ablation potentiates HSCs and enhances inflammation, fibrosis, and HCC progression. p62 directly interacts with VDR and RXR promoting their heterodimerization, which is critical for VDR:RXR target gene recruitment. Loss of p62 in HSCs impairs the repression of fibrosis and inflammation by VDR agonists. This demonstrates that p62 is a negative regulator of liver inflammation and fibrosis through its ability to promote VDR signaling in HSCs, whose activation supports HCC.

title: Metabolic reprogramming of stromal fibroblasts through p62-mTORC1 signaling promotes inflammation and tumorigenesis.

text: The tumor microenvironment plays a critical role in cancer progression, but the precise mechanisms by which stromal cells influence the epithelium are poorly understood. Here we show that p62 levels were reduced in the stroma of several tumors and that its loss in the tumor microenvironment or stromal fibroblasts resulted in increased tumorigenesis of epithelial prostate cancer cells. The mechanism involves the regulation of cellular redox through an mTORC1/c-Myc pathway of stromal glucose and amino acid metabolism, resulting in increased stromal IL-6 production, which is required for tumor promotion in the epithelial compartment. Thus, p62 is an anti-inflammatory tumor suppressor that acts through the modulation of metabolism in the tumor stroma.

title: Promotion of tumorigenesis by heterozygous disruption of the beclin 1 autophagy gene.

text: Malignant cells often display defects in autophagy, an evolutionarily conserved pathway for degrading long-lived proteins and cytoplasmic organelles. However, as yet, there is no genetic evidence for a role of autophagy genes in tumor suppression. The beclin 1 autophagy gene is monoallelically deleted in 40-75% of cases of human sporadic breast, ovarian, and prostate cancer. Therefore, we used a targeted mutant mouse model to test the hypothesis that monoallelic deletion of beclin 1 promotes tumorigenesis. Here we show that heterozygous disruption of beclin 1 increases the frequency of spontaneous malignancies and accelerates the development of hepatitis B virus-induced premalignant lesions. Molecular analyses of tumors in beclin 1 heterozygous mice show that the remaining wild-type allele is neither mutated nor silenced. Furthermore, beclin 1 heterozygous disruption results in increased cellular proliferation and reduced autophagy in vivo. These findings demonstrate that beclin 1 is a haplo-insufficient tumor-suppressor gene and provide genetic evidence that autophagy is a novel mechanism of cell-growth control and tumor suppression. Thus, mutation of beclin 1 or other autophagy genes may contribute to the pathogenesis of human cancers.

title: miR-205 hinders the malignant interplay between prostate cancer cells and associated fibroblasts.

text:
AIMS Tumor microenvironment is a strong determinant for the acquisition of metastatic potential of cancer cells. We have recently demonstrated that cancer-associated fibroblasts (CAFs) elicit a redox-dependent epithelial-mesenchymal transition (EMT) in prostate cancer (PCa) cells, driven by cycloxygenase-2/hypoxia-inducible factor-1 (HIF-1)/nuclear factor-κB pathway and enhancing tumor aggressiveness. Here, we investigated the involvement of microRNAs (miRNAs) in tumor-stroma interplay to identify possible tools to counteract oxidative stress and metastasis dissemination.
 RESULTS We found that miR-205 is the most downmodulated miRNA in PCa cells upon CAF stimulation, due to direct transcriptional repression by HIF-1, a known redox-sensitive transcription factor. Rescue experiments demonstrated that ectopic miR-205 overexpression in PCa cells counteracts CAF-induced EMT, thus impairing enhancement of cell invasion, acquisition of stem cell traits, tumorigenicity, and metastatic dissemination. In addition, miR-205 blocks tumor-driven activation of surrounding fibroblasts by reducing pro-inflammatory cytokine secretion.    INNOVATION Overall, such findings suggest miR-205 as a brake against PCa metastasis by blocking both the afferent and efferent arms of the circuit between tumor cells and associated fibroblasts, thus interrupting the pro-oxidant and pro-inflammatory circuitries engaged by reactive stroma.
 CONCLUSION The evidence that miR-205 replacement in PCa cells is able not only to prevent but also to revert the oxidative/pro-inflammatory axis leading to EMT induced by CAFs sets the rationale for developing miRNA-based approaches to prevent and treat metastatic disease.

title: Oxidative stress in cancer associated fibroblasts drives tumor-stroma co-evolution: A new paradigm for understanding tumor metabolism, the field effect and genomic instability in cancer cells.

text: Loss of stromal fibroblast caveolin-1 (Cav-1) is a powerful single independent predictor of poor prognosis in human breast cancer patients, and is associated with early tumor recurrence, lymph node metastasis and tamoxifen-resistance. We developed a novel co-culture system to understand the mechanism(s) by which a loss of stromal fibroblast Cav-1 induces a "lethal tumor micro-environment. " Here, we propose a new paradigm to explain the powerful prognostic value of stromal Cav-1. In this model, cancer cells induce oxidative stress in cancer-associated fibroblasts, which then acts as a "metabolic" and "mutagenic" motor to drive tumor-stroma co-evolution, DNA damage and aneuploidy in cancer cells. More specifically, we show that an acute loss of Cav-1 expression leads to mitochondrial dysfunction, oxidative stress and aerobic glycolysis in cancer associated fibroblasts. Also, we propose that defective mitochondria are removed from cancer-associated fibroblasts by autophagy/mitophagy that is induced by oxidative stress. As a consequence, cancer associated fibroblasts provide nutrients (such as lactate) to stimulate mitochondrial biogenesis and oxidative metabolism in adjacent cancer cells (the "Reverse Warburg Effect"). We provide evidence that oxidative stress in cancer-associated fibroblasts is sufficient to induce genomic instability in adjacent cancer cells, via a bystander effect, potentially increasing their aggressive behavior. Finally, we directly demonstrate that nitric oxide (NO) over-production, secondary to Cav-1 loss, is the root cause for mitochondrial dysfunction in cancer associated fibroblasts. In support of this notion, treatment with anti-oxidants (such as N-acetyl-cysteine, metformin and quercetin) or NO inhibitors (L-NAME) was sufficient to reverse many of the cancer-associated fibroblast phenotypes that we describe. Thus, cancer cells use "oxidative stress" in adjacent fibroblasts (i) as an "engine" to fuel their own survival via the stromal production of nutrients and (ii) to drive their own mutagenic evolution towards a more aggressive phenotype, by promoting genomic instability. We also present evidence that the "field effect" in cancer biology could also be related to the stromal production of ROS and NO species. eNOS-expressing fibroblasts have the ability to downregulate Cav-1 and induce mitochondrial dysfunction in adjacent fibroblasts that do not express eNOS. As such, the effects of stromal oxidative stress can be laterally propagated, amplified and are effectively "contagious"--spread from cell-to-cell like a virus--creating an "oncogenic/mutagenic" field promoting widespread DNA damage.

title: Defective regulation of adipose tissue autophagy in obesity

text: Objectives: Autophagy is a highly regulated process that has an important role in the control of a wide range of cellular functions, such as organelle recycling, nutrient availability and tissue differentiation. A recent study has shown an increased autophagic activity in the adipose tissue of obese subjects, and a role for autophagy in obesity-associated insulin resistance was proposed. Body mass reduction is the most efficient approach to tackle insulin resistance in over-weight subjects; however, the impact of weight loss in adipose tissue autophagy is unknown. Subjects:Adipose tissue autophagy was evaluated in mice and humans. Results:First, a mouse model of diet-induced obesity and diabetes was maintained on a 15-day, 40% caloric restriction. At baseline, markers of autophagy were increased in obese mice as compared with lean controls. Upon caloric restriction, autophagy increased in the lean mice, whereas it decreased in the obese mice. The reintroduction of ad libitum feeding was sufficient to rapidly reduce autophagy in the lean mice and increase autophagy in the obese mice. In the second part of the study, autophagy was evaluated in the subcutaneous adipose tissue of nine obese-non-diabetic and six obese-diabetic subjects undergoing bariatric surgery for body mass reduction. Specimens were collected during the surgery and approximately 1 year later. Markers of systemic inflammation, such as tumor necrosis factor-1α, interleukin (IL)-6 and IL-1β were evaluated. As in the mouse model, human obesity was associated with increased autophagy, and body mass reduction led to an attenuation of autophagy in the adipose tissue. Conclusion:Obesity and caloric overfeeding are associated with the defective regulation of autophagy in the adipose tissue. The studies in obese-diabetic subjects undergoing improved metabolic control following calorie restriction suggest that autophagy and inflammation are regulated independently.

title: Autophagy and mTORC1 regulate the stochastic phase of somatic cell reprogramming

text: We describe robust induction of autophagy during the reprogramming of mouse fibroblasts to induced pluripotent stem cells by four reprogramming factors (Sox2, Oct4, Klf4 and c-Myc), henceforth 4F. This process occurs independently of p53 activation, and is mediated by the synergistic downregulation of mechanistic target of rapamycin complex 1 (mTORC1) and the induction of autophagy-related genes. The 4F coordinately repress mTORC1, but bifurcate in their regulation of autophagy-related genes, with Klf4 and c-Myc inducing them but Sox2 and Oct4 inhibiting them. On one hand, inhibition of mTORC1 facilitates reprogramming by promoting cell reshaping (mitochondrial remodelling and cell size reduction). On the other hand, mTORC1 paradoxically impairs reprogramming by triggering autophagy. Autophagy does not participate in cell reshaping in reprogramming but instead degrades p62, whose accumulation in autophagy-deficient cells facilitates reprogramming. Our results thus reveal a complex signalling network involving mTORC1 inhibition and autophagy induction in the early phase of reprogramming, whose delicate balance ultimately determines reprogramming efficiency.

title: Stromal-epithelial metabolic coupling in cancer: integrating autophagy and metabolism in the tumor microenvironment.

text: Cancer cells do not exist as pure homogeneous populations in vivo. Instead they are embedded in "cancer cell nests" that are surrounded by stromal cells, especially cancer associated fibroblasts. Thus, it is not unreasonable to suspect that stromal fibroblasts could influence the metabolism of adjacent cancer cells, and visa versa. In accordance with this idea, we have recently proposed that the Warburg effect in cancer cells may be due to culturing cancer cells by themselves, out of their normal stromal context or tumor microenvironment. In fact, when cancer cells are co-cultured with fibroblasts, then cancer cells increase their mitochondrial mass, while fibroblasts lose their mitochondria. An in depth analysis of this phenomenon reveals that aggressive cancer cells are "parasites" that use oxidative stress as a "weapon" to extract nutrients from surrounding stromal cells. Oxidative stress in fibroblasts induces the autophagic destruction of mitochondria, by mitophagy. Then, stromal cells are forced to undergo aerobic glycolysis, and produce energy-rich nutrients (such as lactate and ketones) to "feed" cancer cells. This mechanism would allow cancer cells to seed anywhere, without blood vessels as a food source, as they could simply induce oxidative stress wherever they go, explaining how cancer cells survive during metastasis. We suggest that stromal catabolism, via autophagy and mitophagy, fuels the anabolic growth of tumor cells, promoting tumor progression and metastasis. We have previously termed this new paradigm "The Autophagic Tumor Stroma Model of Cancer Metabolism", or the "Reverse Warburg Effect". We also discuss how glutamine addiction (glutaminolysis) in cancer cells fits well with this new model, by promoting oxidative mitochondrial metabolism in aggressive cancer cells.

title: Epigenetic-induced repression of microRNA-205 is associated with MED1 activation and a poorer prognosis in localized prostate cancer

text: Deregulation of microRNA (miRNA) expression can have a critical role in carcinogenesis. Here we show in prostate cancer that miRNA-205 (miR-205) transcription is commonly repressed and the MIR-205 locus is hypermethylated. LOC642587, the MIR-205 host gene of unknown function, is also concordantly inactivated. We show that miR-205 targets mediator 1 (MED1, also called TRAP220 and PPARBP) for transcriptional silencing in normal prostate cells, leading to reduction in MED1 mRNA levels, and in total and active phospho-MED1 protein. Overexpression of miR-205 in prostate cancer cells negatively affects cell viability, consistent with a tumor suppressor function. We found that hypermethylation of the MIR-205 locus was strongly related with a decrease in miR-205 expression and an increase in MED1 expression in primary tumor samples (n=14), when compared with matched normal prostate (n=7). An expanded patient cohort (tumor n=149, matched normal n=30) also showed significant MIR-205 DNA methylation in tumors compared with normal, and MIR-205 hypermethylation is significantly associated with biochemical recurrence (hazard ratio=2.005, 95% confidence interval (1.109, 3.625), P=0.02), in patients with low preoperative prostate specific antigen. In summary, these results suggest that miR-205 is an epigenetically regulated tumor suppressor that targets MED1 and may provide a potential biomarker in prostate cancer management.

title: Elimination of paternal mitochondria in mouse embryos occurs through autophagic degradation dependent on PARKIN and MUL1

text: A defining feature of mitochondria is their maternal mode of inheritance. However, little is understood about the cellular mechanism through which paternal mitochondria, delivered from sperm, are eliminated from early mammalian embryos. Autophagy has been implicated in nematodes, but whether this mechanism is conserved in mammals has been disputed. Here, we show that cultured mouse fibroblasts and pre-implantation embryos use a common pathway for elimination of mitochondria. Both situations utilize mitophagy, in which mitochondria are sequestered by autophagosomes and delivered to lysosomes for degradation. The E3 ubiquitin ligases PARKIN and MUL1 play redundant roles in elimination of paternal mitochondria. The process is associated with depolarization of paternal mitochondria and additionally requires the mitochondrial outer membrane protein FIS1, the autophagy adaptor P62, and PINK1 kinase. Our results indicate that strict maternal transmission of mitochondria relies on mitophagy and uncover a collaboration between MUL1 and PARKIN in this process.

title: The calcium sensor STIM1 is regulated by androgens in prostate stromal cells.

text:
BACKGROUND Prostate development and maintenance in the adult results from an interaction of stromal and glandular components. Androgens can drive this process by direct action on the stroma. We investigated whether there was a direct link between androgens and another key regulator of stromal cells, intracellular Ca2+ ([Ca2+ ]i ).
 METHODS Prostate stromal cells were freshly obtained and cultures derived from patients with benign prostatic hyperplasia. Gene expression in dihydrotestosterone treated and untreated cells was compared using Affymetrix gene expression arrays and Ca2+ regulated features were identified by Gene Ontology (GO). Changes in [Ca2+]i were determined in Fluo-4 loaded cells. Androgen regulation was confirmed by chromatin immunoprecipitaion.
 RESULTS Stromal cell cultures were sorted for expression of integrin α1 β1 , which enriched for cells expressing the androgen receptor (AR). We identified key functional categories, within the androgen-induced gene expression signature, focusing on genes involved in calcium signaling. From this analysis, stromal interaction molecule-1 (STIM1) was identified as a significantly differentially expressed gene with four relevant associated GO terms. DNA sequence analysis showed that the promoter region of STIM1 contained putative androgen response element sequences in which AR binding ability of STIM1 was confirmed. Androgens directly regulated STIM1 expression and STIM1 effects on store-operated calcium entry were inhibited by STIM1 knock-down. Reduced STIM1 expression in prostate stromal cells led to a reduction in basal Ca2+ levels, the amount of Ca2+ released by thapsigargin and a reduction in store filling following TG-induced store depletion.
 CONCLUSIONS These results indicate that androgens modulate [Ca2+]i through the direct regulation of the STIM1 gene by AR binding to the STIM1 promoter.

title: Rapamycin treatment augments motor neuron degeneration in SOD1(G93A) mouse model of amyotrophic lateral sclerosis.

text: Aberrant protein misfolding may contribute to the pathogenesis of amyotrophic lateral sclerosis (ALS) but the detailed mechanisms are largely unknown. Our previous study has shown that autophagy is altered in the mouse model of ALS. In the present study, we systematically investigated the correlation of the autophagic alteration with the motor neurons (MNs) degeneration in the ALS mice. We have demonstrated that the autophagic protein marker LC3-II is markedly and specifically increased in the spinal cord MNs of the ALS mice. Electron microscopy and immunochemistry studies have shown that autophagic vacuoles are significantly accumulated in the dystrophic axons of spinal cord MNs of the ALS mice. All these changes in the ALS mice appear at the age of 90 d when the ALS mice display modest clinical symptoms; and they become prominent at the age of 120 d. The clinical symptoms are correlated with the progression of MNs degeneration. Moreover, we have found that p62/SQSTM1 is accumulated progressively in the spinal cord, indicating that the possibility of impaired autophagic flux in the SOD1(G93A) mice. Furthermore, to our surprise, we have found that treatment with autophagy enhancer rapamycin accelerates the MNs degeneration, shortens the life span of the ALS mice, and has no obvious effects on the accumulation of SOD1 aggregates. In addition, we have demonstrated that rapamycin treatment in the ALS mice causes more severe mitochondrial impairment, higher Bax levels and greater caspase-3 activation. These findings suggest that selective degeneration of MNs is associated with the impairment of the autophagy pathway and that rapamycin treatment may exacerbate the pathological processing through apoptosis and other mechanisms in the ALS mice.

title: Regulation of E2Fs and senescence by PML nuclear bodies.

text: The tumor suppressor PML (promyelocytic leukemia protein) regulates cellular senescence and terminal differentiation, two processes that implicate a permanent exit from the cell cycle. Here, we show that the mechanism by which PML induces a permanent cell cycle exit and activates p53 and senescence involves a recruitment of E2F transcription factors bound to their promoters and the retinoblastoma (Rb) proteins to PML nuclear bodies enriched in heterochromatin proteins and protein phosphatase 1α. Blocking the functions of the Rb protein family or adding back E2Fs to PML-expressing cells can rescue their defects in E2F-dependent gene expression and cell proliferation, inhibiting the senescent phenotype. In benign prostatic hyperplasia, a neoplastic disease that displays features of senescence, PML was found to be up-regulated and forming nuclear bodies. In contrast, PML bodies were rarely visualized in prostate cancers. The newly defined PML/Rb/E2F pathway may help to distinguish benign tumors from cancers, and suggest E2F target genes as potential targets to induce senescence in human tumors.

title: Antagonistic Controls of Autophagy and Glycogen Accumulation by Snf1p, the Yeast Homolog of AMP-Activated Protein Kinase, and the Cyclin-Dependent

text: In the yeast Saccharomyces cerevisiae, glycogen is accumulated as a carbohydrate reserve when cells are deprived of nutrients. Yeast mutated in SNF1, a gene encoding a protein kinase required for glucose derepression, has diminished glycogen accumulation and concomitant inactivation of glycogen synthase. Restoration of synthesis in an snf1 strain results only in transient glycogen accumulation, implying the existence of other SNF1-dependent controls of glycogen storage. A genetic screen revealed that two genes involved in autophagy, APG1 and APG13, may be regulated by SNF1. Increased autophagic activity was observed in wild-type cells entering the stationary phase, but this induction was impaired in an snf1 strain. Mutants defective for autophagy were able to synthesize glycogen upon approaching the stationary phase, but were unable to maintain their glycogen stores, because subsequent synthesis was impaired and degradation by phosphorylase, Gph1p, was enhanced. Thus, deletion of GPH1 partially reversed the loss of glycogen accumulation in autophagy mutants. Loss of the vacuolar glucosidase, SGA1, also protected glycogen stores, but only very late in the stationary phase. Gph1p and Sga1p may therefore degrade physically distinct pools of glycogen. Pho85p is a cyclin-dependent protein kinase that antagonizes SNF1 control of glycogen synthesis. Induction of autophagy in pho85 mutants entering the stationary phase was exaggerated compared to the level in wild-type cells, but was blocked in apg1 pho85 mutants. We propose that Snf1p and Pho85p are, respectively, positive and negative regulators of autophagy, probably via Apg1 and/or Apg13. Defective glycogen storage in snf1 cells can be attributed to both defective synthesis upon entry into stationary phase and impaired maintenance of glycogen levels caused by the lack of autophagy.

title: Restoration of chaperone-mediated autophagy in aging liver improves cellular maintenance and hepatic function

text: Chaperone-mediated autophagy (CMA), a selective mechanism for degradation of cytosolic proteins in lysosomes, contributes to the removal of altered proteins as part of the cellular quality-control systems. We have previously found that CMA activity declines in aged organisms and have proposed that this failure in cellular clearance could contribute to the accumulation of altered proteins, the abnormal cellular homeostasis and, eventually, the functional loss characteristic of aged organisms. To determine whether these negative features of aging can be prevented by maintaining efficient autophagic activity until late in life, in this work we have corrected the CMA defect in aged rodents. We have generated a double transgenic mouse model in which the amount of the lysosomal receptor for CMA, previously shown to decrease in abundance with age, can be modulated. We have analyzed in this model the consequences of preventing the age-dependent decrease in receptor abundance in aged rodents at the cellular and organ levels. We show here that CMA activity is maintained until advanced ages if the decrease in the receptor abundance is prevented and that preservation of autophagic activity is associated with lower intracellular accumulation of damaged proteins, better ability to handle protein damage and improved organ function.

title: 20-year outcomes following conservative management of clinically localized prostate cancer.

text:
CONTEXT The appropriate therapy for men with clinically localized prostate cancer is uncertain. A recent study suggested an increasing prostate cancer mortality rate for men who are alive more than 15 years following diagnosis.
 OBJECTIVE To estimate 20-year survival based on a competing risk analysis of men who were diagnosed with clinically localized prostate cancer and treated with observation or androgen withdrawal therapy alone, stratified by age at diagnosis and histological findings.
 DESIGN, SETTING, AND PATIENTS A retrospective population-based cohort study using Connecticut Tumor Registry data supplemented by hospital record and histology review of 767 men aged 55 to 74 years with clinically localized prostate cancer diagnosed between January 1, 1971, and December 31, 1984. Patients were treated with either observation or immediate or delayed androgen withdrawal therapy, with a median observation of 24 years.
 MAIN OUTCOME MEASURES Probability of mortality from prostate cancer or other competing medical conditions, given a patient's age at diagnosis and tumor grade.
 RESULTS The prostate cancer mortality rate was 33 per 1000 person-years during the first 15 years of follow-up (95% confidence interval [CI], 28-38) and 18 per 1000 person-years after 15 years of follow-up (95% CI, 10-29). The mortality rates for these 2 follow-up periods were not statistically different, after adjusting for differences in tumor histology (rate ratio, 1.1; 95% CI, 0.6-1.9). Men with low-grade prostate cancers have a minimal risk of dying from prostate cancer during 20 years of follow-up (Gleason score of 2-4, 6 deaths per 1000 person-years; 95% CI, 2-11). Men with high-grade prostate cancers have a high probability of dying from prostate cancer within 10 years of diagnosis (Gleason score of 8-10, 121 deaths per 1000 person-years; 95% CI, 90-156). Men with Gleason score of 5 or 6 tumors have an intermediate risk of prostate cancer death.
 CONCLUSION The annual mortality rate from prostate cancer appears to remain stable after 15 years from diagnosis, which does not support aggressive treatment for localized low-grade prostate cancer.

title: Coactivator SRC-2-dependent metabolic reprogramming mediates prostate cancer survival and metastasis.

text: Metabolic pathway reprogramming is a hallmark of cancer cell growth and survival and supports the anabolic and energetic demands of these rapidly dividing cells. The underlying regulators of the tumor metabolic program are not completely understood; however, these factors have potential as cancer therapy targets. Here, we determined that upregulation of the oncogenic transcriptional coregulator steroid receptor coactivator 2 (SRC-2), also known as NCOA2, drives glutamine-dependent de novo lipogenesis, which supports tumor cell survival and eventual metastasis. SRC-2 was highly elevated in a variety of tumors, especially in prostate cancer, in which SRC-2 was amplified and overexpressed in 37% of the metastatic tumors evaluated. In prostate cancer cells, SRC-2 stimulated reductive carboxylation of α-ketoglutarate to generate citrate via retrograde TCA cycling, promoting lipogenesis and reprogramming of glutamine metabolism. Glutamine-mediated nutrient signaling activated SRC-2 via mTORC1-dependent phosphorylation, which then triggered downstream transcriptional responses by coactivating SREBP-1, which subsequently enhanced lipogenic enzyme expression. Metabolic profiling of human prostate tumors identified a massive increase in the SRC-2-driven metabolic signature in metastatic tumors compared with that seen in localized tumors, further implicating SRC-2 as a prominent metabolic coordinator of cancer metastasis. Moreover, SRC-2 inhibition in murine models severely attenuated the survival, growth, and metastasis of prostate cancer. Together, these results suggest that the SRC-2 pathway has potential as a therapeutic target for prostate cancer.

title: miR-375 inhibits autophagy and reduces viability of hepatocellular carcinoma cells under hypoxic conditions.

text:
BACKGROUND & AIMS Tumor cells survive hypoxic conditions by inducing autophagy. We investigated the roles of microRNAs (miRNAs) in regulating autophagy of hepatocellular carcinoma (HCC) cells under hypoxic conditions.
 METHODS We used gain- and loss-of-function methods to evaluate the effect of miRNAs on autophagy in human HCC cell lines (Huh7 and Hep3B) under hypoxic conditions. Autophagy was quantified by immunoblot, immunofluoresence, and transmission electron microscopy analyses, and after incubation of cells with bafilomycin A1. We used a luciferase reporter assay to confirm associations between miRNAs and their targets. We analyzed growth of HCC xenograft tumors in nude mice.
 RESULTS miR-375 was down-regulated in HCC cells and tissues; it inhibited autophagy under hypoxic conditions by suppressing the conversion of LC3I to LC3II and thereby autophagic flux. The ability of miR-375 to inhibit autophagy was independent of its ability to regulate 3'-phosphoinositide-dependent protein kinase-1-AKT-mammalian target of rapamycin signaling, but instead involved suppression of ATG7, an autophagy-associated gene. miR-375 bound directly to a predicted site in the 3' untranslated region of ATG7. Up-regulating miR-375 or down-regulating ATG7 inhibited mitochondrial autophagy of HCC cells, reduced the elimination of damaged mitochondria under hypoxia, increased release of mitochondrial apoptotic proteins, and reduced viability of HCC cells. In mice, xenograft tumors that expressed miR-375 had fewer autophagic cells, larger areas of necrosis, and grew more slowly than tumors from HCC cells that expressed lower levels of miR-375.
 CONCLUSIONS miR-375 inhibits autophagy by reducing expression of ATG7 and impairs viability of HCC cells under hypoxic conditions in culture and in mice. miRNAs that inhibit autophagy of cancer cells might be developed as therapeutics.

title: Prostate cancer-associated gene expression alterations determined from needle biopsies.

text:
PURPOSE To accurately identify gene expression alterations that differentiate neoplastic from normal prostate epithelium using an approach that avoids contamination by unwanted cellular components and is not compromised by acute gene expression changes associated with tumor devascularization and resulting ischemia.    EXPERIMENTAL DESIGN Approximately 3,000 neoplastic and benign prostate epithelial cells were isolated using laser capture microdissection from snap-frozen prostate biopsy specimens provided by 31 patients who subsequently participated in a clinical trial of preoperative chemotherapy. cDNA synthesized from amplified total RNA was hybridized to custom-made microarrays composed of 6,200 clones derived from the Prostate Expression Database. Expression differences for selected genes were verified using quantitative reverse transcription-PCR.
 RESULTS Comparative analyses identified 954 transcript alterations associated with cancer (q < 0.01%), including 149 differentially expressed genes with no known functional roles. Gene expression changes associated with ischemia and surgical removal of the prostate gland were absent. Genes up-regulated in prostate cancer were statistically enriched in categories related to cellular metabolism, energy use, signal transduction, and molecular transport. Genes down-regulated in prostate cancers were enriched in categories related to immune response, cellular responses to pathogens, and apoptosis. A heterogeneous pattern of androgen receptor expression changes was noted. In exploratory analyses, androgen receptor down-regulation was associated with a lower probability of cancer relapse after neoadjuvant chemotherapy followed by radical prostatectomy.
 CONCLUSIONS Assessments of tumor phenotypes based on gene expression for treatment stratification and drug targeting of oncogenic alterations may best be ascertained using biopsy-based analyses where the effects of ischemia do not complicate interpretation.

title: Changing the course of pancreatic cancer--Focus on recent translational advances.

text: In the past decade, insightful preclinical research has led to important breakthroughs in our understanding of pancreatic cancer. Even though the vast majority of pancreatic cancers are KRAS mutated, not all pancreatic cancer tumors are "KRAS equal"; there seems to be varying dependencies on the KRAS pathway. While KRAS-targeting therapies have been disappointing in the clinic, 'synthetic lethal' approaches hold promise in this setting. The pancreatic cancer stromal microenvironment appears to have contradictory roles. While there is evidence to suggest that stromal barrier prevents drug delivery, in other circumstances, stroma can play a protective role and its disruption enhances tumor dissemination. Clinical trials aimed at manipulating the various stromal components are in progress. BRCA mutation-related pancreatic tumors illustrate a unique subtype with enhanced susceptibility to DNA damaging agents and PARP-inhibition. DNA repair defects in cancer extend beyond germ line BRCA mutation and may extend the indications for DNA repair-targeting agents. Immune strategies are an area of active investigation in pancreatic cancer. Although the initial trials of single-agent checkpoint inhibitors have been negative, combinational approaches using immune-modifying agents and vaccines appear promising and goal is to identify an 'immune-therapy responsive' profile in pancreatic cancer.

title: Tumor engraftment in nude mice and enrichment in stroma- related gene pathways predict poor survival and resistance to gemcitabine in patients with pancreatic cancer.

text:
PURPOSE The goal of this study was to evaluate prospectively the engraftment rate, factors influencing engraftment, and predictability of clinical outcome of low-passage xenografts from patients with resectable pancreatic ductal adenocarcinoma (PDA) and to establish a bank of PDA xenografts.    EXPERIMENTAL DESIGN Patients with resectable PDA scheduled for resection at the Johns Hopkins Hospital were eligible. Representative pieces of tumor were implanted in nude mice. The status of the SMAD4 gene and content of tumor-generating cells were determined by immunohistochemistry. Gene expression was carried out by using a U133 Plus 2.0 array. Patients were followed for progression and survival.
 RESULTS A total of 94 patients with PDA were resected, 69 tumors implanted in nude mice, and 42 (61%) engrafted. Engrafted carcinomas were more often SMAD4 mutant, and had a metastatic gene expression signature and worse prognosis. Tumors from patients resistant to gemcitabine were enriched in stroma-related gene pathways. Tumors sensitive to gemcitabine were enriched in cell cycle and pyrimidine gene pathways. The time to progression for patients who received treatment with gemcitabine for metastatic disease (n = 7) was double in patients with xenografts sensitive to gemcitabine.
 CONCLUSION A successful xenograft was generated in 61% of patients attempted, generating a pool of 42 PDA xenografts with significant biological information and annotated clinical data. Patients with PDA and SMAD4 inactivation have a better engraftment rate. Engraftment is a poor prognosis factor, and engrafted tumors have a metastatic gene expression signature. Tumors from gemcitabine-resistant patients were enriched in stromal pathways.