Comprehensive Molecular Portraits of Invasive Lobular Breast Cancer.

Invasive lobular carcinoma (ILC) is the second most prevalent histologic subtype of invasive breast cancer. Here, we comprehensively profiled 817 breast tumors, including 127 ILC, 490 ductal (IDC), and 88 mixed IDC/ILC. Besides E-cadherin loss, the best known ILC genetic hallmark, we identified mutations targeting PTEN, TBX3, and FOXA1 as ILC enriched features. PTEN loss associated with increased AKT phosphorylation, which was highest in ILC among all breast cancer subtypes. Spatially clustered FOXA1 mutations correlated with increased FOXA1 expression and activity. Conversely, GATA3 mutations and high expression characterized luminal A IDC, suggesting differential modulation of ER activity in ILC and IDC. Proliferation and immune-related signatures determined three ILC transcriptional subtypes associated with survival differences. Mixed IDC/ILC cases were molecularly classified as ILC-like and IDC-like revealing no true hybrid features. This multidimensional molecular atlas sheds new light on the genetic bases of ILC and provides potential clinical options.

Seeking celestial positronium with an OH-suppressed diffraction-limited spectrograph.

Celestially, positronium (Ps) has been observed only through gamma-ray emission produced by its annihilation. However, in its triplet state, a Ps atom has a mean lifetime long enough for electronic transitions to occur between quantum states. This produces a recombination spectrum observable in principle at near IR wavelengths, where angular resolution greatly exceeding that of the gamma-ray observations is possible. However, the background in the near IR is dominated by extremely bright atmospheric hydroxyl (OH) emission lines. In this paper, we present the design of a diffraction-limited spectroscopic system using novel photonic components-a photonic lantern, OH fiber Bragg grating filters, and a photonic TIGER 2D pseudo-slit-to observe the Ps Balmer alpha line at 1.3122 µm for the first time, to our knowledge.

Large-scale profiling of microRNAs for The Cancer Genome Atlas.

The comprehensive multiplatform genomics data generated by The Cancer Genome Atlas (TCGA) Research Network is an enabling resource for cancer research. It includes an unprecedented amount of microRNA sequence data: ~11 000 libraries across 33 cancer types. Combined with initiatives like the National Cancer Institute Genomics Cloud Pilots, such data resources will make intensive analysis of large-scale cancer genomics data widely accessible. To support such initiatives, and to enable comparison of TCGA microRNA data to data from other projects, we describe the process that we developed and used to generate the microRNA sequence data, from library construction through to submission of data to repositories. In the context of this process, we describe the computational pipeline that we used to characterize microRNA expression across large patient cohorts.

Improved white spruce (Picea glauca) genome assemblies and annotation of large gene families of conifer terpenoid and phenolic defense metabolism.

White spruce (Picea glauca), a gymnosperm tree, has been established as one of the models for conifer genomics. We describe the draft genome assemblies of two white spruce genotypes, PG29 and WS77111, innovative tools for the assembly of very large genomes, and the conifer genomics resources developed in this process. The two white spruce genotypes originate from distant geographic regions of western (PG29) and eastern (WS77111) North America, and represent elite trees in two Canadian tree-breeding programs. We present an update (V3 and V4) for a previously reported PG29 V2 draft genome assembly and introduce a second white spruce genome assembly for genotype WS77111. Assemblies of the PG29 and WS77111 genomes confirm the reconstructed white spruce genome size in the 20 Gbp range, and show broad synteny. Using the PG29 V3 assembly and additional white spruce genomics and transcriptomics resources, we performed MAKER-P annotation and meticulous expert annotation of very large gene families of conifer defense metabolism, the terpene synthases and cytochrome P450s. We also comprehensively annotated the white spruce mevalonate, methylerythritol phosphate and phenylpropanoid pathways. These analyses highlighted the large extent of gene and pseudogene duplications in a conifer genome, in particular for genes of secondary (i.e. specialized) metabolism, and the potential for gain and loss of function for defense and adaptation.

Comparative genomics of the pine pathogens and beetle symbionts in the genus Grosmannia.

