Genotyping SNPs in lignin biosynthesis gene (CAD1) and transcription factors (MYB1 and MYB2) exhibits association with wood density in teak (Tectona grandis L.f.)

• Published in: Molecular biology reports Vol. 51; no. 1; p. 169
• Main Authors: Bano, Nuzhat, Mohammad, Naseer, Ansari, Mohammad Israil, Ansari, Shamim Akhtar
• Format: Journal Article
• Language: English
• Published: Dordrecht Springer Netherlands 01.12.2024; Springer Nature B.V
• Subjects: Animal Anatomy; Animal Biochemistry; Biomedical and Life Sciences; Biosynthesis; Decay; economic valuation; exons; forest trees; Gene mapping; Genetic diversity; genetic variation; Genotype; Genotyping; Germplasm; Hardwoods; Histology; India; Lamiaceae; Life Sciences; Lignin; Lignin - genetics; loci; Mathematical models; Morphology; Original Article; outcrossing; phenotypic variation; Phenotypic variations; Polymorphism, Single Nucleotide - genetics; Population genetics; Single-nucleotide polymorphism; Statistical analysis; Tectona grandis; Transcription factors; Transcription Factors - genetics; tropical wood; Wood - genetics; wood density; India; Bonferroni correction; GLM and MLM model; Association mapping; Population structure; Family (genetic) kinship; Linkage disequilibrium
• ISSN: 0301-4851; 1573-4978; 1573-4978
• DOI: 10.1007/s11033-023-09006-y

Background Teak ( Tectona grandis L.f.), an important source of tropical timber with immense economic value, is a highly outcrossing forest tree species. 150 unrelated accessions of teak ( Tectona grandis L.f.) plus trees assembled as clones at National Teak Germplasm Bank, Chandrapur, Maharashtra, India was investigated for association mapping of candidate lignin biosynthesis gene ( CAD1 ) and transcription factors ( MYB1 and MYB2 ). Methods and Results The CAD1 , MYB1 and MYB2 were amplified using specifically designed primers. The amplified sequences were then sequenced and genotyped for 112 SNPs/11 indels. We evaluated the association between SNPs and wood density in teak accessions using GLM and MLM statistical models, with Bonferroni correction applied. The teak accessions recorded an average wood density of 416.69 kg.m −3 (CV 4.97%) and comprised of three loosely structured admixed sub-populations (K = 3), containing 72.05% genetic variation within sub-populations with low intragenic LD (0–21% SNP pairs) at P  < 0.05 and high LD decay (33–934 bp) at R 2  = 0.1. GLM and MLM models discounting systematic biases (Q and K matrices) to avoid false discovery revealed five loci at rare variants (MAF 0.003) and three loci at common variants (MAF 0.05) to be significantly ( P  < 0.05) associated with the wood density. However, the stringent Bonferroni correction (4.06–7.04 × 10 –4 ) yielded only a single associated locus (B1485C/A) from exon of MYB1 transcription factor, contributing to about 10.35% phenotypic variation in wood density trait. Conclusion Scored SNP locus (B1485C/A) can be developed as a molecular probe for selection of improved planting stock with proven wood density trait for a large-scale teak plantation.