Deleterious Substitutions Research Articles

Molecular Analysis of Factor VIII and Factor IX Genes in Hemophilia Patients: Identification of Novel Mutations and Molecular Dynamics Studies.

BackgroundHemophilias A and B are X-linked bleeding disorders caused by mutations in the factor VIII and factor IX genes, respectively. Our objective was to identify the spectrum of mutations of the factor VIII and factor IX genes in Saudi Arabian population and determine the genotype and phenotype correlations by molecular dynamics (MD) simulation.MethodsFor genotyping, blood samples from Saudi Arabian patients were collected, and the genomic DNA was amplified, and then sequenced by Sanger method. For molecular simulations, we have used softwares such as CHARMM (Chemistry at Harvard Macromolecular Mechanics; http://www.charmm-gui.org) and GROMACS. In addition, the secondary structure was determined based on the solvent accessibility for the confirmation of the protein stability at the site of mutation.ResultsSix mutations (three novel and three known) were identified in factor VIII gene, and six mutations (one novel and five known) were identified in factor IX gene. The factor VIII novel mutations identified were c.99G>T, p. (W33C) in exon 1, c.2138 DelA, p. (N713Tfs*9) in eon14, also a novel mutation at splicing acceptor site of exon 23 c.6430 - 1G>A. In factor IX, we found a novel mutation c.855G>C, p. (E285D) in exon 8. These novel mutations were not reported in any factor VIII or factor IX databases previously. The deleterious effects of these novel mutations were confirmed by PolyPhen2 and SIFT programs.ConclusionThe protein functional and structural studies and the models built in this work would be appropriate for predicting the effects of deleterious amino acid substitutions causing these genetic disorders. These findings are useful for genetic counseling in the case of consanguineous marriages which is more common in the Saudi Arabia.

Uniparental Inheritance Promotes Adaptive Evolution in Cytoplasmic Genomes.

Eukaryotes carry numerous asexual cytoplasmic genomes (mitochondria and plastids). Lacking recombination, asexual genomes should theoretically suffer from impaired adaptive evolution. Yet, empirical evidence indicates that cytoplasmic genomes experience higher levels of adaptive evolution than predicted by theory. In this study, we use a computational model to show that the unique biology of cytoplasmic genomes—specifically their organization into host cells and their uniparental (maternal) inheritance—enable them to undergo effective adaptive evolution. Uniparental inheritance of cytoplasmic genomes decreases competition between different beneficial substitutions (clonal interference), promoting the accumulation of beneficial substitutions. Uniparental inheritance also facilitates selection against deleterious cytoplasmic substitutions, slowing Muller’s ratchet. In addition, uniparental inheritance generally reduces genetic hitchhiking of deleterious substitutions during selective sweeps. Overall, uniparental inheritance promotes adaptive evolution by increasing the level of beneficial substitutions relative to deleterious substitutions. When we assume that cytoplasmic genome inheritance is biparental, decreasing the number of genomes transmitted during gametogenesis (bottleneck) aids adaptive evolution. Nevertheless, adaptive evolution is always more efficient when inheritance is uniparental. Our findings explain empirical observations that cytoplasmic genomes—despite their asexual mode of reproduction—can readily undergo adaptive evolution.

Computational prediction of the effects of non-synonymous single nucleotide polymorphisms in the human Quinone Oxidoreductase 1 (NQO1)

The human NAD(P)H Quinone Oxidoreductase 1 (NQO1) gene codes for a flavoenzyme which protects cells from oxidative stress and also regulates ubiquitin-independent proteasomal degradation of several proteins including p53. Polymorphisms in NQO1 are associated with several diseases but of the thousand SNPs reported till date, only a few elevate risk. Thus during an epidemiological investigation, studying a few selected SNPs could be more justified. Computational approaches help in the identification of non-synonymous single nucleotide polymorphisms (nsSNPs) that affect protein functions and impact disease risk. The present study evaluates the structure and functional effect of all known nsSNPs of NQO1. NQO1 SNP information and their corresponding protein sequences were collected from dbSNP and UniProt databases respectively. Genomic analysis of the 75 nsSNPs was initiated using SIFT, Polyphen, SNP & GO, PROVEAN, MutPred, I-Mutant 2.0 and nsSNPAnalyzer. Protein structural analysis of select amino acid variants was performed using Phyre2, RaptorX and Swiss-PDB Viewer, molecular dynamics using GROMOS96 and root mean square deviation (RMSD) values were calculated using UCSF Chimera. We identified eight deleterious substitutions - rs11555215 (R119P), rs368942932 (L7R), rs554822773 (L113P), rs749792989 (W116L), rs764253929 (G160R), rs769655604 (W106R), rs775472154 (F102V) and rs777649241 (D41G). The NQO1 variants -R119P and G160R, were the most energetically unfavorable substitutions and could destabilize the amino acid interactions and hydrogen bond networks. The study concludes that screening for these variants especially rs11555215 and rs764253929 may be useful in understanding disease pathogenesis and developing molecular diagnostics involving the NQO1 pathways.

De novo, systemic, deleterious amino acid substitutions are common in large cytoskeleton-related protein coding regions.

Human mutagenesis is largely random, thus large coding regions, simply on the basis of probability, represent relatively large mutagenesis targets. Thus, we considered the possibility that large cytoskeletal-protein related coding regions (CPCRs), including extra-cellular matrix (ECM) coding regions, would have systemic nucleotide variants that are not present in common SNP databases. Presumably, such variants arose recently in development or in recent, preceding generations. Using matched breast cancer and blood-derived normal datasets from the cancer genome atlas, CPCR single nucleotide variants (SNVs) not present in the All SNPs(142) or 1000 Genomes databases were identified. Using the Protein Variation Effect Analyzer internet-based tool, it was discovered that apparent, systemic mutations (not shared among others in the analysis group) in the CPCRs, represented numerous deleterious amino acid substitutions. However, no such deleterious variants were identified among the (cancer blood-matched) variants shared by other members of the analysis group. These data indicate that private SNVs, which potentially have a medical consequence, occur de novo with significant frequency in the larger, human coding regions that collectively impact the cytoskeleton and ECM.

Extreme genomic erosion after recurrent demographic bottlenecks in the highly endangered Iberian lynx.

BackgroundGenomic studies of endangered species provide insights into their evolution and demographic history, reveal patterns of genomic erosion that might limit their viability, and offer tools for their effective conservation. The Iberian lynx (Lynx pardinus) is the most endangered felid and a unique example of a species on the brink of extinction.ResultsWe generate the first annotated draft of the Iberian lynx genome and carry out genome-based analyses of lynx demography, evolution, and population genetics. We identify a series of severe population bottlenecks in the history of the Iberian lynx that predate its known demographic decline during the 20th century and have greatly impacted its genome evolution. We observe drastically reduced rates of weak-to-strong substitutions associated with GC-biased gene conversion and increased rates of fixation of transposable elements. We also find multiple signatures of genetic erosion in the two remnant Iberian lynx populations, including a high frequency of potentially deleterious variants and substitutions, as well as the lowest genome-wide genetic diversity reported so far in any species.ConclusionsThe genomic features observed in the Iberian lynx genome may hamper short- and long-term viability through reduced fitness and adaptive potential. The knowledge and resources developed in this study will boost the research on felid evolution and conservation genomics and will benefit the ongoing conservation and management of this emblematic species.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-016-1090-1) contains supplementary material, which is available to authorized users.

