A natural substitution of a conserved amino acid in eIF4E confers resistance against multiple potyviruses.

Eukaryotic translation initiation factor 4E (eIF4E), which plays a pivotal role in initiating translation in eukaryotic organisms, is often hijacked by the viral genome-linked protein to facilitate the infection of potyviruses. In this study, we found that the naturally occurring amino acid substitution D71G in eIF4E is widely present in potyvirus-resistant watermelon accessions and disrupts the interaction between watermelon eIF4E and viral genome-linked protein of papaya ringspot virus-watermelon strain, zucchini yellow mosaic virus or watermelon mosaic virus. Multiple sequence alignment and protein modelling showed that the amino acid residue D71 located in the cap-binding pocket of eIF4E is strictly conserved in many plant species. The mutation D71G in watermelon eIF4E conferred resistance against papaya ringspot virus-watermelon strain and zucchini yellow mosaic virus, and the equivalent mutation D55G in tobacco eIF4E conferred resistance to potato virus Y. Therefore, our finding provides a potential precise target for breeding plants resistant to multiple potyviruses.

A Novel Vpu Adaptive Mutation of HIV-1 Degrades Tetherin in Northern Pig-Tailed Macaques (Macaca leonina) Mainly via the Ubiquitin-Proteasome Pathway and Increases Viral Release.

Tetherin prevents viral cross-species transmission by inhibiting the release of multiple enveloped viruses from infected cells. With the evolution of simian immunodeficiency virus of chimpanzees (SIVcpz), a pandemic human immunodeficiency virus type 1 (HIV-1) precursor, its Vpu protein can antagonize human tetherin (hTetherin). Macaca leonina (northern pig-tailed macaque [NPM]) is susceptible to HIV-1, but host-specific restriction factors limit virus replication in vivo. In this study, we isolated the virus from NPMs infected with strain stHIV-1sv (with a macaque-adapted HIV-1 env gene from simian-human immunodeficiency virus SHIV-KB9, a vif gene replaced by SIVmac239, and other genes originating from HIV-1NL4.3) and found that a single acidic amino acid substitution (G53D) in Vpu could increase its ability to degrade the tetherin of macaques (mTetherin) mainly through the proteasome pathway, resulting in an enhanced release and resistance to interferon inhibition of the mutant stHIV-1sv strain, with no influence on the other functions of Vpu. IMPORTANCE HIV-1 has obvious host specificity, which has greatly hindered the construction of animal models and severely restricted the development of HIV-1 vaccines and drugs. To overcome this barrier, we attempted to isolate the virus from NPMs infected with stHIV-1sv, search for a strain with an adaptive mutation in NPMs, and develop a more appropriate nonhuman primate model of HIV-1. This is the first report identifying HIV-1 adaptations in NPMs. It suggests that while tetherin may limit HIV-1 cross-species transmission, the Vpu protein in HIV-1 can overcome this species barrier through adaptive mutation, increasing viral replication in the new host. This finding will be beneficial to building an appropriate animal model for HIV-1 infection and promoting the development of HIV-1 vaccines and drugs.

Constraints on the evolution of toxin-resistant Na,K-ATPases have limited dependence on sequence divergence.

A growing body of theoretical and experimental evidence suggests that intramolecular epistasis is a major determinant of rates and patterns of protein evolution and imposes a substantial constraint on the evolution of novel protein functions. Here, we examine the role of intramolecular epistasis in the recurrent evolution of resistance to cardiotonic steroids (CTS) across tetrapods, which occurs via specific amino acid substitutions to the α-subunit family of Na,K-ATPases (ATP1A). After identifying a series of recurrent substitutions at two key sites of ATP1A that are predicted to confer CTS resistance in diverse tetrapods, we then performed protein engineering experiments to test the functional consequences of introducing these substitutions onto divergent species backgrounds. In line with previous results, we find that substitutions at these sites can have substantial background-dependent effects on CTS resistance. Globally, however, these substitutions also have pleiotropic effects that are consistent with additive rather than background-dependent effects. Moreover, the magnitude of a substitution's effect on activity does not depend on the overall extent of ATP1A sequence divergence between species. Our results suggest that epistatic constraints on the evolution of CTS-resistant forms of Na,K-ATPase likely depend on a small number of sites, with little dependence on overall levels of protein divergence. We propose that dependence on a limited number sites may account for the observation of convergent CTS resistance substitutions observed among taxa with highly divergent Na,K-ATPases (See S1 Text for Spanish translation).

