Search

The Mediterranean lizardPodarcis lilfordiis an emblematic species of the Balearic Islands. The extensive phenotypic diversity among extant isolated populations makes the species a great insular model system for eco-evolutionary studies, as well as a challenging target for conservation management plans. Here we report the first high quality chromosome-level assembly and annotation of theP. lilfordigenome, along with its mitogenome, based on a mixed sequencing strategy (10X Genomics linked reads, Oxford Nanopore Technologies long reads and Hi-C scaffolding) coupled with extensive transcriptomic data (Illumina and PacBio). The genome assembly (1.5 Gb) is highly contiguous (N50 = 90 Mb) and complete, with 99% of the sequence assigned to candidate chromosomal sequences and >97% gene completeness. We annotated a total of 25,663 protein-coding genes, assigning 72% to known functions. Comparison to the genome of the related speciesPodarcis muralisrevealed substantial similarity in genome size, annotation metrics, repeat content, and strong collinearity, despite their evolutionary distance (~18-20 MYA). This genome expands the repertoire of available reptilian genomes and will facilitate the exploration of the molecular and evolutionary processes underlying the extraordinary phenotypic diversity of this insular species, while providing a critical resource for conservation genomics.

The Mediterranean lizard Podarcis lilfordi is an emblematic species of the Balearic Islands. The extensive phenotypic diversity among extant isolated populations makes the species a great insular model system for eco-evolutionary studies, as well as a challenging target for conservation management plans. Here we report the first high-quality chromosome-level assembly and annotation of the P. lilfordi genome, along with its mitogenome, based on a mixed sequencing strategy (10X Genomics linked reads, Oxford Nanopore Technologies long reads and Hi-C scaffolding) coupled with extensive transcriptomic data (Illumina and PacBio). The genome assembly (1.5 Gb) is highly contiguous (N50 = 90 Mb) and complete, with 99% of the sequence assigned to candidate chromosomal sequences and >97% gene completeness. We annotated a total of 25,663 protein-coding genes translating into 38,615 proteins. Comparison to the genome of the related species Podarcis muralis revealed substantial similarity in genome size, annotation metrics, repeat content, and a strong collinearity, despite their evolutionary distance (~18-20 MYA). This genome expands the repertoire of available reptilian genomes and will facilitate the exploration of the molecular and evolutionary processes underlying the extraordinary phenotypic diversity of this insular species, while providing a critical resource for conservation genomics. This study was supported by the Institut d’Estudis Catalans under the Catalan Initiative for the Earth Biogenome Project (PRO2020-S02 to L.B.), the Swedish Research Council (VR 2017-03846 and VR-2021-04656 to T.U. and VR-2020-03650 to N.F.) and Starting Grant from the European Research Council (no. 948126 to N.F.). We also acknowledge the support of the Spanish Ministry of Science and Innovation to the EMBL partnership, the Centro de Excelencia Severo Ochoa, the CERCA Programme/Generalitat de Catalunya, the Spanish Ministry of Science and Innovation through the Instituto de Salud Carlos III and the Generalitat de Catalunya through Departament de Salut and Departament d’Empresa i Coneixement. Co-financing funds were obtained from the European Regional Development Fund by the Spanish Ministry of Science and Innovation corresponding to the Programa Operativo FEDER Plurirregional de España (POPE) 2014-2020 and by the Secretaria d’Universitats i Recerca, Departament d’Empresa i Coneixement of the Generalitat de Catalunya corresponding to the Programa Operatiu FEDER de Catalunya 2014-2020.

Animal coloration is often expressed in periodic patterns that can arise from differential cell migration, yet how these processes are regulated remains elusive. We show that a female-limited polymorphism in dorsal patterning (diamond/chevron) in the brown anole is controlled by a single Mendelian locus. This locus contains the gene

