Chromosome-level genome assembly of Aristolochia contorta provides insights into the biosynthesis of benzylisoquinoline alkaloids and aristolochic acids

Abstract

Aristolochic acids (AAs) and their derivatives are present in multiple Aristolochiaceae species that have been or are being used as medicinal materials. During the past decades, AAs have received increasing attention because of their nephrotoxicity and carcinogenicity. Elimination of AAs from medicinal materials using biotechnological approaches is important for improving medication safety. However, it has not been achieved because of the limited information available on AA biosynthesis. Here, we report a high-quality, reference-grade genome assembly of the AA-containing vine Aristolochia contorta. The total size of the assembly is 209.27 Mb, and it is assembled into 7 pseudochromosomes. Synteny analysis, Ks distribution, and 4DTv suggest an absence of whole-genome duplication (WGD) events in Aristolochia contorta after the angiosperm-wide WGD. Based on genomic, transcriptomic, and metabolic data, pathways and candidate genes were proposed for benzylisoquinoline alkaloid (BIA) and AA biosynthesis in A. contorta. Five O-methyltransferase genes, including AcOMT1–3, AcOMT5, and AcOMT7, were cloned and functionally characterized. The results provide a high-quality reference genome for AA-containing species of Aristolochiaceae. They lay a solid foundation for further elucidation of AA biosynthesis and regulation and for the molecular breeding of Aristolochiaceae medicinal materials.

References

[1] Stiborová M, Arlt VM, Schmeiser HH . DNA adducts formed by aristolochic acid are unique biomarkers of exposure and explain the initiation phase of upper urothelial cancer. Int J Mol Sci. 2017; 18: 2144.

[2] Zhou Q, Pei J, Poon J et al. Worldwide research trends on aristolochic acids (1957-2017): suggestions for researchers. PLoS One. 2019; 14: 1-23.

[3] Balachandran P, Wei F, Lin R et al. Structure activity relationships of aristolochic acid analogues: toxicity in cultured renal epithelial cells. Kidney Int. 2005; 67: 1797-805.

[4] Grosse Y, Baan R, Straif K et al. A review of human carcinogens part a: pharmaceuticals. Lancet Oncol. 2009; 10: 13-4.

[5] Jong T, Lee MR, Hsiao SS et al. Analysis of aristolochic acid in nine sources of Xixin, a traditional Chinese medicine, by liquid chromatography/atmospheric pressure chemical ionization/tandem mass spectrometry. J Pharm Biomed Anal. 2003; 33: 831-7.

[6] Zhao H, Jiang H . Determination of aristolochic acid a in Saruma henryi by HPLC. Guihaia. 2009; 29: 548-51.

[7] Nishida R, Weintraub JD, Feeny P et al. Aristolochic acids from Thottea spp. (Aristolochiaceae) and the osmeterial secretions of Thottea-feeding troidine swallowtail larvae (Papilionidae). J Chem Ecol. 1993; 19: 1587-94.

[8] Luan Y . Thinking and discussion on the carcinogenicity of aristolochic acids, a class of compounds from traditional Chinese medicine: a summary of the presentations in the international symposium on the safety of traditional Chinese medicine. Mod Tradit Chinese Med Mater Medica-World Sci Technol. 2019; 1-7.

[9] Chinese Pharmacopoeia Comission . Pharmacopoeia of the People Republic of China 2020 Volume 1. (China Medical Science Press, 2020).

[10] Spenser DI, Tiwari HP . Biosynthesis of aristolochic acid Chem Commun. 1966; 55-6.

[11] Comer F, Tiwari HP, Spenser ID . Biosynthesis of aristolochic acid. Can J Chem. 1969; 47: 481-7.

[12] Schutte HR, Orban U, Mothes K . Biosynthesis of aristolochic acid. Eur J Biochem. 1967; 1: 70-2.

[13] Sharma V, Bhakuni DS, Kapil RS . Biosynthesis of aristolochic acid JCS Perkin I 40. 1982; 1153-5.