Mutation in the CX3C Motif of G Protein Disrupts Its Interaction with Heparan Sulfate: A Calorimetric, Spectroscopic, and Molecular Docking Study.

Respiratory syncytial virus (RSV) is the leading cause of lower respiratory tract infection in children and infants. To date, there is no effective vaccine available against RSV. Heparan sulfate is a type of glycosaminoglycan that aids in the attachment of the RSV to the host cell membrane via the G protein. In the present study, the effect of amino acid substitution on the structure and stability of the ectodomain G protein was studied. Further, it was investigated whether mutation (K117A) in the CX3C motif of G protein alters the binding with heparan sulfate. The point mutation significantly affects the conformational stability of the G protein. The mutant protein showed a low binding affinity with heparan sulfate as compared to the wild-type G protein, as determined by fluorescence quenching, isothermal titration calorimetry (ITC), and molecular docking studies. The low binding affinity and decreased stability suggested that this mutation may play an important role in prevention of attachment of virion to the host cell receptors. Collectively, this investigation suggests that mutation in the CX3C motif of G protein may likely improve the efficacy and safety of the RSV vaccine.

Amino acid substitutions at the HIV-1 transframe region significantly impair virus infectivity.

A transframe region within HIV-1 Gag-Pol (referred to as p6* or p6pol), directly linked to the protease (PR) N-terminus, plays a pivotal role in modulating PR activation. To identify specific p6* residues involved in PR activation, we created a series of p6* mutants by making substitutions for conserved p6* residues. Our results indicate that some p6* mutants were defective in terms of virus infectivity, despite displaying a wild-type virus particle processing pattern. Mutations at p6* F8 reduced virus infectivity associated with insufficient virus processing, due in part to impaired PR maturation and RT packaging. Our data strongly suggest that conserved Phe (F) residues at position 8 of p6* are involved in the PR maturation process.

A Two-Level, Intramutant Spectrum Haplotype Profile of Hepatitis C Virus Revealed by Self-Organized Maps.

RNA viruses replicate as complex mutant spectra termed viral quasispecies. The frequency of each individual genome in a mutant spectrum depends on its rate of generation and its relative fitness in the replicating population ensemble. The advent of deep sequencing methodologies allows for the first-time quantification of haplotype abundances within mutant spectra. There is no information on the haplotype profile of the resident genomes and how the landscape evolves when a virus replicates in a controlled cell culture environment. Here, we report the construction of intramutant spectrum haplotype landscapes of three amplicons of the NS5A-NS5B coding region of hepatitis C virus (HCV). Two-dimensional (2D) neural networks were constructed for 44 related HCV populations derived from a common clonal ancestor that was passaged up to 210 times in human hepatoma Huh-7.5 cells in the absence of external selective pressures. The haplotype profiles consisted of an extended dense basal platform, from which a lower number of protruding higher peaks emerged. As HCV increased its adaptation to the cells, the number of haplotype peaks within each mutant spectrum expanded, and their distribution shifted in the 2D network. The results show that extensive HCV replication in a monotonous cell culture environment does not limit HCV exploration of sequence space through haplotype peak movements. The landscapes reflect dynamic variation in the intramutant spectrum haplotype profile and may serve as a reference to interpret the modifications produced by external selective pressures or to compare with the landscapes of mutant spectra in complex in vivo environments. IMPORTANCE The study provides for the first time the haplotype profile and its variation in the course of virus adaptation to a cell culture environment in the absence of external selective constraints. The deep sequencing-based self-organized maps document a two-layer haplotype distribution with an ample basal platform and a lower number of protruding peaks. The results suggest an inferred intramutant spectrum fitness landscape structure that offers potential benefits for virus resilience to mutational inputs.

Intermittent macrothrombocytopenia in a novel patient with Takenouchi-Kosaki syndrome and review of literature.

