Gibbons are small arboreal apes that display an accelerated rate of evolutionary chromosomal rearrangement and occupy a key node in the primate phylogeny between Old World monkeys and great apes. Here we present the assembly and analysis of a northern white-cheeked gibbon (Nomascus leucogenys) genome. We describe the propensity for a gibbon-specific retrotransposon (LAVA) to insert into chromosome segregation genes and alter transcription by providing a premature termination site, suggesting a possible molecular mechanism for the genome plasticity of the gibbon lineage. We further show that the gibbon genera (Nomascus, Hylobates, Hoolock and Symphalangus) experienced a near-instantaneous radiation ~5 million years ago, coincident with major geographical changes in Southeast Asia that caused cycles of habitat compression and expansion. Finally, we identify signatures of positive selection in genes important for forelimb development (TBX5) and connective tissues (COL1A1) that may have been involved in the adaptation of gibbons to their arboreal habitat.

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...

Idiopathic scoliosis (IS) is a structural lateral spinal curvature of ≥10° that affects up to 3% of otherwise healthy children and can lead to life-long problems in severe cases. It is well-established that IS is a genetic disorder. Previous studies have identified genes that may contribute to the IS phenotype, but the overall genetic etiology of IS is not well understood. We used exome sequencing to study five multigenerational families with IS. Bioinformatic analyses identified unique and low frequency variants (minor allele frequency ≤5%) that were present in all sequenced members of the family. Across the five families, we identified a total of 270 variants with predicted functional consequences in 246 genes, and found that eight genes were shared by two families. We performed GO term enrichment analyses, with the hypothesis that certain functional annotations or pathways would be enriched in the 246 genes identified in our IS families. Using three complementary programs to complete these analyses, we identified enriched categories that include stereocilia and other actin-based cellular projections, cilia and other microtubule-based cellular projections, and the extracellular matrix (ECM). Our results suggest that there are multiple paths to IS and provide a foundation for future studies of IS pathogenesis.

Idiopathic scoliosis (IS) is a spinal disorder affecting up to 3% of otherwise healthy children. IS has a strong familial genetic component and is believed to be genetically complex due to significant variability in phenotype and heritability. Previous studies identified putative loci and variants possibly contributing to IS susceptibility, including within extracellular matrix, cilia, and actin networks, but the genetic architecture and underlying mechanisms remain unresolved. Here, we used whole‐exome sequencing from three affected individuals in a multigenerational family with IS and identified 19 uncommon variants (minor allele frequency < 0.05). Genotyping of additional family members identified a candidate heterozygous variant (H1115Q, G>C, rs142032413) within the ciliary gene KIF7, a regulator within the hedgehog (Hh) signaling pathway. Resequencing of the second cohort of unrelated IS individuals and controls identified several severe mutations in KIF7 in affected individuals only. Subsequently, we generated a mutant zebrafish model of kif7 using CRISPR‐Cas9. kif7co63/co63 zebrafish displayed severe scoliosis, presenting in juveniles and progressing through adulthood. We observed no deformities in the brain, Reissner fiber, or central canal cilia in kif7co63/co63 embryos, although alterations were seen in Hh pathway gene expression. This study suggests defects in KIF7‐dependent Hh signaling, which may drive pathogenesis in a subset of individuals with IS.