The history and conceptual framework of assays and screens.

Since the 16th century, assays and screens have been essential for scientific investigation. However, most methods could be significantly improved, especially in accuracy, scalability, and often lack adequate comparisons to negative controls. There is a lack of consistency in distinguishing assays, in which accuracy is the main goal, from screens, in which scalability is prioritized over accuracy. We dissected and modernized the original definitions of assays and screens based upon recent developments and the conceptual framework of the original definitions. All methods have three components: design/measurement, performance, and interpretation. We propose a model of method development in which reproducible observations become new methods, initially assessed by sensitivity. Further development can proceed along a path to either screens or assays. The screen path focuses on scalability first, but can later prioritize analysis of negatives. Alternatively, the assay path first compares results to negative controls, assessing specificity and accuracy, later adding scalability. Both pathways converge on a high-accuracy and throughput (HAT) assay, like next generation sequencing, which we suggest should be the ultimate goal of all testing methods. Our model will help scientists better select among available methods, as well as improve existing methods, expanding their impact on science.

Most synonymous allelic variants in HIV tat are not silent.

The genetic code has degenerate codons that produce no change in the translated protein sequence and are generally thought to be silent. However, some synonymous variants are clearly not silent. Herein, we questioned the frequency of non-silent synonymous variants. We tested how random synonymous variants in the HIV Tat transcription factor effect transcription of an LTR-GFP reporter. Our model system has the advantage of directly measuring the function of the gene in human cells. Approximately, 67% of synonymous variants in Tat were non-silent, either having reduced activity or were full loss-of-function alleles. Eight mutant codons had higher codon usage than wild type, accompanied by reduced transcriptional activity. These were clustered on a loop in the Tat structure. We conclude that most synonymous Tat variants are not silent in human cells, and 25% are associated with changes in codon usage, likely effecting protein folding.

GigaAssay - An adaptable high-throughput saturation mutagenesis assay platform.

High-throughput assay systems have had a large impact on understanding the mechanisms of basic cell functions. However, high-throughput assays that directly assess molecular functions are limited. Herein, we describe the "GigaAssay", a modular high-throughput one-pot assay system for measuring molecular functions of thousands of genetic variants at once. In this system, each cell was infected with one virus from a library encoding thousands of Tat mutant proteins, with each viral particle encoding a random unique molecular identifier (UMI). We demonstrate proof of concept by measuring transcription of a GFP reporter in an engineered reporter cell line driven by binding of the HIV Tat transcription factor to the HIV long terminal repeat. Infected cells were flow-sorted into 3 bins based on their GFP fluorescence readout. The transcriptional activity of each Tat mutant was calculated from the ratio of signals from each bin. The use of UMIs in the GigaAssay produced a high average accuracy (95%) and positive predictive value (98%) determined by comparison to literature benchmark data, known C-terminal truncations, and blinded independent mutant tests. Including the substitution tolerance with structure/function analysis shows restricted substitution types spatially concentrated in the Cys-rich region. Tat has abundant intragenic epistasis (10%) when single and double mutants are compared.

Synonymous Variants of Uncertain Silence.

Synonymous variants, traditionally regarded as silent mutations due to their lack of impact on protein sequence, structure and function, have been the subject of increasing scrutiny. This commentary explores the emerging evidence challenging the notion of synonymous variants as functionally inert. Analysis of the activity of 70 synonymous variants in the HIV Tat transcription factor revealed that 50% of the variants exhibited significant deviations from wild-type activity. Our analysis supports previous work and raises important questions about the broader impact of non-silent synonymous variants in human genes. Considering the potential functional implications, the authors propose classifying such variants as "synonymous variants of uncertain silence" (sVUS), highlighting the need for cautious interpretation and further investigations in clinical and genetic testing settings.

Accurate Prediction of Transcriptional Activity of Single Missense Variants in HIV Tat with Deep Learning.

Tat is an essential gene for increasing the transcription of all HIV genes, and affects HIV replication, HIV exit from latency, and AIDS progression. The Tat gene frequently mutates in vivo and produces variants with diverse activities, contributing to HIV viral heterogeneity as well as drug-resistant clones. Thus, identifying the transcriptional activities of Tat variants will help to better understand AIDS pathology and treatment. We recently reported the missense mutation landscape of all single amino acid Tat variants. In these experiments, a fraction of double missense alleles exhibited intragenic epistasis. However, it is too time-consuming and costly to determine the effect of the variants for all double mutant alleles through experiments. Therefore, we propose a combined GigaAssay/deep learning approach. As a first step to determine activity landscapes for complex variants, we evaluated a deep learning framework using previously reported GigaAssay experiments to predict how transcription activity is affected by Tat variants with single missense substitutions. Our approach achieved a 0.94 Pearson correlation coefficient when comparing the predicted to experimental activities. This hybrid approach can be extensible to more complex Tat alleles for a better understanding of the genetic control of HIV genome transcription.