Takenouchi-Kosaki syndrome (TKS) is a recently delineated syndromic form of thrombocytopenia strictly related to an hot-spot missense variant, p.Tyr64Cys, in CDC42 (Cell Division Control protein 42). Herein we report an additional patient with the p.Tyr64Cys aminoacidic substitution who showed the well-defined phenotypical TKS features and an intermittent, very mild, macrothrombocytopenia at 10.7 years of age (93,000/mL), that was only retrospectively valorized. Outside of this value the PLT count had always been higher than 100,000/mL. We also review literature data from patients carrying this recurrent variant. Our female patient presented with prenatal onset of short stature and microcephaly, camptodactyly, heart defects, typical facial gestalt, developmental delay, and not specific brain abnormalities. After several genetic investigations (karyotype, CGH-Array, targeted NGS analysis for short stature genes), by whole exome sequencing we identified the p.Tyr64Cys in CDC42, occurring de novo. The case presented here provides further evidence that macrothrombocytopenia can be intermittent and thus it might escape attention of clinicians. Without this key feature, TKS clinical presentation can overlap other syndromic forms of short stature. Immunodeficiency, autoimmunity, and malignancies were recently reported in patients with the p.Tyr64Cys substitution, making imperative an early diagnosis of Takenouchi-Kosaki syndrome to organize the most proper follow-up of these pediatric patients. The whole exome sequencing can be a solving tool in the challenge to the rare diseases.

Compartmentalization of Resistance-Associated Substitutions in HIV/HCV-Infected Patients: Possible Correlation with Infecting HCV Genotype.

Resistance-associated substitutions (RASs) may exist prior to treatment and contribute to the failure of treatment with direct-acting antivirals (DAAs). As the major site of HCV replication, naturally occurring variants with RASs may segregate into the liver. In the present study, we performed viral population sequencing to retrospectively investigate the NS3 and NS5A RAS profiles in 34 HIV/HCV coinfected patients naïve to anti-HCV treatment who underwent diagnostic liver biopsy between 2000 and 2006 and had liver and plasma samples available. Sixteen were infected by HCV genotype (GT) 1a, 11 by GT3a, and 7 by GT4d. The analysis of the NS3 domain in GT1a showed a difference in strain between the liver and plasma in three cases, with a preponderance of specific RASs in the liver compartment. In GT4d samples, 6/7 coupled liver and plasma samples were concordant with no RASs. Sequence analysis of the NS5A domain showed the presence of RASs in the livers of 2/16 patients harboring GT1a but not in the corresponding plasma. In GT4d, NS5A RASs were detected in 7/7 liver tissues and 5/7 plasma samples. NS3 domain and NS5A domain were found to be conserved in plasma and livers of patients infected with GT3a. Thus, RASs within GT1a and GT4d more likely segregate into the liver and may explain the emergence of resistant strains during DAA treatment.

Overexpression of ERAP2N in Human Trophoblast Cells Promotes Cell Death.

The genes involved in implantation and placentation are tightly regulated to ensure a healthy pregnancy. The endoplasmic reticulum aminopeptidase 2 (ERAP2) gene is associated with preeclampsia (PE). Our studies have determined that an isoform of ERAP2-arginine (N), expressed in trophoblast cells (TC), significantly activates immune cells, and ERAP2N-expressing TCs are preferentially killed by both cytotoxic T lymphocytes (CTLs) and Natural Killer cells (NKCs). To understand the cause of this phenomenon, we surveyed differentially expressed genes (DEGs) between ERAP2N expressing and non-expressing TCs. Our RNAseq data revealed 581 total DEGs between the two groups. 289 genes were up-regulated, and 292 genes were down-regulated. Interestingly, most of the down-regulated genes of significance were pro-survival genes that play a crucial role in cell survival (LDHA, EGLN1, HLA-C, ITGB5, WNT7A, FN1). However, the down-regulation of these genes in ERAP2N-expressing TCs translates into a propensity for cell death. The Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that 64 DEGs were significantly enriched in nine pathways, including "Protein processing in endoplasmic reticulum" and "Antigen processing and presentation", suggesting that the genes may be associated with peptide processes involved in immune recognition during the reproductive cycle.

Neuroprotective effects of bone marrow Sca-1+ cells against age-related retinal degeneration in OPTN E50K mice.

