The disclosure relates to transcription cassettes and vectors comprising a nucleotide sequence encoding a transcription factor and wherein the transcription factor controls expression of genes that encode polypeptides involved in the production of plant benzylisoquinoline (BIA) alkaloids such as for example, morphine, codeine, oripavine and thebaine, the disclosure further relates to plants, plant cells/cultures and seeds that are genetically engineered to include a recombinant copy or copies of the nucleic acid encoding said transcription factor to increase expression of said transcription factor thereby modifying BIA manufacture processes for the cultivation of the genetically engineered plants or plant cell cultures and for the modulation of BIAs and processes for the extraction of BIAs from the genetically engineered plants or plant cell cultures.
The opium poppy, Papaver somniferum is an important source of a variety of BIAs. Due to their narcotic and analgesic properties, some BIAs, and their derivatives, are desired for use in therapy. P. somniferum is a source of clinically useful alkaloids such as morphine, codeine, thebaine, noscapine and papaverine. BIAs are extracted from latex harvested from the green seed pods of opium poppy or from the poppy straw which is the dried mature plant. The pathway to produce alkaloids and the various genes involved in the pathway are known and are disclosed in U.S. application Ser. Nos. 15/182,761 and 15/304,455 the contents of which are hereby incorporated by reference in their entirety. Morphinan alkaloids such as codeine and morphine are known to derive from the intermediate (R)-reticuline. (R)-reticuline is thought to be formed by its enantiomer (S)-reticuline in a two-step isomerization process. (R)-reticuline is then further transformed to thebaine. Thebaine is transformed either to oripavine to morphinone and morphine or via an alternative route to codeinone and codeine which is then subsequently transformed to morphine.
Various attempts to increase the content of alkaloids in planta have been attempted by natural breeding methods or reverse genetics shifting the alkaloid pathway to increase intermediate and end-products such as for example obtaining plants with increased thebaine or codeine content; see WO2016/207643; WO2017/112011; EP3398430. Plants comprising knock out mutations are disclosed in patent application No U.S. Ser. No. 14/375,120 the content of which is hereby incorporated by reference. Alternative methods to modulate plant metabolism are described in WO2021/026119 by overexpressing a transcription factor that positively regulates nicotine biosynthesis. Modified Nicotiana tabacum plants comprised increased nicotine levels in tobacco leaves.
Fast Neutron Mutagenesis (FNM) was carried out on seed of a morphine & noscapine producing P. sominferum cultivar. The mutagenized M1 seed was sown and the M1 generation self-pollinated to generate M2 seed. The M2 generation was then screened for any unusual metabolite profiles compared to the non-mutagenized parental variety. The screen identified two independent mutants that no longer accumulate any BIAs and no longer express most BIA biosynthesis genes. Molecular characterisation of the mutants at the DNA level revealed that each mutant carries a unique deletion spanning a shared region encoding a transcription factor. Because of the loss of BIA synthesis gene expression and lack of BIA accumulation in the mutants, the transcription factor was named REGULATOR OF BENZYLISOQUINOLINE ALKALOIDS (RBA) and the mutants retrospectively labelled rba-1 and rba-2 mutants.
This disclosure characterises a transcription factor that controls the expression of genes involved in the production of BIAs such as morphinan alkaloids.
According to an aspect of the invention there is provided an expression cassette adapted for plant expression comprising a nucleic acid molecule selected from:
Hybridization of a nucleic acid molecule occurs when two complementary nucleic acid molecules undergo an amount of hydrogen bonding to each other. Isolated nucleic acid molecules as referred herein include genomic DNA, cDNA molecules and RNA molecules. The stringency of hybridization can vary according to the environmental conditions surrounding the nucleic acids, the nature of the hybridization method, and the composition and length of the nucleic acid molecules used. Calculations regarding hybridization conditions required for attaining particular degrees of stringency are discussed in Sambrook et al., Molecular Cloning: A Laboratory Manual (Cold Spring Harbor Laboratory Press, Cold Spring Harbor, N Y, 2001); and Tijssen, Laboratory Techniques in Biochemistry and Molecular Biology-Hybridization with Nucleic Acid Probes Part I, Chapter 2 (Elsevier, New York, 1993). The Tm is the temperature at which 50% of a given strand of a nucleic acid molecule is hybridized to its complementary strand. The following is an exemplary set of hybridization conditions and is not limiting:
In a preferred embodiment of the invention said transcription factor is a helix loop helix transcription factor.
In a preferred embodiment of the invention said isolated nucleic acid molecule is at least 50% identical to the nucleotide sequence set forth in SEQ ID NO: 1.
Preferably, the isolated nucleic acid molecule is at least 55%, 60%, 65%, 70%, 75%, 80%, 90%, 91%, 92%, 93%, 94%, 95%, 96%, 97%, 98% or 99% identical to the nucleotide sequence set forth in SEQ ID NO: 1 over the full-length nucleotide sequence.
In a preferred embodiment of the invention said nucleic acid molecule comprises or consists of sequence ID NO 1, or polymorphic sequence variant thereof.
In a preferred embodiment said nucleic acid encodes a polypeptide comprising or consisting of an amino acid sequence set forth in SEQ ID NO 2, or polymorphic sequence variant thereof.
A modified polypeptide as herein disclosed may differ in amino acid sequence by one or more substitutions, additions, deletions, truncations that may be present in any combination and includes polymorphic sequence variants that the skilled person would expect to exist in nature. Among preferred variants are those that vary from a reference polypeptide by conservative amino acid substitutions. Such substitutions are those that substitute a given amino acid by another amino acid of like characteristics. The following non-limiting list of amino acids are considered conservative replacements (similar): a) alanine, serine, and threonine; b) glutamic acid and aspartic acid; c) asparagine and glutamine d) arginine and lysine; e) isoleucine, leucine, methionine and valine and f) phenylalanine, tyrosine and tryptophan. Most highly preferred are variants that retain or enhance the same biological function and activity as the reference polypeptide from which it varies.
In a preferred embodiment of the invention said modified polypeptide is a variant and is at least 50%, 55%, 59%, 60%, 65%, 70%, 75%, 80%, 85%, 90%, 95%, 96%, 97%, 98% or 99% similar to the amino acid sequence set forth in SEQ ID NO: 2 over the full amino acid sequence.
In a preferred embodiment of the invention said modified polypeptide is a variant and is at least 50%, 55%, 59%, 60%, 65%, 70%, 75%, 80%, 85%, 90%, 95%, 96%, 97%, 98% or 99% identical to the amino acid sequence set forth in SEQ ID NO: 2 over the full amino acid sequence.
In a preferred embodiment said one or more genes involved in the synthesis of benzylisoquinoline alkaloids are selected from: thebaine synthase, salutaridine reductase, salutaridinol 7-O-acetyltransferase, salutaridine synthase, codeine O-demethylase, codeinone reductase, thebaine 6-O-demethylase, (S)-to-(R)-reticuline P450-oxidoreductase, O-methyltransferase 1, canadine synthase, O-methyltransferase 3, cytochrome P450 CYP82Y1, O-methyltransferase 2, acetyltransferase 1, cytochrome P450 CYP82X2, cytochrome P450 CYP82X1, carboxylesterase 1, short-chain dehydrogenase/reductase, (S)-tetrahydroprotoberberine-cis-N-methyltransferase berberine bridge enzyme, 3′-hydroxy-N-methylcoclaurine 4′-O-methyltransferase 2, 3′-hydroxy-N-methylcoclaurine 4′-O-methyltransferase 1, (S)—N-methylcoclaurine 3′-hydroxylase, coclaurine N-methyltransferase, norcoclaurine 6-O-methyltransferase, S-norcoclaurine synthase 1, protopine 6-hydroxylase, norreticuline-7-O-methyltransferase, 3′-hydroxy-N-methylcoclaurine 4′-O-methyltransferase 2,3′-hydroxy-N-methylcoclaurine 4′-O-methyltransferase 1, (S)—N-methylcoclaurine 3′-hydroxylase, coclaurine N-methyltransferase, norcoclaurine 6-O-methyltransferase, S-norcoclaurine synthase 1, (S)-tetrahydroprotoberberine-cis-N-methyltransferase, norreticuline-7-O-methyltransferase, scoulerine O-demethylase.
In a preferred embodiment said one or more genes involved in the synthesis of benzylisoquinoline alkaloids are selected from the genes listed in Table 2.
In a preferred embodiment of the invention said nucleic acid molecule is operably linked to a transcription promoter.
In an alternative preferred embodiment of the invention said transcription promoter is a heterologous promoter.
In an alternative preferred embodiment said transcription promoter is an endogenous promoter that naturally controls transcription of said transcription factor.
In a preferred embodiment of the invention said transcription promoter is an inducible promoter.
In an alternative preferred embodiment of the invention said transcription promotor is a constitutive promoter.
In an alternative preferred embodiment of the invention said promoter is a tissue specific promoter.
In a further preferred method of the invention said promotor is selected from the group consisting of CaMV 35S promoter and Glycine Max Ubiquitin 3 promoter.
By “transcription promoter” is meant a nucleotide sequence upstream from the transcriptional initiation site and which contains all the regulatory regions required for transcription. Said promoters include viral, fungal, bacterial, animal and plant-derived promoters capable of functioning in plant cells. Constitutive promoters include, for example CaMV 35S promoter (Odell et al. (1985) Nature 313, 9810-812); rice actin (McElroy et al. (1990) Plant Cell 2: 163-171); ubiquitin (Christian et al. (1989) Plant Mol. Biol. 18 (675-689); pEMU (Last et al. (1991) Theor Appl. Genet. 81: 581-588); MAS (Velten et al. (1984) EMBO J. 3. 2723-2730); ALS promoter (U.S. application Ser. No. 08/409,297), and the like. Other constitutive promoters include those in U.S. Pat. Nos. 5,608,149; 5,608,144; 5,604,121; 5,569,597; 5,466,785; 5,399,680, 5,268,463; and 5,608,142, each of which is incorporated by reference.
Chemical-regulated promoters can be used to modulate the expression of a gene in a plant through the application of an exogenous chemical regulator. Depending upon the objective, the promoter may be a chemical-inducible promoter, where application of the chemical induced gene expression, or a chemical-repressible promoter, where application of the chemical represses gene expression. Chemical-inducible promoters are known in the art and include, but are not limited to, the maize In2-2 promoter, which is activated by benzene sulphonamide herbicide safeners, the maize GST promoter, which is activated by hydrophobic electrophilic compounds that are used as pre-emergent herbicides, and the tobacco PR-1a promoter, which is activated by salicylic acid. Other chemical-regulated promoters of interest include steroid-responsive promoters (see, for example, the glucocorticoid-inducible promoter in Schena et al. (1991) Proc. Natl. Acad. Sci. USA 88: 10421-10425 and McNellis et al. (1998) Plant J. 14(2): 247-257) and tetracycline-inducible and tetracycline-repressible promoters (see, for example, Gatz et al. (1991) Mol. Gen. Genet. 227: 229-237, and U.S. Pat. Nos. 5,814,618 and 5,789,156, herein incorporated by reference.
Where enhanced expression in particular tissues is desired, tissue-specific promoters can be utilised. Tissue-specific promoters include those described by Yamamoto et al. (1997) Plant J. 12(2): 255-265; Kawamata et al. (1997) Plant Cell Physiol. 38(7): 792-803; Hansen et al. (1997) Mol. Gen. Genet. 254(3): 337-343; Russell et al. (1997) Transgenic Res. 6(2): 157-168; Rinehart et al. (1996) Plant Physiol. 112(3): 1331-1341; Van Camp et al. (1996) Plant Physiol. 112(2): 525-535; Canevascni et al. (1996) Plant Physiol. 112(2): 513-524; Yamamoto et al. (1994) Plant Cell Physiol. 35(5): 773-778; Lam (1994) Results Probl. Cell Differ. 20: 181-196; Orozco et al. (1993) Plant Mol. Biol. 23(6): 1129-1138; Mutsuoka et al. (1993) Proc. Natl. Acad. Sci. USA 90 (20): 9586-9590; and Guevara-Garcia et al (1993) Plant J. 4(3): 495-50.
“Operably linked” means functionally associated as part of the same nucleic acid molecule, suitably positioned and oriented for transcription to be initiated from the promoter. DNA operably linked to a promoter is “under transcriptional initiation regulation” of the promoter.
According to an aspect of the invention there is provided a vector comprising an expression cassette according to the invention.
In a preferred embodiment of the invention said expression vector is a transposon.
In an alternative embodiment of the invention said vector is a viral based vector.
Of interest in the present context are nucleic acid constructs which operate as plant vectors. Specific procedures and vectors previously used with wide success in plants are described by Guerineau and Mullineaux (1993) (Plant transformation and expression vectors. In: Plant Molecular Biology Labfax (Croy RRD ed) Oxford, BIOS Scientific Publishers, pp 121-148. Suitable vectors may include plant viral-derived vectors (see e.g. EP194809). If desired, selectable genetic markers may be included in the construct, such as those that confer selectable phenotypes such as resistance to herbicides (e.g. kanamycin, hygromycin, phosphinotricin, chlorsulfuron, methotrexate, gentamycin, spectinomycin, imidazolinones and glyphosate).
According to a further aspect of the invention there is provided a plant wherein said plant is transformed or transfected with a transcription cassette or expression vector according to the invention.
Preferably, said plant is of the genus Papaver spp.
Genes involved in the pathway of producing benzylisoquinoline alkaloids are known in the art and encode polypeptides with cytochrome p450 and methyltransferase activity. Deletion of one or more genes in this pathway can alter the levels of one or more benzylisoquinoline alkaloids such as the deletion of codeine 3-O-demethylase results in Papaver somniferum plants with increased content of codeine, whereas deletion of thebaine 6-O-demethylase results in increased content of thebaine and oripavine. Plants with altered benzylisoquinoline alkaloids pathways are known in the art and disclosed in, for example, PCT/GB2017/050068 and EP3398430.
In a preferred embodiment of the invention said plant is selected from: Papaver somniferum, P. setigerum, P. bracteatum, P. orientale, P. pseudo-orientale, P. lasiothrix, P. cylindricum, P. fugax, P. triniifolium.
In a preferred embodiment of the invention said plant is a Papaver somniferum plant.
Benzylisoquinoline alkaloids include oripavine, codeine, thebaine, morphine, noscapine, morphinone, codeinone, reticuline, scoulerine, papaverine, cryptopine, laudanosine, protopine, O-methylsomniferin.
In a preferred embodiment of the invention expression of said transcription factor from said transcription cassette or expression vector encoding is increased when compared to a non-transformed plant of the same species.
In the context of this invention the expression “increased expression levels” relates to increased transcription subsequently resulting in higher nucleic acid transcript levels in said transformed plant when compared to an untransformed plant. Said increased transcript levels can be measured by, for example, quantitative PCR. Increased expression can be achieved by increasing the rate of transcription and/or by increasing the copy number of genes encoding the basic helix-loop-helix transcription factor according to the invention.
In a preferred embodiment of the invention said expression levels are increased by at least 2, 3, 4, 5, or 10-fold.
In a preferred embodiment of the invention said transformed plant or plant material thereof comprises at least two copies of said nucleic acid molecule.
In a further preferred embodiment of the invention said transformed plant or plant material thereof comprises 3 or more copies of said nucleic acid molecule encoding the basic helix-loop-helix transcription factor according to the invention.
In a preferred embodiment of the invention said transformed plant or plant cell comprises increased content of one or more BIAs when compared to an untransformed Papaver plant or plant cell.
The total sum of benzylisoquinoline alkaloids weight is the sum of the weight of oripavine, codeine, thebaine, morphine, noscapine, morphinone, codeinone, reticuline, scoulerine, papaverine, cryptopine, laudanosine, protopine and O-methylsomniferine, or more preferably, the sum of the weight of codeine, morphine, thebaine, noscapine and oripavine, or even more preferably, the sum of the weight of codeine, morphine, thebaine and oripavine.
In a preferred embodiment of the invention said levels are increased by 2, 3, 4, 5 or 10-fold
According to a further aspect of the invention there is provided plant material obtained from the plant according to the invention.
Plant material in the context of this invention refers to leaves, capsules or seeds.
According to a further aspect of the invention there is provide a plant cell transformed or transfected with a transcription cassette or an expression vector according to the invention.
According to a further aspect of the invention there is provide a plant cell culture comprising a transformed plant cell according to the invention.
In a preferred embodiment of the invention said plant cell is a Papaver spp plant cell.
In a preferred embodiment of the invention said plant cell is selected from the group consisting of Papaver somniferum, P. setigerum, P. bracteatum, P. orientale, P. pseudo-orientale, P. lasiothrix, P. cylindricum, P. fugax, P. triniifolium.
In a preferred embodiment of the invention said plant cell is a Papaver somniferum cell.
According to a further aspect of the invention there is provided a bioreactor comprising a plant cell culture according to the invention.
According to a further aspect of the invention there is provided a method to produce one or more benzylisoquinoline alkaloids comprising:
According to an aspect of the invention there is provided a process for the extraction of from a Papaver plant comprising the steps:
In a preferred method of the invention said plant material comprises poppy capsule, poppy straw and/or poppy latex.
According to a further aspect of the invention there is provided a method to produce a plant that has altered expression of a polypeptide according to the invention comprising the steps of:
In a preferred embodiment of the invention said plant is of the genus Papaver spp.
In a preferred method of the invention said Papaver plant is selected from: Papaver somniferum, P. setigerum, P. bracteatum, P. orientale, P. pseudo-orientale, P. lasiothrix, P. cylindricum, P. fugax, P. triniifolium.
According to a further aspect of the invention there is provided a plant obtained by the method according to the invention.
According to a further aspect of the invention there is provided a Papaver plant, or plant part thereof, comprising a gene encoding a transcription factor comprising a nucleotide sequence set forth in SEQ ID NO: 1, or a polymorphic sequence variant of SEQ ID NO: 1 wherein said nucleotide sequence is modified wherein the modified nucleotide sequence encodes a transcription factor with reduced or undetectable transcription factor activity.
In a preferred embodiment of the invention said nucleotide sequence is modified by addition, deletion or substitution of at least one nucleotide base.
In a preferred embodiment of the invention said Papaver plant is modified wherein said modification is the deletion of all or part of the nucleotide sequence set forth in SEQ ID NO: 1, or a polymorphic sequence variant of SEQ ID NO:1.
In a preferred embodiment of the invention said plant part thereof is a Papaver capsule.
In an alternative preferred embodiment said plant part thereof is a Papaver seed.
In a preferred embodiment of the invention said Papaver plant, capsule or seed substantially lacks benzylisoquinoline alkaloids or benzylisoquinoline precursors.
In a preferred embodiment of the invention said Papaver plant, or plant part thereof, is selected from: Papaver somniferum, P. setigerum, P. bracteatum, P. orientale, P. pseudo-orientale, P. lasiothrix, P. cylindricum, P. fugax, P. triniifolium.
In a preferred embodiment of the invention said Papaver plant is Papaver somniferum.
Throughout the description and claims of this specification, the words “comprise” and “contain” and variations of the words, for example “comprising” and “comprises”, means “including but not limited to”, and is not intended to (and does not) exclude other moieties, additives, components, integers or steps. “Consisting essentially” means having the essential integers but including integers which do not materially affect the function of the essential integers.
Throughout the description and claims of this specification, the singular encompasses the plural unless the context otherwise requires. Where the indefinite article is used, the specification is to be understood as contemplating plurality as well as singularity, unless the context requires otherwise.
Features, integers, characteristics, compounds, chemical moieties or groups described in conjunction with a particular aspect, embodiment or example of the invention are to be understood to be applicable to any other aspect, embodiment or example described herein unless incompatible therewith.
An embodiment of the invention will now be described by example only and with reference to the following figures and tables:
All plant material was grown under glass in 16 hour days at the University of York horticulture facilities. The growth substrate consisted of 4 parts John Innes No. 2, 1 part Perlite and 2 parts Vermiculite.
Mother plants used for crosses were emasculated at the hook stage of flower development. At the time of anthesis emasculated flowers were fertilised with pollen from synchronously developing father plants. To prevent contaminating pollen from reaching the receptive stigmas, emasculated flowers were covered with micro perforated clear bakery bags until four days after fertilisation.
Self-pollination was ensured by covering the flowers shortly before onset of anthesis with micro perforated clear bakery bags and manually pollinating the flower at anthesis.
Fast neutron mutagenesis (FNM) was carried out by the HAS Centre for Energy Research, Radiation Protection Department, H-1121 Budapest, Konkoly Thege út 29-33, X. epulet, Hungary. About 40 g of seed of the Sun Pharmaceutical Industries (Australia) cultivar High Noscapine 4 (HN4) were exposed to 20 Gy (beam time 50 min). The M1 generation was self-pollinated. One M2 plant per M2 seed family was screened. Latex was collected from 9 weeks old M2 seedlings and analysed. M2 plants displaying an altered alkaloid composition in the latex compared to HN4 controls were grown to maturity to confirm the phenotype in dry M2 capsule material. M3 seed was collected from M2 plants with confirmed phenotypes. Five M3 plants were then grown per M2 mutant plant to confirm the heritability of the phenotype in mature dry M3 capsule material. The rba-1 and rba-2 mutants were further propagated through self-pollinations to the M5 generation with each of the successive generations being phenotyped to confirm the ‘loss of all BIAs” phenotype. All generations were grown in the glasshouse.
Latex was collected from 9 weeks old plants from cut petioles, with a single drop dispersed into 500 μL of 10% acetic acid. This was diluted 10fold in 1% acetic acid to give an alkaloid solution in 2% acetic acid for further analysis. Capsules were harvested from the same plants used for latex analysis and single capsules were ground to a fine powder in a ball mill (Model MM04, Retsch, Haan, Germany). Samples of ground poppy straw were then weighed accurately to 10±0.1 mg and extracted in 0.5 ml of a 10% acetic acid solution with gentle shaking for 1 h at room temperature. Samples were then clarified by centrifugation and a 50 μl subsample diluted 10fold in 1% acetic acid to give an alkaloid solution in 2% acetic acid for further analysis.
All solutions were analysed using a Waters Acquity UPLC system (Waters Ltd., Elstree, UK) fitted with a Waters Acquity BEH C18 column, 2.1 mm×100 mm with 1.7 micron packing. The mobile phase used a gradient profile with eluent A consisting of 10 mM ammonium bicarbonate of pH 10.2 and eluent B methanol. The mobile phase gradient conditions used are shown in the table below with a linear gradient. The flow rate was 0.5 ml per minute. The column temperature was maintained at 60° C. The injection volume was 2 μl. The eluted peaks were ionised in positive APCI mode and detected within 5 ppm mass accuracy using a Thermo LTQ-Orbitrap. The runs were controlled by Thermo Xcalibur software (Thermo Fisher Scientific Inc., Hemel Hempstead, UK).
All data analysis was carried out using the R programming language and custom scripts. Standards for morphine, oripavine, codeine, thebaine and noscapine were used to identify and quantify these alkaloids using exact mass, retention time and peak areas for the pseudomolecular ion. Other putative alkaloid peaks were quantified by their pseudomolecular ion areas relative to the thebaine standard. For putative alkaloids, the Bioconductor rcdk package (Guha, (2007) J. Stat. Software 6(18)) was used to generate pseudomolecular formulae from exact masses within elemental constraints C=1-100, H=1-200, O=0-200, N=0-3 and mass accuracy <5 ppm. The hit with the lowest ppm error within these constraints was used to assign a putative formula.
Three weeks old rba-1 (M5 generation) and rba-2 (seM4 generation) mutant seedlings were used for DNA extraction. The seedling material was grown in complete darkness for 24 hours prior to harvest. Seedlings were severed from their primary root, flash frozen in liquid nitrogen and stored at −80° C. until shipment on dry ice to Amplicon Express, Inc. (Pullman, WA, USA) for DNA extraction.
High quality genomic DNA was prepared from the seedling material by Amplicon Express using their inhouse NGS-Grade genomic DNA preparation protocol optimised for long reads on various NextGenerationSequencing (NGS) platforms such as Chromium. The DNA preparations were shipped directly to HudsonAlphas (HudsonAlpha Institute for Biotechnology, 601 Genome Way, Huntsville, AL 35806) where Chromium libraries were constructed for each sample for the 10× Genomics Platform and sequenced with Illumina NovaSeq technology after quality control of the DNA samples. In total, 180 billion bases (2×91,141,870,748 and 2×91,161,893,801 respectively) were generated for each of the rba-1 and rba-2 mutant samples (named TC06 and TC07, respectively) in the resulting datasets that consist of 2×600 million (2×603,588,548 and 2×603,721,151 respectively) pair-ended 150 bases long sequence reads. Based on the 2.7 Gb size of the HN1 reference genome (ASM357369v1), this represents approximately 60 times sequencing depth for each base position for the two mutants.
The draft genomes TC06 and TC07, respectively, for the rba-1 and rba-2 mutants were assembled with the Supernova (version2.1.1), a 10× Genomics Linked-Read Diploid De Novo Assembler (https://support.10×genomics.com/de-novo-assembly/software/overview/latest/performance) on the Viking high performance computing cluster at the University of York. The TC06 (rba-1) mutant assembly has a total size of 2.6 Gb (2,581,235,059) with 67,648 scaffolds and a scaffold N50 of 836 kb (836,342), and the TC07 (rba-2) mutant assembly has a total size of 2.5Gb (2,459,139,609) with 87,225 scaffolds and a scaffold N50 of 210 kb (209,519).
The whole 51,213 annotated protein coding gene transcripts of the HN1 reference genome were used to perform BLASTN searches against both TC06 (rba-1) and TC07 (rba-2) mutant assemblies. The top scaffold hits were ordered for each gene based on the positions of each gene in the reference genome. It is expected that each scaffold would appear as top matches for a stretch of consecutive reference genes as the mutant genome is near identical to the reference genome. When a scaffold matches two such stretches but is interrupted by one or more genes in between with low matching top hits scores it was identified as a potential candidate for carrying a deleted region in the mutant genome. Two such scaffolds, TC06_scaffold_251445 in the TC06 (rba-1) assembly and TC07_scaffold_186274 in the TC07 (rba-2) assembly, were identified as potentially carrying overlapping deletions corresponding to the same region on chromosome 9 in the reference genome.
Firstly, dotplot analysis was used to confirm and locate the putative deletions for the above two scaffolds. A dotplot analysis is a graphic representation of the comparison of two sequences which identifies regions of close similarity after sequence alignment. A dot plot is a two-dimensional similarity matrix that have the two sequences being compared along the vertical and horizontal axes. For a simple visual representation of the similarity between two sequences, individual dots in the matrix are painted black if a defined length of bases are identical, so that matching sequence segments appear as runs of diagonal lines across the matrix. As a result, when two DNA sequences are completely identical, a solid diagonal line with a downward slope will appear. If two DNA sequences are in reverse complement to each other, a solid diagonal line with an upward slope will appear. When two sequences are random, no continuous line pattern would appear. Thus, a dotplot result would identify potential deletions as well as inversions between two sequences. Corresponding regions to both scaffolds (TC06_scaffold_251445 and TC07_scaffold_186274) were retrieved from the HN1 reference genome and compared to the respective scaffolds with Gepard dotplot tool, resulting in the detection of the respective deletions in the scaffolds as well as an inversion in TC06_scaffold_251445.
Cloning and Sanger-Sequencing of the RBA Gene from HN4 Genomic DNA and cDNA
Genomic DNA was extracted from the first true leaves of HN4 seedlings using he BioSprint 96 Plant kit on the BioSprint 96 Workstation (Qiagen) according to the manufacturer's protocol. Extracted DNA was quantified on the NanoDrop™ 8000 Spectrophotometer (Fisher Scientific) and normalized to 10 ng/μL.
Total RNA was extracted from HN4 stems harvested at anthesis and flash frozen liquid nitrogen using the Direct-Zol RNA miniprep kit (Zymo) according to the manufacturer's instructions. RNA concentration and purity were tested by NanoDrop (Thermo Fisher Scientific) and agarose gel electrophoresis. cDNA was synthesised using the SuperScript IV First-Strand cDNA Synthesis kit (Thermo Fisher Scientific) using Oligo d(T)20 primers, according to the manufacturer's instructions.
Full-length RBA gene sequences and coding sequences were amplified from HN4 genomic DNA and cDNA respectively using Q5 High-Fidelity DNA Polymerase (New England Biolabs). PCR reactions containing 1× Q5 High-Fidelity buffer, 2 mM dNTPs, 0.5 μM forward primer GAAGGGGTAGTGGAGTGGTAGTTG, 0.5 μM reverse primer GTGTTATACGATCATCGTTTTGC and 0.5 U Q5 DNA polymerase in a total volume of 25 μl were run on a Tetrad thermocycler (Bio-Rad) under the following conditions: 98° C. for 30 seconds, followed by 10 cycles of 98° C. for 10 seconds, 65° C. decreasing to 55° C. by 1° C. per cycle for 20 seconds, 72° C. for 1 minute, followed by 25 cycles of 98° C. for 10 seconds, 55° C. for 20 seconds, 72° C. for 1 minute, followed by 72° C. for 2 minutes, then a hold at 7° C. PCR products were checked on 1% agarose gel, extracted by Wizard SV PCR extraction kit (Promega) according to the manufacturer's instructions, and the concentration determined using the Qubit™ dsDNA BR Assay Kit (Thermo Fisher Scientific).
Purified RBA PCR product was cloned using the CloneJET PCR Cloning Kit (Thermo Fisher Scientific) according to the manufacturer's protocol and transformed into Subcloning Efficiency™ DH5a Competent Cells (Invitrogen) according to the manufacturer's protocol.
Between 2 and 12 single colonies positive for each insert were selected by colony PCR. Plasmids were extracted using the QIAprep Spin Miniprep kit (Qiagen) and used for Sanger-sequencing.
Upper stem samples (defined as the 2 cm stem section immediately underneath the flower head) were collected from rba-1 and rba-2 mutant plants of the M4 generation as well as HN4 wild type at the time of anthesis. The samples were flash frozen in liquid nitrogen and stored at −80° C. until extraction. Grinding was performed in jars chilled in liquid nitrogen on the Qiagen TissueLyser as follows: 15 seconds at 20 Hz, followed by re-chilling in liquid nitrogen, then a further 15 seconds at 20 Hz. 100 mg of ground material was used per RNA extraction. RNA was isolated from the powder using a CTAB-based extraction method (Chang et al. (1993) Plant Mol. Biol. Rep.11: 113-116) with small modifications: (i) three sequential extractions with chloroform:isoamylalcohol (24:1) were performed and (ii) the RNA was precipitated overnight with lithium chloride at −20° C. Following extraction the samples were treated with DNAse I using Ambion's DNA-free kit according to the manufacturers' protocol. After spectrophotometric quantification, equal amounts of RNA were pooled from three plants per mutant and HN4 wild type.
RNA sequencing library preparation was carried out by the Genomics Laboratory of the Bioscience Technology Facility of the Department of Biology, University of York. RNA quality was assessed by running 1 ul from each RNA pool on an Agilent Bioanalyzer RNA Nano chip. 900 ng good quality total RNA was then used per pool for mRNA sequencing library preparation using the NEB RNA Ultra Library preparation kit for Illumina in conjunction with the NEBNext Poly(A) mRNA Magnetic Isolation Module (New England BioLabs Inc.), and NEB Next single 6 bp indexing primers, according to the manufacturer's instructions. First, mRNA was purified from total RNA using two rounds of sample binding to oligo d(T)-coupled paramagnetic beads and washing. Purified mRNA was eluted from the beads into a first strand synthesis reaction buffer plus random primer mix, incubating at 94° C. for 12 minutes to fragment RNA. After addition of RNase inhibitors and ProtoScript II Reverse Transcriptase, first strand cDNA synthesis was performed by incubating at 10 minutes at 25° C., 50 minutes at 42° C. then 15 minutes at 70° C. Second strand synthesis and sample clean up were performed according to the manufacturer's guidelines, as were subsequent end preparation and adapter ligation steps. Libraries were amplified and complementary indices added in a 12 cycle PCR reaction. Following a final clean up step, the yield and size distribution of each amplified cDNA library was assessed using the Agilent High Sensitivity DNA kit with the Agilent 2100 Bioanalyzer and quantified using the Qubit with a HS dsDNA kit (Thermo Fisher Scientific). Libraries were then sent for 2×150 base paired end sequencing on a HiSeq 3000 at the University of Leeds Next Generation Sequencing Facility.
Three samples were sequenced, the number of reads were obtained as the following: rba-1, 2×28,154,370; rba-2, 2×21,681,680; and HN4 wild type,2×29,529,431. FASTQC method was used for QC analysis and Ribosomal RNA was filtered with mapping to rRNA_115_tax_silva_v1.0 downloaded from SILVA database (https://www.arb-silva.de/) using BOWTIE2 mapping software. The remaining RNA-seq reads were mapped to the reference transcript dataset of the 51,213 annotated proteins coding genes of the HN1 reference genome (ASM357369v1), using BWA mapping software with default parameters. Mapped reads were counted and retrieved using SAMTOOLS software package and the expression matrix were normalised to RPKM (Reads Per Kilobase of transcript, per Million mapped reads) values for subsequent comparative analyses.
Production of Stable Opium Poppy Transformants Constitutively Expressing RBA Under the Control of the CaMV 35S Promoter from Cauliflower Mosaic Virus
The open reading frame of RBA was cloned into binary vector pRI 201-AN (Takara Bio Inc.) that includes a 5′-UTR untranslated region (5′UTR enhancer region) from the Arabidopsis alcohol dehydrogenase downstream of the CaMV 35S promoter to drive constitutive expression of candidate genes in plants.
Plasmid pRI 201-AN was linearised by a double digest with SalI and NdeI (NEB) in rCutSmart™ buffer (NEB) for 1 hour at 37° C. according to the manufacturer's instructions. The digest was purified using the NucleoSpin™ Gel and PCR Clean-up kit (Macherey-Nagel™) according to the manufacturer's instructions to remove a small fragment from between the restriction sites.
The open reading frame (ORF) of RBA was synthesised by GeneArt Gene Synthesis (Invitrogen/Thermo Fisher Scientific UK) and further amplified using PCRBIO HiFi Polymerase kit (PCRbiosystems) according to the manufacturer's instructions using primers 470 & 471 (Table 4), which contained 5′extensions to facilitate homology-based cloning into pRI 201-AN linearised with SalI and NdeI. PCR reactions were carried out in 5×20 μL reactions on a Thermocycler as follows: initial denaturing at 95° C. for 1 minute followed by 35 cycles of 95° C. for 15 seconds, 65° C. for 15 seconds and 72° C. for 30 seconds, followed by a final extension for 1 minute at 72° C. The PCR reactions were pooled and resolved on a 0.6% agarose gel. The PCR product of the expected size was cut out from the gel and purified using NucleoSpin™ Gel and PCR Clean-up kit (Macherey-Nagel™) according to the manufacturer's instructions. This purified RBA ORF PCR product was used for homology-based cloning into pRI 201-AN.
Homology-Based Cloning of RBA into pRI 201-AN Linearised with SalI and NdeI
The pRI 201-AN_35 S::RBA expression construct (
Preparation of Competent Cells of Agrobacterium tumefaciens Strain GV3101
Agrobacterium tumefaciens strain GV3101 cells were streaked out onto YEB plates containing Gentamicin (50 μg/ml) and Rifampicin (10 μg/ml) for antibiotic selection and grown at 28° C. with shaking_for 2-3 days. A single colony was selected to inoculate a starter culture of 5 ml YEB liquid media with Gentamicin and Rifampicin selection. This was grown at 28° C. overnight, with shaking. The starter culture was used to inoculate a larger 100 ml culture with Gentamicin and Rifampicin selection that was grown for a further 3 hours at 28° C. with shaking to an OD600 of 0.5-1.0. The culture was then chilled on ice for 10 mins before the cells were collected by centrifugation. The pelleted cells were resuspended in 2 ml of 20 mM CaCl2) plus 10% glycerol and aliquoted prior to being immediately frozen in liquid nitrogen and stored at −80° C.
Transformation of A. tumefaciens Strain GV3101 with pRI 201-AN_35 S::RBA
One 50 μl aliquot of competent A. tumefaciens strain GV3101cells was defrosted on ice, and approximately 200 ng of the pRI 201-AN_35 S::RBA expression vector was added and swirled to mix. The cells were flash-frozen in liquid nitrogen, then defrosted at 37° C. in a water-bath for 5 minutes to heat shock followed by an incubation on ice for 30 minutes. Then, 250 μL of liquid LB medium was added, and the tube placed in a shaking incubator at 28° C. for 3 hours. The transformed cells were plated onto selective agar plates using Kanamycin (50 μg/ml) and Gentamicin (50 μg/ml) selection and incubated at 28° C. for 2 days to obtain single colonies.
A single colony was used to inoculate 5 ml LB plus Kanamycin and Gentamicin at 28° C. until dense culture was obtained to make a glycerol stock with 15% glycerol. Aliquots of the glycerol stocks were kept at −80° C.
Transformation of Opium Poppy with pRI 201-An_35 S::RBA
Wild type line HN4 and homozygous rba-1 deletion mutant material were used for stable transformation and expression of RBA using the pRI 201-An_35 S::RBA expression construct.
Prior to sowing poppy seeds were sterilised using a chlorine vapour method. A glass petri dish containing the seed was placed lid-off in a sealable box in the fume hood. The petri dish lid was also placed in the box along with a beaker containing a 3% hydrochloric acid solution. To release chlorine vapour, a 2.5 g Presept® Disinfectant Tablet (Johnson & Johnson) was added to the 3% hydrochloric acid solution and the box was sealed with parafilm and placed in a fume hood. After 90 minutes the seal was broken and the lid removed allowing the petri dish containing the sterilised seed to be carefully covered by the glass lid and removed from the sealable box to the laminar flow hood.
The sterilised seed was sown onto B50 agar medium in sterile tissue culture pots. B50 media composition is—B5 micro and macro elements (Duchefa Biochemie), B5 vitamins (Duchefa Biochemie), 20% sucrose, 2-(N-morpholino)ethanesulfonic acid (MES) buffer and 0.8% plant cell culture agar. The sown seed was stratified at 4° C. for 48 hours before being transferred to a growth cabinet to germinate and grow at 24° C. for 7-9 days.
A. tumefaciens pre-cultures were set up using glycerol stocks prepared from A. tumefaciens GV3101 cells containing the pRI 201-AN_35 S:RBA construct. A single 50 μl aliquot was added to 5 ml LB containing Kanamycin (50 μg/ml) and Gentamicin (50 μg/ml). The cultures were grown in a shaking incubator at 28° C. overnight. About 100 μl of the pre-culture was used to inoculate a 50 ml main culture of LB medium containing Kanamycin (50 μg/ml) and Gentamicin (50 μg/ml) which was grown in a shaking incubator at 28° C. overnight. Once an OD600 of 0.3-0.8 was reached cells were harvested by centrifugation at 4000 rpm at 4° C. for 20 minutes. Inoculation culture was prepared by resuspending the cell pellet in inoculation medium to an OD600 of between 0.2-0.3 and incubated for 1-3 hours at 28° C. on a shaker. The inoculation medium consisted of a liquid B50 medium without MES, but including Acetosyringone (100 μM).
Hypocotyl explants were carefully dissected from the roots and cotyledons of 7-9 day old poppy seedlings and immediately transferred into petri dishes containing the Agrobacterium inoculation culture. After gently shaking on an orbital shaker at room temperature for 15 minutes the explants were blotted on sterile filter paper to remove excess culture medium and transferred to solid co-cultivation medium for 2 to 3 days at 24° C. Co-cultivation medium consisted of the B50 medium plus agar with the addition of 2,4-Dichlorophenoxyacetic acid (1 μg/ml) and Acetosyringone (100 μM).
After the co-cultivation period the explants were immersed in Timentin®-containing solution (150 μg/mL) in 50 mL Falcon tubes, agitated for 60 seconds and washed three times in sterile water. After blotting on sterile filter paper, explants were transferred to callus induction (CM) plates. CM medium consisted of B50 medium plus agar with Timentin® (150 μg/mL) and Paromomycin (30 μg/mL) in addition to 2,4-Dichlorophenoxyacetic acid (1 μg/ml). The CM plates were incubated in a growth cabinet at 24° C. Explants were transferred to fresh CM plates of the same composition at three-weekly intervals.
Two types of calli formed on the explants (Chitty et al., 2003): Type I and Type II. Type I is a colourless loose callus that becomes brown over time. Type II forms as small regions of white compact embryogenic callus that appear after about 12 weeks on transformed explants.
Type II callus was removed from the explants and transferred to B50 media containing Timentin® (150 μg/mL) and Paromomycin (30 μg/mL) and incubated in a growth cabinet at 18-20° C. Type II calli were transferred to fresh B50 plates of the same composition at three-weekly intervals until somatic embryos started to form (typically after 2-3 culture periods).
Development of Somatic Embryos into Plantlets
Once plantlets were starting to develop from the somatic embryos, they were transferred to B50 medium with Timentin® (150 μg/mL) and Paromomycin (30 μg/mL) selection in baby jars to allow them to grow taller and to develop roots. Any basal callus as well as brown or dead leaves or shoots were removed before transfer.
Plantlets with fully established root system were transferred to soil and initially kept in a humid environment and slowly hardened off by decreasing humidity over time. Once acclimatised to normal humidity levels the plants were transferred to plant growth cabinets or the glasshouse.
Exactly the same procedure was followed to transform HN4 wild type and homozygous rba-1 deletion mutant material with pRI 201-AN empty vector, which is identical to the vector pRI 201-AN_35 S::RBA shown in
A forward screen was carried out on the M2 generation of fast-neutron-mutagenized material of the HN4 variety. The screen was carried out in two rounds whereby in the first round latex of 12 week-old plants from the entire M2 mutant population was analysed by high throughput screening using Direct Injection MS without prior separation on the UPLC to monitor relative proportions of non-isobaric alkaloid groups. Any material showing latex phenotypes substantially different to HN4 control material where selected for the second round of screening by UPLC-MS against reference standards to obtain a full benzylisoquinoline alkaloid profile of the mature dry capsules.
Two independent M2 mutant plants belonging to different M2 seed families were identified in the latex screen as lacking any BIAs. The complete absence of measurable quantities of BIAs in these plants was confirmed by analysing mature dry M2 capsule material (Table 1). Following the characterisation of the molecular basis for the loss of BIAs in the mutant plants (Example 3 and 4) the mutant alleles were named ‘regulator of benzylisoquinoline alkaloids-1’ and ‘-2’ (rba-1 and rba-2), respectively.
The rba-1 and rba-2 M2 mutant plants were self-pollinated and their respective progenies further propagated through several generations by self-pollinations to obtain a larger amount of material for testing the heritability and stability of the ‘loss of all BIAs’ phenotype.
As was the case for the M2 capsules, dry capsule material of the M3, M4 and M5 generations were devoid of the major alkaloids, morphine and noscapine, found in substantial amounts in the HN4 parental line. They also did not contain any other alkaloids in measurable quantities (Table 1). The analysis thus confirmed the ‘loss of all BIAs’ phenotype and showed that it is stably inherited.
The rba-1 mutant was crossed with two different morphine varieties, HM1 and HM8. Capsules of the F1 progeny, in which the rba-1 mutant allele is in the heterozygous state, did contain morphine and other morphinan alkaloids in similar amounts as the parental morphine varieties demonstrating that the rba-1 mutant allele is recessive (Table 1). Given that the rba-1 and rba-2 mutants carry different mutant alleles of the same genetic locus (Examples 3 and 4) the rba-2 allele can also be considered recessive.
A cross was carried out between the rba-1 and rba-2 mutants to test if they would complement each other with respect to BIA production. Since the HN4 wildtype phenotype was not rescued in capsules of the resulting F1 generation (Table 1) the rba-1 and rba-2 mutant alleles must belong to same complementation group with the respective mutations affecting same genetic locus. This was later confirmed by the characterisation of the molecular basis of the rba-1 and rba-2 mutants (Example 3).
To further characterise the mutants with respect to gene expressions levels of BIA biosynthesis genes, RNA sequencing was carried out on stem material of the M3 generation of the rba-1 and rba-2 as well as of the parental variety HN4. RPKM-normalised expression analysis revealed that compared to non-mutagenised HN4 wild type material the vast majority of BIA biosynthesis genes are strongly down-regulated in rba-1 and rba-2 mutant material compared to HN4 controls (Table 2). One example is NORCOCLAURINE SYNTHASE (NCS), which catalyses the first committed step in the biosynthesis of all BIAs in opium poppy (Samanani et al. (2004) Plant J. 40: 302-313). There are two NCS genes annotated in the opium poppy reference genome: PSUN64910.1 on scaffold ups_21 and PSUN03620.1 on scaffold ups_10. The RPKM-values of the rba-1 and rba-2 mutants are just 2 compared to 300 in the HN4 parental variety. Likewise, the respective RPKM-values of PSUN03620.1 are 8 and 2 in the rba-1 and -2 mutants compared to 130 in HN4.
Determination of the Molecular Basis for of the ‘Loss of all BIAs’ Phenotype in the rba-1 and rba-2 Mutants
To determine the molecular basis for the ‘loss of all BIAs’ phenotype, the genomes of rba-1 and rba-2 mutants were re-sequenced using the 10× Genomics platform and assembled inhouse using the Supernova software. The respective deletions in each mutant causing the ‘loss of all BIAs’ phenotype were identified using bioinformatic approaches as described in materials and methods.
With respect to the HN1 genome reference assembly 989 kb of genomic sequence on chromosome 9 containing 12 genes are deleted in the rba-1 mutant (
In the rba-2 mutant, 780 kb of genomic sequence on chromosome 9 containing 11 genes were found to be deleted (
The fact the rba-1 and rba-2 deletion alleles identified in the respective mutants are unique in terms of their overall size and precise boundaries is consistent with them having arisen through independent mutation events in different M1 plants.
Consistent with the genetic complementation analysis (Example 1), the rba-1 and rba-2 deletion alleles overlap for 506 kb of genomic sequence (
The deletions identified in the rba-1 and rba-2 mutants by means of genome re-sequencing and assembly were further corroborated by gene expression analysis using the RNA sequencing data set from upper stems described above. The RPKM values of those genes residing in the respective deleted regions on chromosome 9 were compared between the rba-1 and rba-2 mutants and HN4 wild type (Table 3).
RPKM values of deleted genes would be expected to be low in the respective mutants compared to HN4 wild type—provided they expressed to reasonable levels in stem tissue of the HN4 wild type in the first place (RPKM values of deleted genes may not be absolutely zero in the respective mutants: low RPKM values may reflect background noise inherent to the methodology such as non-specific mapping of short reads from other genes with some degree of sequence similarity rather than genuinely low expression).
The RPKM values of the genes residing in the overlapping region deleted in both mutants are generally zero or low in both mutants compared to HN4 wild type. In contrast, RPKM-values of genes outside the overlapping region (and expressed to reasonable levels in HN4) show low expression only in the respective mutant in which they are deleted but not the other.
The RNA sequencing gene expression analysis thus confirms the respective deletions identified in the rba-1 and rba-2 genome assemblies as bona fide deletions.
Of the six genes residing in the overlapping region deleted in both mutants only one shows homology with genes known to be involved in transcriptional regulation: PS0917990.1 encodes a bHLH transcription factor.
Since its deletion results in loss of expression of most genes known to be involved in BIA biosynthesis and, consequently, the ‘loss of all BIAs’ phenotype PS0917990.1 is a master regulator of BIA synthesis and accordingly named REGULATOR OF BENZYLISOQUINOLINE ALKALOIDS (RBA).
Cloning and Sanger-sequencing from genomic and cDNA confirmed the RBA sequence from HN4 to be identical to that from the HN1 reference genome assembly.
Number | Date | Country | Kind |
---|---|---|---|
2107438.0 | May 2021 | GB | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/IB2022/054862 | 5/24/2022 | WO |