The order Chiroptera, commonly known as bats, is the only group of mammals to have evolved the capability of flight. They are estimated to have diverged from their arboreal ancestors ~51 million years ago 1 . Their adaptions for flight include substantial specialization of the forelimb, characterized by the notable extension of digits II-V, a decrease in wing bone mineralization along the proximal-distal axis, and the retention and expansion of interdigit webbing, which is controlled by a novel complex of muscles 2,3 . Bat hindlimbs are comparatively short, with free, symmetrical digits, providing an informative contrast that can be used to highlight the genetic processes involved in bat wing formation. Previous studies that examined gene expression in developing bat forelimbs and hindlimbs reported differential expression of several genes, including Tbx3, Brinp3, Meis2, the 5′ HoxD genes and components of the Shh-Fgf signaling loop, suggesting that multiple genes and processes are involved in generating these morphological innovations [4][5][6][7][8] . Gene regulatory elements are thought to be important drivers of these changes: for example, replacement of the mouse Prx1 limb enhancer with the equivalent bat sequence resulted in elongated forelimbs 9 . However, an integrated understanding of how changes in regulatory elements, various genes and signaling pathways combine to collectively shape the bat wing remains largely elusive.To characterize the genetic differences that underlie divergence in bat forelimb and hindlimb development, we used a comprehensive, genome-wide strategy. We generated a de novo whole-genome assembly for the vesper bat, M. natalensis, for which a well-characterized stage-by-stage morphological comparison between developing bat and mouse limbs is available 10 . In this species, the developing forelimb noticeably diverges from the hindlimb from developmental stages CS15 and CS16, with clear morphological differences seen at a subsequent stage, CS17 (ref. 10). This developmental window is equivalent to embryonic day (E) 12.0 to E13.5 in mouse 4,10 . M. natalensis embryos were obtained and transcriptomic (RNA-seq) data and ChIP-seq data for both an active (acetylation of histone H3 at lysine 27, H3K27ac; refs. 11,12) and a repressive (trimethylation of histone H3 at lysine 27, H3K27me3; ref. 13) mark were generated for these three developmental stages (Fig. 1). 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 and H3K27me3) analyses 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 long noncoding RNAs, including Tbx5-as1 and Hottip, were differentially expressed between forelimb and hindlimb, and across different stages. ChIP-seq analysis identified thousands of regions that are differentiall...

SummaryBacteria differ from eukaryotes by having the enzyme DNA gyrase, which catalyses the ATP-dependent negative supercoiling of DNA. Negative supercoils are essential for condensing chromosomes into an interwound (plectonemic) and branched structure known as the nucleoid. Topo-1 removes excess supercoiling in an ATP-independent reaction and works with gyrase to establish a topological equilibrium where supercoils move within 10 kb domains bounded by stochastic barriers along the sequence. However, transcription changes the stochastic pattern by generating supercoil diffusion barriers near the sites of gene expression. Using supercoildependent Tn3 and gd resolution assays, we studied DNA topology upstream, downstream and across highly transcribed operons. Whenever two Res sites flanked efficiently transcribed genes, resolution was inhibited and the loss in recombination efficiency was proportional to transcription level. Ribosomal RNA operons have the highest transcription rates, and resolution assays at the rrnG and rrnH operons showed inhibitory levels 40-100 times those measured in low-transcription zones. Yet, immediately upstream and downstream of RNA polymerase (RNAP) initiation and termination sites, supercoiling characteristics were similar to poorly transcribed zones. We present a model that explains why RNAP blocks plectonemic supercoil movement in the transcribed track and suggests how gyrase and TopA control upstream and downstream transcriptiondriven supercoiling.

The molecular events leading to the development of the bat wing remain largely unknown, and are thought to be caused, in part, by changes in gene expression during limb development. These expression changes could be instigated by variations in gene regulatory enhancers. Here, we used a comparative genomics approach to identify regions that evolved rapidly in the bat ancestor, but are highly conserved in other vertebrates. We discovered 166 bat accelerated regions (BARs) that overlap H3K27ac and p300 ChIP-seq peaks in developing mouse limbs. Using a mouse enhancer assay, we show that five Myotis lucifugus BARs drive gene expression in the developing mouse limb, with the majority showing differential enhancer activity compared to the mouse orthologous BAR sequences. These include BAR116, which is located telomeric to the HoxD cluster and had robust forelimb expression for the M. lucifugus sequence and no activity for the mouse sequence at embryonic day 12.5. Developing limb expression analysis of Hoxd10-Hoxd13 in Miniopterus natalensis bats showed a high-forelimb weak-hindlimb expression for Hoxd10-Hoxd11, similar to the expression trend observed for M. lucifugus BAR116 in mice, suggesting that it could be involved in the regulation of the bat HoxD complex. Combined, our results highlight novel regulatory regions that could be instrumental for the morphological differences leading to the development of the bat wing.

No abstract

The molecular events leading to the development of the bat wing remain largely unknown, and are thought to be caused, in part, by changes in gene expression during limb development. These expression changes could be instigated by variations in gene regulatory enhancers. Here, we used a comparative genomics approach to identify regions that evolved rapidly in the bat ancestor, but are highly conserved in other vertebrates. We discovered 166 bat accelerated regions (BARs) that overlap H3K27ac and p300 ChIP-seq peaks in developing mouse limbs. Using a mouse enhancer assay, we show that five Myotis lucifugus BARs drive gene expression in the developing mouse limb, with the majority showing differential enhancer activity compared to the mouse orthologous BAR sequences. These include BAR116, which is located telomeric to the HoxD cluster and had robust forelimb expression for the M. lucifugus sequence and no activity for the mouse sequence at embryonic day 12.5. Developing limb expression analysis of Hoxd10-Hoxd13 in Miniopterus natalensis bats showed a high-forelimb weak-hindlimb expression for Hoxd10-Hoxd11, similar to the expression trend observed for M. lucifugus BAR116 in mice, suggesting that it could be involved in the regulation of the bat HoxD complex. Combined, our results highlight novel regulatory regions that could be instrumental for the morphological differences leading to the development of the bat wing. Author SummaryThe limb is a classic example of vertebrate homology and is represented by a large range of morphological structures such as fins, legs and wings. The evolution of these structures could be driven by alterations in gene regulatory elements that have critical roles during development. To identify elements that may contribute to bat wing development, we characterized sequences that are conserved between vertebrates, but changed significantly in the bat lineage. We then overlapped these sequences with predicted developing limb enhancers as determined by ChIP-seq, finding 166 bat accelerated sequences (BARs). Five BARs that were tested for enhancer activity in mice all drove expression in the limb. Testing the mouse orthologous sequence showed that three had differences in their limb enhancer activity as compared to the bat sequence. Of these, BAR116 was of particular interest as it is located near the HoxD locus, an essential gene complex required for proper spatiotemporal patterning of the developing limb. The bat BAR116 sequence drove robust forelimb expression but the mouse BAR116 sequence did not show enhancer activity. These experiments correspond to analyses of HoxD gene expressions in developing bat limbs, which had strong forelimb versus weak hindlimb expression for Hoxd10-11. Combined, our studies highlight specific genomic regions that could be important in shaping the morphological differences that led to the development of the bat wing.