Genome of the early spider-orchid Ophrys sphegodes provides insights into sexual deception and pollinator adaptation.

Pollinator-driven evolution of floral traits is thought to be a major driver of angiosperm speciation and diversification. Ophrys orchids mimic female insects to lure male pollinators into pseudocopulation. This strategy, called sexual deception, is species-specific, thereby providing strong premating reproductive isolation. Identifying the genomic architecture underlying pollinator adaptation and speciation may shed light on the mechanisms of angiosperm diversification. Here, we report the 5.2 Gb chromosome-scale genome sequence of Ophrys sphegodes. We find evidence for transposable element expansion that preceded the radiation of the O. sphegodes group, and for gene duplication having contributed to the evolution of chemical mimicry. We report a highly differentiated genomic candidate region for pollinator-mediated evolution on chromosome 2. The Ophrys genome will prove useful for investigations into the repeated evolution of sexual deception, pollinator adaptation and the genomic architectures that facilitate evolutionary radiations.

Gene flow throughout the evolutionary history of a colour polymorphic and generalist clownfish.

Even seemingly homogeneous on the surface, the oceans display high environmental heterogeneity across space and time. Indeed, different soft barriers structure the marine environment, which offers an appealing opportunity to study various evolutionary processes such as population differentiation and speciation. Here, we focus on Amphiprion clarkii (Actinopterygii; Perciformes), the most widespread of clownfishes that exhibits the highest colour polymorphism. Clownfishes can only disperse during a short pelagic larval phase before their sedentary adult lifestyle, which might limit connectivity among populations, thus facilitating speciation events. Consequently, the taxonomic status of A. clarkii has been under debate. We used whole-genome resequencing data of 67 A. clarkii specimens spread across the Indian and Pacific Oceans to characterize the species' population structure, demographic history and colour polymorphism. We found that A. clarkii spread from the Indo-Pacific Ocean to the Pacific and Indian Oceans following a stepping-stone dispersal and that gene flow was pervasive throughout its demographic history. Interestingly, colour patterns differed noticeably among the Indonesian populations and the two populations at the extreme of the sampling distribution (i.e. Maldives and New Caledonia), which exhibited more comparable colour patterns despite their geographic and genetic distances. Our study emphasizes how whole-genome studies can uncover the intricate evolutionary past of wide-ranging species with diverse phenotypes, shedding light on the complex nature of the species concept paradigm.

Recurrent gene flow events occurred during the diversification of clownfishes of the skunk complex.

Clownfish (subfamily Amphiprioninae) are an iconic group of coral reef fish that evolved a mutualistic interaction with sea anemones, which triggered the adaptive radiation of the clade. Within clownfishes, the "skunk complex" is particularly interesting. Besides ecological speciation, interspecific gene flow and hybrid speciation are thought to have shaped the evolution of the group. We investigated the mechanisms characterizing the diversification of this complex. By taking advantage of their disjunct geographical distribution, we obtained whole-genome data of sympatric and allopatric populations of the three main species of the complex (Amphiprion akallopisos, A. perideraion and A. sandaracinos). We examined population structure, genomic divergence and introgression signals and performed demographic modelling to identify the most realistic diversification scenario. We excluded scenarios of strict isolation or hybrid origin of A. sandaracinos. We discovered moderate gene flow from A. perideraion to the ancestor of A. akallopisos + A. sandaracinos and weak gene flow between the species in the Indo-Australian Archipelago throughout the diversification of the group. We identified introgressed regions in A. sandaracinos and detected in A. perideraion two large regions of high divergence from the two other species. While we found that gene flow has occurred throughout the species' diversification, we also observed that recent admixture was less pervasive than initially thought, suggesting a role of host repartition or behavioural barriers in maintaining the genetic identity of the species in sympatry.

Characterization of the complete plastome of Delphinium montanum, a polyploid, endemic and endangered Pyrenean Larkspur.