[14] Wang X, Hui F, Yang Y et al. Deep sequencing and transcriptome analysis to identify genes related to biosynthesis of aristolochic acid in Asarum heterotropoides. Sci Rep. 2018; 8: 1-14.

[15] Liu X, Liu Y, Huang P et al. The genome of medicinal plant Macleaya cordata provides new insights into benzylisoquinoline alkaloids metabolism. Mol Plant. 2017; 10: 975-89.

[16] Singh A, Menéndez-Perdomo IM, Facchini PJ . Benzylisoquinoline alkaloid biosynthesis in opium poppy: an update. Phytochem Rev. 2019; 18: 1457-82.

[17] Dembitsky VM, Gloriozova TA, Poroikov VV . Naturally occurring plant isoquinoline N-oxide alkaloids: their pharmacological and SAR activities. Phytomedicine. 2015; 22: 183-202.

[18] Comission CP . Pharmacopoeia of the People Republic of China Vol. 1. China Medical Science Press; 2015.

[19] Chen M, Zhu Z . Progress on study of pharmacological functions of Aristolochia L. Plant. 2007; 3: 59-62.

[20] Li K . The pharmacology of Aristolochia species. Acta Acad Med CPAPF. 2000; 9: 230-2.

[21] Zhao H, Liu X . Review on the studies of Chinese Aristolochia L. herbs. J Henan Univ Natural Sci. 2003; 33: 73-7.

[22] Tian H, Cheng X, Wang C et al. Simultaneous determination of four constituents in Aristolochia debilis by HPLC. Chinese Tradit Pat Med. 2016; 38: 560-5.

[23] Liu J, Dai Z, Cheng X et al. Research progress on detection and analysis of aristolochic acids. Drug Eval Res. 2019; 1644-9.

[24] Lou F, Linsheng D, Li L . Chemical studies on Atristolochia contorta Bunge. Chinese Tradit Herb Drugs. 1986; 17: 390-1.

[25] Zhang C, Yu J, Wang B et al. Analysis of volatile oil constituents in the fruits of Aristolochia contorta by GC-MS. Chin J Nat Med. 2004; 2: 126-8.

[26] Sun R, Zhang L, Luo Z et al. Overview of utilization value and cultivation technique of the medicinal plants Aristolochia contorta Bunge. J Anhui Agri Sci. 2011; 39: 14620-1.

[27] Chen Y, Yu L . Studies on the chemical constituents of Aristolochia Contorta. J Yunnan Norm Univ. 2005; 25: 41-4.

[28] Wickett NJ, Mirarab S, Nguyen N et al. Phylotranscriptomic analysis of the origin and early diversification of land plants. Proc Natl Acad Sci. 2014; 111: E4859-68.

[29] Zeng L, Zhang Q, Sun R et al. Resolution of deep angiosperm phylogeny using conserved nuclear genes and estimates of early divergence times. Nat Commun. 2014; 5: 4956.

[30] Chaw S, Liu YC, Wu YW et al. Stout camphor tree genome fills gaps in understanding of flowering plant genome evolution. Nat Plants. 2019; 5: 63-73.

[31] One Thousand Plant Transcriptomes Initiative . One thousand plant transcriptomes and the phylogenomics of green plants. Nature. 2019; 574: 679-85.

[32] Dong S, Liu M, Liu Y et al. The genome of Magnolia biondii Pamp. Provides insights into the evolution of Magnoliales and biosynthesis of terpenoids. Hortic Res. 2021; 8.

[33] Chen J, Hao Z, Guang X et al. Liriodendron genome sheds light on angiosperm phylogeny and species-pair differentiation. Nat Plants. 2019; 5: 18-25.

[34] Chen SP, Sun WH, Xiong YF et al. The Phoebe genome sheds light on the evolution of magnoliids. Hortic Res. 2020; 7: 146.

[35] Qin L, Hu Y, Wang J et al. Insights into angiosperm evolution, floral development and chemical biosynthesis from the Aristolochia fimbriata genome. Nat Plants. 2021; 7: 1239-53.

