Experimental and Field Data Support Range Expansion in an Allopolyploid Arabidopsis Owing to Parental Legacy of Heavy Metal Hyperaccumulation.

Empirical evidence is limited on whether allopolyploid species combine or merge parental adaptations to broaden habitats. The allopolyploid Arabidopsis kamchatica is a hybrid of the two diploid parents Arabidopsis halleri and Arabidopsis lyrata. A. halleri is a facultative heavy metal hyperaccumulator, and may be found in cadmium (Cd) and zinc (Zn) contaminated environments, as well as non-contaminated environments. A. lyrata is considered non-tolerant to these metals, but can be found in serpentine habitats. Therefore, the parents have adaptation to different environments. Here, we measured heavy metals in soils from native populations of A. kamchatica. We found that soil Zn concentration of nearly half of the sampled 40 sites was higher than the critical toxicity level. Many of the sites were near human construction, suggesting adaptation of A. kamchatica to artificially contaminated soils. Over half of the A. kamchatica populations had >1,000 µg g-1 Zn in leaf tissues. Using hydroponic treatments, most genotypes accumulated >3,000 µg g-1 Zn, with high variability among them, indicating substantial genetic variation in heavy metal accumulation. Genes involved in heavy metal hyperaccumulation showed an expression bias in the A. halleri-derived homeolog in widely distributed plant genotypes. We also found that two populations were found growing on serpentine soils. These data suggest that A. kamchatica can inhabit a range of both natural and artificial soil environments with high levels of ions that either of the parents specializes and that it can accumulate varying amount of heavy metals. Our field and experimental data provide a compelling example of combining genetic toolkits for soil adaptations to expand the habitat of an allopolyploid species.

Parental legacy and regulatory novelty in Brachypodium diurnal transcriptomes accompanying their polyploidy.

Polyploidy is a widespread phenomenon in eukaryotes that can lead to phenotypic novelty and has important implications for evolution and diversification. The modification of phenotypes in polyploids relative to their diploid progenitors may be associated with altered gene expression. However, it is largely unknown how interactions between duplicated genes affect their diurnal expression in allopolyploid species. In this study, we explored parental legacy and hybrid novelty in the transcriptomes of an allopolyploid species and its diploid progenitors. We compared the diurnal transcriptomes of representative Brachypodium cytotypes, including the allotetraploid Brachypodium hybridum and its diploid progenitors Brachypodium distachyon and Brachypodium stacei. We also artificially induced an autotetraploid B. distachyon. We identified patterns of homoeolog expression bias (HEB) across Brachypodium cytotypes and time-dependent gain and loss of HEB in B. hybridum. Furthermore, we established that many genes with diurnal expression experienced HEB, while their expression patterns and peak times were correlated between homoeologs in B. hybridum relative to B. distachyon and B. stacei, suggesting diurnal synchronization of homoeolog expression in B. hybridum. Our findings provide insight into the parental legacy and hybrid novelty associated with polyploidy in Brachypodium, and highlight the evolutionary consequences of diurnal transcriptional regulation that accompanied allopolyploidy.

Genetic Variation for Seed Metabolite Levels in Brachypodium distachyon.

Metabolite composition and concentrations in seed grains are important traits of cereals. To identify the variation in the seed metabolotypes of a model grass, namely Brachypodium distachyon, we applied a widely targeted metabolome analysis to forty inbred lines of B. distachyon and examined the accumulation patterns of 183 compounds in the seeds. By comparing the metabolotypes with the population structure of these lines, we found signature metabolites that represent different accumulation patterns for each of the three B. distachyon subpopulations. Moreover, we found that thirty-seven metabolites exhibited significant differences in their accumulation between the lines Bd21 and Bd3-1. Using a recombinant inbred line (RIL) population from a cross between Bd3-1 and Bd21, we identified the quantitative trait loci (QTLs) linked with this variation in the accumulation of thirteen metabolites. Our metabolite QTL analysis illustrated that different genetic factors may presumably regulate the accumulation of 4-pyridoxate and pyridoxamine in vitamin B6 metabolism. Moreover, we found two QTLs on chromosomes 1 and 4 that affect the accumulation of an anthocyanin, chrysanthemin. These QTLs genetically interacted to regulate the accumulation of this compound. This study demonstrates the potential for metabolite QTL mapping in B. distachyon and provides new insights into the genetic dissection of metabolomic traits in temperate grasses.