Glaucoma is characterized by retinal ganglion cell (RGC) death, the underlying mechanisms of which are still largely unknown. An E50K mutation in the Optineurin (OPTN) gene is a leading cause of normal-tension glaucoma (NTG), which directly affects RGCs in the absence of high intraocular pressure and causes severe glaucomatous symptoms in patients. Bone marrow (BM) stem cells have been demonstrated to play a key role in regenerating damaged tissue during ageing and disease through their trophic effects and homing capability. Here, we separated BM stem cells into Sca-1+ and Sca-1- cells and transplanted them into lethally irradiated aged OPTN E50K mice to generate Sca-1+ and Sca-1- chimaeras, respectively. After 3 months of BM repopulation, we investigated whether Sca-1+ cells maximized the regenerative effects in the retinas of NTG model mice with the OPTN E50K mutation. We found that the OPTN E50K mutation aggravated age-related deficiency of neurotrophic factors in both retinas and BM during NTG development, leading to retinal degeneration and BM dysfunction. Sca-1+ cells from young healthy mice had greater paracrine trophic effects than Sca-1- cells and Sca-1+ cells from young OPTN E50K mice. In addition, Sca-1+ chimaeras demonstrated better visual functions than Sca-1- chimaeras and untreated OPTN E50K mice. More Sca-1+ cells than Sca-1- cells were recruited to repair damaged retinas and reverse visual impairment in NTG resulting from high expression levels of neurotrophic factors. These findings indicated that the Sca-1+ cells from young, healthy mice may have exhibited an enhanced ability to repair retinal degeneration in NTG because of their excellent neurotrophic capability.

Fitness selection of hyperfusogenic measles virus F proteins associated with neuropathogenic phenotypes.

Measles virus (MeV) is resurgent and caused >200,000 deaths in 2019. MeV infection can establish a chronic latent infection of the brain that can recrudesce months to years after recovery from the primary infection. Recrudescent MeV leads to fatal subacute sclerosing panencephalitis (SSPE) or measles inclusion body encephalitis (MIBE) as the virus spreads across multiple brain regions. Most clinical isolates of SSPE/MIBE strains show mutations in the fusion (F) gene that result in a hyperfusogenic phenotype in vitro and allow for efficient spread in primary human neurons. Wild-type MeV receptor-binding protein is indispensable for manifesting these mutant F phenotypes, even though neurons lack canonical MeV receptors (CD150/SLAMF1 or nectin-4). How such hyperfusogenic F mutants are selected and whether they confer a fitness advantage for efficient neuronal spread is unresolved. To better understand the fitness landscape that allows for the selection of such hyperfusogenic F mutants, we conducted a screen of ≥3.1 × 105 MeV-F point mutants in their genomic context. We rescued and amplified our genomic MeV-F mutant libraries in BSR-T7 cells under conditions in which MeV-F-T461I (a known SSPE mutant), but not wild-type MeV, can spread. We recovered known SSPE mutants but also characterized at least 15 hyperfusogenic F mutations with an SSPE phenotype. Structural mapping of these mutants onto the prefusion MeV-F trimer confirm and extend our understanding of the F regulatory domains in MeV-F. Our list of hyperfusogenic F mutants is a valuable resource for future studies into MeV neuropathogenesis and the regulation of paramyxovirus F.

Protein kinase A controls the hexosamine pathway by tuning the feedback inhibition of GFAT-1.

The hexosamine pathway (HP) is a key anabolic pathway whose product uridine 5'-diphospho-N-acetyl-D-glucosamine (UDP-GlcNAc) is an essential precursor for glycosylation processes in mammals. It modulates the ER stress response and HP activation extends lifespan in Caenorhabditis elegans. The highly conserved glutamine fructose-6-phosphate amidotransferase 1 (GFAT-1) is the rate-limiting HP enzyme. GFAT-1 activity is modulated by UDP-GlcNAc feedback inhibition and via phosphorylation by protein kinase A (PKA). Molecular consequences of GFAT-1 phosphorylation, however, remain poorly understood. Here, we identify the GFAT-1 R203H substitution that elevates UDP-GlcNAc levels in C. elegans. In human GFAT-1, the R203H substitution interferes with UDP-GlcNAc inhibition and with PKA-mediated Ser205 phosphorylation. Our data indicate that phosphorylation affects the interactions of the two GFAT-1 domains to control catalytic activity. Notably, Ser205 phosphorylation has two discernible effects: it lowers baseline GFAT-1 activity and abolishes UDP-GlcNAc feedback inhibition. PKA controls the HP by uncoupling the metabolic feedback loop of GFAT-1.