Delphinium montanum DC. 1815, is an endangered larkspur endemic to the Eastern Pyrenees. For biogeographic and conservation purpose, a hybrid assembly approach based on long- and short-read genomic data allowed us to successfully assemble whole plastid genome of Delphinium montanum. The complete plastome is 154,185 bp in length, consisting of a pair of inverted repeats (IRs) of 26,559 bp, a large single-copy (LSC) region and a small single-copy region (SSC) of 84,746 and 16,320 bp, respectively. It was found to contain 136 genes, including 84 protein-coding genes, 44 trRNA genes and 8 rRNA genes. The overall GC content of the plastid genome is 38.3%. Phylogenetic inference supports the polyphyly of the Delphinium genus.

Little hope for the polyploid endemic Pyrenean Larkspur (Delphinium montanum): Evidences from population genomics and Ecological Niche Modeling.

Species endemic to restricted geographical ranges represent a particular conservation issue, be it for their heritage interest. In a context of global change, this is particularly the case for plants which belong to high-mountain ecosystems and, because of their ecological requirements, are doomed to survive or disappear on their "sky islands". The Pyrenean Larkspur (Delphinium montanum, Ranunculaceae) is endemic to the Eastern part of the Pyrenees (France and Spain). It is now only observable at a dozen of localities and some populations show signs of decline, such as a recurrent lack of flowering. Implementing population genomics approach (e.g., RAD-seq like) is particularly useful to understand genomic patterns of diversity and differentiation in order to provide recommendations in term of conservation. However, it remains challenging for species such as D. montanum that are autotetraploid with a large genome size (1C-value >10 pg) as most methods currently available were developed for diploid species. A Bayesian framework able to call genotypes with uncertainty allowed us to assess genetic diversity and population structure in this system. Our results show evidence for inbreeding (mean G IS = 0.361) within all the populations and substantial population structure (mean G ST = 0.403) at the metapopulation level. In addition to a lack of connectivity between populations, spatial projections of Ecological Niche Modeling (ENM) analyses under different climatic scenarios predict a dramatic decrease of suitable habitat for D. montanum in the future. Based on these results, we discuss the relevance and feasibility of different conservation measures.

Why are there so many bee-orchid species? Adaptive radiation by intra-specific competition for mnesic pollinators.

Adaptive radiations occur mostly in response to environmental variation through the evolution of key innovations that allow emerging species to occupy new ecological niches. Such biological innovations may play a major role in niche divergence when emerging species are engaged in reciprocal ecological interactions. To demonstrate coevolution is a difficult task; only a few studies have confirmed coevolution as driver of speciation and diversification. Herein we review current knowledge about bee orchid (Ophrys spp.) reproductive biology. We propose that the adaptive radiation of the Mediterranean orchid genus Ophrys, comprising several hundred species, is due to coevolutionary dynamics between these plants and their pollinators. We suggest that pollination by sexual swindling used by Ophrys orchids is the main driver of this coevolution. Flowers of each Ophrys species mimic a sexually receptive female of one particular insect species, mainly bees. Male bees are first attracted by pseudo-pheromones emitted by Ophrys flowers that are similar to the sexual pheromones of their females. Males then are lured by the flower shape, colour and hairiness, and attempt to copulate with the flower, which glues pollen onto their bodies. Pollen is later transferred to the stigma of another flower of the same Ophrys species during similar copulation attempts. In contrast to rewarding pollination strategies, Ophrys pollinators appear to be parasitized. Here we propose that this apparent parasitism is in fact a coevolutionary relationship between Ophrys and their pollinators. For plants, pollination by sexual swindling could ensure pollination efficiency and specificity, and gene flow among populations. For pollinators, pollination by sexual swindling could allow habitat matching and inbreeding avoidance. Pollinators might use the pseudo-pheromones emitted by Ophrys to locate suitable habitats from a distance within complex landscapes. In small populations, male pollinators would disperse once they have memorized the local diversity of sexual pseudo-pheromone bouquets or if all Ophrys flowers are fertilized and thus repel pollinators via production of repulsive pheromones that mimic those produced by fertilized female bees. We propose the following evolutionary scenario: Ophrys radiation is driven by strong intra-specific competition among Ophrys individuals for the attraction of species-specific pollinators, which is a consequence of the high cognitive abilities of pollinators. Male bees record the pheromone signatures of kin or of previously courted partners to avoid further copulation attempts, thereby inducing strong selection on Ophrys for variation in odour bouquets emitted by individual flowers. The resulting odour bouquets could by chance correspond to pseudo-pheromones of the females of another bee species, and thus attract a new pollinator. If such pollinator shifts occur simultaneously in several indivuals, pollen exchanges might occur and initiate speciation. To reinforce the attraction of the new pollinator and secure prezygotic isolation, the following step is directional selection on flower phenotypes (shape, colour and hairiness) towards a better match with the body of the pollinator's female. Pollinator shift and the resulting prezygotic isolation is adaptive for new Ophrys species because they may benefit from competitor-free space for limited pollinators. We end our review by proritizing several critical research avenues.

