Unique expression patterns of multiple key genes associated with the evolution of mammalian flight

Analysis of changes in the action potential morphology of the mouse sinoatrial node true pacemaker cells during ontogenetic development in vitro and in silico

BACKGROUND

Maturation of the mouse is accompanied by the increase in heart rate. However, the mechanisms underlying this process remain unclear. We performed an action potentials (APs) recordings in mouse sinoatrial node (SAN) true pacemaker cells and in silico analysis to clarify the mechanisms underlying pre-postnatal period heart rate changes.

RESULTS

The APs of true pacemaker cells at different stages had similar configurations and dV/dtmax values. The cycle length, action potential duration (APD90 ), maximal diastolic potential (MDP), and AP amplitude decreased, meanwhile the velocity of diastolic depolarization (DDR) increased from E12.5 stage to adult. Using a pharmacological approach we found that in SAN true pacemaker cells ivabradine reduces the DDR and the cycle length significantly stronger in E12.5 than in newborn and adult mice, whereas the effects of Ni2+ and nifedipine were significantly stronger in adult mice. Computer simulations further suggested that the density of the hyperpolarization-activated pacemaker сurrent (If ) decreased during development, whereas transmembrane and intracellular Ca2+ flows increased.

CONCLUSIONS

The ontogenetic decrease in IK1 density from E12.5 to adult leads to depolarization of MDP to the voltage range in which calcium currents are activated, thereby shifting the balance from the "membrane-clock" to the "calcium-clock."

Novel Heredity Basis of the Four-Horn Phenotype in Sheep Using Genome-Wide Sequence Data

Horns are an important breeding trait for sheep. However, no widely recognized viewpoint on the regulatory genes and mechanisms of horns is available, and the genetic basis of the four-horn phenotype (FHP) is unclear. This work conducted a genome-wide association study with 100 sheep genomes from multiple breeds to investigate the genetic basis of the FHP. The results revealed three significant associations (corrected as p < 1.64 × 10-8) of the InDels (CHR2: g.133,742,709delA, g.133,743,215insC, and g.133,743,940delT) for FHP in the intergenic sequence (IGS) between the MTX2 and the LOC105609047 of CHR2. Moreover, 14 significant associations (corrected as p < 1.42 × 10-9) of SNPs with the FHP phenotype were identified in CHR2 and CHR16, including five (e.g., CHR16: g.40,351,378G > A and g.40,352,577G > A) located in the intron of the ADAMTS12 gene, eight (e.g., CHR2: g.133,727,513C > T and g.133,732,145T > G) in the IGS between MTX2 and LOC105609047, and only one (CHR2: g.133,930,761A > G) in the IGS between HOXD1 and MTX2. Obvious divergence was also observed in genotype patterns between the FHP and others (two horns and hornless) in the HOXD1 and ADAMTS12 gene regions. An extremely significant linkage also occurred between Loci I and Loci II within 100 individuals (LD = -156.02186, p < 0.00001). In summary, our study indicated that the genomic sequences from CHR2 and CHR16 contributed to the FHP in sheep, specifically the key candidate genes HOXD1 and ADAMTS12. These results improved our understanding of the Mendelian genetic basis of the FHP in sheep.

(Non)Parallel developmental mechanisms in vertebrate appendage reduction and loss

Appendages have been reduced or lost hundreds of times during vertebrate evolution. This phenotypic convergence may be underlain by shared or different molecular mechanisms in distantly related vertebrate clades. To investigate, we reviewed the developmental and evolutionary literature of appendage reduction and loss in more than a dozen vertebrate genera from fish to mammals. We found that appendage reduction and loss was nearly always driven by modified gene expression as opposed to changes in coding sequences. Moreover, expression of the same genes was repeatedly modified across vertebrate taxa. However, the specific mechanisms by which expression was modified were rarely shared. The multiple routes to appendage reduction and loss suggest that adaptive loss of function phenotypes might arise routinely through changes in expression of key developmental genes.

Making a bat: The developmental basis of bat evolution

Bats are incredibly diverse, both morphologically and taxonomically. Bats are the only mammalian group to have achieved powered flight, an adaptation that is hypothesized to have allowed them to colonize various and diverse ecological niches. However, the lack of fossils capturing the transition from terrestrial mammal to volant chiropteran has obscured much of our understanding of bat evolution. Over the last 20 years, the emergence of evo-devo in non-model species has started to fill this gap by uncovering some developmental mechanisms at the origin of bat diversification. In this review, we highlight key aspects of studies that have used bats as a model for morphological adaptations, diversification during adaptive radiations, and morphological novelty. To do so, we review current and ongoing studies on bat evolution. We first investigate morphological specialization by reviewing current knowledge about wing and face evolution. Then, we explore the mechanisms behind adaptive diversification in various ecological contexts using vision and dentition. Finally, we highlight the emerging work into morphological novelties using bat wing membranes.