Evolutionary switches between two serine codon sets are driven by selection

Serine is the only amino acid that is encoded by two disjoint codon sets so that a tandem substitution of two nucleotides is required to switch between the two sets. Previously published evidence suggests that, for the most evolutionarily conserved serines, the codon set switch occurs by simultaneous substitution of two nucleotides. Here we report a genome-wide reconstruction of the evolution of serine codons in triplets of closely related species from diverse prokaryotes and eukaryotes. The results indicate that the great majority of codon set switches proceed by two consecutive nucleotide substitutions, via a threonine or cysteine intermediate, and are driven by selection. These findings imply a strong pressure of purifying selection in protein evolution, which in the case of serine codon set switches occurs via an initial deleterious substitution quickly followed by a second, compensatory substitution. The result is frequent reversal of amino acid replacements and, at short evolutionary distances, pervasive homoplasy.

Extensively Parameterized Mutation-Selection Models Reliably Capture Site-Specific Selective Constraint.

The mutation-selection model of coding sequence evolution has received renewed attention for its use in estimating site-specific amino acid propensities and selection coefficient distributions. Two computationally tractable mutation-selection inference frameworks have been introduced: One framework employs a fixed-effects, highly parameterized maximum likelihood approach, whereas the other employs a random-effects Bayesian Dirichlet Process approach. While both implementations follow the same model, they appear to make distinct predictions about the distribution of selection coefficients. The fixed-effects framework estimates a large proportion of highly deleterious substitutions, whereas the random-effects framework estimates that all substitutions are either nearly neutral or weakly deleterious. It remains unknown, however, how accurately each method infers evolutionary constraints at individual sites. Indeed, selection coefficient distributions pool all site-specific inferences, thereby obscuring a precise assessment of site-specific estimates. Therefore, in this study, we use a simulation-based strategy to determine how accurately each approach recapitulates the selective constraint at individual sites. We find that the fixed-effects approach, despite its extensive parameterization, consistently and accurately estimates site-specific evolutionary constraint. By contrast, the random-effects Bayesian approach systematically underestimates the strength of natural selection, particularly for slowly evolving sites. We also find that, despite the strong differences between their inferred selection coefficient distributions, the fixed- and random-effects approaches yield surprisingly similar inferences of site-specific selective constraint. We conclude that the fixed-effects mutation-selection framework provides the more reliable software platform for model application and future development.