Differential divergence in autosomes and sex chromosomes is associated with intra-island diversification at a very small spatial scale in a songbird lineage.

Recently diverged taxa showing marked phenotypic and ecological diversity provide optimal systems to understand the genetic processes underlying speciation. We used genome-wide markers to investigate the diversification of the Reunion grey white-eye (Zosterops borbonicus) on the small volcanic island of Reunion (Mascarene archipelago), where this species complex exhibits four geographical forms that are parapatrically distributed across the island and differ strikingly in plumage colour. One form restricted to the highlands is separated by a steep ecological gradient from three distinct lowland forms which meet at narrow hybrid zones that are not associated with environmental variables. Analyses of genomic variation based on single nucleotide polymorphism data from genotyping-by-sequencing and pooled RAD-seq approaches show that signatures of selection associated with elevation can be found at multiple regions across the genome, whereas most loci associated with the lowland forms are located on the Z sex chromosome. We identified TYRP1, a Z-linked colour gene, as a likely candidate locus underlying colour variation among lowland forms. Tests of demographic models revealed that highland and lowland forms diverged in the presence of gene flow, and divergence has progressed as gene flow was restricted by selection at loci across the genome. This system holds promise for investigating how adaptation and reproductive isolation shape the genomic landscape of divergence at multiple stages of the speciation process.

Characterization of the complete plastome of Ophrys aveyronensis, a Euro-Mediterranean orchid with an intriguing disjunct geographic distribution.

Ophrys aveyronensis is an orchid with disjunct geographic distribution. For biogeographic and conservation purpose, we sequenced its complete plastome using Illumina data. The complete plastome is 146,816 bp in length, consisting of a pair of inverted repeats (IRs) of 25,006 bp, a large single-copy (LSC) region and a small single-copy region (SSC) of 80,495 and 16,309 bp, respectively. It was found to contain 133 genes, including 86 protein-coding genes, up to 39 trRNA genes and 8 rRNA genes. The overall GC content of the plastid genome is 36.9%. Phylogenetic inference confirms that O. aveyronensis is very close to O. sphegodes.

First draft genome of an iconic clownfish species (Amphiprion frenatus).

Clownfishes (or anemonefishes) form an iconic group of coral reef fishes, principally known for their mutualistic interaction with sea anemones. They are characterized by particular life history traits, such as a complex social structure and mating system involving sequential hermaphroditism, coupled with an exceptionally long lifespan. Additionally, clownfishes are considered to be one of the rare groups to have experienced an adaptive radiation in the marine environment. Here, we assembled and annotated the first genome of a clownfish species, the tomato clownfish (Amphiprion frenatus). We obtained 17,801 assembled scaffolds, containing a total of 26,917 genes. The completeness of the assembly and annotation was satisfying, with 96.5% of the Actinopterygii Benchmarking Universal Single-Copy Orthologs (BUSCOs) being retrieved in A. frenatus assembly. The quality of the resulting assembly is comparable to other bony fish assemblies. This resource is valuable for advancing studies of the particular life history traits of clownfishes, as well as being useful for population genetic studies and the development of new phylogenetic markers. It will also open the way to comparative genomics. Indeed, future genomic comparison among closely related fishes may provide means to identify genes related to the unique adaptations to different sea anemone hosts, as well as better characterize the genomic signatures of an adaptive radiation.

Molecular systematics of the anchovy genus Encrasicholina in the Northwest Pacific.

