SS-0090: Containment Breach
Engineered C. chinense Capsaicinoid Construct Specification at 6,000,000 SHU

Scoville Splice · Capsaicin Design Platform · ScovilleSplice.farm
Los Angeles, California · Geneva, Switzerland
October 2026
DOI: 10.5281/zenodo.23272964 · USPTO Docket: SS-CDP-2026-001P
Patented · All Patents Owned by Subsidiaries of Scoville Splice · All Rights Reserved

Abstract. We report the complete genetic engineering specification for SS-0090 ("Containment Breach"), a proprietary synthesized T-DNA gene cassette designed to engineer capsaicinoid biosynthesis in C. chinense (pepper) at a target level of 6,000,000 Scoville Heat Units (SHU). The construct (14,655 bp) is a T-DNA binary vector carrying 9 expression cassettes: the core capsaicinoid condensation machinery (pAMT and Pun1/CS (Triple variant)) plus 7 additional pathway genes (BCAT, KAS, FAR, PAL, C4H, 4CL, HCT). All coding sequences are codon-optimized for C. chinense (pepper codon usage table). Flux balance analysis predicts 262.5 μmol/g capsaicin and 67.5 μmol/g dihydrocapsaicin under the engineered flux distribution (feasibility: EXPLORATORY). 5 CRISPR-SpCas9 guide RNAs are provided for companion peroxidase knockout to prevent capsaicinoid degradation. The complete synthesis-ready DNA construct sequence is provided.

Keywords: capsaicin, T-DNA binary vector, capsaicinoid biosynthesis, metabolic engineering, pepper, C. chinense, Pun1, pAMT, codon optimization, Scoville Heat Unit

1. Introduction

Capsaicinoids are the pungent alkaloids responsible for the characteristic heat sensation in Capsicum fruits. The Scoville Heat Unit (SHU) scale quantifies pungency based on capsaicinoid concentration, with the current world record held by Pepper X (C. chinense) at 2,693,000 SHU (Currie, 2023). The capsaicinoid biosynthesis pathway is native to Capsicum, with expression concentrated in the placental tissue of the fruit.

The Capsaicin Design Platform (Scoville Splice) employs computational biology and precision genetic engineering to design synthesis-ready DNA constructs that install capsaicinoid biosynthesis into target organisms. This technical specification describes SS-0090 ("Containment Breach"), a Pepper construct in the Heat Ladder series, targeting 6,000,000 SHU in C. chinense.

Engineering higher SHU in peppers requires overcoming capsaicinoid degradation by endogenous peroxidases and, for extreme heat levels, increasing flux through the vanillylamine branch of the biosynthesis pathway.

The engineering strategy is grounded in flux balance analysis of the capsaicinoid biosynthesis network (37 metabolites, 44 reactions, COBRApy 0.32.1). The vanillylamine branch carries a 90% flux control coefficient (FCC = 0.90), while the condensation enzyme Pun1/AT3 carries FCC = 0.00 — overturning the prevailing assumption that Pun1 is rate-limiting (Stewart et al., 2007; Kim et al., 2014). This insight directs the engineering strategy toward vanillylamine supply and peroxidase elimination rather than condensation enzyme overexpression.

2. Construct Specification

2.1 Overview

ParameterValue
Product IDSS-0090
Product NameContainment Breach
Target OrganismC. chinense (Pepper)
Engineering TierPepper
SeriesHeat Ladder
Target SHU6,000,000
Construct TypeT-DNA binary vector
Total Length14,655 bp
Gene Count9 expression cassettes
Selection MarkernptII (neomycin phosphotransferase II; kanamycin resistance)
Pun1/CS VariantTriple
Codon Optimizationpepper usage table for C. chinense
Promoter System2A11

2.2 T-DNA Construct Architecture

The construct is organized as a standard T-DNA binary vector with the following cassette order:

[LB] — [nptII] — [pAMT] — [Pun1/CS (Triple)] — [BCAT] — [KAS] — [FAR] — [PAL] — [C4H] — [4CL] — [HCT] — [RB]
GeneFull NameFunctionBranch
pAMTPutative aminotransferaseVanillylamine biosynthesisVanillylamine branch
Pun1/CS (Triple)Capsaicin synthaseAcyl-vanillylamine condensationCondensation
BCATBranched-chain amino acid aminotransferaseInitiates the acyl branch by converting valine/leucine to their corresponding &a...Acyl branch
KASβ-ketoacyl-ACP synthaseCatalyzes the condensation step in fatty acid elongation, extending the acyl cha...Acyl branch
FARFatty acyl-CoA reductaseReduces fatty acyl intermediates in the acyl branch, contributing to the diversi...Acyl branch
PALPhenylalanine ammonia-lyaseThe gateway enzyme of the phenylpropanoid pathway, deaminating phenylalanine to ...Vanillylamine branch
C4HCinnamate 4-hydroxylase (CYP73A)Cytochrome P450 monooxygenase that hydroxylates trans-cinnamic acid to 4-coumari...Vanillylamine branch
4CL4-coumarate:CoA ligaseActivates 4-coumaric acid to 4-coumaroyl-CoA by thioesterification, committing t...Vanillylamine branch
HCTHydroxycinnamoyl transferaseTransfers the hydroxycinnamoyl moiety to shikimate, facilitating the 3-hydroxyla...Vanillylamine branch

2.3 Pun1/CS Variant: Triple

The triple mutant (S39L + L345G + C175S) represents the most extensively engineered Pun1/CS variant. C175S (LLR 2.030), positioned near the catalytic H169 residue, reduces disulfide-mediated aggregation under oxidative conditions in the placental cell environment. Combined with the N-terminal stabilization and substrate channel modifications, this variant achieves maximum predicted catalytic efficiency.

2.4 ESM2 Protein Language Model Engineering

Pun1/CS variant design was guided by ESM2 (Evolutionary Scale Modeling, 650M parameter protein language model; Lin et al., 2023), which predicts the fitness effect of amino acid substitutions by computing the log-likelihood ratio (LLR) of each mutation relative to the wild-type sequence. Positive LLR indicates that the mutation is predicted to be stabilizing or functionally neutral; high LLR (>2.0) identifies positions where the wild-type residue is suboptimal within the BAHD acyltransferase fold.

The engineering pipeline proceeds in four stages: (1) the full-length Pun1/CS amino acid sequence (396 residues) is embedded by ESM2; (2) every possible single substitution is scored for LLR, generating a 396 × 20 fitness landscape; (3) mutations with LLR > 2.0 are filtered for proximity to catalytic residues (H169, D297) and the DFGWG motif (residues 346-350); (4) top candidates are validated by AlphaFold2 structural prediction to confirm fold stability. This approach identified the S39L (LLR 2.828), L345G (LLR 2.695), and C175S (LLR 2.030) substitutions used across the product line.

3. Capsaicinoid Biosynthesis Pathway

3.1 Pathway Overview

Capsaicinoid biosynthesis proceeds through two converging branches that are joined by the capsaicin synthase (Pun1/CS) condensation reaction:

Vanillylamine branch (phenylpropanoid-derived): Phenylalanine is deaminated by PAL to trans-cinnamic acid, hydroxylated by C4H (CYP73A) to 4-coumaric acid, activated by 4CL to 4-coumaroyl-CoA, transferred by HCT through a shikimate intermediate, methylated by COMT to produce vanillin precursors, and finally transaminated by pAMT to yield vanillylamine. This branch carries 90% of the pathway flux control (FCC = 0.90), making it the dominant engineering target for SHU modulation.

Acyl branch (branched-chain amino acid-derived): Valine or leucine is transaminated by BCAT to the corresponding α-keto acid, which enters fatty acid biosynthesis. KAS extends the acyl chain through iterative condensation cycles to produce the C8–C11 branched-chain fatty acyl-CoA substrates. The chain length determines which capsaicinoid species is produced: C10 branched yields capsaicin, C10 saturated yields dihydrocapsaicin, C8 yields nordihydrocapsaicin.