Benzothiadiazole, a plant defense inducer, negatively regulates sheath blight resistance in Brachypodium distachyon.

Plant defense inducers that mimic functions of the plant immune hormone salicylic acid (SA) often affect plant growth. Although benzothiadiazole (BTH), a synthetic analog of SA, has been widely used to protect crops from diseases by inducing plant defense responses, we recently demonstrated that SA, but not BTH, confers resistance against Rhizoctonia solani, the causal agent of sheath blight disease, in Brachypodium distachyon. Here, we demonstrated that BTH compromised the resistance of Bd3-1 and Gaz4, the two sheath blight-resistant accessions of B. distachyon, which activate SA-dependent signaling following challenge by R. solani. Moreover, upon analyzing our published RNA-seq data from B. distachyon treated with SA or BTH, we found that BTH specifically induces expression of genes related to chloroplast function and jasmonic acid (JA) signaling, suggesting that BTH attenuates R. solani resistance by perturbing growth-defense trade-offs and/or by inducing a JA response that may increase susceptibility to R. solani. Our findings demonstrated that BTH does not work as a simple mimic of SA in B. distachyon, and consequently may presumably cause unfavorable side effects through the transcriptional alteration, particularly with respect to R. solani resistance.

Homoeolog-specific activation of genes for heat acclimation in the allopolyploid grass Brachypodium hybridum.

Background: Allopolyploid plants often show wider environmental tolerances than their ancestors; this is expected to be due to the merger of multiple distinct genomes with a fixed heterozygosity. The complex homoeologous gene expression could have been evolutionarily advantageous for the adaptation of allopolyploid plants. Despite multiple previous studies reporting homoeolog-specific gene expression in allopolyploid species, there are no clear examples of homoeolog-specific function in acclimation to a long-term stress condition. Results: We found that the allopolyploid grass Brachypodium hybridum and its ancestor Brachypodium stacei show long-term heat stress tolerance, unlike its other ancestor, Brachypodium distachyon. To understand the physiological traits of B. hybridum, we compared the transcriptome of the 3 Brachypodium species grown under normal and heat stress conditions. We found that the expression patterns of approximately 26% and approximately 38% of the homoeolog groups in B. hybridum changed toward nonadditive expression and nonancestral expression, respectively, under normal condition. Moreover, we found that B. distachyon showed similar expression patterns between normal and heat stress conditions, whereas B. hybridum and B. stacei significantly altered their transcriptome in response to heat after 3 days of stress exposure, and homoeologs that were inherited from B. stacei may have contributed to the transcriptional stress response to heat in B. hybridum. After 15 days of heat exposure, B. hybridum and B. stacei maintained transcriptional states similar to those under normal conditions. These results suggest that an earlier response to heat that was specific to homoeologs originating from B. stacei contributed to cellular homeostasis under long-term heat stress in B. hybridum. Conclusions: Our results provide insights into different regulatory events of the homoeo-transcriptome that are associated with stress acclimation in allopolyploid plants.

Multiplex PCR Targeted Amplicon Sequencing (MTA-Seq): Simple, Flexible, and Versatile SNP Genotyping by Highly Multiplexed PCR Amplicon Sequencing.