The anchovy genus Encrasicholina is an important coastal marine resource of the tropical Indo-West Pacific (IWP) region for which insufficient comparative data are available to evaluate the effects of current exploitation levels on the sustainability of its species and populations. Encrasicholina currently comprises nine valid species that are morphologically very similar. Only three, Encrasicholina punctifer, E. heteroloba, and E. pseudoheteroloba, occur in the Northwest Pacific subregion of the northeastern part of the IWP region. These species are otherwise broadly distributed and abundant in the IWP region, making them the most important anchovy species for local fisheries. In this study, we reconstructed the phylogeny of these three species of Encrasicholina within the Engraulidae. We sequenced 10 complete mitochondrial genomes (using high-throughput and Sanger DNA sequencing technologies) and compared those sequences to 21 previously published mitochondrial genomes from various engraulid taxa. The phylogenetic results showed that the genus Encrasicholina is monophyletic, and it is the sister group to the more-diverse "New World anchovy" clade. The mitogenome-based dating results indicated that the crown group Encrasicholina originated about 33.7 million years ago (nearby the limit Eocene/Oligocene), and each species of Encrasicholina has been reproductively isolated from the others for more than 20 million years, despite their morphological similarities. In contrast, preliminary population genetic analyses across the Northwest Pacific region using four mitogenomic sequences revealed very low levels of genetic differentiation within Encrasicholina punctifer. These molecular results combined with recent taxonomic revisions are important for designing further studies on the population structure and phylogeography of these anchovies.

Phylogenetic position of the rainbow sardine Dussumieria (Dussumieriidae) and its bearing on the early evolution of the Clupeoidei.

The fish family Dussumieriidae (suborder Clupeoidei), commonly called round herrings, is traditionally considered to be a key taxon for understanding the evolution of the Clupeoidei because some of its morphological characteristics have been interpreted as being either derived or primitive, such as the nearly complete absence of abdominal scutes. Recent molecule-based studies showed that the Dussumieriidae is likely not a monophyletic group. None of those studies, however, included the genus Dussumieria (rainbow sardines) which is the type genus of the family Dussumieriidae. Herein, we investigated the phylogenetic position of Dussumieria within the Clupeoidei, using a dataset of complete mitogenomic sequences, including five newly determined using high-throughput sequencing technology. In the inferred phylogenetic reconstructions, the Dussumieriidae was never recovered as monophyletic, and Dussumieria was not exclusively related to any other of the three dussumieriid genera. Although the position of Dussumieria is not fully resolved, this genus represents a major evolutionary lineage within the Clupeoidei, along with the Engraulidae, Pristigasteridae, Clupeinae, Etrumeus, and two unnamed clades, one containing the Chirocentridae and Spratelloidinae and the other containing the Ehiravinae, Dorosomatinae, and Alosinae. Our results allow some comments regarding the early evolution of the Clupeoidei. In particular, they strongly support the hypothesis that the W-shaped pelvic scute does not represent a good phylogenetic character within the Clupeoidei as either it is primitive or, alternatively, it has independently evolved several times.

A novel locus on chromosome 1 underlies the evolution of a melanic plumage polymorphism in a wild songbird.

Understanding the mechanisms responsible for phenotypic diversification within and among species ultimately rests with linking naturally occurring mutations to functionally and ecologically significant traits. Colour polymorphisms are of great interest in this context because discrete colour patterns within a population are often controlled by just a few genes in a common environment. We investigated how and why phenotypic diversity arose and persists in the Zosterops borbonicus white-eye of Reunion (Mascarene archipelago), a colour polymorphic songbird in which all highland populations contain individuals belonging to either a brown or a grey plumage morph. Using extensive phenotypic and genomic data, we demonstrate that this melanin-based colour polymorphism is controlled by a single locus on chromosome 1 with two large-effect alleles, which was not previously described as affecting hair or feather colour. Differences between colour morphs appear to rely upon complex cis-regulatory variation that either prevents the synthesis of pheomelanin in grey feathers, or increases its production in brown ones. We used coalescent analyses to show that, from a 'brown' ancestral population, the dominant 'grey' allele spread quickly once it arose from a new mutation. Since colour morphs are always found in mixture, this implies that the selected allele does not go to fixation, but instead reaches an intermediate frequency, as would be expected under balancing selection.

