An estimated 15% or more of the cancer burden worldwide is attributable to known infectious agents. We screened colorectal carcinoma and matched normal tissue specimens using RNA-seq followed by host sequence subtraction and found marked over-representation of Fusobacterium nucleatum sequences in tumors relative to control specimens. F. nucleatum is an invasive anaerobe that has been linked previously to periodontitis and appendicitis, but not to cancer. Fusobacteria are rare constituents of the fecal microbiota, but have been cultured previously from biopsies of inflamed gut mucosa. We obtained a Fusobacterium isolate from a frozen tumor specimen; this showed highest sequence similarity to a known gut mucosa isolate and was confirmed to be invasive. We verified overabundance of Fusobacterium sequences in tumor versus matched normal control tissue by quantitative PCR analysis from a total of 99 subjects (p = 2.5 × 10⁻⁶), and we observed a positive association with lymph node metastasis.

Next-generation sequencing approaches have been used to investigate the genomes and transcriptomes of an oestrogen-receptor-α-positive metastatic lobular breast cancer from a patient — rather than from a cell line or xenograft — over a 9-year period between the diagnosis of the primary tumour and the appearance of metastasis. Comparison of the somatic non-synonymous coding mutations in the metastasis and the primary tumour of the same patient and the combined analysis of genome and transcriptome data provided insights into the mutational evolution that can occur with disease progression. The cover shows sequence elements of the HAUS3 locus, one of the genes found to be mutated in the tissue (shown in the background) from the primary lobular cancer used for this work. Advances in next generation sequencing have made it possible to precisely characterize the coding mutations that occur during the development and progression of individual cancers. Here, this technique is used to sequence the genomes and transcriptomes of an oestrogen-receptor-α-positive metastatic lobular breast cancer; significant evolution is found to occur with disease progression. Recent advances in next generation sequencing1,2,3,4 have made it possible to precisely characterize all somatic coding mutations that occur during the development and progression of individual cancers. Here we used these approaches to sequence the genomes (>43-fold coverage) and transcriptomes of an oestrogen-receptor-α-positive metastatic lobular breast cancer at depth. We found 32 somatic non-synonymous coding mutations present in the metastasis, and measured the frequency of these somatic mutations in DNA from the primary tumour of the same patient, which arose 9 years earlier. Five of the 32 mutations (in ABCB11, HAUS3, SLC24A4, SNX4 and PALB2) were prevalent in the DNA of the primary tumour removed at diagnosis 9 years earlier, six (in KIF1C, USP28, MYH8, MORC1, KIAA1468 and RNASEH2A) were present at lower frequencies (1–13%), 19 were not detected in the primary tumour, and two were undetermined. The combined analysis of genome and transcriptome data revealed two new RNA-editing events that recode the amino acid sequence of SRP9 and COG3. Taken together, our data show that single nucleotide mutational heterogeneity can be a property of low or intermediate grade primary breast cancers and that significant evolution can occur with disease progression.

The assembly of DNA sequences de novo is fundamental to genomics research. It is the first of many steps toward elucidating and characterizing whole genomes. Downstream applications, including analysis of genomic variation between species, between or within individuals critically depend on robustly assembled sequences. In the span of a single decade, the sequence throughput of leading DNA sequencing instruments has increased drastically, and coupled with established and planned large-scale, personalized medicine initiatives to sequence genomes in the thousands and even millions, the development of efficient, scalable and accurate bioinformatics tools for producing high-quality reference draft genomes is timely. With ABySS 1.0, we originally showed that assembling the human genome using short 50-bp sequencing reads was possible by aggregating the half terabyte of compute memory needed over several computers using a standardized message-passing system (MPI). We present here its redesign, which departs from MPI and instead implements algorithms that employ a Bloom filter, a probabilistic data structure, to represent a de Bruijn graph and reduce memory requirements. We benchmarked ABySS 2.0 human genome assembly using a Genome in a Bottle data set of 250-bp Illumina paired-end and 6-kbp mate-pair libraries from a single individual. Our assembly yielded a NG50 (NGA50) scaffold contiguity of 3.5 (3.0) Mbp using <35 GB of RAM. This is a modest memory requirement by today's standards and is often available on a single computer. We also investigate the use of BioNano Genomics and 10x Genomics’ Chromium data to further improve the scaffold NG50 (NGA50) of this assembly to 42 (15) Mbp.