Next-generation sequencing (NGS) technologies have enabled genome re-sequencing for exploring genome-wide polymorphisms among individuals, as well as targeted re-sequencing for the rapid and simultaneous detection of polymorphisms in genes associated with various biological functions. Therefore, a simple and robust method for targeted re-sequencing should facilitate genotyping in a wide range of biological fields. In this study, we developed a simple, custom, targeted re-sequencing method, designated "multiplex PCR targeted amplicon sequencing (MTA-seq)," and applied it to the genotyping of the model grass Brachypodium distachyon. To assess the practical usability of MTA-seq, we applied it to the genotyping of genome-wide single-nucleotide polymorphisms (SNPs) identified in natural accessions (Bd1-1, Bd3-1, Bd21-3, Bd30-1, Koz-1, Koz-3, and Koz-4) by comparing the re-sequencing data with that of reference accession Bd21. Examination of SNP-genotyping accuracy in 443 amplicons from eight parental accessions and an F1 progeny derived by crossing of Bd21 and Bd3-1 revealed that ~95% of the SNPs were correctly called. The assessment suggested that the method provided an efficient framework for accurate and robust SNP genotyping. The method described here enables easy design of custom target SNP-marker panels in various organisms, facilitating a wide range of high-throughput genetic applications, such as genetic mapping, population analysis, molecular breeding, and genomic diagnostics.

Salicylic acid-dependent immunity contributes to resistance against Rhizoctonia solani, a necrotrophic fungal agent of sheath blight, in rice and Brachypodium distachyon.

Rhizoctonia solani is a soil-borne fungus causing sheath blight. In consistent with its necrotrophic life style, no rice cultivars fully resistant to R. solani are known, and agrochemical plant defense activators used for rice blast, which upregulate a phytohormonal salicylic acid (SA)-dependent pathway, are ineffective towards this pathogen. As a result of the unavailability of genetics, the infection process of R. solani remains unclear. We used the model monocotyledonous plants Brachypodium distachyon and rice, and evaluated the effects of phytohormone-induced resistance to R. solani by pharmacological, genetic and microscopic approaches to understand fungal pathogenicity. Pretreatment with SA, but not with plant defense activators used in agriculture, can unexpectedly induce sheath blight resistance in plants. SA treatment inhibits the advancement of R. solani to the point in the infection process in which fungal biomass shows remarkable expansion and specific infection machinery is developed. The involvement of SA in R. solani resistance is demonstrated by SA-deficient NahG transgenic rice and the sheath blight-resistant B. distachyon accessions, Bd3-1 and Gaz-4, which activate SA-dependent signaling on inoculation. Our findings suggest a hemi-biotrophic nature of R. solani, which can be targeted by SA-dependent plant immunity. Furthermore, B. distachyon provides a genetic resource that can confer disease resistance against R. solani to plants.

Diurnal Transcriptome and Gene Network Represented through Sparse Modeling in Brachypodium distachyon.

We report the comprehensive identification of periodic genes and their network inference, based on a gene co-expression analysis and an Auto-Regressive eXogenous (ARX) model with a group smoothly clipped absolute deviation (SCAD) method using a time-series transcriptome dataset in a model grass, Brachypodium distachyon. To reveal the diurnal changes in the transcriptome in B. distachyon, we performed RNA-seq analysis of its leaves sampled through a diurnal cycle of over 48 h at 4 h intervals using three biological replications, and identified 3,621 periodic genes through our wavelet analysis. The expression data are feasible to infer network sparsity based on ARX models. We found that genes involved in biological processes such as transcriptional regulation, protein degradation, and post-transcriptional modification and photosynthesis are significantly enriched in the periodic genes, suggesting that these processes might be regulated by circadian rhythm in B. distachyon. On the basis of the time-series expression patterns of the periodic genes, we constructed a chronological gene co-expression network and identified putative transcription factors encoding genes that might be involved in the time-specific regulatory transcriptional network. Moreover, we inferred a transcriptional network composed of the periodic genes in B. distachyon, aiming to identify genes associated with other genes through variable selection by grouping time points for each gene. Based on the ARX model with the group SCAD regularization using our time-series expression datasets of the periodic genes, we constructed gene networks and found that the networks represent typical scale-free structure. Our findings demonstrate that the diurnal changes in the transcriptome in B. distachyon leaves have a sparse network structure, demonstrating the spatiotemporal gene regulatory network over the cyclic phase transitions in B. distachyon diurnal growth.

