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.

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.

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.

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.

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.

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.

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.

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.

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.

Single-cell RNA sequencing deconvolutes the in vivo heterogeneity of human bone marrow-derived mesenchymal stem cells.

Bone marrow-derived mesenchymal stem cells (BM-MSCs) are multipotent stromal cells that have a critical role in the maintenance of skeletal tissues such as bone, cartilage, and the fat in bone marrow. In addition to providing microenvironmental support for hematopoietic processes, BM-MSCs can differentiate into various mesodermal lineages including osteoblast/osteocyte, chondrocyte, and adipocyte that are crucial for bone metabolism. While BM-MSCs have high cell-to-cell heterogeneity in gene expression, the cell subtypes that contribute to this heterogeneity in vivo in humans have not been characterized. To investigate the transcriptional diversity of BM-MSCs, we applied single-cell RNA sequencing (scRNA-seq) on freshly isolated CD271+ BM-derived mononuclear cells (BM-MNCs) from two human subjects. We successfully identified LEPRhiCD45low BM-MSCs within the CD271+ BM-MNC population, and further codified the BM-MSCs into distinct subpopulations corresponding to the osteogenic, chondrogenic, and adipogenic differentiation trajectories, as well as terminal-stage quiescent cells. Biological functional annotations of the transcriptomes suggest that osteoblast precursors induce angiogenesis coupled with osteogenesis, and chondrocyte precursors have the potential to differentiate into myocytes. We also discovered transcripts for several clusters of differentiation (CD) markers that were either highly expressed (e.g., CD167b, CD91, CD130 and CD118) or absent (e.g., CD74, CD217, CD148 and CD68) in BM-MSCs, representing potential novel markers for human BM-MSC purification. This study is the first systematic in vivo dissection of human BM-MSCs cell subtypes at the single-cell resolution, revealing an insight into the extent of their cellular heterogeneity and roles in maintaining bone homeostasis.

Severe COVID-19 in Alzheimer's disease: APOE4's fault again?

Challenges have been recognized in healthcare of patients with Alzheimer's disease (AD) in the COVID-19 pandemic, given a high infection and mortality rate of COVID-19 in these patients. This situation urges the identification of underlying risks and preferably biomarkers for evidence-based, more effective healthcare. Towards this goal, current literature review and network analysis synthesize available information on the AD-related gene APOE into four lines of mechanistic evidence. At a cellular level, the risk isoform APOE4 confers high infectivity by the underlying coronavirus SARS-CoV-2; at a genetic level, APOE4 is associated with severe COVID-19; at a pathway level, networking connects APOE with COVID-19 risk factors such as ACE2, TMPRSS2, NRP1, and LZTFL1; at a behavioral level, APOE4-associated dementia may increase the exposure to coronavirus infection which causes COVID-19. Thus, APOE4 could exert multiple actions for high infection and mortality rates of the patients, or generally, with COVID-19.

Identification and Functional Characterization of Metabolites for Bone Mass in Peri- and Postmenopausal Chinese Women.

CONTEXT: Although metabolic profiles appear to play an important role in menopausal bone loss, the functional mechanisms by which metabolites influence bone mineral density (BMD) during menopause are largely unknown. OBJECTIVE: We aimed to systematically identify metabolites associated with BMD variation and their potential functional mechanisms in peri- and postmenopausal women. DESIGN AND METHODS: We performed serum metabolomic profiling and whole-genome sequencing for 517 perimenopausal (16%) and early postmenopausal (84%) women aged 41 to 64 years in this cross-sectional study. Partial least squares regression and general linear regression analysis were applied to identify BMD-associated metabolites, and weighted gene co-expression network analysis was performed to construct co-functional metabolite modules. Furthermore, we performed Mendelian randomization analysis to identify causal relationships between BMD-associated metabolites and BMD variation. Finally, we explored the effects of a novel prominent BMD-associated metabolite on bone metabolism through both in vivo/in vitro experiments. RESULTS: Twenty metabolites and a co-functional metabolite module (consisting of fatty acids) were significantly associated with BMD variation. We found dodecanoic acid (DA), within the identified module causally decreased total hip BMD. Subsequently, the in vivo experiments might support that dietary supplementation with DA could promote bone loss, as well as increase the osteoblast and osteoclast numbers in normal/ovariectomized mice. Dodecanoic acid treatment differentially promoted osteoblast and osteoclast differentiation, especially for osteoclast differentiation at higher concentrations in vitro (eg,10, 100 µM). CONCLUSIONS: This study sheds light on metabolomic profiles associated with postmenopausal osteoporosis risk, highlighting the potential importance of fatty acids, as exemplified by DA, in regulating BMD.