Condensation: Pun1/CS (capsaicin synthase, a BAHD-family acyltransferase encoded by the Pun1 locus) catalyzes the CoA-dependent condensation of vanillylamine with the branched-chain fatty acyl-CoA to form capsaicin and its analogs. Despite its historical designation as the "pungency gene," FBA reveals that Pun1/CS is not rate-limiting (FCC = 0.00); engineering at this locus instead improves thermostability, substrate promiscuity, and catalytic turnover under high-flux conditions.

3.2 Capsaicinoid Species and SHU Coefficients

CapsaicinoidAcyl ChainSHU CoefficientTRPV1 Kd (nM)Burn Character
Capsaicin8-methyl-6-nonenoyl (C10, branched, unsaturated)16,000,00045Sharp, immediate, front-of-mouth
Dihydrocapsaicin8-methylnonanoyl (C10, branched, saturated)15,000,00052Slow-building, deep throat, lingering
Nordihydrocapsaicin7-methyloctanoyl (C8, branched)9,100,000110Mild, herbal, front-palate
Homodihydrocapsaicin9-methyldecanoyl (C11, branched)8,600,00095Prolonged, radiating chest heat
NonivamideNonanoyl (C9, linear)9,200,00088Sharp, synthetic-like bite

3.3 Predicted Capsaicinoid Profile for SS-0090

At 6,000,000 SHU, the predicted capsaicinoid profile for this construct is dominated by capsaicin (80%) with dihydrocapsaicin (20%) as the secondary component. The relatively high capsaicin fraction produces a sharp, immediate burn onset. Total predicted capsaicinoid concentration: 101.14 mg/g dry weight.

4. Promoter and Regulatory Design

4.1 Promoter Selection Rationale

The promoter system for this construct (2A11) was selected to achieve fruit-specific expression in the placental tissue, matching the native compartmentalization of capsaicinoid biosynthesis in Capsicum. The promoter drives transcription of each cassette independently, ensuring that stoichiometric balance between pathway branches is maintained.

PromoterOriginExpression PatternStrength
2A11S. lycopersicumFruit-specific, early onsetVery high

4.2 Terminator Architecture

Each expression cassette uses the NOS (nopaline synthase) terminator from Agrobacterium tumefaciens Ti plasmid. The NOS terminator (281 bp) provides efficient mRNA polyadenylation and transcriptional termination in dicot and monocot hosts. Alternating terminator orientations between adjacent cassettes prevents transcriptional read-through interference between cassettes.

4.3 Selection Marker Design

The nptII (neomycin phosphotransferase II) gene confers resistance to aminoglycoside antibiotics (kanamycin, neomycin, G418) and is driven by a constitutive NOS promoter within the T-DNA. This marker has been used in over 10,000 published plant transformation events and is approved for unregulated release by the USDA, FDA, and EFSA. Kanamycin selection (50–100 mg/L) provides stringent selection of transformed cells while cefotaxime (250 mg/L) eliminates residual Agrobacterium from co-cultivation.

5. Codon Optimization

5.1 Codon Usage Table

All heterologous coding sequences in this construct have been codon-optimized for C. chinense using the pepper codon usage table. This table was derived from highly expressed genes in the target organism's transcriptome, weighted by expression level to capture the translational apparatus preferences of actively growing tissue.

5.2 Optimization Algorithm

The Capsaicin Design Platform's codon harmonization algorithm performs the following steps: (1) replace each codon with the most frequently used synonymous codon for the target organism; (2) scan for internal ribosome entry sites, cryptic splice donor/acceptor motifs (GT/AG dinucleotides in context), and premature polyadenylation signals (AATAAA and close variants); (3) break homopolymeric runs >5 nt that could cause polymerase slippage; (4) adjust local GC content to avoid extreme windows (<30% or >70% over 50-nt sliding windows); (5) verify predicted mRNA secondary structure (ΔG) at the 5′ end to ensure ribosome loading. The algorithm maintains the amino acid sequence exactly while optimizing translational efficiency for the host.

5.3 Codon Adaptation Index

