Research Article |
|
Corresponding author: Justin T. Scholten ( jts329@cornell.edu ) Academic editor: Myriam Heuertz
© 2026 Justin T. Scholten, Chelsea D. Specht.
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.
Citation:
Scholten JT, Specht CD (2026) Does size matter? The integrated roles of light, adaptive sex expression, and hybridization in a widespread Arisaema (Araceae) species from western Japan. Plant Ecology and Evolution 159(2): 197-207. https://doi.org/10.5091/plecevo.177783
|
Background and aims – 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. Arisaema serratum (Araceae), 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 A. serratum 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.
Materials and methods – We conducted field-based demographic surveys across five populations of Arisaema 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.
Key results – 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 A. serratum and sympatric species, indicating ongoing hybridization and introgression.
Conclusions – 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 A. serratum.
Araceae, colonization, diphasic sexual expression, niche expansion, plasticity, Shikoku Island
Species of Arisaema 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 (
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 (
More recent field-based work by
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 (
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 (
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 (
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 Arisaema 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.
In this study, we test how environmental and morphological factors interact to shape reproductive phase transitions in Arisaema serratum, a widespread Japanese species commonly found in post-disturbance evergreen plantations and secondary forests (
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 Arisaema serratum and closely related congeners occurring in managed plantation and secondary-forest habitats. All populations were situated within mixed evergreen conifer plantations dominated by Cryptomeria japonica (Thunb. ex L.f.) D.Don and Chamaecyparis obtusa (Siebold & 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.
At each site, all flowering Arisaema individuals encountered along 100–250 m transects were recorded and identified to species following
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.
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.
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.
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 (
sex ~ height + pseudostem diameter + light intensity
Sex was coded as a binary response variable (male = 0, female = 1). Model performance was assessed using Akaike’s Information Criterion (AIC) and pseudo-R2. 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.
To complement the ecological analyses, we generated exploratory genomic data to assess potential genetic relationships and gene flow among sympatric Arisaema species. Fresh leaf tissue samples were collected from individuals of A. ternatipartitum Makino (n = 1), A. ringens (n = 1), A. ishizuchiense Murata (n = 1), A. sikokianum Franch. & Sav. (n = 1), A. iyoanum (n = 1), A. serratum (n = 3), A. ehimense (n = 2), and A. tosaense (n = 1) and dried in silica. DNA was extracted using a CTAB-based method (
Genotyping-by-sequencing (GBS) 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 (
Sequence data from all individuals were pre-processed using fastp (
Variant calling and genotyping were performed using GATK v.4.2.0 (
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 (
Phylogenetic structuring was used to quantify admixture and assess the influences of hybridization and gene flow between Arisaema serratum 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 (
To assess patterns of introgression and excess allele sharing among taxa, we used Dsuite (
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.
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 > 0.5, P < 0.05, P < 0.005, and P > 0005, respectively. A. Arisaema serratum. B. A. ehimense. C. A. tosaense. D. A. ringens. E. A. iyoanum. 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.
In A. serratum, by contrast, there is no detected significant correlation between sex expression and plant size (Fig.
The maximum likelihood (ML) 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, A. serratum does not form a clade, instead appearing as a grade with respect to other sampled species (Fig.
Assessment of admixture in Arisaema species on Shikoku Island, Japan. A. Sample locations on map of Shikoku Island. B. Maximum likelihood phylogeny of Arisaema individuals. Support values represented by bootstrap values. Branch lengths are proportional to the tree scale. Outgroup (not shown): Pinellia tripartita. Localities of collections indicated following the species name in parentheses with letters corresponding to abbreviations indicated on the map. “GF” after A. ehimense refers to the green appendix morphotype. C. 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 fb 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 A. serratum.
Exploratory f-branch (fb) 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.
Our results demonstrate that Arisaema serratum 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 Arisaema species, including the four congeners analysed here, and challenges the long-held assumption that sex expression in this genus is strictly allometric (
While similar light–sex relationships have been observed in A. triphyllum populations in high-light habitats of Ontario (
Our findings contradict early foundational work by Maekawa (
The departure from strict size dependence in A. serratum also provides an explanation for regional variation in reproductive strategies across Japan.
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 A. serratum with congeners such as A. ehimense, A. tosaense, and A. ringens 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 A. serratum and sympatric species, consistent with recurrent introgression. These results align with previous reports of hybrid individuals and morphologically intermediate populations throughout western Japan (
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 A. serratum is decoupled from plant size and instead strongly associated with light availability (Fig.
Conceptual model illustrating the ecological and evolutionary sequence linking light availability, reproductive plasticity, and hybridization in Arisaema serratum. Under a closed canopy, populations of A. serratum 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 A. serratum 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.
An intriguing question emerging from these findings is why sympatric species such as A. ehimense, A. tosaense, and A. ringens 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 A. serratum 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.
The genomic results presented here, though limited in resolution, underscore the evolutionary fluidity within this species complex. The paraphyly of A. serratum relative to A. ehimense and A. tosaense, 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 A. serratum—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.
The unexpected polyphyly of A. ehimense in the phylogeny is particularly significant, as this study represents the first phylogenetic resolution of Arisaema species on Shikoku Island. The original description of A. ehimense by
The increasing number of documented natural hybrids between A. serratum and its congeners (
Taken together, our results demonstrate that light availability, rather than plant size, is the dominant predictor of sex expression in Arisaema serratum 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 A. serratum 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 A. serratum 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 Arisaema.
All sample sequence data, the multisample VCF, and the logistic regression matrix are openly available for download from Zenodo (https://doi.org/10.5281/zenodo.18163207). Additional supporting information may be found online in Suppl. material
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.
Maximum likelihood phylogeny of Arisaema species. Support values are represented by bootstrap values. Branch lengths are proportional to the tree scale. Outgroup: Pinellia tripartita. “GF” refers to the green appendix morphotype of A. ehimense.