The centromere landscapes of four karyotypically diverse Papaver species provide insights into chromosome evolution and speciation.

Understanding the roles played by centromeres in chromosome evolution and speciation is complicated by the fact that centromeres comprise large arrays of tandemly repeated satellite DNA, which hinders high-quality assembly. Here, we used long-read sequencing to generate nearly complete genome assemblies for four karyotypically diverse Papaver species, P. setigerum (2n = 44), P. somniferum (2n = 22), P. rhoeas (2n = 14), and P. bracteatum (2n = 14), collectively representing 45 gapless centromeres. We identified four centromere satellite (cenSat) families and experimentally validated two representatives. For the two allopolyploid genomes (P. somniferum and P. setigerum), we characterized the subgenomic distribution of each satellite and identified a "homogenizing" phase of centromere evolution in the aftermath of hybridization. An interspecies comparison of the peri-centromeric regions further revealed extensive centromere-mediated chromosome rearrangements. Taking these results together, we propose a model for studying cenSat competition after hybridization and shed further light on the complex role of the centromere in speciation.

Developing benzylisoquinoline alkaloid-enriched opium poppy via CRISPR-directed genome editing: A review.

Among plant-derived secondary metabolites are benzylisoquinoline alkaloids (BIAs) that play a vital role in medicine. The most conspicuous BIAs frequently found in opium poppy are morphine, codeine, thebaine, papaverine, sanguinarine, and noscapine. BIAs have provided abundant clinically useful drugs used in the treatment of various diseases and ailments With an increasing demand for these herbal remedies, genetic improvement of poppy plants appears to be essential to live up to the expectations of the pharmaceutical industry. With the advent of clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated9 (Cas9), the field of metabolic engineering has undergone a paradigm shift in its approach due to its appealing attributes, such as the transgene-free editing capability, precision, selectivity, robustness, and versatility. The potentiality of the CRISPR system for manipulating metabolic pathways in opium poppy was demonstrated, but further investigations regarding the use of CRISPR in BIA pathway engineering should be undertaken to develop opium poppy into a bioreactor synthesizing BIAs at the industrial-scale levels. In this regard, the recruitment of RNA-guided genome editing for knocking out miRNAs, flower responsible genes, genes involved in competitive pathways, and base editing are described. The approaches presented here have never been suggested or applied in opium poppy so far.

Active-Site Oxygen Accessibility and Catalytic Loop Dynamics of Plant Aromatic Amino Acid Decarboxylases from Molecular Simulations.

Aromatic amino acid decarboxylases (AAADs) are pyridoxal-5'-phosphate (PLP)-dependent enzymes that catalyze the decarboxylation of aromatic amino acid l-amino acids. In plants, apart from canonical AAADs that catalyze the straightforward decarboxylation reaction, other members of the AAAD family function as aromatic acetaldehyde synthases (AASs) and catalyze more complex decarboxylation-dependent oxidative deamination. The interconversion between a canonical AAAD and an AAS can be achieved by a single tyrosine-phenylalanine mutation in the large catalytic loop of the enzymes. In this work, we report implicit ligand sampling (ILS) calculations of the canonical l-tyrosine decarboxylase from Papaver somniferum (PsTyDC) that catalyzes l-tyrosine decarboxylation and its Y350F mutant that instead catalyzes the decarboxylation-dependent oxidative deamination of the same substrate. Through comparative analysis of the resulting three-dimensional (3D) O2 free energy profiles, we evaluate the impact of the key tyrosine/phenylalanine mutation on oxygen accessibility to both the wild type and Y350F mutant of PsTyDC. Additionally, using molecular dynamics (MD) simulations of the l-tryptophan decarboxylase from Catharanthus roseus (CrTDC), we further investigate the dynamics of a large catalytic loop known to be indispensable to all AAADs. Results of our ILS and MD calculations shed new light on how key structural elements and loop conformational dynamics underlie the enzymatic functions of different members of the plant AAAD family.

Phylogeography and genetic structure of Papaver bracteatum populations in Iran based on genotyping-by-sequencing (GBS).