XRCC4 and MRE11 Roles and Transcriptional Response to Repair of TALEN-Induced Double-Strand DNA Breaks.

Double-strand breaks (DSB) are one of the most lethal forms of DNA damage that, if left unrepaired, can lead to genomic instability, cellular transformation, and cell death. In this work, we examined how repair of transcription activator-like effector nuclease (TALEN)-induced DNA damage was altered when knocking out, or inhibiting a function of, two DNA repair proteins, XRCC4 and MRE11, respectively. We developed a fluorescent reporter assay that uses TALENs to introduce DSB and detected repair by the presence of GFP fluorescence. We observed repair of TALEN-induced breaks in the XRCC4 knockout cells treated with mirin (a pharmacological inhibitor of MRE11 exonuclease activity), albeit with ~40% reduced efficiency compared to normal cells. Editing in the absence of XRCC4 or MRE11 exonuclease was robust, with little difference between the indel profiles amongst any of the groups. Reviewing the transcriptional profiles of the mirin-treated XRCC4 knockout cells showed 307 uniquely differentially expressed genes, a number far greater than for either of the other cell lines (the HeLa XRCC4 knockout sample had 83 genes, and the mirin-treated HeLa cells had 30 genes uniquely differentially expressed). Pathways unique to the XRCC4 knockout+mirin group included differential expression of p53 downstream pathways, and metabolic pathways indicating cell adaptation for energy regulation and stress response. In conclusion, our study showed that TALEN-induced DSBs are repaired, even when a key DSB repair protein or protein function is not operational, without a change in indel profiles. However, transcriptional profiles indicate the induction of unique cellular responses dependent upon the DNA repair protein(s) hampered.

The carboxy-terminus, a key regulator of protein function.

All proteins end with a carboxyl terminus that has unique biophysical properties and is often disordered. Although there are examples of important C-termini functions, a more global role for the C-terminus is not yet established. In this review, we summarize research on C-termini, a unique region in proteins that cells exploit. Alternative splicing and proteolysis increase the diversity of proteins and peptides in cells with unique C-termini. The C-termini of proteins contain minimotifs, short peptides with an encoded function generally characterized as binding, posttranslational modifications, and trafficking. Many of these activities are specific to minimotifs on the C-terminus. Approximately 13% of C-termini in the human proteome have a known minimotif, and the majority, if not all of the remaining termini have conserved motifs inferring a function that remains to be discovered. C-termini, their predictions, and their functions are collated in the C-terminome, Proteus, and Terminus Oriented Protein Function INferred Database (TopFIND) database/web systems. Many C-termini are well conserved, and some have a known role in health and disease. We envision that this summary of C-termini will guide future investigation of their biochemical and physiological significance.

Can Designer Indels Be Tailored by Gene Editing?: Can Indels Be Customized?

Genome editing with engineered nucleases (GEENs) introduce site-specific DNA double-strand breaks (DSBs) and repairs DSBs via nonhomologous end-joining (NHEJ) pathways that eventually create indels (insertions/deletions) in a genome. Whether the features of indels resulting from gene editing could be customized is asked. A review of the literature reveals how gene editing technologies via NHEJ pathways impact gene editing. The survey consolidates a body of literature that suggests that the type (insertion, deletion, and complex) and the approximate length of indel edits can be somewhat customized with different GEENs and by manipulating the expression of key NHEJ genes. Structural data suggest that binding of GEENs to DNA may interfere with binding of key components of DNA repair complexes, favoring either classical- or alternative-NHEJ. The hypotheses have some limitations, but if validated, will enable scientists to better control indel makeup, holding promise for basic science and clinical applications of gene editing. Also see the video abstract here https://youtu.be/vTkJtUsLi3w.

Minimotif Miner 4: a million peptide minimotifs and counting.

