Isolation, expression and single nucleotide polymorphisms (SNPs) analysis of LACCASE gene (LkLAC8) from Japanese larch (Larix kaempferi)

Abstract

Nucleotide diversity (π) and linkage disequilibrium (LD) analysis based on SNP marker could provide a sound basis for choosing an association analysis method. Japanese larch (Larix kaempferi) is an important timber coniferous tree species for pulping and papermaking, but its high lignin content has significantly restricted it application potential. In this study, the LACCASE gene, that plays an important regulatory role for lignin biosynthesis, was selected as research target. The full-length cDNA and genomic sequences of the encoding LkLAC8 gene were isolated from the LACCASE expressed sequence tags of the Japanese larch transcriptome database using the rapid amplification of cDNA ends-polymerase chain reaction (RACE-PCR). The cDNA was determined to be 1940 bp, with an open reading frame (ORF, 1734 bp) that encoded a protein of 577 AA. This protein contains four highly specific Cu²⁺ binding sites and 11 glycosylation sites, thus belonging to the LACCASE family. The deduced protein sequence shared an 89% identity with the PtaLAC from Pinus taeda. A real-time PCR analysis showed that the LkLAC8 transcript was expressed predominantly in mature xylem, with moderate levels in the immature xylem, cambium and mature leaves, the lowest in the roots. Lastly, the genomic sequences of LkLAC8 in 40 individuals from six naturally distributed populations of Japanese larch were amplified, and a total of 201 SNPs (103 and 98 mutation types of transition and transversion, respectively) were detected; the frequency of the SNPs was 1/19 bp. Nucleotide diversity among the six populations ranged from 0.0034 to 0.0053, which suggested that there were no significant differences among the populations. The LD analysis showed that the LD level decayed rapidly within the increasing length of the LkLAC8 gene. These results implied that LD mapping and association analysis based on candidate gene may be feasible for the marker-assisted breeding of new germplasms with low lignin in Japanese larch.

Keywords

Gene cloning / LACCASE / Larix kaempferi / Linkage disequilibrium / Real-time PCR / Single nucleotide polymorphisms

References

[1] Beaulieu J, Doerksen T, Boyle B, Clément S, Deslauriers M, Beauseigle S, Blais S, Poulin PL, Lenz P, Caron S, Rigault P, Bicho P, Bousquet J, MacKay J. Association genetics of wood physical traits in the conifer white spruce and relationships with gene expression. Genetics, 2011, 2: 197-212.

[2] Berthet S, Demont-Caulet N, Pollet B, Bidzinsk P, Cézard L, Le Bris P, Borrega N, Hervé J, Blondet E, Balzergue S, Lapierre C, Jouanin L (2011) Disruption of LACCASE4 and 17 results in tissue-specific alterations to lignification of Arabidopsis thaliana stems. Plant Cell 23:1124–1137

[3] Brown GR, Gill GP, Kuntz RJ, Langley CH, Neale DB. Nucleotide diversity and linkage disequilibrium in loblolly pine. Proc Natl Acad Sci USA, 2004, 101: 15255-15260.

[4] Caicedo AL, Williamson S, Hernandez RD, Boyko A, Fledel-Alon A, York TL, Polato NR, Olsen KM, Nielsen R, McCouch SR, Bustamante CD, Purugganan MD. Genome-wide patterns of nucleotide polymorphism in domesticated rice. PLoS Genet, 2007, 3: 1745-1756.

[5] Cheng R, Ma JZ, Elston RC, Li MD. Fine mapping functional sites or regions from case–control data using haplotypes of multiple linked SNPs. Ann Hum Genet, 2005, 69: 102-112.

[6] Chu YG, Su XH. Research progress of single nucleotide polymorphisms in forest trees. Heredites, 2008, 30: 1272-1278.

[7] Du QZ, Pan W, Tian JX, Li BL, Zhang DQ. The UDP-glucuronate decarboxylase gene family in Populus: structure, expression, and association genetics. PLoS ONE, 2013, 8: 1-14.

[8] Dvornyk V, Sirvio A, Mikkonen M. Low nucleotide diversity at the pall locus in the widely distributed Pinus sylvestris. Mol Biol Evol, 2002, 19: 179-188.

[9] Flint-Garcia SA, Thomsberry JM, Buckler ES. Structure of linkage disequilibrium in plants. Annu Rev Plant Bio, 2003, 54: 357-374.

[10] Fu YX, Li WS. Statistical tests of neutrality of mutations. Genetics, 1993, 133: 693-709.