Papaver bracteatum, known for its high thebaine content and absence of morphine, has emerged as a promising alternative to opium poppy for codeine production. In this study, our objective was to create a diverse panel representing the natural variation of this species in Iran. To achieve this, we employed genotyping-by-sequencing to obtain genome-wide distributed single-nucleotide polymorphisms (SNPs) for phylogeographic analysis, population structure assessment, and evaluation of genetic diversity within P. bracteatum populations. A total of 244 P. bracteatum individuals from 13 distinct populations formed seven genetic groups, along with one highly admixed population. We observed a clear split between the populations inhabiting the Alborz Mts. in the east and Zagros Mts. in the west. In between these mountain ranges, the population of Kachal Mangan exhibited a high degree of genetic admixture between both genetic groups. At or after the end of the last glacial maximum, when climate conditions rapidly changed, all P. bracteatum populations experienced a strong demographic bottleneck reducing the already small effective population sizes further before they increased to their recent strengths. Our results suggest that the ongoing climate change together with human pressure on the species' habitats and limited seed-dispersal ability are potential factors contributing today to rising genetic isolation of P. bracteatum populations. Our results provide genetic data that can be used for conservation measures to safeguard the species' genetic diversity as a resource for future breeding approaches in this medicinally important species.

Development of a two-layer machine learning model for the forensic application of legal and illegal poppy classification based on sequence data.

Poppies are beneficial plants with a variety of applications, including medicinal, edible, ornamental, and industrial purposes. Some Papaver species are forensically significant plants because they contain opium, a narcotic substance. Internationally trafficked species of illegal poppies are being identified by DNA barcoding employing multiple markers in response to their forensic value. However, effective markers for precise species identification of legal and illegal poppies are still under discussion, with research on illegal poppies focusing on Papaver somniferum L., and species identification studies of Papaver bracteatum and Papaver setigerum DC. still lacking. As a result, in order to evaluate the performance of genetic markers and classify their DNA sequences in the genus Papaver, this study developed the first machine learning-based two-layer model, in which the first layer classifies legal and illegal poppies from the given sequence and the second layer identifies species of illegal poppies using their sequences. We constructed the dataset and investigated biological features from four markers, internal transcribed spacer 1 (ITS1), internal transcribed spacer 2 (ITS2), transfer RNA Leucine (trnL), transfer RNA Leucine - transfer RNA Phenylalanine intergenic spacer (trnL-trnF intergenic spacer) and their combination, using four machine learning algorithms, K-nearest neighbor (KNN), Naïve Bayes (NB), extreme gradient boost (XGBoost) and Random Forest (RF). According to our findings, for Layer 1 to classify legal and illegal poppies, KNN-based models using combined ITS region achieved the greatest performance of accuracy 0.846 and 0.889 using training and test sets, respectively. Additionally, for Layer 2 to identify illegal poppy species, KNN-based models using combined ITS region achieved the best performance of 0.833 and 1.000 for using training and test sets, respectively. To validate the model, the combined ITS region, which includes ITS 1 and 2 sequences, from blind poppy samples were used as a case study, with the Layer 1 correctly classifying legal and illegal poppies with over 0.830 accuracy. Layer 2 correctly identified P. setigerum DC., however, only one of the three P. somniferum L. species was accurately identified. Nevertheless, our research shows that machine learning can be used to classify and identify legal and illegal poppy species using DNA barcodes which can then be used as an efficient and effective forensic tool for improved law enforcement and a safer society.

The effects of green and chemically-synthesized copper oxide nanoparticles on the production and gene expression of morphinan alkaloids in Oriental poppy.

