Multisite Substrate Recognition in Asf1-Dependent Acetylation of Histone H3 K56 by Rtt109.

Rtt109 is a unique histone acetyltransferase acetylating histone H3 lysine 56 (H3K56), a modification critical for DNA replication-coupled nucleosome assembly and genome stability. In cells, histone chaperone Asf1 is essential for H3K56 acetylation, yet the mechanisms for H3K56 specificity and Asf1 requirement remain unknown. We have determined the crystal structure of the Rtt109-Asf1-H3-H4 complex and found that unwinding of histone H3 αN, where K56 is normally located, and stabilization of the very C-terminal ß strand of histone H4 by Asf1 are prerequisites for H3K56 acetylation. Unexpectedly, an interaction between Rtt109 and the central helix of histone H3 is also required. The observed multiprotein, multisite substrate recognition mechanism among histone modification enzymes provides mechanistic understandings of Rtt109 and Asf1 in H3K56 acetylation, as well as valuable insights into substrate recognition by histone modification enzymes in general.

Origin and evolutionary malleability of T cell receptor α diversity.

Lymphocytes of vertebrate adaptive immune systems acquired the capability to assemble, from split genes in the germline, billions of functional antigen receptors1-3. These receptors show specificity; unlike the broadly tuned receptors of the innate system, antibodies (Ig) expressed by B cells, for instance, can accurately distinguish between the two enantiomers of organic acids4, whereas T cell receptors (TCRs) reliably recognize single amino acid replacements in their peptide antigens5. In developing lymphocytes, antigen receptor genes are assembled from a comparatively small set of germline-encoded genetic elements in a process referred to as V(D)J recombination6,7. Potential self-reactivity of some antigen receptors arising from the quasi-random somatic diversification is suppressed by several robust control mechanisms8-12. For decades, scientists have puzzled over the evolutionary origin of somatically diversifying antigen receptors13-16. It has remained unclear how, at the inception of this mechanism, immunologically beneficial expanded receptor diversity was traded against the emerging risk of destructive self-recognition. Here we explore the hypothesis that in early vertebrates, sequence microhomologies marking the ends of recombining elements became the crucial targets of selection determining the outcome of non-homologous end joining-based repair of DNA double-strand breaks generated during RAG-mediated recombination. We find that, across the main clades of jawed vertebrates, TCRα repertoire diversity is best explained by species-specific extents of such sequence microhomologies. Thus, selection of germline sequence composition of rearranging elements emerges as a major factor determining the degree of diversity of somatically generated antigen receptors.

ZAKα Recognizes Stalled Ribosomes through Partially Redundant Sensor Domains.

Impairment of ribosome function activates the MAPKKK ZAK, leading to activation of mitogen-activated protein (MAP) kinases p38 and JNK and inflammatory signaling. The mechanistic basis for activation of this ribotoxic stress response (RSR) remains completely obscure. We show that the long isoform of ZAK (ZAKα) directly associates with ribosomes by inserting its flexible C terminus into the ribosomal intersubunit space. Here, ZAKα binds helix 14 of 18S ribosomal RNA (rRNA). An adjacent domain in ZAKα also probes the ribosome, and together, these sensor domains are critically required for RSR activation after inhibition of both the E-site, the peptidyl transferase center (PTC), and ribotoxin action. Finally, we show that ablation of the RSR response leads to organismal phenotypes and decreased lifespan in the nematode Caenorhabditis elegans (C. elegans). Our findings yield mechanistic insight into how cells detect ribotoxic stress and provide experimental in vivo evidence for its physiological importance.

CnoX Is a Chaperedoxin: A Holdase that Protects Its Substrates from Irreversible Oxidation.