Environmental stress can result in epigenetic modifications that are passed down several generations. Such epigenetic inheritance can have significant impact on eco-evolutionary dynamics, but the phenomenon remains controversial in ecological model systems. Here, we used whole-genome bisulfite sequencing on individual water fleas (Daphnia magna) to assess whether environmentally-induced DNA methylation can persist for up to four generations. Genetically identical females were exposed to a control treatment, one of three natural stressors (high temperature, zinc, microcystin), or the methylation-inhibitor 5-azacytidine. After exposure, lines were propagated clonally for four generations under control conditions. We identified between 70 and 225 differentially methylated CpG positions (DMPs) between controls and F1 individuals whose mothers (and therefore they themselves as germ cells) were exposed to one of the three natural stressors. Between 46% and 58% of these environmentally-induced DMPs persisted until generation F4 without attenuation in their magnitude of differential methylation. DMPs were enriched in exons and largely stressor-specific, suggesting a possible role in environment-dependent gene regulation. In contrast, treatment with the compound 5-azacytidine demonstrated that pervasive hypo-methylation upon exposure is reset almost completely after a single generation. These results suggest that environmentally-induced DNA methylation is non-random and stably inherited across generations in Daphnia, making epigenetic inheritance a putative factor in the eco-evolutionary dynamics of fresh-water communities.Author summaryWater fleas are important keystone species mediating eco-evolutionary dynamics in lakes and ponds. It is currently an open question in how far epigenetic inheritance contributes to the ability of Daphnia populations to adapt to environmental stress. Using a range of naturally occurring stressors and a multi-generational design, we show that environmentally-induced DNA methylation variants are stably inherited for at least four generations in Daphnia magna. The induced variation in DNA methylation are stressor-specific and almost exclusively found in exons, bearing the signatures of functional adaptations. Our findings imply that ecological adaptations of Daphnia to seasonal fluctuations can be underpinned by epigenetic inheritance of DNA methylation without changes in gene frequencies.

Rapid technological improvements are democratizing access to high quality, chromosome-scale genome assemblies. No longer the domain of only the most highly studied model organisms, now non-traditional and emerging model species can be genome-enabled using a combination of sequencing technologies and assembly software. Consequently, old ideas built on sparse sampling across the tree of life have recently been amended in the face of genomic data drawn from a growing number of high-quality reference genomes. Arguably the most valuable are those long-studied species for which much is already known about their biology; what many term emerging model species. Here, we report a highly complete chromosome-scale genome assembly for the brown anole, Anolis sagrei - a lizard species widely studied across a variety of disciplines and for which a high-quality reference genome was long overdue. This assembly exceeds the vast majority of existing reptile and snake genomes in contiguity (N50 = 253.6 Mb) and annotation completeness. Through the analysis of this genome and population resequence data, we examine the history of repetitive element accumulation, identify the X chromosome, and propose a hypothesis for the evolutionary history of fusions between autosomes and the X that led to the sex chromosomes of A. sagrei.

Rapid technological improvements are democratizing access to high quality, chromosome-scale genome assemblies. No longer the domain of only the most highly studied model organisms, now non-traditional and emerging model species can be genome-enabled using a combination of sequencing technologies and assembly software. Consequently, old ideas built on sparse sampling across the tree of life have recently been amended in the face of genomic data drawn from a growing number of high-quality reference genomes. Arguably the most valuable are those long-studied species for which much is already known about their biology; what many term emerging model species. Here, we report a new, highly complete chromosome-scale genome assembly for the brown anole, Anolis sagrei – a lizard species widely studied across a variety of disciplines and for which a high-quality reference genome was long overdue.

