MaveRegistry: a collaboration platform for multiplexed assays of variant effect

Abstract Summary Multiplexed assays of variant effect (MAVEs) are capable of experimentally testing all possible single nucleotide or amino acid variants in selected genomic regions, generating ‘variant effect maps’, which provide biochemical insight and functional evidence to enable more rapid and accurate clinical interpretation of human variation. Because the international community applying MAVE approaches is growing rapidly, we developed the online MaveRegistry platform to catalyze collaboration, reduce redundant efforts, allow stakeholders to nominate targets and enable tracking and sharing of progress on ongoing MAVE projects. Availability and implementation MaveRegistry service: https://registry.varianteffect.org. MaveRegistry source code: https://github.com/kvnkuang/maveregistry-front-end.

AbstractThe effect of single nucleotide variants (SNVs) in coding and non-coding regions is of great interest in genetics. Although many computational methods aim to elucidate the effects of SNVs on cellular mechanisms, it is not straightforward to comprehensively cover different molecular effects. To address this we compiled and benchmarked sequence and structure-based variant effect predictors and we analyzed the impact of nearly all possible amino acid and nucleotide variants in the reference genomes of H. sapiens, S. cerevisiae and E. coli. Studied mechanisms include protein stability, interaction interfaces, post-translational modifications and transcription factor binding sites. We apply this resource to the study of natural and disease coding variants. We also show how variant effects can be aggregated to generate protein complex burden scores that uncover protein complex to phenotype associations based on a set of newly generated growth profiles of 93 sequenced S. cerevisiae strains in 43 conditions. This resource is available through mutfunc, a tool by which users can query precomputed predictions by providing amino acid or nucleotide-level variants.

AbstractSummaryquasitools is a collection of newly-developed, open-source tools for analyzing viral quasispcies data. The application suite includes tools with the ability to create consensus sequences, call nucleotide, codon, and amino acid variants, calculate the complexity of a quasispecies, and measure the genetic distance between two similar quasispecies. These tools may be run independently or in user-created workflows.AvailabilityThe quasitools suite is a freely available application licensed under the Apache License, Version 2.0. The source code, documentation, and file specifications are available at: https://phac-nml.github.io/quasitools/[email protected]

Previous studies have identified variants in genes encoding proteins associated with the degree of addiction, smoking onset, and cessation. We aimed to describe thirty-one single nucleotide polymorphisms (SNPs) in seven candidate genomic regions spanning six genes associated with tobacco-smoking in a cross-sectional study from two different interventions for quitting smoking: (1) thirty-eight smokers were recruited via multimedia to participate in e-Decídete! program (e-Dec) and (2) ninety-four attended an institutional smoking cessation program on-site. SNPs genotyping was done by real-time PCR using TaqMan probes. The analysis of alleles and genotypes was carried out using the EpiInfo v7. on-site subjects had more years smoking and tobacco index than e-Dec smokers (p < 0.05, both); in CYP2A6 we found differences in the rs28399433 (p < 0.01), the e-Dec group had a higher frequency of TT genotype (0.78 vs. 0.35), and TG genotype frequency was higher in the on-site group (0.63 vs. 0.18), same as GG genotype (0.03 vs. 0.02). Moreover, three SNPs in NRXN1, two in CHRNA3, and two in CHRNA5 had differences in genotype frequencies (p < 0.01). Cigarettes per day were different (p < 0.05) in the metabolizer classification by CYP2A6 alleles. In conclusion, subjects attending a mobile smoking cessation intervention smoked fewer cigarettes per day, by fewer years, and by fewer cumulative pack-years. There were differences in the genotype frequencies of SNPs in genes related to nicotine metabolism and nicotine dependence. Slow metabolizers smoked more cigarettes per day than intermediate and normal metabolizers.

Mountain plants, challenged by vegetation time contractions and dynamic changes in environmental conditions, developed adaptations that help them to balance their growth, reproduction, survival, and regeneration. However, knowledge regarding the genetic basis of species adaptation to higher altitudes remain scarce for most plant species. Here, we attempted to identify such corresponding genomic regions of high evolutionary importance in two closely related European pines, Pinus mugo and P. uncinata, contrasting them with a reference lowland relative—P. sylvestris. We genotyped 438 samples at thousands of single nucleotide polymorphism (SNP) markers, tested their genetic differentiation and population structure followed by outlier detection and gene ontology annotations. Markers clearly differentiated the species and uncovered patterns of population structure in two of them. In P. uncinata three Pyrenean sites were grouped together, while two outlying populations constituted a separate cluster. In P. sylvestris, Spanish population appeared distinct from the remaining four European sites. Between mountain pines and the reference species, 35 candidate genes for altitude-dependent selection were identified, including such encoding proteins responsible for photosynthesis, photorespiration and cell redox homeostasis, regulation of transcription, and mRNA processing. In comparison between two mountain pines, 75 outlier SNPs were found in proteins involved mainly in the gene expression and metabolism.