Oriental poppy (Papaver orientale L.) belonging to the Papaveraceae family, has the capacity to synthesize a wide range of benzylisoquinoline alkaloids (BIAs). This experiment was conducted to investigate the effects of green and chemical copper oxide nanoparticles (CuO NPs) elicitors on oxidative stress and the BIAs biosynthesis pathway in the cell suspension culture of P. orientale. This research shows that both green and chemical CuO NPs at concentrations of 20 mg/L and 40 mg/L, induce oxidative stress in the cell suspension of P. orientale by increasing the production of H2O2 and the activity of antioxidant enzymes. The comparison of treatments revealed that utilizing a lower concentration of CuO NPs (20 mg/L) and extending the duration of cell suspension incubation (up to 48 h) play a more influential role in inducing the expression of the BIAs biosynthesis pathway genes (PsWRKY, TYDC, SalSyn, SalR, SalAT, T6ODM, COR and CODM) and increasing the production of morphinan alkaloids (thebaine, codeine, and morphine). The overarching results indicate that the concentration of CuO NPs and the duration of cell treatment have a more significant impact than the nature of CuO NPs in inducing oxidative stress and stimulating the expression of the BIAs pathway genes.

Prevalence, spatial structure and evolution of resistance to acetolactate-synthase (ALS) inhibitors and 2,4-D in the major weed Papaver rhoeas (L.) assessed using a massive, country-wide sampling.

BACKGROUND: Corn poppy (Papaver rhoeas) is the most damaging broadleaf weed in France. Massively parallel amplicon sequencing was used to investigate the prevalence, mode of evolution and spread of resistance-endowing ALS alleles in 422 populations randomly sampled throughout poppy's range in France. Bioassays were used to detect resistance to the synthetic auxin 2,4-D in 43 of these populations. RESULTS: A total of 21 100 plants were analysed and 24 mutant ALS alleles carrying an amino-acid substitution involved or potentially involved in resistance were identified. The vast majority (97.6%) of the substitutions occurred at codon Pro197, where all six possible single-nucleotide non-synonymous substitutions plus four double-nucleotide substitutions were identified. Changes observed in the enzymatic properties of the mutant ALS isoforms could not explain the differences in prevalence among the corresponding alleles. Sequence read analysis showed that mutant ALS alleles had multiple, independent evolutionary origins, and could have evolved several times independently within an area of a few kilometres. Finally, 2,4-D resistance was associated with mutant ALS alleles in individual plants in one third of the populations assayed. CONCLUSION: The intricate geographical mosaic of mutant ALS alleles observed is the likely result of the combination of huge population sizes, multiple independent mutation events and human-mediated spread of resistance. Our work highlights the ability of poppy populations and individual plants to accumulate different ALS alleles and as yet unknown mechanisms conferring resistance to synthetic auxins. This does not bode well for the continued use of chemical herbicides to control poppy. © 2023 Society of Chemical Industry.

Engineering of CYP82Y1, a cytochrome P450 monooxygenase: a key enzyme in noscapine biosynthesis in opium poppy.

Protein engineering provides a powerful base for the circumvention of challenges tied with characteristics accountable for enzyme functions. CYP82Y1 introduces a hydroxyl group (-OH) into C1 of N-methylcanadine as the substrate to yield 1-hydroxy-N-methylcanadine. This chemical process has been found to be the gateway to noscapine biosynthesis. Owning to the importance of CYP82Y1 in this biosynthetic pathway, it has been selected as a target for enzyme engineering. The insertion of tags to the N- and C-terminal of CYP82Y1 was assessed for their efficiencies for improvement of the physiological performances of CYP82Y1. Although these attempts achieved some positive results, further strategies are required to dramatically enhance the CYP82Y1 activity. Here methods that have been adopted to achieve a functionally improved CYP82Y1 will be reviewed. In addition, the possibility of recruitment of other techniques having not yet been implemented in CYP82Y1 engineering, including the substitution of the residues located in the substrate recognition site, formation of the synthetic fusion proteins, and construction of the artificial lipid-based scaffold will be discussed. Given the fact that the pace of noscapine synthesis is constrained by the CYP82Y1-catalyzing step, the methods proposed here are capable of accelerating the rate of reaction performed by CYP82Y1 through improving its properties, resulting in the enhancement of noscapine accumulation.

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.

Contrasting self-recognition rejection systems for self-incompatibility in Brassica and Papaver.

