Protein Sequence Variation Research Articles

Campylobacter coli Clade 3 Isolates Induce Rapid Cell Death In Vitro.

Campylobacter bacteria are major human enteropathogens. Campylobacter coli shows less genetic diversity than C. jejuni and clusters into three clades, of which clade 1 includes most human and farm animal isolates, while environmental C. coli isolates mainly belong to clades 2 and 3. Recently, we sequenced the whole genomes of eight C. coli clade 2 and 3 isolates cultivated from water, and here we studied their interaction with human HT-29 colon cancer cells compared to that of clinical clade 1 isolates. All C. coli clade 3 isolates already caused cell necrosis 1 to 2 h after inoculation, whereas none of the clade 1 and 2 isolates analyzed induced cell death. Isolates from clades 2 and 3 adhered to epithelial cells better than clade 1 isolates, but all isolates induced similar levels of interleukin-8 (IL-8). Alignment and phylogenetic analysis of the translated putative virulence genes cadF, flpA, iamA, ciaB, and ceuE revealed clade-specific protein sequence variations, with clade 1 and 2 sequences being more closely related and clade 3 sequences being further apart, in general. Moreover, when RNA levels were measured, clade 3 isolates showed significantly lower levels of expression of cadF, iamA, and ceuE than clade 2 isolates, while flpA expression levels were higher in clade 3 isolates. The cytolethal distending toxin genes were also expressed in clades 2 and 3, although there was no difference between clades. Our findings demonstrate differences between the effects of C. coli clade 1, 2, and 3 isolates on human cells and suggest that C. coli clade 3 might be more virulent than clade 2 due to the observed cytotoxicity.IMPORTANCECampylobacter coli is a common zoonotic cause of gastroenteritis in humans worldwide. The majority of infections are caused by C. coli clade 1 isolates, whereas infections due to clade 2 and 3 isolates are rare. Whether this depends on a low prevalence of clade 2 and 3 isolates in reservoirs important for human infections or their lower ability to cause human disease is unknown. Here, we studied the effects of C. coli clade 2 and 3 isolates on a human cell line. These isolates adhered to human cells to a higher degree than clinical clade 1 isolates. Furthermore, we could show that C. coli clade 3 isolates rapidly induced cell death, suggesting differences in the virulence of C. coli The exact mechanism of cell death remains to be revealed, but selected genes showed interesting clade-specific expression patterns.

Contrasting patterns of molecular evolution in metazoan germ line genes

BackgroundGerm lines are the cell lineages that give rise to the sperm and eggs in animals. The germ lines first arise from primordial germ cells (PGCs) during embryogenesis: these form from either a presumed derived mode of preformed germ plasm (inheritance) or from an ancestral mechanism of inductive cell-cell signalling (induction). Numerous genes involved in germ line specification and development have been identified and functionally studied. However, little is known about the molecular evolutionary dynamics of germ line genes in metazoan model systems.ResultsHere, we studied the molecular evolution of germ line genes within three metazoan model systems. These include the genus Drosophila (N=34 genes, inheritance), the fellow insect Apis (N=30, induction), and their more distant relative Caenorhabditis (N=23, inheritance). Using multiple species and established phylogenies in each genus, we report that germ line genes exhibited marked variation in the constraint on protein sequence divergence (dN/dS) and codon usage bias (CUB) within each genus. Importantly, we found that de novo lineage-specific inheritance (LSI) genes in Drosophila (osk, pgc) and in Caenorhabditis (pie-1, pgl-1), which are essential to germ plasm functions under the derived inheritance mode, displayed rapid protein sequence divergence relative to the other germ line genes within each respective genus. We show this may reflect the evolution of specialized germ plasm functions and/or low pleiotropy of LSI genes, features not shared with other germ line genes. In addition, we observed that the relative ranking of dN/dS and of CUB between genera were each more strongly correlated between Drosophila and Caenorhabditis, from different phyla, than between Drosophila and its insect relative Apis, suggesting taxonomic differences in how germ line genes have evolved.ConclusionsTaken together, the present results advance our understanding of the evolution of animal germ line genes within three well-known metazoan models. Further, the findings provide insights to the molecular evolution of germ line genes with respect to LSI status, pleiotropy, adaptive evolution as well as PGC-specification mode.

Pnpla3 silencing with antisense oligonucleotides ameliorates nonalcoholic steatohepatitis and fibrosis in Pnpla3 I148M knock-in mice

