Published February 27, 2026 | Version v1

The different routes of parallel evolution in epiarenic growth in a hyperarid desert environment

Authors/Creators

  • 1. ROR icon Heidelberg University

Description

The Atacama Desert, one of the driest and oldest regions on Earth, represents an extreme environment that has driven remarkable adaptive evolution over millions of years. Within this setting, the genus Tillandsia comprises specialists capable of surviving at the dry limit of plant life. These species exhibit crassulacean acid metabolism (CAM), lack functional roots, and possess trichomes adapted for water and nutrient absorption. Of the more than 650 known Tillandsia species, nine occur exclusively in the Chilean–Peruvian Atacama Desert, where they colonize bare sand surfaces (epiarenic growth) under hyperarid conditions without significant rainfall, relying solely on nocturnal fog for moisture. Despite their striking adaptations, the evolutionary mechanisms and timing underlying the emergence of epiarenic growth remain poorly understood. Here, we reconstructed a maximum-likelihood phylogeny based on 278 plastome sequences of Tillandsioideae and other Bromeliaceae subfamilies to study respective sister species relationships and further divergence times estimates for the origin of epiarenic Tillandsia were estimated using bayesian inference in BEAST2. In addition, we analyzed variation in orthologous copies of the nuclear-encoded Agt1 gene to test for interspecific and interploidal gene flow among epiarenic taxa and their closest relatives. This gene has been previously established as a barcoding marker in bromeliads. The results indicate that epiarenic Tillandsia evolved multiple times independently from the Late Pliocene through the Pleistocene, consistent with a long-term hyperarid evolutionary arena promoting extreme adaptations. Moreover, evidence for frequent interspecific hybridization and gene flow suggests that hybridization may have also contributed to the long-term success of epiarenic Tillandsia species and is reflecting also the spatio-temporal dynamics of the Atacama’s hyperarid landscapes during the Pleistocene. 

Table of contents

SUPPLEMENTARY MATERIAL 1

Accession list, and metadata, Genebank codes

SUPPLEMENTARY MATERIAL 2

Taxon sampling of epiarenic sister species

SUPPLEMENTARY MATERIAL 3

Plastome alignment

SUPPLEMENTARY MATERIAL 4

ML substitution models

SUPPLEMENTARY MATERIAL 5

Agt1 alignment

SUPPLEMENTARY MATERIAL 6

Genome size estimates

SUPPLEMENTARY MATERIAL 7

Plastome ML tree

SUPPLEMENTARY MATERIAL 8

Plastome BEAST2 tree

SUPPLEMENTARY MATERIAL 9

Agt1 ML tree

Files

Supplementary Material 007new.pdf

Files (13.3 MB)

Name Size Download all
md5:cde880afdd61a729c8cdb5437d5a6ea4
140.8 kB Download
md5:ca7349d1d97bc0342702e90881362e23
15.7 kB Download
md5:06250a9c3e668ac0a460619938d55648
9.9 MB Download
md5:2cf241263efdcf90bcc0770f27dad3e8
18.7 kB Download
md5:808bcc83d05e8a8bf0442140326736ef
23.4 kB Download
md5:64da840491cb83f369d940ddd2c4591a
20.8 kB Download
md5:b3184a3238571b554794d88dc817c508
1.2 MB Preview Download
md5:55721cfef5258888568de89206986c30
1.6 MB Preview Download
md5:feeac1cb7d4c733457c9b895abc4b714
315.5 kB Preview Download