Identifying Pleiotropic SNPs Associated With Femoral Neck and Heel Bone Mineral Density.

BACKGROUND: Genome-wide association studies (GWASs) routinely identify loci associated with risk factors for osteoporosis. However, GWASs with relatively small sample sizes still lack sufficient power to ascertain the majority of genetic variants with small to modest effect size, which may together truly influence the phenotype. The loci identified only account for a small percentage of the heritability of osteoporosis. This study aims to identify novel genetic loci associated with DXA-derived femoral neck (FNK) bone mineral density (BMD) and quantitative ultrasound of the heel calcaneus estimated BMD (eBMD), and to detect shared/causal variants for the two traits, to assess whether the SNPs or putative causal SNPs associated with eBMD were also associated with FNK-BMD. METHODS: Novel loci associated with eBMD and FNK-BMD were identified by the genetic pleiotropic conditional false discovery rate (cFDR) method. Shared putative causal variants between eBMD and FNK-BMD and putative causal SNPs for each trait were identified by the colocalization method. Mendelian randomization analysis addresses the causal relationship between eBMD/FNK-BMD and fracture. RESULTS: We identified 9,500 (cFDR < 9.8E-6), 137 (cFDR < 8.9E-4) and 124 SNPs associated with eBMD, FNK-BMD, and both eBMD and FNK-BMD, respectively, with 37 genomic regions where there was a SNP that influences both eBMD and FNK-BMD. Most genomic regions only contained putative causal SNPs associated with eBMD and 3 regions contained two distinct putative causal SNPs influenced both traits, respectively. We demonstrated a causal effect of FNK-BMD/eBMD on fracture. CONCLUSION: Most of SNPs or putative causal SNPs associated with FNK-BMD were also associated with eBMD. However, most of SNPs or putative causal SNPs associated with eBMD were not associated with FNK-BMD. The novel variants we identified may help to account for the additional proportion of variance of each trait and advance our understanding of the genetic mechanisms underlying osteoporotic fracture.

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.

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.

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.

Minimotifs dysfunction is pervasive in neurodegenerative disorders.

Minimotifs are modular contiguous peptide sequences in proteins that are important for posttranslational modifications, binding to other molecules, and trafficking to specific subcellular compartments. Some molecular functions of proteins in cellular pathways can be predicted from minimotif consensus sequences identified through experimentation. While a role for minimotifs in regulating signal transduction and gene regulation during disease pathogenesis (such as infectious diseases and cancer) is established, the therapeutic use of minimotif mimetic drugs is limited. In this review, we discuss a general theme identifying a pervasive role of minimotifs in the pathomechanism of neurodegenerative diseases. Beyond their longstanding history in the genetics of familial neurodegeneration, minimotifs are also major players in neurotoxic protein aggregation, aberrant protein trafficking, and epigenetic regulation. Generalizing the importance of minimotifs in neurodegenerative diseases offers a new perspective for the future study of neurodegenerative mechanisms and the investigation of new therapeutics.

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.

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 Functional Human C-Terminome.

All translated proteins end with a carboxylic acid commonly called the C-terminus. Many short functional sequences (minimotifs) are located on or immediately proximal to the C-terminus. However, information about the function of protein C-termini has not been consolidated into a single source. Here, we built a new "C-terminome" database and web system focused on human proteins. Approximately 3,600 C-termini in the human proteome have a minimotif with an established molecular function. To help evaluate the function of the remaining C-termini in the human proteome, we inferred minimotifs identified by experimentation in rodent cells, predicted minimotifs based upon consensus sequence matches, and predicted novel highly repetitive sequences in C-termini. Predictions can be ranked by enrichment scores or Gene Evolutionary Rate Profiling (GERP) scores, a measurement of evolutionary constraint. By searching for new anchored sequences on the last 10 amino acids of proteins in the human proteome with lengths between 3-10 residues and up to 5 degenerate positions in the consensus sequences, we have identified new consensus sequences that predict instances in the majority of human genes. All of this information is consolidated into a database that can be accessed through a C-terminome web system with search and browse functions for minimotifs and human proteins. A known consensus sequence-based predicted function is assigned to nearly half the proteins in the human proteome. Weblink: http://cterminome.bio-toolkit.com.