Studies on beetle/tree fungal symbionts typically characterize the ecological and geographic distributions of the fungal populations. There is limited understanding of the genome-wide evolutionary processes that act within and between species as such fungi adapt to different environments, leading to physiological differences and reproductive isolation. Here, we assess genomic evidence for such evolutionary processes by extending our recent work on Grosmannia clavigera, which is vectored by the mountain pine beetle and jeffrey pine beetle. We report the genome sequences of an additional 11 G. clavigera (Gc) sensu lato strains from the two known sibling species, Grosmannia sp. (Gs) and Gc. The 12 fungal genomes are structurally similar, showing large-scale synteny within and between species. We identified 103,430 single-nucleotide variations that separated the Grosmannia strains into divergent Gs and Gc clades, and further divided each of these clades into two subclades, one of which may represent an additional species. Comparing variable genes between these lineages, we identified truncated genes and potential pseudogenes, as well as seven genes that show evidence of positive selection. As these variable genes are involved in secondary metabolism and in detoxifying or utilizing host-tree defense chemicals (e.g., polyketide synthases, oxidoreductases, and mono-oxygenases), their variants may reflect adaptation to the specific chemistries of the host trees Pinus contorta, P. ponderosa, and P. jeffreyi. This work provides a comprehensive resource for developing informative markers for landscape population genomics of these ecologically and economically important fungi, and an approach that could be extended to other beetle-tree-associated fungi.

A Notch-dependent transcriptional hierarchy promotes mesenchymal transdifferentiation in the cardiac cushion.

BACKGROUND: Valvuloseptal defects are the most common congenital heart defects. Notch signaling-induced endothelial-to-mesenchymal transition (EMT) in the atrioventricular canal (AVC) cushions at murine embryonic day (E)9.5 is a required step during early valve development. Insights to the transcriptional network that is activated in endocardial cells (EC) during EMT and how these pathways direct valve maturation are lacking. RESULTS: We show that at E11.5, AVC-EC retain the ability to undergo Notch-dependent EMT when explanted on collagen. EC-Notch inhibition at E10.5 blocks expression of known mesenchymal genes in E11.5 AVC-EC. To understand the genetic network and AVC development downstream of Notch signaling beyond E9.5, we constructed Tag-Seq libraries corresponding to different cell types of the E11.5 AVC and atrium in wild-type mice and in EC-Notch inhibited mice. We identified 1,400 potential Notch targets in the AVC-EC, of which 124 are transcription factors (TF). From the 124 TFs, we constructed a transcriptional hierarchy and identify 10 upstream TFs within the network. CONCLUSIONS: We validated 4 of the upstream TFs as Notch targets that are enriched in AVC-EC. Functionally, we show these 4 TFs regulate EMT in AVC explant assays. These novel signaling pathways downstream of Notch are potentially relevant to valve development.

Identification and characterization of Hoxa9 binding sites in hematopoietic cells.

The clustered homeobox proteins play crucial roles in development, hematopoiesis, and leukemia, yet the targets they regulate and their mechanisms of action are poorly understood. Here, we identified the binding sites for Hoxa9 and the Hox cofactor Meis1 on a genome-wide level and profiled their associated epigenetic modifications and transcriptional targets. Hoxa9 and the Hox cofactor Meis1 cobind at hundreds of highly evolutionarily conserved sites, most of which are distant from transcription start sites. These sites show high levels of histone H3K4 monomethylation and CBP/P300 binding characteristic of enhancers. Furthermore, a subset of these sites shows enhancer activity in transient transfection assays. Many Hoxa9 and Meis1 binding sites are also bound by PU.1 and other lineage-restricted transcription factors previously implicated in establishment of myeloid enhancers. Conditional Hoxa9 activation is associated with CBP/P300 recruitment, histone acetylation, and transcriptional activation of a network of proto-oncogenes, including Erg, Flt3, Lmo2, Myb, and Sox4. Collectively, this work suggests that Hoxa9 regulates transcription by interacting with enhancers of genes important for hematopoiesis and leukemia.