[36] Chen YC, Li Z, Zhao YX et al. The Litsea genome and the evolution of the laurel family. Nat Commun. 2020; 11: 1675.

[37] Susko E, Roger AJ . Long Branch attraction biases in Phylogenetics. Syst Biol. 2021; 70: 838-43.

[38] Jost M, Samain MS, Marques I et al. Discordant Phylogenomic placement of Hydnoraceae and Lactoridaceae within Piperales using data from all three genomes. Front Plant Sci. 2021; 12.

[39] Ren R, Wang H, Guo C et al. Widespread whole genome duplications contribute to genome complexity and species diversity in angiosperms. Mol Plant. 2018; 11: 414-28.

[40] Guo L, Winzer T, Yang X et al. The opium poppy genome and morphinan production. Science. 2018; 362: 343-7.

[41] Hu L, Xu Z, Wang M et al. The chromosome-scale reference genome of black pepper provides insight into piperine biosynthesis. Nat Commun. 2019; 10: 4702.

[42] Jaillon O, Aury J, Noel B et al. The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007; 449: 463-7.

[43] Salse J . Ancestors of modern plant crops. Curr Opin Plant Biol. 2016; 30: 134-42.

[44] Zhang L, Chen F, Zhang X et al. The water lily genome and the early evolution of flowering plants. Nature. 2020; 577: 79-84.

[45] Shi T, Rahmani RS, Gugger PF et al. Distinct expression and methylation patterns for genes with different fates following a single whole-genome duplication in flowering plants. Mol Biol Evol. 2020; 37: 2394-413.

[46] Badouin H, Gouzy J, Grassa CJ et al. The sunflower genome provides insights into oil metabolism, flowering and Asterid evolution. Nature. 2017; 546: 148-52.

[47] Amborella Genome Project . The Amborella genome and the evolution of flowering plants. Science. 2013; 342: 1241089-9.

[48] Ghelardini C, Galeotti N, Di Cesare Mannelli L et al. A. Local anaesthetic activity of β-caryophyllene. ChemMedChem. 2001; 56: 387-9.

[49] Liu X, Luo Y, Wu H et al. Systematic analysis of O-methyltransferase gene family and identification of potential members involved in the formation of O-methylated flavonoids in citrus. Gene. 2016; 575: 458-72.

[50] Morris JS, Facchini PJ . Molecular origins of functional diversity in benzylisoquinoline alkaloid methyltransferases. Front Plant Sci. 2019; 10.

[51] Sato F, Tsujita T, Katagiri Y et al. Purification and characterization of S-adenosyl-L-methionine:norcoclaurine 6-O-methyltransferase from cultured Coptis japonica cells. Eur J Biochem. 1994; 225: 125-31.

[52] Robin AY, Giustini C, Graindorge M et al. Crystal structure of norcoclaurine-6-O-methyltransferase, a key rate-limiting step in the synthesis of benzylisoquinoline alkaloids. Plant J. 2016; 87: 641-53.

[53] Ounaroon A, Decker G, Schmidt J et al. (R,S)-Reticuline 7-O-methyltransferase and (R,S)-norcoclaurine 6-O-methyltransferase of Papaver somniferum - cDNA cloning and characterization of methyl transfer enzymes of alkaloid biosynthesis in opium poppy. Plant J. 2003; 36: 808-19.

[54] Ziegler J, Diaz-Chávez ML, Kramell R et al. Comparative macroarray analysis of morphine containing Papaver somniferum and eight morphine free Papaver species identifies an O-methyltransferase involved in benzylisoquinoline biosynthesis. Planta. 2005; 222: 458-71.

[55] Morishige T, Tsujita T, Yamada Y et al. Molecular characterization of the S-Adenosyl-l-methionine:3′-Hydroxy-N-methylcoclaurine 4′-O-Methyltransferase involved in Isoquinoline alkaloid biosynthesis in Coptis japonica. J Biol Chem. 2000; 275: 23398-405.

[56] Mizutani M, Ohta D . Diversification of P450 genes during land plant evolution. Annu Rev Plant Biol. 2010; 61: 291-315.