Hox genes are a remarkable example of conservation in animal development and their nested expression along the head-to-tail axis orchestrates embryonic patterning. Early in vertebrate history, two duplications led to the emergence of four Hox clusters (A-D) and redundancy within paralog groups has been partially accommodated with gene losses. Here we conduct an inventory of squamate Hox genes using the genomes of 10 lizard and 7 snake species. Although the HoxC1 gene has been hypothesized to be lost in the amniote ancestor, we reveal that it is retained in lizards. In contrast, all snakes lack functional HoxC1 and -D12 genes. Varying levels of degradation suggest differences in the process of gene loss between the two genes. The vertebrate HoxC1 gene is prone to gene loss and its functional domains are more variable than those of other Hox1 genes. We describe for the first time the HoxC1 expression patterns in tetrapods. HoxC1 is broadly expressed during development in the diencephalon, the neural tube, dorsal root ganglia, and limb buds in two lizard species. Our study emphasizes the value of revisiting Hox gene repertoires by densely sampling taxonomic groups and its feasibility owing to growing sequence resources in evaluating gene repertoires across taxa.

Diversifications often proceed along highly conserved, evolutionary trajectories. These patterns of covariation arise in ontogeny, which raises the possibility that adaptive morphologies are biased towards trait covariations that resemble growth trajectories. Here, we test this prediction in the diverse clade of Anolis lizards by investigating the covariation of embryonic growth of 13 fore- and hindlimb bones in 15 species, and compare these to the evolutionary covariation of these limb bones across 267 Anolis species. Our results demonstrate that species differences in relative limb length are established already at hatching, and are resulting from both differential growth and differential sizes of cartilaginous anlagen. Multivariate analysis revealed that Antillean Anolis share a common ontogenetic allometry that is characterized by positive allometric growth of the long bones relative to metapodial and phalangeal bones. This major axis of ontogenetic allometry in limb bones deviated from the major axis of evolutionary allometry of the Antillean Anolis and the two clades of mainland Anolis lizards. These results demonstrate that the remarkable diversification of locomotor specialists in Anolis lizards are accessible through changes that are largely independent from ontogenetic growth trajectories, and therefore likely to be the result of modifications that manifest at the earliest stages of limb development.

The Mediterranean Basin has experienced extensive change in geology and climate over the past six million years. Yet, the relative importance of key geological events for the distribution and genetic structure of the Mediterranean fauna remains poorly understood. Here, we use population genomic and phylogenomic analyses to establish the evolutionary history and genetic structure of common wall lizards (Podarcis muralis). This species is particularly informative because, in contrast to other Mediterranean lizards, it is widespread across the Iberian, Italian, and Balkan peninsulas, and in extra-Mediterranean regions. We found strong support for six major lineages within P. muralis, which were largely discordant with the phylogenetic relationship of mitochondrial DNA. The most recent common ancestor of extant P. muralis was likely distributed in the Italian Peninsula, and experienced an “Out-of-Italy” expansion following the Messinian salinity crisis (~5 Mya), resulting in the differentiation into the extant lineages on the Iberian, Italian and Balkan peninsulas. Introgression analysis revealed that both inter- and intraspecific gene flow have been pervasive throughout the evolutionary history of P. muralis. For example, the Southern Italy lineage has a hybrid origin, formed through admixture between the Central Italy lineage and an ancient lineage that was the sister to all other P. muralis. More recent genetic differentiation is associated with the onset of the Quaternary glaciations, which influenced population dynamics and genetic diversity of contemporary lineages. These results demonstrate the pervasive role of Mediterranean geology and climate for the evolutionary history and population genetic structure of extant species.

The survival prospects of threatened species or populations can sometimes be improved by adaptive change. Such evolutionary rescue is particularly relevant when the threat comes from changing environments, or when long-term population persistence requires range expansion into new habitats. Conservation biologists are therefore often interested in whether or not populations or lineages show a disposition for adaptive evolution, that is, if they are evolvable. Here, we discuss four alternative perspectives that target different causes of evolvability and outline some of the key challenges those perspectives are designed to address. Standing genetic variation provides one familiar estimate of evolvability. Yet, the mere presence of genetic variation is often insufficient to predict if a population will adapt, or how it will adapt. The reason is that adaptive change not only depends on genetic variation, but also on the extent to which this genetic variation can be realized as adaptive phenotypic variation. This requires attention to developmental systems and how plasticity influences evolutionary potential. Finally, we discuss how a better understanding of the different factors that contribute to evolvability can be exploited in conservation practice.