Genome and transcriptome analyses of the mountain pine beetle-fungal symbiont Grosmannia clavigera, a lodgepole pine pathogen.

In western North America, the current outbreak of the mountain pine beetle (MPB) and its microbial associates has destroyed wide areas of lodgepole pine forest, including more than 16 million hectares in British Columbia. Grosmannia clavigera (Gc), a critical component of the outbreak, is a symbiont of the MPB and a pathogen of pine trees. To better understand the interactions between Gc, MPB, and lodgepole pine hosts, we sequenced the ∼30-Mb Gc genome and assembled it into 18 supercontigs. We predict 8,314 protein-coding genes, and support the gene models with proteome, expressed sequence tag, and RNA-seq data. We establish that Gc is heterothallic, and report evidence for repeat-induced point mutation. We report insights, from genome and transcriptome analyses, into how Gc tolerates conifer-defense chemicals, including oleoresin terpenoids, as they colonize a host tree. RNA-seq data indicate that terpenoids induce a substantial antimicrobial stress in Gc, and suggest that the fungus may detoxify these chemicals by using them as a carbon source. Terpenoid treatment strongly activated a ∼100-kb region of the Gc genome that contains a set of genes that may be important for detoxification of these host-defense chemicals. This work is a major step toward understanding the biological interactions between the tripartite MPB/fungus/forest system.

Improved Productivity with Multilaser Rotary Powder Bed Fusion Additive Manufacturing.

Laser powder bed fusion (LPBF) enables the fabrication of intricate, geometrically complex structures with a sufficiently fine surface finish for many engineering applications with a diversity of available feedstock metals. However, the production rate of LPBF systems is not well suited for mass production in comparison to traditional manufacturing methods. LPBF systems measure their deposition rates in 100's of grams per hour, while other processes measure in kilograms per hour or even in the case of processes such as forming, stamping, and casting, 100's of kilograms per hour. To be widely adopted in industry for mass production, LPBF requires a new scalable architecture that enables many orders of magnitude improvement in deposition rate, while maintaining the geometry freedom of additive manufacturing. This article explores concepts that could achieve as much as four orders of magnitude increase in the production rate through the application of (1) rotary table kinematic arrangements; (2) a dramatic number of simultaneously operating lasers; (3) reductions of laser optic size; (4) improved scanning techniques; and (5) an optimization of toroidal build plate size. To theoretically demonstrate the possibilities of production improvements, a productivity analysis is proposed for synchronous reluctance motors with relevance to the electric vehicle industry, given the recent increase in the diversity of printable soft magnetic alloys. The analysis provides insights into the impact of the architecture and process parameters necessary to optimize rotary powder bed fusion for mass production.

The genome and transcriptome of the pine saprophyte Ophiostoma piceae, and a comparison with the bark beetle-associated pine pathogen Grosmannia clavigera.

BACKGROUND: Ophiostoma piceae is a wood-staining fungus that grows in the sapwood of conifer logs and lumber. We sequenced its genome and analyzed its transcriptomes under a range of growth conditions. A comparison with the genome and transcriptomes of the mountain pine beetle-associated pathogen Grosmannia clavigera highlights differences between a pathogen that colonizes and kills living pine trees and a saprophyte that colonizes wood and the inner bark of dead trees. RESULTS: We assembled a 33 Mbp genome in 45 scaffolds, and predicted approximately 8,884 genes. The genome size and gene content were similar to those of other ascomycetes. Despite having similar ecological niches, O. piceae and G. clavigera showed no large-scale synteny. We identified O. piceae genes involved in the biosynthesis of melanin, which causes wood discoloration and reduces the commercial value of wood products. We also identified genes and pathways involved in growth on simple carbon sources and in sapwood, O. piceae's natural substrate. Like the pathogen, the saprophyte is able to tolerate terpenes, which are a major class of pine tree defense compounds; unlike the pathogen, it cannot utilize monoterpenes as a carbon source. CONCLUSIONS: This work makes available the second annotated genome of a softwood ophiostomatoid fungus, and suggests that O. piceae's tolerance to terpenes may be due in part to these chemicals being removed from the cells by an ABC transporter that is highly induced by terpenes. The data generated will provide the research community with resources for work on host-vector-fungus interactions for wood-inhabiting, beetle-associated saprophytes and pathogens.