Charged Residues in the Membrane Anchor of the Pestiviral Erns Protein Are Important for Processing and Secretion of Erns and Recovery of Infectious Viruses.

The pestivirus envelope protein Erns is anchored in membranes via a long amphipathic helix. Despite the unusual membrane topology of the Erns membrane anchor, it is cleaved from the following glycoprotein E1 by cellular signal peptidase. This was proposed to be enabled by a salt bridge-stabilized hairpin structure (so-called charge zipper) formed by conserved charged residues in the membrane anchor. We show here that the exchange of one or several of these charged residues reduces processing at the Erns carboxy-terminus to a variable extend, but reciprocal mutations restoring the possibility to form salt bridges did not necessarily restore processing efficiency. When introduced into an Erns-only expression construct, these mutations enhanced the naturally occurring Erns secretion significantly, but again to varying extents that did not correlate with the number of possible salt bridges. Equivalent effects on both processing and secretion were also observed when the proteins were expressed in avian cells, which points at phylogenetic conservation of the underlying principles. In the viral genome, some of the mutations prevented recovery of infectious viruses or immediately (pseudo)reverted, while others were stable and neutral with regard to virus growth.

Systematic investigation of PRMT6 substrate recognition reveals broad specificity with a preference for an RG motif or basic and bulky residues.