ObjectiveNonalcoholic fatty liver disease (NAFLD) is becoming a leading cause of advanced chronic liver disease. The progression of NAFLD, including nonalcoholic steatohepatitis (NASH), has a strong genetic component, and the most robust contributor is the patatin-like phospholipase domain-containing 3 (PNPLA3) rs738409 encoding the 148M protein sequence variant. We hypothesized that suppressing the expression of the PNPLA3 148M mutant protein would exert a beneficial effect on the entire spectrum of NAFLD. MethodsWe examined the effects of liver-targeted GalNAc3-conjugated antisense oligonucleotide (ASO)-mediated silencing of Pnpla3 in a knock-in mouse model in which we introduced the human PNPLA3 I148M mutation. ResultsASO-mediated silencing of Pnpla3 reduced liver steatosis (p = 0.038) in homozygous Pnpla3 148M/M knock-in mutant mice but not in wild-type littermates fed a steatogenic high-sucrose diet. In mice fed a NASH-inducing diet, ASO-mediated silencing of Pnpla3 reduced liver steatosis score and NAFLD activity score independent of the Pnpla3 genotype, while reductions in liver inflammation score (p = 0.018) and fibrosis stage (p = 0.031) were observed only in the Pnpla3 knock-in 148M/M mutant mice. These responses were accompanied by reduced liver levels of Mcp1 (p = 0.026) and Timp2 (p = 0.007) specifically in the mutant knock-in mice. This may reduce levels of chemokine attracting inflammatory cells and increase the collagenolytic activity during tissue regeneration. ConclusionThis study provides the first evidence that a Pnpla3 ASO therapy can improve all features of NAFLD, including liver fibrosis, and suppress the expression of a strong innate genetic risk factor, Pnpla3 148M, which may open up a precision medicine approach in NASH.

Variation in the repulsive guidance molecule family in human populations.

Repulsive guidance molecules, RGMA, RGMB, and RGMC, are related proteins discovered independently through different experimental paradigms. They are encoded by single copy genes in mammalian and other vertebrate genomes, and are ~50% identical in amino acid sequence. The importance of RGM actions in human physiology has not been realized, as most research has focused on non‐human models, although mutations in RGMC are the cause of the severe iron storage disorder, juvenile hemochromatosis. Here I show that repositories of human genomic and population genetic data can be used as starting points for discovery and for developing new testable hypotheses about each of these paralogs in human biology and disease susceptibility. Information was extracted, aggregated, and analyzed from the Ensembl and UCSC Genome Browsers, the Exome Aggregation Consortium, the Genotype‐Tissue Expression project portal, the cBio portal for Cancer Genomics, and the National Cancer Institute Genomic Data Commons data site. Results identify extensive variation in gene expression patterns, substantial alternative RNA splicing, and possible missense alterations and other modifications in the coding regions of each of the three genes, with many putative mutations being detected in individuals with different types of cancers. Moreover, selected amino acid substitutions are highly prevalent in the world population, with minor allele frequencies of up to 37% for RGMA and up to 8% for RGMB. These results indicate that protein sequence variation is common in the human RGM family, and raises the possibility that individual variants will have a significant population impact on human physiology and/or disease predisposition.

Challenges of developing protein assays by mass spectrometry – method validation and quality issues

Protein assays by liquid chromatography mass spectrometry (LCMS) are gaining momentum as an alternative to immunoassays. The high selectivity of mass spectrometry allows for multiple proteins to be targeted in a single assay without cross reactivity and increases the breadth of data for the identification of clinically significant protein sequence variants and their post translational modifications (PTMs).

Designing antibody against highly conserved region of dengue envelope protein by in silico screening of scFv mutant library.

Dengue being one of the deadliest diseases of tropical regions, enforces to put continuous efforts for the development of vaccine and effective therapeutics. Most of the antibodies generated during dengue infection are non-neutralizing and cause antibody dependent enhancement. Hence, making a potent neutralizing antibody against all four dengue serotypes could be very effective for the treatment. However, designing a single antibody for all serotypes is difficult due to variation in protein sequences. Therefore, the objective is to identify conserved region of dengue envelope protein and then develop an antibody against that conserved region. Before advancing to the development of such an antibody, it is desirable to validate the interactions between antibody and dengue envelope protein. In silico analysis of such interactions provides a good platform to find out a suitable region to design and construct an antibody against it by analyzing antigen-antibody interaction before synthesizing the antibody. In this study, two highly conserved regions of dengue envelope protein were identified and an scFv was constructed against it. Both scFv and FuBc proteins were expressed in bacterial expression system and binding efficiency was analyzed by SPR analysis with KD value 2.3 μM. In order to improve binding efficiency, an in silico scFv mutant library was created which was virtually screened for higher binding efficiency. Six mutants with high binding efficiency were selected for further analysis. The binding ability of these mutants were predicted using simulation analysis which shows these mutations were stabilizing scFv-FuBc complex.

