**摘要**
开花、受精和结果的时间组织是被子植物生活史进化的一个基本轴心。虽然大多数物种在一个生长季节内完成果实发育（“一年结果”），但许多壳斗科物种将受精和果实成熟推迟到开花后的第二年（“两年结果”）。尽管这种策略在生态上非常普遍，但其进化起源及其与其他功能性状的协调机制仍不甚明了。在这项研究中，我们通过重建包含88个壳斗科物种的系统发育树来探讨两年结果策略的进化起源，这些物种代表了壳斗科的所有八个属。鉴于两年结果与动物传粉属的关联以及叶习性可能限制繁殖时间的因素，传粉方式和叶习性都可能影响两年结果的进化。为了验证这些可能性，我们进一步采用了系统发育比较分析，以确定结果性状是否受到传粉方式（动物传粉与风传粉）或叶习性（常绿与落叶）的进化限制。我们的结果支持两年结果起源于排除栎属（Fagus）和三角叶栎属（Trigonobalanus）的主要分支的共同祖先，并随后在栗属（Castanea）、栎属（Quercus）和栲属（Castanopsis）等谱系中多次独立恢复为一年结果。对传粉方式和叶习性的祖先状态重建强烈表明，昆虫传粉和常绿叶习性是栎亚科（Quercoideae）的祖先特征，早于两年结果策略的出现。尽管昆虫传粉和常绿习性是栎亚科的祖先特征，但我们的相关性分析表明，结果策略的转变并不依赖于这些特征的变化。这些结果提供了宏观进化证据，表明传粉和受精之间的时间间隔的极端延长在进化上可能是稳定的，但在漫长的进化时间尺度上很少重新出现。
**1 引言**
关键繁殖事件——开花、受精和果实成熟的时间组织——在塑造生物适应性方面起着基础性作用，并构成了被子植物物候和生活史多样化的基础（Satake等人2022；Schwartz 2003）。大多数被子植物在24-48小时内完成从传粉到受精的转变（Williams 2008；Williams和Reese 2019）。然而，有一部分分类单元进化出了受精显著延迟的繁殖策略，从几天到开花后一年以上不等（Sogo和Tobe 2006；Williams 2008；Williams和Reese 2019）。这些策略统称为延迟受精（Sogo和Tobe 2006；Williams 2008；Williams和Reese 2019），它们代表了与典型被子植物繁殖物候的显著偏离，并为探索延长繁殖周期的适应性意义提供了独特的机会（Deng等人2022；Satake和Kelly 2021；Satake等人2023；Schoonderwoerd和Friedman 2016；Shagawa等人2025；Yao等人2023）。壳斗科，包括橡树和山毛榉，就是这种偏离的典型例子。在许多壳斗科物种中，受精被推迟到开花后的第二年，从而导致“两年结果”策略，即受精和随后的果实成熟发生在开花后的那一年。相比之下，在“一年结果”策略下，从开花到受精再到果实成熟的整个繁殖过程都在一个生长季节内完成（图1；Benson 1894；Deng等人2022；Sogo和Tobe 2006）。在整个科中，一年结果和两年结果策略共存，不同属之间的普遍程度差异显著：所有栎属物种都表现为一年结果，而大多数栲属和石栎属物种表现出两年结果，栎属则呈现中间模式（Araye等人2023；Satake和Kelly 2021）。尽管延迟受精的两年结果已经为人所知一个多世纪，但其进化起源、一年结果和两年结果之间的转换频率和方向仍不清楚（Satake和Kelly 2021；Satake等人2023；Shagawa等人2025）。
**2 材料与方法**
2.1 特征数据收集
我们遵循了Zhou等人（2022）的分类单元采样方案，该方案涵盖了八个属的90个壳斗科物种，并提供了一个基因组规模的核基因系统发育树。现存的壳斗科物种大约有900-1000种（Deng等人2022；Kremer等人2012；Satake和Kelly 2021）。尽管这个样本只代表了总物种多样性的一部分，但它涵盖了科内的所有主要谱系，是目前可用于所有属之间比较分析的最全面的数据集。为了与最近的分子系统发育框架保持一致（Hipp等人2020；Liu等人2025；Shen等人2025；Yang等人2025；Zhou等人2022），我们采用了更新的属定义：在中国植物志（eFlora of China）中列为Cyclobalanopsis和Formanodendron的物种分别被归类为栎属（Quercus）和三角叶栎属（Trigonobalanus）（Brach和Song 2006），而在北美植物志（eFlora of North America）中列为Lithocarpus的物种被归类为Notholithocarpus（Nixon 1997）。我们排除了Quercus liaotungensis Koidz.，因为它在中国植物志中被视为Q. mongolica Fisch. ex Ledeb.的同义词（Brach和Song 2006），而后者已经包含在数据集中。我们从中国植物志（Brach和Song 2006）和北美植物志（Nixon 1997）中编译了88个物种的果实特征数据——一年结果或两年结果（表S1）。对于在中国植物志或北美植物志中关于这三个特征的信息不可用或可能不准确的19个物种，我们参考了额外的文献（见表S1）。对于文献中信息缺失的8个物种，我们检查了表S2中列出的在线标本库标本。值得注意的是，Lithocarpus tephrocarpus（Drake）A.Camus根据野外观察被归类为两年结果物种。传粉方式（风媒或虫媒）和叶习性（落叶或常绿）的特征也是根据表S1中列出的参考文献编译的。在这项研究中，叶特征被分为两类：落叶和常绿。然而，有四个物种被报道表现出中间状态，通常称为“亚常绿”（也记录为短落叶和半常绿）。在这种情况下，将三个状态分配给单一特征将阻止在BayesTraits（Pagel等人2004；Pagel和Meade 2006）中应用离散模型。在Hipp等人（2018）对橡树叶物候的详细评估中，Quercus engelmanii Greene、Q. fusiformis Small和Q. virginiana Mill被归类为短落叶，定义为“叶子全年存在或仅有短暂的落叶期”。鉴于叶子大部分时间都存在，我们将这些物种视为与落叶不同，并将其归类为常绿。对于Q. castaneifolia C.A.Mey，我们参考了额外的文献（见表S1），所有文献都一致将其归类为落叶物种；因此，我们将其视为落叶物种。利用在线标本库数据库中的标本图像，我们首先根据诊断性形态特征确认了物种身份，然后确定每个物种是一年结果还是两年结果。对于每个标本，我们检查了带有雌性生殖器官的枝条。使用叶痕作为形态标记来区分当年和前一年的枝条。为了分类果实特征，我们将来自同一采集事件的所有重复标本视为一个单一标本组。如果至少有一个标本组包含两个不同的雌性生殖器官发育阶段，则将该物种归类为两年结果：（i）前一年的花朵或果实；（ii）采集当年产生的花朵。不符合这一标准的物种被归类为一年结果（见表S3中的注释）。由于个别标本页通常只包含一个发育阶段，评估标本组而不是单个标本使我们能够检测到同一采集事件中的两个阶段。
2.2 系统发育重建
我们使用MrBayes v3.2.7（Ronquist等人2012）进行了贝叶斯系统发育分析，使用了Zhou等人（2022）提供的2124个核基因的核苷酸比对序列。不属于我们研究对象的序列被排除在外。我们进行了四次独立的马尔可夫链蒙特卡洛（MCMC）分析，每次分析运行1,000,000代，每500代采样一次。每次运行中前25%的采样树被丢弃作为燃烧期。剩余的燃烧后树被合并生成一个50%多数规则共识树，并计算节点的后验概率（图A1）。模型设置遵循Zhou等人（2022），采用PartitionFinder 2 v1.1（Lanfear等人2017）在赤池信息量准则（AIC）（Akaike 1974）下确定的最佳分割方案和核苷酸替换模型。考虑到后续祖先特征状态重建的拓扑不确定性，我们从燃烧后的树中随机重新采样了500棵树。
2.3 重建祖先状态和测试相关进化
我们应用BayesTraits v3.0.5（Pagel等人2004；Pagel和Meade 2006）来推断结果类型、传粉方式和叶习性的祖先状态，并使用多状态模型来测试这些特征之间的相关进化。多状态模型允许重建在系统发育树上存在的有限离散条件下的特征进化变化，从而便于祖先状态重建和进化模型的假设检验（Pagel等人2004）。离散模型用于通过比较依赖和独立的进化情景来测试两个离散二元特征是否以相关方式进化（Pagel 1994；Pagel和Meade 2006）。为了满足离散模型的要求，我们按照以下方式编码特征：结果特征（一年结果=0，两年结果=1），传粉方式（风媒=0，虫媒=1），以及叶特征（落叶=0，常绿=1）（见表S1）。在多状态模型框架内，我们比较了全速率不同模型（ARD），该模型估计所有状态之间的不同转换率，以及等速率模型（ER），后者将所有转换限制为相同，使用BayesTraits中的Restrict函数。对于结果特征，我们还拟合了群体异质性变体，允许主要谱系之间的转换率有所不同，考虑到果实特征在属水平上的显著变异（Satake和Kelly 2021）。为了验证这一点，我们使用了BayesTraits中的AddPattern函数，为四个物种丰富的支系（每个支系至少包含10个物种）分配了特定于该支系的转变率：Quercus subg. Quercus、Quercus subg. Cerris、Lithocarpus以及Castanea + Castanopsis（图1B）。这导致了52种不同的模型规格（图A2）。此外，一些Quercus物种在果实类型上表现出种内甚至个体内的多态性（例如，Q. suber L.：Díaz-Fernández等人2004年；Elena-Rossello等人1993年）。因此，我们对Q. suber在两种不同的假设下进行了祖先状态重建，分别将其视为1年或2年结果实的物种；这些数据被作为单独的数据集进行分析（“Fruiting_A”和“Fruiting_B”；表S4）。对于授粉方式和叶型，我们仅比较了等速率（ER）和所有速率不同（ARD）模型（即没有支系异质性变异）。在离散模型中（Pagel 1994），我们比较了四种不同的性状进化模型，以测试果实性状与其他性状（包括授粉和叶型性状）之间的相关性进化：（i）独立的ARD模型（4个参数），（ii）独立的ER模型（2个参数），（iii）依赖的ARD模型（8个参数），以及（iv）依赖的ER模型（4个参数）（图A3）。在独立模型中，这两个二元性状沿着独立的连续时间马尔可夫链独立进化，表明一个性状的变化不受另一个性状状态的影响。在依赖模型中，转变率是状态依赖的，即一个性状的变化概率取决于另一个性状的状态。这反映了一种情况，即一个性状的进化可能影响、限制或与另一个性状的进化协调。模型选择是使用从竞争模型的边际似然值计算出的贝叶斯因子（BF）来进行的（Kass和Raftery 1995）。我们计算了，定义为，其中表示边际似然值。所得到的值按照Kass和Raftery（1995）的方法进行解释：正值表示强烈支持更复杂模型的证据；非常强的证据（> 10）表示极强支持。所有BayesTraits分析都进行了1,250,000次迭代，前250,000次迭代被丢弃作为燃烧期，每1000次迭代进行一次采样。在MCMC运行的每一步中，从500个燃烧后重采样树的后验样本中随机选择一棵树，以整合拓扑结构和分支长度的不确定性。边际似然值是使用步石采样器（Xie等人2011）估计的，每个采样器有100次尝试和10,000次迭代；模型特定的值在表S4中报告。我们对转变率参数放置了一个指数超先验，其均值来自0到100之间的均匀分布，并使用可逆跳跃MCMC来遍历模型空间（Pagel和Meade 2006；Plummer等人2006）。使用R包“coda”评估了每个模型的MCMC链的收敛性，所有波动参数的有效样本量（ESS）超过了200。然后将后验状态概率映射到一个50%多数规则共识树上。我们使用R包ape中的chronos函数（Paradis等人2004；Paradis和Schliep 2019）和八个化石校准（Zhou等人2022年的年龄范围）对时间进行了缩放，以添加简化的时间轴。
3 结果
3.1 祖先特征状态重建
3.1.1 果实性状的估计进化历史
最适合果实性状祖先状态重建的模型是单转变率模型（表S4），在该模型中，1年和2年结果的转变在整个树中以相同的速率发生（0.757 ± 0.066；图2）。根据这个模型，节点N3——连接了除Fagus和Trigonobalanus之外的六个支系——被强烈支持为2年结果实的（86.8% ± 2.56%；图2；表S5）。相比之下，更深层次的节点N1和N2的祖先状态，即涉及Fagus和Trigonobalanus的早期分化，无法得到明确的解析（图2），反映了谱系基部的分辨率有限。图2显示了基于多状态模型分析中最佳支持结果的果实性状的估计进化历史，当Quercus suber被视为2年结果实时。粉色和绿色分别代表1年和2年结果实。系统发育树中每个节点的饼图表示1年和2年结果实的相对概率。系统发育树下方的值表示1年和2年结果实之间的相对转变率。在Lithocarpus、Chrysolepis和Notholithocarpus在节点N4-6分化后，2年结果实状态一直得到保持，每个节点的后验概率都超过了95%（图2；表S5）。在Castanea和Castanopsis支系内，祖先节点N9的结果是不确定的，后验概率为64.3% ± 4.76%，略微支持1年结果实。相比之下，所有其他包含2年结果实物种的共同祖先被推断为2年结果实，具有强烈的后验支持（Castanopsis：节点N10为80.2% ± 1.52%；Lithocarpus：“N12”为98.7% ± 1.06%；Chrysolepis：节点N13为99.9% ± 0.0702%；Quercus：节点N14为92.6% ± 4.45%；表S5）。相反，仅包含1年结果实物种的祖先被推断为1年结果实（Fagus：节点N7为98.3% ± 1.12%；Trigonobalanus：节点N8为79.6% ± 8.19%；Castanea：节点N11为99.2% ± 0.602%；表S5）。我们的分析还显示，从2年结果实到1年结果实的转变在多个谱系中独立发生了多次。1年结果实的物种主要来自Castanopsis、Quercus subg. Cerris和subg. Quercus的2年结果实支系。值得注意的是，Sects. Ponticae、Virentes和Quercus subg. Quercus的共同祖先被重建为1年结果实，后验概率为99.2% ± 0.644%，与这些部分的共同状态一致。此外，将Quercus suber视为1年结果实物种的分析得出的结果与上述结果一致（图A4）。
3.1.2 授粉性状的估计进化历史
多状态分析显示，选择了一个限制转变率的模型（表S4）。估计的相对转变率为0.200 ± 0.116（图3A）。对于授粉方式，包含Trigonobalanus及其姐妹支系的节点N2被强烈支持为虫媒传粉，后验概率为96.3% ± 2.46%，尽管Fagaceae的祖先状态是模糊的（图3A）。这种状态一直保持到Lithocarpus分化之后，之后Notholithocarpus和Quercus在节点N6的共同祖先被强烈支持为风媒传粉（88.9% ± 4.13%）。Fagus、Trigonobalanus doichangensis和Quercus的风媒传粉独立起源得到了进一步的支持（图3A）。
3.1.3 叶型性状的估计进化历史
对于叶型，选择了一个限制转变率的模型（表S4），估计的相对转变率为0.572 ± 0.226（图3B）。Fagaceae的祖先状态同样模糊，如节点N1所示，其落叶和常绿状态的可能性相当（图3B）。然而，包含Trigonobalanus及其姐妹支系的节点N2被强烈支持为常绿（94.5% ± 4.70%）。常绿属的共同祖先被支持为常绿，后验概率超过85.9%（Trigonobalanus：节点N8为85.9% ± 7.59%；Castanopsis：节点N10为98.7% ± 1.07%；Lithocarpus：节点N12为99.3% ± 0.716%；Chrysolepis：节点N13为99.9% ± 0.0507%；Quercus：节点N14为97.5% ± 1.68%；图3B），表明Castanea和某些Quercus谱系中的落叶性是一种衍生特征。由于Fagaceae的祖先状态仍未解决，因此无法确定Fagus的落叶习性是祖先状态还是衍生状态。
3.2 测试两种性状之间的相关进化
最后，我们检查了果实性状是否与授粉或叶型性状相关进化。基于似然比测试，独立ER模型被选为两种性状组合的最佳拟合模型（图4；表S4），表明没有证据支持相关性状进化。在两种分析中，果实、授粉和叶型性状的转变都是独立发生的，估计的转变率与分别为每个性状重建祖先状态时获得的转变率非常接近（图2-4；图A4）。图4显示了基于离散模型分析中最佳支持结果的果实性状与（A）授粉性状和（B）叶型性状之间相关进化的估计进化历史。系统发育树中每个节点的饼图表示四种状态的比例概率。系统发育树下方的值表示状态之间的相对转变率。（A）授粉性状的祖先状态重建。浅蓝色和黄色分别代表风媒传粉（AN）和虫媒传粉（EN）。（B）叶型性状的祖先状态重建。浅棕色和深绿色分别代表落叶（DEC）和常绿（EVE）。
4 讨论
我们的祖先状态重建表明，2年结果实特征起源于今天具有这一特征的属的最近期共同祖先（节点N3，连接了除Fagus和Trigonobalanus之外的六个支系：见图2；图A4）。先前的分子年代学研究将这一谱系的出现时间定在晚白垩世到早古新世（大约85–60百万年前；Liu等人2025年；Yang等人2025年；Zhou等人2022年）。由于这一时间框架与全球气候变化和季节性的出现相吻合（Toumoulin等人2022年），人们很容易推测这些环境变化可能促进了或驱动了延迟受精和由此产生的2年结果实的进化，因为在这些时期不适合繁殖。要更详细地理解这一进化背景，需要进行基于化石的生物地理学分析，明确将环境动态与延迟受精的起源和多样化联系起来。在最初获得2年结果实之后，整个家族中多次独立地恢复了1年结果实。这些恢复发生在Castanea的共同祖先、Castanopsis内部以及Quercus的几个支系中，表明尽管恢复到1年结果实的情况相对罕见，但它们在不同属中反复出现。在Lithocarpus中也可能存在类似的案例，其中一些未包含在我们的分析中的物种（例如，L. dodonaeifolius（Hayata）Hayata，L. formosanus（Skan）Hayata）已知具有1年结果实（Brach和Song 2006）。值得注意的是，在任何恢复到1年结果实的谱系中都没有检测到2年结果实的恢复案例。一个显著的例子涉及Quercus Sects. Ponticae、Virentes和Quercus的共同祖先，它出现在始新世（56–33.9百万年前）。尽管这个祖先被推断为经历了回到1年结果实的转变，但其随后在欧亚大陆和北美的分布扩展并没有伴随着2年结果实的二次获得（图2；图A4；Denk等人2017年；Hipp等人2018年，2020年；Zhou等人2022年）。同样，在Castanea中，所有现存物种都是1年结果实，没有恢复到2年结果实的案例。在Quercus suber中，已经报告了1年结果实和2年结果实在同一物种中共存（Elena-Rossello等人1993年；Díaz-Fernández等人2004年）。Díaz-Fernández等人（2004年）显示，在具有不同环境条件的多个地点，晚开花特征的个体倾向于表现出更高比例的2年结果实雌花。他们还报告说，1年结果实的个体比例随着纬度的增加而增加。这些发现表明，当延迟开花限制了在一个季节内完成果实发育的时间时，2年结果实可能更受青睐。这种模式与理论预测和观察到的开花和结果实物候一致（Satake和Kelly 2021年；Araye等人2023年），这表明在不利季节开花的物种更有可能表现出2年结果实。祖先状态重建强烈支持虫媒传粉和常绿叶型作为Quercoideae的祖先特征（节点N2；图3），早于最近共同祖先（节点N3）出现2年结果实策略。这种模式表明，2年结果实是在虫媒传粉和常绿叶型的背景下进化的。昆虫授粉是现存属如Castanea、Castanopsis、Lithocarpus、Chrysolepis、Notholithocarpus和Trigonobalanus verticillata的特征，这些属具有复花序，并能吸引广食性的昆虫传粉者（Chen等人2025；Crepet和Nixon 1989；Larue和Petit 2024；Larue等人2021；Kaul和Abbe 1984；Kaul 1986；Nixon和Crepet 1989；Petit和Larue 2022；Sadowski等人2020；Soepadmo 1972；Wright和Dodd 2013；Yuan等人2024），而Quercus和Trigonobalanus doichangensis的风力授粉则代表了衍生状态（Manos等人2001；Nixon和Crepet 1989；Zhou等人2022）；对于Fagus来说，其风力授粉的祖先状态或衍生状态仍然没有定论，尽管之前有研究表明Fagales和Fagaceae中的风力授粉是祖先特征（Stephens等人2023；Yao等人2023）。常绿叶习性可能是在温暖的白垩纪-古新世气候下获得的，并在低纬度的亚洲谱系中得到系统发育保守，如Castanopsis、Lithocarpus以及某些Quercus的亚属（Sects. Cyclobalanopsis和Ilex），而Castanea和几个Quercus谱系（Sects. Lobatae和Quercus）则独立地转变为落叶性，并伴随着向北的分布扩展，这可能是由于更冷和更季节性的气候（Barrón等人2017；Cavender-Bares 2018；Deng等人2018；Fontes等人2025；Hipp等人2018；Jiang等人2019；Kikuzawa 1991；Sancho-Knapik等人2021；Soepadmo 1972；Zanne等人2014）。从落叶到常绿状态的多次次级转变进一步表明，叶习性的进化既反映了系统发育的保守性，也反映了反复的适应性变化（Furze等人2021；Hipp等人2018；Sancho-Knapik等人2021）。尽管存在这些明显的进化动态，我们的离散性状分析并未发现果实结实策略与授粉方式或叶习性之间存在相关进化的证据。虽然转变次数本身限制了统计功效，但这一结果表明这些性状之间没有强烈的耦合。除了授粉方式和叶习性之外，与果实发育相关的资源投入成本可能是影响两年一次结实进化的重要因素。对果实组织的分配可能会施加重大的生理限制，从而可能影响繁殖时间。先前的理论研究（Satake和Kelly 2021）预测，当种子成熟所需的时间——与投入成本密切相关——延长时，延迟受精更有可能发生，这表明发育持续时间与两年一次结实进化之间存在机制上的联系。然而，来自现存物种的实证证据尚无定论。现有数据表明，1年结实和2年结实物种在杯果和果实中投入的碳、氮和磷的数量并没有一致差异（Wang等人2022）。此外，一些产生相对较大果实的物种仍然表现出1年结实（Wang等人2022），这表明仅靠投入成本不足以解释两年一次结实的进化。这些观察结果强调了需要采取更综合的方法来全面理解两年一次结实的进化驱动因素。同时，将这些根本驱动因素与直接的分子机制结合起来，对于全面理解延迟受精的机制至关重要。新兴的基因组学和转录组学方法现在使得能够在自然环境条件下研究发育过程是如何被调控的（Satake等人2022, 2023；Kudo等人2025；Lara-De La Cruz和Chávez-Vergara 2026）。例如，控制胚珠发育时间的基因以及涉及激素控制休眠及其释放的调控途径可能提供了季节性环境与繁殖时间之间的关键机制联系。因此，将比较系统发育分析与分子和发育研究相结合的未来研究对于阐明Fagaceae中延迟受精的进化过程至关重要。
作者贡献：
Shagawa Takenori：数据管理（负责人），正式分析（平等），验证（平等），可视化（平等），撰写——原始草稿（负责人），撰写——审阅和编辑（平等）。
Myotoishi Chihiro：数据管理（支持），正式分析（平等），方法学（平等），验证（平等），可视化（平等）。
Yahara Tetsukazu：数据管理（支持），监督（支持），撰写——审阅和编辑（平等）。
Imai Ryosuke：方法学（平等）。
Deng Min：数据管理（支持）。
Satake Akiko：概念化（负责人），资金获取（负责人），方法学（平等），项目管理（负责人），监督（负责人），验证（平等），撰写——审阅和编辑（平等）。
致谢：
我们感谢以下标本馆的管理员提供在线标本数据：A、ANDA、AU、BM、BR、CAF、CDBI、CSFI、E、GAC、GFS、GH、GZAC、HITBC、HK、IBK、IBSC、JXAU、K、KUN、L、MEL、MICH、MO、NAS、NY、P、PE、SHM、SN、SYS、SZ、US和WU。我们衷心感谢Nagahama Ai对本研究的宝贵讨论和深刻评论。我们还要感谢九州大学信息技术研究所提供的用于系统发育重建的计算资源，这些资源属于一般项目类别。本研究得到了日本学术振兴会（JSPS KAKENHI）的资助（JP23H04965和JP23H04966）。
利益冲突：
作者声明没有利益冲突。
附录A：
图A1：使用MrBayes v.3.2.7从Zhou等人（2022）的核基因组数据集中重建的多数规则共识树，在去除排除样本后。每个节点显示了后验概率（PP）。刻度条表示每个位点的替换次数。
图A2：多状态模型分析中考虑的52种替代模型设置的概览。模型ID显示在每棵树的左上角，分配不同转换率的分支用不同的颜色表示。
图A3：离散模型分析中四种不同模型配置的示意图。在独立模型（左侧）中，每个性状分别以不同的转换率进化，导致四个参数（不受限制）或两个参数（受限制）。在依赖模型（右侧）中，两个性状共同进化，导致八个参数（不受限制）或四个参数（受限制）。框表示两个性状可能的组合状态，箭头代表允许的转换及其相关的转换率参数。
图A4：基于多状态模型分析中最佳支持结果的果实性状估计进化历史，其中Quercus suber被视为1年结实。粉色和绿色分别代表1年和2年结实。系统发育树中的每个节点处的饼图表示1年和2年结实的比例概率。系统发育树下方的值表示1年和2年结实之间的相对转换率。
数据可用性声明：
本研究使用的性状数据可从中国eFlora（Brach和Song 2006）和北美eFlora（Nixon 1997）公开获取。用于性状数据汇编的额外文献来源列在表S1中。用于验证和细化性状分配的数字化标本信息提供在表S2和S3中。