RNA-Seq using bulked recombinant inbred line populations uncovers the importance of brassinosteroid for seed longevity after priming treatments.

Seed priming is a commercially used technique for improving seed performance including germination. However, the treatment sometimes reduces seed longevity as a side effect, limiting the storable period or longevity of the seeds. To overcome this problem, molecular mechanisms involved in the loss of seed longevity during priming were analyzed using natural variations of Arabidopsis thaliana. We found that the Est-1 accession retained longevity for longer after priming compared to the reference accession Col-0. QTL analysis using 279 recombinant inbred lines (RILs) derived from the Est-1 × Col-0 detected three QTL regions associated with the loss of seed longevity during priming. Bulked transcriptome analysis (RNA-Seq with bulked RIL populations) revealed that genes related to brassinosteroid (BR) biosynthesis/signaling and cell wall modification were highly expressed in primed seeds with shorter longevity. After priming, BR-deficient mutants cyp85a1/a2 and det2 showed significantly longer longevity than the wild type (WT). Moreover, tetrazolium staining indicated that mutant seed coats were less permeable after priming than those of WT. We suggest that the loss of seed longevity in primed seed is due to increased seed coat permeability, which is positively regulated, at least partly, via BR signaling.

Exploring Genetic Diversity in Plants Using High-Throughput Sequencing Techniques.

Food security has emerged as an urgent concern because of the rising world population. To meet the food demands of the near future, it is required to improve the productivity of various crops, not just of staple food crops. The genetic diversity among plant populations in a given species allows the plants to adapt to various environmental conditions. Such diversity could therefore yield valuable traits that could overcome the food-security challenges. To explore genetic diversity comprehensively and to rapidly identify useful genes and/or allele, advanced high-throughput sequencing techniques, also called next-generation sequencing (NGS) technologies, have been developed. These provide practical solutions to the challenges in crop genomics. Here, we review various sources of genetic diversity in plants, newly developed genetic diversity-mining tools synergized with NGS techniques, and related genetic approaches such as quantitative trait locus analysis and genome-wide association study.

Conserved but Attenuated Parental Gene Expression in Allopolyploids: Constitutive Zinc Hyperaccumulation in the Allotetraploid Arabidopsis kamchatica.

Allopolyploidization combines parental genomes and often confers broader species distribution. However, little is known about parentally transmitted gene expression underlying quantitative traits following allopolyploidization because of the complexity of polyploid genomes. The allopolyploid species Arabidopsis kamchatica is a natural hybrid of the zinc hyperaccumulator Arabidopsis halleri and of the nonaccumulator Arabidopsis lyrata We found that A. kamchatica retained the ability to hyperaccumulate zinc from A. halleri and grows in soils with both low and high metal content. Hyperaccumulation of zinc by A. kamchatica was reduced to about half of A. halleri, but is 10-fold greater than A. lyrata Homeologs derived from A. halleri had significantly higher levels of expression of genes such as HEAVY METAL ATPASE4 (HMA4), METAL TRANSPORTER PROTEIN1 and other metal ion transporters than those derived from A. lyrata, which suggests cis-regulatory differences. A. kamchatica has on average about half the expression of these genes compared with A. halleri due to fixed heterozygosity inherent in allopolyploids. Zinc treatment significantly changed the ratios of expression of 1% of homeologous pairs, including genes putatively involved in metal homeostasis. Resequencing data showed a significant reduction in genetic diversity over a large genomic region (290 kb) surrounding the HMA4 locus derived from the A. halleri parent compared with the syntenic A. lyrata-derived region, which suggests different evolutionary histories. We also estimated that three A. halleri-derived HMA4 copies are present in A. kamchatica Our findings support a transcriptomic model in which environment-related transcriptional patterns of both parents are conserved but attenuated in the allopolyploids.