Crystal structure of a Thermus aquaticus diversity-generating retroelement variable protein.

Diversity-generating retroelements (DGRs) are widely distributed in bacteria, archaea, and microbial viruses, and bring about unparalleled levels of sequence variation in target proteins. While DGR variable proteins share low sequence identity, the structures of several such proteins have revealed the C-type lectin (CLec)-fold as a conserved scaffold for accommodating massive sequence variation. This conservation has led to the suggestion that the CLec-fold may be useful in molecular surface display applications. Thermostability is an attractive feature in such applications, and thus we studied the variable protein of a DGR encoded by a prophage of the thermophile Thermus aquaticus. We report here the 2.8 Å resolution crystal structure of the variable protein from the T. aquaticus DGR, called TaqVP, and confirm that it has a CLec-fold. Remarkably, its variable region is nearly identical in structure to those of several other CLec-fold DGR variable proteins despite low sequence identity among these. TaqVP was found to be thermostable, which appears to be a property shared by several CLec-fold DGR variable proteins. These results provide impetus for the pursuit of the DGR variable protein CLec-fold in molecular display applications.

An enhanced workflow for variant interpretation in UniProtKB/Swiss-Prot improves consistency and reuse in ClinVar

Personalized genomic medicine depends on integrated analyses that combine genetic and phenotypic data from individual patients with reference knowledge of the functional and clinical significance of sequence variants. Sources of this reference knowledge include the ClinVar repository of human genetic variants, a community resource that accepts submissions from external groups, and UniProtKB/Swiss-Prot, an expert-curated resource of protein sequences and functional annotation. UniProtKB/Swiss-Prot provides knowledge on the functional impact and clinical significance of over 30 000 human protein-coding sequence variants, curated from peer-reviewed literature reports. Here we present a pilot study that lays the groundwork for the integration of curated knowledge of protein sequence variation from UniProtKB/Swiss-Prot with ClinVar. We show that existing interpretations of variant pathogenicity in UniProtKB/Swiss-Prot and ClinVar are highly concordant, with 88% of variants that are common to the two resources having interpretations of clinical significance that agree. Re-curation of a subset of UniProtKB/Swiss-Prot variants according to American College of Medical Genetics and Genomics (ACMG) guidelines using ClinGen tools further increases this level of agreement, mainly due to the reclassification of supposedly pathogenic variants as benign, based on newly available population frequency data. We have now incorporated ACMG guidelines and ClinGen tools into the UniProt Knowledgebase (UniProtKB) curation workflow and routinely submit variant data from UniProtKB/Swiss-Prot to ClinVar. These efforts will increase the usability and utilization of UniProtKB variant data and will facilitate the continuing (re-)evaluation of clinical variant interpretations as data sets and knowledge evolve.

Regulation of organohalide respiration.

Organohalide respiration (OHR) is an anaerobic metabolism by which bacteria conserve energy with the use of halogenated compounds as terminal electron acceptors. Genes involved in OHR are organized in reductive dehalogenase (rdh) gene clusters and can be found in relatively high copy numbers in the genomes of organohalide-respiring bacteria (OHRB). The minimal rdh gene set is composed by rdhA and rdhB, encoding the catalytic enzyme involved in reductive dehalogenation and its putative membrane anchor, respectively. In this chapter, we present the major findings concerning the regulatory strategies developed by OHRB to control the expression of the rdh gene clusters. The first section focuses on the description of regulation patterns obtained from targeted transcriptional analyses, and from transcriptomic and proteomic studies, while the second section offers a detailed overview of the biochemically characterized OHR regulatory proteins identified so far. Depending on OHRB, transcriptional regulators belonging to three different protein families are found in the direct vicinity of rdh gene clusters, suggesting that they activate the transcription of their cognate gene cluster. In this chapter, strong emphasis was laid on the family of CRP/FNR-type RdhK regulators which belong to members of the genera Dehalobacter and Desulfitobacterium. Whereas only chlorophenols have been identified as effectors for RdhK regulators, the protein sequence diversity suggests a broader organohalide spectrum. Thus, effector identification of new regulators offers a promising alternative to elucidate the substrates of yet uncharacterized reductive dehalogenases. Future work investigating the possible cross-talk between OHR regulators and their possible use as biosensors is discussed.

