RESUMO
BACKGROUND AND AIMS: Heterotrophic plants have long been a challenge for systematists, exemplified by the base of the orchid subfamily Epidendroideae, which contains numerous mycoheterotrophic species. METHODS: Here we address the utility of organellar genomes in resolving relationships at the epidendroid base, specifically employing models of heterotachy, or lineage-specific rate variation over time. We further conduct comparative analyses of plastid genome evolution in heterotrophs and structural variation in matK. KEY RESULTS: We present the first complete plastid genomes (plastomes) of Wullschlaegelia, the sole genus of the tribe Wullschlaegelieae, revealing a highly reduced genome of 37 kilobases, which retains a fraction of the genes present in related autotrophs. Plastid phylogenomic analyses recovered a strongly supported clade composed exclusively of mycoheterotrophic species with long branches. We further analyzed mitochondrial gene sets, which recovered similar relationships to those in other studies using nuclear data, but the placement of Wullschlaegelia remains uncertain. We conducted comparative plastome analyses among Wullschlaegelia and other heterotrophic orchids, revealing a suite of correlated substitutional and structural changes relative to autotrophic species. Lastly, we investigated evolutionary and structural variation in matK, which is retained in Wullschlaegelia and a few other 'late stage' heterotrophs and found evidence for structural conservation despite rapid substitution rates in both Wullschlaegelia and the leafless Gastrodia. CONCLUSIONS: Our analyses reveal the limits of what the plastid genome can tell us on orchid relationships in this part of the tree, even when applying parameter-rich heterotachy models. Our study underscores the need for increased taxon sampling across all three genomes at the epidendroid base, and illustrates the need for further research on addressing heterotachy in phylogenomic analyses.