The tissue-specific chromatin accessibility landscape of Papaver somniferum

Genomic and cell-specific regulation of benzylisoquinoline alkaloid biosynthesis in opium poppy

Opium poppy is a crop of great commercial value as a source of several opium alkaloids for the pharmaceutical industries including morphine, codeine, thebaine, noscapine, and papaverine. Most enzymes involved in benzylisoquinoline alkaloid (BIA) biosynthesis in opium poppy have been functionally characterized, and opium poppy currently serves as a model system to study BIA metabolism in plants. BIA biosynthesis in opium poppy involves two biosynthetic gene clusters associated respectively with the morphine and noscapine branches. Recent reports have shown that genes in the same cluster are co-expressed, suggesting they might also be co-regulated. However, the transcriptional regulation of opium poppy BIA biosynthesis is not well studied. Opium poppy BIA biosynthesis involves three cell types associated with the phloem system: companion cells, sieve elements, and laticifers. The transcripts and enzymes associated with BIA biosynthesis are distributed across cell types, requiring the translocation of key enzymes and pathway intermediates between cell types. Together, these suggest that the regulation of BIA biosynthesis in opium poppy is multilayered and complex, involving biochemical, genomic, and physiological mechanisms. In this review, we highlight recent advances in genome sequencing and single cell and spatial transcriptomics with a focus on how these efforts can improve our understanding of the genomic and cell-specific regulation of BIA biosynthesis. Such knowledge is vital for opium poppy genetic improvement and metabolic engineering efforts targeting the modulation of alkaloid yield and composition.

Two ubiquitous aldo-keto reductases in the genus Papaver support a patchwork model for morphine pathway evolution

The evolution of morphinan alkaloid biosynthesis in plants of the genus Papaver includes permutation of several processes including gene duplication, fusion, neofunctionalization, and deletion resulting in the present chemotaxonomy. A critical gene fusion event resulting in the key bifunctional enzyme reticuline epimerase (REPI), which catalyzes the stereochemical inversion of (S)-reticuline, was suggested to precede neofunctionalization of downstream enzymes leading to morphine biosynthesis in opium poppy (Papaver somniferum). The ancestrally related aldo-keto reductases 1,2-dehydroreticuline reductase (DRR), which occurs in some species as a component of REPI, and codeinone reductase (COR) catalyze the second and penultimate steps, respectively, in the pathway converting (S)-reticuline to morphine. Orthologs for each enzyme isolated from the transcriptomes of 12 Papaver species were shown to catalyze their respective reactions in species that capture states of the metabolic pathway prior to key evolutionary events, including the gene fusion event leading to REPI, thus suggesting a patchwork model for pathway evolution. Analysis of the structure and substrate preferences of DRR orthologs in comparison with COR orthologs revealed structure-function relationships underpinning the functional latency of DRR and COR orthologs in the genus Papaver, thus providing insights into the molecular events leading to the evolution of the pathway.

The Role of Liriodendron Dof Gene Family in Abiotic Stress Response

The DOF (DNA-binding with one finger) transcription factors are exclusive to plants and play crucial roles in plant growth, development, and environmental adaptation. Although extensive research has been conducted on the Dof gene family in Arabidopsis, maize, and Solanum, investigations concerning the role of this gene family in Liriodendron remain unreported, leaving its biological function largely unknown. In this study, we performed a comprehensive genome-wide identification of the Dof gene family based on the Liriodendron genome, resulting in the discovery of a total of 17 LcDof gene members. Based on the results of phylogenetic analysis, the 17 LcDof proteins were classified into eight subfamilies. The motif analysis revealed the diverse nature of motifs within the D1 subfamily, which includes a distinct type of Dof transcription factor known as CDF (Cycling Dof Factor). We further characterized the chromosomal distribution, gene structure, conserved protein motifs, and cis-elements in the promoter regions. Additionally, utilizing transcriptome data from Liriodendron hybrids and conducting RT-qPCR experiments, we investigated the expression patterns of LhDofs under various abiotic stresses such as drought, cold, and heat stress. Notably, we found that several LhDofs, particularly LhDof4 and LhDof6, were significantly upregulated in response to abiotic stress. Furthermore, we cloned LhDof4 and LhDof6 genes and found that its encoding protein was mainly located in the nucleus by transient transformation in Liriodendron hybrids protoplast. Subsequently, we used LhDof6-overexpressing Liriodendron hybrid seedlings. We found that overexpression of LhDof6 enhanced the cold tolerance of the plants, increasing their survival rate at -20 °C. This result was further validated by changes in physiological indicators.

Evolutionary analysis of conserved non-coding elements subsequent to whole-genome duplication in opium poppy

Whole-genome duplication (WGD) leads to the duplication of both coding and non-coding sequences within an organism's genome, providing an abundant supply of genetic material that can drive evolution, ultimately contributing to plant adaptation and speciation. Although non-coding sequences contain numerous regulatory elements, they have been understudied compared to coding sequences. In order to address this gap, we explored the evolutionary patterns of regulatory sequences, coding sequences and transcriptomes using conserved non-coding elements (CNEs) as regulatory element proxies following the recent WGD event in opium poppy (Papaver somniferum). Our results showed similar evolutionary patterns in subgenomes of regulatory and coding sequences. Specifically, the biased or unbiased retention of coding sequences reflected the same pattern as retention levels in regulatory sequences. Further, the divergence of gene expression patterns mediated by regulatory element variations occurred at a more rapid pace than that of gene coding sequences. However, gene losses were purportedly dependent on relaxed selection pressure in coding sequences. Specifically, the rapid evolution of tissue-specific benzylisoquinoline alkaloid production in P. somniferum was associated with regulatory element changes. The origin of a novel stem-specific ACR, which utilized ancestral cis-elements as templates, is likely to be linked to the evolutionary trajectory behind the transition of the PSMT1-CYP719A21 cluster from high levels of expression solely in P. rhoeas root tissue to its elevated expression in P. somniferum stem tissue. Our findings demonstrate that rapid regulatory element evolution can contribute to the emergence of new phenotypes and provide valuable insights into the high evolvability of regulatory elements.