[11] Gehrig HH, Winter K, Cushman J, Borland A, Taybi T. An improved RNA isolation method for succulent plant species rich in polyphenols and polysaccharides. Plant Mol Biol Report, 2000, 18: 369-376.

[12] González-Martínez SC, Wheeler NC, Ersoz E, Nelson DC, Neale DB. Association genetics in Pinus taeda L.I. wood property traits. Genetics, 2007, 175(1): 399-409.

[13] Guerra FP, Wegrzyn JL, Sykes R, Davis MF, Stanton BJ, Neale DB. Association genetics of chemical wood properties in black poplar (Populus nigra). New Phytol, 2013, 197: 162-176.

[14] Guo M, Rupe MA, Zinselneier C, Habben J, Bewen BA, Smith OS (2004) Allelic variation of gene expression in maize hybrids. Plant Cell 16:1707–1716

[15] Guo Q, Wang BL, Wang BW, Li BL, Zhang DQ. Isolation, expression and single nucleotide polymorphisms analysis of PtDREB2A in Populus tomentosa. Sci silvae sin, 2011, 47: 49-56.

[16] Gupta K. Linkage disequilibrium and association studies in higher plants: present status and future prospects. Plant Mol Biol, 2005, 57: 46-485.

[17] Harald C (2004) Laccases: structure, reactions, distribution. Micron 35(1–2):93–96

[18] Heuertz M, De Paoli E, Källman T, Larsson H, Jurman I, Morgante M, Lascoux M, Gyllenstrand N. Multi-locus patterns of nucleotide diversity, linkage disequilibrium and demographic history of Norway spruce (Picea abies(L.) Karst). Genetics, 2006, 174: 2095-2105.

[19] Hudson RR, Kreitman M, Aguade M. A test of neutral molecular evolution based on nucleotide data. Genetics, 1987, 116: 153-159.

[20] Ingvarsson PK. Nucleotide polymorphism and linkage disequilibrium within and among natural populations of European aspen (Populus tremula L. Salicaceae). Genetics, 2005, 169: 945-953.

[21] Jin L, Bao JS. Progress on the trait-marker association analysis in plants. Mol Plant Breed, 2009, 7: 1048-1063.

[22] Kado T, Yoshimaru H, Tsumura Y. DNA variation in a conifer Cryptomeria japonica (Copressaceae sensu lato). Genetics, 2014, 164: 1547-1599.

[23] Krutovsky KV, Neale DB. Nucleotide diversity and linkage disequilibrium in cold-hardiness and wood quality-related candidate genes in douglas fir. Genetics, 2005, 171: 2029-2041.

[24] Kuittinen H, Aguage M. Nucleotide variation at the chalcone isomerase locus in Arabidopsis thaliana. Genetics, 2000, 155: 863-872.

[25] LaFayette PR, Eriksson KE, Dean JF. Characterization and heterologous expression of laccase cDNAs from xylem tisssues yellow poplar (Liriodendron tulipifera). Plant Mol Biol, 1999, 40: 23-35.

[26] Lu SF, Li QZ, Wei HR, Chang MJ. Ptr-miR397a is a negative regulator of laccase genes affecting lignin content in Populus trichocarpa. Plant Biol, 2013, 110: 10848-10853.

[27] Ma C, Sun X. Larch genetic improvement and its future development in China. World For Res, 2008, 21: 58-63.

[28] McCaig BC, Meagher RB, Dean JF. Gene structure and molecular analysis of the Laccase-like multicopper oxidase (LMCO) gene family in Arabidopsis thaliana. Planta, 2005, 221: 619-636.

[29] Mohamad SB, Ong AL, Ripen AM. Evolutionary trace analysis at the ligand binding site of laccase. Bioinformation, 2008, 2: 369-372.

[30] Neale DB, Savolainen O. Association genetics of complex traits in conifers. Trends Plant Sci, 2004, 9: 325-330.

[31] Nei M (1987) Molecular evolutionary genetics. Columbia University Press, New York

[32] Nordborg M, Hu TT, Ishino Y, Jhaveri J, Toomajian C, Zheng H, Bakker E. The pattern of polymorphism in Arabidopsis thaliana. PLoS Biol, 2005, 3: 1289-1299.

[33] Nowotny P, Kwon JM, Goate AM. SNP analysis to dissect human traits. Curr Opin Neurobiol, 2001, 11: 637-641.

[34] Olivier M (2004) From SNPs to function: the effect of sequence variation on gene expression. Physiol Genomics 16:182–183

