Prediction of deleterious non-synonymous single-nucleotide polymorphisms of human uridine diphosphate glucuronosyltransferase genes

Evolutionary history of MEK1 illuminates the nature of deleterious mutations

Mutations in signal transduction pathways lead to various diseases including cancers. MEK1 kinase, encoded by the human MAP2K1 gene, is one of the central components of the MAPK pathway and more than a hundred somatic mutations in the MAP2K1 gene were identified in various tumors. Germline mutations deregulating MEK1 also lead to congenital abnormalities, such as the cardiofaciocutaneous syndrome and arteriovenous malformation. Evaluating variants associated with a disease is a challenge, and computational genomic approaches aid in this process. Establishing evolutionary history of a gene improves computational prediction of disease-causing mutations; however, the evolutionary history of MEK1 is not well understood. Here, by revealing a precise evolutionary history of MEK1, we construct a well-defined dataset of MEK1 metazoan orthologs, which provides sufficient depth to distinguish between conserved and variable amino acid positions. We matched known and predicted disease-causing and benign mutations to evolutionary changes observed in corresponding amino acid positions and found that all known and many suspected disease-causing mutations are evolutionarily intolerable. We selected several variants that cannot be unambiguously assessed by automated prediction tools but that are confidently identified as "damaging" by our approach, for experimental validation in Drosophila. In all cases, evolutionary intolerant variants caused increased mortality and severe defects in fruit fly embryos confirming their damaging nature. We anticipate that our analysis will serve as a blueprint to help evaluate known and novel missense variants in MEK1 and that our approach will contribute to improving automated tools for disease-associated variant interpretation.

Concurrence of novel mutations causing Gilbert's and Dubin-Johnson syndrome with poor clinical outcomes in a Han Chinese family

