<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//TaxonX//DTD Taxonomic Treatment Publishing DTD v0 20100105//EN" "https://plecevo.eu/nlm/tax-treatment-NS0.dtd">
<article xmlns:tp="http://www.plazi.org/taxpub" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:mml="http://www.w3.org/1998/Math/MathML" article-type="research-article" xml:lang="en">
  <front>
    <journal-meta>
      <journal-id journal-id-type="publisher-id">118</journal-id>
      <journal-id journal-id-type="index">urn:lsid:arphahub.com:pub:71cc5dc6-a767-5334-951f-ef6ae8936459</journal-id>
      <journal-title-group>
        <journal-title xml:lang="en">Plant Ecology and Evolution</journal-title>
        <abbrev-journal-title xml:lang="en">plecevo</abbrev-journal-title>
      </journal-title-group>
      <issn pub-type="ppub">2032-3913</issn>
      <issn pub-type="epub">2032-3921</issn>
      <publisher>
        <publisher-name>Meise Botanic Garden and Royal Botanical Society of Belgium</publisher-name>
      </publisher>
    </journal-meta>
    <article-meta>
      <article-id pub-id-type="doi">10.5091/plecevo.177783</article-id>
      <article-id pub-id-type="publisher-id">177783</article-id>
      <article-categories>
        <subj-group subj-group-type="heading">
          <subject>Research Article</subject>
        </subj-group>
        <subj-group subj-group-type="biological_taxon">
          <subject>Araceae</subject>
          <subject>Basal Angiosperms: Monocots</subject>
        </subj-group>
        <subj-group subj-group-type="scientific_subject">
          <subject>Ecology</subject>
          <subject>Genetics</subject>
          <subject>Morphology &amp; Anatomy</subject>
        </subj-group>
        <subj-group subj-group-type="geographical_area">
          <subject>Asia</subject>
        </subj-group>
      </article-categories>
      <title-group>
        <article-title>Does size matter? The integrated roles of light, adaptive sex expression, and hybridization in a widespread <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part></tp:taxon-name></italic> (<tp:taxon-name><tp:taxon-name-part taxon-name-part-type="family" reg="Araceae">Araceae</tp:taxon-name-part></tp:taxon-name>) species from western Japan</article-title>
      </title-group>
      <contrib-group content-type="authors">
        <contrib contrib-type="author" corresp="yes">
          <name name-style="western">
            <surname>Scholten</surname>
            <given-names>Justin T.</given-names>
          </name>
          <email xlink:type="simple">jts329@cornell.edu</email>
          <xref ref-type="aff" rid="A1">1</xref>
          <role content-type="http://credit.niso.org/contributor-roles/conceptualization/">Conceptualization</role>
          <role content-type="http://credit.niso.org/contributor-roles/writing-original-draft/">Writing - original draft</role>
          <role content-type="http://credit.niso.org/contributor-roles/formal-analysis/">Formal analysis</role>
        </contrib>
        <contrib contrib-type="author" corresp="no">
          <name name-style="western">
            <surname>Specht</surname>
            <given-names>Chelsea D.</given-names>
          </name>
          <xref ref-type="aff" rid="A1">1</xref>
        </contrib>
      </contrib-group>
      <aff id="A1">
        <label>1</label>
        <addr-line content-type="verbatim">Section of Plant Biology and the LH Bailey Hortorium, Cornell University, Ithaca, New York, United States of America</addr-line>
        <institution>Section of Plant Biology and the LH Bailey Hortorium, Cornell University</institution>
        <addr-line content-type="city">Ithaca</addr-line>
        <country>United States of America</country>
        <uri content-type="ror">https://ror.org/05bnh6r87</uri>
      </aff>
      <author-notes>
        <fn fn-type="corresp">
          <p>Corresponding author: Justin T. Scholten (<email xlink:type="simple">jts329@cornell.edu</email>)</p>
        </fn>
        <fn fn-type="edited-by">
          <p><bold>Academic editor</bold>: Myriam Heuertz</p>
        </fn>
      </author-notes>
      <pub-date pub-type="collection">
        <year>2026</year>
      </pub-date>
      <pub-date pub-type="epub">
        <day>14</day>
        <month>04</month>
        <year>2026</year>
      </pub-date>
      <volume>159</volume>
      <issue>2</issue>
      <fpage>197</fpage>
      <lpage>207</lpage>
      <uri content-type="arpha" xlink:href="http://openbiodiv.net/AED69A8F-A563-5B18-ACE3-C62A21DA478B">AED69A8F-A563-5B18-ACE3-C62A21DA478B</uri>
      <history>
        <date date-type="received">
          <day>11</day>
          <month>11</month>
          <year>2025</year>
        </date>
        <date date-type="accepted">
          <day>06</day>
          <month>01</month>
          <year>2026</year>
        </date>
      </history>
      <permissions>
        <copyright-statement>Justin T. Scholten, Chelsea D. Specht</copyright-statement>
        <license license-type="creative-commons-attribution" xlink:href="http://creativecommons.org/licenses/by/4.0/" xlink:type="simple">
          <license-p>This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</license-p>
        </license>
      </permissions>
      <abstract>
        <label>Abstract</label>
        <p><bold>Background and aims</bold> – Reproductive plasticity and hybridization are known to play key roles in local adaptation and range expansion in plants, yet the environmental conditions under which canonical size-dependent models of sex expression break down remain poorly understood. <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> (<tp:taxon-name><tp:taxon-name-part taxon-name-part-type="family" reg="Araceae">Araceae</tp:taxon-name-part></tp:taxon-name>), a widespread species in Japan, has long been characterized as exhibiting size-dependent sex determination, whereby individuals transition from producing staminate to pistillate inflorescences upon reaching a size threshold. Field observations from western Japan further suggest that sex expression in <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> may be decoupled from plant size, raising the possibility that environmental cues such as light availability override intrinsic size thresholds in some populations. Because such ecological flexibility can increase spatial and temporal overlap among sympatric congeners, we additionally evaluate whether relaxed sex determination may coincide with increased introgression within this species complex.</p>
        <p><bold>Materials and methods</bold> – We conducted field-based demographic surveys across five populations of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part></tp:taxon-name></italic> on Shikoku Island, Japan, quantifying plant size, sex expression, and light intensity to evaluate environmental and morphological predictors of reproductive phase. To assess hybridization, we generated genome-wide SNP data through genotyping-by-sequencing, constructed a maximum likelihood phylogeny, and performed f-branch analyses to detect patterns of introgression among sympatric species.</p>
        <p><bold>Key results</bold> – Our data indicate that size in fact may not be the key factor in sex determination, with light availability emerging as a stronger predictor of sex expression than size alone for populations on Shikoku Island. Preliminary genomic analyses support evidence of gene flow between <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> and sympatric species, indicating ongoing hybridization and introgression.</p>
        <p><bold>Conclusions</bold> – Our findings indicate that the shift away from a strict allometric model of sex determination on Shikoku Island to a photosensitive model allowing for diphasic sexual expression likely enhances reproductive flexibility and leads to an increase in seed production in open-canopy and disturbed environments. Together, these findings provide ecological and demographic insight contributing to the broad distribution, sex determination, and morphological diversity of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic>.</p>
      </abstract>
      <kwd-group>
        <label>Keywords</label>
        <kwd>
          <tp:taxon-name>
            <tp:taxon-name-part taxon-name-part-type="family" reg="Araceae">Araceae</tp:taxon-name-part>
          </tp:taxon-name>
        </kwd>
        <kwd>colonization</kwd>
        <kwd>diphasic sexual expression</kwd>
        <kwd>niche expansion</kwd>
        <kwd>plasticity</kwd>
        <kwd>Shikoku Island</kwd>
      </kwd-group>
      <funding-group>
        <funding-statement>This work was supported by funds provided by the Lewis and Clark Fund for Exploration and Research, the Torrey Botanical Society, the Society of Systematic Biologists, the Mario Einaudi Center of International Studies, and the Moore Fund of the L.H. Bailey Hortorium Herbarium. Additionally, the stipend of the first author was supported in part by the Tara Atluri Memorial Fund.</funding-statement>
      </funding-group>
    </article-meta>
  </front>
  <body>
    <sec sec-type="Introduction" id="sec1">
      <title>Introduction</title>
      <p>Species of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part></tp:taxon-name></italic> Mart. are perennial, tuberous geophytes that occur primarily as forest understory herbs, persisting across years via a subterranean corm but relying largely on sexual reproduction rather than clonal spread (<xref ref-type="bibr" rid="B4">Bierzychudek 1982</xref>; <xref ref-type="bibr" rid="B32">Murata 2011</xref>; <xref ref-type="bibr" rid="B35">Murata et al. 2018</xref>). Although typically associated with shaded forest interiors, many taxa, including <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> (Thunb.) Schott, readily establish in canopy gaps, forest edges, and managed plantation forests where light availability is spatially heterogeneous (<xref ref-type="bibr" rid="B18">Kinoshita 1987</xref>; <xref ref-type="bibr" rid="B29">Masaki et al. 2017</xref>; <xref ref-type="bibr" rid="B30">Matsumoto et al. 2018</xref>). Pollination is mediated mainly by small dipterans attracted and temporarily trapped within the spathe, a system that can facilitate heterospecific pollen transfer when congeners flower synchronously, while seeds are dispersed locally by gravity and small animals into newly opened microsites (<xref ref-type="bibr" rid="B1">Atkinson 1898</xref>; <xref ref-type="bibr" rid="B4">Bierzychudek 1982</xref>; <xref ref-type="bibr" rid="B33">Murata and Ohno 1989</xref>; <xref ref-type="bibr" rid="B22">Kudo et al. 2008</xref>).</p>
      <p>The development of sexual expression in plants, here defined as the production of male or female flowers by an individual, is often constrained by intrinsic factors such as plant size and resource availability (<xref ref-type="bibr" rid="B40">Policansky 1981</xref>). In sexually diphasic species, these constraints can influence or trigger transitions between reproductive phases. In <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part></tp:taxon-name></italic> (<tp:taxon-name><tp:taxon-name-part taxon-name-part-type="family" reg="Araceae">Araceae</tp:taxon-name-part></tp:taxon-name>), sex expression has traditionally been considered size-dependent, with smaller individuals producing staminate (male) flowers and larger individuals producing “pistillate” or ovulate (female) flowers, due to the higher energetic demands associated with female function (<xref ref-type="bibr" rid="B40">Policansky 1981</xref>; <xref ref-type="bibr" rid="B4">Bierzychudek 1982</xref>). This pattern is supported by early morphological and demographic studies in North American congeners such as <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="triphyllum">triphyllum</tp:taxon-name-part></tp:taxon-name></italic> (L.) Schott and <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="dracontium">dracontium</tp:taxon-name-part></tp:taxon-name></italic> (L.) Schott, where sex changes are linked to plant size (<xref ref-type="bibr" rid="B1">Atkinson 1898</xref>; <xref ref-type="bibr" rid="B5">Camp 1932</xref>; <xref ref-type="bibr" rid="B13">Doust and Cavers 1982</xref>; <xref ref-type="bibr" rid="B16">Ewing and Klein 1982</xref>; <xref ref-type="bibr" rid="B9">Clay 1993</xref>), nutrient status (<xref ref-type="bibr" rid="B44">Schaffner 1922</xref>), and age (<xref ref-type="bibr" rid="B17">Gow 1913</xref>; <xref ref-type="bibr" rid="B5">Camp 1932</xref>). Similar trends have been documented in several Japanese <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part></tp:taxon-name></italic> species. For example, individuals of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="urashima">urashima</tp:taxon-name-part></tp:taxon-name></italic> H.Hara are reported to transition from male to female flowering with increasing plant size (<xref ref-type="bibr" rid="B48">Takasu 1987</xref>). The same size-dependent model has been applied to <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="japonicum">japonicum</tp:taxon-name-part></tp:taxon-name></italic> Blume (now considered a synonym of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic>), with early studies noting a general male-to-female transition in flowering accompanied by increases in plant age and plant size under controlled settings (<xref ref-type="bibr" rid="B25">Maekawa 1924</xref>, <xref ref-type="bibr" rid="B26">1927</xref>). However, these studies were conducted under homogeneous greenhouse conditions, limiting their ability to assess the influence of natural environmental gradients.</p>
      <p>More recent field-based work by <xref ref-type="bibr" rid="B18">Kinoshita (1987)</xref> provides explicit support for a size–sex relationship in <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> populations from Nagano Prefecture in north-central Japan. This study documented that larger plants were more likely to produce female flowers, consistent with the canonical size-dependent model. Yet, this pattern may not be universal. Preliminary field observations of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> populations on Shikoku Island in western Japan suggested an absence or weakening of this correlation, with female-flowering individuals frequently observed across a broad range of plant sizes, including individuals below traditionally inferred size thresholds. These observations prompted a re-evaluation of the relative roles of intrinsic (size-based) and extrinsic (environmentally induced) cues in determining sex expression in this variable species.</p>
      <p>The size-dependent model of sexual expression has been widely documented across plant families and is thought to optimize reproductive success under stable ecological conditions (<xref ref-type="bibr" rid="B7">Charnov and Bull 1977</xref>). However, in disturbed environments, where rapid colonization is favoured, strict size-dependent sex expression may be maladaptive (<xref ref-type="bibr" rid="B12">Dorken and Barrett 2004</xref>). In these situations, environmental parameters might select for, or result in the selection of, greater plasticity in reproductive strategies (<xref ref-type="bibr" rid="B2">Barot et al. 2005</xref>) including the timing of transition between the production of male vs female flowers.</p>
      <p>Disturbance regimes, particularly those affecting canopy structure and light availability, can substantially alter reproductive dynamics in forest understory plants by modifying patterns of resource allocation and growth (<xref ref-type="bibr" rid="B41">Protchazka et al. 2024</xref>). In species capable of facultative or environmentally sensitive sex expression, external cues such as light intensity may override or interact with intrinsic, size-dependent factors to determine reproductive phase transitions (<xref ref-type="bibr" rid="B36">Obeso 2002</xref>). Light availability is especially critical in shaping fitness within early successional or plantation forests, where fluctuating microsite conditions favour rapid growth and reproductive responsiveness (<xref ref-type="bibr" rid="B24">Lichstein et al. 2021</xref>). In such contexts, sexually diphasic plants that relax their dependence on intrinsic size thresholds and respond directly to environmental conditions may gain a selective advantage, particularly in shaded or periodically disturbed landscapes common throughout western Japan.</p>
      <p>Across many plant lineages, ecological expansion into novel or transitional habitats can also increase opportunities for interspecific hybridization. Species capable of colonizing light-variable, post-disturbance environments often come into secondary contact with congeners occupying adjacent native forests, creating conditions conducive to gene flow and introgression (<xref ref-type="bibr" rid="B15">Ellstrand and Schierenbeck 2000</xref>; <xref ref-type="bibr" rid="B43">Rieseberg et al. 2007</xref>). When such contact zones coincide with individuals exhibiting flexible or asynchronous sex expression, the likelihood of heterospecific pollen transfer may increase, facilitating admixture and the emergence of hybrid lineages.</p>
      <p>The interaction between environmentally responsive sex expression and ecological disturbance may therefore have broader evolutionary implications. Morphological observations alone cannot determine whether apparent intermediacy among sympatric <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part></tp:taxon-name></italic> taxa reflects phenotypic plasticity, shared ancestral variation, or genetic introgression, limiting inferences based solely on form. Because environmentally responsive sex expression may increase flowering synchrony and microsite overlap among congeners under variable light regimes, we evaluate introgression not as a predictor of sex expression itself, but as a potential evolutionary outcome of this reproductive flexibility in the plantation and secondary-forest landscapes of western Japan.</p>
      <p>In this study, we test how environmental and morphological factors interact to shape reproductive phase transitions in <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic>, a widespread Japanese species commonly found in post-disturbance evergreen plantations and secondary forests (<xref ref-type="bibr" rid="B30">Matsumoto et al. 2018</xref>). Specifically, we evaluate whether light availability serves as a stronger predictor of sex expression than plant size—challenging the traditional size-dependent model—and whether this ecological flexibility may contribute to the species’ success in dynamic, human-modified landscapes. To explore potential evolutionary consequences of this flexibility, we also analyse genome-wide SNP data as an exploratory framework for detecting possible gene flow between <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> and sympatric congeners. Together, these approaches aim to clarify how environmentally responsive reproductive strategies may facilitate both local adaptation and the maintenance of porous species boundaries in forest understory taxa.</p>
    </sec>
    <sec sec-type="materials|methods" id="sec2">
      <title>Material and methods</title>
      <sec sec-type="Sampling" id="sec3">
        <title>Sampling</title>
        <p>Field sampling was conducted from 2–20 May 2024 across five populations on north-western Shikoku Island, Japan (Kugawa (n = 71), Ishizuchi (n = 24), Kumakogen (n = 95), Kamigawa (n = 84), and Ozu (n = 31)). These sites were selected to encompass the regional range of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> and closely related congeners occurring in managed plantation and secondary-forest habitats. All populations were situated within mixed evergreen conifer plantations dominated by <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Cryptomeria">Cryptomeria</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="japonica">japonica</tp:taxon-name-part></tp:taxon-name></italic> (Thunb. ex L.f.) D.Don and <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Chamaecyparis">Chamaecyparis</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="obtusa">obtusa</tp:taxon-name-part></tp:taxon-name></italic> (Siebold &amp; Zucc.) Endl., landscapes shaped by post-World War II reforestation and periodic thinning. These management practices create heterogeneous light environments with canopy gaps and intermittent ground disturbance, conditions representative of much of western Japan’s low- to mid-elevation forests.</p>
        <p>At each site, all flowering <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part></tp:taxon-name></italic> individuals encountered along 100–250 m transects were recorded and identified to species following <xref ref-type="bibr" rid="B32">Murata (2011)</xref> and <xref ref-type="bibr" rid="B35">Murata et al. (2018)</xref>. Sampling captured the full range of vegetative and reproductive morphologies present within each population. The five study sites spanned elevations from approximately 250 m (Ozu) to 1050 m (Ishizuchi) and included both interior forest and edge microhabitats. All sites occurred within mixed evergreen conifer plantations dominated by <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Cryptomeria">Cryptomeria</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="japonica">japonica</tp:taxon-name-part></tp:taxon-name></italic> and <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Chamaecyparis">Chamaecyparis</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="obtusa">obtusa</tp:taxon-name-part></tp:taxon-name></italic>, with occasional secondary broadleaf elements such as <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Castanopsis">Castanopsis</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="cuspidata">cuspidata</tp:taxon-name-part></tp:taxon-name></italic> Schottky and <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Acer">Acer</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="sieboldianum">sieboldianum</tp:taxon-name-part></tp:taxon-name></italic> Miq. Across all localities, we measured and recorded a total of 305 flowering individuals (<italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> (n = 114), <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ehimense">ehimense</tp:taxon-name-part></tp:taxon-name></italic> J.Murata &amp; J.Ohno (n = 56), <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="tosaense">tosaense</tp:taxon-name-part></tp:taxon-name></italic> Makino (n = 86), <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ringens">ringens</tp:taxon-name-part></tp:taxon-name></italic> (Thunb.) Schott. (n = 30), and <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="iyoanum">iyoanum</tp:taxon-name-part></tp:taxon-name></italic> Makino (n = 19).</p>
      </sec>
      <sec sec-type="Morphological and light availability measurements" id="sec4">
        <title>Morphological and light availability measurements</title>
        <p>For each flowering individual, plant height (cm; measured from soil level to the apex of the tallest leaflet) and pseudostem base diameter (cm; measured at the midpoint of the pseudostem formed by petiole sheaths) were recorded as indicators of vegetative size from live plants. Within each species and site, we aimed to capture a wide range of morphological variation, and sampling included individuals spanning the full observed spectrum of plant height and pseudostem diameter present at the time of survey, given that an inflorescence was present for sexual classification. This approach minimized sampling bias by size class by ensuring that both smaller and larger sexually mature individuals were proportionally represented within each population.</p>
        <p>To estimate microsite light availability, photosynthetically relevant light intensity (kilolux) was measured at each individual’s position using a URCERI MT-912 digital illuminance meter. Measurements were taken on clear days between 10:00 and 14:00 local time to minimize temporal variation due to sun angle or cloud cover. The sensor was positioned 30 cm above the soil surface and oriented parallel to the plant lamina to approximate the ambient light environment experienced by each individual.</p>
        <p>Individuals were categorized as staminate (male) or pistillate (female) based on the presence of respective floral morphs within the inflorescence. Only individuals bearing mature inflorescences were included in the analyses.</p>
      </sec>
      <sec sec-type="Statistical analyses" id="sec5">
        <title>Statistical analyses</title>
        <p>To evaluate the relative influence of size and light on sex expression, logistic regression models were implemented for each species using the glm() function in R v.4.2.1 (<xref ref-type="bibr" rid="B42">R Core Team 2021</xref>) with a binomial error distribution and logit link function:</p>
        <p>sex ~ height + pseudostem diameter + light intensity</p>
        <p>Sex was coded as a binary response variable (male = 0, female = 1). Model performance was assessed using Akaike’s Information Criterion (<abbrev xlink:title="Akaike’s Information Criterion">AIC</abbrev>) and pseudo-R<sup>2</sup>. Pairwise correlations among predictor variables (height × light, diameter × light) were tested using Pearson’s r to ensure predictors were not collinear. Light availability was not strongly correlated with either plant height or pseudostem diameter (Pearson’s r ≤ 0.11; p ≥ 0.28 across species), indicating low collinearity and supporting the inclusion of these predictors within the same model. All tests were evaluated at α = 0.05. Models were specified a priori to include all biologically relevant predictors (plant height, pseudostem diameter, and light intensity), and inference focused on evaluating the relative contribution of each variable within a single multivariate framework rather than on stepwise model building or sequential model comparison.</p>
      </sec>
      <sec sec-type="Genomic data generation" id="sec6">
        <title>Genomic data generation</title>
        <p>To complement the ecological analyses, we generated exploratory genomic data to assess potential genetic relationships and gene flow among sympatric <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part></tp:taxon-name></italic> species. Fresh leaf tissue samples were collected from individuals of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ternatipartitum">ternatipartitum</tp:taxon-name-part></tp:taxon-name></italic> Makino (n = 1), <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ringens">ringens</tp:taxon-name-part></tp:taxon-name></italic> (n = 1), <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ishizuchiense">ishizuchiense</tp:taxon-name-part></tp:taxon-name></italic> Murata (n = 1), <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="sikokianum">sikokianum</tp:taxon-name-part></tp:taxon-name></italic> Franch. &amp; Sav. (n = 1), <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="iyoanum">iyoanum</tp:taxon-name-part></tp:taxon-name></italic> (n = 1), <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> (n = 3), <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ehimense">ehimense</tp:taxon-name-part></tp:taxon-name></italic> (n = 2), and <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="tosaense">tosaense</tp:taxon-name-part></tp:taxon-name></italic> (n = 1) and dried in silica. DNA was extracted using a CTAB-based method (<xref ref-type="bibr" rid="B14">Doyle and Doyle 1987</xref>), and fragments within the 5–10 kb range were selected using SPRIselect beads (<xref ref-type="bibr" rid="B3">Beckman Coulter 2012</xref>) through Solid Phase Reversible Immobilization (<abbrev xlink:title="Solid Phase Reversible Immobilization">SPRI</abbrev>). Because of limited sample sizes, these genomic data are treated as exploratory rather than population-level estimates, intended to reveal qualitative patterns of relatedness and introgression rather than to infer demographic parameters. Likewise, we did not perform PCA or model-based genetic structure analyses, which require broader population-level sampling to yield stable and interpretable results.</p>
        <p>Genotyping-by-sequencing (<abbrev xlink:title="Genotyping-by-sequencing">GBS</abbrev>) libraries were generated with the ApeKI restriction enzyme. Library preparation and sequencing were conducted at the University of Wisconsin Biotechnology Center DNA Sequencing Facility on an Illumina NovaSeq 6000 platform, producing 150 bp paired-end reads. Raw sequencing data were demultiplexed and analysed using Stacks (<xref ref-type="bibr" rid="B6">Catchen et al. 2013</xref>). All generated sequence data are publicly available on Zenodo (<xref ref-type="bibr" rid="B46">Scholten and Specht 2026</xref>).</p>
        <p>Sequence data from all individuals were pre-processed using fastp (<xref ref-type="bibr" rid="B8">Chen 2023</xref>) before being aligned to the draft genome of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="tosaense">tosaense</tp:taxon-name-part></tp:taxon-name></italic> (<xref ref-type="bibr" rid="B45">Scholten and Specht 2025</xref>) with the Burrows-Wheeler Aligner (<abbrev xlink:title="Burrows-Wheeler Aligner">BWA</abbrev> v.0.07.18) (<xref ref-type="bibr" rid="B23">Li and Durbin 2009</xref>). The alignment output in SAM format was converted to BAM, sorted, and indexed using SAMtools v.1.21 (<xref ref-type="bibr" rid="B10">Danecek et al. 2021</xref>).</p>
        <p>Variant calling and genotyping were performed using GATK v.4.2.0 (<xref ref-type="bibr" rid="B31">McKenna et al. 2010</xref>). Genotypes were called in diploid mode, and both homozygous and heterozygous genotypes were retained in the resulting VCF. Because sequencing depth varied among loci, particularly across taxa represented by few individuals, variants were filtered to retain sites with a minimum genotype depth (<abbrev xlink:title="genotype depth">DP</abbrev>) of 3 and minimum genotype quality of 50%. This threshold was selected to balance retention of informative heterozygous sites with locus recovery in a low-coverage <abbrev xlink:title="Genotyping-by-sequencing">GBS</abbrev> dataset, and to maximize the number of loci available for allele-sharing analyses.</p>
        <p>The resulting VCF file was filtered to retain variants with minimum genotype quality scores of 50% and a sequencing depth of at least 3 using VCFR (<xref ref-type="bibr" rid="B19">Knaus and Grunwald 2017</xref>) and SNPfiltR (<xref ref-type="bibr" rid="B11">DeRaad 2022</xref>). While higher depth thresholds can increase confidence in heterozygous genotype calls, exploratory analyses with stricter filtering (<abbrev xlink:title="genotype depth">DP</abbrev> ≥ 5) substantially reduced locus retention without qualitatively altering patterns of allele sharing; therefore, the more inclusive threshold was retained for downstream analyses.</p>
      </sec>
      <sec sec-type="Phylogenetic structure of genetic diversity" id="sec7">
        <title>Phylogenetic structure of genetic diversity</title>
        <p>Phylogenetic structuring was used to quantify admixture and assess the influences of hybridization and gene flow between <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> and sympatric species. The VCF file generated from genomic sequencing was converted into phylip format for maximum likelihood phylogenetic inference using vcf2phylip v.2.8 (<xref ref-type="bibr" rid="B37">Ortiz 2019</xref>). The resulting phylip alignment contained 165,325 SNPs after filtering across 11 taxa, 64.37% of which are invariant sites. A maximum likelihood (<abbrev xlink:title="maximum likelihood">ML</abbrev>) phylogeny was generated for this dataset with RAxML-NG (<xref ref-type="bibr" rid="B21">Kozlov et al. 2019</xref>) using a GTR+G model to account for both substitution rate heterogeneity across sites and unequal base frequencies, 50 starting trees to ensure broad exploration of tree space and avoid suboptimal local maxima, and 200 bootstrap replicates to provide robust support values for inferred relationships and assess the stability of the phylogenetic topology (<xref ref-type="bibr" rid="B39">Pattengale et al. 2010</xref>). <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Pinellia">Pinellia</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="tripartita">tripartita</tp:taxon-name-part></tp:taxon-name></italic> Schott (<tp:taxon-name><tp:taxon-name-part taxon-name-part-type="family" reg="Araceae">Araceae</tp:taxon-name-part></tp:taxon-name>) was included as an outgroup taxon for genomic library construction and downstream phylogenetic analyses; leaf material was collected concurrently with <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part></tp:taxon-name></italic> samples on Shikoku Island and prepared using identical silica-drying and DNA extraction protocols, and the species was selected based on its well-supported placement outside <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part></tp:taxon-name></italic> within <tp:taxon-name><tp:taxon-name-part taxon-name-part-type="family" reg="Araceae">Araceae</tp:taxon-name-part></tp:taxon-name> (<xref ref-type="bibr" rid="B35">Murata et al. 2018</xref>). Ascertainment bias correction was not applied, as invariant sites were retained in the alignment and such correction is deemed unnecessary when the data include both variant and invariant positions (<xref ref-type="bibr" rid="B21">Kozlov et al. 2019</xref>).</p>
        <p>To assess patterns of introgression and excess allele sharing among taxa, we used Dsuite (<xref ref-type="bibr" rid="B28">Malinsky et al. 2021</xref>), which implements Patterson’s D-statistics and related f-statistics directly from VCF files. Using the maximum likelihood phylogeny as a guide topology, we first calculated D-statistics for all possible taxon quartets using the Dtrios module, which tests for asymmetries in allele sharing inconsistent with the inferred species tree. These statistics were then summarized using the f-branch (<abbrev xlink:title="f-branch">fb</abbrev>) statistic, which estimates the proportion of the genome exhibiting excess allele sharing between pairs of lineages and assigns this signal to specific internal branches of the phylogeny. This approach allows introgression to be visualized across multiple species simultaneously while accounting for shared ancestry and phylogenetic structure. Because of the small number of individuals per taxon, these analyses are exploratory and interpreted as indicative of potential admixture rather than definitive evidence of ongoing hybridization.</p>
      </sec>
    </sec>
    <sec sec-type="Results" id="sec8">
      <title>Results</title>
      <sec sec-type="Sex expression in relation to plant size and light availability" id="sec9">
        <title>Sex expression in relation to plant size and light availability</title>
        <p>Logistic regression analyses comparing plant sex (male vs female) and two size variables (plant height and pseudostem diameter) for each species demonstrates that both plant height (Fig. <xref ref-type="fig" rid="F1">1</xref>, column 1) and pseudostem diameter (Fig. <xref ref-type="fig" rid="F1">1</xref>, column 2) are significant predictors of sex expression in four of the five sampled species, with <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> being the only exception (Fig. <xref ref-type="fig" rid="F1">1A</xref>). In <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ehimense">ehimense</tp:taxon-name-part></tp:taxon-name></italic>, <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="tosaense">tosaense</tp:taxon-name-part></tp:taxon-name></italic>, <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ringens">ringens</tp:taxon-name-part></tp:taxon-name></italic>, and <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="iyoanum">iyoanum</tp:taxon-name-part></tp:taxon-name></italic>, plant size is shown to be a strong predictor of sex expression (Fig. <xref ref-type="fig" rid="F1">1B–E</xref>), with shorter individuals bearing a thin pseudostem producing male flowers (left/top of graphs) and taller individuals exhibiting wider a pseudostem more likely to produce female flowers (bottom/right of graphs). The logistic regression between the individuals shows a steep curve (P &lt; 0.05), indicating a significant correlation between sex (male vs female) and plant size (columns 1, 2).</p>
        <fig id="F1">
          <object-id content-type="doi">10.5091/plecevo.177783.figure1</object-id>
          <object-id content-type="arpha">C537E0D5-AAF8-56F0-965D-2C1202B413FB</object-id>
          <label>Figure 1.</label>
          <caption>
            <p>Pairwise logistic regression between plant sex and size/light intensity variables. N.S., *P, **P, and ***P indicate statistical significance at the levels of P &gt; 0.5, P &lt; 0.05, P &lt; 0.005, and P &gt; 0005, respectively. <bold>A</bold>. <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic>. <bold>B</bold>. <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ehimense">ehimense</tp:taxon-name-part></tp:taxon-name></italic>. <bold>C</bold>. <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="tosaense">tosaense</tp:taxon-name-part></tp:taxon-name></italic>. <bold>D</bold>. <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ringens">ringens</tp:taxon-name-part></tp:taxon-name></italic>. <bold>E</bold>. <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="iyoanum">iyoanum</tp:taxon-name-part></tp:taxon-name></italic>. Open circles in each plot indicate individuals, either Male (top) or Female (bottom), recorded at each height in cm (column 1), pseudostem diameter in cm (column 2), or light intensity in klx (column 3). The red line indicates the logistic regression curve with the pitch indicating the strength of differentiation between measurements pertaining to individuals producing either male or female inflorescences.</p>
          </caption>
          <graphic xlink:href="plecevo-159-197-g001.jpg" id="oo_1586660.jpg">
            <uri content-type="original_file">https://binary.pensoft.net/fig/1586660</uri>
          </graphic>
        </fig>
        <p>In <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic>, by contrast, there is no detected significant correlation between sex expression and plant size (Fig. <xref ref-type="fig" rid="F1">1A</xref>, columns 1, 2). Instead, logistic regression analyses for this species indicate that light intensity rather than plant size is the primary determinant of sex expression (Fig. <xref ref-type="fig" rid="F1">1A</xref>, column 3). The logistic regression reveals that in <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic>, individuals exposed to higher light intensity are significantly more likely to express female reproductive structures (P &lt; 0.05), whereas no significant correlation is found between sex and plant height (column 1) or pseudostem diameter (column 2). Correlation tests indicated that light availability was not strongly correlated with plant height (r = 0.11, p = 0.28) or pseudostem diameter (r = 0.09, p = 0.34), confirming that the observed relationship between light and sex was not simply an indirect effect of larger plants occupying brighter sites.</p>
      </sec>
      <sec sec-type="Phylogenetic structure of genetic diversity" id="sec10">
        <title>Phylogenetic structure of genetic diversity</title>
        <p>The maximum likelihood (<abbrev xlink:title="maximum likelihood">ML</abbrev>) phylogenetic tree constructed from genome-wide SNP data allows us to assess overall relationships among sampled individuals while also detecting levels of admixture that represent hybridization and/or ongoing gene flow. In our reconstruction, <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> does not form a clade, instead appearing as a grade with respect to other sampled species (Fig. <xref ref-type="fig" rid="F2">2B</xref>; Suppl. material <xref ref-type="supplementary-material" rid="S1">1</xref>). Bootstrap support values are high (&gt; 85%) for deeper nodes in the tree but lower (&lt; 50%) for relationships among closely related species, particularly among <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic>, <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ehimense">ehimense</tp:taxon-name-part></tp:taxon-name></italic>, and <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="tosaense">tosaense</tp:taxon-name-part></tp:taxon-name></italic>. The analysis also indicates that <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ehimense">ehimense</tp:taxon-name-part></tp:taxon-name></italic> is polyphyletic, with its two defined morphotypes (red-appendix and green-appendix) appearing separately in the phylogeny. The red form is recovered as sister to <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="tosaense">tosaense</tp:taxon-name-part></tp:taxon-name></italic> while the green form is sister to the large clade formed by the <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> grade plus the <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="tosaense">tosaense</tp:taxon-name-part></tp:taxon-name>–<tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ehimense">ehimense</tp:taxon-name-part></tp:taxon-name></italic> red form sister pair. Neither of these relationships for <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ehimense">ehimense</tp:taxon-name-part></tp:taxon-name></italic> is well supported (bootstrap &lt; 55), indicating either incomplete lineage sorting or ongoing hybridization among the populations represented by these individuals.</p>
        <fig id="F2">
          <object-id content-type="doi">10.5091/plecevo.177783.figure2</object-id>
          <object-id content-type="arpha">7C21AF12-269B-58FB-B2C4-E79A5C1F6B26</object-id>
          <label>Figure 2.</label>
          <caption>
            <p>Assessment of admixture in <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part></tp:taxon-name></italic> species on Shikoku Island, Japan. <bold>A</bold>. Sample locations on map of Shikoku Island. <bold>B</bold>. Maximum likelihood phylogeny of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part></tp:taxon-name></italic> individuals. Support values represented by bootstrap values. Branch lengths are proportional to the tree scale. Outgroup (not shown): <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Pinellia">Pinellia</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="tripartita">tripartita</tp:taxon-name-part></tp:taxon-name></italic>. Localities of collections indicated following the species name in parentheses with letters corresponding to abbreviations indicated on the map. “GF” after <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ehimense">ehimense</tp:taxon-name-part></tp:taxon-name></italic> refers to the green appendix morphotype. <bold>C</bold>. Results from the f-branch analysis of excessive allele sharing. Grey dotted lines indicate internal (ancestral) branches of the phylogeny to which introgression signals are assigned in the f-branch framework. These ancestral branches are shown only to aid interpretation of the f-branch analysis and do not represent additional taxa or alternative phylogenetic relationships beyond those depicted in panel B. Shaded cells represent <abbrev xlink:title="f-branch">fb</abbrev> values, with darker red indicating a higher proportion of the genome exhibiting excess allele sharing between a focal lineage and a given internal branch. Cells outlined in green highlight allele sharing involving <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic>.</p>
          </caption>
          <graphic xlink:href="plecevo-159-197-g002.jpg" id="oo_1586658.jpg">
            <uri content-type="original_file">https://binary.pensoft.net/fig/1586658</uri>
          </graphic>
        </fig>
        <p>Exploratory f-branch (<abbrev xlink:title="f-branch">fb</abbrev>) analyses were conducted to identify putative instances of gene flow among sampled individuals and to assess whether admixture patterns correspond with the observed phylogenetic relationships. The strength of admixture between pairs of taxa is visualized using shaded boxes (Fig. <xref ref-type="fig" rid="F2">2C</xref>), with darker hues indicating higher levels of excess allele sharing between species pairs. The dotted lines correspond to inferred ancestral branches, helping to contextualize observed admixture patterns within the context of the broader phylogeny. Our results reveal elevated allele sharing (<abbrev xlink:title="f-branch">fb</abbrev> &gt; 0.025) between <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> and <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ternatipartitum">ternatipartitum</tp:taxon-name-part></tp:taxon-name></italic>, <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ringens">ringens</tp:taxon-name-part></tp:taxon-name></italic>, <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ishizuchiense">ishizuchiense</tp:taxon-name-part></tp:taxon-name></italic>, <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="sikokianum">sikokianum</tp:taxon-name-part></tp:taxon-name></italic>, <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="iyoanum">iyoanum</tp:taxon-name-part></tp:taxon-name></italic>, <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ehimense">ehimense</tp:taxon-name-part></tp:taxon-name></italic>, and <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="tosaense">tosaense</tp:taxon-name-part></tp:taxon-name></italic>, suggesting a history of introgression within this regional complex (Fig. <xref ref-type="fig" rid="F2">2C</xref>). This suggests that hybridization is likely not limited to a single species pair but instead may be widespread across the sampled taxa, with <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> acting as a central contributor to gene flow. Although individual <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> samples showed qualitative differences in allele-sharing signals, all sampled populations occurred in similar plantation or secondary-forest habitats and exhibited broadly comparable sex-expression patterns, with no obvious population-level differences detected. Notably, <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="tosaense">tosaense</tp:taxon-name-part></tp:taxon-name></italic> also exhibits strong signals of excess allele sharing with multiple species, further suggesting that gene flow may be a recurrent process within this lineage and that there is widespread historical and/or ongoing introgression among taxa on Shikoku Island. Both morphotypes of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ehimense">ehimense</tp:taxon-name-part></tp:taxon-name></italic> show strong evidence of admixture with <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> (Fig. <xref ref-type="fig" rid="F2">2C</xref>), despite their polyphyletic position in the phylogeny.</p>
      </sec>
    </sec>
    <sec sec-type="Discussion" id="sec11">
      <title>Discussion</title>
      <p>Our results demonstrate that <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> exhibits a relaxed dependence on plant size for determining sex expression, with light availability emerging as a more consistent predictor of reproductive phase among western Japanese populations. This contrasts sharply with the canonical size–sex relationship documented in other <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part></tp:taxon-name></italic> species, including the four congeners analysed here, and challenges the long-held assumption that sex expression in this genus is strictly allometric (<xref ref-type="bibr" rid="B1">Atkinson 1898</xref>; <xref ref-type="bibr" rid="B5">Camp 1932</xref>; <xref ref-type="bibr" rid="B13">Doust and Cavers 1982</xref>; <xref ref-type="bibr" rid="B16">Ewing and Klein 1982</xref>; <xref ref-type="bibr" rid="B48">Takasu 1987</xref>; <xref ref-type="bibr" rid="B9">Clay 1993</xref>). In <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic>, the probability of producing pistillate inflorescences increases under brighter microsite conditions, regardless of individual size, implying that light availability may act as an environmental cue overriding the intrinsic thresholds governing sex expression (<xref ref-type="bibr" rid="B36">Obeso 2002</xref>).</p>
      <p>While similar light–sex relationships have been observed in <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="triphyllum">triphyllum</tp:taxon-name-part></tp:taxon-name></italic> populations in high-light habitats of Ontario (<xref ref-type="bibr" rid="B13">Doust and Cavers 1982</xref>), those systems showed a strong correlation between light and plant size, making the causal driver ambiguous. In contrast, our analyses revealed no significant correlation between light intensity and either plant height or pseudostem diameter in <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic>, suggesting that growth traits are not the proximate mediators of the light effect. This decoupling may reflect a shift toward direct physiological sensitivity to light in determining reproductive investment, a response that could confer adaptive advantages in patchy or periodically illuminated forest floors. For an understory geophyte such as <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic>, this plasticity allows individuals to opportunistically reproduce when canopy openings increase resource availability, rather than waiting to achieve a fixed vegetative size threshold.</p>
      <p>Our findings contradict early foundational work by Maekawa (<xref ref-type="bibr" rid="B25">1924</xref>, <xref ref-type="bibr" rid="B26">1927</xref>), who reported that individuals of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="japonicum">japonicum</tp:taxon-name-part></tp:taxon-name></italic> (now synonymized under <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> (<xref ref-type="bibr" rid="B47">Scholten et al. 2025</xref>)) generally transitioned from male to female with increasing size and age. However, Maekawa’s studies were conducted under homogeneous greenhouse conditions and did not assess the influence of light on reproductive expression. In fact, Maekawa himself acknowledged that energy assimilation and nutrient status likely played a role in sex transitions, hinting at complexity beyond mere size thresholds.</p>
      <p>The departure from strict size dependence in <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> also provides an explanation for regional variation in reproductive strategies across Japan. <xref ref-type="bibr" rid="B18">Kinoshita (1987)</xref> found a clear size–sex correlation in <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> populations from Nagano Prefecture, where deciduous broadleaf forests dominate. In contrast, the evergreen conifer plantations and secondary mixed forests of Shikoku Island—characterized by reduced understory light and irregular canopy gaps (<xref ref-type="bibr" rid="B38">Ovington 1983</xref>)—present a markedly different selective environment. In such low-light systems, a reproductive strategy responsive to transient illumination may increase both individual fitness and population persistence. These findings suggest that forest composition and management history can exert a strong influence over the evolution of reproductive plasticity in perennial understory plants.</p>
      <p>The ecological consequences of such flexibility extend beyond individual fitness. Light-mediated sex expression likely facilitates colonization of plantation landscapes, where disturbance and canopy thinning create transient windows of high light. Increased reproductive output in these microsites can elevate local population density and the probability of contact among species. The co-occurrence of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> with congeners such as <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ehimense">ehimense</tp:taxon-name-part></tp:taxon-name></italic>, <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="tosaense">tosaense</tp:taxon-name-part></tp:taxon-name></italic>, and <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ringens">ringens</tp:taxon-name-part></tp:taxon-name></italic> within the same plantation mosaics provides numerous opportunities for heterospecific pollen transfer and hybridization. Our genomic analyses—although exploratory due to limited sampling—revealed extensive allele sharing among <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> and sympatric species, consistent with recurrent introgression. These results align with previous reports of hybrid individuals and morphologically intermediate populations throughout western Japan (<xref ref-type="bibr" rid="B33">Murata and Ohno 1989</xref>; <xref ref-type="bibr" rid="B34">Murata and Ohashi 2009</xref>) and indicate that anthropogenic habitats may serve as arenas for both ecological convergence and genetic exchange.</p>
      <p>This pattern of admixture supports a causal sequence in which light-driven reproductive plasticity promotes colonization of plantation forests, thereby increasing opportunities for hybridization. Specifically, our demographic analyses show that sex expression in <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> is decoupled from plant size and instead strongly associated with light availability (Fig. <xref ref-type="fig" rid="F1">1A</xref>), a pattern not observed in sympatric congeners, indicating that ecological flexibility is a species-specific trait rather than a consequence of hybridization. In parallel, phylogenetic and f-branch analyses reveal widespread allele sharing between <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> and multiple sympatric taxa (Fig. <xref ref-type="fig" rid="F2">2B–C</xref>), consistent with recurrent gene flow occurring in the same disturbed, light-variable habitats where relaxed sex expression is observed. Together, these results favour an interpretation in which ecological flexibility enables range expansion and increased interspecific contact, after which gene flow among sympatric taxa reinforces variation and potentially facilitates adaptation, rather than hybridization acting as the primary driver of sex-expression plasticity. The frequent occurrence of morphologically intermediate individuals in plantation understories further supports this feedback, in which environmental heterogeneity, reproductive plasticity, and introgression operate jointly to maintain demographic and genetic diversity within the <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> complex (Fig. <xref ref-type="fig" rid="F3">3</xref>). While our data do not allow us to test whether introgression quantitatively modifies sex-expression strategies, the concurrence of relaxed size dependence and widespread allele sharing suggests that reproductive plasticity may facilitate opportunities for gene flow rather than result from it.</p>
      <fig id="F3">
        <object-id content-type="doi">10.5091/plecevo.177783.figure3</object-id>
        <object-id content-type="arpha">485B2A38-6C74-542D-8E76-A17565CBB833</object-id>
        <label>Figure 3.</label>
        <caption>
          <p>Conceptual model illustrating the ecological and evolutionary sequence linking light availability, reproductive plasticity, and hybridization in <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic>. Under a closed canopy, populations of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> exhibit primarily staminate (male) expression, and interspecific contact with congeners is limited (Box 1). Following canopy disturbance (e.g. through treefall) (Box 2), increased light availability promotes a shift toward pistillate (female) expression (Box 3), enhancing local reproductive output and seed production. As population density increases in canopy gaps and forest edges (Box 4), spatial and temporal overlap with congeners expands, elevating opportunities for heterospecific pollen transfer (Box 5). Gene flow between <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> and sympatric species leads to introgression, generating morphological intermediates and increased genetic diversity (Box 6). Orange arrows represent gene flow between individuals. Staminate and pistillate individuals are noted by (♂) and (♀), respectively.</p>
        </caption>
        <graphic xlink:href="plecevo-159-197-g003.jpg" id="oo_1586659.jpg">
          <uri content-type="original_file">https://binary.pensoft.net/fig/1586659</uri>
        </graphic>
      </fig>
      <p>An intriguing question emerging from these findings is why sympatric species such as <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ehimense">ehimense</tp:taxon-name-part></tp:taxon-name></italic>, <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="tosaense">tosaense</tp:taxon-name-part></tp:taxon-name></italic>, and <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ringens">ringens</tp:taxon-name-part></tp:taxon-name></italic> retain strict size-dependent sex determination despite occupying similar habitats. One possibility is that their divergent phenology or differing energy reserves limit responsiveness to transient light changes, constraining the expression of plasticity. Alternatively, underlying hormonal or gene-regulatory mechanisms may differ among taxa, with <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> having a greater physiological sensitivity to photic cues. Future work should experimentally test these hypotheses through controlled shading and light-manipulation experiments, coupled with hormone assays to determine whether endogenous regulators such as gibberellins or cytokinins mediate the observed responses. Comparative transcriptomic or methylation studies could also identify candidate pathways underlying light responsiveness and reveal whether these mechanisms are shared or species-specific.</p>
      <p>The genomic results presented here, though limited in resolution, underscore the evolutionary fluidity within this species complex. The paraphyly of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> relative to <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ehimense">ehimense</tp:taxon-name-part></tp:taxon-name></italic> and <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="tosaense">tosaense</tp:taxon-name-part></tp:taxon-name></italic>, combined with widespread allele sharing, points to weak reproductive barriers and frequent historical gene flow. While the small number of individuals analysed prevents robust population-genetic inference, the congruence of morphological and molecular patterns strongly suggests that introgression is ongoing. Broader sampling across the geographic range of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic>—including both plantation and native-forest populations—will be essential to test whether hybridization rates are indeed higher in managed habitats and to disentangle incomplete lineage sorting from contemporary gene exchange. Such work could also clarify whether adaptive alleles associated with light sensitivity have experienced introgression across species boundaries, potentially linking hybridization to the evolution of such plastic sexual expression.</p>
      <p>The unexpected polyphyly of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ehimense">ehimense</tp:taxon-name-part></tp:taxon-name></italic> in the phylogeny is particularly significant, as this study represents the first phylogenetic resolution of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part></tp:taxon-name></italic> species on Shikoku Island. The original description of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ehimense">ehimense</tp:taxon-name-part></tp:taxon-name></italic> by <xref ref-type="bibr" rid="B33">Murata and Ohno (1989)</xref> identifies two distinct morphotypes, one with a red appendix and another with a green appendix. Our analysis indicates that individuals with a green appendix do not form a monophyletic group with the more widely recognized red-appendix form. This suggests that the green-appendix morphotype, while superficially matching the morphological criteria of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ehimense">ehimense</tp:taxon-name-part></tp:taxon-name></italic>, may represent a hybrid between <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> and another species outside the <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic>–<italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="tosaense">tosaense</tp:taxon-name-part></tp:taxon-name></italic> hybrid complex. The f-branch analysis (Fig. <xref ref-type="fig" rid="F2">2C</xref>) supports this hypothesis, indicating one-way allele sharing with <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> and no other sampled species, suggesting that the second parental species was not included in our analysis. Further genetic sampling across Japan, targeting more individuals and taxa not included in this study, could help provide insight into the taxonomy of this species.</p>
      <p>The increasing number of documented natural hybrids between <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> and its congeners (<xref ref-type="bibr" rid="B33">Murata and Ohno 1989</xref>; <xref ref-type="bibr" rid="B34">Murata and Ohashi 2009</xref>) further supports the idea that hybridization and introgression significantly influence the genetic structure of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part></tp:taxon-name></italic> populations on Shikoku Island. Collectively, our findings suggest that the demographic success of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> across western Japan arises from the interaction of environmental heterogeneity, reproductive plasticity, and gene flow. By decoupling the reproductive phase from plant size and coupling it instead to light availability, individuals of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> can exploit brief periods of high irradiance within plantation forests, ensuring reproductive output under otherwise limiting conditions. This ecological strategy promotes persistence and expansion in anthropogenically altered habitats, while also facilitating contact and hybridization with congeners. The resulting feedback between ecological opportunity and genetic exchange likely contributes to the morphological diversity and phylogenetic complexity characteristic of this species. Future studies integrating demographic monitoring, pollinator observations, and expanded genomic datasets will be crucial for quantifying the fitness consequences of this flexibility and for understanding how human-modified environments continue to shape the evolutionary trajectories of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part></tp:taxon-name></italic> and other clonal, sexually plastic forest herbs.</p>
    </sec>
    <sec sec-type="Conclusions" id="sec12">
      <title>Conclusions</title>
      <p>Taken together, our results demonstrate that light availability, rather than plant size, is the dominant predictor of sex expression in <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> populations on Shikoku Island. This shift from a size- to a light-dependent reproductive strategy likely reflects adaptation to the shaded, spatially heterogeneous conditions characteristic of evergreen conifer plantations. The decoupling of size and light effects suggests that reproductive transitions in <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> are triggered directly by environmental cues rather than by accumulated biomass, allowing individuals to reproduce under the variable light regimes typical of managed forests. These ecological patterns coincide with genomic evidence of allele sharing between <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="serratum">serratum</tp:taxon-name-part></tp:taxon-name></italic> and several sympatric congeners, consistent with ongoing or historical introgression within plantation and secondary-forest mosaics. We therefore interpret reproductive plasticity as the initial driver of colonization in these modified habitats, with hybridization emerging as a secondary outcome of increased spatial overlap and flowering synchrony among species. Future work combining broader population sampling, experimental manipulations of light intensity, and higher-resolution genomic data will be necessary to test whether the observed plasticity is heritable, to quantify the direction and extent of gene flow, and to clarify how anthropogenic forest structure continues to influence reproductive and evolutionary dynamics in <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part></tp:taxon-name></italic>.</p>
    </sec>
    <sec sec-type="Data availability" id="sec13">
      <title>Data availability</title>
      <p>All sample sequence data, the multisample VCF, and the logistic regression matrix are openly available for download from Zenodo (<ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.5281/zenodo.18163207">https://doi.org/10.5281/zenodo.18163207</ext-link>). Additional supporting information may be found online in Suppl. material <xref ref-type="supplementary-material" rid="S1">1</xref>.</p>
    </sec>
  </body>
  <back>
    <ack>
      <title>Acknowledgments</title>
      <p>The authors thank the University of Wisconsin Biotechnology Center DNA Sequencing Facility (Research Resource Identifier – RRID: SCR_017759) for providing sequencing services, Jacob B. Landis for assistance in genome data processing, and Shun Tanaka and Olivia Hullihen for field assistance during sample collection. This work was supported by funds provided by the Lewis and Clark Fund for Exploration and Research, the Torrey Botanical Society, the Society of Systematic Biologists, the Mario Einaudi Center of International Studies, and the Moore Fund of the L.H. Bailey Hortorium Herbarium. Additionally, the stipend of the first author was supported in part by the Tara Atluri Memorial Fund.</p>
    </ack>
    <ref-list>
      <title>References</title>
      <ref id="B1">
        <mixed-citation>Atkinson G (1898) Experiments on the morphology of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus">Arisaema</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species">triphyllum</tp:taxon-name-part></tp:taxon-name></italic>. Botanical Gazette 25: 114.</mixed-citation>
      </ref>
      <ref id="B2">
        <mixed-citation>Barot S, Mitja D., Miranda I, Meija G, Grimaldi M (2005) Reproductive plasticity in an Amazonian palm. Evolutionary Ecology Research 7: 1051 –1075.</mixed-citation>
      </ref>
      <ref id="B3">
        <mixed-citation>Beckman Coulter (2012) SPRI-based size selection. Brea, California, 1–14.</mixed-citation>
      </ref>
      <ref id="B4">
        <mixed-citation>Bierzychudek P (1982) The demography of Jack-in-the-Pulpit, a forest perennial that changes sex. Ecological Monographs 82: 335–351. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.2307/2937350">https://doi.org/10.2307/2937350</ext-link></mixed-citation>
      </ref>
      <ref id="B5">
        <mixed-citation>Camp W (1932) Sex in <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus">Arisaema</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species">triphyllum</tp:taxon-name-part></tp:taxon-name></italic>. Ohio Journal of Science 32: 147–151.</mixed-citation>
      </ref>
      <ref id="B6">
        <mixed-citation>Catchen J, Hohenlohe PA, Bassham S, Amores A, Cresko WA (2013) Stacks: an analysis tool set for population genomics. Molecular Ecology 22: 3124–3140. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1111/mec.12354">https://doi.org/10.1111/mec.12354</ext-link></mixed-citation>
      </ref>
      <ref id="B7">
        <mixed-citation>Charnov EL, Bull J (1977) When is sex environmentally determined? Nature 266: 828–830. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1038/266828a0">https://doi.org/10.1038/266828a0</ext-link></mixed-citation>
      </ref>
      <ref id="B8">
        <mixed-citation>Chen S (2023) Ultrafast one-pass FASTQ data preprocessing, quality control, and deduplication using fastp. iMeta 2: e107. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1002/imt2.107">https://doi.org/10.1002/imt2.107</ext-link></mixed-citation>
      </ref>
      <ref id="B9">
        <mixed-citation>Clay K (1993) Size-dependent gender change in green dragon (<italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus">Arisaema</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species">dracontium</tp:taxon-name-part></tp:taxon-name></italic>; <tp:taxon-name><tp:taxon-name-part taxon-name-part-type="family">Araceae</tp:taxon-name-part></tp:taxon-name>). American Journal of Botany 80: 769–777. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1002/j.1537-2197.1993.tb15293.x">https://doi.org/10.1002/j.1537-2197.1993.tb15293.x</ext-link></mixed-citation>
      </ref>
      <ref id="B10">
        <mixed-citation>Danecek P, Bonfield J, Liddle J, Marshall J, Ohan V, Pollard M, Whitwham A, Keane T, McCarthy SA, Davies RM, Li H (2021) Twelve years of SAMtools and BCFtools. GigaScience 10: giab008. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1093/gigascience/giab008">https://doi.org/10.1093/gigascience/giab008</ext-link></mixed-citation>
      </ref>
      <ref id="B11">
        <mixed-citation>DeRaad DA (2022) SNPfiltR: an R package for interactive and reproducible SNP filtering. Molecular Ecology Resources 22: 2443–2453. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1111/1755-0998.13618">https://doi.org/10.1111/1755-0998.13618</ext-link></mixed-citation>
      </ref>
      <ref id="B12">
        <mixed-citation>Dorken ME, Barrett SCH (2004) Phenotypic plasticity of vegetative and reproductive traits in monoecious and dioecious populations of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus">Sagittaria</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species">latifolia</tp:taxon-name-part></tp:taxon-name></italic> (<tp:taxon-name><tp:taxon-name-part taxon-name-part-type="family">Alismataceae</tp:taxon-name-part></tp:taxon-name>): a clonal aquatic plant. Journal of Ecology 92: 32–44. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1111/j.1365-2745.2004.00857.x">https://doi.org/10.1111/j.1365-2745.2004.00857.x</ext-link></mixed-citation>
      </ref>
      <ref id="B13">
        <mixed-citation>Doust JL, Cavers PB (1982) Sex and gender dynamics in Jack-in-the-Pulpit, <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus">Arisaema</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species">triphyllum</tp:taxon-name-part></tp:taxon-name></italic> (<tp:taxon-name><tp:taxon-name-part taxon-name-part-type="family">Araceae</tp:taxon-name-part></tp:taxon-name>). Ecology 63: 797–808. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.2307/1936800">https://doi.org/10.2307/1936800</ext-link></mixed-citation>
      </ref>
      <ref id="B14">
        <mixed-citation>Doyle JJ, Doyle JL (1987) A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochemical Bulletin 19: 11–15.</mixed-citation>
      </ref>
      <ref id="B15">
        <mixed-citation>Ellstrand NC, Schierenbeck KA (2000) Hybridization as a stimulus for the evolution of invasiveness in plants? PNAS 97: 7043–7050. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1073/pnas.97.13.7043">https://doi.org/10.1073/pnas.97.13.7043</ext-link></mixed-citation>
      </ref>
      <ref id="B16">
        <mixed-citation>Ewing JW, Klein RM (1982) Sex expression in Jack-in-the-Pulpit. Bulletin of the Torrey Botanical Club 109: 47–50. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.2307/2484467">https://doi.org/10.2307/2484467</ext-link></mixed-citation>
      </ref>
      <ref id="B17">
        <mixed-citation>Gow JE (1913) Observations on the morphology of the Aroids. Botanical Gazette 56: 127–142. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1086/331128">https://doi.org/10.1086/331128</ext-link></mixed-citation>
      </ref>
      <ref id="B18">
        <mixed-citation>Kinoshita E (1987) Sex change and population dynamics in <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus">Arisaema</tp:taxon-name-part></tp:taxon-name></italic> (<tp:taxon-name><tp:taxon-name-part taxon-name-part-type="family">Araceae</tp:taxon-name-part></tp:taxon-name>) I. <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus">Arisaema</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species">serratum</tp:taxon-name-part></tp:taxon-name></italic> (Thunb.) Schott. Plant Species Biology 2: 12–28. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1111/j.1442-1984.1987.tb00030.x">https://doi.org/10.1111/j.1442-1984.1987.tb00030.x</ext-link></mixed-citation>
      </ref>
      <ref id="B19">
        <mixed-citation>Knaus BJ, Grunwald NJ (2017) VCFR: a package to manipulate and visualize variant call format data in R. Molecular Ecology Resources 17: 44–53. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1111/1755-0998.12549">https://doi.org/10.1111/1755-0998.12549</ext-link></mixed-citation>
      </ref>
      <ref id="B20">
        <mixed-citation>Korneliussen TS, Albrechtsen A, Nielsen R (2014) ANGSD: analysis of next generation sequencing data. BMC Bioinformatics 15: 356. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1186/s12859-014-0356-4">https://doi.org/10.1186/s12859-014-0356-4</ext-link></mixed-citation>
      </ref>
      <ref id="B21">
        <mixed-citation>Kozlov AM, Darriba D, Flouri T, Morel B, Stamatakis A (2019) RAxML-NG: a fast, scalable, and user-friendly tool for maximum likelihood phylogenetic inference. Bioinformatics 35: 4453–4455. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1093/bioinformatics/btz305">https://doi.org/10.1093/bioinformatics/btz305</ext-link></mixed-citation>
      </ref>
      <ref id="B22">
        <mixed-citation>Kudo G, Ida TY, Tani T (2008) Linkages between phenology, pollination, photosynthesis, and reproduction in deciduous forest understory plants. Ecology 89: 321–331. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1890/06-2131.1">https://doi.org/10.1890/06-2131.1</ext-link></mixed-citation>
      </ref>
      <ref id="B23">
        <mixed-citation>Li H, Durbin R (2009) Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25: 1754–1760. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1093/bioinformatics/btp324">https://doi.org/10.1093/bioinformatics/btp324</ext-link></mixed-citation>
      </ref>
      <ref id="B24">
        <mixed-citation>Lichstein JW, Peterson BT, Langebreak J, McKinley SA (2021) Leaf economics of early- and late-successional plants. The American Naturalist 198: 347–359. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1086/715453">https://doi.org/10.1086/715453</ext-link></mixed-citation>
      </ref>
      <ref id="B25">
        <mixed-citation>Maekawa T (1924) On the phenomena of sex transition in <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus">Arisaema</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species">japonica</tp:taxon-name-part></tp:taxon-name></italic> Bl. Journal of the College of Agriculture, Hakkaido Imperial University 13: 217–305.</mixed-citation>
      </ref>
      <ref id="B26">
        <mixed-citation>Maekawa T (1927) On intersexualism in <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus">Arisaema</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species">japonica</tp:taxon-name-part></tp:taxon-name></italic> Bl. Japanese Journal of Botany 3: 204–216.</mixed-citation>
      </ref>
      <ref id="B27">
        <mixed-citation>Mahall BE, Bormann FH (1978) A quantitative description of the vegetative phenology of herbs in a northern hardwood forest. Botanical Gazette 139: 467–481. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1086/337022">https://doi.org/10.1086/337022</ext-link></mixed-citation>
      </ref>
      <ref id="B28">
        <mixed-citation>Malinsky M, Matschiner M, Svardal H (2021) Dsuite - Fast D-statistics and related admixture evidence from VCF files. Molecular Ecology Resources 21: 584–595. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1111/1755-0998.13265">https://doi.org/10.1111/1755-0998.13265</ext-link></mixed-citation>
      </ref>
      <ref id="B29">
        <mixed-citation>Masaki T, Oguro M, Yamashita N, Otani T, Utsugi H (2017) Reforestation following harvesting of conifer plantations in Japan: current issues from silvicultural and ecological perspectives. Reforesta 3: 125–142. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.21750/REFOR.3.11.35">https://doi.org/10.21750/REFOR.3.11.35</ext-link></mixed-citation>
      </ref>
      <ref id="B30">
        <mixed-citation>Matsumoto T, Sakuwa S, Murata J (2018) A taxonomic study of the <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus">Arisaema</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species">serratum</tp:taxon-name-part></tp:taxon-name></italic> group (<tp:taxon-name><tp:taxon-name-part taxon-name-part-type="family">Araceae</tp:taxon-name-part></tp:taxon-name>) in northern Okayama Prefecture, western Japan. Journal of Japanese Botany 93: 253–268.</mixed-citation>
      </ref>
      <ref id="B31">
        <mixed-citation>McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Ternytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, DePristo MA (2010) The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Research 20: 1297–1303. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1101/gr.107524.110">https://doi.org/10.1101/gr.107524.110</ext-link></mixed-citation>
      </ref>
      <ref id="B32">
        <mixed-citation>Murata J (2011) The Picture Book of Plant Systematics in Color: <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus">Arisaema</tp:taxon-name-part></tp:taxon-name></italic> in Japan. Hokuryukan Publishing Co., Tokyo.</mixed-citation>
      </ref>
      <ref id="B33">
        <mixed-citation>Murata J, Ohno J (1989) <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus">Arisaema</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species">ehimense</tp:taxon-name-part></tp:taxon-name></italic> J. Murata et Ohno (<tp:taxon-name><tp:taxon-name-part taxon-name-part-type="family">Araceae</tp:taxon-name-part></tp:taxon-name>), a new species from Shikoku, Japan, of putative hybrid origin. Journal of Japanese Botany 64: 341–351.</mixed-citation>
      </ref>
      <ref id="B34">
        <mixed-citation>Murata J, Ohashi H (2009) Taxonomic history of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus">Arisaema</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species">serratum</tp:taxon-name-part></tp:taxon-name></italic> and <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species">japonicum</tp:taxon-name-part></tp:taxon-name></italic>. Bunrui 9: 37–45.</mixed-citation>
      </ref>
      <ref id="B35">
        <mixed-citation>Murata J, Ohno J, Koeayashi T, Ohi-Toma T (2018) The Genus <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus">Arisaema</tp:taxon-name-part></tp:taxon-name></italic> in Japan. Hokuryukan Publishing Co., Tokyo.</mixed-citation>
      </ref>
      <ref id="B36">
        <mixed-citation>Obeso JR (2002) The cost of reproduction in plants. New Phytologist 155: 321–348. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1046/j.1469-8137.2002.00477.x">https://doi.org/10.1046/j.1469-8137.2002.00477.x</ext-link></mixed-citation>
      </ref>
      <ref id="B37">
        <mixed-citation>Ortiz E (2019) vcf2phylip v2.0: convert a VCF matrix into several matrix formats for phylogenetic analysis. (v2.0). Zenodo. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.5281/zenodo.1257057">https://doi.org/10.5281/zenodo.1257057</ext-link></mixed-citation>
      </ref>
      <ref id="B38">
        <mixed-citation>Ovington JD (1983) Temperate Broad-Leaved Evergreen Forests (Ecosystems of the World). Elsevier, New York.</mixed-citation>
      </ref>
      <ref id="B39">
        <mixed-citation>Pattengale ND, Alipour M, Bininda-Emonds ORP, Moret BME, Stamatakis A (2010) How many bootstrap replicates are necessary? Journal of Computational Biology 17: 337–354. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1089/cmb.2009.0179">https://doi.org/10.1089/cmb.2009.0179</ext-link></mixed-citation>
      </ref>
      <ref id="B40">
        <mixed-citation>Policansky D (1981) Sex choice and the size advantage model in jack-in-the-pulpit (<italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus">Arisaema</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species">triphyllum</tp:taxon-name-part></tp:taxon-name></italic>). Population Biology 78: 1306–1308. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1073/pnas.78.2.1306">https://doi.org/10.1073/pnas.78.2.1306</ext-link></mixed-citation>
      </ref>
      <ref id="B41">
        <mixed-citation>Protchazka LS, Alcantara A, Rando JG, Vasconcelos T, Pizzardo RC, Nogueira A (2024) Resource availability and disturbance frequency shape evolution of plant life forms in Neotropical habitats. New Phytologist 242: 760–773. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1111/nph.19601">https://doi.org/10.1111/nph.19601</ext-link></mixed-citation>
      </ref>
      <ref id="B42">
        <mixed-citation>R Core Team (2021) R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna. <ext-link ext-link-type="uri" xlink:href="https://www.r-project.org">https://www.r-project.org</ext-link> [accessed 26.01.2026]</mixed-citation>
      </ref>
      <ref id="B43">
        <mixed-citation>Rieseberg LH, Kim S-C, Randell RA, Whitney KD, Gross BL, Lexer C, Clay K (2007) Hybridization and the colonization of novel habitats by annual sunflowers. Genetica 129: 149–165. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1007/s10709-006-9011-y">https://doi.org/10.1007/s10709-006-9011-y</ext-link></mixed-citation>
      </ref>
      <ref id="B44">
        <mixed-citation>Schaffner JH (1922) Control of the sexual state in <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus">Arisaema</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species">triphyllum</tp:taxon-name-part></tp:taxon-name></italic> and <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus">Arisaema</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species">dracontium</tp:taxon-name-part></tp:taxon-name></italic>. American Journal of Botany 9: 72–78. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.2307/2435278">https://doi.org/10.2307/2435278</ext-link></mixed-citation>
      </ref>
      <ref id="B45">
        <mixed-citation>Scholten J, Specht CD (2025) <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus">Arisaema</tp:taxon-name-part></tp:taxon-name></italic> species from Shikoku Island, Japan [dataset]. Dryad. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.5061/dryad.7pvmcvf4q">https://doi.org/10.5061/dryad.7pvmcvf4q</ext-link></mixed-citation>
      </ref>
      <ref id="B46">
        <mixed-citation>Scholten J, Specht CD (2026) Genomic sequences and ecological data for <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus">Arisaema</tp:taxon-name-part></tp:taxon-name></italic> species on Shikoku Island, Japan [dataset]. Zenodo. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.5281/zenodo.18163207">https://doi.org/10.5281/zenodo.18163207</ext-link></mixed-citation>
      </ref>
      <ref id="B47">
        <mixed-citation>Scholten J, Sprenger A, Hullihen O, Specht CD (2025) A bookend in the <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus">Arisaema</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species">japonicum</tp:taxon-name-part></tp:taxon-name></italic> (<tp:taxon-name><tp:taxon-name-part taxon-name-part-type="family">Araceae</tp:taxon-name-part></tp:taxon-name>) taxonomic debate: morphological and genetic evidence for synonymization. Phytotaxa 751: 1–10. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.11646/phytotaxa.725.1.1">https://doi.org/10.11646/phytotaxa.725.1.1</ext-link></mixed-citation>
      </ref>
      <ref id="B48">
        <mixed-citation>Takasu H (1987) Life history studies on <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus">Arisaema</tp:taxon-name-part></tp:taxon-name></italic> (<tp:taxon-name><tp:taxon-name-part taxon-name-part-type="family">Araceae</tp:taxon-name-part></tp:taxon-name>) I. Growth and reproductive biology of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus">Arisaema</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species">urashima</tp:taxon-name-part></tp:taxon-name></italic> Hara. Plant Species Biology 2: 29–56. <ext-link xlink:type="simple" ext-link-type="doi" xlink:href="10.1111/j.1442-1984.1987.tb00031.x">https://doi.org/10.1111/j.1442-1984.1987.tb00031.x</ext-link></mixed-citation>
      </ref>
    </ref-list>
    <sec sec-type="supplementary-material">
      <title>Supplementary materials</title>
      <supplementary-material id="S1" position="float" orientation="portrait" xlink:type="simple">
        <object-id content-type="doi">10.5091/plecevo.177783.suppl1</object-id>
        <object-id content-type="arpha">276CE52B-3786-5905-9BBE-165B58735AC1</object-id>
        <label>Supplementary material 1</label>
        <statement content-type="notes">
          <p>Maximum likelihood phylogeny of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">Arisaema</tp:taxon-name-part></tp:taxon-name></italic> species. Support values are represented by bootstrap values. Branch lengths are proportional to the tree scale. Outgroup: <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Pinellia">Pinellia</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="tripartita">tripartita</tp:taxon-name-part></tp:taxon-name></italic>. “GF” refers to the green appendix morphotype of <italic><tp:taxon-name><tp:taxon-name-part taxon-name-part-type="genus" reg="Arisaema">A.</tp:taxon-name-part> <tp:taxon-name-part taxon-name-part-type="species" reg="ehimense">ehimense</tp:taxon-name-part></tp:taxon-name></italic>.</p>
        </statement>
        <media xlink:href="plecevo-159-197-s001.pdf" mimetype="application" mime-subtype="pdf" position="float" orientation="portrait" id="oo_1586661.pdf">
          <uri content-type="original_file">https://binary.pensoft.net/file/1586661</uri>
        </media>
      </supplementary-material>
    </sec>
  </back>
</article>
