In eukaryotes, the highly conserved U3 small nucleolar RNA (snoRNA) base-pairs to multiple sites in the pre-ribosomal RNA (pre-rRNA) to promote early cleavage and folding events. Binding of the U3 box A region to the pre-rRNA is mutually exclusive with folding of the central pseudoknot (CPK), a universally conserved rRNA structure of the small ribosomal subunit essential for protein synthesis. Here, we report that the DEAH-box helicase Dhr1 (Ecm16) is responsible for displacing U3. An active site mutant of Dhr1 blocked release of U3 from the pre-ribosome, thereby trapping a pre-40S particle. This particle had not yet achieved its mature structure because it contained U3, pre-rRNA, and a number of early-acting ribosome synthesis factors but noticeably lacked ribosomal proteins (r-proteins) that surround the CPK. Dhr1 was cross-linked in vivo to the pre-rRNA and to U3 sequences flanking regions that base-pair to the pre-rRNA including those that form the CPK. Point mutations in the box A region of U3 suppressed a cold-sensitive mutation of Dhr1, strongly indicating that U3 is an in vivo substrate of Dhr1. To support the conclusions derived from in vivo analysis we showed that Dhr1 unwinds U3-18S duplexes in vitro by using a mechanism reminiscent of DEAD box proteins.

Motile ciliopathies are characterized by specific defects in cilia beating that result in chronic airway disease, subfertility, ectopic pregnancy, and hydrocephalus. While many patients harbor mutations in the dynein motors that drive cilia beating, the disease also results from mutations in so-called dynein axonemal assembly factors (DNAAFs) that act in the cytoplasm. The mechanisms of DNAAF action remain poorly defined. Here, we show that DNAAFs concentrate together with axonemal dyneins and chaperones into organelles that form specifically in multiciliated cells, which we term DynAPs, for dynein axonemal particles. These organelles display hallmarks of biomolecular condensates, and remarkably, DynAPs are enriched for the stress granule protein G3bp1, but not for other stress granule proteins or P-body proteins. Finally, we show that both the formation and the liquid-like behaviors of DynAPs are disrupted in a model of motile ciliopathy. These findings provide a unifying cell biological framework for a poorly understood class of human disease genes and add motile ciliopathy to the growing roster of human diseases associated with disrupted biological phase separation.

The gram-negative opportunistic pathogen Pseudomonas aeruginosa is the primary cause of chronic respiratory infections in individuals with the heritable disease cystic fibrosis (CF). These infections can last for decades, during which time P. aeruginosa has been proposed to acquire beneficial traits via adaptive evolution. Because CF lacks an animal model that can acquire chronic P. aeruginosa infections, identifying genes important for long-term in vivo fitness remains difficult. However, since clonal, chronological samples can be obtained from chronically infected individuals, traits undergoing adaptive evolution can be identified. Recently we identified 24 P. aeruginosa gene expression traits undergoing parallel evolution in vivo in multiple individuals, suggesting they are beneficial to the bacterium. The goal of this study was to determine if these genes impact P. aeruginosa phenotypes important for survival in the CF lung. By using a gain-of-function genetic screen, we found that 4 genes and 2 operons undergoing parallel evolution in vivo promote P. aeruginosa biofilm formation. These genes/operons promote biofilm formation by increasing levels of the non-alginate exopolysaccharide Psl. One of these genes, phaF, enhances Psl production via a post-transcriptional mechanism, while the other 5 genes/operons do not act on either psl transcription or translation. Together, these data demonstrate that P. aeruginosa has evolved at least two pathways to over-produce a non-alginate exopolysaccharide during long-term colonization of the CF lung. More broadly, this approach allowed us to attribute a biological significance to genes with unknown function, demonstrating the power of using evolution as a guide for targeted genetic studies.