Bleach (HOCl) is a powerful oxidant that kills bacteria in part by causing protein aggregation. It inactivates ATP-dependent chaperones, rendering cellular proteins mostly dependent on holdases. Here we identified Escherichia coli CnoX (YbbN) as a folding factor that, when activated by bleach via chlorination, functions as an efficient holdase, protecting the substrates of the major folding systems GroEL/ES and DnaK/J/GrpE. Remarkably, CnoX uniquely combines this function with the ability to prevent the irreversible oxidation of its substrates. This dual activity makes CnoX the founding member of a family of proteins, the "chaperedoxins." Because CnoX displays a thioredoxin fold and a tetratricopeptide (TPR) domain, two structural motifs conserved in all organisms, this investigation sets the stage for the discovery of additional chaperedoxins in bacteria and eukaryotes that could cooperate with proteins from both the Hsp60 and Hsp70 families.

Evolutionary analyses of intrinsically disordered regions reveal widespread signals of conservation.

Intrinsically disordered regions (IDRs) are segments of proteins without stable three-dimensional structures. As this flexibility allows them to interact with diverse binding partners, IDRs play key roles in cell signaling and gene expression. Despite the prevalence and importance of IDRs in eukaryotic proteomes and various biological processes, associating them with specific molecular functions remains a significant challenge due to their high rates of sequence evolution. However, by comparing the observed values of various IDR-associated properties against those generated under a simulated model of evolution, a recent study found most IDRs across the entire yeast proteome contain conserved features. Furthermore, it showed clusters of IDRs with common "evolutionary signatures," i.e. patterns of conserved features, were associated with specific biological functions. To determine if similar patterns of conservation are found in the IDRs of other systems, in this work we applied a series of phylogenetic models to over 7,500 orthologous IDRs identified in the Drosophila genome to dissect the forces driving their evolution. By comparing models of constrained and unconstrained continuous trait evolution using the Brownian motion and Ornstein-Uhlenbeck models, respectively, we identified signals of widespread constraint, indicating conservation of distributed features is mechanism of IDR evolution common to multiple biological systems. In contrast to the previous study in yeast, however, we observed limited evidence of IDR clusters with specific biological functions, which suggests a more complex relationship between evolutionary constraints and function in the IDRs of multicellular organisms.

p38γ is essential for cell cycle progression and liver tumorigenesis.

The cell cycle is a tightly regulated process that is controlled by the conserved cyclin-dependent kinase (CDK)-cyclin protein complex1. However, control of the G0-to-G1 transition is not completely understood. Here we demonstrate that p38 MAPK gamma (p38γ) acts as a CDK-like kinase and thus cooperates with CDKs, regulating entry into the cell cycle. p38γ shares high sequence homology, inhibition sensitivity and substrate specificity with CDK family members. In mouse hepatocytes, p38γ induces proliferation after partial hepatectomy by promoting the phosphorylation of retinoblastoma tumour suppressor protein at known CDK target residues. Lack of p38γ or treatment with the p38γ inhibitor pirfenidone protects against the chemically induced formation of liver tumours. Furthermore, biopsies of human hepatocellular carcinoma show high expression of p38γ, suggesting that p38γ could be a therapeutic target in the treatment of this disease.

Genome manipulation by guide-directed Argonaute cleavage.

Many prokaryotic argonautes (pAgos) mediate DNA interference by using small DNA guides to cleave target DNA. A recent study shows that CbAgo, a pAgo from Clostridium butyricum, induces DNA interference between homologous sequences and generates double-stranded breaks (DSBs) in target DNAs. This mechanism enables the host to defend against invading DNAs such as plasmids and viruses. However, whether such a CbAgo-mediated DNA cleavage is mutagenic remains unexplored. Here we demonstrate that CbAgo, directed by plasmid-encoded guide sequences, can cleave genome target sites and induce chromosome recombination between downstream homologous sequences in Escherichia coli. The recombination rate correlates well with pAgo DNA cleavage activity and the mechanistic study suggests the recombination involves DSBs and RecBCD processing. In RecA-deficient E. coli strain, guide-directed CbAgo cleavage on chromosomes severely impairs cell growth, which can be utilized as counter-selection to assist Lambda-Red recombineering. These findings demonstrate the guide-directed cleavage of pAgo on the host genome is mutagenic and can lead to different outcomes according to the function of the host DNA repair machinery. We anticipate this novel DNA-guided interference to be useful in broader genetic manipulation. Our study also provides an in vivo assay to characterize or engineer pAgo DNA cleavage activity.