Developmental and Evolutionary Allometry of the Mammalian Limb Skeleton

The variety of limb skeletal proportions enables a remarkable diversity of behaviors that include powered flight in bats and flipper-propelled swimming in whales using extremes of a range of homologous limb architectures. Even within human limbs, bone lengths span more than an order of magnitude from the short finger and toe bones to the long arm and leg bones. Yet all of this diversity arises from embryonic skeletal elements that are each a very similar size at formation. In this review article, I survey what is and is not yet known of the development and evolution of skeletal proportion at multiple hierarchical levels of biological organization. These include the cellular parameters of skeletal elongation in the cartilage growth plate, genes associated with differential growth, and putative gene regulatory mechanisms that would allow both covariant and independent evolution of the forelimbs and hindlimbs and of individual limb segments. Although the genetic mechanisms that shape skeletal proportion are still largely unknown, and most of what is known is limited to mammals, it is becoming increasingly apparent that the diversity of bone lengths is an emergent property of a complex system that controls elongation of individual skeletal elements using a genetic toolkit shared by all.

Accelerated Evolution of Limb-Related Gene Hoxd11 in the Common Ancestor of Cetaceans and Ruminants (Cetruminantia)

Reduced numbers of carpal and tarsal bones (wrist and ankle joints) are extensively observed in the clade of Cetacea and Ruminantia (Cetruminantia). Homebox D11 (Hoxd11) is one of the important genes required for limb development in mammals. Mutations in Hoxd11 can lead to defects in particular bones of limbs, including carpus and tarsus. To test whether evolutionary changes in Hoxd11 underlie the loss of these bones in Cetruminantia, we sequenced and analyzed Hoxd11 coding sequences and compared them with other 5′ HoxA and HoxD genes in a taxonomic coverage of Cetacea, Ruminantia and other mammalian relatives. Statistical tests on the Hoxd11 sequences found an accelerated evolution in the common ancestor of cetaceans and ruminants, which coincided with the reduction of carpal and tarsal bones in this clade. Five amino acid substitutions (G222S, G227A, G229S, A240T and G261V) and one amino acid deletion (G254Del) occurred in this lineage. In contrast, other 5′ HoxA and HoxD genes do not show this same evolutionary pattern, but instead display a highly conserved pattern of evolution in this lineage. Accelerated evolution of Hoxd11, but not other 5′ HoxA and HoxD genes, is probably related to the reduction of the carpal and tarsal bones in Cetruminantia. Moreover, we found two amino acid substitutions (G110S and D223N) in Hoxd11 that are unique to the lineage of Cetacea, which coincided with hindlimb loss in the common ancestor of cetaceans. Our results give molecular evidence of Hoxd11 adaptive evolution in cetaceans and ruminants, which could be correlated with limb morphological adaptation.

A metaanalysis of bat phylogenetics and positive selection based on genomes and transcriptomes from 18 species

This work represents a large, order-wide evolutionary analysis of the order Chiroptera (bats). Our pipeline for assembling sequence data and curating orthologous multiple sequence alignments includes methods for improving results when combining genomic and transcriptomic data sources. The resulting phylogenetic tree divides the order Chiroptera into Yinpterochiroptera and Yangochiroptera, in disagreement with the previous division into Megachiroptera and Microchiroptera and in agreement with some other recent molecular studies, and also provides evidence for other contested branch placements. We also performed a genome-wide analysis of positive selection and found 181 genes with signatures of positive selection. Enrichment analysis shows these positively selected genes to be primarily related to immune responses but also, surprisingly, collagen formation.

Historically, the evolution of bats has been analyzed using a small number of genetic loci for many species or many genetic loci for a few species. Here we present a phylogeny of 18 bat species, each of which is represented in 1,107 orthologous gene alignments used to build the tree. We generated a transcriptome sequence of Hypsignathus monstrosus, the African hammer-headed bat, and additional transcriptome sequence for Rousettus aegyptiacus, the Egyptian fruit bat. We then combined these data with existing genomic and transcriptomic data from 16 other bat species. In the analysis of such datasets, there is no clear consensus on the most reliable computational methods for the curation of quality multiple sequence alignments since these public datasets represent multiple investigators and methods, including different source materials (chromosomal DNA or expressed RNA). Here we lay out a systematic analysis of parameters and produce an advanced pipeline for curating orthologous gene alignments from combined transcriptomic and genomic data, including a software package: the Mismatching Isoform eXon Remover (MIXR). Using this method, we created alignments of 11,677 bat genes, 1,107 of which contain orthologs from all 18 species. Using the orthologous gene alignments created, we assessed bat phylogeny and also performed a holistic analysis of positive selection acting in bat genomes. We found that 181 genes have been subject to positive natural selection. This list is dominated by genes involved in immune responses and genes involved in the production of collagens.