Barnacle: detecting and characterizing tandem duplications and fusions in transcriptome assemblies.

BACKGROUND: Chimeric transcripts, including partial and internal tandem duplications (PTDs, ITDs) and gene fusions, are important in the detection, prognosis, and treatment of human cancers. RESULTS: We describe Barnacle, a production-grade analysis tool that detects such chimeras in de novo assemblies of RNA-seq data, and supports prioritizing them for review and validation by reporting the relative coverage of co-occurring chimeric and wild-type transcripts. We demonstrate applications in large-scale disease studies, by identifying PTDs in MLL, ITDs in FLT3, and reciprocal fusions between PML and RARA, in two deeply sequenced acute myeloid leukemia (AML) RNA-seq datasets. CONCLUSIONS: Our analyses of real and simulated data sets show that, with appropriate filter settings, Barnacle makes highly specific predictions for three types of chimeric transcripts that are important in a range of cancers: PTDs, ITDs, and fusions. High specificity makes manual review and validation efficient, which is necessary in large-scale disease studies. Characterizing an extended range of chimera types will help generate insights into progression, treatment, and outcomes for complex diseases.

Locus co-occupancy, nucleosome positioning, and H3K4me1 regulate the functionality of FOXA2-, HNF4A-, and PDX1-bound loci in islets and liver.

The liver and pancreas share a common origin and coexpress several transcription factors. To gain insight into the transcriptional networks regulating the function of these tissues, we globally identify binding sites for FOXA2 in adult mouse islets and liver, PDX1 in islets, and HNF4A in liver. Because most eukaryotic transcription factors bind thousands of loci, many of which are thought to be inactive, methods that can discriminate functionally active binding events are essential for the interpretation of genome-wide transcription factor binding data. To develop such a method, we also generated genome-wide H3K4me1 and H3K4me3 localization data in these tissues. By analyzing our binding and histone methylation data in combination with comprehensive gene expression data, we show that H3K4me1 enrichment profiles discriminate transcription factor occupied loci into three classes: those that are functionally active, those that are poised for activation, and those that reflect pioneer-like transcription factor activity. Furthermore, we demonstrate that the regulated presence of H3K4me1-marked nucleosomes at transcription factor occupied promoters and enhancers controls their activity, implicating both tissue-specific transcription factor binding and nucleosome remodeling complex recruitment in determining tissue-specific gene expression. Finally, we apply these approaches to generate novel insights into how FOXA2, PDX1, and HNF4A cooperate to drive islet- and liver-specific gene expression.

Alternative expression analysis by RNA sequencing.

In alternative expression analysis by sequencing (ALEXA-seq), we developed a method to analyze massively parallel RNA sequence data to catalog transcripts and assess differential and alternative expression of known and predicted mRNA isoforms in cells and tissues. As proof of principle, we used the approach to compare fluorouracil-resistant and -nonresistant human colorectal cancer cell lines. We assessed the sensitivity and specificity of the approach by comparison to exon tiling and splicing microarrays and validated the results with reverse transcription-PCR, quantitative PCR and Sanger sequencing. We observed global disruption of splicing in fluorouracil-resistant cells characterized by expression of new mRNA isoforms resulting from exon skipping, alternative splice site usage and intron retention. Alternative expression annotation databases, source code, a data viewer and other resources to facilitate analysis are available at http://www.alexaplatform.org/alexa_seq/.

De novo assembly and analysis of RNA-seq data.

We describe Trans-ABySS, a de novo short-read transcriptome assembly and analysis pipeline that addresses variation in local read densities by assembling read substrings with varying stringencies and then merging the resulting contigs before analysis. Analyzing 7.4 gigabases of 50-base-pair paired-end Illumina reads from an adult mouse liver poly(A) RNA library, we identified known, new and alternative structures in expressed transcripts, and achieved high sensitivity and specificity relative to reference-based assembly methods.