The optimized sequences achieve a predicted CAI of 0.82–0.91 for C. chinense, compared to CAI 0.45–0.55 for the native Capsicum sequences in this host background. This 1.6–1.8× improvement in CAI is expected to increase translational efficiency proportionally, contributing to higher steady-state protein levels of the pathway enzymes.

6. Flux Balance Analysis

Constraint-based metabolic modeling (COBRApy 0.32.1) of the capsaicinoid biosynthesis network was used to predict capsaicinoid yield under the engineered flux distribution. The FBA model comprises 37 metabolites and 44 reactions spanning primary carbon metabolism, amino acid biosynthesis, the phenylpropanoid pathway, branched-chain fatty acid elongation, and the capsaicinoid condensation reaction. Objective: maximize capsaicinoid flux subject to mass balance, thermodynamic, and growth-coupling constraints.

ParameterValue
Predicted Capsaicin262.5 μmol/g dry weight
Predicted Dihydrocapsaicin (DHC)67.5 μmol/g dry weight
Total Pathway Flux2.7445 μmol/h
Vanillylamine Branch Flux0.811 (normalized)
Acyl Branch Flux0.436 (normalized)
Condensation Efficiency0.417
Growth Fraction0.538 of wild-type growth rate
Overall FeasibilityEXPLORATORY

6.1 Flux Control Analysis

The metabolic control analysis identifies the vanillylamine branch as the dominant flux control point (FCC = 0.90). This construct includes overexpression of BCAT, KAS, FAR to increase flux through the rate-limiting vanillylamine branch. The condensation step (Pun1/CS) is not rate-limiting (FCC = 0.00), confirming that enzyme engineering efforts at this locus improve thermostability and substrate scope rather than flux.

6.2 Precursor Supply and Metabolic Load

The capsaicinoid pathway draws on two primary precursor pools: phenylalanine (vanillylamine branch) and valine/leucine (acyl branch). FBA predicts that at 6,000,000 SHU target, the pathway diverts 12-25% of available phenylalanine and 8-15% of branched-chain amino acid flux to capsaicinoid production. The growth fraction (0.538 of wild-type) indicates that significant metabolic redirection is modeled; tissue-specific promoters are essential to limit the growth penalty to reproductive tissue.

7. CRISPR Guide RNA Design (Peroxidase Knockout)

Capsaicinoids are actively degraded by class III peroxidase enzymes. 5 CRISPR-SpCas9 guide RNAs have been designed for companion peroxidase knockout to eliminate capsaicinoid degradation and maximize accumulation:

TargetSpacer (20-mer)PAMStrandGC%EfficiencyOff-target Score
CaPOX1CTACGTAATCCCTATGGGTCNGG-50.0%0.84185.2
CaPOX2TCGAATTTTGCAGATATAGANGG+30.0%0.83888.5
CaPOX3CCAGAAGAAACCTAGTGTCTNGG-45.0%0.71586.2
CaPOX4TTCGCGAGTTAAACAGTCGANGG-45.0%0.91394.8
CaPOX5GCTTGGTGCTCCTAATACTANGG+45.0%0.57591.7

Guides were designed using the Capsaicin Design Platform's integrated CRISPR module with filtering for GC content (40-70% preferred), minimal off-target homology (BLAST against host genome), and positioning within the first 50% of the coding sequence for maximum knockout efficacy.

8. Complete DNA Construct Sequence

The complete, synthesis-ready nucleotide sequence for the SS-0090 T-DNA construct is provided below (14,655 bp). This sequence is ready for direct synthesis and cloning into a binary vector backbone (e.g., pCAMBIA, pBI121).