Three-Dimensional Model of Human Nicotinamide Nucleotide Transhydrogenase (NNT) and Sequence-Structure Analysis of its Disease-Causing Variations.

ABSTRACTDefective mitochondrial proteins are emerging as major contributors to human disease. Nicotinamide nucleotide transhydrogenase (NNT), a widely expressed mitochondrial protein, has a crucial role in the defence against oxidative stress. NNT variations have recently been reported in patients with familial glucocorticoid deficiency (FGD) and in patients with heart failure. Moreover, knockout animal models suggest that NNT has a major role in diabetes mellitus and obesity. In this study, we used experimental structures of bacterial transhydrogenases to generate a structural model of human NNT (H‐NNT). Structure‐based analysis allowed the identification of H‐NNT residues forming the NAD binding site, the proton canal and the large interaction site on the H‐NNT dimer. In addition, we were able to identify key motifs that allow conformational changes adopted by domain III in relation to its functional status, such as the flexible linker between domains II and III and the salt bridge formed by H‐NNT Arg882 and Asp830. Moreover, integration of sequence and structure data allowed us to study the structural and functional effect of deleterious amino acid substitutions causing FGD and left ventricular non‐compaction cardiomyopathy. In conclusion, interpretation of the function–structure relationship of H‐NNT contributes to our understanding of mitochondrial disorders.

Systematic TOR1A non-c.907_909delGAG variant analysis in isolated dystonia and controls

BackgroundAn increasing number of rare, functionally relevant non-c.907_909delGAG (non-ΔGAG) variants in TOR1A have been recognized, associated with phenotypic expressions different from classic DYT1 childhood-onset generalized dystonia. Only recently, DYT1 genotype-phenotype correlations have been proposed, awaiting further elucidation in independent cohorts. MethodsWe screened the entire coding sequence and the 5′-UTR region of TOR1A for rare non-ΔGAG sequence variants in a large series of 940 individuals with various forms of isolated dystonia as well as in 376 ancestry-matched controls. The frequency of rare, predicted deleterious non-ΔGAG TOR1A variants was assessed in the European sample of the Exome Aggregation Consortium (ExAC) dataset. ResultsIn the case cohort, we identified a rare 5′-UTR variant (c.-39G > T), a rare splice-region variant (c.445-8T > C), as well as one novel (p.Ile231Asn) and two rare (p.Ala163Val, p.Thr321Met) missense variants, each in a single patient with adult-onset focal/segmental isolated dystonia. Of these variants, only p.Thr321Met qualified as possibly disease-related according to variant interpretation criteria. One novel, predicted deleterious missense substitution (p.Asn208Ser) was detected in the control cohort. Among European ExAC individuals, the carrier rate of rare, predicted deleterious non-ΔGAG variants was 0.4%. ConclusionsOur study does not allow the establishment of genotype-specific clinical correlations for DYT1. Further large-scale genetic screening accompanied by comprehensive segregation and functional studies is required to conclusively define the contribution of TOR1A whole-gene variation to the pathogenesis of isolated dystonia.

