Comparison of tumor-informed and tumor-naïve sequencing assays for ctDNA detection in breast cancer.

Analysis of circulating tumor DNA (ctDNA) to monitor cancer dynamics and detect minimal residual disease has been an area of increasing interest. Multiple methods have been proposed but few studies have compared the performance of different approaches. Here, we compare detection of ctDNA in serial plasma samples from patients with breast cancer using different tumor-informed and tumor-naïve assays designed to detect structural variants (SVs), single nucleotide variants (SNVs), and/or somatic copy-number aberrations, by multiplex PCR, hybrid capture, and different depths of whole-genome sequencing. Our results demonstrate that the ctDNA dynamics and allele fractions (AFs) were highly concordant when analyzing the same patient samples using different assays. Tumor-informed assays showed the highest sensitivity for detection of ctDNA at low concentrations. Hybrid capture sequencing targeting between 1,347 and 7,491 tumor-identified mutations at high depth was the most sensitive assay, detecting ctDNA down to an AF of 0.00024% (2.4 parts per million, ppm). Multiplex PCR targeting 21-47 tumor-identified SVs per patient detected ctDNA down to 0.00047% AF (4.7 ppm) and has potential as a clinical assay.

NRG1 fusions in breast cancer.

BACKGROUND: NRG1 gene fusions may be clinically actionable, since cancers carrying the fusion transcripts can be sensitive to tyrosine kinase inhibitors. The NRG1 gene encodes ligands for the HER2(ERBB2)-ERBB3 heterodimeric receptor tyrosine kinase, and the gene fusions are thought to lead to autocrine stimulation of the receptor. The NRG1 fusion expressed in the breast cancer cell line MDA-MB-175 serves as a model example of such fusions, showing the proposed autocrine loop and exceptional drug sensitivity. However, its structure has not been properly characterised, its oncogenic activity has not been fully explained, and there is limited data on such fusions in breast cancer. METHODS: We analysed genomic rearrangements and transcripts of NRG1 in MDA-MB-175 and a panel of 571 breast cancers. RESULTS: We found that the MDA-MB-175 fusion-originally reported as a DOC4(TENM4)-NRG1 fusion, lacking the cytoplasmic tail of NRG1-is in reality a double fusion, PPP6R3-TENM4-NRG1, producing multiple transcripts, some of which include the cytoplasmic tail. We hypothesise that many NRG1 fusions may be oncogenic not for lacking the cytoplasmic domain but because they do not encode NRG1's nuclear-localised form. The fusion in MDA-MB-175 is the result of a very complex genomic rearrangement, which we partially characterised, that creates additional expressed gene fusions, RSF1-TENM4, TPCN2-RSF1, and MRPL48-GAB2. We searched for NRG1 rearrangements in 571 breast cancers subjected to genome sequencing and transcriptome sequencing and found four cases (0.7%) with fusions, WRN-NRG1, FAM91A1-NRG1, ARHGEF39-NRG1, and ZNF704-NRG1, all splicing into NRG1 at the same exon as in MDA-MB-175. However, the WRN-NRG1 and ARHGEF39-NRG1 fusions were out of frame. We identified rearrangements of NRG1 in many more (8% of) cases that seemed more likely to inactivate than to create activating fusions, or whose outcome could not be predicted because they were complex, or both. This is not surprising because NRG1 can be pro-apoptotic and is inactivated in some breast cancers. CONCLUSIONS: Our results highlight the complexity of rearrangements of NRG1 in breast cancers and confirm that some do not activate but inactivate. Careful interpretation of NRG1 rearrangements will therefore be necessary for appropriate patient management.

Imaging breast cancer using hyperpolarized carbon-13 MRI.

Our purpose is to investigate the feasibility of imaging tumor metabolism in breast cancer patients using 13C magnetic resonance spectroscopic imaging (MRSI) of hyperpolarized 13C label exchange between injected [1-13C]pyruvate and the endogenous tumor lactate pool. Treatment-naïve breast cancer patients were recruited: four triple-negative grade 3 cancers; two invasive ductal carcinomas that were estrogen and progesterone receptor-positive (ER/PR+) and HER2/neu-negative (HER2-), one grade 2 and one grade 3; and one grade 2 ER/PR+ HER2- invasive lobular carcinoma (ILC). Dynamic 13C MRSI was performed following injection of hyperpolarized [1-13C]pyruvate. Expression of lactate dehydrogenase A (LDHA), which catalyzes 13C label exchange between pyruvate and lactate, hypoxia-inducible factor-1 (HIF1α), and the monocarboxylate transporters MCT1 and MCT4 were quantified using immunohistochemistry and RNA sequencing. We have demonstrated the feasibility and safety of hyperpolarized 13C MRI in early breast cancer. Both intertumoral and intratumoral heterogeneity of the hyperpolarized pyruvate and lactate signals were observed. The lactate-to-pyruvate signal ratio (LAC/PYR) ranged from 0.021 to 0.473 across the tumor subtypes (mean ± SD: 0.145 ± 0.164), and a lactate signal was observed in all of the grade 3 tumors. The LAC/PYR was significantly correlated with tumor volume (R = 0.903, P = 0.005) and MCT 1 (R = 0.85, P = 0.032) and HIF1α expression (R = 0.83, P = 0.043). Imaging of hyperpolarized [1-13C]pyruvate metabolism in breast cancer is feasible and demonstrated significant intertumoral and intratumoral metabolic heterogeneity, where lactate labeling correlated with MCT1 expression and hypoxia.