Minimotif Miner (MnM) is a database and web system for analyzing short functional peptide motifs, termed minimotifs. We present an update to MnM growing the database from ∼300 000 to >1 000 000 minimotif consensus sequences and instances. This growth comes largely from updating data from existing databases and annotation of articles with high-throughput approaches analyzing different types of post-translational modifications. Another update is mapping human proteins and their minimotifs to know human variants from the dbSNP, build 150. Now MnM 4 can be used to generate mechanistic hypotheses about how human genetic variation affect minimotifs and outcomes. One example of the utility of the combined minimotif/SNP tool identifies a loss of function missense SNP in a ubiquitylation minimotif encoded in the excision repair cross-complementing 2 (ERCC2) nucleotide excision repair gene. This SNP reaches genome wide significance for many types of cancer and the variant identified with MnM 4 reveals a more detailed mechanistic hypothesis concerning the role of ERCC2 in cancer. Other updates to the web system include a new architecture with migration of the web system and database to Docker containers for better performance and management. Weblinks:minimotifminer.org and mnm.engr.uconn.edu.

Prioritization of Variants for Investigation of Genotype-Directed Nutrition in Human Superpopulations.

Dietary guidelines recommended by key health agencies are generally designed for a global population. However, ethnicity affects human disease and environment-gene interactions, including nutrient intake. Historically, isolated human populations with different genetic backgrounds have adapted to distinct environments with varying food sources. Ethnicity is relevant to the interaction of food intake with genes and disease susceptibility; yet major health agencies generally do not recommend food and nutrients codified by population genotypes and their frequencies. In this paper, we have consolidated published nutrigenetic variants and examine their frequencies in human superpopulations to prioritize these variants for future investigation of population-specific genotype-directed nutrition. The nutrients consumed by individuals interact with their genome and may alter disease risk. Herein, we searched the literature, designed a data model, and manually curated hundreds of papers. The resulting database houses 101 variants that reached significance (p < 0.05), from 35 population studies. Nutrigenetic variants associated with modified nutrient intake have the potential to reduce the risk of colorectal cancer, obesity, metabolic syndrome, type 2 diabetes, and several other diseases. Since many nutrigenetic studies have identified a major variant in some populations, we suggest that superpopulation-specific genotype-directed nutrition modifications be prioritized for future study and evaluation. Genotype-directed nutrition approaches to dietary modification have the potential to reduce disease risk in select human populations.

Natural variability of minimotifs in 1092 people indicates that minimotifs are targets of evolution.

Since the function of a short contiguous peptide minimotif can be introduced or eliminated by a single point mutation, these functional elements may be a source of human variation and a target of selection. We analyzed the variability of ∼300 000 minimotifs in 1092 human genomes from the 1000 Genomes Project. Most minimotifs have been purified by selection, with a 94% invariance, which supports important functional roles for minimotifs. Minimotifs are generally under negative selection, possessing high genomic evolutionary rate profiling (GERP) and sitewise likelihood-ratio (SLR) scores. Some are subject to neutral drift or positive selection, similar to coding regions. Most SNPs in minimotif were common variants, but with minor allele frequencies generally <10%. This was supported by low substation rates and few newly derived minimotifs. Several minimotif alleles showed different intercontinental and regional geographic distributions, strongly suggesting a role for minimotifs in adaptive evolution. We also note that 4% of PTM minimotif sites in histone tails were common variants, which has the potential to differentially affect DNA packaging among individuals. In conclusion, minimotifs are a source of functional genetic variation in the human population; thus, they are likely to be an important target of selection and evolution.

TALEN gene editing takes aim on HIV.

Transcription activator-like effector nucleases (TALENs) are one of several types of programmable, engineered nucleases that bind and cleave specific DNA sequences. Cellular machinery repairs the cleaved DNA by introducing indels. In this review, we emphasize the potential, explore progress, and identify challenges in using TALENs as a therapeutic tool to treat HIV infection. TALENs have less off-target editing and can be more effective at tolerating HIV escape mutations than CRISPR/Cas-9. Scientists have explored TALEN-mediated editing of host genes such as viral entry receptors (CCR5 and CXCR4) and a protein involved in proviral integration (LEDGF/p75). Viral targets include the proviral DNA, particularly focused on the long terminal repeats. Major challenges with translating gene therapy from bench to bedside are improving cleavage efficiency and delivery, while minimizing off-target editing, cytotoxicity, and immunogenicity. However, rapid improvements in TALEN technology are enhancing cleavage efficiency and specificity. Therapeutic testing in animal models of HIV infection will help determine whether TALENs are a viable HIV treatment therapy. TALENs or other engineered nucleases could shift the therapeutic paradigm from life-long antiretroviral therapy toward eradication of HIV infection.