Characterizing the in vivo dynamics of the polyclonal antibody repertoire in serum, such as that which might arise in response to stimulation with an antigen, is difficult due to the presence of many highly similar immunoglobulin proteins, each specified by distinct B lymphocytes. These challenges have precluded the use of conventional mass spectrometry for antibody identification based on peptide mass spectral matches to a genomic reference database. Recently, progress has been made using bottom-up analysis of serum antibodies by nanoflow liquid chromatography/high-resolution tandem mass spectrometry combined with a sample-specific antibody sequence database generated by high-throughput sequencing of individual B cell immunoglobulin variable domains (V genes). Here, we describe how intrinsic features of antibody primary structure, most notably the interspersed segments of variable and conserved amino acid sequences, generate recurring patterns in the corresponding peptide mass spectra of V gene peptides, greatly complicating the assignment of correct sequences to mass spectral data. We show that the standard method of decoy-based error modeling fails to account for the error introduced by these highly similar sequences, leading to a significant underestimation of the false discovery rate. Because of these effects, antibody-derived peptide mass spectra require increased stringency in their interpretation. The use of filters based on the mean precursor ion mass accuracy of peptide-spectrum matches is shown to be particularly effective in distinguishing between "true" and "false" identifications. These findings highlight important caveats associated with the use of standard database search and error-modeling methods with nonstandard data sets and custom sequence databases.

Despite recent advances in human genetics, model organisms are indispensable for human disease research. Most human disease pathways are evolutionally conserved among other species, where they may phenocopy the human condition or be associated with seemingly unrelated phenotypes. Much of the known gene-to-phenotype association information is distributed across diverse databases, growing rapidly due to new experimental techniques. Accessible bioinformatics tools will therefore facilitate translation of discoveries from model organisms into human disease biology. Here, we present a web-based discovery tool for human disease studies, MORPHIN (model organisms projected on a human integrated gene network), which prioritizes the most relevant human diseases for a given set of model organism genes, potentially highlighting new model systems for human diseases and providing context to model organism studies. Conceptually, MORPHIN investigates human diseases by an orthology-based projection of a set of model organism genes onto a genome-scale human gene network. MORPHIN then prioritizes human diseases by relevance to the projected model organism genes using two distinct methods: a conventional overlap-based gene set enrichment analysis and a network-based measure of closeness between the query and disease gene sets capable of detecting associations undetectable by the conventional overlap-based methods. MORPHIN is freely accessible at http://www.inetbio.org/morphin.

To explore the origins and consequences of tetraploidy in the African clawed frog, we sequenced the Xenopus laevis genome and compared it to the related diploid X. tropicalis genome. We characterize the allotetraploid origin of X. laevis by partitioning its genome into two homoeologous subgenomes, marked by distinct families of 'fossil' transposable elements. On the basis of the activity of these elements and the age of hundreds of unitary pseudogenes, we estimate that the two diploid progenitor species diverged around 34 million years ago (Ma) and combined to form an allotetraploid around 17-18 Ma. More than 56% of all genes were retained in two homoeologous copies. Protein function, gene expression, and the amount of conserved flanking sequence all correlate with retention rates. The subgenomes have evolved asymmetrically, with one chromosome set more often preserving the ancestral state and the other experiencing more gene loss, deletion, rearrangement, and reduced gene expression.

GPR151 is a G-protein coupled receptor for which the endogenous ligand remains unknown. In the nervous system of vertebrates, its expression is enriched in specific diencephalic structures, where the highest levels are observed in the habenular area. The habenula has been implicated in a range of different functions including behavioral flexibility, decision making, inhibitory control, and pain processing, which makes it a promising target for treating psychiatric and neurological disease. This study aimed to further characterize neurons expressing the Gpr151 gene, by tracing the afferent connectivity of this diencephalic cell population. Using pseudotyped rabies virus in a transgenic Gpr151-Cre mouse line, monosynaptic afferents of habenular and thalamic Gpr151-expressing neuronal populations could be visualized. The habenular and thalamic Gpr151 systems displayed both shared and distinct connectivity patterns. The habenular neurons primarily received input from basal forebrain structures, the bed nucleus of stria terminalis, the lateral preoptic area, the entopeduncular nucleus, and the lateral hypothalamic area. The Gpr151-expressing neurons in the paraventricular nucleus of the thalamus was primarily contacted by medial hypothalamic areas as well as the zona incerta and projected to specific forebrain areas such as the prelimbic cortex and the accumbens nucleus. Gpr151 mRNA was also detected at low levels in the lateral posterior thalamic nucleus which received input from areas associated with visual processing, including the superior colliculus, zona incerta, and the visual and retrosplenial cortices. Knowledge about the connectivity of Gpr151-expressing neurons will facilitate the interpretation of future functional studies of this receptor.