Expression profiling of marker genes responsive to the defence-associated phytohormones salicylic acid, jasmonic acid and ethylene in Brachypodium distachyon.

BACKGROUND: Brachypodium distachyon is a promising model plants for grasses. Infections of Brachypodium by various pathogens that severely impair crop production have been reported, and the species accordingly provides an alternative platform for investigating molecular mechanisms of pathogen virulence and plant disease resistance. To date, we have a broad picture of plant immunity only in Arabidopsis and rice; therefore, Brachypodium may constitute a counterpart that displays the commonality and uniqueness of defence systems among plant species. Phytohormones play key roles in plant biotic stress responses, and hormone-responsive genes are used to qualitatively and quantitatively evaluate disease resistance responses during pathogen infection. For these purposes, defence-related phytohormone marker genes expressed at time points suitable for defence-response monitoring are needed. Information about their expression profiles over time as well as their response specificity is also helpful. However, useful marker genes are still rare in Brachypodium. RESULTS: We selected 34 candidates for Brachypodium marker genes on the basis of protein-sequence similarity to known marker genes used in Arabidopsis and rice. Brachypodium plants were treated with the defence-related phytohormones salicylic acid, jasmonic acid and ethylene, and their transcription levels were measured 24 and 48 h after treatment. Two genes for salicylic acid, 7 for jasmonic acid and 2 for ethylene were significantly induced at either or both time points. We then focused on 11 genes encoding pathogenesis-related (PR) 1 protein and compared their expression patterns with those of Arabidopsis and rice. Phylogenetic analysis suggested that Brachypodium contains several PR1-family genes similar to rice genes. Our expression profiling revealed that regulation patterns of some PR1 genes as well as of markers identified for defence-related phytohormones are closely related to those in rice. CONCLUSION: We propose that the Brachypodium immune hormone marker genes identified in this study will be useful to plant pathologists who use Brachypodium as a model pathosystem, because the timing of their transcriptional activation matches that of the disease resistance response. Our results using Brachypodium also suggest that monocots share a characteristic immune system, defined as the common defence system, that is different from that of dicots.

Determination of growth stages and metabolic profiles in Brachypodium distachyon for comparison of developmental context with Triticeae crops.

Brachypodium distachyon is an emerging model plant for studying biological phenomena in temperate grasses. Study of the growth scale is essential to analyse spatio-temporal changes in molecular factors throughout the life cycle. For sensitive and robust staging based on morphology in B. distachyon, we demonstrated the utility of the BBCH (Biologische Bundesanstalt, Bundessortenamt and CHemical industry) scale, which is comparable to the Zadoks scale conventionally used for Triticeae crops. We compared the chronological progression of B. distachyon accessions Bd21 and Bd3-1, in addition to the progression of Chinese Spring wheat. The comparison of growth stages illustrates the morphological similarities and differences in the timing of life cycle events. Furthermore, we compared metabolite accumulation patterns across different growth stages and across different stress conditions using a widely targeted metabolome analysis. Metabolic profiling determined commonalities and specificities in chemical properties that were dependent on organisms, growth stages and/or stress conditions. Most metabolites accumulated equivalently in B. distachyon and wheat. This qualitative similarity indicated the superiority of B. distachyon as a model for Triticeae crops. The growth scale of B. distachyon should provide a conceptual framework for comparative analysis and for knowledge integration between this model grass and crops in the Pooideae subfamily.

Transcriptome analysis of thermogenic Arum concinnatum reveals the molecular components of floral scent production.