The SIAH E3 ubiquitin ligases promote Wnt/ß-catenin signaling through mediating Wnt-induced Axin degradation.

The Wnt/ß-catenin signaling pathway plays essential roles in embryonic development and adult tissue homeostasis. Axin is a concentration-limiting factor responsible for the formation of the ß-catenin destruction complex. Wnt signaling itself promotes the degradation of Axin. However, the underlying molecular mechanism and biological relevance of this targeting of Axin have not been elucidated. Here, we identify SIAH1/2 (SIAH) as the E3 ligase mediating Wnt-induced Axin degradation. SIAH proteins promote the ubiquitination and proteasomal degradation of Axin through interacting with a VxP motif in the GSK3-binding domain of Axin, and this function of SIAH is counteracted by GSK3 binding to Axin. Structural analysis reveals that the Axin segment responsible for SIAH binding is also involved in GSK3 binding but adopts distinct conformations in Axin/SIAH and Axin/GSK3 complexes. Knockout of SIAH1 blocks Wnt-induced Axin ubiquitination and attenuates Wnt-induced ß-catenin stabilization. Our data suggest that Wnt-induced dissociation of the Axin/GSK3 complex allows SIAH to interact with Axin not associated with GSK3 and promote its degradation and that SIAH-mediated Axin degradation represents an important feed-forward mechanism to achieve sustained Wnt/ß-catenin signaling.

An FPGA-based hardware accelerator supporting sensitive sequence homology filtering with profile hidden Markov models.

BACKGROUND: Sequence alignment lies at the heart of genome sequence annotation. While the BLAST suite of alignment tools has long held an important role in alignment-based sequence database search, greater sensitivity is achieved through the use of profile hidden Markov models (pHMMs). Here, we describe an FPGA hardware accelerator, called HAVAC, that targets a key bottleneck step (SSV) in the analysis pipeline of the popular pHMM alignment tool, HMMER. RESULTS: The HAVAC kernel calculates the SSV matrix at 1739 GCUPS on a ∼  $3000 Xilinx Alveo U50 FPGA accelerator card, ∼  227× faster than the optimized SSV implementation in nhmmer. Accounting for PCI-e data transfer data processing, HAVAC is 65× faster than nhmmer's SSV with one thread and 35× faster than nhmmer with four threads, and uses ∼  31% the energy of a traditional high end Intel CPU. CONCLUSIONS: HAVAC demonstrates the potential offered by FPGA hardware accelerators to produce dramatic speed gains in sequence annotation and related bioinformatics applications. Because these computations are performed on a co-processor, the host CPU remains free to simultaneously compute other aspects of the analysis pipeline.

The Discovery of Chicken Foxp3 Demands Redefinition of Avian Regulatory T Cells.

Since the publication of the first chicken genome sequence, we have encountered genes playing key roles in mammalian immunology, but being seemingly absent in birds. One of those was, until recently, Foxp3, the master transcription factor of regulatory T cells in mammals. Therefore, avian regulatory T cell research is still poorly standardized. In this study we identify a chicken ortholog of Foxp3 We prove sequence homology with known mammalian and sauropsid sequences, but also reveal differences in major domains. Expression profiling shows an association of Foxp3 and CD25 expression levels in CD4+CD25+ peripheral T cells and identifies a CD4-CD25+Foxp3high subset of thymic lymphocytes that likely represents yet undescribed avian regulatory T precursor cells. We conclude that Foxp3 is existent in chickens and that it shares certain functional characteristics with its mammalian ortholog. Nevertheless, pathways for regulatory T cell development and Foxp3 function are likely to differ between mammals and birds. The identification and characterization of chicken Foxp3 will help to define avian regulatory T cells and to analyze their functional properties and thereby advance the field of avian immunology.

High frequency of an otherwise rare phenotype in a small and isolated tiger population.