A specialized ABC efflux transporter GcABC-G1 confers monoterpene resistance to Grosmannia clavigera, a bark beetle-associated fungal pathogen of pine trees.

Grosmannia clavigera is a bark beetle-vectored pine pathogen in the mountain pine beetle epidemic in western North America. Grosmannia clavigera colonizes pines despite the trees' massive oleoresin terpenoid defences. We are using a functional genomics approach to identify G. clavigera's mechanisms of adaptation to pine defences. We annotated the ABC transporters in the G. clavigera genome and generated RNA-seq transcriptomes from G. clavigera grown with a range of terpenes. We functionally characterized GcABC-G1, a pleiotropic drug resistance (PDR) transporter that was highly induced by terpenes, using qRT-PCR, gene knock-out and heterologous expression in yeast. Deleting GcABC-G1 increased G. clavigera's sensitivity to monoterpenes and delayed development of symptoms in inoculated young lodgepole pine trees. Heterologous expression of GcABC-G1 in yeast increased tolerance to monoterpenes. G. clavigera but not the deletion mutant, can use (+)-limonene as a carbon source. Phylogenetic analysis placed GcABC-G1 outside the ascomycete PDR transporter clades. G. clavigera appears to have evolved two mechanisms to survive and grow when exposed to monoterpenes: GcABC-G1 controls monoterpene levels within the fungal cells and G. clavigera uses monoterpenes as a carbon source. This work has implications for understanding adaptation to host defences in an important forest insect-fungal system, and potentially for metabolic engineering of terpenoid production in yeast.

Hexabundles: imaging fiber arrays for low-light astronomical applications.

We demonstrate a novel imaging fiber bundle ("hexabundle") that is suitable for low-light applications in astronomy. The most successful survey instruments at optical-infrared wavelengths use hundreds to thousands of multimode fibers fed to one or more spectrographs. Since most celestial sources are spatially extended on the celestial sphere, a hexabundle provides spectroscopic information at many distinct locations across the source. We discuss two varieties of hexabundles: (i) lightly fused, closely packed, circular cores; (ii) heavily fused non-circular cores with higher fill fractions. In both cases, we find the important result that the cladding can be reduced to ~2 µm over the short fuse length, well below the conventional ~10λ thickness employed more generally, with a consequent gain in fill factor. Over the coming decade, it is to be expected that fiber-based instruments will be upgraded with hexabundles in order to increase the spatial multiplex capability by two or more orders of magnitude.

PICS: probabilistic inference for ChIP-seq.

ChIP-seq combines chromatin immunoprecipitation with massively parallel short-read sequencing. While it can profile genome-wide in vivo transcription factor-DNA association with higher sensitivity, specificity, and spatial resolution than ChIP-chip, it poses new challenges for statistical analysis that derive from the complexity of the biological systems characterized and from variability and biases in its sequence data. We propose a method called PICS (Probabilistic Inference for ChIP-seq) for identifying regions bound by transcription factors from aligned reads. PICS identifies binding event locations by modeling local concentrations of directional reads, and uses DNA fragment length prior information to discriminate closely adjacent binding events via a Bayesian hierarchical t-mixture model. It uses precalculated, whole-genome read mappability profiles and a truncated t-distribution to adjust binding event models for reads that are missing due to local genome repetitiveness. It estimates uncertainties in model parameters that can be used to define confidence regions on binding event locations and to filter estimates. Finally, PICS calculates a per-event enrichment score relative to a control sample, and can use a control sample to estimate a false discovery rate. Using published GABP and FOXA1 data from human cell lines, we show that PICS' predicted binding sites were more consistent with computationally predicted binding motifs than the alternative methods MACS, QuEST, CisGenome, and USeq. We then use a simulation study to confirm that PICS compares favorably to these methods and is robust to model misspecification.

Caenorhabditis elegans cisRED: a catalogue of conserved genomic elements.