Several plant species can generate enough heat to increase their internal floral temperature above ambient temperature. Among thermogenic plants, Arum concinnatum shows the highest respiration activity during thermogenesis. However, an overall understanding of the genes related to plant thermogenesis has not yet been achieved. In this study, we performed de novo transcriptome analysis of flower organs in A. concinnatum. The de novo transcriptome assembly represented, in total, 158,490 non-redundant transcripts, and 53,315 of those showed significant homology with known genes. To explore genes associated with thermogenesis, we filtered 1266 transcripts that showed a significant correlation between expression pattern and the temperature trend of each sample. We confirmed five putative alternative oxidase transcripts were included in filtered transcripts as expected. An enrichment analysis of the Gene Ontology terms for the filtered transcripts suggested over-representation of genes involved in 1-deoxy-D-xylulose-5-phosphate synthase (DXS) activity. The expression profiles of DXS transcripts in the methyl-D-erythritol 4-phosphate (MEP) pathway were significantly correlated with thermogenic levels. Our results suggest that the MEP pathway is the main biosynthesis route for producing scent monoterpenes. To our knowledge, this is the first report describing the candidate pathway and the key enzyme for floral scent production in thermogenic plants.

Effects of floral thermogenesis on pollen function in Asian skunk cabbage Symplocarpus renifolius.

The effects of temperature on pollen germination and pollen tube growth rate were measured in vitro in thermogenic skunk cabbage, Symplocarpus renifolius Schott ex Tzvelev, and related to floral temperatures in the field. This species has physiologically thermoregulatory spadices that maintain temperatures near 23 degrees C, even in sub-freezing air. Tests at 8, 13, 18, 23, 28 and 33 degrees C showed sharp optima at 23 degrees C for both variables, and practically no development at 8 degrees C. Thermogenesis is therefore a requirement for fertilization in early spring. The narrow temperature tolerance is probably related to a long period of evolution in flowers that thermoregulate within a narrow range.

Two cys or not two cys? That is the question; alternative oxidase in the thermogenic plant sacred Lotus.

Sacred lotus (Nelumbo nucifera) regulates temperature in its floral chamber to 32 degrees C to 35 degrees C across ambient temperatures of 8 degrees C to 40 degrees C with heating achieved through high alternative pathway fluxes. In most alternative oxidase (AOX) isoforms, two cysteine residues, Cys(1) and Cys(2), are highly conserved and play a role in posttranslational regulation of AOX. Further control occurs via interaction of reduced Cys(1) with alpha-keto acids, such as pyruvate. Here, we report on the in vitro regulation of AOX isolated from thermogenic receptacle tissues of sacred lotus. AOX protein was mostly present in the reduced form, and only a small fraction could be oxidized with diamide. Cyanide-resistant respiration in isolated mitochondria was stimulated 4-fold by succinate but not pyruvate or glyoxylate. Insensitivity of the alternative pathway of respiration to pyruvate and the inability of AOX protein to be oxidized by diamide suggested that AOX in these tissues may lack Cys(1). Subsequently, we isolated two novel cDNAs for AOX from thermogenic tissues of sacred lotus, designated as NnAOX1a and NnAOX1b. Deduced amino acid sequences of both confirmed that Cys(1) had been replaced by serine; however, Cys(2) was present. This contrasts with AOXs from thermogenic Aroids, which contain both Cys(1) and Cys(2). An additional cysteine was present at position 193 in NnAOX1b. The significance of the sequence data for regulation of the AOX protein in thermogenic sacred lotus is discussed and compared with AOXs from other thermogenic and nonthermogenic species.

Functional coexpression of the mitochondrial alternative oxidase and uncoupling protein underlies thermoregulation in the thermogenic florets of skunk cabbage.