Anolis lizards originated in continental America but have colonized the Greater Antillean islands and recolonized the mainland, resulting in three major groups (Primary and Secondary Mainland and Greater Antillean). The adaptive radiation in the Greater Antilles has famously resulted in the repeated evolution of ecomorphs. Yet, it remains poorly understood to what extent this island radiation differs from diversification on the mainland. Here, we demonstrate that the evolutionary modularity between girdles and limbs is fundamentally different in the Greater Antillean and Primary Mainland Anolis. This is consistent with ecological opportunities on islands driving the adaptive radiation along distinct evolutionary trajectories. However, Greater Antillean Anolis share evolutionary modularity with the group that recolonized the mainland, demonstrating a persistent phylogenetic inertia. A comparison of these two groups support an increased morphological diversity and faster and more variable evolutionary rates on islands. These macroevolutionary trends of the locomotor skeleton in Anolis illustrate that ecological opportunities on islands can have lasting effects on morphological diversification., Both ecological opportunity and phenotypic modularity have been suggested to facilitate adaptive radiations. Feiner et al. show that Anolis lizards evolved a new modularity structure in their island adaptive radiation, but that this modularity did not produce the same extreme diversification when Anolis returned to the mainland.

International audience; Parasites have evolved a diversity of lifestyles that exploit the biology of their hosts. Some nematodes that parasitize mammals pass via the placenta or milk from one host to another. Similar cases of vertical transmission have never been reported in avian and nonavian reptiles, suggesting that egg laying may constrain the means of parasite transmission. However, here we report the first incidence of transovarial transmission of a previously undescribed nematode in an egg-laying amniote, the common wall lizard (Podarcis muralis). Nematodes enter the developing brain from the female ovary early in embryonic development. Infected lizard embryos develop normally and hatch with nematodes residing in their braincase. We present a morphological and molecular phylogenetic characterization of the nematode and suggest that particular features of lizard biology that are absent from birds and turtles facilitated the evolutionary origin of this novel life history.

Many features of the development of reptiles are affected by temperature, but very little is known about how incubation temperature affects gene expression. Here, we provide a detailed case study of gene expression profiles in common wall lizard (Podarcis muralis) embryos developing at stressfully low (15°C) versus benign (24°C) temperature. For maximum comparability between the two temperature regimes, we selected a precise developmental stage early in embryogenesis defined by the number of somites. We used a split-clutch design and lizards from four different populations to evaluate the robustness of temperature-responsive gene expression profiles. Embryos incubated at stressfully low incubation temperature expressed on average 20% less total RNA than those incubated at benign temperatures, presumably reflecting lower rates of transcription at cool temperature. After normalizing for differences in total amounts of input RNA, we find that approximately 50% of all transcripts show significant expression differences between the two incubation temperatures. Transcripts with the most extreme changes in expression profiles are associated with transcriptional and translational regulation and chromatin remodeling, suggesting possible epigenetic mechanisms underlying acclimation of early embryos to cool temperature. We discuss our findings in light of current advances in the use of transcriptomic data to study how individuals acclimatize and populations adapt to thermal stress.

Populations adapting independently to the same environment provide important insights into the repeatability of evolution at different levels of biological organization. In the 20th century, common wall lizards (Podarcis muralis) from southern and western Europe were introduced to England, north of their native range. Nonnative populations of both lineages have adapted to the shorter season and lower egg incubation temperature by increasing the absolute rate of embryonic development. Here, we tested if this adaptation is accompanied by signatures of directional selection in the transcriptomes of early embryos and, if so, if nonnative populations show adaptive convergence. Embryos from nonnative populations exhibited gene expression profiles consistent with directional selection following introduction, but different genes were affected in the two lineages. Despite this, the functional enrichment of genes that changed their expression following introduction showed substantial similarity between lineages, and was consistent with mechanisms that should promote developmental rate. Moreover, the divergence between nonnative and native populations was enriched for genes that were temperature-responsive in native populations. These results indicate that small populations are able to adapt to new climatic regimes, but the means by which they do so may largely be determined by founder effects and other sources of genetic drift.

