Proteins modulate signaling pathways through dynamic protein-protein interactions, which provide critical insights into the underlying mechanisms of cellular processes leading to human diseases. Our previous study revealed the biochemical and genetic interactions between FUZ and GPR161 that regulate sonic hedgehog signaling during spinal neural tube development. This study aimed to identify novel co-interacting proteins of FUZ and GPR161 and explore their biochemical and functional associations. We performed affinity-based liquid chromatography tandem mass spectrometry on immunoprecipitated FUZ and GPR161 proteins. Using this approach, we identified 159 co-interacting proteins, of which 289 proteins that interacted exclusively with FUZ and 617 proteins that interacted exclusively with GPR161. Gene Ontology (GO) analysis of the FUZ and GPR161 co-interactome revealed an enrichment of proteins associated with proteasomal catabolic processes and trafficking. GO analysis of the exclusive FUZ and GPR161 interactomes identified cell cycle progression, mitochondrial membrane, RNA metabolism, receptor complex, and Endoplasmic Reticulum-Golgi transport. These findings were further validated using STRING network analysis. Among the identified proteins, we prioritized FKBP8 and confirmed its biochemical interactions with both FUZ and GPR161 exclusively. In summary, our proteomic profiling uncovered the protein network of FUZ and GPR161, revealing their individual and cooperative functions in multiple cellular processes.