The tissue‐restricted GATA‐4 transcription factor and Nkx2‐5 homeodomain protein are two early markers of precardiac cells. Both are essential for heart formation, but neither can initiate cardiogenesis. Overexpression of GATA‐4 or Nkx2‐5 enhances cardiac development in committed precursors, suggesting each interacts with a cardiac cofactor. We tested whether GATA‐4 and Nkx2‐5 are cofactors for each other by using transcription and binding assays with the cardiac atrial natriuretic factor (ANF) promoter—the only known target for Nkx2‐5. Co‐expression of GATA‐4 and Nkx2‐5 resulted in synergistic activation of the ANF promoter in heterologous cells. The synergy involves physical Nkx2‐5–GATA‐4 interaction, seen in vitro and in vivo , which maps to the C‐terminal zinc finger of GATA‐4 and a C‐terminus extension; similarly, a C‐terminally extended homeodomain of Nkx2‐5 is required for GATA‐4 binding. Structure/function studies suggest that binding of GATA‐4 to the C‐terminus autorepressive domain of Nkx2‐5 may induce a conformational change that unmasks Nkx2‐5 activation domains. GATA‐6 cannot substitute for GATA‐4 for interaction with Nkx2‐5. This interaction may impart functional specificity to GATA factors and provide cooperative crosstalk between two pathways critical for early cardiogenesis. Given the co‐expression of GATA proteins and NK2 class members in other tissues, the GATA/Nkx partnership may represent a paradigm for transcription factor interaction during organogenesis.

Rhodococcus sp. RHA1 (RHA1) is a potent polychlorinated biphenyl-degrading soil actinomycete that catabolizes a wide range of compounds and represents a genus of considerable industrial interest. RHA1 has one of the largest bacterial genomes sequenced to date, comprising 9,702,737 bp (67% G+C) arranged in a linear chromosome and three linear plasmids. A targeted insertion methodology was developed to determine the telomeric sequences. RHA1's 9,145 predicted protein-encoding genes are exceptionally rich in oxygenases (203) and ligases (192). Many of the oxygenases occur in the numerous pathways predicted to degrade aromatic compounds ( 30 ) or steroids ( 4 ). RHA1 also contains 24 nonribosomal peptide synthase genes, six of which exceed 25 kbp, and seven polyketide synthase genes, providing evidence that rhodococci harbor an extensive secondary metabolism. Among sequenced genomes, RHA1 is most similar to those of nocardial and mycobacterial strains. The genome contains few recent gene duplications. Moreover, three different analyses indicate that RHA1 has acquired fewer genes by recent horizontal transfer than most bacteria characterized to date and far fewer than Burkholderia xenovorans LB400, whose genome size and catabolic versatility rival those of RHA1. RHA1 and LB400 thus appear to demonstrate that ecologically similar bacteria can evolve large genomes by different means. Overall, RHA1 appears to have evolved to simultaneously catabolize a diverse range of plant-derived compounds in an O ₂ -rich environment. In addition to establishing RHA1 as an important model for studying actinomycete physiology, this study provides critical insights that facilitate the exploitation of these industrially important microorganisms.

Somatic missense mutations can initiate tumorogenesis and, conversely, anti-tumor cytotoxic T cell (CTL) responses. Tumor genome analysis has revealed extreme heterogeneity among tumor missense mutation profiles, but their relevance to tumor immunology and patient outcomes has awaited comprehensive evaluation. Here, for 515 patients from six tumor sites, we used RNA-seq data from The Cancer Genome Atlas to identify mutations that are predicted to be immunogenic in that they yielded mutational epitopes presented by the MHC proteins encoded by each patient’s autologousHLA-Aalleles. Mutational epitopes were associated with increased patient survival. Moreover, the corresponding tumors had higher CTL content, inferred fromCD8Agene expression, and elevated expression of the CTL exhaustion markersPDCD1andCTLA4. Mutational epitopes were very scarce in tumors without evidence of CTL infiltration. These findings suggest that the abundance of predicted immunogenic mutations may be useful for identifying patients likely to benefit from checkpoint blockade and related immunotherapies.