Hybrid zones provide insights into the evolution of reproductive isolation. Sexual selection can contribute to the evolution of reproductive barriers, but it remains poorly understood how sexual traits impact gene flow in secondary contact. Here, we show that a recently evolved suite of sexual traits that function in male-male competition mediates gene flow between two lineages of wall lizards (Podarcis muralis). Gene flow was relatively low and asymmetric in the presence of exaggerated male morphology and coloration compared to when the lineages share the ancestral phenotype. Putative barrier loci were enriched in genomic regions that were highly differentiated between the two lineages and showed low concordance between the transects. The exception was a consistently low genetic exchange around ATXN1, a gene that modulates social behavior. We suggest that this gene may contribute to the male mate preferences that are known to cause lineage-assortative mating in this species. Although female choice modulates the degree of reproductive isolation in a variety of taxa, wall lizards demonstrate that both male-male competition and male mate choice can contribute to the extent of gene flow between lineages.

The radiation of Anolis lizards in the Caribbean is associated with a diversification of the functional match between morphology, habitat use, and locomotor performance. It has been hypothesized that the microhabitat a lizard is reared in can achieve a similar fit of form and function within a species. This predicts that plasticity in the locomotor apparatus is accompanied by changes in perching behavior or improved locomotor performance. To test this, we raised juveniles of two species (Anolis sagrei and Anolis carolinensis) on either broad or narrow surfaces and examined perching behavior and locomotor performance as well as the shape of the pectoral and pelvic girdles, limb length, and thickness of the long bones. Perching behavior was not affected by the habitat surface experienced during ontogeny. However, individuals raised on broad surfaces showed better locomotor performance on broad surfaces, and the magnitude of the effect was as large as the difference between the two species. Both species showed modifications of pectoral and pelvic shape, but only A. carolinensis developed longer limbs on broad surfaces. However, these morphological adjustments induced by physical activity did not explain why lizards raised on broad surfaces performed better. Thus, it appears that early-life experiences can affect both the morphology of the locomotor apparatus and locomotor performance in Anolis lizards, without the two being functionally connected.

Developmental plasticity looks like a promising bridge between ecological and developmental perspectives on evolution. Yet, there is no consensus on whether plasticity is part of the explanation for adaptive evolution or an optional "add-on" to genes and natural selection. Here, we suggest that these differences in opinion are caused by differences in the simplifying assumptions, and particular idealizations, that enable evolutionary explanation. We outline why idealizations designed to explain evolution through natural selection prevent an understanding of the role of development, and vice versa. We show that representing plasticity as a reaction norm conforms with the idealizations of selective explanations, which can give the false impression that plasticity has no explanatory power for adaptive evolution. Finally, we use examples to illustrate why evolutionary explanations that include developmental plasticity may in fact be more satisfactory than explanations that solely refer to genes and natural selection.

Hox genes orchestrate development by patterning the embryonic axis. Vertebrate Hox genes are arranged in four compact clusters, and the spacing between genes is assumed to be crucial for their function. The genomes of squamate reptiles are unusually rich and variable in transposable elements (TEs), and it has been suggested that TE invasion is responsible for the Hox cluster expansion seen in snakes and lizards. Using de novo TE prediction on 17 genomes of lizards and snakes, I show that TE content of Hox clusters are generally 50% lower than genome‐wide TE levels. However, two distantly related lizards of the species‐rich genus Anolis have Hox clusters with a TE content that approaches genomic levels. The age distribution of TEs in Anolis lizards revealed that peaks of TE activity broadly coincide with speciation events. In accordance with theoretical models of Hox cluster regulation, I find that Anolis species with many TEs in their Hox clusters show aberrant Hox gene expression patterns, suggesting a functional link between TE accumulation and embryonic development. These results are consistent with the hypothesis that TEs play a role in developmental processes as well as in evolutionary diversifications.