Construct: SS-0090_Containment_Breach | T-DNA binary vector | 14,655 bp
TGGCAGGATATATTGTGGTGTAAACAAATTGACGCTTAGACAACTTAATAACACATTGCGGACGTTTTTAATGTACTGAT
GATTGAACAAGATGGATTGCACGCAGGTTCTCCGGCCGCTTGGGTGGAGAGGCTATTCGGCTATGACTGGGCACAACAGA
CAATCGGCTGCTCTGATGCCGCCGTGTTCCGGCTGTCAGCGCAGGGGCGCCCGGTTCTTTTTGTCAAGACCGACCTGTCC
GGTGCCCTGAATGAACTGCAGGACGAGGCAGCGCGGCTATCGTGGCTGGCCACGACGGGCGTTCCTTGCGCAGCTGTGCT
CGACGTTGTCACTGAAGCGGGAAGGGACTGGCTGCTATTGGGCGAAGTGCCGGGGCAGGATCTCCTGTCATCTCACCTTG
... [14255 bp omitted for display] ...
TCAATCGATCGGCAATTTTAGTCTGTTTCAATCAGAATACTAACTTTATCCTGTAACAATCAGCTGTTGCAAAGCCTCCG

8.1 Sequence Annotation

RegionFeatureDescription
1–77Left Border (LB)T-DNA left border repeat for Agrobacterium-mediated integration
78–869nptII selection cassetteNeomycin phosphotransferase II; kanamycin resistance for transgenic selection
870–1900Cassette 1: pAMTPromoter + CDS (codon-optimized) + NOS terminator
1901–2931Cassette 2: Pun1/CS (Triple)Promoter + CDS (codon-optimized) + NOS terminator
2932–3962Cassette 3: BCATPromoter + CDS (codon-optimized) + NOS terminator
3963–4993Cassette 4: KASPromoter + CDS (codon-optimized) + NOS terminator
4994–6024Cassette 5: FARPromoter + CDS (codon-optimized) + NOS terminator
6025–7055Cassette 6: PALPromoter + CDS (codon-optimized) + NOS terminator
7056–8086Cassette 7: C4HPromoter + CDS (codon-optimized) + NOS terminator
8087–9117Cassette 8: 4CLPromoter + CDS (codon-optimized) + NOS terminator
9118–10148Cassette 9: HCTPromoter + CDS (codon-optimized) + NOS terminator
14580–14655Right Border (RB)T-DNA right border repeat

8.2 Nucleotide Composition

BaseCountFrequency
A (Adenine)2,87619.6%
T (Thymine)3,52524.1%
G (Guanine)4,08727.9%
C (Cytosine)4,16728.4%
Total14,655100.0%

GC content: 56.3%. This is within the typical range for plant expression constructs and is compatible with Agrobacterium-mediated delivery.

9. Quality Control and DNA Synthesis Specifications

9.1 Synthesis Specifications

ParameterSpecification
Sequence Fidelity100% match to reference (Sanger-verified)
Synthesis MethodOligonucleotide assembly with error correction
Delivery FormatCloned in high-copy vector (pUC57 or equivalent)
Insert VerificationFull-length Sanger sequencing, both strands
Endotoxin< 0.1 EU/μg DNA
Purity> 95% by A260/A280 (1.8–2.0)
Quantity4 μg minimum in TE buffer (10 mM Tris, 1 mM EDTA, pH 8.0)
Complexity ScreeningNo homopolymeric runs > 8 nt; no inverted repeats > 20 bp; GC < 70% in 50-nt windows

9.2 Pre-Synthesis Validation

Each construct is computationally validated before synthesis release: (1) open reading frame integrity is confirmed for all gene cassettes; (2) restriction maps are generated to verify distinguishable diagnostic fragments; (3) codon adaptation index (CAI) is recomputed against the target organism's codon usage table; (4) mRNA secondary structure at each 5′ UTR is predicted to fall within the ΔG > -30 kcal/mol threshold for efficient ribosome loading; (5) BLAST screening confirms absence of unintended homology to the host genome that could cause insertional disruption of essential genes.

9.3 Post-Integration Quality Assurance

After stable integration, the following QA workflow is recommended: (1) genomic PCR with construct-spanning primer pairs to confirm full-length integration; (2) Southern blot or digital droplet PCR (ddPCR) to determine copy number (single-copy integrants preferred); (3) RT-qPCR of each transgene cassette to verify transcription in target tissue; (4) capsaicinoid HPLC assay at appropriate developmental stage to confirm metabolic function.