AbstractWe introduce a new statistical model for patterns of linkage disequilibrium (LD) among multiple SNPs in a population sample. The model overcomes limitations of existing approaches to understanding, summarizing, and interpreting LD by (i) relating patterns of LD directly to the underlying recombination process; (ii) considering all loci simultaneously, rather than pairwise; (iii) avoiding the assumption that LD necessarily has a “block-like” structure; and (iv) being computationally tractable for huge genomic regions (up to complete chromosomes). We examine in detail one natural application of the model: estimation of underlying recombination rates from population data. Using simulation, we show that in the case where recombination is assumed constant across the region of interest, recombination rate estimates based on our model are competitive with the very best of current available methods. More importantly, we demonstrate, on real and simulated data, the potential of the model to help identify and quantify fine-scale variation in recombination rate from population data. We also outline how the model could be useful in other contexts, such as in the development of more efficient haplotype-based methods for LD mapping.

A highly conserved arginine residue is close to the catalytic center of PPM/PP2C-type protein phosphatases. Different crystal structures of PPM/PP2C homologues revealed that the guanidinium side chain of this arginine residue can adopt variable conformations and may bind ligands, suggesting an important role of this residue during catalysis. In this paper, we randomly mutated Arginine 13 of tPphA, a PPM/PP2C-type phosphatase from Thermosynechococcus elongatus, and obtained 18 different amino acid variants. The generated variants were tested towards p-nitrophenyl phosphate and various phosphopeptides. Towards p-nitrophenyl phosphate as substrate, twelve variants showed 3–7 times higher Km values than wild-type tPphA and four variants (R13D, R13F, R13L, and R13W) completely lost activity. Strikingly, these variants were still able to dephosphorylate phosphopeptides, although with strongly reduced activity. The specific inability of some Arg-13 variants to hydrolyze p-nitrophenyl phosphate highlights the importance of additional substrate interactions apart from the substrate phosphate for catalysis. The properties of the R13 variants indicate that this residue assists in substrate binding.

The common approaches for haplotype inference from genotype data are targeted toward phasing short genomic regions. Longer regions are often tackled in a heuristic manner, due to the high computational cost. Here, we describe a novel approach for phasing genotypes over long regions, which is based on combining information from local predictions on short, overlapping regions. The phasing is done in a way, which maximizes a natural maximum likelihood criterion. Among other things, this criterion takes into account the physical length between neighboring single nucleotide polymorphisms. The approach is very efficient and is applied to several large scale datasets and is shown to be successful in two recent benchmarking studies (Zaitlen et al., in press; Marchini et al., in preparation). Our method is publicly available via a webserver at .

AbstractThere have been five waves of influenza A (H7N9) epidemics in Zhejiang Province between 2013 and 2017. Although the epidemiological characteristics of the five waves have been reported, the molecular genetics aspects, including the phylogeny, evolution, and mutation of hemagglutinin (HA), have not been systematically investigated. A total of 154 H7N9 samples from Zhejiang Province were collected between 2013 and 2017 and sequenced using an Ion Torrent Personal Genome Machine. The starting dates of the waves were 16 March 2013, 1 July 2013, 1 July 2014, 1 July 2015, and 1 July 2016. Single-nucleotide polymorphisms (SNPs) and amino acid mutations were counted after the HA sequences were aligned. The evolution of H7N9 matched the temporal order of the five waves, among which wave 3 played an important role. The 55 SNPs and 14 amino acid mutations with high frequency identified among the five waves revealed the dynamic occurrence of mutation in the process of viral dissemination. Wave 3 contributed greatly to the subsequent epidemic of waves 4 and 5 of H7N9. Compared with wave 1, wave 5 was characterized by more mutations, including A143V and R148K, two mutations that have been reported to weaken the immune response. In addition, some amino acid mutations were observed in wave 5 that led to more lineages. It is necessary to strengthen the surveillance of subsequent H7N9 influenza outbreaks.