Complexity of gene paralogues resolved in biosynthetic pathway of hepatoprotective iridoid glycosides in a medicinal herb, Picrorhiza kurroa through differential NGS transcriptomes

Picrorhiza kurroa is a medicinal herb with diverse pharmacological applications due to the presence of iridoid glycosides, picroside-I (P-I), and picroside-II (P-II), among others. Any genetic improvement in this medicinal herb can only be undertaken if the biosynthetic pathway genes are correctly identified. Our previous studies have deciphered biosynthetic pathways for P-I and P-II, however, the occurrence of multiple copies of genes has been a stumbling block in their usage. Therefore, a methodological strategy was designed to identify and prioritize paralogues of pathway genes associated with contents of P-I and P-II. We used differential transcriptomes varying for P-I and P-II contents in different tissues of P. kurroa. All transcripts for a particular pathway gene were identified, clustered based on multiple sequence alignment to notify as a representative of the same gene (>= 99% sequence identity) or a paralogue of the same gene. Further, individual paralogues were tested for their expression level via qRT-PCR in tissue-specific manner. In total 44 paralogues in 14 key genes have been identified out of which 19 gene paralogues showed the highest expression pattern via qRT-PCR. Overall analysis shortlisted 6 gene paralogues, PKHMGR3, PKPAL2, PKDXPS1, PK4CL2, PKG10H2 and PKIS2 that might be playing role in the biosynthesis of P-I and P-II, however, their functional analysis need to be further validated either through gene silencing or over-expression. The usefulness of this approach can be expanded to other non-model plant species for which transcriptome resources have been generated.