Phylogeography of the sergeants Abudefduf sexfasciatus and A. vaigiensis reveals complex introgression patterns between two widespread and sympatric Indo-West Pacific reef fishes.

On evolutionary timescales, sea level oscillations lead to recurrent spatio-temporal variation in species distribution and population connectivity. In this situation, applying classical concepts of biogeography is challenging yet necessary to understand the mechanisms underlying biodiversity in highly diverse marine ecosystems such as coral reefs. We aimed at studying the outcomes of such complex biogeographic dynamics on reproductive isolation by sampling populations across a wide spatial range of a species-rich fish genus: the sergeants (Pomacentridae: Abudefduf). We generated a mutlilocus data set that included ten morpho-species from 32 Indo-West Pacific localities. We observed a pattern of mito-nuclear discordance in two common and widely distributed species: Abudefduf sexfasciatus and Abudefduf vaigiensis. The results showed three regional sublineages (Indian Ocean, Coral Triangle region, western Pacific) in A. sexfasciatus (0.6-1.5% divergence at cytb). The other species, A. vaigiensis, is polyphyletic and consists of three distinct genetic lineages (A, B and C) (9% divergence at cytb) whose geographic ranges overlap. Although A. vaigiensis A and A. sexfasciatus were found to be distinct based on nuclear information, A. vaigiensis A was found to be nested within A. sexfasciatus in the mitochondrial gene tree. A. sexfasciatus from the Coral Triangle region and A. vaigiensis A were not differentiated from each other at the mitochondrial locus. We then used coalescent-based simulation to characterize a spatially widespread but weak gene flow between the two species. We showed that these fishes are good candidates to investigate the evolutionary complexity of the discrepancies between phenotypic and genetic similarity in closely related species.

Candidate Gene Analysis Suggests Untapped Genetic Complexity in Melanin-Based Pigmentation in Birds.

Studies on melanin-based color variation in a context of natural selection have provided a wealth of information on the link between phenotypic and genetic variation. Here, we evaluated associations between melanic plumage patterns and genetic polymorphism in the Réunion grey white-eye (Zosterops borbonicus), a species in which mutations on MC1R do not seem to play any role in explaining melanic variation. This species exhibits 5 plumage color variants that can be grouped into 3 color forms which occupy discrete geographic regions in the lowlands of Réunion, and a fourth high-elevation form which comprises 2 color morphs (grey and brown) and represents a true color polymorphism. We conducted a comprehensive survey of sequence variation in 96 individuals at a series of 7 candidate genes other than MC1R that have been previously shown to influence melanin-based color patterns in vertebrates, including genes that have rarely been studied in a wild bird species before: POMC, Agouti, TYR, TYRP1, DCT, Corin, and SLC24A5 Of these 7 genes, 2 (Corin and TYRP1) displayed an interesting shift in allele frequencies between lowland and highland forms and a departure from mutation-drift equilibrium consistent with balancing selection in the polymorphic highland form only. Sequence variation at Agouti, a gene frequently involved in melanin-based pigmentation patterning, was not associated with color forms or morphs. Thus, we suggest that functionally important changes in loci other than those classically studied are involved in the color polymorphism exhibited by the Réunion grey white-eye and possibly many other nonmodel species.

A Novel Female-Specific and Sexual Reproduction-Associated Dmrt Gene Discovered in the Stony Coral, Euphyllia ancora.