[57] Weitzel C, Simonsen HT . Cytochrome P450-enzymes involved in the biosynthesis of mono- and sesquiterpenes. Phytochem Rev. 2015; 14: 7-24.

[58] Gou M, Ran X, Martin DW et al. The scaffold proteins of lignin biosynthetic cytochrome P450 enzymes. Nat Plants. 2018; 4: 299-310.

[59] Schoenbohm C, Martens S, Eder C et al. Identification of the Arabidopsis thaliana flavonoid 3′-hydroxylase gene and functional expression of the encoded P450 enzyme. Biol Chem. 2000; 381: 749-53.

[60] Mizutani M, Sato F . Unusual P450 reactions in plant secondary metabolism. Arch Biochem Biophys. 2011; 507: 194-203.

[61] Ikezawa N, Iwasa K, Sato F . Molecular cloning and characterization of CYP80G2, a cytochrome P450 that catalyzes an intramolecular C-C phenol coupling of (S)-reticuline in magnoflorine biosynthesis, from cultured Coptis japonica cells. J Biol Chem. 2008; 283: 8810-21.

[62] Facchini PJ, Hagel JM, Liscombe DK et al. Opium poppy: blueprint for an alkaloid factory. Phytochem Rev. 2007; 6: 97-124.

[63] Hagel JM, Facchini PJ . Benzylisoquinoline alkaloid metabolism: a century of discovery and a brave new world. Plant Cell Physiol. 2013; 54: 647-72.

[64] Beaudoin GAW, Facchini PJ . Benzylisoquinoline alkaloid biosynthesis in opium poppy. Planta. 2014; 240: 19-32.

[65] He S, Liang YL, Cong K et al. Identification and characterization of genes involved in benzylisoquinoline alkaloid biosynthesis in Coptis species. Front Plant Sci. 2018; 9: 731.

[66] Xu T, Kuang T, Du H et al. Magnoflorine: a review of its pharmacology, pharmacokinetics and toxicity. Pharmacol Res. 2020; 152: 104632-10.

[67] Choi K-B, Morishige T, Shitan N et al. Molecular cloning and characterization of coclaurine N-methyltransferase from cultured cells of Coptis japonica. J Biol Chem. 2002; 277: 830-5.

[68] Canedo-Téxon A, Ramon-Farias F, Monribot-Villanueva JL et al. Novel findings to the biosynthetic pathway of magnoflorine and taspine through transcriptomic and metabolomic analysis of Croton draco (Euphorbiaceae). BMC Plant Biol. 2019; 19: 1-21.

[69] Zhong F, Huang L, Qi L et al. Full-length transcriptome analysis of Coptis deltoidea and identification of putative genes involved in benzylisoquinoline alkaloids biosynthesis based on combined sequencing platforms. Plant Mol Biol. 2020; 102: 477-99.

[70] Minami H, Dubouzet E, Iwasa K et al. Functional analysis of norcoclaurine synthase in Coptis japonica. J Biol Chem. 2007; 282: 6274-82.

[71] Li J, Lee E, Chang L et al. Genes encoding norcoclaurine synthase occur as tandem fusions in the Papaveraceae. Sci Rep. 2016; 6: 39256.

[72] Alcantara J, Bird DA, Franceschi VR et al. Sanguinarine biosynthesis is associated with the endoplasmic reticulum in cultured opium poppy cells after elicitor treatment. Plant Physiol. 2005; 138: 173-83.

[73] Pauli HH, Kutchan TM . Molecular cloning and functional heterologous expression of two alleles encoding (S)-N-methylcoclaurine 3′-hydroxylase (CYP80B1), a new methyl jasmonate-inducible cytochrome P450-dependent monooxygenase of benzylisoquinoline alkaloid biosynthesis. Plant J. 1998; 13: 793-801.

[74] Frick S, Kramell R, Kutchan TM . Metabolic engineering with a morphine biosynthetic P450 in opium poppy surpasses breeding. Metab Eng. 2007; 9: 169-76.