The human phosphotyrosine signaling network: evolution and hotspots of hijacking in cancer.

Phosphotyrosine (pTyr) signaling, which plays a central role in cell-cell and cell-environment interactions, has been considered to be an evolutionary innovation in multicellular metazoans. However, neither the emergence nor the evolution of the human pTyr signaling system is currently understood. Tyrosine kinase (TK) circuits, each of which consists of a TK writer, a kinase substrate, and a related reader, such as Src homology (SH) 2 domains and pTyr-binding (PTB) domains, comprise the core machinery of the pTyr signaling network. In this study, we analyzed the evolutionary trajectories of 583 literature-derived and 50,000 computationally predicted human TK circuits in 19 representative eukaryotic species and assigned their evolutionary origins. We found that human TK circuits for intracellular pTyr signaling originated largely from primitive organisms, whereas the inter- or extracellular signaling circuits experienced significant expansion in the bilaterian lineage through the "back-wiring" of newly evolved kinases to primitive substrates and SH2/PTB domains. Conversely, the TK circuits that are involved in tissue-specific signaling evolved mainly in vertebrates by the back-wiring of vertebrate substrates to primitive kinases and SH2/PTB domains. Importantly, we found that cancer signaling preferentially employs the pTyr sites, which are linked to more TK circuits. Our work provides insights into the evolutionary paths of the human pTyr signaling circuits and suggests the use of a network approach for cancer intervention through the targeting of key pTyr sites and their associated signaling hubs in the network.

Loss-of-function mutations of ILDR1 cause autosomal-recessive hearing impairment DFNB42.

By using homozygosity mapping in a consanguineous Pakistani family, we detected linkage of nonsyndromic hearing loss to a 7.6 Mb region on chromosome 3q13.31-q21.1 within the previously reported DFNB42 locus. Subsequent candidate gene sequencing identified a homozygous nonsense mutation (c.1135G>T [p.Glu379X]) in ILDR1 as the cause of hearing impairment. By analyzing additional consanguineous families with homozygosity at this locus, we detected ILDR1 mutations in the affected individuals of 10 more families from Pakistan and Iran. The identified ILDR1 variants include missense, nonsense, frameshift, and splice-site mutations as well as a start codon mutation in the family that originally defined the DFNB42 locus. ILDR1 encodes the evolutionarily conserved immunoglobulin-like domain containing receptor 1, a putative transmembrane receptor of unknown function. In situ hybridization detected expression of Ildr1, the murine ortholog, early in development in the vestibule and in hair cells and supporting cells of the cochlea. Expression in hair cell- and supporting cell-containing neurosensory organs is conserved in the zebrafish, in which the ildr1 ortholog is prominently expressed in the developing ear and neuromasts of the lateral line. These data identify loss-of-function mutations of ILDR1, a gene with a conserved expression pattern pointing to a conserved function in hearing in vertebrates, as underlying nonsyndromic prelingual sensorineural hearing impairment.

Minimotif Miner 3.0: database expansion and significantly improved reduction of false-positive predictions from consensus sequences.