Cystinuria Associated with Different SLC7A9 Gene Variants in the Cat.

Cystinuria is a classical inborn error of metabolism characterized by a selective proximal renal tubular defect affecting cystine, ornithine, lysine, and arginine (COLA) reabsorption, which can lead to uroliths and urinary obstruction. In humans, dogs and mice, cystinuria is caused by variants in one of two genes, SLC3A1 and SLC7A9, which encode the rBAT and bo,+AT subunits of the bo,+ basic amino acid transporter system, respectively. In this study, exons and flanking regions of the SLC3A1 and SLC7A9 genes were sequenced from genomic DNA of cats (Felis catus) with COLAuria and cystine calculi. Relative to the Felis catus-6.2 reference genome sequence, DNA sequences from these affected cats revealed 3 unique homozygous SLC7A9 missense variants: one in exon 5 (p.Asp236Asn) from a non-purpose-bred medium-haired cat, one in exon 7 (p.Val294Glu) in a Maine Coon and a Sphinx cat, and one in exon 10 (p.Thr392Met) from a non-purpose-bred long-haired cat. A genotyping assay subsequently identified another cystinuric domestic medium-haired cat that was homozygous for the variant originally identified in the purebred cats. These missense variants result in deleterious amino acid substitutions of highly conserved residues in the bo,+AT protein. A limited population survey supported that the variants found were likely causative. The remaining 2 sequenced domestic short-haired cats had a heterozygous variant at a splice donor site in intron 10 and a homozygous single nucleotide variant at a branchpoint in intron 11 of SLC7A9, respectively. This study identifies the first SLC7A9 variants causing feline cystinuria and reveals that, as in humans and dogs, this disease is genetically heterogeneous in cats.

The Role of Deleterious Substitutions in Crop Genomes.

Populations continually incur new mutations with fitness effects ranging from lethal to adaptive. While the distribution of fitness effects of new mutations is not directly observable, many mutations likely either have no effect on organismal fitness or are deleterious. Historically, it has been hypothesized that a population may carry many mildly deleterious variants as segregating variation, which reduces the mean absolute fitness of the population. Recent advances in sequencing technology and sequence conservation-based metrics for inferring the functional effect of a variant permit examination of the persistence of deleterious variants in populations. The issue of segregating deleterious variation is particularly important for crop improvement, because the demographic history of domestication and breeding allows deleterious variants to persist and reach moderate frequency, potentially reducing crop productivity. In this study, we use exome resequencing of 15 barley accessions and genome resequencing of 8 soybean accessions to investigate the prevalence of deleterious single nucleotide polymorphisms (SNPs) in the protein-coding regions of the genomes of two crops. We conclude that individual cultivars carry hundreds of deleterious SNPs on average, and that nonsense variants make up a minority of deleterious SNPs. Our approach identifies known phenotype-altering variants as deleterious more frequently than the genome-wide average, suggesting that putatively deleterious variants are likely to affect phenotypic variation. We also report the implementation of a SNP annotation tool BAD_Mutations that makes use of a likelihood ratio test based on alignment of all currently publicly available Angiosperm genomes.

Evolutionary Consequences of DNA Methylation in a Basal Metazoan.