Protein arginine methyltransferase 6 (PRMT6) catalyses the asymmetric dimethylation of arginines on numerous substrate proteins within the human cell. In particular, PRMT6 methylates histone H3 arginine 2 (H3R2) which affects both gene repression and activation. However, the substrate specificity of PRMT6 has not been comprehensively analysed. Here, we systematically characterise the substrate recognition motif of PRMT6, finding that it has broad specificity and recognises the RG motif. Working with a H3 tail peptide as a template, on which we made 204 amino acid substitutions, we use targeted mass spectrometry to measure their effect on PRMT6 in vitro activity. We first show that PRMT6 methylates R2 and R8 in the H3 peptide, although H3R8 is methylated with lower efficiency and is not an in vivo PRMT6 substrate. We then quantify the effect of 194 of these amino acid substitutions on methylation at both H3R2 and H3R8. In both cases, we find that PRMT6 tolerates essentially any amino acid substitution in the H3 peptide, but that positively charged and bulky residues are preferred near the target arginine. We show that PRMT6 also has preference for glycine, but only in the position immediately following the target arginine. This indicates that PRMT6 recognises the RG motif rather than the RGG motif. We further confirm this preference for the RG motif on another PRMT6 substrate, histone H4R3. This broad specificity and recognition of RG rather than RGG are distinctive among the PRMT family and has implications for the development of drugs to selectively target PRMT6. DATABASES: Panorama Public (https://panoramaweb.org/PRMT6motif.url); ProteomeXchange (PXD016711).

A Glycoprotein Mutation That Emerged during the 2013-2016 Ebola Virus Epidemic Alters Proteolysis and Accelerates Membrane Fusion.

Genomic surveillance of viral isolates during the 2013-2016 Ebola virus epidemic in Western Africa, the largest and most devastating filovirus outbreak on record, revealed several novel mutations. The responsible strain, named Makona, carries an A-to-V substitution at position 82 (A82V) in the glycoprotein (GP), which is associated with enhanced infectivity in vitro Here, we investigated the mechanistic basis for this enhancement as well as the interplay between A82V and a T-to-I substitution at residue 544 of GP, which also modulates infectivity in cell culture. We found that both 82V and 544I destabilize GP, with the residue at position 544 impacting overall stability, while 82V specifically destabilizes proteolytically cleaved GP. Both residues also promote faster kinetics of lipid mixing of the viral and host membranes in live cells, individually and in tandem, which correlates with faster times to fusion following colocalization with the viral receptor Niemann-Pick C1 (NPC1). Furthermore, GPs bearing 82V are more sensitive to proteolysis by cathepsin L (CatL), a key host factor for viral entry. Intriguingly, CatL processed 82V variant GPs to a novel product with a molecular weight of approximately 12,000 (12K), which we hypothesize corresponds to a form of GP that is pre-triggered for fusion. We thus propose a model in which 82V promotes more efficient GP processing by CatL, leading to faster viral fusion kinetics and higher levels of infectivity.IMPORTANCE The 2013-2016 outbreak of Ebola virus disease in West Africa demonstrated the potential for previously localized outbreaks to turn into regional, or even global, health emergencies. With over 28,000 cases and 11,000 confirmed deaths, this outbreak was over 50 times as large as any previously recorded. This outbreak also afforded the largest-ever collection of Ebola virus genomic sequence data, allowing new insights into viral transmission and evolution. Viral mutants arising during the outbreak have attracted attention for their potentially altered patterns of infectivity in cell culture, with potential, if unclear, implications for increased viral spread and/or virulence. Here, we report the properties of one such mutation in the viral glycoprotein, A82V, and its interplay with a previously described polymorphism at position 544. We show that mutations at both residues promote infection and fusion activation in cells but that A82V additionally leads to increased infectivity under cathepsin-limited conditions and the generation of a novel glycoprotein cleavage product.

Molecular basis of accessible plasma membrane cholesterol recognition by the GRAM domain of GRAMD1b.

Cholesterol is essential for cell physiology. Transport of the "accessible" pool of cholesterol from the plasma membrane (PM) to the endoplasmic reticulum (ER) by ER-localized GRAMD1 proteins (GRAMD1a/1b/1c) contributes to cholesterol homeostasis. However, how cells detect accessible cholesterol within the PM remains unclear. We show that the GRAM domain of GRAMD1b, a coincidence detector for anionic lipids, including phosphatidylserine (PS), and cholesterol, possesses distinct but synergistic sites for sensing accessible cholesterol and anionic lipids. We find that a mutation within the GRAM domain of GRAMD1b that is associated with intellectual disability in humans specifically impairs cholesterol sensing. In addition, we identified another point mutation within this domain that enhances cholesterol sensitivity without altering its PS sensitivity. Cell-free reconstitution and cell-based assays revealed that the ability of the GRAM domain to sense accessible cholesterol regulates membrane tethering and determines the rate of cholesterol transport by GRAMD1b. Thus, cells detect the codistribution of accessible cholesterol and anionic lipids in the PM and fine-tune the non-vesicular transport of PM cholesterol to the ER via GRAMD1s.

Temperature-dependent autoactivation associated with clinical variability of PDGFRB Asn666 substitutions.

Ocular pterygium-digital keloid dysplasia (OPDKD) presents in childhood with ingrowth of vascularized connective tissue on the cornea leading to severely reduced vision. Later the patients develop keloids on digits but are otherwise healthy. The overgrowth in OPDKD affects body parts that typically have lower temperature than 37°C. We present evidence that OPDKD is associated with a temperature sensitive, activating substitution, p.(Asn666Tyr), in PDGFRB. Phosphorylation levels of PDGFRB and downstream targets were higher in OPDKD fibroblasts at 37°C but were further greatly increased at the average corneal temperature of 32°C. This suggests that the substitution cause significant constitutive autoactivation mainly at lower temperature. In contrast, a different substitution in the same codon, p.(Asn666Ser), is associated with Penttinen type of premature aging syndrome. This devastating condition is characterized by widespread tissue degeneration, including pronounced chronic ulcers and osteolytic resorption in distal limbs. In Penttinen syndrome fibroblasts, equal and high levels of phosphorylated PDGFRB was present at both 32°C and 37°C. This indicates that this substitution causes severe constitutive autoactivation of PDGFRB regardless of temperature. In line with this, most downstream targets were not affected by lower temperature. However, STAT1, important for tissue wasting, did show further increased phosphorylation at 32°C. Temperature-dependent autoactivation offers an explanation to the strikingly different clinical outcomes of substitutions in the Asn666 codon of PDGFRB.

Alignment-free method for functional annotation of amino acid substitutions: Application on epigenetic factors involved in hematologic malignancies.

For the last couple of decades, there has been a significant growth in sequencing data, leading to an extraordinary increase in the number of gene variants. This places a challenge on the bioinformatics research community to develop and improve computational tools for functional annotation of new variants. Genes coding for epigenetic regulators have important roles in cancer pathogenesis and mutations in these genes show great potential as clinical biomarkers, especially in hematologic malignancies. Therefore, we developed a model that specifically focuses on these genes, with an assumption that it would outperform general models in predicting the functional effects of amino acid substitutions. EpiMut is a standalone software that implements a sequence based alignment-free method. We applied a two-step approach for generating sequence based features, relying on the biophysical and biochemical indices of amino acids and the Fourier Transform as a sequence transformation method. For each gene in the dataset, the machine learning algorithm-Naïve Bayes was used for building a model for prediction of the neutral or disease-related status of variants. EpiMut outperformed state-of-the-art tools used for comparison, PolyPhen-2, SIFT and SNAP2. Additionally, EpiMut showed the highest performance on the subset of variants positioned outside conserved functional domains of analysed proteins, which represents an important group of cancer-related variants. These results imply that EpiMut can be applied as a first choice tool in research of the impact of gene variants in epigenetic regulators, especially in the light of the biomarker role in hematologic malignancies. EpiMut is freely available at https://www.vin.bg.ac.rs/180/tools/epimut.php.

Missense3D-DB web catalogue: an atom-based analysis and repository of 4M human protein-coding genetic variants.

The interpretation of human genetic variation is one of the greatest challenges of modern genetics. New approaches are urgently needed to prioritize variants, especially those that are rare or lack a definitive clinical interpretation. We examined 10,136,597 human missense genetic variants from GnomAD, ClinVar and UniProt. We were able to perform large-scale atom-based mapping and phenotype interpretation of 3,960,015 of these variants onto 18,874 experimental and 84,818 in house predicted three-dimensional coordinates of the human proteome. We demonstrate that 14% of amino acid substitutions from the GnomAD database that could be structurally analysed are predicted to affect protein structure (n = 568,548, of which 566,439 rare or extremely rare) and may, therefore, have a yet unknown disease-causing effect. The same is true for 19.0% (n = 6266) of variants of unknown clinical significance or conflicting interpretation reported in the ClinVar database. The results of the structural analysis are available in the dedicated web catalogue Missense3D-DB ( http://missense3d.bc.ic.ac.uk/ ). For each of the 4 M variants, the results of the structural analysis are presented in a friendly concise format that can be included in clinical genetic reports. A detailed report of the structural analysis is also available for the non-experts in structural biology. Population frequency and predictions from SIFT and PolyPhen are included for a more comprehensive variant interpretation. This is the first large-scale atom-based structural interpretation of human genetic variation and offers geneticists and the biomedical community a new approach to genetic variant interpretation.

Role of core protein mutations in the development of occult HBV infection.

BACKGROUND & AIMS: Occult HBV infection (OBI) is associated with transfusion-transmitted HBV infection and hepatocellular carcinoma. Studies on OBI genesis have concentrated on mutations in the S region and the regulatory elements. Herein, we aimed to determine the role of mutations in the core region on OBIs. METHODS: An OBI strain (SZA) carrying 9 amino acid (aa) substitutions in the core protein/capsid (Cp) was selected by sequence alignment and Western blot analysis from 26 genotype B OBI samples to extensively explore the impact of Cp mutations on viral antigen production in vitro and in vivo. RESULTS: A large panel of 30 Cp replicons were generated by a replication-competent pHBV1.3 carrying SZA or wild-type (WT) Cp in a 1.3-fold over-length of HBV genome, in which the various Cp mutants were individually introduced by repairing site mutations of SZA-Cp or creating site mutations of WT-Cp by site-directed mutagenesis. The expression of HBcAg, HBeAg, and HBsAg and viral RNA was quantified from individual SZA and WT Cp mutant replicons in transfected Huh7 cells or infected mice, respectively. An analysis of the effect of Cp mutants on intracellular or extracellular viral protein production indicated that the W62R mutation in Cp had a critical impact on the reduction of HBcAg and HBeAg production during HBV replication, whereas P50H and/or S74G mutations played a limited role in influencing viral protein production invivo. CONCLUSIONS: W62R and its combination mutations in HBV Cp might massively affect HBcAg and HBeAg production during viral replication, which, in turn, might contribute to the occurrence of OBI. LAY SUMMARY: Occult hepatitis B virus infections (OBIs) have been found to be associated with amino acid mutations in the S region of the HBV, but the role of mutations in the core protein (Cp) remains unclear. In this study, an OBI strain (SZA) carrying 9 amino acid substitutions in Cp has been examined comprehensively in vitro and in vivo. The W62R mutation in Cp majorly reduces HBcAg and HBeAg production during HBV replication, potentially contributing to the occurrence of OBI.