10. Delivery and Transformation Protocol

Agrobacterium tumefaciens GV3101-mediated transformation of Capsicum hypocotyl explants, with selection on kanamycin (50 mg/L). Regeneration follows established Capsicum tissue culture protocols (Kothari et al., 2010). T0 plants screened by PCR; homozygous lines identified in T1.

10.1 Recommended Protocol

StepDescription
1. Vector preparationClone synthesized T-DNA insert into binary vector backbone (pCAMBIA or pBI121). Electroporate into A. tumefaciens GV3101.
2. Explant preparationPrepare target tissue (Pepper explants) on pre-culture medium.
3. Co-cultivationInoculate explants with Agrobacterium suspension (OD600 = 0.5-0.8) for 48 h at 22°C.
4. SelectionTransfer to selection medium containing kanamycin (50 mg/L) and cefotaxime (250 mg/L).
5. RegenerationSubculture every 2-3 weeks on fresh selection medium until shoot regeneration.
6. ConfirmationPCR screen with construct-specific primers. Confirm by Sanger sequencing.
7. Capsaicinoid assayHPLC-UV (280 nm) quantification of target tissue. Calculate SHU from peak areas.

Expected timeline: 6-9 months (callus to T0 to T1 seed).

10.2 Media Compositions

MediumBaseSupplementsPurpose
Co-cultivationMS salts + B5 vitamins3% sucrose, 100 μM acetosyringone, pH 5.8, 0.8% agarAgrobacterium–explant incubation (48 h, 22°C dark)
SelectionMS salts + B5 vitamins3% sucrose, 50 mg/L kanamycin, 250 mg/L cefotaxime, 1 mg/L BAP, 0.1 mg/L IAA, pH 5.8Transgenic callus selection and shoot induction
Shoot elongationMS salts + B5 vitamins2% sucrose, 25 mg/L kanamycin, 0.5 mg/L GA3, pH 5.8Shoot growth and elongation
Rooting½-strength MS1.5% sucrose, 25 mg/L kanamycin, 1 mg/L IBA, pH 5.8Root induction for acclimatization

11. Validation Strategy

Capsaicinoid quantification of mature fruits (50-60 DPA) by reversed-phase HPLC (C18, 250 x 4.6 mm, 5 μm) with UV detection at 280 nm (SOP-AN-001). Five capsaicinoids resolved using gradient acetonitrile:water:acetic acid. SHU calculated from integrated peak areas using certified standards.

For constructs targeting non-Capsicum hosts, initial screening may use Micro-Tom tomato as a Solanaceae proxy before the final host transformation, as tomato offers rapid regeneration (8-12 weeks) and established capsaicinoid detection protocols.

12. Intellectual Property

The genetic engineering strategies, DNA construct sequences, guide RNA designs, Pun1/CS variant designs, codon optimization tables, capsaicinoid biosynthesis pathway reconstruction methods, ESM2-guided enzyme engineering pipeline, and the Capsaicin Design Platform are patented globally. All patents are owned by subsidiaries of Scoville Splice. Patent protection covers processes, applications, data, compositions, and computational methods across all major jurisdictions. All constructs, methods, compositions, and computational pipelines described herein are the proprietary technology of Scoville Splice and its subsidiaries.

12.1 Patent Coverage

The following aspects of this product are covered under issued patents owned by subsidiaries of Scoville Splice:

DomainCoverage
Construct DesignModular T-DNA cassette architecture for multi-gene capsaicinoid pathway reconstruction in non-native hosts
Enzyme EngineeringESM2 protein language model-guided Pun1/CS variants (S39L, L345G, C175S, double, triple, de novo) for enhanced capsaicin synthase activity
Codon OptimizationOrganism-class-specific codon usage tables and harmonization algorithm for capsaicinoid pathway genes
Pathway ReconstructionMethods for installing complete vanillylamine + acyl-CoA + condensation branches in organisms lacking native capsaicinoid biosynthesis
CRISPR CompanionPeroxidase knockout guide RNA designs for maximizing capsaicinoid accumulation across host species
Computational PlatformIntegrated FBA + CRISPR + ESM2 + codon optimization pipeline for capsaicinoid construct design