[35] Paques LE. Roles of European and Japanese larch in the genetic control of growth, architecture and wood quality traits in interspecific hybrids (Larix x eurolepis Henry). Ann Forest Sci, 2004, 61: 25-33.

[36] Pot D, McMillan L, Echt C, Le Provost G, Garnier-Géré P, Cato S, Plomion C. Nucleotide variation in genes involved in wood formation in two pine species. New Phytol, 2005, 167: 101-112.

[37] Ranocha P, Dougall MC, Hawkins G, Sterjiades S. Biochemical characterization, molecular cloning and expression of laccases-a divergent gene family in poplar. Eur J Biochem, 1999, 159: 1-2.

[38] Reinhammar B. Lontie R. Laccase. Copper proteins and copper enzymes, 1984, Boca Raton: CRC Press Inc 1 35

[39] Richardson A, Duncan J, Mcdougall GJ. Oxidase activity in lignifying xylem of a taxonomically diverse range of trees: identification of a conifer Laccase. Tree Physiol, 2000, 20: 1039-1047.

[40] Rozas J, Sánchez-DelBarrio JC, Messeguer X, Rozas R (2003) DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics 19:2496–2497

[41] Salisbury BA, Pungliya M, Choi JY, Jiang RH, Sun XJ, Stephens JC. SNP and haplotype variation in the human genome. Mutat Res, 2003, 526: 53-61.

[42] Sato Y, Wuli B, Sederoff R. Molecular cloning and expression of eight laccase cDNAs in loblolly pine (Pinus taeda). J Plant Res, 2001, 114: 147-155.

[43] Schuetz M, Smith R, Ellis B. Xylem tissue specification, patterning, and differentiation mechanisms. Exp Botany, 2013, 64: 11-31.

[44] Sun XM, Zhang SG, Li SY, Hou YM. Multi traits selection of open pollinated Larix kaempferi families for pulpwood purpose. Sci silvae sin, 2005, 41: 48-54.

[45] Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics, 1989, 123: 585-595.

[46] Takata K, Kurinobu S, Akio K, Yasue K, Tamai Y, Kisanuki M. Bibliography on Larix kaempferi (Larix kaempferi (Lamb.) Carr.). Eurasian J For Res, 2005, 8: 111-126.

[47] Tenaillon M, Sawkins MC, Long AD, Gaut RL, Doebley JF, Gaut BS. Patterns of DNA sequence polymorphism along chromosome1 of maize (Zea mays ssp. mays L.). Proc Natl Acad Sci USA, 2001, 98: 9161-9166.

[48] Tian JX, Du QZ. Allelic variation in PtGA20Ox associates with growth and wood properties in Populus spp. PLoS ONE, 2012, 12: 1-11.

[49] Watterson GA (1975) On the number of segregating sites in genetical models without recombination. Theor Popul Biol 7:256–276

[50] Witayakran S (2008) Laccase in organic synthesis and its applications. Doctor Degree, Georgia Institute of Technology United States

[51] Yasushi S, Ross W, Whetten (2006) Characterization of two laccases of loblolly pine (Pinus taeda) expressed in tobacco BY-2 cells. J Plant Res 119:581–588

[52] Yi M (2014) Polymorphic single nucleotide polymorphisms (SNP) loci within wood quality-related candidate genes are associated with growth and wood properties in Larix kaempferi. Chinese Academy of Forestry, Beijing

[53] Yi M, Zhang SG, Xie YH, Sun XM. Isolation and single nucleotide polymorphisms analysis of caffeic acid O-methyltransferase gene (LkCOMT) from Larix kaempferi. For Res, 2013, 26: 52-59.

[54] Yoshida KT, Naito S, Takada G. cDNA cloning of regeneration-specific genes in rice by differential screening of randomly amplified cDNAs using RAPD primers. Plant Cell Physiol, 1994, 35: 1003-1009.

[55] Zhang DQ, Zhang ZY. Single nucleotide polymorphisms (SNPs) discovery and linkage disequilibrium (LD) in forest trees. For Stud China, 2005, 7: 1-14.

[56] Zhang DQ, Xu BH, Yang XH, Zhang ZY, Li BL. The sucrose synthase gene family in Populus: structure, expression and evolution. Tree Genet Genomes, 2011, 7: 443-456.

[57] Zhao Q, Nakashima Jin, Fang Chen F. Laccase is necessary and non-redundant with peroxidase for lignin polymerization during vascular development in Arabidopsis. Plant Cell, 2013, 25: 3976-3987.

[58] Zollner S, Wen X, Pritchard JK. Association mapping and fine mapping with tree LD. Bioinformatics, 2005, 21: 3168-3170.