Self-incompatibility (SI) plays a pivotal role in whether self-pollen is accepted or rejected. Most SI systems employ two tightly linked loci encoding highly polymorphic pollen (male) and pistil (female) S-determinants that control whether self-pollination is successful or not. In recent years our knowledge of the signalling networks and cellular mechanisms involved has improved considerably, providing an important contribution to our understanding of the diverse mechanisms used by plant cells to recognise each other and elicit responses. Here, we compare and contrast two important SI systems employed in the Brassicaceae and Papaveraceae. Both use 'self-recognition' systems, but their genetic control and S-determinants are quite different. We describe the current knowledge about the receptors and ligands, and the downstream signals and responses utilized to prevent self-seed set. What emerges is a common theme involving the initiation of destructive pathways that block the key processes that are required for compatible pollen-pistil interactions.

Rational Engineering of (S)-Norcoclaurine Synthase for Efficient Benzylisoquinoline Alkaloids Biosynthesis.

(S)-Norcoclaurine is synthesized in vivo through a metabolic pathway that ends with (S)-norcoclaurine synthase (NCS). The former constitutes the scaffold for the biosynthesis of all benzylisoquinoline alkaloids (BIAs), including many drugs such as the opiates morphine and codeine and the semi-synthetic opioids oxycodone, hydrocodone, and hydromorphone. Unfortunately, the only source of complex BIAs is the opium poppy, leaving the drug supply dependent on poppy crops. Therefore, the bioproduction of (S)-norcoclaurine in heterologous hosts, such as bacteria or yeast, is an intense area of research nowadays. The efficiency of (S)-norcoclaurine biosynthesis is strongly dependent on the catalytic efficiency of NCS. Therefore, we identified vital NCS rate-enhancing mutations through the rational transition-state macrodipole stabilization method at the Quantum Mechanics/Molecular Mechanics (QM/MM) level. The results are a step forward for obtaining NCS variants able to biosynthesize (S)-norcoclaurine on a large scale.

Structural Characterization and Molecular Dynamics Study of the REPI Fusion Protein from Papaver somniferum L.

REPI is a pivotal point enzyme in plant benzylisoquinoline alkaloid metabolism as it promotes the evolution of the biosynthetic branch of morphinan alkaloids. Experimental studies of its activity led to the identification of two modules (DRS and DRR) that catalyze two sequential steps of the epimerization of (S)- to (R)-reticuline. Recently, special attention has been paid to its genetic characterization and evolutionary history, but no structural analyses of the REPI protein have been conducted to date. We present here a computational structural characterization of REPI with heme and NADP cofactors in the apo state and in three complexes with substrate (S)-reticuline in DRS and intermediate 1,2-dehydroreticuline in DRS and in DRR. Since no experimental structure exists for REPI, we used its AlphaFold model as a scaffold to build up these four systems, which were submitted to all-atom molecular dynamics (MD) simulations. A comparison of MD results for the four systems revealed key dynamic changes associated with cofactor and ligand binding and provided a dynamic picture of the evolution of their structures and interactions. We also explored the possible dynamic occurrence of tunnels and electrostatic highways potentially involved in alternative mechanisms for channeling the intermediate from DRS to DRR.

Alkaloid binding to opium poppy major latex proteins triggers structural modification and functional aggregation.

Opium poppy accumulates copious amounts of several benzylisoquinoline alkaloids including morphine, noscapine, and papaverine, in the specialized cytoplasm of laticifers, which compose an internal secretory system associated with phloem throughout the plant. The contiguous latex includes an abundance of related proteins belonging to the pathogenesis-related (PR)10 family known collectively as major latex proteins (MLPs) and representing at least 35% of the total cellular protein content. Two latex MLP/PR10 proteins, thebaine synthase and neopione isomerase, have recently been shown to catalyze late steps in morphine biosynthesis previously assigned as spontaneous reactions. Using a combination of sucrose density-gradient fractionation-coupled proteomics, differential scanning fluorimetry, isothermal titration calorimetry, and X-ray crystallography, we show that the major latex proteins are a family of alkaloid-binding proteins that display altered conformation in the presence of certain ligands. Addition of MLP/PR10 proteins to yeast strains engineered with morphine biosynthetic genes from the plant significantly enhanced the conversion of salutaridine to morphinan alkaloids.

