Search

The spike (S) glycoprotein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the binding to the permissive cells. The receptor-binding domain (RBD) of SARS-CoV-2 S protein directly interacts with the human angiotensin-converting enzyme 2 (ACE2) on the host cell membrane. In this study, we used computational saturation mutagenesis approaches, including structure-based energy calculations and sequence-based pathogenicity predictions, to quantify the systemic effects of missense mutations on SARS-CoV-2 S protein structure and function. A total of 18 354 mutations in S protein were analyzed, and we discovered that most of these mutations could destabilize the entire S protein and its RBD. Specifically, residues G431 and S514 in SARS-CoV-2 RBD are important for S protein stability. We analyzed 384 experimentally verified S missense variations and revealed that the dominant pandemic form, D614G, can stabilize the entire S protein. Moreover, many mutations in N-linked glycosylation sites can increase the stability of the S protein. In addition, we investigated 3705 mutations in SARS-CoV-2 RBD and 11 324 mutations in human ACE2 and found that SARS-CoV-2 neighbor residues G496 and F497 and ACE2 residues D355 and Y41 are critical for the RBD–ACE2 interaction. The findings comprehensively provide potential target sites in the development of drugs and vaccines against COVID-19.

ACE2 plays a critical role in SARS-CoV-2 infection to cause COVID-19 and SARS-CoV-2 spike protein binds to ACE2 and probably functionally inhibits ACE2 to aggravate the underlying diseases of COVID-19. The important factors that affect the severity and fatality of COVID-19 include patients' underlying diseases and ages. Therefore, particular care to the patients with underlying diseases is needed during the treatment of COVID-19 patients.

Coronavirus possesses the largest RNA genome among all the RNA viruses. Its genome encodes about 29 proteins. Most of the viral proteins are non-structural proteins (NSP) except envelop (E), membrane (M), nucleocapsid (N) and Spike (S) proteins that constitute the viral nucleocapsid, envelop and surface. We have recently cloned all the 29 SARS-CoV-2 genes into vectors for their expressions in mammalian cells except NSP11 that has only 14 amino acids (aa). We are able to express all the 28 cloned SARS-CoV-2 genes in human cells to characterize their subcellular distributions. The proteins of SARS-CoV-2 are mostly cytoplasmic but some are both cytoplasmic and nuclear. Those punctate staining proteins were further investigated by immunofluorescent assay (IFA) using specific antibodies or by co-transfection with an organelle marker-expressing plasmid. As a result, we found that NSP15, ORF6, M and ORF7a are related to Golgi apparatus, and that ORF7b, ORF8 and ORF10 colocalize with endoplasmic reticulum (ER). Interestingly, ORF3a distributes in cell membrane, early endosome, endosome, late endosome and lysosome, which suggests that ORF3a might help the infected virus to usurp endosome and lysosome for viral use. Furthermore, we revealed that NSP13 colocalized with SC35, a protein standing for splicing compartments in the nucleus. Our studies for the first time visualized the subcellular locations of SARS-CoV-2 proteins and might provide novel insights into the viral proteins’ biological functions.

Although the genetic contribution is under debate, biological studies in multiple mouse models have suggested that the Disrupted-in-Schizophrenia-1 (DISC1) protein may contribute to susceptibility to psychiatric disorders. In the present study, we took the advantages of the Drosophila model to dissect the molecular pathways that can be affected by DISC1 in the context of pathology-related phenotypes. We found that three pathways that include the homologs of Drosophila Dys, Trio, and Shot were downregulated by introducing a C-terminal truncated mutant DISC1. Consistently, these three molecules were downregulated in the induced pluripotent stem cell-derived forebrain neurons from the subjects carrying a frameshift deletion in DISC1 C-terminus. Importantly, the three pathways were underscored in the pathophysiology of psychiatric disorders in bioinformatics analysis. Taken together, our findings are in line with the polygenic theory of psychiatric disorders.

BackgroundLymphotoxin-α (LTα), located in the Major Histocompatibility Complex (MHC) class III region on chromosome 6, encodes a cytotoxic protein that mediates a variety of antiviral responses among other biological functions. Furthermore, several genotypes at this gene have been implicated in the onset of a number of complex diseases, including myocardial infarction, autoimmunity, and various types of cancer. However, little is known about levels of nucleotide variation and linkage disequilibrium (LD) in and nearLTα, which could also influence phenotypic variance. To address this gap in knowledge, we examined sequence variation across ~ 10 kilobases (kbs), encompassingLTαand the upstream region, in 2039 individuals from the 1000 Genomes Project originating from 21 global populations.ResultsHere, we observed striking patterns of diversity, including an excess of intermediate-frequency alleles, the maintenance of multiple common haplotypes and a deep coalescence time for variation (dating > 1.0 million years ago), in global populations. While these results are generally consistent with a model of balancing selection, we also uncovered a signature of positive selection in the form of long-range LD on chromosomes with derived alleles primarily in Eurasian populations. To reconcile these findings, which appear to support different models of selection, we argue that selective sweeps (particularly, soft sweeps) of multiple derived alleles in and/or nearLTαoccurred in non-Africans after their ancestors left Africa. Furthermore, these targets of selection were predicted to alter transcription factor binding site affinity and protein stability, suggesting they play a role in gene function. Additionally, our data also showed that a subset of these functional adaptive variants are present in archaic hominin genomes.ConclusionsOverall, this study identified candidate functional alleles in a biologically-relevant genomic region, and offers new insights into the evolutionary origins of these loci in modern human populations.

