Search

SummaryTraits shared among distantly related lineages are indicators of common evolutionary constraints, at the ecological, physiological or molecular level. The vertebral stripe is a color pattern that is widespread across the anuran phylogeny. Despite its prevalence in the order, surprisingly little is known about the genetic basis and evolutionary dynamic of this color pattern. Here we combine histology, genome- and transcriptome-wide analyses with order-scale phylogenetic comparative analyses to investigate this common phenotype. We show that the vertebral stripe has evolved hundreds of times in the evolutionary history of anurans and is selected for in terrestrial habitats. Using the Ethiopian Ptychadena radiation as a model system, we demonstrate that variation at the ASIP gene is responsible for the different vertebral stripe phenotypes. Alleles associated to these phenotypes are younger than the split between closely related Ptychadena species, thus indicating that the vertebral stripe results from parallel evolution within the group. Our findings demonstrate that this widespread color pattern evolves rapidly and recurrently in terrestrial anurans, and therefore constitute an ideal system to study repeated evolution.

Gaining a better understanding on how selection and neutral processes affect genomic diversity is essential to gain better insights into the mechanisms driving adaptation and speciation. However, the evolutionary processes affecting variation at a genomic scale have not been investigated in most vertebrate lineages. Previous studies have been limited to a small number of model species, mostly mammals, and no studies have investigated genomic variation in non avian reptiles. Here we present the first population genomics survey using whole genome re sequencing in the green anole (Anolis carolinensis). This species has emerged as a model for the study of genomic evolution in squamates. We quantified how demography, recombination and selection have led to the current genetic diversity of the green anole by using whole-genome resequencing of five genetic clusters covering the entire species range. The differentiation of green anole’s populations is consistent with a northward expansion from South Florida followed by genetic isolation and subsequent gene flow among adjacent genetic clusters. Dispersal out-of-Florida was accompanied by a drastic population bottleneck followed by a rapid population expansion. This event was accompanied by male-biased dispersal and/or selective sweeps on the X chromosome. We show that the combined effect of background selection and recombination rates is the main contributor to the genomic landscape of differentiation in the anole genome. We further demonstrate that recombination rates are positively correlated with GC content at third codon position (GC3) and confirm the importance of biased gene conversion in shaping genome wide patterns of diversity in reptiles.

We investigated the evolution of the families of LINE-1 (L1) retrotransposons that have amplified in the human lineage since the origin of primates. We identified two phases in the evolution of L1. From approximately 70 million years ago (Mya) until approximately 40 Mya, three distinct L1 lineages were simultaneously active in the genome of ancestral primates. In contrast, during the last 40 million years (Myr), i.e., during the evolution of anthropoid primates, a single lineage of families has evolved and amplified. We found that novel (i.e., unrelated) regulatory regions (5'UTR) have been frequently recruited during the evolution of L1, whereas the two open-reading frames (ORF1 and ORF2) have remained relatively conserved. We found that L1 families coexisted and formed independently evolving L1 lineages only when they had different 5'UTRs. We propose that L1 families with different 5'UTR can coexist because they don't rely on the same host-encoded factors for their transcription and therefore do not compete with each other. The most prolific L1 families (families L1PA8 to L1PA3) amplified between 40 and 12 Mya. This period of high activity corresponds to an episode of adaptive evolution in a segment of ORF1. The correlation between the high activity of L1 families and adaptive evolution could result from the coevolution of L1 and a host-encoded repressor of L1 activity.