Minimotif Miner (MnM available at http://minimotifminer.org or http://mnm.engr.uconn.edu) is an online database for identifying new minimotifs in protein queries. Minimotifs are short contiguous peptide sequences that have a known function in at least one protein. Here we report the third release of the MnM database which has now grown 60-fold to approximately 300,000 minimotifs. Since short minimotifs are by their nature not very complex we also summarize a new set of false-positive filters and linear regression scoring that vastly enhance minimotif prediction accuracy on a test data set. This online database can be used to predict new functions in proteins and causes of disease.

Systematic Assessment of Protein C-Termini Mutated in Human Disorders.

All proteins have a carboxyl terminus, and we previously summarized eight mutations in binding and trafficking sequence determinants in the C-terminus that, when disrupted, cause human diseases. These sequence elements for binding and trafficking sites, as well as post-translational modifications (PTMs), are called minimotifs or short linear motifs. We wanted to determine how frequently mutations in minimotifs in the C-terminus cause disease. We searched specifically for PTMs because mutation of a modified amino acid almost always changes the chemistry of the side chain and can be interpreted as loss-of-function. We analyzed data from ClinVar for disease variants, Minimotif Miner and the C-terminome for PTMs, and RefSeq for protein sequences, yielding 20 such potential disease-causing variants. After additional screening, they include six with a previously reported PTM disruption mechanism and nine with new hypotheses for mutated minimotifs in C-termini that may cause disease. These mutations were generally for different genes, with four different PTM types and several different diseases. Our study helps to identify new molecular mechanisms for nine separate variants that cause disease, and this type of analysis could be extended as databases grow and to binding and trafficking motifs. We conclude that mutated motifs in C-termini are an infrequent cause of disease.

Single-cell RNA sequencing reveals in vivo osteoimmunology interactions between the immune and skeletal systems.

Background: While osteoimmunology interactions between the immune and skeletal systems are known to play an important role in osteoblast development, differentiation and bone metabolism related disease like osteoporosis, such interactions in either bone microenvironment or peripheral circulation in vivo at the single-cell resolution have not yet been characterized. Methods: We explored the osteoimmunology communications between immune cells and osteoblastic lineage cells (OBCs) by performing CellphoneDB and CellChat analyses with single-cell RNA sequencing (scRNA-seq) data from human femoral head. We also explored the osteoimmunology effects of immune cells in peripheral circulation on skeletal phenotypes. We used a scRNA-seq dataset of peripheral blood monocytes (PBMs) to perform deconvolution analysis. Then weighted gene co-expression network analysis (WGCNA) was used to identify monocyte subtype-specific subnetworks. We next used cell-specific network (CSN) and the least absolute shrinkage and selection operator (LASSO) to analyze the correlation of a gene subnetwork identified by WGCNA with bone mineral density (BMD). Results: We constructed immune cell and OBC communication networks and further identified L-R genes, such as JAG1 and NOTCH1/2, with ossification related functions. We also found a Mono4 related subnetwork that may relate to BMD variation in both older males and postmenopausal female subjects. Conclusions: This is the first study to identify numerous ligand-receptor pairs that likely mediate signals between immune cells and osteoblastic lineage cells. This establishes a foundation to reveal advanced and in-depth osteoimmunology interactions to better understand the relationship between local bone microenvironment and immune cells in peripheral blood and the impact on bone phenotypes.

Data supporting a saturation mutagenesis assay for Tat-driven transcription with the GigaAssay.

The data in this article are associated with the research paper "GigaAssay - an adaptable high-throughput saturation mutagenesis assay" [1]. The raw data are sequence reads of HIV-1 Tat cDNA amplified from cellular genomic DNA in a new single-pot saturation mutagenesis assay designated the "GigaAssay". A bioinformatic pipeline and parameters used to analyze the data. Raw, processed, analyzed, and filtered data are reported. The data is processed to calculate the Tat-driven transcription activity for cells with each possible single amino acid substitution in Tat. This data can be reused to interpret Tat intermolecular interactions and HIV latency. This is one of the largest and most complete datasets regarding the impact of amino acid substitutions within a single protein on a molecular function.

A computational tool for identifying minimotifs in protein-protein interactions and improving the accuracy of minimotif predictions.

Protein-protein interactions are important to understanding cell functions; however, our theoretical understanding is limited. There is a general discontinuity between the well-accepted physical and chemical forces that drive protein-protein interactions and the large collections of identified protein-protein interactions in various databases. Minimotifs are short functional peptide sequences that provide a basis to bridge this gap in knowledge. However, there is no systematic way to study minimotifs in the context of protein-protein interactions or vice versa. Here we have engineered a set of algorithms that can be used to identify minimotifs in known protein-protein interactions and implemented this for use by scientists in Minimotif Miner. By globally testing these algorithms on verified data and on 100 individual proteins as test cases, we demonstrate the utility of these new computation tools. This tool also can be used to reduce false-positive predictions in the discovery of novel minimotifs. The statistical significance of these algorithms is demonstrated by an ROC analysis (P = 0.001).

VENN, a tool for titrating sequence conservation onto protein structures.

Residue conservation is an important, established method for inferring protein function, modularity and specificity. It is important to recognize that it is the 3D spatial orientation of residues that drives sequence conservation. Considering this, we have built a new computational tool, VENN that allows researchers to interactively and graphically titrate sequence homology onto surface representations of protein structures. Our proposed titration strategies reveal critical details that are not readily identified using other existing tools. Analyses of a bZIP transcription factor and receptor recognition of Fibroblast Growth Factor using VENN revealed key specificity determinants. Weblink: http://sbtools.uchc.edu/venn/.