Most endangered species exist today in small populations, many of which are isolated. Evolution in such populations is largely governed by genetic drift. Empirical evidence for drift affecting striking phenotypes based on substantial genetic data are rare. Approximately 37% of tigers (Panthera tigris) in the Similipal Tiger Reserve (in eastern India) are pseudomelanistic, characterized by wide, merged stripes. Camera trap data across the tiger range revealed the presence of pseudomelanistic tigers only in Similipal. We investigated the genetic basis for pseudomelanism and examined the role of drift in driving this phenotype's frequency. Whole-genome data and pedigree-based association analyses from captive tigers revealed that pseudomelanism cosegregates with a conserved and functionally important coding alteration in Transmembrane Aminopeptidase Q (Taqpep), a gene responsible for similar traits in other felid species. Noninvasive sampling of tigers revealed a high frequency of the Taqpep p.H454Y mutation in Similipal (12 individuals, allele frequency = 0.58) and absence from all other tiger populations (395 individuals). Population genetic analyses confirmed few (minimal number) tigers in Similipal, and its genetic isolation, with poor geneflow. Pairwise FST (0.33) at the mutation site was high but not an outlier. Similipal tigers had low diversity at 81 single nucleotide polymorphisms (mean heterozygosity = 0.28, SD = 0.27). Simulations were consistent with founding events and drift as possible drivers for the observed stark difference of allele frequency. Our results highlight the role of stochastic processes in the evolution of rare phenotypes. We highlight an unusual evolutionary trajectory in a small and isolated population of an endangered species.

Decoding the link of microbiome niches with homologous sequences enables accurately targeted protein structure prediction.

Information derived from metagenome sequences through deep-learning techniques has significantly improved the accuracy of template free protein structure modeling. However, most of the deep learning-based modeling studies are based on blind sequence database searches and suffer from low efficiency in computational resource utilization and model construction, especially when the sequence library becomes prohibitively large. We proposed a MetaSource model built on 4.25 billion microbiome sequences from four major biomes (Gut, Lake, Soil, and Fermentor) to decode the inherent linkage of microbial niches with protein homologous families. Large-scale protein family folding experiments on 8,700 unknown Pfam families showed that a microbiome targeted approach with multiple sequence alignment constructed from individual MetaSource biomes requires more than threefold less computer memory and CPU (central processing unit) time but generates contact-map and three-dimensional structure models with a significantly higher accuracy, compared with that using combined metagenome datasets. These results demonstrate an avenue to bridge the gap between the rapidly increasing metagenome databases and the limited computing resources for efficient genome-wide database mining, which provides a useful bluebook to guide future microbiome sequence database and modeling development for high-accuracy protein structure and function prediction.

Highly conserved and cis-acting lncRNAs produced from paralogous regions in the center of HOXA and HOXB clusters in the endoderm lineage.

Long noncoding RNAs (lncRNAs) have been shown to play important roles in gene regulatory networks acting in early development. There has been rapid turnover of lncRNA loci during vertebrate evolution, with few human lncRNAs conserved beyond mammals. The sequences of these rare deeply conserved lncRNAs are typically not similar to each other. Here, we characterize HOXA-AS3 and HOXB-AS3, lncRNAs produced from the central regions of the HOXA and HOXB clusters. Sequence-similar orthologs of both lncRNAs are found in multiple vertebrate species and there is evident sequence similarity between their promoters, suggesting that the production of these lncRNAs predates the duplication of the HOX clusters at the root of the vertebrate lineage. This conservation extends to similar expression patterns of the two lncRNAs, in particular in cells transiently arising during early development or in the adult colon. Functionally, the RNA products of HOXA-AS3 and HOXB-AS3 regulate the expression of their overlapping HOX5-7 genes both in HT-29 cells and during differentiation of human embryonic stem cells. Beyond production of paralogous protein-coding and microRNA genes, the regulatory program in the HOX clusters therefore also relies on paralogous lncRNAs acting in restricted spatial and temporal windows of embryonic development and cell differentiation.

Bacterial translation machinery for deliberate mistranslation of the genetic code.