The growing availability of large-scale functional networks has promoted the development of many successful techniques for predicting functions of genes. Here we extend these network-based principles and techniques to functionally characterize whole sets of genes. We present RIDDLE (Reflective Diffusion and Local Extension), which uses well developed guilt-by-association principles upon a human gene network to identify associations of gene sets. RIDDLE is particularly adept at characterizing sets with no annotations, a major challenge where most traditional set analyses fail. Notably, RIDDLE found microRNA-450a to be strongly implicated in ocular diseases and development. A web application is available at http://www.functionalnet.org/RIDDLE.

Myelination of the central nervous system (CNS) is critical to vertebrate nervous systems for efficient neural signaling. CNS myelination occurs as oligodendrocytes terminally differentiate, a process regulated in part by the myelin regulatory factor, MYRF. Using bioinformatics and extensive biochemical and functional assays, we find that MYRF is generated as an integral membrane protein that must be processed to release its transcription factor domain from the membrane. In contrast to most membrane-bound transcription factors, MYRF proteolysis seems constitutive and independent of cell- and tissue-type, as we demonstrate by reconstitution in E. coli and yeast. The apparent absence of physiological cues raises the question as to how and why MYRF is processed. By using computational methods capable of recognizing extremely divergent sequence homology, we identified a MYRF protein domain distantly related to bacteriophage tailspike proteins. Although occurring in otherwise unrelated proteins, the phage domains are known to chaperone the tailspike proteins' trimerization and auto-cleavage, raising the hypothesis that the MYRF domain might contribute to a novel activation method for a membrane-bound transcription factor. We find that the MYRF domain indeed serves as an intramolecular chaperone that facilitates MYRF trimerization and proteolysis. Functional assays confirm that the chaperone domain-mediated auto-proteolysis is essential both for MYRF's transcriptional activity and its ability to promote oligodendrocyte maturation. This work thus reveals a previously unknown key step in CNS myelination. These data also reconcile conflicting observations of this protein family, different members of which have been identified as transmembrane or nuclear proteins. Finally, our data illustrate a remarkable evolutionary repurposing between bacteriophages and eukaryotes, with a chaperone domain capable of catalyzing trimerization-dependent auto-proteolysis in two entirely distinct protein and cellular contexts, in one case participating in bacteriophage tailspike maturation and in the other activating a key transcription factor for CNS myelination.

Phenotypes and diseases may be related to seemingly dissimilar phenotypes in other species by means of the orthology of underlying genes. Such "orthologous phenotypes," or "phenologs," are examples of deep homology, and may be used to predict additional candidate disease genes.

In eukaryotic cells the final maturation of ribosomes occurs in the cytoplasm, where trans-acting factors are removed and critical ribosomal proteins are added for functionality. Here, we have carried out a comprehensive analysis of cytoplasmic maturation, ordering the known steps into a coherent pathway. Maturation is initiated by the ATPase Drg1. Downstream, assembly of the ribosome stalk is essential for the release of Tif6. The stalk recruits GTPases during translation. Because the GTPase Efl1, which is required for the release of Tif6, resembles the translation elongation factor eEF2, we suggest that assembly of the stalk recruits Efl1, triggering a step in 60S biogenesis that mimics aspects of translocation. Efl1 could thereby provide a mechanism to functionally check the nascent subunit. Finally, the release of Tif6 is a prerequisite for the release of the nuclear export adaptor Nmd3. Establishing this pathway provides an important conceptual framework for understanding ribosome maturation.

Several years after sequencing the human genome and the mouse genome, much remains to be discovered about the functions of most human and mouse genes. Computational prediction of gene function promises to help focus limited experimental resources on the most likely hypotheses. Several algorithms using diverse genomic data have been applied to this task in model organisms; however, the performance of such approaches in mammals has not yet been evaluated.

Post-translational modification (PTM) of proteins is central to many cellular processes across all domains of life, but despite decades of study and a wealth of genomic and proteomic data the biological function of many PTMs remains unknown. This is especially true for prokaryotic PTM systems, many of which have only recently been recognized and studied in depth. It is increasingly apparent that a deep sampling of abundance across a wide range of environmental stresses, growth conditions, and PTM types, rather than simply cataloging targets for a handful of modifications, is critical to understanding the complex pathways that govern PTM deposition and downstream effects.