The vertebrate gene repertoire is characterized by “cryptic” genes whose identification has been hampered by their absence from the genomes of well-studied species. One example is the Bmp16 gene, a paralog of the developmental key genes Bmp2 and -4. We focus on the Bmp2/4/16 group of genes to study the evolutionary dynamics following gen(om)e duplications with special emphasis on the poorly studied Bmp16 gene. We reveal the presence of Bmp16 in chondrichthyans in addition to previously reported teleost fishes and reptiles. Using comprehensive, vertebrate-wide gene sampling, our phylogenetic analysis complemented with synteny analyses suggests that Bmp2, -4 and -16 are remnants of a gene quartet that originated during the two rounds of whole-genome duplication (2R-WGD) early in vertebrate evolution. We confirm that Bmp16 genes were lost independently in at least three lineages (mammals, archelosaurs and amphibians) and report that they have elevated rates of sequence evolution. This finding agrees with their more “flexible” deployment during development; while Bmp16 has limited embryonic expression domains in the cloudy catshark, it is broadly expressed in the green anole lizard. Our study illustrates the dynamics of gene family evolution by integrating insights from sequence diversification, gene repertoire changes, and shuffling of expression domains.

Reptiles use pterin and carotenoid pigments to produce yellow, orange, and red colors. These conspicuous colors serve a diversity of signaling functions, but their molecular basis remains unresolved. Here, we show that the genomes of sympatric color morphs of the European common wall lizard, which differ in orange and yellow pigmentation and in their ecology and behavior, are virtually undifferentiated. Genetic differences are restricted to two small regulatory regions, near genes associated with pterin (SPR) and carotenoid metabolism (BCO2), demonstrating that a core gene in the housekeeping pathway of pterin biosynthesis has been co-opted for bright coloration in reptiles and indicating that these loci exert pleiotropic effects on other aspects of physiology. Pigmentation differences are explained by extremely divergent alleles and haplotype analysis revealed abundant trans-specific allele sharing with other lacertids exhibiting color polymorphisms. The evolution of these conspicuous color ornaments is the result of ancient genetic variation and cross-species hybridization.

Phylogenetic approaches are indispensable in any comparative molecular study involving multiple species. These approaches are in increasing demand as the amount and availability of DNA sequence information continues to increase exponentially, even for organisms that were previously not extensively studied. Without the sound application of phylogenetic concepts and knowledge, one can be misled when attempting to infer ancestral character states as well as the timing and order of evolutionary events, both of which are frequently exerted in evolutionary developmental biology. The ignorance of phylogenetic approaches can also impact non-evolutionary studies and cause misidentification of the target gene or protein to be examined in functional characterization. This review aims to promote tree-thinking in evolutionary conjecture and stress the importance of a sense of time scale in cross-species comparisons, in order to enhance the understanding of phylogenetics in all biological fields including developmental biology. To this end, molecular phylogenies of several developmental regulatory genes, including those denoted as "cryptic pan-vertebrate genes", are introduced as examples.