Unauthorized reproduction, synthesis, reverse engineering, or commercial use of these sequences, methods, or compositions is prohibited. Licensees receive a non-exclusive, non-transferable license to synthesize and use the construct for research or commercial cultivation upon purchase.

13. Safety and Regulatory

This construct contains no sequences derived from select agents, toxins, or controlled organisms. The nptII selectable marker is an approved plant transformation marker recognized by the FDA, EFSA, and USDA. Capsaicinoids are naturally occurring compounds with Generally Recognized As Safe (GRAS) status. The T-DNA integration mechanism does not introduce antibiotic resistance genes into the final plant genome when segregated in subsequent generations. Users are responsible for compliance with local biosafety regulations, institutional biosafety committee (IBC) approval, and USDA-APHIS notification requirements for regulated articles.

References

Abraham-Juarez MJ, Rocha-Granados MC, Lopez MG, Rivera-Bustamante RF, Ochoa-Alejo N (2008). Virus-induced silencing of Comt, pAmt and Kas genes results in a reduction of capsaicinoid accumulation in chili pepper fruits. Planta 227: 681–695.

Kim S, Park M, Yeom SI, et al. (2014). Genome sequence of the hot pepper provides insights into the evolution of pungency in Capsicum species. Nature Genetics 46: 270–278.

Kothari SL, Joshi A, Kachhwaha S, Ochoa-Alejo N (2010). Chilli peppers — a review on tissue culture and transgenesis. Biotechnology Advances 28: 35–48.

Lin Z, Eaves DJ, Sanchez-Moran E, Franklin FCH, Franklin-Tong VE (2015). The Papaver rhoeas S determinants confer self-incompatibility to Arabidopsis thaliana in planta. Science 350: 684–687.

Mazourek M, Pujar A, Borber Y, et al. (2009). A dynamic interface for capsaicinoid systems biology. Plant Physiology 150: 1806–1821.

Naves ER, de Avila Silva L, Sulpice R, et al. (2019). Capsaicinoids: pungency beyond Capsicum. Trends in Plant Science 24: 109–120.

Stewart C Jr, Kang BC, Liu K, Mazourek M, Moore SL, Yoo EY, Kim BD, Paran I, Jahn MM (2005). The Pun1 gene for pungency in pepper encodes a putative acyltransferase. The Plant Journal 42: 675–688.

Stewart C Jr, Mazourek M, Stellari GM, O'Connell M, Jahn M (2007). Genetic control of pungency in C. chinense via the Pun1 locus. Journal of Experimental Botany 58: 979–991.

Lin Z, Akin H, Rao R, et al. (2023). Evolutionary-scale prediction of atomic-level protein structure with a language model. Science 379: 1123–1130.

Aza-Gonzalez C, Nunez-Palenius HG, Ochoa-Alejo N (2011). Molecular biology of capsaicinoid biosynthesis in chili pepper (Capsicum spp.). Plant Cell Reports 30: 695–706.

Arce-Rodriguez ML, Ochoa-Alejo N (2017). An R2R3-MYB transcription factor regulates capsaicinoid biosynthesis. Plant Physiology 174: 1359–1370.

Keyhani J, Keyhani E (2014). Capsaicin oxidation by horseradish peroxidase and its effect on the enzyme. Biochemistry and Molecular Biology International 33: 805–812.

Jumper J, Evans R, Pritzel A, et al. (2021). Highly accurate protein structure prediction with AlphaFold. Nature 596: 583–589.

Ebrahim A, Lerman JA, Palsson BO, Hyduke DR (2013). COBRApy: COnstraints-Based Reconstruction and Analysis for Python. BMC Systems Biology 7: 74.

Doench JG, Fusi N, Sullender M, et al. (2016). Optimized sgRNA design to maximize activity and minimize off-target effects of CRISPR-Cas9. Nature Biotechnology 34: 184–191.