Computational insights into diverse aspects of glutathione S-transferase gene family in Papaver somniferum.

Plant glutathione S-transferases are an ancient protein superfamily having antioxidant activity. These proteins are primarily involved in diverse plant functions such as plant growth and development, secondary metabolism, signaling pathways and defense against biotic and abiotic stresses. The current study aimed to comprehensively identify and characterize the GST gene family in the medicinally important crop Papaver somniferum. A total of 93 GST proteins were identified belonging to eight GST classes and found to be majorly localized in the cytoplasm. All GST genes were found on eleven opium chromosomes. Gene duplication analysis showed segmental duplication as a key factor for opium GST gene family expansion under strong purifying selection. Phylogenetic analysis with gymnosperm, angiosperm and bryophyte revealed the evolution of GSTs earlier than their division into separate groups and also prior to the divergence of monocot and dicot. The secondary structure prediction showed the dominance of α-helices indicative of PsomGSTs as structurally stable and elastic proteins. Gene architecture showed the conservation of number of exons across the classes. MEME analysis revealed only a few class specific and many across class conserved motifs. Ser was found to be the active site residue of tau, phi, theta and zeta class and Cys was catalytic residue of DHAR, lambda and GHR class. Promoter analyses identified many cis-acting regulatory elements related to hormonal, cellular, stress and light response functions. Ser was the key phosphorylation site. Only three glycosylation sites were found across the 93 PsomGSTs. 3D structure prediction was also performed and was validated. Interactome analyses revealed the correlation of PsomGSTs with glutathione metabolizing proteins. Gene enrichment analysis and KEGG pathway analyzed the involvement of PsomGSTs in three major pathways i.e. glutathione metabolism, tyrosine metabolism and ascorbate metabolism. The outcome revealed high model quality of PsomGSTs. The results of the current study will be of potential significance to understand the functional and structural importance of the GST gene family in opium, a medicinally important crop.

Development of diagnostic SNP markers and a novel SNP genotyping assay for distinguishing opium poppies.