Inaccurate expression of the genetic code, also known as mistranslation, is an emerging paradigm in microbial studies. Growing evidence suggests that many microbial pathogens can deliberately mistranslate their genetic code to help invade a host or evade host immune responses. However, discovering different capacities for deliberate mistranslation remains a challenge because each group of pathogens typically employs a unique mistranslation mechanism. In this study, we address this problem by studying duplicated genes of aminoacyl-transfer RNA (tRNA) synthetases. Using bacterial prolyl-tRNA synthetase (ProRS) genes as an example, we identify an anomalous ProRS isoform, ProRSx, and a corresponding tRNA, tRNAProA, that are predominately found in plant pathogens from Streptomyces species. We then show that tRNAProA has an unusual hybrid structure that allows this tRNA to mistranslate alanine codons as proline. Finally, we provide biochemical, genetic, and mass spectrometric evidence that cells which express ProRSx and tRNAProA can translate GCU alanine codons as both alanine and proline. This dual use of alanine codons creates a hidden proteome diversity due to stochastic Ala→Pro mutations in protein sequences. Thus, we show that important plant pathogens are equipped with a tool to alter the identity of their sense codons. This finding reveals the initial example of a natural tRNA synthetase/tRNA pair for dedicated mistranslation of sense codons.

HyPep: An Open-Source Software for Identification and Discovery of Neuropeptides Using Sequence Homology Search.

Neuropeptides are a class of endogenous peptides that have key regulatory roles in biochemical, physiological, and behavioral processes. Mass spectrometry analyses of neuropeptides often rely on protein informatics tools for database searching and peptide identification. As neuropeptide databases are typically experimentally built and comprised of short sequences with high sequence similarity to each other, we developed a novel database searching tool, HyPep, which utilizes sequence homology searching for peptide identification. HyPep aligns de novo sequenced peptides, generated through PEAKS software, with neuropeptide database sequences and identifies neuropeptides based on the alignment score. HyPep performance was optimized using LC-MS/MS measurements of peptide extracts from various Callinectes sapidus neuronal tissue types and compared with a commercial database searching software, PEAKS DB. HyPep identified more neuropeptides from each tissue type than PEAKS DB at 1% false discovery rate, and the false match rate from both programs was 2%. In addition to identification, this report describes how HyPep can aid in the discovery of novel neuropeptides.

Prediction of NB-LRR resistance genes based on full-length sequence homology.

The activation of plant immunity is mediated by resistance (R)-gene receptors, also known as nucleotide-binding leucine-rich repeat (NB-LRR) genes, which in turn trigger the authentic defense response. R-gene identification is a crucial goal for both classic and modern plant breeding strategies for disease resistance. The conventional method identifies NB-LRR genes using a protein motif/domain-based search (PDS) within an automatically predicted gene set of the respective genome assembly. PDS proved to be imprecise since repeat masking prior to automatic genome annotation unwittingly prevented comprehensive NB-LRR gene detection. Furthermore, R-genes have diversified in a species-specific manner, so that NB-LRR gene identification cannot be universally standardized. Here, we present the full-length Homology-based R-gene Prediction (HRP) method for the comprehensive identification and annotation of a genome's R-gene repertoire. Our method has substantially addressed the complex genomic organization of tomato (Solanum lycopersicum) NB-LRR gene loci, proving to be more performant than the well-established RenSeq approach. HRP efficiency was also tested on three differently assembled and annotated Beta sp. genomes. Indeed, HRP identified up to 45% more full-length NB-LRR genes compared to previous approaches. HRP also turned out to be a more refined strategy for R-gene allele mining, testified by the identification of hitherto undiscovered Fom-2 homologs in five Cucurbita sp. genomes. In summary, our high-performance method for full-length NB-LRR gene discovery will propel the identification of novel R-genes towards development of improved cultivars.

CRISPR/Cas9 searches for a protospacer adjacent motif by lateral diffusion.