Lymphotoxin-α (LTα), located in the Major Histocompatibility Complex (MHC) class III region on chromosome 6, encodes a cytotoxic protein that mediates a variety of antiviral responses among other biological functions. Furthermore, several genotypes at this gene have been implicated in the onset of a number of complex diseases, including myocardial infarction, autoimmunity, and various types of cancer. However, little is known about levels of nucleotide variation and linkage disequilibrium (LD) in and near LTα, which could also influence phenotypic variance. To address this gap in knowledge, we examined sequence variation across ~ 10 kilobases (kbs), encompassing LTα and the upstream region, in 2039 individuals from the 1000 Genomes Project originating from 21 global populations.Here, we observed striking patterns of diversity, including an excess of intermediate-frequency alleles, the maintenance of multiple common haplotypes and a deep coalescence time for variation (dating 1.0 million years ago), in global populations. While these results are generally consistent with a model of balancing selection, we also uncovered a signature of positive selection in the form of long-range LD on chromosomes with derived alleles primarily in Eurasian populations. To reconcile these findings, which appear to support different models of selection, we argue that selective sweeps (particularly, soft sweeps) of multiple derived alleles in and/or near LTα occurred in non-Africans after their ancestors left Africa. Furthermore, these targets of selection were predicted to alter transcription factor binding site affinity and protein stability, suggesting they play a role in gene function. Additionally, our data also showed that a subset of these functional adaptive variants are present in archaic hominin genomes.Overall, this study identified candidate functional alleles in a biologically-relevant genomic region, and offers new insights into the evolutionary origins of these loci in modern human populations.

Hereditary spastic paraplegias (HSPs) are a genetically heterogeneous collection of neurodegenerative disorders categorized by progressive lower-limb spasticity and frailty. The complex HSP forms are characterized by various neurological features including progressive spastic weakness, urinary sphincter dysfunction, extra pyramidal signs and intellectual disability (ID). The kinesin superfamily proteins (KIFs) are microtubule-dependent molecular motors involved in intracellular transport. Kinesins directionally transport membrane vesicles, protein complexes, and mRNAs along neurites, thus playing important roles in neuronal development and function. Recent genetic studies have identified kinesin mutations in patients with HSPs. In this study, we used the computational approaches to investigate the 40 missense mutations associated with HSP and ID in KIF1A and KIF5A. We performed homology modeling to construct the structures of kinesin–microtubule binding domain and kinesin–tubulin complex. We applied structure-based energy calculation methods to determine the effects of missense mutations on protein stability and protein–protein interaction. The results revealed that the most of disease-causing mutations could change the folding free energy of kinesin motor domain and the binding free energy of kinesin–tubulin complex. We found that E253K associated with ID in KIF1A decrease the protein stability of kinesin motor domains. We showed that the HSP mutations located in kinesin–tubulin complex interface, such as K253N and R280C in KIF5A, can destabilize the kinesin–tubulin complex. The computational analysis provides useful information for understanding the roles of kinesin mutations in the development of ID and HSPs.

The ability to predict if a given mutation is disease-causing or not has enormous potential to impact human health. Typically, these predictions are made by assessing the effects of mutation on macromolecular stability and amino acid conservation. Here we report a novel feature: the electrostatic component of the force acting between a kinesin motor domain and tubulin. We demonstrate that changes in the electrostatic component of the binding force are able to discriminate between disease-causing and non-disease-causing mutations found in human kinesin motor domains using the receiver operating characteristic (ROC). Because diseases may originate from multiple effects not related to kinesin-microtubule binding, the prediction rate of 0.843 area under the ROC plot due to the change in magnitude of the electrostatic force alone is remarkable. These results reflect the dependence of kinesin’s function on motility along the microtubule, which suggests a precise balance of microtubule binding forces is required.