Shared characteristics underpinning C4 leaf maturation derived from analysis of multiple C3 and C4 species of Flaveria.

Most terrestrial plants use C3 photosynthesis to fix carbon. In multiple plant lineages a modified system known as C4 photosynthesis has evolved. To better understand the molecular patterns associated with induction of C4 photosynthesis, the genus Flaveria that contains C3 and C4 species was used. A base to tip maturation gradient of leaf anatomy was defined, and RNA sequencing was undertaken along this gradient for two C3 and two C4 Flaveria species. Key C4 traits including vein density, mesophyll and bundle sheath cross-sectional area, chloroplast ultrastructure, and abundance of transcripts encoding proteins of C4 photosynthesis were quantified. Candidate genes underlying each of these C4 characteristics were identified. Principal components analysis indicated that leaf maturation and the photosynthetic pathway were responsible for the greatest amount of variation in transcript abundance. Photosynthesis genes were over-represented for a prolonged period in the C4 species. Through comparison with publicly available data sets, we identify a small number of transcriptional regulators that have been up-regulated in diverse C4 species. The analysis identifies similar patterns of expression in independent C4 lineages and so indicates that the complex C4 pathway is associated with parallel as well as convergent evolution.

Ancestral light and chloroplast regulation form the foundations for C4 gene expression.

C4 photosynthesis acts as a carbon concentrating mechanism that leads to large increases in photosynthetic efficiency. The C4 pathway is found in more than 60 plant lineages1 but the molecular enablers of this evolution are poorly understood. In particular, it is unclear how non-photosynthetic proteins in the ancestral C3 system have repeatedly become strongly expressed and integrated into photosynthesis gene regulatory networks in C4 leaves. Here, we provide clear evidence that in C3 leaves, genes encoding key enzymes of the C4 pathway are already co-regulated with photosynthesis genes and are controlled by both light and chloroplast-to-nucleus signalling. In C4 leaves this regulation becomes increasingly dependent on the chloroplast. We propose that regulation of C4 cycle genes by light and the chloroplast in the ancestral C3 state has facilitated the repeated evolution of the complex and convergent C4 trait.

TransRate: reference-free quality assessment of de novo transcriptome assemblies.

TransRate is a tool for reference-free quality assessment of de novo transcriptome assemblies. Using only the sequenced reads and the assembly as input, we show that multiple common artifacts of de novo transcriptome assembly can be readily detected. These include chimeras, structural errors, incomplete assembly, and base errors. TransRate evaluates these errors to produce a diagnostic quality score for each contig, and these contig scores are integrated to evaluate whole assemblies. Thus, TransRate can be used for de novo assembly filtering and optimization as well as comparison of assemblies generated using different methods from the same input reads. Applying the method to a data set of 155 published de novo transcriptome assemblies, we deconstruct the contribution that assembly method, read length, read quantity, and read quality make to the accuracy of de novo transcriptome assemblies and reveal that variance in the quality of the input data explains 43% of the variance in the quality of published de novo transcriptome assemblies. Because TransRate is reference-free, it is suitable for assessment of assemblies of all types of RNA, including assemblies of long noncoding RNA, rRNA, mRNA, and mixed RNA samples.

Transcript residency on ribosomes reveals a key role for the Arabidopsis thaliana bundle sheath in sulfur and glucosinolate metabolism.

Leaves of angiosperms are made up of multiple distinct cell types. While the function of mesophyll cells, guard cells, phloem companion cells and sieve elements are clearly described, this is not the case for the bundle sheath (BS). To provide insight into the role of the BS in the C3 species Arabidopsis thaliana, we labelled ribosomes in this cell type with a FLAG tag. We then used immunocapture to isolate these ribosomes, followed by sequencing of resident mRNAs. This showed that 5% of genes showed specific splice forms in the BS, and that 15% of genes were preferentially expressed in these cells. The BS translatome strongly implies that the BS plays specific roles in sulfur transport and metabolism, glucosinolate biosynthesis and trehalose metabolism. Much of the C4 cycle is differentially expressed between the C3 BS and the rest of the leaf. Furthermore, the global patterns of transcript residency on BS ribosomes overlap to a greater extent with cells of the root pericycle than any other cell type. This analysis provides the first insight into the molecular function of this cell type in C3 species, and also identifies characteristics of BS cells that are probably ancestral to both C3 and C4 plants.

The Pfam protein families database.

Pfam is a widely used database of protein families, currently containing more than 13,000 manually curated protein families as of release 26.0. Pfam is available via servers in the UK (http://pfam.sanger.ac.uk/), the USA (http://pfam.janelia.org/) and Sweden (http://pfam.sbc.su.se/). Here, we report on changes that have occurred since our 2010 NAR paper (release 24.0). Over the last 2 years, we have generated 1840 new families and increased coverage of the UniProt Knowledgebase (UniProtKB) to nearly 80%. Notably, we have taken the step of opening up the annotation of our families to the Wikipedia community, by linking Pfam families to relevant Wikipedia pages and encouraging the Pfam and Wikipedia communities to improve and expand those pages. We continue to improve the Pfam website and add new visualizations, such as the 'sunburst' representation of taxonomic distribution of families. In this work we additionally address two topics that will be of particular interest to the Pfam community. First, we explain the definition and use of family-specific, manually curated gathering thresholds. Second, we discuss some of the features of domains of unknown function (also known as DUFs), which constitute a rapidly growing class of families within Pfam.