Gene body methylation (gbM) is an ancestral and widespread feature in Eukarya, yet its adaptive value and evolutionary implications remain unresolved. The occurrence of gbM within protein-coding sequences is particularly puzzling, because methylation causes cytosine hypermutability and hence is likely to produce deleterious amino acid substitutions. We investigate this enigma using an evolutionarily basal group of Metazoa, the stony corals (order Scleractinia, class Anthozoa, phylum Cnidaria). We show that patterns of coral gbM are similar to other invertebrate species, predicting wide and active transcription and slower sequence evolution. We also find a strong correlation between gbM and codon bias, resulting from systematic replacement of CpG bearing codons. We conclude that gbM has strong effects on codon evolution and speculate that this may influence establishment of optimal codons.

Coupling high-throughput genetics with phylogenetic information reveals an epistatic interaction on the influenza A virus M segment.

BackgroundEpistasis is one of the central themes in viral evolution due to its importance in drug resistance, immune escape, and interspecies transmission. However, there is a lack of experimental approach to systematically probe for epistatic residues.ResultsBy utilizing the information from natural occurring sequences and high-throughput genetics, this study established a novel strategy to identify epistatic residues. The rationale is that a substitution that is deleterious in one strain may be prevalent in nature due to the presence of a naturally occurring compensatory substitution. Here, high-throughput genetics was applied to influenza A virus M segment to systematically identify deleterious substitutions. Comparison with natural sequence variation showed that a deleterious substitution M1 Q214H was prevalent in circulating strains. A coevolution analysis was then performed and indicated that M1 residues 121, 207, 209, and 214 naturally coevolved as a group. Subsequently, we experimentally validated that M1 A209T was a compensatory substitution for M1 Q214H.ConclusionsThis work provided a proof-of-concept to identify epistatic residues by coupling high-throughput genetics with phylogenetic information. In particular, we were able to identify an epistatic interaction between M1 substitutions A209T and Q214H. This analytic strategy can potentially be adapted to study any protein of interest, provided that the information on natural sequence variants is available.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-2358-7) contains supplementary material, which is available to authorized users.

Apolipoprotein A-V gene therapy for disease prevention / treatment: a critical analysis.

Apolipoprotein (apo) A-V is a novel member of the class of exchangeable apo's involved in triacylglycerol (TG) homeostasis. Whereas a portion of hepatic-derived apoA-V is secreted into plasma and functions to facilitate lipo­protein lipase-mediated TG hydrolysis, another portion is recovered intracellularly, in association with cytosolic lipid droplets. Loss of apoA-V function is positively correlated with elevated plasma TG and increased risk of car­diovascular disease. Single nucleotide polymorphisms (SNP) in the APOA5 locus can affect transcription efficiency or introduce deleterious amino acid substitutions. Likewise, rare mutations in APOA5 that compromise functionality are associated with increased plasma TG and premature myocardial infarction. Genetically engineered mouse mod­els and human population studies suggest that, in certain instances, supplementation with wild type (WT) apoA-V may have therapeutic benefit. It is hypothesized that individuals that manifest elevated plasma TG owing to deleter­ious APOA5 SNPs or rare mutations would respond to WT apoA-V supplementation with improved plasma TG clearance. On the other hand, subjects with hypertriglyceridemia of independent origin (unrelated to apoA-V func­tion) may not respond to apoA-V augmentation in this manner. Improvement in the ability to identify individuals predicted to benefit, advances in gene transfer technology and the strong connection between HTG and heart disease, point to apoA-V supplementation as a viable disease prevention / therapeutic strategy. Candidates would include individuals that manifest chronic TG elevation, have low plasma apoA-V due to an APOA5 mutation/polymorphism and not have deleterious mutations/polymorphisms in other genes known to influence plasma TG levels.

Developmental, transcriptome, and genetic alterations associated with parthenocarpy in the grapevine seedless somatic variant Corinto bianco.