Acquisition and storage of Latimeria chalumnae samples was supported by grants from the African Coelacanth Ecosystem Programme of the South African National Department of Science and Technology. Generation of the Latimeria chalumnae and Protopterus annectens sequences by the Broad Institute of the Massachusetts Institute of Technology (MIT) and Harvard University was supported by grants from the National Human Genome Research Institute (NHGRI). K.L.T. is the recipient of a EURYI award from the European Science Foundation. We would also like to thank the Genomics Sequencing Platform of the Broad Institute for sequencing the L. chalumnae genome and L. chalumnae and P. annectens transcriptomes, S. Ahamada, R. Stobbs and the Association pour le Protection de Gombesa (APG) for their help in obtaining coelacanth samples, Y. Zhao for the use of data from Rana chensinensis, and L. Gaffney, C. Hamilton and J. Westlund for assistance with figure preparation. 10; International audience; The discovery of a living coelacanth specimen in 1938 was remarkable, as this lineage of lobe-finned fish was thought to have become extinct 70 million years ago. The modern coelacanth looks remarkably similar to many of its ancient relatives, and its evolutionary proximity to our own fish ancestors provides a glimpse of the fish that first walked on land. Here we report the genome sequence of the African coelacanth, Latimeria chalumnae. Through a phylogenomic analysis, we conclude that the lungfish, and not the coelacanth, is the closest living relative of tetrapods. Coelacanth protein-coding genes are significantly more slowly evolving than those of tetrapods, unlike other genomic features. Analyses of changes in genes and regulatory elements during the vertebrate adaptation to land highlight genes involved in immunity, nitrogen excretion and the development of fins, tail, ear, eye, brain and olfaction. Functional assays of enhancers involved in the fin-to-limb transition and in the emergence of extra-embryonic tissues show the importance of the coelacanth genome as a blueprint for understanding tetrapod evolution.

Populations adapting independently to the same environment provide important insights into the repeatability of evolution at different levels of biological organization. In the 20th century, common wall lizards (Podarcis muralis) from southern and western Europe were introduced to England, north of their native range. Nonnative populations of both lineages have adapted to the shorter season and lower egg incubation temperature by increasing the absolute rate of embryonic development. Here, we tested if this adaptation is accompanied by signatures of directional selection in the transcriptomes of early embryos and, if so, if nonnative populations show adaptive convergence. Embryos from nonnative populations exhibited gene expression profiles consistent with directional selection following introduction, but different genes were affected in the two lineages. Despite this, the functional enrichment of genes that changed their expression following introduction showed substantial similarity between lineages, and was consistent with mechanisms that should promote developmental rate. Moreover, the divergence between nonnative and native populations was enriched for genes that were temperature-responsive in native populations. These results indicate that small populations are able to adapt to new climatic regimes, but the means by which they do so may largely be determined by founder effects and other sources of genetic drift.

We have recently shown that the human Nuclear pore-associated protein (NPAP1)/C15orf2 gene encodes a nuclear pore-associated protein. This gene is one of several paternally expressed imprinted genes in the genomic region 15q11q13. Because the Prader-Willi syndrome is known to be caused by the loss of function of paternally expressed genes in 15q11q13, a phenotypic contribution of NPAP1 cannot be excluded. NPAP1 appears to be under strong positive Darwinian selection in primates, suggesting an important function in primatebiology. Interestingly, however, incontrast to all other protein-coding genes in15q11q13, NPAP1 has no ortholog in the mouse. Our investigation of the evolutionary origin of NPAP1 showed that the gene is specific to primate species and absent from the 15q11q13-orthologous regions in all non primatemammals. However, we identified a group of paralogousgenes, which we call NPAP1L, in all placental mammals except rodents. Phylogenetic analysis revealed that NPAP1, NPAP1L, and another group of genes (UPF0607), which is also restricted to primates, are closely related to the vertebrate transmembrane nucleoporin gene POM121, although they lack the transmembrane domain. These three newly identified groups of genes all lack conserved introns, and hence, are likely retrogenes. We hypothesize that, in the common ancestor of placentals, the POM121 gene retro transposed and gave rise to an NPAP1-ancestral retrogene NPAP1L/NPAP1/UPF0607. Our results suggest that the nuclear pore-associated gene NPAP1 originates from the vertebrate nucleoporingene POM121 and - after several steps of retrotransposition and duplication - has been subjected to genomic imprinting and positive selection after integration into the imprinted SNRPN-UBE3A chromosomal domain. \copyright The Author(s) 2014.