The availability of completely sequenced genomes from eight species of nematodes has provided an opportunity to identify novel cis-regulatory elements in the promoter regions of Caenorhabditis elegans transcripts using comparative genomics. We determined orthologues for C. elegans transcripts in C. briggsae, C. remanei, C. brenneri, C. japonica, Pristionchus pacificus, Brugia malayi and Trichinella spiralis using the WABA alignment algorithm. We pooled the upstream region of each transcript in C. elegans with the upstream regions of its orthologues and identified conserved DNA sequence elements by de novo motif discovery. In total, we discovered 158 017 novel conserved motifs upstream of 3847 C. elegans transcripts for which three or more orthologues were available, and identified 82% of 44 experimentally proven regulatory elements from ORegAnno. We annotated 26% of the motifs as similar to known binding sequences of transcription factors from ORegAnno, TRANSFAC and JASPAR. This is the first catalogue of annotated conserved upstream elements for nematodes and can be used to find putative regulatory elements, improve gene models, discover novel RNA genes, and understand the evolution of transcription factors and their binding sites in phylum Nematoda. The annotated motifs provide novel binding site candidates for both characterized transcription factors and orthologues of characterized mammalian transcription factors.

Gene discovery for the bark beetle-vectored fungal tree pathogen Grosmannia clavigera.

BACKGROUND: Grosmannia clavigera is a bark beetle-vectored fungal pathogen of pines that causes wood discoloration and may kill trees by disrupting nutrient and water transport. Trees respond to attacks from beetles and associated fungi by releasing terpenoid and phenolic defense compounds. It is unclear which genes are important for G. clavigera's ability to overcome antifungal pine terpenoids and phenolics. RESULTS: We constructed seven cDNA libraries from eight G. clavigera isolates grown under various culture conditions, and Sanger sequenced the 5' and 3' ends of 25,000 cDNA clones, resulting in 44,288 high quality ESTs. The assembled dataset of unique transcripts (unigenes) consists of 6,265 contigs and 2,459 singletons that mapped to 6,467 locations on the G. clavigera reference genome, representing ~70% of the predicted G. clavigera genes. Although only 54% of the unigenes matched characterized proteins at the NCBI database, this dataset extensively covers major metabolic pathways, cellular processes, and genes necessary for response to environmental stimuli and genetic information processing. Furthermore, we identified genes expressed in spores prior to germination, and genes involved in response to treatment with lodgepole pine phloem extract (LPPE). CONCLUSIONS: We provide a comprehensively annotated EST dataset for G. clavigera that represents a rich resource for gene characterization in this and other ophiostomatoid fungi. Genes expressed in response to LPPE treatment are indicative of fungal oxidative stress response. We identified two clusters of potentially functionally related genes responsive to LPPE treatment. Furthermore, we report a simple method for identifying contig misassemblies in de novo assembled EST collections caused by gene overlap on the genome.

Global analysis of in vivo Foxa2-binding sites in mouse adult liver using massively parallel sequencing.

Foxa2 (HNF3 beta) is a one of three, closely related transcription factors that are critical to the development and function of the mouse liver. We have used chromatin immunoprecipitation and massively parallel Illumina 1G sequencing (ChIP-Seq) to create a genome-wide profile of in vivo Foxa2-binding sites in the adult liver. More than 65% of the approximately 11.5 k genomic sites associated with Foxa2 binding, mapped to extended gene regions of annotated genes, while more than 30% of intragenic sites were located within first introns. 20.5% of all sites were further than 50 kb from any annotated gene, suggesting an association with novel gene regions. QPCR analysis demonstrated a strong positive correlation between peak height and fold enrichment for Foxa2-binding sites. We measured the relationship between Foxa2 and liver gene expression by overlapping Foxa2-binding sites with a SAGE transcriptome profile, and found that 43.5% of genes expressed in the liver were also associated with Foxa2 binding. We also identified potential Foxa2-interacting transcription factors whose motifs were enriched near Foxa2-binding sites. Our comprehensive results for in vivo Foxa2-binding sites in the mouse liver will contribute to resolving transcriptional regulatory networks that are important for adult liver function.