Transcriptomic and epigenomic characterization of the developing bat wing

Bats are the only mammals capable of powered flight, but little is known about the genetic determinants that shape their wings. Here, we generated a genome for Miniopterus natalensis and performed RNA-seq and ChIP-seq (H3K27ac, H3K27me3) on its developing forelimb and hindlimb autopods at sequential embryonic stages to decipher the molecular events that underlie bat wing development. Over 7,000 genes and several lncRNAs, including Tbx5-as1 and Hottip, were differentially expressed between forelimb, hindlimb and different stages. ChIP-seq identified thousands of regions that are differentially modified in forelimb versus hindlimb. Comparative genomics found 2,796 bat-accelerated regions within H3K27ac peaks, several of which cluster near limb-associated genes. Pathway analyses revealed multiple ribosomal proteins and known limb patterning signaling pathways as differentially regulated, and implicated increased forelimb mesenchymal condensations with differential growth. Combined, our work outlines multiple genetic components that contribute to bat wing formation, providing a genomic blueprint for this morphological innovation.