Transposable elements (TEs) are DNA sequences that can insert elsewhere in the genome and modify genome structure and gene regulation. The role of TEs in evolution is contentious. One hypothesis posits that TE activity generates genomic incompatibilities that can cause reproductive isolation between incipient species. This predicts that TEs will accumulate during speciation events. Here, I tested the prediction that extant lineages with a relatively high rate of speciation have a high number of TEs in their genomes. I sequenced and analysed the TE content of a marker genomic region (Hoxclusters) inAnolislizards, a classic case of an adaptive radiation. Unlike other vertebrates, including closely related lizards,Anolislizards have high numbers of TEs in theirHoxclusters, genomic regions that regulate development of the morphological adaptations that characterize habitat specialists in these lizards. Following a burst of TE activity in the lineage leading to extantAnolis, TEs have continued to accumulate during or after speciation events, resulting in a positive relationship between TE density and lineage speciation rate. These results are consistent with the prediction that TE activity contributes to adaptive radiation by promoting speciation. Although there was no evidence that TE densityper seis associated with ecological morphology, the activity of TEs inHoxclusters could have been a rich source for phenotypic variation that may have facilitated the rapid parallel morphological adaptation to microhabitats seen in extantAnolislizards.

Bilaterian Hox genes play pivotal roles in the specification of positional identities along the anteroposterior axis. Particularly in vertebrates, their regulation is tightly coordinated by tandem arrays of genes [paralogy groups (PGs)] in four gene clusters (HoxA-D). Traditionally, the uninterrupted Hox cluster (Hox1-14) of the invertebrate chordate amphioxus was regarded as an archetype of the vertebrate Hox clusters. In contrast to Hox1-13 that are globally regulated by the "Hox code" and are often phylogenetically conserved, vertebrate Hox14 members were only recently revealed to be present in an African lungfish, a coelacanth, chondrichthyans and a lamprey, and decoupled from the Hox code. In this study we performed a PCR-based search of Hox14 members from diverse vertebrates, and identified one in the Australian lungfish, Neoceratodus forsteri. Based on a molecular phylogenetic analysis, this gene was designated NfHoxA14. Our real-time RT-PCR suggested its hindgut-associated expression, previously observed also in cloudy catshark HoxD14 and lamprey Hox14Alpha. It is likely that this altered expression scheme was established before the Hox cluster quadruplication, probably at the base of extant vertebrates. To investigate the origin of vertebrate Hox14, by including this sarcopterygian Hox14 member, we performed focused phylogenetic analyses on its relationship with other vertebrate posterior Hox PGs (Hox9-13) as well as amphioxus posterior Hox genes. Our results confirmed the hypotheses previously proposed by other studies that vertebrate Hox14 does not have any amphioxus ortholog, and that none of 1-to-1 pairs of vertebrate and amphioxus posterior Hox genes, based on their relative location in the clusters, is orthologous.

Background Whole genome sequences have allowed us to have an overview of the evolution of gene repertoires. The target of the present study, the TGFβ superfamily, contains many genes involved in vertebrate development, and provides an ideal system to explore the relationships between evolution of gene repertoires and that of developmental programs. Results As a result of a bioinformatic survey of sequenced vertebrate genomes, we identified an uncharacterized member of the TGFβ superfamily, designated bmp16, which is confined to teleost fish species. Our molecular phylogenetic study revealed a high affinity of bmp16 to the Bmp2/4 subfamily. Importantly, further analyses based on the maximum-likelihood method unambiguously ruled out the possibility that this teleost-specific gene is a product of teleost-specific genome duplication. This suggests that the absence of a bmp16 ortholog in tetrapods is due to a secondary loss. In situ hybridization showed embryonic expression of the zebrafish bmp16 in the developing swim bladder, heart, tail bud, and ectoderm of pectoral and median fin folds in pharyngula stages, as well as gut-associated expression in 5-day embryos. Conclusion Comparisons of expression patterns revealed (1) the redundancy of bmp16 expression with its homologs in presumably plesiomorphic expression domains, such as the fin fold, heart, and tail bud, which might have permitted its loss in the tetrapod lineage, and (2) the loss of craniofacial expression and gain of swim bladder expression of bmp16 after the gene duplication between Bmp2, -4 and -16. Our findings highlight the importance of documenting secondary changes of gene repertoires and expression patterns in other gene families.