The Streptococcus pyogenes CRISPR/Cas9 (SpCas9) nuclease has been widely applied in genetic engineering. Despite its importance in genome editing, aspects of the precise molecular mechanism of Cas9 activity remain ambiguous. In particular, because of the lack of a method with high spatio-temporal resolution, transient interactions between Cas9 and DNA could not be reliably investigated. It therefore remains controversial how Cas9 searches for protospacer adjacent motif (PAM) sequences. We have developed single-molecule Förster resonance energy transfer (smFRET) assays to monitor transient interactions of Cas9 and DNA in real time. Our study shows that Cas9 interacts with the PAM sequence weakly, yet probing neighboring sequences via facilitated diffusion. This dynamic mode of interactions leads to translocation of Cas9 to another PAM nearby and consequently an on-target sequence. We propose a model in which lateral diffusion competes with three-dimensional diffusion and thus is involved in PAM finding and consequently on-target binding. Our results imply that the neighboring sequences can be very important when choosing a target in genetic engineering applications.

Structures of the heart specific SERCA2a Ca2+-ATPase.

The sarcoplasmic/endoplasmic reticulum Ca2+-ATPase 2a (SERCA2a) performs active reuptake of cytoplasmic Ca2+ and is a major regulator of cardiac muscle contractility. Dysfunction or dysregulation of SERCA2a is associated with heart failure, while restoring its function is considered as a therapeutic strategy to restore cardiac performance. However, its structure has not yet been determined. Based on native, active protein purified from pig ventricular muscle, we present the first crystal structures of SERCA2a, determined in the CPA-stabilized E2-AlF4- form (3.3 Å) and the Ca2+-occluded [Ca2]E1-AMPPCP form (4.0 Å). The structures are similar to the skeletal muscle isoform SERCA1a pointing to a conserved mechanism. We seek to explain the kinetic differences between SERCA1a and SERCA2a. We find that several isoform-specific residues are acceptor sites for post-translational modifications. In addition, molecular dynamics simulations predict that isoform-specific residues support distinct intramolecular interactions in SERCA2a and SERCA1a. Our experimental observations further indicate that isoform-specific intramolecular interactions are functionally relevant, and may explain the kinetic differences between SERCA2a and SERCA1a.

Ssu72 phosphatase is a conserved telomere replication terminator.

Telomeres, the protective ends of eukaryotic chromosomes, are replicated through concerted actions of conventional DNA polymerases and elongated by telomerase, but the regulation of this process is not fully understood. Telomere replication requires (Ctc1/Cdc13)-Stn1-Ten1, a telomeric ssDNA-binding complex homologous to RPA Here, we show that the evolutionarily conserved phosphatase Ssu72 is responsible for terminating the cycle of telomere replication in fission yeast. Ssu72 controls the recruitment of Stn1 to telomeres by regulating Stn1 phosphorylation at Ser74, a residue located within its conserved OB-fold domain. Consequently, ssu72∆ mutants are defective in telomere replication and exhibit long 3'-ssDNA overhangs, indicative of defective lagging-strand DNA synthesis. We also show that hSSU72 regulates telomerase activation in human cells by controlling recruitment of hSTN1 to telomeres. These results reveal a previously unknown yet conserved role for the phosphatase SSU72, whereby this enzyme controls telomere homeostasis by activating lagging-strand DNA synthesis, thus terminating the cycle of telomere replication.

Structures of neurexophilin-neurexin complexes reveal a regulatory mechanism of alternative splicing.

Neurexins are presynaptic, cell-adhesion molecules that specify the functional properties of synapses via interactions with trans-synaptic ligands. Neurexins are extensively alternatively spliced at six canonical sites that regulate multifarious ligand interactions, but the structural mechanisms underlying alternative splicing-dependent neurexin regulation are largely unknown. Here, we determined high-resolution structures of the complex of neurexophilin-1 and the second laminin/neurexin/sex-hormone-binding globulin domain (LNS2) of neurexin-1 and examined how alternative splicing at splice site #2 (SS2) regulates the complex. Our data reveal a unique, extensive, neurexophilin-neurexin binding interface that extends the jelly-roll ß-sandwich of LNS2 of neurexin-1 into neurexophilin-1. The SS2A insert of LNS2 augments this interface, increasing the binding affinity of LNS2 for neurexophilin-1. Taken together, our data reveal an unexpected architecture of neurexophilin-neurexin complexes that accounts for the modulation of binding by alternative splicing, which in turn regulates the competition of neurexophilin for neurexin binding with other ligands.