[75] Ikezawa N, Tanaka M, Nagayoshi M et al. Molecular cloning and characterization of CYP719, a methylenedioxy bridge-forming enzyme that belongs to a novel P450 family, from cultured Coptis japonica cells. J Biol Chem. 2003; 278: 38557-65.

[76] Gurkok T, Ozhuner E, Parmaksiz I et al. Functional characterization of 4′OMT and 7OMT genes in BIA biosynthesis. Front Plant Sci. 2016; 7.

[77] Xu D, Lin H, Tang Y et al. Integration of full-length transcriptomics and targeted metabolomics to identify benzylisoquinoline alkaloid biosynthetic genes in corydalis yanhusuo. Hortic Res. 2021; 8: 16.

[78] Meelaph T, Kobtrakul K, Chansilpa NN et al. Coregulation of biosynthetic genes and transcription factors for aporphine-type alkaloid production in wounded lotus provides insight into the biosynthetic pathway of nuciferine. ACS Omega. 2018; 3: 8794-802.

[79] Mao WW, Gao W, Liang ZT et al. Characterization and quantitation of aristolochic acid analogs in different parts of Aristolochiae fructus, using UHPLC-Q/TOF-MS and UHPLC-QqQ-MS. Chin J Nat Med. 2017; 15: 392-400.

[80] Debelle FD, Vanherweghem JL, Nortier JL . Aristolochic acid nephropathy: a worldwide problem. Kidney Int. 2008; 74: 158-69.

[81] Koren S, Walenz BP, Berlin K et al. Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res. 2017; 27: 722-36.

[82] Vaser R, Sović I, Nagarajan N, Šikić M . Fast and accurate de novo genome assembly from long uncorrected reads. Genome Res. 2017; 27: 737-46.

[83] Walker BJ, Abeel T, Shea T et al. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS One. 2014; 9: e112963.

[84] Parra G, Bradnam K, Korf I . CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes. Bioinformatics. 2007; 23: 1061-7.

[85] Li H, Durbin R . Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009; 25: 1754-60.

[86] Burton JN, Adey A, Patwadhan RP et al. Chromosome-scale scaffolding of de novo genome assemblies based on chromatin interactions. Nat Biotechnol. 2013; 31: 1119-25.

[87] Servant N, Varoquaux N, Lajoie BR et al. HiC-pro: an optimized and flexible pipeline for Hi-C data processing. Genome Biol. 2015; 16: 259.

[88] Xu Z, Wang H . LTR_FINDER: an efficient tool for the prediction of full-length LTR retrotransposons. Nucleic Acids Res. 2007; 35: W265-8.

[89] Price AL, Jones NC, Pevzner PA . De novo identification of repeat families in large genomes. Bioinformatics. 2005; 21: i351-8.

[90] Hoede C, Arnoux S, Moisset M et al. PASTEC: an automatic transposable element classification tool. PLoS One. 2014; 9: e91929.

[91] Jurka J, Kapitonov VV, Pavlicek A et al. Repbase update, a database of eukaryotic repetitive elements. Cytogenet Genome Res. 2005; 110: 462-7.

[92] Tarailo-Graovac M, Chen N . Using repeatMasker to identify repetitive elements in genomic sequencesV. Curr Protoc Bioinformatics. 2009; 25.

[93] Griffiths-Jones S . Rfam: annotating non-coding RNAs in complete genomes. Nucleic Acids Res. 2004; 33: D121-4.

[94] Lowe TM, Eddy SR . tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997; 25: 955-64.

[95] Birney E. GeneWise and genomewise. Genome Res. 2004; 14: 988-95.

[96] Burge C, Karlin S . Prediction of complete gene structures in human genomic DNA. J Mol Biol. 1997; 268: 78-94.

[97] Stanke M, Waack S . Gene prediction with a hidden Markov model and a new intron submodel. Bioinformatics. 2003; 19: ii215-25.

[98] Majoros WH, Pertea M, Salzberg SL . TigrScan and GlimmerHMM: two open source ab initio eukaryotic gene-finders. Bioinformatics. 2004; 20: 2878-9.

