Functional convergence of regulatory regions provides vital insights into mammalian gliding adaptation

Uncovering the key genetic basis of complex phenotypic convergence in distantly related species has been a long-standing focus in evolutionary biology and genetics, and the convergent evolution of gliding in mammals offers a valuable opportunity to address this question. Here, we investigated the genomic basis of convergent evolution of gliding in mammals by analyzing both protein-coding genes and conserved non-coding elements (CNEs). We first de novo assembled and annotated two chromosome-level genomes of gliding mammals, the red and white giant flying squirrel (Petaurista alborufus) and sugar gliders (Petaurus breviceps), and conducted comprehensive comparative genomic analysis combined with another gliding mammal, the Sunda flying lemur (Galeopterus variegatus) and 14 background species. We found that the convergent evolution of protein-coding genes provided relatively limited but functionally relevant evidence linked to gliding phenotypes. By contrast, we found that gliding-accelerated CNEs (GACNEs) cluster near functionally equivalent genes and frequently aggregate into highly diverged yet functionally convergent hotspot regions. Across the three gliding lineages, both GACNEs and hotspot GACNEs show strong convergence in their functional enrichment profiles, suggesting a broad genetic basis underlying the convergent gliding phenotype. Furthermore, we identified 72 core transcription factors underpinning the genetic basis of gliding convergence, including EMX2 and ZFHX3, potentially involved in multiple aspects of gliding adaptation. Our study highlights the role of functional convergence in regulatory regions as a key mechanism in mammalian gliding convergence, offering valuable insights and strategies for uncovering the genetic basis of complex convergent traits, thereby advancing understanding of the molecular basis of convergent traits.