Two distinct mitochondrial energy dissipating systems, alternative oxidase (AOX) and uncoupling protein (UCP), have been implicated as crucial components of thermogenesis in plants and animals, respectively. To further clarify the physiological roles of AOX and UCP during homeothermic heat production in the thermogenic skunk cabbage (Symplocarpus renifolius), we identified the thermogenic cells and performed expression and functional analyses of these genes in this organism. Thermographic analysis combined with in situ hybridization revealed that the putative thermogenic cells surround the stamens in the florets of skunk cabbage and coexpress transcripts for SrAOX, encoding Symplocarpus AOX, and SrUCPb, encoding a novel UCP that lacks a fifth transmembrane segment. Mitochondria isolated from the thermogenic florets exhibited substantial linoleic acid (LA)-inducible uncoupling activities. Moreover, our results demonstrate that LA is capable of inhibiting the mitochondrial AOX pathway, whereas the proportion of pyruvate-stimulated AOX capacity was not significantly affected by LA. Intriguingly, the protein expression levels for SrAOX and SrUCPb were unaffected even when the ambient air temperatures increased from 10.3 degrees C to 23.1 degrees C or from 8.3 degrees C to 24.9 degrees C. Thus, our results suggest that functional coexpression of AOX and UCP underlies the molecular basis of heat production, and that posttranslational modifications of these proteins play a crucial role in regulating homeothermic heat production under conditions of natural ambient temperature fluctuations in skunk cabbage.

Pyruvate-sensitive AOX exists as a non-covalently associated dimer in the homeothermic spadix of the skunk cabbage, Symplocarpus renifolius.

The cyanide-resistant alternative oxidase (AOX) is a homodimeric protein whose activity can be regulated by the oxidation/reduction state and by alpha-keto acids. To further clarify the role of AOX in the skunk cabbage, Symplocarpus renifolius, we have performed expression and functional analyses of the encoding gene. Among the various tissues in the skunk cabbage, SrAOX transcripts were found to be specifically expressed in the thermogenic spadix. Moreover, our data demonstrate that the SrAOX protein exists as a non-covalently associated dimer in the thermogenic spadix, and is more sensitive to pyruvate than to other carboxylic acids. Our results suggest that the pyruvate-mediated modification of SrAOX activity plays a significant role in thermoregulation in the skunk cabbage.

Modeling of the thermoregulation system in the skunk cabbage: Symplocarpus foetidus.

This paper presents a model of the thermoregulation system of the spadix of skunk cabbage Symplocarpus foetidus which regulates its internal temperature at around 20 degrees C during flowering even when the ambient air temperature drops below freezing. From the temperature responses of the spadix to changing ambient air temperature, we assumed that the thermoregulation system of the spadix is probably one of negative feedback control. The feedback signals are based on the rate of temperature change of the spadix over time. A signal is factored into the biochemical energy generator, and becomes biochemical energy, some of which becomes heat. Comparing our proposed model temperature responses and those of the living spadix, we found good agreement. In the process of engineering the model, the existence of two regulatory pathways in the thermoregulation system was simulated, and our proposed model appears to provide the necessary elements to explain the fundamental mechanism of the thermoregulation system of S. foetidus.

Changes in the composition of xylem sap during development of the spadix of skunk cabbage (Symplocarpus foetidus).

The spadix of skunk cabbage, Symplocarpus foetidus, is thermogenic and maintains an internal temperature of around 20 degrees C even when the ambient air temperature drops below freezing. This homeothermic heat production is observed only during the stigma stage, and thereafter ceases at the male stage when pollen is shed. To clarify the regulatory mechanism by which the stigma stage-specific heat production occurs in the spadix, sugars, organic acids, and amino acids in xylem sap were analyzed and compared with those of post-thermogenic plants. Interestingly, no significant difference was observed in the total volume of xylem sap per fresh weight of the spadix between thermogenic (31.2+/-24.7 microl h(-1) g(-1)) and post-thermogenic (50.5+/-30.4 microl h(-1) g(-1)) plants. However, concentrations of sugars (sucrose, glucose, and fructose), organic acids (malate and succinate), and amino acids (Asp, Asn, Glu, Gln, Gly, and Ala) in xylem sap decreased remarkably in post-thermogenic plants. Our results indicate that the composition of the xylem sap differs during the development of the spadix of S. foetidus.