[99] Blanco, E., Parra, G. & Guigó, R. Using geneid to identify genes. in Current Protocols in Bioinformatics (John Wiley & Sons, Inc., 2007).

[100] Keilwagen J, Hartung F, Paulini M et al. Combining RNA-seq data and homology-based gene prediction for plants, animals and fungi. BMC Bioinformatics. 2018; 19: 189.

[101] Kim D, Langmead B, Salzberg SL . HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015; 12: 357-60.

[102] Pertea M, Pertea GM, Antonescu CM et al. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol. 2015; 33: 290-5.

[103] Tang S, Lomsadze A, Borodovsky M . Identification of protein coding regions in RNA transcripts. Nucleic Acids Res. 2015; 43: e78-8.

[104] Campbell MA, Haas BJ, Hamilton JP et al. Comprehensive analysis of alternative splicing in rice and comparative analyses with Arabidopsis. BMC Genomics. 2006; 7: 327.

[105] Altschul SF, Gish W, Miller W et al. Basic local alignment search tool. J Mol Biol. 1990; 215: 403-10.

[106] Marchler-Bauer A, Lu S, Anderson JB et al. CDD: a conserved domain database for the functional annotation of proteins. Nucleic Acids Res. 2011; 39: D225-9.

[107] Koonin EV, Fedorova ND, Jackson JD et al. A comprehensive evolutionary classification of proteins encoded in complete eukaryotic genomes. Genome Biol. 2004; 5: R7.

[108] Kanehisa M . KEGG: Kyoto Encyclopedia of genes and genomes. Nucleic Acids Res. 2000; 28: 27-30.

[109] Boeckmann B. The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003. Nucleic Acids Res. 2003; 31: 365-70.

[110] Nguyen LT, Schmidt HA, Haeseler A et al. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2015; 32: 268-74.

[111] Katoh K, Asimenos G, Toh H . Multiple alignment of DNA sequences with MAFFT. 2009; 39-64.

[112] Suyama M, Torrents D, Bork P . PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments. Nucleic Acids Res. 2006; 34: W609-12.

[113] Talavera G, Castresana J . Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst Biol. 2007; 56: 564-77.

[114] Kalyaanamoorthy S, Minh BQ, Wong TKF et al. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017; 14: 587-9.

[115] Yang Z . PAML: a program package for phylogenetic analysis by maximum likelihood. Bioinformatics. 1997; 13: 555-6.

[116] Emms DM, Kelly S . OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biol. 2019; 20: 238.

[117] Mi H, Muruganujan A, Huang X et al. Protocol update for large-scale genome and gene function analysis with the PANTHER classification system (v.14.0). Nat Protoc. 2019; 14: 703-21.

[118] Yu G, Wang L, Han Y et al. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012; 16: 284-7.

[119] Han MV, Thomas GWC, Lugo-Martinez J et al. Estimating gene gain and loss rates in the presence of error in genome assembly and annotation using CAFE 3. Mol Biol Evol. 2013; 30: 1987-97.

[120] Buchfink B, Xie C, Huson DH . Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015; 12: 59-60.

[121] Wang Y, Tang H, Debarry JD et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012; 40: e49-9.

[122] Zwaenepoel A, Van de Peer Y . Wgd-simple command line tools for the analysis of ancient whole-genome duplications. Bioinformatics. 2019; 35: 2153-5.

[123] Nelson DR, Koymans L, Kamataki T et al. P450 superfamily: update on new sequences, gene mapping, accession numbers and nomenclature. Pharmacogenetics. 1996; 6: 1-42.

[124] Florea L, Song L, Salzberg SL . Thousands of exon skipping events differentiate among splicing patterns in sixteen human tissues. F1000Research. 2013; 2: 188.

[125] Smith CA, Want EJ, O’Maille G et al. XCMS: processing mass spectrometry data for metabolite profiling using nonlinear peak alignment, matching, and identification. Anal Chem. 2006; 78: 779-87.

[126] Chen C, Chen H, Zhang Y et al. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020; 13: 1194-202.1