Studying the genetic and molecular characteristics of brewing yeast strains is crucial for understanding their domestication history and adaptations accumulated over time in fermentation environments, and for guiding optimizations to the brewing process itself. Saccharomyces cerevisiae (brewing yeast) is among the most profiled organisms on the planet, yet the temporal molecular changes that underlie industrial fermentation and beer brewing remain understudied. Here, we characterized the genomic makeup of a Saccharomyces cerevisiae ale yeast widely used in the production of Hefeweizen beers, and applied shotgun mass spectrometry to systematically measure the proteomic changes throughout 2 fermentation cycles which were separated by 14 rounds of serial repitching. The resulting brewing yeast proteomics resource includes 64,740 protein abundance measurements. We found that this strain possesses typical genetic characteristics of Saccharomyces cerevisiae ale strains and displayed progressive shifts in molecular processes during fermentation based on protein abundance changes. We observed protein abundance differences between early fermentation batches compared to those separated by 14 rounds of serial repitching. The observed abundance differences occurred mainly in proteins involved in the metabolism of ergosterol and isobutyraldehyde. Our systematic profiling serves as a starting point for deeper characterization of how the yeast proteome changes during commercial fermentations and additionally serves as a resource to guide fermentation protocols, strain handling, and engineering practices in commercial brewing and fermentation environments. Finally, we created a web interface (https://brewing-yeast-proteomics.ccbb.utexas.edu/) to serve as a valuable resource for yeast geneticists, brewers, and biochemists to provide insights into the global trends underlying commercial beer production.

All eukaryotes share a common ancestor from roughly 1.5 - 1.8 billion years ago, a single-celled, swimming microbe known as LECA, the Last Eukaryotic Common Ancestor. Nearly half of the genes in modern eukaryotes were present in LECA, and many current genetic diseases and traits stem from these ancient molecular systems. To better understand these systems, we compared genes across modern organisms and identified a core set of 10,092 shared protein-coding gene families likely present in LECA, a quarter of which are uncharacterized. We then integrated >26,000 mass spectrometry proteomics analyses from 31 species to infer how these proteins interact in higher-order complexes. The resulting interactome describes the biochemical organization of LECA, revealing both known and new assemblies. We analyzed these ancient protein interactions to find new human gene-disease relationships for bone density and congenital birth defects, demonstrating the value of ancestral protein interactions for guiding functional genetics today.

Replication of the complex retrovirus mouse mammary tumor virus (MMTV) is antagonized by murine Apobec3 (mA3), a member of the Apobec family of cytidine deaminases. We have shown that MMTV-encoded Rem protein inhibits proviral mutagenesis by the Apobec enzyme, activation-induced cytidine deaminase (AID) during viral replication in BALB/c mice. To further study the role of Rem in vivo, we have infected C57BL/6 (B6) mice with a superantigen-independent lymphomagenic strain of MMTV (TBLV-WT) or a mutant strain that is defective in Rem and its cleavage product Rem-CT (TBLV-SD). Compared to BALB/c, B6 mice were more susceptible to TBLV infection and tumorigenesis. Furthermore, unlike MMTV, TBLV induced T-cell tumors in B6 μMT mice, which lack membrane-bound IgM and conventional B-2 cells. At limiting viral doses, loss of Rem expression in TBLV-SD-infected B6 mice accelerated tumorigenesis compared to TBLV-WT in either wild-type B6 or AID-knockout mice. Unlike BALB/c results, high-throughput sequencing indicated that proviral G-to-A or C-to-T mutations were unchanged regardless of Rem expression in B6 tumors. However, knockout of both AID and mA3 reduced G-to-A mutations. Ex vivo stimulation showed higher levels of mA3 relative to AID in B6 compared to BALB/c splenocytes, and effects of agonists differed in the two strains. RNA-Seq revealed increased transcripts related to growth factor and cytokine signaling in TBLV-SD-induced tumors relative to TBLV-WT-induced tumors, consistent with another Rem function. Thus, Rem-mediated effects on tumorigenesis in B6 mice are independent of Apobec-mediated proviral hypermutation.

Southern flounder skin pigmentation is a critical phenotypic characteristic for this species’ survival in the natural environment. Normal pigmentation allows rapid changes of color for concealment to capture prey and UV light protection. In contrast, highly visible hypopigmented pseudo-albinos exhibit a compromised immune system and are vulnerable to predation, sensitive to UV exposure, and likely have poor survival in the wild. Skin and brain tissue samples from normally pigmented and hypopigmented individuals were analyzed with next-generation RNA sequencing. A total of 1,589,613 transcripts were used to identify 952,825 genes to assemble a de novo transcriptome, with 99.43% of genes mapped to the assembly. Differential gene expression and gene enrichment analysis of contrasting tissues and phenotypes revealed that pseudo-albino individuals appeared more susceptible to environmental stress, UV light exposure, hypoxia, and osmotic stress. The pseudo-albinos’ restricted immune response showed upregulated genes linked to cancer development, signaling and response, skin tissue formation, regeneration, and healing. The data indicate that a modified skin collagen structure likely affects melanocyte differentiation and distribution, generating the pseudo-albino phenotype. In addition, the comparison of the brain transcriptome revealed changes in myelination and melanocyte stem cell activity, which may indicate modified brain function, reduced melanocyte migration, and impaired vision.

AbstractReplication of the complex retrovirus mouse mammary tumor virus (MMTV) is antagonized by murine Apobec3 (mA3), a member of the Apobec family of cytidine deaminases. We have shown that MMTV-encoded Rem protein inhibits proviral mutagenesis by the Apobec enzyme, activation-induced cytidine deaminase (AID) during viral replication in BALB/c mice. To further study the role of Remin vivo, we have infected C57BL/6 (B6) mice with a superantigen-independent lymphomagenic strain of MMTV (TBLV-WT) or a mutant strain (TBLV-SD) that is defective in Rem and its cleavage product Rem-CT. Unlike MMTV, TBLV induced T-cell tumors in µMT mice, indicating that mature B cells, which express the highest AID levels, are not required for TBLV replication. Compared to BALB/c, B6 mice were more susceptible to TBLV infection and tumorigenesis. The lack of Rem expression accelerated B6 tumorigenesis at limiting doses compared to TBLV-WT in either wild-type B6 or AID-deficient mice. However, unlike proviruses from BALB/c mice, high-throughput sequencing indicated that proviral G-to-A or C-to-T changes did not significantly differ in the presence and absence of Rem expression.Ex vivostimulation showed higher levels of mA3 relative to AID in B6 compared to BALB/c splenocytes, but effects of agonists differed in the two strains. RNA-Seq revealed increased transcripts related to growth factor and cytokine signaling in TBLV-SD-induced tumors relative to those from TBLV-WT, consistent with a third Rem function. Thus, Rem-mediated effects on tumorigenesis in B6 mice are independent of Apobec-mediated proviral hypermutation.

Merged Tuxedo Suite output for tumor data, combining the cuffdiff and cuffnorm tables by unique ENSEMBL identifier

<p>Bivalent normalized expression in FPKM of the 840 genes identified by DAVID in Fig. 6A.</p>

Many cancers, including melanoma, have a higher requirement for l-methionine in comparison with noncancerous cells. In this study, we show that administration of an engineered human methionine-& gamma;-lyase (hMGL) significantly reduced the survival of both human and mouse melanoma cells in vitro. A multiomics approach was utilized to identify global changes in gene expression and in metabolite levels with hMGL treatment in melanoma cells. There was considerable overlap in the perturbed pathways identified in the two data sets. Common pathways were flagged for further investigation to understand their mechanistic importance. In this regard, hMGL treatment induced S and G2 phase cell cycle arrest, decreased nucleotide levels, and increased DNA double-strand breaks suggesting an important role for replication stress in the mechanism of hMGL effects on melanoma cells. Further, hMGL treatment resulted in increased cellular reactive oxygen species levels and increased apoptosis as well as uncharged transfer RNA pathway upregulation. Finally, treatment with hMGL significantly inhibited the growth of both mouse and human melanoma cells in orthotopic tumor models in vivo. Overall, the results of this study provide a strong rationale for further mechanistic evaluation and clinical development of hMGL for the treatment of melanoma skin cancer and other cancers.

Table S1: Output of ChromHMM and color key. Table S2: Enhancer counts. Table S3: Enhancer distribution. Table S4: Bivalent region counts.

Figure S1: Copy number alterations in GBM tumor genomes. Figure S2: Clustering of primary GBM tumors relative to GBM-derived cell lines. Figure S3: Gene expression over all 21 chromatin states in GBM tumors. Figure S4: Methylation over all identified chromatin states, over the 11 tumors analyzed in this study. Figure S5: 5hmC signal over enhancers, enhancer distance to genes, average expression of enhancer-associated genes by tumor. Figure S6: Sequence motifs in enhancers. Figure S7: Gene expression and chromatin signal for frequently bivalent domains and associated genes. Figure S8: Histone mark profiles over the HOXB locus, comparing GBM data with Roadmap Epigenomics brain data. Figure S9: Bivalent regions in T98G cells and their overlap with bivalent regions identified in at least 8 tumors.

<p>sequencing statistics for all datasets in this study.</p>

Enhancer normalized expression in FPKM of the 307 genes identified by DAVID in Fig. 5A.

The ongoing evolution of SARS-CoV-2 into more easily transmissible and infectious variants has sparked concern over the continued effectiveness of existing therapeutic antibodies and vaccines. Hence, together with increased genomic surveillance, methods to rapidly develop and assess effective interventions are critically needed. Here we report the discovery of SARS-CoV-2 neutralizing antibodies isolated from COVID-19 patients using a high-throughput platform. Antibodies were identified from unpaired donor B-cell and serum repertoires using yeast surface display, proteomics, and public light chain screening. Cryo-EM and functional characterization of the antibodies identified N3-1, an antibody that binds avidly (Kd,app = 68 pM) to the receptor binding domain (RBD) of the spike protein and robustly neutralizes the virus in vitro. This antibody likely binds all three RBDs of the trimeric spike protein with a single IgG. Importantly, N3-1 equivalently binds spike proteins from emerging SARS-CoV-2 variants of concern, neutralizes UK variant B.1.1.7, and binds SARS-CoV spike with nanomolar affinity. Taken together, the strategies described herein will prove broadly applicable in interrogating adaptive immunity and developing rapid response biological countermeasures to emerging pathogens.

Female fertility in mammals requires iterative remodeling of the entire adult female reproductive tract across the menstrual/estrous cycle. However, while transcriptome dynamics across the estrous cycle have been reported in human and bovine models, no global analysis of gene expression across the estrous cycle has yet been reported for the mouse. Here, we examined the cellular composition and global transcriptional dynamics of the mouse oviduct along the anteroposterior axis and across the estrous cycle. We observed robust patterns of differential gene expression along the anteroposterior axis, but we found surprisingly few changes in gene expression across the estrous cycle. Notable gene expression differences along the anteroposterior axis included a surprising enrichment for genes related to embryonic development, such as Hox and Wnt genes. The relatively stable transcriptional dynamics across the estrous cycle differ markedly from other mammals, leading us to speculate that this is an evolutionarily derived state that may reflect the extremely rapid five-day mouse estrous cycle. This dataset fills a critical gap by providing an important genomic resource for a highly tractable genetic model of mammalian female reproduction.

Deciphering the mechanisms that regulate the sensitivity of pathogen recognition receptors is imperative to understanding infection and inflammation. Here we demonstrate that the RNA triphosphatase dual-specificity phosphatase 11 (DUSP11) acts on both host and virus-derived 5'-triphosphate RNAs rendering them less active in inducing a RIG-I-mediated immune response. Reducing DUSP11 levels alters host triphosphate RNA packaged in extracellular vesicles and induces enhanced RIG-I activation in cells exposed to extracellular vesicles. Virus infection of cells lacking DUSP11 results in a higher proportion of triphosphorylated viral transcripts and attenuated virus replication, which is rescued by reducing RIG-I expression. Consistent with the activity of DUSP11 in the cellular RIG-I response, mice lacking DUSP11 display lower viral loads, greater sensitivity to triphosphorylated RNA, and a signature of enhanced interferon activity in select tissues. Our results reveal the importance of controlling 5'-triphosphate RNA levels to prevent aberrant RIG-I signaling and demonstrate DUSP11 as a key effector of this mechanism.

Transcription start sites (TSS) in eukaryotes are characterized by a nucleosome-depleted region (NDR), which appears to be flanked upstream and downstream by strongly positioned nucleosomes incorporating the histone variant H2A.Z. H2A.Z associates with both active and repressed TSS and is important for priming genes for rapid transcriptional activation. However, the determinants of H2A.Z occupancy at specific nucleosomes and its relationship to transcription initiation remain unclear. To further elucidate the specificity of H2A.Z, we determined its genomic localization at single nucleosome resolution, as well as the localization of its chromatin remodelers Swr1 and Ino80. By analyzing H2A.Z occupancy in conjunction with RNA expression data that captures promoter-derived antisense initiation, we find that H2A.Z's bimodal incorporation on either side of the NDR is not a general feature of TSS, but is specifically a marker for bidirectional transcription, such that the upstream flanking -1 H2A.Z-containing nucleosome is more appropriately considered as a +1 H2A.Z nucleosome for antisense transcription. The localization of H2A.Z almost exclusively at the +1 nucleosome suggests that a transcription-initiation dependent process could contribute to its specific incorporation.

Bromodomain proteins (BRD) are key chromatin regulators of genome function and stability as well as therapeutic targets in cancer. Here, we systematically delineate the contribution of human BRD proteins for genome stability and DNA double-strand break (DSB) repair using several cell-based assays and proteomic interaction network analysis. Applying these approaches, we identify 24 of the 42 BRD proteins as promoters of DNA repair and/or genome integrity. We identified a BRD-reader function of PCAF that bound TIP60-mediated histone acetylations at DSBs to recruit a DUB complex to deubiquitylate histone H2BK120, to allowing direct acetylation by PCAF, and repair of DSBs by homologous recombination. We also discovered the bromo-and-extra-terminal (BET) BRD proteins, BRD2 and BRD4, as negative regulators of transcription-associated RNA-DNA hybrids (R-loops) as inhibition of BRD2 or BRD4 increased R-loop formation, which generated DSBs. These breaks were reliant on topoisomerase II, and BRD2 directly bound and activated topoisomerase I, a known restrainer of R-loops. Thus, comprehensive interactome and functional profiling of BRD proteins revealed new homologous recombination and genome stability pathways, providing a framework to understand genome maintenance by BRD proteins and the effects of their pharmacological inhibition.

RNA-binding proteins (RBPs) play essential roles in biology and are frequently associated with human disease. Although recent studies have systematically identified individual RNA-binding proteins, their higher-order assembly into ribonucleoprotein (RNP) complexes has not been systematically investigated. Here, we describe a proteomics method for systematic identification of RNP complexes in human cells. We identify 1,428 protein complexes that associate with RNA, indicating that more than 20% of known human protein complexes contain RNA. To explore the role of RNA in the assembly of each complex, we identify complexes that dissociate, change composition, or form stable protein-only complexes in the absence of RNA. We use our method to systematically identify cell-type-specific RNA-associated proteins in mouse embryonic stem cells and finally, distribute our resource, rna.MAP, in an easy-to-use online interface (rna.proteincomplexes.org). Our system thus provides a methodology for explorations across human tissues, disease states, and throughout all domains of life.

To detect functional somatic mutations in tumor samples, whole-exome sequencing (WES) is often used for its reliability and relative low cost.RNA-seq, while generally used to measure gene expression, can potentially also be used for identification of somatic mutations.However there has been little systematic evaluation of the utility of RNA-seq for identifying somatic mutations.Here, we develop and evaluate a pipeline for processing RNA-seq data from glioblastoma multiforme (GBM) tumors in order to identify somatic mutations.The pipeline entails the use of the STAR aligner 2-pass procedure jointly with MuTect2 from GATK to detect somatic variants.Variants identified from RNA-seq data were evaluated by comparison against the COSMIC and dbSNP databases, and also compared to somatic variants identified by exome sequencing.We also estimated the putative functional impact of coding variants in the most frequently mutated genes in GBM.Interestingly, variants identified by RNA-seq alone showed better representation of GBMrelated mutations cataloged by COSMIC.RNA-seq-only data substantially outperformed the ability of WES to reveal potentially new somatic mutations in known GBM-related pathways, and allowed us to build a high-quality set of somatic mutations common to exome and RNA-seq calls.Using RNA-seq data in parallel with WES data to detect somatic mutations in cancer genomes can thus broaden the scope of discoveries and lend additional support to somatic variants identified by exome sequencing alone.