Creating diversity in mammalian facial morphology: a review of potential developmental mechanisms

Mammals (class Mammalia) have evolved diverse craniofacial morphology to adapt to a wide range of ecological niches. However, the genetic and developmental mechanisms underlying the diversification of mammalian craniofacial morphology remain largely unknown. In this paper, we focus on the facial length and orofacial clefts of mammals and deduce potential mechanisms that produced diversity in mammalian facial morphology. Small-scale changes in facial morphology from the common ancestor, such as slight changes in facial length and the evolution of the midline cleft in some lineages of bats, could be attributed to heterochrony in facial bone ossification. In contrast, large-scale changes of facial morphology from the common ancestor, such as a truncated, widened face as well as the evolution of the bilateral cleft possessed by some bat species, could be brought about by changes in growth and patterning of the facial primordium (the facial processes) at the early stages of embryogenesis.

Timing the developmental origins of mammalian limb diversity

Mammals have highly diverse limbs that have contributed to their occupation of almost every niche. Researchers have long been investigating the development of these diverse limbs, with the goals of identifying developmental processes and potential biases that shape mammalian limb diversity. To date, researchers have used techniques ranging from the genomic to the anatomic to investigate the developmental processes shaping the limb morphology of mammals from five orders (Marsupialia, Chiroptera, Rodentia, Cetartiodactyla, and Perissodactyla). Results of these studies suggest that the differential expression of genes controlling diverse cellular processes underlies mammalian limb diversity. Results also suggest that the earliest development of the limb tends to be conserved among mammalian species, while later limb development tends to be more variable. This research has established the mammalian limb as a model system for evolutionary developmental biology, and set the stage for more in-depth, cross-disciplinary research into the genetic controls, tissue-level cellular behaviors, and selective pressures that have driven the developmental evolution of mammalian limbs. Ideally, these studies will be performed in a diverse suite of mammalian species within a comparative, phylogenetic framework.

Transcriptomic insights into the genetic basis of mammalian limb diversity

Background

From bat wings to whale flippers, limb diversification has been crucial to the evolutionary success of mammals. We performed the first transcriptome-wide study of limb development in multiple species to explore the hypothesis that mammalian limb diversification has proceeded through the differential expression of conserved shared genes, rather than by major changes to limb patterning. Specifically, we investigated the manner in which the expression of shared genes has evolved within and among mammalian species.

Results

We assembled and compared transcriptomes of bat, mouse, opossum, and pig fore- and hind limbs at the ridge, bud, and paddle stages of development. Results suggest that gene expression patterns exhibit larger variation among species during later than earlier stages of limb development, while within species results are more mixed. Consistent with the former, results also suggest that genes expressed at later developmental stages tend to have a younger evolutionary age than genes expressed at earlier stages. A suite of key limb-patterning genes was identified as being differentially expressed among the homologous limbs of all species. However, only a small subset of shared genes is differentially expressed in the fore- and hind limbs of all examined species. Similarly, a small subset of shared genes is differentially expressed within the fore- and hind limb of a single species and among the forelimbs of different species.

Conclusions

Taken together, results of this study do not support the existence of a phylotypic period of limb development ending at chondrogenesis, but do support the hypothesis that the hierarchical nature of development translates into increasing variation among species as development progresses.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0902-6) contains supplementary material, which is available to authorized users.

Limb development: a paradigm of gene regulation

The limb is a commonly used model system for developmental biology. Given the need for precise control of complex signalling pathways to achieve proper patterning, the limb is also becoming a model system for gene regulation studies. Recent developments in genomic technologies have enabled the genome-wide identification of regulatory elements that control limb development, yielding insights into the determination of limb morphology and forelimb versus hindlimb identity. The modulation of regulatory interactions - for example, through the modification of regulatory sequences or chromatin architecture - can lead to morphological evolution, acquired regeneration capacity or limb malformations in diverse species, including humans.

The origins, scaling and loss of tetrapod digits