Seedlessness is a relevant trait in grapevine cultivars intended for fresh consumption or raisin production. Previous DNA marker analysis indicated that Corinto bianco (CB) is a parthenocarpic somatic variant of the seeded cultivar Pedro Ximenes (PX). This study compared both variant lines to determine the basis of this parthenocarpic phenotype. At maturity, CB seedless berries were 6-fold smaller than PX berries. The macrogametophyte was absent from CB ovules, and CB was also pollen sterile. Occasionally, one seed developed in 1.6% of CB berries. Microsatellite genotyping and flow cytometry analyses of seedlings generated from these seeds showed that most CB viable seeds were formed by fertilization of unreduced gametes generated by meiotic diplospory, a process that has not been described previously in grapevine. Microarray and RNA-sequencing analyses identified 1958 genes that were differentially expressed between CB and PX developing flowers. Genes downregulated in CB were enriched in gametophyte-preferentially expressed transcripts, indicating the absence of regular post-meiotic germline development in CB. RNA-sequencing was also used for genetic variant calling and 14 single-nucleotide polymorphisms distinguishing the CB and PX variant lines were detected. Among these, CB-specific polymorphisms were considered as candidate parthenocarpy-responsible mutations, including a putative deleterious substitution in a HAL2-like protein. Collectively, these results revealed that the absence of a mature macrogametophyte, probably due to meiosis arrest, coupled with a process of fertilization-independent fruit growth, caused parthenocarpy in CB. This study provides a number of grapevine parthenocarpy-responsible candidate genes and shows how genomic approaches can shed light on the genetic origin of woody crop somatic variants.

The Roles of Compensatory Evolution and Constraint in Aminoacyl tRNA Synthetase Evolution.

Mitochondrial protein translation requires interactions between transfer RNAs encoded by the mitochondrial genome (mt-tRNAs) and mitochondrial aminoacyl tRNA synthetase proteins (mt-aaRS) encoded by the nuclear genome. It has been argued that animal mt-tRNAs have higher deleterious substitution rates relative to their nuclear-encoded counterparts, the cytoplasmic tRNAs (cyt-tRNAs). This dynamic predicts elevated rates of compensatory evolution of mt-aaRS that interact with mt-tRNAs, relative to aaRS that interact with cyt-tRNAs (cyt-aaRS). We find that mt-aaRS do evolve at significantly higher rates (exemplified by higher dN and dN/dS) relative to cyt-aaRS, across mammals, birds, and Drosophila. While this pattern supports a model of compensatory evolution, the level at which a gene is expressed is a more general predictor of protein evolutionary rate. We find that gene expression level explains 10–56% of the variance in aaRS dN/dS, and that cyt-aaRS are more highly expressed in addition to having lower dN/dS values relative to mt-aaRS, consistent with more highly expressed genes being more evolutionarily constrained. Furthermore, we find no evidence of positive selection acting on either class of aaRS protein, as would be expected under a model of compensatory evolution. Nevertheless, the signature of faster mt-aaRS evolution persists in mammalian, but not bird or Drosophila, lineages after controlling for gene expression, suggesting some additional effect of compensatory evolution for mammalian mt-aaRS. We conclude that gene expression is the strongest factor governing differential amino acid substitution rates in proteins interacting with mitochondrial versus cytoplasmic factors, with important differences in mt-aaRS molecular evolution among taxonomic groups.

A Missense Mutation in SLC45A2 Is Associated with Albinism in Several Small Long Haired Dog Breeds.

Homozygosity for a large deletion in the solute carrier family 45, member 2 (SLC45A2) gene causes oculocutaneous albinism (OCA) in the Doberman Pinscher breed. An albino Lhasa Apso did not have this g.27141_31223del (CanFam2) deletion in her SLC45A2 sequence. Therefore, SLC45A2 was investigated in this female Lhasa Apso to search for other possible variants that caused her albinism. The albino Lhasa Apso was homozygous for a nonsynonymous substitution in the seventh exon, a c.1478G>A base change that resulted in a glycine to aspartic acid substitution (p.G493D). This mutation was not found in a wolf, a coyote, or any of the 15 other Lhasa Apso dogs or 32 other dogs of breeds related to the Lhasa Apso. However, an albino Pekingese, 2 albino Pomeranians, and an albino mixed breed dog that was small and long haired were also homozygous for the 493D allele. The colored puppies of the albino Lhasa Apso and the colored dam of the 2 albino Pomeranians were heterozygous for this allele. However, an albino Pug was homozygous for the 493G allele and therefore although we suggest the 493D allele causes albinism when homozygous in several small, long haired dog breeds, it does not explain all albinism in dogs. A variant effect prediction for the albino Lhasa Apso confirms that p.G493D is a deleterious substitution, and a topology prediction for SLC45A2 suggests that the 11th transmembrane domain where the 493rd amino acid was located, has an altered structure.