Transcription factors encoded by the Dmrt gene family regulate multiple aspects of animal reproduction. Most studies investigating the Dmrt gene family were conducted in model organisms from bilateral species, with a particular emphasis on gene function in male sex determination. It is still unclear whether the E. ancora Dmrt (EaDmrt) genes found in basal metazoans such as cnidarians share similar characteristics with orthologs in other metazoans. In this study, seven full Dmrt gene transcript sequences for a gonochoric coral, Euphyllia ancora (phylum: Cnidaria; class: Anthozoa), were obtained through transcriptome data mining, RT-PCR analysis, rapid amplification of cDNA ends, and sequencing. These EaDmrts were subjected to quantitative assays measuring temporal and tissue-specific expression. Results demonstrated a unique gene expression pattern for EaDmrtE, which is enriched in female germ cells during the spawning season. Based on the phylogenetic analyses performed across the homologous Dmrt genes in metazoans, we found that the female-specific EaDmrtE gene is not related to the DM1 gene of Acropora spp. coral nor to Dmrt1 of vertebrates, which are involved in sexual reproduction, especially in sex determination (vertebrate Dmrt1). Additionally, high levels of EaDmrtE transcripts detected in unfertilized mature eggs are retained in newly formed zygotes but decrease during embryonic development. We suggest that the newly discovered gene may play a role in oogenesis and early embryogenesis as a maternal factor in corals. Therefore, the sexual reproduction-associated Dmrt gene(s) should have arisen in cnidarians and might have evolved multiple times in metazoans.

Timing and number of colonizations but not diversification rates affect diversity patterns in hemosporidian lineages on a remote oceanic archipelago.

Parasite diversity on remote oceanic archipelagos is determined by the number and timing of colonizations and by in situ diversification rate. In this study, we compare intra-archipelago diversity of two hemosporidian parasite genera, Plasmodium and Leucocytozoon, infecting birds of the Mascarene archipelago. Despite the generally higher vagility of Plasmodium parasites, we report a diversity of Plasmodium cytochrome b haplotypes in the archipelago much lower than that of Leucocytozoon. Using phylogenetic data, we find that this difference in diversity is consistent with differences in the timing and number of colonizations, while rates of diversification do not vary significantly between the two genera. The prominence of immigration history in explaining current diversity patterns highlights the importance of historical contingencies in driving disparate biogeographic patterns in potentially harmful blood parasites infecting island birds.

Mass production of SNP markers in a nonmodel passerine bird through RAD sequencing and contig mapping to the zebra finch genome.

Here, we present an adaptation of restriction-site-associated DNA sequencing (RAD-seq) to the Illumina HiSeq2000 technology that we used to produce SNP markers in very large quantities at low cost per unit in the Réunion grey white-eye (Zosterops borbonicus), a nonmodel passerine bird species with no reference genome. We sequenced a set of six pools of 18-25 individuals using a single sequencing lane. This allowed us to build around 600 000 contigs, among which at least 386 000 could be mapped to the zebra finch (Taeniopygia guttata) genome. This yielded more than 80 000 SNPs that could be mapped unambiguously and are evenly distributed across the genome. Thus, our approach provides a good illustration of the high potential of paired-end RAD sequencing of pooled DNA samples combined with comparative assembly to the zebra finch genome to build large contigs and characterize vast numbers of informative SNPs in nonmodel passerine bird species in a very efficient and cost-effective way.

Investigating the role of the melanocortin-1 receptor gene in an extreme case of microgeographical variation in the pattern of melanin-based plumage pigmentation.

UNLABELLED: The Réunion grey white-eye (Zosterops borbonicus) is a single-island endemic passerine bird that exhibits striking geographically structured melanic polymorphism at a very small spatial scale. We investigated the genetic basis of this color polymorphism by testing whether the melanocortin-1 receptor (MC1R), a gene often involved in natural melanic polymorphism in birds, was associated with the observed plumage variation. Although we found three non-synonymous mutations, we detected no association between MC1R variants and color morphs, and the main amino-acid variant found in the Réunion grey white-eye was also present at high frequency in the Mauritius grey white-eye (Zosterops mauritianus), its sister species which shows no melanic polymorphism. In addition, neutrality tests and analysis of population structure did not reveal any obvious pattern of positive or balancing selection acting on MC1R. Altogether these results indicate that MC1R does not play a role in explaining the melanic variation observed in the Réunion grey white-eye. We propose that other genes such as POMC, Agouti or any other genes involved in pigment synthesis will need to be investigated in future studies if we are to understand how selection shapes complex patterns of melanin-based plumage pigmentation. TRIAL REGISTRATION: All sequences submitted to Genbank. Accession number: JX914505 to JX914564.