Plants are foundational for global ecological and economic systems, but most plant proteins remain uncharacterized. Protein interaction networks often suggest protein functions and open new avenues to characterize genes and proteins. We therefore systematically determined protein complexes from 13 plant species of scientific and agricultural importance, greatly expanding the known repertoire of stable protein complexes in plants. By using co-fractionation mass spectrometry, we recovered known complexes, confirmed complexes predicted to occur in plants, and identified previously unknown interactions conserved over 1.1 billion years of green plant evolution. Several novel complexes are involved in vernalization and pathogen defense, traits critical for agriculture. We also observed plant analogs of animal complexes with distinct molecular assemblies, including a megadalton-scale tRNA multi-synthetase complex. The resulting map offers a cross-species view of conserved, stable protein assemblies shared across plant cells and provides a mechanistic, biochemical framework for interpreting plant genetics and mutant phenotypes.

Many gene families have been expanded by gene duplications along the human lineage, relative to ancestral opisthokonts, but the extent to which the duplicated genes function similarly is understudied. Here, we focused on structural cytoskeletal genes involved in critical cellular processes, including chromosome segregation, macromolecular transport, and cell shape maintenance. To determine functional redundancy and divergence of duplicated human genes, we systematically humanized the yeast actin, myosin, tubulin, and septin genes, testing ∼81% of human cytoskeletal genes across seven gene families for their ability to complement a growth defect induced by inactivation or deletion of the corresponding yeast ortholog. In five of seven families-all but α-tubulin and light myosin, we found at least one human gene capable of complementing loss of the yeast gene. Despite rescuing growth defects, we observed differential abilities of human genes to rescue cell morphology, meiosis, and mating defects. By comparing phenotypes of humanized strains with deletion phenotypes of their interaction partners, we identify instances of human genes in the actin and septin families capable of carrying out essential functions, but failing to fully complement the cytoskeletal roles of their yeast orthologs, thus leading to abnormal cell morphologies. Overall, we show that duplicated human cytoskeletal genes appear to have diverged such that only a few human genes within each family are capable of replacing the essential roles of their yeast orthologs. The resulting yeast strains with humanized cytoskeletal components now provide surrogate platforms to characterize human genes in simplified eukaryotic contexts.

The molecular and cellular basis of novelty is an active area of research in evolutionary biology. Until very recently, the vast majority of cellular phenomena were so difficult to sample that cross-species studies of biochemistry were rare and comparative analysis at the level of biochemical systems was almost impossible. Recent advances in systems biology are changing what is possible, however, and comparative phylogenetic methods that can handle this new data are wanted. Here, we introduce the term "phylogenetic latent variable models" (PLVMs, pronounced "plums") for a class of models that has recently been used to infer the evolution of cellular states from systems-level molecular data, and develop a new parameterization and fitting strategy that is useful for comparative inference of biochemical networks. We deploy this new framework to infer the ancestral states and evolutionary dynamics of protein-interaction networks by analyzing >16,000 predominantly metazoan co-fractionation and affinity-purification mass spectrometry experiments. Based on these data, we estimate ancestral interactions across unikonts, broadly recovering protein complexes involved in translation, transcription, proteostasis, transport, and membrane trafficking. Using these results, we predict an ancient core of the Commander complex made up of CCDC22, CCDC93, C16orf62, and DSCR3, with more recent additions of COMMD-containing proteins in tetrapods. We also use simulations to develop model fitting strategies and discuss future model developments.

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.

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.

Cilia are hairlike protrusions that project from the surface of eukaryotic cells and play key roles in cell signaling and motility. Ciliary motility is regulated by the conserved nexin-dynein regulatory complex (N-DRC), which links adjacent doublet microtubules and regulates and coordinates the activity of outer doublet complexes. Despite its critical role in cilia motility, the assembly and molecular basis of the regulatory mechanism are poorly understood. Here, using cryo-electron microscopy in conjunction with biochemical cross-linking and integrative modeling, we localize 12 DRC subunits in the N-DRC structure of Tetrahymena thermophila. We also find that the CCDC96/113 complex is in close contact with the DRC9/10 in the linker region. In addition, we reveal that the N-DRC is associated with a network of coiled-coil proteins that most likely mediates N-DRC regulatory activity.