Many of the great morphologists of the nineteenth century marvelled at similarities between the limbs of diverse species, and Charles Darwin noted these homologies as significant supporting evidence for descent with modification from a common ancestor. Sir Richard Owen also took great care to highlight each of the elements of the forelimb and hindlimb in a multitude of species with focused attention on the homology between the hoof of the horse and the middle digit of man. The ensuing decades brought about a convergence of palaeontology, experimental embryology and molecular biology to lend further support to the homologies of tetrapod limbs and their developmental origins. However, for all that we now understand about the conserved mechanisms of limb development and the development of gross morphological disturbances, little of what is presented in the experimental or medical literature reflects the remarkable diversity resulting from the 450 million year experiment of natural selection. An understanding of conserved and divergent limb morphologies in this new age of genomics and genome engineering promises to reveal more of the developmental potential residing in all limbs and to unravel the mechanisms of evolutionary variation in limb size and shape. In this review, we present the current state of our rapidly advancing understanding of the evolutionary origin of hands and feet and highlight what is known about the mechanisms that shape diverse limbs.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

Bat Accelerated Regions Identify a Bat Forelimb Specific Enhancer in the HoxD Locus

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.

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 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 differentially modified in forelimb and hindlimb. Comparative genomics found 2,796 bat-accelerated regions within H3K27ac peaks, several of which cluster near limb-associated genes. Pathway analyses highlighted multiple ribosomal proteins and known limb patterning signaling pathways as differentially regulated and implicated increased forelimb mesenchymal condensation in differential growth. In combination, our work outlines multiple genetic components that likely contribute to bat wing formation, providing insights into this morphological innovation.

Geometric Morphometrics on Gene Expression Patterns Within Phenotypes: A Case Example on Limb Development

How the genotype translates into the phenotype through development is critical to fully understand the evolution of phenotypes. We propose a novel approach to directly assess how changes in gene expression patterns are associated with changes in morphology using the limb as a case example. Our method combines molecular biology techniques, such as whole-mount in situ hybridization, with image and shape analysis, extending the use of Geometric Morphometrics to the analysis of nonanatomical shapes, such as gene expression domains. Elliptical Fourier and Procrustes-based semilandmark analyses were used to analyze the variation and covariation patterns of the limb bud shape with the expression patterns of two relevant genes for limb morphogenesis, Hoxa11 and Hoxa13. We devised a multiple thresholding method to semiautomatically segment gene domains at several expression levels in large samples of limb buds from C57Bl6 mouse embryos between 10 and 12 postfertilization days. Besides providing an accurate phenotyping tool to quantify the spatiotemporal dynamics of gene expression patterns within developing structures, our morphometric analyses revealed high, non-random, and gene-specific variation undergoing canalization during limb development. Our results demonstrate that Hoxa11 and Hoxa13, despite being paralogs with analogous functions in limb patterning, show clearly distinct dynamic patterns, both in shape and size, and are associated differently with the limb bud shape. The correspondence between our results and already well-established molecular processes underlying limb development confirms that this morphometric approach is a powerful tool to extract features of development regulating morphogenesis. Such multilevel analyses are promising in systems where not so much molecular information is available and will advance our understanding of the genotype–phenotype map. In systematics, this knowledge will increase our ability to infer how evolution modified a common developmental pattern to generate a wide diversity of morphologies, as in the vertebrate limb.

Retinoic acid-independent expression of Meis2 during autopod patterning in the developing bat and mouse limb

BACKGROUND

The bat has strikingly divergent forelimbs (long digits supporting wing membranes) and hindlimbs (short, typically free digits) due to the distinct requirements of both aerial and terrestrial locomotion. During embryonic development, the morphology of the bat forelimb deviates dramatically from the mouse and chick, offering an alternative paradigm for identifying genes that play an important role in limb patterning.

RESULTS

Using transcriptome analysis of developing Natal long-fingered bat (Miniopterus natalensis) fore- and hindlimbs, we demonstrate that the transcription factor Meis2 has a significantly higher expression in bat forelimb autopods compared to hindlimbs. Validation by reverse transcriptase and quantitative polymerase chain reaction (RT-qPCR) and whole mount in situ hybridisation shows that Meis2, conventionally known as a marker of the early proximal limb bud, is upregulated in the bat forelimb autopod from CS16. Meis2 expression is localised to the expanding interdigital webbing and the membranes linking the wing to the hindlimb and tail. In mice, Meis2 is also expressed in the interdigital region prior to tissue regression. This interdigital Meis2 expression is not activated by retinoic acid (RA) signalling as it is present in the retained interdigital tissue of Rdh10 (trex/trex) mice, which lack RA. Additionally, genes encoding RA-synthesising enzymes, Rdh10 and Aldh1a2, and the RA nuclear receptor Rarβ are robustly expressed in bat fore- and hindlimb interdigital tissues indicating that the mechanism that retains interdigital tissue in bats also occurs independently of RA signalling.