In Vitro Neutralization Assay Using Cultured Hepatitis C Virus.

The method outlined here enables evaluation of the neutralization potency of monoclonal and polyclonal antibodies against in vitro cultured hepatitis C virus (HCV). The high variation in envelope protein sequence among HCV isolates necessitates the inclusion of several isolates, spanning the major genotypes of HCV, in order to make strong conclusions concerning the cross-reactive neutralization potential of a given antibody. This would be particularly relevant for any neutralization experiments aimed at uncovering novel therapeutic- or vaccine-relevant antibodies. In addition, these assays can also be used to compare neutralization sensitivity of novel cultured HCV to that of previously characterized isolates.

Keratin Gene Expression Differences in Wool Follicles and Sequence Diversity of High Glycine-Tyrosine Keratin-Associated Proteins (Kaps) in Magra Sheep of India

ABSTRACT Wool fibers consist of proteins with majority belongs to keratin family. Longitudinally arrayed keratin intermediate filaments constitute a matrix surrounded by keratin-associated proteins (KAPs). These KAPs play essential role in providing the physico-mechanical properties of the wool fiber. Real time PCR assay revealed that type 1 hair cortex keratin K33A (K1.2) was significantly (p ≤ 0.05) upregulated in lustrous Magra wool follicles however other keratin and keratin associated protein (KAP) genes were found down regulated. Furthermore, nucleotide sequence analysis revealed 4, 30, and 7 single nucleotide polymorphisms (SNPs) in KAP6, KAP7, and KAP8 genes, respectively. The presence of differing amount of type I and type II keratin and hydrophobic high glycine-tyrosine proteins in wool might be partially responsible for different glass transitions exhibited by wool and consequently imparts different physical properties like wool luster.

Scop3D: Online Visualization of Mutation Rates on Protein Structure.

Scop3D is a tool that automatically annotates protein structure with sequence conservation starting from a set of protein sequence variants. We present a complete upgrade and rewrite of Scop3D. We have included a DNA module that allows the analysis of single nucleotide polymorphisms in relation to the structural context of the protein. Scop3D therefore forms a bridge between genomics and protein structure. Moreover, Scop3D is now also available through an intuitive web-interface that makes the tool highly user-friendly.

PISTILLATA paralogs in Tarenaya hassleriana have diverged in interaction specificity

BackgroundFloral organs are specified by MADS-domain transcription factors that act in a combinatorial manner, as summarized in the (A)BCE model. However, this evolutionarily conserved model is in contrast to a remarkable amount of morphological diversity in flowers. One of the mechanisms suggested to contribute to this diversity is duplication of floral MADS-domain transcription factors. Although gene duplication is often followed by loss of one of the copies, sometimes both copies are retained. If both copies are retained they will initially be redundant, providing freedom for one of the paralogs to change function. Here, we examine the evolutionary fate and functional consequences of a transposition event at the base of the Brassicales that resulted in the duplication of the floral regulator PISTILLATA (PI), using Tarenaya hassleriana (Cleomaceae) as a model system.ResultsThe transposition of a genomic region containing a PI gene led to two paralogs which are located at different positions in the genome. The original PI copy is syntenic in position with most angiosperms, whereas the transposed copy is syntenic with the PI genes in Brassicaceae. The two PI paralogs of T. hassleriana have very similar expression patterns. However, they may have diverged in function, as only one of these PI proteins was able to act heterologously in the first whorl of A. thaliana flowers. We also observed differences in protein complex formation between the two paralogs, and the two paralogs exhibit subtle differences in DNA-binding specificity. Sequence analysis indicates that most of the protein sequence divergence between the two T. hassleriana paralogs emerged in a common ancestor of the Cleomaceae and the Brassicaceae.ConclusionsWe found that the PI paralogs in T. hassleriana have similar expression patterns, but may have diverged at the level of protein function. Data suggest that most protein sequence divergence occurred rapidly, prior to the origin of the Brassicaceae and Cleomaceae. It is tempting to speculate that the interaction specificities of the Brassicaceae-specific PI proteins are different compared to the PI found in other angiosperms. This could lead to PI regulating partly different genes in the Brassicaceae, and ultimately might result in change floral in morphology.