The opium poppy acts as an important natural pain reliever but is also responsible for increased rates of severe drug abuse and addiction owing to its characteristic psychoactive effect. Non-medical illicit use of the poppy plant is markedly increasing worldwide, thereby highlighting the need for a robust species identification strategy. In this study, we identified SNPs within the region of two universal DNA barcodes, matK (maturase K) and the trnL-trnF (tRNA-Leu [3'exon]-tRNA-Phe [exon] intergenic spacer, that are forensically applicable for distinguishing opium poppy species based on a genetic analysis of 164 samples of family Papaveraceae obtained from locations spanning Jeolla-do and Jeju Island, Republic of Korea. A comparative analysis of the DNA barcode sequences for two narcotic types of the Papaver species (Papaver somniferum, Papaver somniferum subs. setigerum) to eight non-narcotic species revealed three unique nucleotide substitution events. Newly identified SNPs were located at position 255 of matK and at positions 305 and 306 of trnL-trnF; the narcotic species contained C, A, and T, whereas non-narcotic species contained T, G, and C at these positions. Phylogenetic analysis demonstrated that newly identified SNPs, which we named PsMAT255 and PsLF305/306, could be used to clearly differentiate between the narcotic and non-narcotic types of Papaver species based on the patterns of nucleotide variation. These results indicate that the nucleotide differences between the narcotic and non-narcotic species may influence genetic markers. We, therefore, developed a novel SNP-based allelic genotyping assay using the RT-PCR system that can reliably differentiate the narcotic type of the Papaver species. In summary, our findings suggest that the newly identified species-specific SNPs of both matK and trnL-trnF can be used as identification markers of narcotic Papaver species. Furthermore, a newly developed TaqMan allelic discrimination assay may be used as a practically applicable diagnostic method to survey several illicit narcotic specimens carrying the type-specific SNP.

A functionally conserved STORR gene fusion in Papaver species that diverged 16.8 million years ago.

The STORR gene fusion event is considered essential for the evolution of the promorphinan/morphinan subclass of benzylisoquinoline alkaloids (BIAs) in opium poppy as the resulting bi-modular protein performs the isomerization of (S)- to (R)-reticuline essential for their biosynthesis. Here, we show that of the 12 Papaver species analysed those containing the STORR gene fusion also contain promorphinans/morphinans with one important exception. P. californicum encodes a functionally conserved STORR but does not produce promorphinans/morphinans. We also show that the gene fusion event occurred only once, between 16.8-24.1 million years ago before the separation of P. californicum from other Clade 2 Papaver species. The most abundant BIA in P. californicum is (R)-glaucine, a member of the aporphine subclass of BIAs, raising the possibility that STORR, once evolved, contributes to the biosynthesis of more than just the promorphinan/morphinan subclass of BIAs in the Papaveraceae.

Development of SSR and SNP markers for identifying opium poppy.

Papaver somniferum L. (opium poppy) is the original plant of heroin, which is a major narcotic drug, and this plant has brought great harm to human health. However, the ban on opium poppy cultivation and trafficking is facing great challenges because of abnormal profits. Therefore, rapid and accurate identification is important to address the abovementioned problems. In this study, eleven simple sequence repeats (SSR) markers and two single nucleotide polymorphism (SNP) markers were mined to distinguish opium poppy from other six Papaver species. These molecular markers were further verified through a large number of plant materials of these seven Papaver species. An excellent multiplex polymerase chain reaction (PCR) system that simultaneously amplifies the three of eleven SSR markers was developed, which effectively improves the efficiency and speed of identification. The present research is of great implication for identifying and investigating the illegal cultivation and trafficking of opium poppy.

Self-incompatibility requires GPI anchor remodeling by the poppy PGAP1 ortholog HLD1.

Glycosylphosphatidylinositol-anchored proteins (GPI-APs) are tethered to the outer leaflet of the plasma membrane where they function as key regulators of a plethora of biological processes in eukaryotes. Self-incompatibility (SI) plays a pivotal role regulating fertilization in higher plants through recognition and rejection of "self" pollen. Here, we used Arabidopsis thaliana lines that were engineered to be self-incompatible by expression of Papaver rhoeas SI determinants for an SI suppressor screen. We identify HLD1/AtPGAP1, an ortholog of the human GPI-inositol deacylase PGAP1, as a critical component required for the SI response. Besides a delay in flowering time, no developmental defects were observed in HLD1/AtPGAP1 knockout plants, but SI was completely abolished. We demonstrate that HLD1/AtPGAP1 functions as a GPI-inositol deacylase and that this GPI-remodeling activity is essential for SI. Using GFP-SKU5 as a representative GPI-AP, we show that the HLD1/AtPGAP1 mutation does not affect GPI-AP production and targeting but affects their cleavage and release from membranes in vivo. Our data not only implicate GPI-APs in SI, providing new directions to investigate SI mechanisms, but also identify a key functional role for GPI-AP remodeling by inositol deacylation in planta.

Interindividual variability in diabetic patients' response to opium poppy: an overview of impressive factors.

Diabetic patients always seek alternative treatments to lower their blood glucose level efficiently, because antidiabetic drugs produce adverse effects and many patients experience reduced response after a treatment period. Opium poppy (Papaver somniferum) is frequently consumed by diabetic patients for reduction of blood glucose level. Scientific studies found controversial results in the investigation of the blood glucose-lowering effects of opium poppy. In this regard, we explored the antidiabetic effect of opium poppy more closely. The antidiabetic or antihyperglycemic effect of P. somniferum alkaloids were reviewed. Next, opioid receptors and their role in diabetes were explored. In the final part origins of interindividual variabilities in opioid receptors and metabolizing enzymes' functions including genetic and epigenetic factors were reviewed.