Schizophrenia (SCZ), bipolar disorder (BD) and recurrent major depressive disorder (rMDD) are common psychiatric illnesses. All have been associated with lower cognitive ability, and show evidence of genetic overlap and substantial evidence of pleiotropy with cognitive function and neuroticism. Disrupted in schizophrenia 1 (DISC1) protein directly interacts with a large set of proteins (DISC1 Interactome) that are involved in brain development and signaling. Modulation of DISC1 expression alters the expression of a circumscribed set of genes (DISC1 Regulome) that are also implicated in brain biology and disorder. Here we report targeted sequencing of 59 DISC1 Interactome genes and 154 Regulome genes in 654 psychiatric patients and 889 cognitively-phenotyped control subjects, on whom we previously reported evidence for trait association from complete sequencing of the DISC1 locus. Burden analyses of rare and singleton variants predicted to be damaging were performed for psychiatric disorders, cognitive variables and personality traits. The DISC1 Interactome and Regulome showed differential association across the phenotypes tested. After family-wise error correction across all traits (FWERacross), an increased burden of singleton disruptive variants in the Regulome was associated with SCZ (FWERacross P=0.0339). The burden of singleton disruptive variants in the DISC1 Interactome was associated with low cognitive ability at age 11 (FWERacross P=0.0043). These results identify altered regulation of schizophrenia candidate genes by DISC1 and its core Interactome as an alternate pathway for schizophrenia risk, consistent with the emerging effects of rare copy number variants associated with intellectual disability.

Schizophrenia (SCZ), bipolar disorder (BD) and recurrent major depressive disorder (rMDD) are common psychiatric illnesses. All have been associated with lower cognitive ability, and show evidence of genetic overlap and substantial evidence of pleiotropy with cognitive function and neuroticism. Disrupted in schizophrenia 1 (DISC1) protein directly interacts with a large set of proteins (DISC1 Interactome) that are involved in brain development and signaling. Modulation of DISC1 expression alters the expression of a circumscribed set of genes (DISC1 Regulome) that are also implicated in brain biology and disorder. Here, we report targeted sequencing of 59 DISC1 Interactome genes and 154 Regulome genes in 654 psychiatric patients and 889 cognitively-phenotyped control subjects, on whom we previously reported evidence for trait association from complete sequencing of the DISC1 locus. Burden analyses of rare and singleton variants predicted to be damaging were performed for psychiatric disorders, cognitive variables and personality traits. The DISC1 Interactome and Regulome showed differential association across the phenotypes tested. After family-wise error correction across all traits (FWERacross), an increased burden of singleton disruptive variants in the Regulome was associated with SCZ (FWERacrossP=0.0339). The burden of singleton disruptive variants in the DISC1 Interactome was associated with low cognitive ability at age 11 (FWERacrossP=0.0043). These results suggest that variants in the DISC1 Interactome effect the risk of psychiatric illness through altered expression of schizophrenia-associated genes. The biological impact of rare variants highlighted here merit further study.

Recent technological advances in next-generation sequencing (NGS) provide unprecedented power to sequence personal genomes, characterize genomic landscapes, and detect a large number of sequence variants. The discovery of disease-causing variants in patients’ genomes has dramatically changed our perspective on precision medicine. This chapter provides an overview of sequence variant detection and analysis in NGS study. We outline the general methods for identifying different types of sequence variants from NGS data. We summarize the common approaches for analyzing and visualizing casual variants associated with complex diseases on precision medicine informatics.

Schizophrenia (SCZ) is a serious psychiatric disorder that affects 1% of general population and places a heavy burden worldwide. The underlying genetic mechanism of SCZ remains unknown, but studies indicate that the disease is associated with a global gene expression disturbance across many genes. Next-generation sequencing, particularly of RNA sequencing (RNA-Seq), provides a powerful genome-scale technology to investigate the pathological processes of SCZ. RNA-Seq has been used to analyze the gene expressions and identify the novel splice isoforms and rare transcripts associated with SCZ. This paper provides an overview on the genetics of SCZ, the advantages of RNA-Seq for transcriptome analysis, the accomplishments of RNA-Seq in SCZ cohorts, and the applications of induced pluripotent stem cells and RNA-Seq in SCZ research.

Many studies have shown that missense mutations might play an important role in carcinogenesis. However, the extent to which cancer mutations might affect biomolecular interactions remains unclear. Here, we map glioblastoma missense mutations on the human protein interactome, model the structures of affected protein complexes and decipher the effect of mutations on protein-protein, protein-nucleic acid and protein-ion binding interfaces. Although some missense mutations over-stabilize protein complexes, we found that the overall effect of mutations is destabilizing, mostly affecting the electrostatic component of binding energy. We also showed that mutations on interfaces resulted in more drastic changes of amino acid physico-chemical properties than mutations occurring outside the interfaces. Analysis of glioblastoma mutations on interfaces allowed us to stratify cancer-related interactions, identify potential driver genes, and propose two dozen additional cancer biomarkers, including those specific to functions of the nervous system. Such an analysis also offered insight into the molecular mechanism of the phenotypic outcomes of mutations, including effects on complex stability, activity, binding and turnover rate. As a result of mutated protein and gene network analysis, we observed that interactions of proteins with mutations mapped on interfaces had higher bottleneck properties compared to interactions with mutations elsewhere on the protein or unaffected interactions. Such observations suggest that genes with mutations directly affecting protein binding properties are preferably located in central network positions and may influence critical nodes and edges in signal transduction networks.