In silico screening, genotyping, molecular dynamics simulation and activity studies of SNPs in pyruvate kinase M2.

Role of, 29-non-synonymous, 15-intronic, 3-close to UTR, single nucleotide polymorphisms (SNPs) and 2 mutations of Human Pyruvate Kinase (PK) M2 were investigated by in-silico and in-vitro functional studies. Prediction of deleterious substitutions based on sequence homology and structure based servers, SIFT, PANTHER, SNPs&GO, PhD-SNP, SNAP and PolyPhen, depicted that 19% emerged common between all the mentioned programs. SNPeffect and HOPE showed three substitutions (C31F, Q310P and S437Y) in-silico as deleterious and functionally important. In-vitro activity assays showed C31F and S437Y variants of PKM2 with reduced activity, while Q310P variant was catalytically inactive. The allosteric activation due to binding of fructose 1-6 bisphosphate (FBP) was compromised in case of S437Y nsSNP variant protein. This was corroborated through molecular dynamics (MD) simulation study, which was also carried out in other two variant proteins. The 5 intronic SNPs of PKM2, associated with sporadic breast cancer in a case-control study, when subjected to different computational analyses, indicated that 3 SNPs (rs2856929, rs8192381 and rs8192431) could generate an alternative transcript by influencing splicing factor binding to PKM2. We propose that these, potentially functional and important variations, both within exons and introns, could have a bearing on cancer metabolism, since PKM2 has been implicated in cancer in the recent past.

The Evolutionarily Stable Distribution of Fitness Effects

The distribution of fitness effects (DFE) of new mutations is a key parameter in determining the course of evolution. This fact has motivated extensive efforts to measure the DFE or to predict it from first principles. However, just as the DFE determines the course of evolution, the evolutionary process itself constrains the DFE. Here, we analyze a simple model of genome evolution in a constant environment in which natural selection drives the population toward a dynamic steady state where beneficial and deleterious substitutions balance. The distribution of fitness effects at this steady state is stable under further evolution and provides a natural null expectation for the DFE in a population that has evolved in a constant environment for a long time. We calculate how the shape of the evolutionarily stable DFE depends on the underlying population genetic parameters. We show that, in the absence of epistasis, the ratio of beneficial to deleterious mutations of a given fitness effect obeys a simple relationship independent of population genetic details. Finally, we analyze how the stable DFE changes in the presence of a simple form of diminishing-returns epistasis.

Carbapenem resistance in cystic fibrosis strains of Pseudomonas aeruginosa as a result of amino acid substitutions in porin OprD

The aim of this work was to investigate the impact of single amino acid substitutions occurring in specific porin OprD on carbapenem resistance of cystic fibrosis (CF) strains of Pseudomonas aeruginosa. A PAO1ΔoprD mutant was complemented with the oprD genes from five carbapenem-resistant CF strains exhibiting very low amounts of mutated OprD porins in their outer membrane despite wild-type levels of oprD transcripts. Compared with wild-type porin from strain PAO1, single amino acid substitutions S403P (in periplasmic loop 8), Y242H, S278P and L345P (in β-sheets 10, 12 and 14, respectively) were found to result in reduced amounts of OprD in the outer membrane, increased carbapenem resistance, and slower growth in minimal medium containing gluconate, an OprD substrate, as the sole source of carbon and energy. This indicates that in CF strains of P. aeruginosa, loss of porin OprD may not only result from mutations downregulating the expression of or disrupting the oprD gene, but also from mutations generating deleterious amino acid substitutions in the porin structure.