Toxin-antitoxin systems shows variability among Mycobacterium tuberculosis lineages.

The toxin-antitoxin (TA) systems are operons involved in the formation of persistent cells and in stress situations in microorganism. TA systems are widely distributed in Mycobacterium tuberculosis (MTB). The objective of this study was to determine the distribution and variability of protein sequences of TA systems in seven MTB lineages. Protein prediction on 73 genomes of different lineage was made using Prodigal, and profile hidden Markov models (PHMMs) of 225 reference TA proteins were constructed with HMMER. An homology search was made comparing the predicted proteins to PHMMs. Consecutive proteins that matched PHMMs (forming an operon) were selected. Thereafter, clustering of orthologous genes was made for further mutation scanning through multiple alignments. All proteins found belong to TA types II and IV, and 45 proteins were found completely conserved. Six uncharacterized antitoxins homologous to VapB11, VapB5, VapB45, VapB13, ParD1 and RelB were found. Multiple alignments revealed differences among lineages with specific mutations suitable for genotyping. Significant changes in amino acid sequences caused by frameshift mutations were found in 46 TA proteins.

Adaptive evolution of rbcL in Reaumuria soongarica (Tamaricaceae)

In the field of phylogenetic analyses, the rbcL gene encoded large subunit Ribulose-1, 5-biphosphate carboxylase/oxygenase (Rubisco, EC4.1.1.39), which plays a crucial role in the process of photosynthesis for most terrestrial plants, has been considered to be conserved; however, recent controversy regarding rbcL conservation has appeared since it was proposed to be under natural selection within all principal lineages of land plants. In this study, by examining the variation of DNA and protein sequences among 17 species in the family Tamaricaceae, three nonsynonymous mutations were identified to be under positive selection. The favored sites were located in the alph-helix domains of Rubisco, with decreased hydrophobicity and increased entropy, which could facilitate CO 2 penetration into the active site of Rubisco. We also found that the expression level of rbcL in different genotypes of Reaumuria soongarica shifted in response to various stresses such as drought, temperature, salt, and light. This study not only sheds light on the functional/structural features of Rubisco in the evolution scenarios from C 3 -like into C 4 in Tamaricaceae but also provides useful information on directing genetic performance to enhance photosynthesis efficiency of desert plants for sustaining fragile desert ecosystems; furthermore, it promotes the ability to cope with desert aridification and global warming.

ProteomeGenerator: A Framework for Comprehensive Proteomics Based on de Novo Transcriptome Assembly and High-Accuracy Peptide Mass Spectral Matching.

Modern mass spectrometry now permits genome-scale and quantitative measurements of biological proteomes. However, analysis of specific specimens is currently hindered by the incomplete representation of biological variability of protein sequences in canonical reference proteomes and the technical demands for their construction. Here, we report ProteomeGenerator, a framework for de novo and reference-assisted proteogenomic database construction and analysis based on sample-specific transcriptome sequencing and high-accuracy mass spectrometry proteomics. This enables the assembly of proteomes encoded by actively transcribed genes, including sample-specific protein isoforms resulting from non-canonical mRNA transcription, splicing, or editing. To improve the accuracy of protein isoform identification in non-canonical proteomes, ProteomeGenerator relies on statistical target-decoy database matching calibrated using sample-specific controls. Its current implementation includes automatic integration with MaxQuant mass spectrometry proteomics algorithms. We applied this method for the proteogenomic analysis of splicing factor SRSF2 mutant leukemia cells, demonstrating high-confidence identification of non-canonical protein isoforms arising from alternative transcriptional start sites, intron retention, and cryptic exon splicing as well as improved accuracy of genome-scale proteome discovery. Additionally, we report proteogenomic performance metrics for current state-of-the-art implementations of SEQUEST HT, MaxQuant, Byonic, and PEAKS mass spectral analysis algorithms. Finally, ProteomeGenerator is implemented as a Snakemake workflow within a Singularity container for one-step installation in diverse computing environments, thereby enabling open, scalable, and facile discovery of sample-specific, non-canonical, and neomorphic biological proteomes.

Haplosaurus computes protein haplotypes for use in precision drug design