Dual-hereditary jaundice (Dubin-Johnson syndrome (DJS) and Gilbert's syndrome (GS)) is a rare clinical entity resulting from defects of the ATP binding cassette subfamily C member 2 (ABCC2) and UDP glucuronosyltransferase family 1 member A1 (UGT1A1) genes with autosomal recessive inheritance. In this study, we aimed to investigate the mutation profiles and characterize the phenotypes in a Han Chinese family with DJS and GS. Genetic screening for variants in the ABCC2 and UGT1A1, immunohistochemistry for expression of ABCC2, and histopathological examination were carried out. The proband and his brother had unconjugated and conjugated hyperbilirubinemia after birth. The proband's sister had only conjugated hyperbilirubinemia after birth. The proband developed into pleural effusions and ascites, pericardial thickening, intrahepatic and extrahepatic biliary duct dilatation, and enlarged gallbladder at age 50. Hepatocellular carcinoma occurred in the proband's brother at age 46. Seven compound defects of the ABCC2 gene [c.2414delG, p.(Ile1489Gly), p.(Thr1490Pro), and p.(Ile1491Gln)] and the UGT1A1 gene (c.-3279T>G, p.(Gly71Arg), and p.(Pro451Leu)) were identified in family members. Accumulation of pigment in hepatocytes characteristic of that in DJS was present in the proband and his brother. Expression of ABCC2 protein was markedly diminished in the patient's liver. Our results show a different genetic profile of DJS and GS in a Han Chinese family, indicating a more complex pattern of dual-hereditary jaundice among different populations. The present study illuminates the underpinnings of DJS and GS and extends the mutation profiles and phenotypes of these two syndromes in dual-hereditary jaundice.

Assessment of susceptibility to phthalate and DINCH exposure through CYP and UGT single nucleotide polymorphisms

Single nucleotide polymorphisms (SNPs) of cytochrome P450 (CYPs) and UDP-glucuronosyltransferase (UGTs) genes have been proposed to influence phthalates and 1,2-cyclo-hexanedicarboxylic acid diisononyl ester (DINCH) biotransformation but have not been investigated on a populational level. We investigated the role of SNPs in CYP2C9, CYP2C19, CYP2D6, UGT2B15, and UGT1A7 genes in the biotransformation of phthalates (DEHP, DEP, DiBP, DnBP, BBzP, DiNP, DidP) and DINCH by determining their urine metabolites. From the Slovenian study population of 274 men and 289 lactating primiparous women we obtained data on phthalate and DINCH urine metabolite levels (MEHP, 5OH-MEHP, 5oxo-MEHP, 5cx-MEPP, MEP, MiBP, MnBP, MBzP, cx-MINP, OH-MiDP, MCHP, MnPeP, MnOP, 5OH-MINCH, 5oxo-MINCH), SNP genotypes (rs1057910 = CYP2C9*3, rs1799853 = CYP2C9*2, rs4244285 = CYP2C19*2, rs12248560 = CYP2C19*17, rs3892097 = CYP2D6*4, rs1902023 = UGT2B15*2, and rs11692021 = UGT1A7*3) and questionnaires. Associations of SNPs with levels of metabolites and their ratios were assessed by multiple linear regression and ordinary logistic regression analyses. Significant associations were observed for CYP2C9*2, CYP2C9*3, CYP2C19*17, and UGT1A7*3 SNPs. The most pronounced was the influence of CYP2C9*2 and *3 on the reduced DEHP biotransformation, with lower levels of metabolites and their ratios in men and women. In contrast, carriers of CYP2C19*17 showed higher urine levels of DEHP metabolites in both genders, and in women also in higher DiNP, DiDP, and DINCH metabolite levels. The presence of UGT1A7*3 was associated with increased metabolite levels of DINCH in men and of DiBP and DBzP in women. Statistical models explained up to 27% of variability in metabolite levels or their ratios. Our observations confirm the effect of CYP2C9*2 and *3 SNPs towards reduced DEHP biotransformation. We show that CYP2C9*2, CYP2C9*3, CYP2C19*17, and UGT1A7*3 SNPs might represent biomarkers of susceptibility or resilience in phthalates and DINCH exposure that have been so far unrecognised.

Prediction of deleterious mutations in coding regions of mammals with transfer learning

The genomes of mammals contain thousands of deleterious mutations. It is important to be able to recognize them with high precision. In conservation biology, the small size of fragmented populations results in accumulation of damaging variants. Preserving animals with less damaged genomes could optimize conservation efforts. In breeding of farm animals, trade-offs between farm performance versus general fitness might be better avoided if deleterious mutations are well classified. In humans, the problem of such a precise classification has been successfully solved, in large part due to large databases of disease-causing mutations. However, this kind of information is very limited for other mammals. Here, we propose to better use information available on human mutations to enable classification of damaging mutations in other mammalian species. Specifically, we apply transfer learning-machine learning methods-improving small dataset for solving a focal problem (recognizing damaging mutations in our companion and farm animals) due to the use of much large datasets available for solving a related problem (recognizing damaging mutations in humans). We validate our tools using mouse and dog annotated datasets and obtain significantly better results in companion to the SIFT classifier. Then, we apply them to predict deleterious mutations in cattle genomewide dataset.

Clinical Significance of UGT1A1 Genetic Analysis in Chinese Neonates with Severe Hyperbilirubinemia

BACKGROUND

Neonatal hyperbilirubinemia is common in Asia, and the importance of genetically determined conditions has been recently recognized. The aim of this study was to assess the clinical utility of genetic testing in Chinese neonates with severe hyperbilirubinemia.

METHODS

Fifty-eight term infants with bilirubin level ≥ 20 mg/dL (342 μmol/L), and 65 controls were enrolled in the study. Variation status of UGT1A1, G6PD, and thalassemia genes in our study cohort was determined by direct sequencing or genotype assays.

RESULTS

Among these case infants, seven were confirmed with G6PD deficiency, four were heterozygous for α- or β-thalassemia, and forty-four were detected with at least one heterozygous UGT1A1 functional variant, including nine homozygous for UGT1A1 variation. As well as the predominant c.211G>A (Gly71Arg) variant, three UGT1A1 coding variants [c.1091C>T (Pro364Leu), c.1352C>T (pro451leu), and c.1456C>T (Tyr486Asp)] were observed in our case neonates. The results of multivariate logistic regressions, adjusted for covariates, revealed odds ratios for neonates who carried heterozygous, homozygous variation at nucleotide 211 of UGT1A1, and G6PD deficiency of 3.47 (1.26-9.55), 12.46 (1.09-142.7) ,and 12.87 (1.32-135.87) compared with those having the wild genotype and normal G6PD activity, respectively.

CONCLUSION

Besides G6PD-deficiency screening, UGT1A1 genetic analysis, and especially the UGT1A1*6(c.211G>A, p.Arg71Gly) polymorphism detection, may be taken into consideration for early diagnosis and treatment of severe hyperbilirubinemic newborns in southern China.

Performance of In Silico Tools for the Evaluation of UGT1A1 Missense Variants

Variations in the gene encoding uridine diphosphate glucuronosyltransferase 1A1 (UGT1A1) are particularly important because they have been associated with hyperbilirubinemia in Gilbert's and Crigler-Najjar syndromes as well as with changes in drug metabolism. Several variants associated with these phenotypes are nonsynonymous single-nucleotide polymorphisms (nsSNPs). Bioinformatics approaches have gained increasing importance in predicting the functional significance of these variants. This study was focused on the predictive ability of bioinformatics approaches to determine the pathogenicity of human UGT1A1 nsSNPs, which were previously characterized at the protein level by in vivo and in vitro studies. Using 16 Web algorithms, we evaluated 48 nsSNPs described in the literature and databases. Eight of these algorithms reached or exceeded 90% sensitivity and six presented a Matthews correlation coefficient above 0.46. The best-performing method was MutPred, followed by Sorting Intolerant from Tolerant (SIFT). The prediction measures varied significantly when predictors such us SIFT, polyphen-2, and Prediction of Pathological Mutations on Proteins were run with their native alignment generated by the tool, or with an input alignment that was strictly built with UGT1A1 orthologs and manually curated. Our results showed that the prediction performance of some methods based on sequence conservation analysis can be negatively affected when nsSNPs are positioned at the hypervariable or constant regions of UGT1A1 ortholog sequences.

Pharmacokinetic considerations for therapies used to treat interstitial cystitis

INTRODUCTION

Interstitial cystitis (IC) or bladder pain syndrome (BPS) is defined as supra-pubic pain related to bladder filling. IC is characterized by a particular symptom complex with no identifiable causes; as with bladder hypersensitivity it is usually associated with urinary frequency and urgency with bladder pain. No current treatments have a significant impact on symptoms over time.

AREAS COVERED

This systematic review examines the pharmacokinetic aspects and adverse event of present IC therapy to highlight appropriate treatment to improve the symptoms of IC. This article reviews material obtained via Medline, PubMed, and EMBASE literature searches up to October 2013.

EXPERT OPINION

The correct approach to IC should consider a multidisciplinary team of specialists and a multimodal treatment package that include psychotherapy, behavior change, physical activation, and analgesic treatment. Unfortunately, a single therapeutic target for IC is not yet known. With regard to pathophysiology and therapy, there is more to discover. The first insult damages the bladder urothelium, hence vehicles that lead the drug to penetrate the wall of the bladder might be a novel strategic approach.

Mutation inactivation of Nijmegen breakage syndrome gene (NBS1) in hepatocellular carcinoma and intrahepatic cholangiocarcinoma

Nijmegen breakage syndrome (NBS) with NBS1 germ-line mutation is a human autosomal recessive disease characterized by genomic instability and enhanced cancer predisposition. The NBS1 gene codes for a protein, Nbs1(p95/Nibrin), involved in the processing/repair of DNA double-strand breaks. Hepatocellular carcinoma (HCC) is a complex and heterogeneous tumor with several genomic alterations. Recent studies have shown that heterozygous NBS1 mice exhibited a higher incidence of HCC than did wild-type mice. The objective of the present study is to assess whether NBS1 mutations play a role in the pathogenesis of human primary liver cancer, including HBV-associated HCC and intrahepatic cholangiocarcinoma (ICC). Eight missense NBS1 mutations were identified in six of 64 (9.4%) HCCs and two of 18 (11.1%) ICCs, whereas only one synonymous mutation was found in 89 control cases of cirrhosis and chronic hepatitis B. Analysis of the functional consequences of the identified NBS1 mutations in Mre11-binding domain showed loss of nuclear localization of Nbs1 partner Mre11, one of the hallmarks for Nbs1 deficiency, in one HCC and two ICCs with NBS1 mutations. Moreover, seven of the eight tumors with NBS1 mutations had at least one genetic alteration in the TP53 pathway, including TP53 mutation, MDM2 amplification, p14ARF homozygous deletion and promoter methylation, implying a synergistic effect of Nbs1 disruption and p53 inactivation. Our findings provide novel insight on the molecular pathogenesis of primary liver cancer characterized by mutation inactivation of NBS1, a DNA repair associated gene.

UGT1A1 genetic analysis as a diagnostic aid for individuals with unconjugated hyperbilirubinemia

OBJECTIVE

To assess the clinical utility of UGT1A1 genetic testing and describe the spectrum and prevalence of UGT1A1 variations identified in pediatric unconjugated hyperbilirubinemia (UCH), and to characterize specific genotype-phenotype relationships in suspected Gilbert and Crigler-Najjar syndromes.

STUDY DESIGN

A retrospective study was conducted to review clinical information and UGT1A1 genotyping data from 181 pediatric patients referred for UCH. In silico analyses were performed to aid in the assessment of novel UGT1A1 variants.

RESULTS

Overall, 146/181 pediatric patients had at least one heterozygous UGT1A1 functional variant. Identified UGT1A1 variants included 17 novel variants, 7 rare star alleles, and 1 rare variant. There were 129 individuals who possessed the TA7 (*28) promoter repeat and 15 individuals who possessed the *6 (c.211G > A) variation. Out of the 104 individuals with accompanying bilirubin levels, 41 individuals did not have identifiable UGT1A1 variants that explained their UCH, although glucose-6-phosphate dehydrogenase deficiency and other causes of UCH could not be ruled out.

CONCLUSION

Much of the observed UCH could be attributed to variation at the UGT1A1 locus, and UGT1A1 testing helped to substantiate a genetic diagnosis, thereby aiding in individual and family disease management. Although UGT1A1 variation plays a large role in UCH, genetic assessment of UGT1A1 alone may not be comprehensive. Assessment of additional genes may also be useful to evaluate genetic causes for UCH.

Performance of computational tools in evaluating the functional impact of laboratory-induced amino acid mutations

Summary: Site-directed mutagenesis is frequently used by scientists to investigate the functional impact of amino acid mutations in the laboratory. Over 10 000 such laboratory-induced mutations have been reported in the UniProt database along with the outcomes of functional assays. Here, we explore the performance of state-of-the-art computational tools (Condel, PolyPhen-2 and SIFT) in correctly annotating the function-altering potential of 10 913 laboratory-induced mutations from 2372 proteins. We find that computational tools are very successful in diagnosing laboratory-induced mutations that elicit significant functional change in the laboratory (up to 92% accuracy). But, these tools consistently fail in correctly annotating laboratory-induced mutations that show no functional impact in the laboratory assays. Therefore, the overall accuracy of computational tools for laboratory-induced mutations is much lower than that observed for the naturally occurring human variants. We tested and rejected the possibilities that the preponderance of changes to alanine and the presence of multiple base-pair mutations in the laboratory were the reasons for the observed discordance between the performance of computational tools for natural and laboratory mutations. Instead, we discover that the laboratory-induced mutations occur predominately at the highly conserved positions in proteins, where the computational tools have the lowest accuracy of correct prediction for variants that do not impact function (neutral). Therefore, the comparisons of experimental-profiling results with those from computational predictions need to be sensitive to the evolutionary conservation of the positions harboring the amino acid change.

Contact:s.kumar@asu.edu

Morphine metabolism, transport and brain disposition

The chemical structures of morphine and its metabolites are closely related to the clinical effects of drugs (analgesia and side-effects) and to their capability to cross the Blood Brain Barrier (BBB). Morphine-6-glucuronide (M6G) and Morphine-3-glucuronide (M3G) are both highly hydrophilic, but only M6G can penetrate the BBB; accordingly, M6G is considered a more attractive analgesic than the parent drug and the M3G. Several hypotheses have been made to explain these differences. In this review we will discuss recent advances in the field, considering brain disposition of M6G, UDP-glucoronosyltransferases (UGT) involved in morphine metabolism, UGT interindividual variability and transport proteins.

Impact of UGT1A1 gene variants on total bilirubin levels in Gilbert syndrome patients and in healthy subjects

The Gilbert syndrome is a benign form of unconjugated hyperbilirubinemia, mainly associated with alterations in UGT1A1 gene. This work investigated the effect of UGT1A1 variants on total bilirubin levels in Gilbert patients (n=45) and healthy controls (n=161). Total bilirubin levels were determined using a colorimetric method; molecular analysis of exons 1-5 and two UGT1A1 promoter regions were performed by direct sequencing and automatic analysis of fragments. Five in silico methods predicted the effect of new identified variants. A significant different allelic distribution, in Gilbert patients and in controls, was found for two promoter polymorphisms. Among patients, 82.2% were homozygous and 17.8% heterozygous for the c.-41_-40dupTA allele; in control group, 9.9% were homozygous and 43.5% heterozygous for this promoter variant, while 46.6% (n=75) presented the [A(TA)6TAA]. For the T>G transition at c.-3279 promoter region, in patients, 86.7% were homozygous and 13.3% heterozygous; in control group, 33.5% were homozygous for the wild type allele, 44.1% were heterozygous and 22.4% homozygous for the mutated allele. The two polymorphisms were in Hardy-Weinberg equilibrium in both groups. Sequencing of UGT1A1 coding region identified nine novel variants, five in patients and four in controls. In silico analysis of these amino acids replacements predicted four of them as benign and three as damaging. In conclusion, we demonstrated that total bilirubin levels are mainly determined by the TA duplication in the TATA-box promoter and by the c.-3279T>G variant. Alterations in the UGT1A1 coding region seem to be associated with increased bilirubin levels, and, therefore, with Gilbert syndrome.

Phenotype prediction of non-synonymous single-nucleotide polymorphisms in human ATP-binding cassette transporter genes

A large number of non-synonymous single-nucleotide polymorphisms (nsSNPs) have been found in human genome, but there is poor knowledge on the relationship between the genotype and phenotype of these nsSNPs. Human ATP-binding cassette (ABC) transporters are able to transport a number of important substrates including endogenous and exogenous compounds. This study aimed to predict the phenotypical impact of nsSNPs of human ABC transporter genes, and the predicted results were further validated by reported phenotypical data from site-directed mutagenesis and clinical genetic studies. One thousand and six hundred thirty-two nsSNPs were found from 49 human ABC transporter genes. Using the PolyPhen and SIFT algorithms, 41.8-53.6% of nsSNPs in ABC transporter genes were predicted to have an impact on protein function. The prediction accuracy was up to 63-85% when compared with known phenotypical data from in vivo and in vitro studies. There was a significant concordance between the prediction results using SIFT and PolyPhen. Of nsSNPs predicted as deleterious, the prediction scores by SIFT and PolyPhen were significantly related to the number of nsSNPs with known phenotypes confirmed by experimental and human studies. The amino acid substitution variants are supposed to be the pathogenetic basis of increased susceptibility to certain diseases with Mendelian or complex inheritance, altered drug resistance and altered drug clearance and response. Predicting the phenotypic consequence of nsSNPs using computational algorithms may provide a better understanding of genetic differences in susceptibility to diseases and drug response. The prediction of nsSNPs in human ABC transporter genes would be useful hints for further genotype-phenotype studies.