CONCLUSIONS

Mammalian interdigital Meis2 expression, and upregulation in the interdigital webbing of bat wings, suggests an important role for Meis2 in autopod development. Interdigital Meis2 expression is RA-independent, and retention of interdigital webbing in bat wings is not due to the suppression of RA-induced cell death. Rather, RA signalling may play a role in the thinning (rather than complete loss) of the interdigital tissue in the bat forelimb, while Meis2 may interact with other factors during both bat and mouse autopod development to maintain a pool of interdigital cells that contribute to digit patterning and growth.

How the pterosaur got its wings

Throughout the evolutionary history of life, only three vertebrate lineages took to the air by acquiring a body plan suitable for powered flight: birds, bats, and pterosaurs. Because pterosaurs were the earliest vertebrate lineage capable of powered flight and included the largest volant animal in the history of the earth, understanding how they evolved their flight apparatus, the wing, is an important issue in evolutionary biology. Herein, I speculate on the potential basis of pterosaur wing evolution using recent advances in the developmental biology of flying and non-flying vertebrates. The most significant morphological features of pterosaur wings are: (i) a disproportionately elongated fourth finger, and (ii) a wing membrane called the brachiopatagium, which stretches from the posterior surface of the arm and elongated fourth finger to the anterior surface of the leg. At limb-forming stages of pterosaur embryos, the zone of polarizing activity (ZPA) cells, from which the fourth finger eventually differentiates, could up-regulate, restrict, and prolong expression of 5'-located Homeobox D (Hoxd) genes (e.g. Hoxd11, Hoxd12, and Hoxd13) around the ZPA through pterosaur-specific exploitation of sonic hedgehog (SHH) signalling. 5'Hoxd genes could then influence downstream bone morphogenetic protein (BMP) signalling to facilitate chondrocyte proliferation in long bones. Potential expression of Fgf10 and Tbx3 in the primordium of the brachiopatagium formed posterior to the forelimb bud might also facilitate elongation of the phalanges of the fourth finger. To establish the flight-adapted musculoskeletal morphology shared by all volant vertebrates, pterosaurs probably underwent regulatory changes in the expression of genes controlling forelimb and pectoral girdle musculoskeletal development (e.g. Tbx5), as well as certain changes in the mode of cell-cell interactions between muscular and connective tissues in the early phase of their evolution. Developmental data now accumulating for extant vertebrate taxa could be helpful in understanding the cellular and molecular mechanisms of body-plan evolution in extinct vertebrates as well as extant vertebrates with unique morphology whose embryonic materials are hard to obtain.

Differential expression of Meis2, Mab21l2 and Tbx3 during limb development associated with diversification of limb morphology in mammals

Bats are the only mammals capable of self-powered flight using wings. Differing from mouse or human limbs, four elongated digits within a broad wing membrane support the bat wing, and the foot of the bat has evolved a long calcar that spread the interfemoral membrane. Our recent mRNA sequencing (mRNA-Seq) study found unique expression patterns for genes at the 5' end of the Hoxd gene cluster and for Tbx3 that are associated with digit elongation and wing membrane growth in bats. In this study, we focused on two additional genes, Meis2 and Mab21l2, identified from the mRNA-Seq data. Using whole-mount in situ hybridization (WISH) we validated the mRNA-Seq results for differences in the expression patterns of Meis2 and Mab21l2 between bat and mouse limbs, and further characterize the timing and location of the expression of these two genes. These analyses suggest that Meis2 may function in wing membrane growth and Mab21l2 may have a role in AP and DV axial patterning. In addition, we found that Tbx3 is uniquely expressed in the unique calcar structure found in the bat hindlimb, suggesting a role for this gene in calcar growth and elongation. Moreover, analysis of the coding sequences for Meis2, Mab21l2 and Tbx3 showed that Meis2 and Mab21l2 have high sequence identity, consistent with the functions of genes being conserved, but that Tbx3 showed accelerated evolution in bats. However, evidence for positive selection in Tbx3 was not found, which would suggest that the function of this gene has not been changed. Together, our findings support the hypothesis that the modulation of the spatiotemporal expression patterns of multiple functional conserved genes control limb morphology and drive morphological change in the diversification of mammalian limbs.