Selecting the most appropriate protein sequences is critical for precision drug design. Here we describe Haplosaurus, a bioinformatic tool for computation of protein haplotypes. Haplosaurus computes protein haplotypes from pre-existing chromosomally-phased genomic variation data. Integration into the Ensembl resource provides rapid and detailed protein haplotypes retrieval. Using Haplosaurus, we build a database of unique protein haplotypes from the 1000 Genomes dataset reflecting real-world protein sequence variability and their prevalence. For one in seven genes, their most common protein haplotype differs from the reference sequence and a similar number differs on their most common haplotype between human populations. Three case studies show how knowledge of the range of commonly encountered protein forms predicted in populations leads to insights into therapeutic efficacy. Haplosaurus and its associated database is expected to find broad applications in many disciplines using protein sequences and particularly impactful for therapeutics design.

The core promoter controls basal and inducible expression of duck retinoic acid inducible gene-I (RIG-I)

Retinoic acid inducible gene-I (RIG-I) is a cytoplasmic RNA sensor for detecting a variety of RNA viruses including influenza A viruses. Detection ultimately produces Type I interferon (IFN), which stimulates expression of interferon stimulated genes (ISGs), including RIG-I itself in a positive feedback loop. The structure and function of RIG-I is conserved across phylogeny, despite significant protein sequence divergence, however, the promoter sequences do not show the expected phylogenetic relationships and it is not known whether they are similarly regulated. We previously cloned duck RIG-I and showed it is highly induced during influenza A infection consistent with induction by the interferon produced. Here, we identified the Pekin duck RIG-I promoter and constructed promoter reporter vectors, which we transfected into duck embryonic fibroblasts or chicken DF-1 cells and tested in dual luciferase assays. We showed that activation of the Mitochondrial Antiviral Signalling (MAVS) pathway using the constitutively active N-terminal region of RIG-I or polyinosinic-polycytidylic acid (poly I:C) led to stimulation of duck RIG-I promoter activity. Using deletion constructs we showed the core promoter lies in the proximal 250 basepairs, and we identified essential cis-regulatory elements, a GC-box and an interferon-sensitive response element (ISRE), responsible for basal and inducible expression, respectively. Using mCherry-tagged interferon regulatory factors (IRFs) cloned from chickens and ducks, we show overexpression of chIRF7 induced the duck RIG-I promoter, and this required the ISRE site. Finally, we also demonstrated that overexpressed chIRF7 translocated to the nucleus, which was augmented by MAVS activation using RIG-I 2CARD. Our findings demonstrate that RIG-I expression is induced by chIRF7, in a positive regulatory loop. These studies show that the duck RIG-I promoter is appropriately regulated in chicken cells, necessary for the potential generation of transgenic chickens expressing RIG-I.

RNA editing derived epitopes function as cancer antigens to elicit immune responses

In addition to genomic mutations, RNA editing is another major mechanism creating sequence variations in proteins by introducing nucleotide changes in mRNA sequences. Deregulated RNA editing contributes to different types of human diseases, including cancers. Here we report that peptides generated as a consequence of RNA editing are indeed naturally presented by human leukocyte antigen (HLA) molecules. We provide evidence that effector CD8+ T cells specific for edited peptides derived from cyclin I are present in human tumours and attack tumour cells that are presenting these epitopes. We show that subpopulations of cancer patients have increased peptide levels and that levels of edited RNA correlate with peptide copy numbers. These findings demonstrate that RNA editing extends the classes of HLA presented self-antigens and that these antigens can be recognised by the immune system.

Hgvs: A Python package for manipulating sequence variants using HGVS nomenclature: 2018 Update.

The Human Genome Variation Society (HGVS) nomenclature guidelines encourage the accurate and standard description of DNA, RNA, and protein sequence variants in public variant databases and the scientific literature. Inconsistent application of the HGVS guidelines can lead to misinterpretation of variants in clinical settings. Reliable software tools are essential to ensure consistent application of the HGVS guidelines when reporting and interpreting variants. We present the hgvs Python package, a comprehensive tool for manipulating sequence variants according to the HGVS nomenclature guidelines. Distinguishing features of the hgvs package include: (1) parsing, formatting, validating, and normalizing variants on genome, transcript, and protein sequences; (2) projecting variants between aligned sequences, including those with gapped alignments; (3) flexible installation using remote or local data (fully local installations eliminate network dependencies); (4) extensive automated tests; and (5) open source development by a community from eight organizations worldwide. This report summarizes recent and significant updates to the hgvs package since its original release in 2014, and presents results of extensive validation using clinical relevant variants from ClinVar and HGMD.