Whole genome based identification of BAHD acyltransferase gene involved in piperine biosynthetic pathway in black pepper

Black pepper (Piper nigrum L.), a crop of the genus Piper, is an important spice that has both economic and ecological significance. It is widely regarded as the "King of Spices" because of its pungency, attributed to the presence of piperine. BAHD acyl transferase, the crucial enzyme involved in the final step in piperine biosynthesis was the focus of our study and the aim was to identify the candidate isoform involved in biosynthesis of piperine. Reference genome-based analysis of black pepper identified six BAHD-AT isoforms and mapping of these sequences revealed that the isoforms were situated on six distinct chromosomes. By using specific primers for each of these transcripts, qPCR analysis was done in different tissues as well as berry stages to obtain detectable amplification products. Expression profiles of isoforms from chromosome 6 correlated well with piperine content compared to other five isoforms, across tissues and was therefore assumed to be involved in biosynthesis of piperine. In addition to this, we could also identify the binding sites of MYB transcription factor in the cis-regulatory regions of the isoforms. We also used in-silico docking and molecular dynamics simulation to calculate the binding free energy of the ligand and confirmed that among all the isoforms, BAHD-AT from chromosome 6 had the lowest free binding energy and highest affinity towards the ligand. Our findings are expected to aid the identification of new genes connecting enzymes involved in the biosynthetic pathway of piperine, which will have major implications for future research in metabolic engineering.