Molecular phylogeny and evolution of inflorescence types in Eperua.

PREMISE: The Amazonian hyperdominant genus Eperua (Fabaceae) currently holds 20 described species and has two strongly different inflorescence and flower types, with corresponding different pollination syndrome. The evolution of these vastly different inflorescence types within this genus was unknown and the main topic in this study. METHODS: We constructed a molecular phylogeny, based on the full nuclear ribosomal DNA and partial plastome, using Bayesian inference and maximum likelihood methods, to test whether the genus is monophyletic, whether all species are monophyletic and if the shift from bat to bee pollination (or vice versa) occurred once in this genus. RESULTS: All but two species are well supported by the nuclear ribosomal phylogeny. The plastome phylogeny, however, shows a strong geographic signal suggesting strong local hybridization or chloroplast capture, rendering chloroplast barcodes meaningless in this genus. CONCLUSIONS: With our data, we cannot fully resolve the backbone of the tree to clarify sister genera relationships and confirm monophyly of the genus Eperua. Within the genus, the shift from bat to bee and bee to bat pollination has occurred several times but, with the bee to bat not always leading to a pendant inflorescence.

Interpreting phylogenetic conflict: Hybridization in the most speciose genus of lichen-forming fungi.

While advances in sequencing technologies have been invaluable for understanding evolutionary relationships, increasingly large genomic data sets may result in conflicting evolutionary signals that are often caused by biological processes, including hybridization. Hybridization has been detected in a variety of organisms, influencing evolutionary processes such as generating reproductive barriers and mixing standing genetic variation. Here, we investigate the potential role of hybridization in the diversification of the most speciose genus of lichen-forming fungi, Xanthoparmelia. As Xanthoparmelia is projected to have gone through recent, rapid diversification, this genus is particularly suitable for investigating and interpreting the origins of phylogenomic conflict. Focusing on a clade of Xanthoparmelia largely restricted to the Holarctic region, we used a genome skimming approach to generate 962 single-copy gene regions representing over 2 Mbp of the mycobiont genome. From this genome-scale dataset, we inferred evolutionary relationships using both concatenation and coalescent-based species tree approaches. We also used three independent tests for hybridization. Although different species tree reconstruction methods recovered largely consistent and well-supported trees, there was widespread incongruence among individual gene trees. Despite challenges in differentiating hybridization from ILS in situations of recent rapid radiations, our genome-wide analyses detected multiple potential hybridization events in the Holarctic clade, suggesting one possible source of trait variability in this hyperdiverse genus. This study highlights the value in using a pluralistic approach for characterizing genome-scale conflict, even in groups with well-resolved phylogenies, while highlighting current challenges in detecting the specific impacts of hybridization.

The genomic basis of army ant chemosensory adaptations.

The evolution of mass raiding has allowed army ants to become dominant arthropod predators in the tropics. Although a century of research has led to many discoveries about behavioural, morphological and physiological adaptations in army ants, almost nothing is known about the molecular basis of army ant biology. Here we report the genome of the iconic New World army ant Eciton burchellii, and show that it is unusually compact, with a reduced gene complement relative to other ants. In contrast to this overall reduction, a particular gene subfamily (9-exon ORs) expressed predominantly in female antennae is expanded. This subfamily has previously been linked to the recognition of hydrocarbons, key olfactory cues used in insect communication and prey discrimination. Confocal microscopy of the brain showed a corresponding expansion in a putative hydrocarbon response centre within the antennal lobe, while scanning electron microscopy of the antenna revealed a particularly high density of hydrocarbon-sensitive sensory hairs. E. burchellii shares these features with its predatory and more cryptic relative, the clonal raider ant. By integrating genomic, transcriptomic and anatomical analyses in a comparative context, our work thus provides evidence that army ants and their relatives possess a suite of modifications in the chemosensory system that may be involved in behavioural coordination and prey selection during social predation. It also lays the groundwork for future studies of army ant biology at the molecular level.

Population genomics of a reindeer lichen species from North American lichen woodlands.

PREMISE: Lichens are one of the main structural components of plant communities in the North American boreal biome. They play a pivotal role in lichen woodlands, a large ecosystem situated north of the closed-crown forest zone, and south of the forest-tundra zone. In Eastern Canada (Quebec), there is a remnant LW found 500 km south of its usual distribution range, in the Parc National des Grands-Jardins, originated mainly because of wildfires. We inferred the origin of the lichen Cladonia stellaris from this LW and assessed its genetic diversity in a postfire succession. METHODS: We genotyped 122 individuals collected across a latitudinal gradient in Quebec. Using the software Stacks, we compared four different approaches of locus selection and single-nucleotide polymorphism calling. We identified the best fitting approach to investigate population structure and estimate genetic diversity of C. stellaris. RESULTS: Populations in southern Quebec are not genetically different from those of northern LWs. The species consists of at least four phylogenetic lineages with elevated levels of genetic diversity and low co-ancestry. In Parc National des Grands-Jardins, we reported high values of genetic diversity not related with time since fire disturbance and low genetic differentiation among populations with different fire histories. CONCLUSIONS: This first population genomic study of C. stellaris is an important step forward to understand the origin and biogeographic patterns of lichen woodlands in North America. Our findings also contribute to the understanding of the effect of postfire succession on the genetic structure of the species.

Museum genomics reveals the Xerces blue butterfly (Glaucopsyche xerces) was a distinct species driven to extinction.

The last Xerces blue butterfly was seen in the early 1940s, and its extinction is credited to human urban development. This butterfly has become a North American icon for insect conservation, but some have questioned whether it was truly a distinct species, or simply an isolated population of another living species. To address this question, we leveraged next-generation sequencing using a 93-year-old museum specimen. We applied a genome skimming strategy that aimed for the organellar genome and high-copy fractions of the nuclear genome by a shallow sequencing approach. From these data, we were able to recover over 200 million nucleotides, which assembled into several phylogenetically informative markers and the near-complete mitochondrial genome. From our phylogenetic analyses and haplotype network analysis we conclude that the Xerces blue butterfly was a distinct species driven to extinction.

Characterizing the ribosomal tandem repeat and its utility as a DNA barcode in lichen-forming fungi.

BACKGROUND: Regions within the nuclear ribosomal operon are a major tool for inferring evolutionary relationships and investigating diversity in fungi. In spite of the prevalent use of ribosomal markers in fungal research, central features of nuclear ribosomal DNA (nrDNA) evolution are poorly characterized for fungi in general, including lichenized fungi. The internal transcribed spacer (ITS) region of the nrDNA has been adopted as the primary DNA barcode identification marker for fungi. However, little is known about intragenomic variation in the nrDNA in symbiotic fungi. In order to better understand evolution of nrDNA and the utility of the ITS region for barcode identification of lichen-forming fungal species, we generated nearly complete nuclear ribosomal operon sequences from nine species in the Rhizoplaca melanophthalma species complex using short reads from high-throughput sequencing. RESULTS: We estimated copy numbers for the nrDNA operon, ranging from nine to 48 copies for members of this complex, and found low levels of intragenomic variation in the standard barcode region (ITS). Monophyly of currently described species in this complex was supported in phylogenetic inferences based on the ITS, 28S, intergenic spacer region, and some intronic regions, independently; however, a phylogenetic inference based on the 18S provided much lower resolution. Phylogenetic analysis of concatenated ITS and intergenic spacer sequence data generated from 496 specimens collected worldwide revealed previously unrecognized lineages in the nrDNA phylogeny. CONCLUSIONS: The results from our study support the general assumption that the ITS region of the nrDNA is an effective barcoding marker for fungi. For the R. melanophthalma group, the limited amount of potential intragenomic variability in the ITS region did not correspond to fixed diagnostic nucleotide position characters separating taxa within this species complex. Previously unrecognized lineages inferred from ITS sequence data may represent undescribed species-level lineages or reflect uncharacterized aspects of nrDNA evolution in the R. melanophthalma species complex.

Lycophyte plastid genomics: extreme variation in GC, gene and intron content and multiple inversions between a direct and inverted orientation of the rRNA repeat.

Lycophytes are a key group for understanding vascular plant evolution. Lycophyte plastomes are highly distinct, indicating a dynamic evolutionary history, but detailed evaluation is hindered by the limited availability of sequences. Eight diverse plastomes were sequenced to assess variation in structure and functional content across lycophytes. Lycopodiaceae plastomes have remained largely unchanged compared with the common ancestor of land plants, whereas plastome evolution in Isoetes and especially Selaginella is highly dynamic. Selaginella plastomes have the highest GC content and fewest genes and introns of any photosynthetic land plant. Uniquely, the canonical inverted repeat was converted into a direct repeat (DR) via large-scale inversion in some Selaginella species. Ancestral reconstruction identified additional putative transitions between an inverted and DR orientation in Selaginella and Isoetes plastomes. A DR orientation does not disrupt the activity of copy-dependent repair to suppress substitution rates within repeats. Lycophyte plastomes include the most archaic examples among vascular plants and the most reconfigured among land plants. These evolutionary trends correlate with the mitochondrial genome, suggesting shared underlying mechanisms. Copy-dependent repair for DR-localized genes indicates that recombination and gene conversion are not inhibited by the DR orientation. Gene relocation in lycophyte plastomes occurs via overlapping inversions rather than transposase/recombinase-mediated processes.

The Reverse Transcriptase/RNA Maturase Protein MatR Is Required for the Splicing of Various Group II Introns in Brassicaceae Mitochondria.

Group II introns are large catalytic RNAs that are ancestrally related to nuclear spliceosomal introns. Sequences corresponding to group II RNAs are found in many prokaryotes and are particularly prevalent within plants organellar genomes. Proteins encoded within the introns themselves (maturases) facilitate the splicing of their own host pre-RNAs. Mitochondrial introns in plants have diverged considerably in sequence and have lost their maturases. In angiosperms, only a single maturase has been retained in the mitochondrial DNA: the matR gene found within NADH dehydrogenase 1 (nad1) intron 4. Its conservation across land plants and RNA editing events, which restore conserved amino acids, indicates that matR encodes a functional protein. However, the biological role of MatR remains unclear. Here, we performed an in vivo investigation of the roles of MatR in Brassicaceae. Directed knockdown of matR expression via synthetically designed ribozymes altered the processing of various introns, including nad1 i4. Pull-down experiments further indicated that MatR is associated with nad1 i4 and several other intron-containing pre-mRNAs. MatR may thus represent an intermediate link in the gradual evolutionary transition from the intron-specific maturases in bacteria into their versatile spliceosomal descendants in the nucleus. The similarity between maturases and the core spliceosomal Prp8 protein further supports this intriguing theory.

Monilophyte mitochondrial rps1 genes carry a unique group II intron that likely originated from an ancient paralog in rpl2.

Intron patterns in plant mitochondrial genomes differ significantly between the major land plant clades. We here report on a new, clade-specific group II intron in the rps1 gene of monilophytes (ferns). This intron, rps1i25g2, is strikingly similar to rpl2i846g2 previously identified in the mitochondrial rpl2 gene of seed plants, ferns, and the lycophyte Phlegmariurus squarrosus Although mitochondrial ribosomal protein genes are frequently subject to endosymbiotic gene transfer among plants, we could retrieve the mitochondrial rps1 gene in a taxonomically wide sampling of 44 monilophyte taxa including basal lineages such as the Ophioglossales, Psilotales, and Marattiales with the only exception being the Equisetales (horsetails). Introns rps1i25g2 and rpl2i846g2 were likewise consistently present with only two exceptions: Intron rps1i25g2 is lost in the genus Ophioglossum and intron rpl2i846g2 is lost in Equisetum bogotense Both intron sequences are moderately affected by RNA editing. The unprecedented primary and secondary structure similarity of rps1i25g2 and rpl2i846g2 suggests an ancient retrotransposition event copying rpl2i846g2 into rps1, for which we suggest a model. Our phylogenetic analysis adding the new rps1 locus to a previous data set is fully congruent with recent insights on monilophyte phylogeny and further supports a sister relationship of Gleicheniales and Hymenophyllales.

Ginkgo and Welwitschia Mitogenomes Reveal Extreme Contrasts in Gymnosperm Mitochondrial Evolution.

Mitochondrial genomes (mitogenomes) of flowering plants are well known for their extreme diversity in size, structure, gene content, and rates of sequence evolution and recombination. In contrast, little is known about mitogenomic diversity and evolution within gymnosperms. Only a single complete genome sequence is available, from the cycad Cycas taitungensis, while limited information is available for the one draft sequence, from Norway spruce (Picea abies). To examine mitogenomic evolution in gymnosperms, we generated complete genome sequences for the ginkgo tree (Ginkgo biloba) and a gnetophyte (Welwitschia mirabilis). There is great disparity in size, sequence conservation, levels of shared DNA, and functional content among gymnosperm mitogenomes. The Cycas and Ginkgo mitogenomes are relatively small, have low substitution rates, and possess numerous genes, introns, and edit sites; we infer that these properties were present in the ancestral seed plant. By contrast, the Welwitschia mitogenome has an expanded size coupled with accelerated substitution rates and extensive loss of these functional features. The Picea genome has expanded further, to more than 4 Mb. With regard to structural evolution, the Cycas and Ginkgo mitogenomes share a remarkable amount of intergenic DNA, which may be related to the limited recombinational activity detected at repeats in Ginkgo Conversely, the Welwitschia mitogenome shares almost no intergenic DNA with any other seed plant. By conducting the first measurements of rates of DNA turnover in seed plant mitogenomes, we discovered that turnover rates vary by orders of magnitude among species.

Reverse U-to-C editing exceeds C-to-U RNA editing in some ferns - a monilophyte-wide comparison of chloroplast and mitochondrial RNA editing suggests independent evolution of the two processes in both organelles.

BACKGROUND: RNA editing by C-to-U conversions is nearly omnipresent in land plant chloroplasts and mitochondria, where it mainly serves to reconstitute conserved codon identities in the organelle mRNAs. Reverse U-to-C RNA editing in contrast appears to be restricted to hornworts, some lycophytes, and ferns (monilophytes). A well-resolved monilophyte phylogeny has recently emerged and now allows to trace the side-by-side evolution of both types of pyrimidine exchange editing in the two endosymbiotic organelles. RESULTS: Our study of RNA editing in four selected mitochondrial genes show a wide spectrum of divergent RNA editing frequencies including a dominance of U-to-C over the canonical C-to-U editing in some taxa like the order Schizaeales. We find that silent RNA editing leaving encoded amino acids unchanged is highly biased with more than ten-fold amounts of silent C-to-U over U-to-C edits. In full contrast to flowering plants, RNA editing frequencies are low in early-branching monilophyte lineages but increase in later emerging clades. Moreover, while editing rates in the two organelles are usually correlated, we observe uncoupled evolution of editing frequencies in fern mitochondria and chloroplasts. Most mitochondrial RNA editing sites are shared between the recently emerging fern orders whereas chloroplast editing sites are mostly clade-specific. Finally, we observe that chloroplast RNA editing appears to be completely absent in horsetails (Equisetales), the sister clade of all other monilophytes. CONCLUSIONS: C-to-U and U-to-C RNA editing in fern chloroplasts and mitochondria follow disinct evolutionary pathways that are surprisingly different from what has previously been found in flowering plants. The results call for careful differentiation of the two types of RNA editing in the two endosymbiotic organelles in comparative evolutionary studies.

Laser capture microdissection microscopy and genome sequencing of the avian malaria parasite, Plasmodium relictum.

Acquiring genomic material from avian malaria parasites for genome sequencing has proven problematic due to the nucleation of avian erythrocytes, which produces a large ratio of host to parasite DNA (∼1 million to 1 bp). We tested the ability of laser capture microdissection microscopy to isolate parasite cells from individual avian erythrocytes for four avian Plasmodium species, and subsequently applied whole genome amplification and Illumina sequencing methods to Plasmodium relictum (lineage pSGS1) to produce sequence reads of the P. relictum genome. We assembled ∼335 kbp of parasite DNA from this species, but were unable to completely avoid contamination by host DNA and other sources. However, it is clear that laser capture microdissection holds promise for the isolation of genomic material from haemosporidian parasites in intracellular life stages. In particular, laser capture microdissection may prove useful for isolating individual parasite species from co-infected hosts. Although not explicitly tested in this study, laser capture microdissection may also have important applications for isolation of rare parasite lineages and museum specimens for which no fresh material exists.

Dynamic evolution of Geranium mitochondrial genomes through multiple horizontal and intracellular gene transfers.

The exchange of genetic material between cellular organelles through intracellular gene transfer (IGT) or between species by horizontal gene transfer (HGT) has played an important role in plant mitochondrial genome evolution. The mitochondrial genomes of Geraniaceae display a number of unusual phenomena including highly accelerated rates of synonymous substitutions, extensive gene loss and reduction in RNA editing. Mitochondrial DNA sequences assembled for 17 species of Geranium revealed substantial reduction in gene and intron content relative to the ancestor of the Geranium lineage. Comparative analyses of nuclear transcriptome data suggest that a number of these sequences have been functionally relocated to the nucleus via IGT. Evidence for rampant HGT was detected in several Geranium species containing foreign organellar DNA from diverse eudicots, including many transfers from parasitic plants. One lineage has experienced multiple, independent HGT episodes, many of which occurred within the past 5.5 Myr. Both duplicative and recapture HGT were documented in Geranium lineages. The mitochondrial genome of Geranium brycei contains at least four independent HGT tracts that are absent in its nearest relative. Furthermore, G. brycei mitochondria carry two copies of the cox1 gene that differ in intron content, providing insight into contrasting hypotheses on cox1 intron evolution.

Horsetails are the sister group to all other monilophytes and Marattiales are sister to leptosporangiate ferns.

The "Monilophyte" clade comprising ferns, horsetails and whisk ferns receives unequivocal support from molecular data as the sister clade to seed plants. However, the branching order of its earliest emerging lineages, the Equisetales (horsetails), the Marattiales, the Ophioglossales/Psilotales and the large group of leptosporangiate ferns has remained dubious. We investigated the mitochondrial nad2 and rpl2 genes as two new, intron-containing loci for a wide sampling of taxa. We found that both group II introns - nad2i542g2 and rpl2i846g2 - are universally present among monilophytes. Both introns have orthologues in seed plants where nad2i542g2 has evolved into a trans-arrangement. In contrast and despite substantial size extensions to more than 5kb in Psilotum, nad2i542g2 remains cis-arranged in the monilophytes. For phylogenetic analyses, we filled taxonomic gaps in previously investigated mitochondrial (atp1, nad5) and chloroplast (atpA, atpB, matK, rbcL, rps4) loci and created a 9-gene matrix that also included the new mitochondrial nad2 and rpl2 loci. We extended the taxon sampling with two taxa each for all land plant outgroups (liverworts, mosses, hornworts, lycophytes and seed plants) to minimize the risk of phylogenetic artefacts. We ultimately obtained a well-supported molecular phylogeny placing Marattiales as sister to leptosporangiate ferns and horsetails as sister to all remaining monilophytes. In addition, an indel in an exon of the here introduced rpl2 locus independently supports the placement of horsetails. We conclude that under dense taxon sampling, phylogenetic information from a prudent choice of loci is currently superior to character-rich phylogenomic approaches at low taxon sampling. As here shown the selective choice of loci and taxa enabled us to resolve the long-enigmatic diversifications of the earliest monilophyte lineages.

Complete plastid genomes from Ophioglossum californicum, Psilotum nudum, and Equisetum hyemale reveal an ancestral land plant genome structure and resolve the position of Equisetales among monilophytes.

BACKGROUND: Plastid genome structure and content is remarkably conserved in land plants. This widespread conservation has facilitated taxon-rich phylogenetic analyses that have resolved organismal relationships among many land plant groups. However, the relationships among major fern lineages, especially the placement of Equisetales, remain enigmatic. RESULTS: In order to understand the evolution of plastid genomes and to establish phylogenetic relationships among ferns, we sequenced the plastid genomes from three early diverging species: Equisetum hyemale (Equisetales), Ophioglossum californicum (Ophioglossales), and Psilotum nudum (Psilotales). A comparison of fern plastid genomes showed that some lineages have retained inverted repeat (IR) boundaries originating from the common ancestor of land plants, while other lineages have experienced multiple IR changes including expansions and inversions. Genome content has remained stable throughout ferns, except for a few lineage-specific losses of genes and introns. Notably, the losses of the rps16 gene and the rps12i346 intron are shared among Psilotales, Ophioglossales, and Equisetales, while the gain of a mitochondrial atp1 intron is shared between Marattiales and Polypodiopsida. These genomic structural changes support the placement of Equisetales as sister to Ophioglossales + Psilotales and Marattiales as sister to Polypodiopsida. This result is augmented by some molecular phylogenetic analyses that recover the same relationships, whereas others suggest a relationship between Equisetales and Polypodiopsida. CONCLUSIONS: Although molecular analyses were inconsistent with respect to the position of Marattiales and Equisetales, several genomic structural changes have for the first time provided a clear placement of these lineages within the ferns. These results further demonstrate the power of using rare genomic structural changes in cases where molecular data fail to provide strong phylogenetic resolution.

A unique transcriptome: 1782 positions of RNA editing alter 1406 codon identities in mitochondrial mRNAs of the lycophyte Isoetes engelmannii.

The analysis of the mitochondrial DNA of Isoetes engelmannii as a first representative of the lycophytes recently revealed very small introns and indications for extremely frequent RNA editing. To analyze functionality of intron splicing and the extent of RNA editing in I. engelmannii, we performed a comprehensive analysis of its mitochondrial transcriptome. All 30 groups I and II introns were found to be correctly removed, showing that intron size reduction does not impede splicing. We find that mRNA editing affects 1782 sites, which lead to a total of 1406 changes in codon meanings. This includes the removal of stop codons from 23 of the 25 mitochondrial protein encoding genes. Comprehensive sequence analysis of multiple cDNAs per locus allowed classification of partially edited sites as either inefficiently edited but relevant or as non-specifically edited at mostly low frequencies. Abundant RNA editing was also found to affect tRNAs in hitherto unseen frequency, taking place at 41 positions in tRNA-precursors, including the first identification of U-to-C exchanges in two tRNA species. We finally investigated the four group II introns of the nad7 gene and could identify 27 sites of editing, most of which improve base pairing for proper secondary structure formation.

Reference-Based RADseq Unravels the Evolutionary History of Polar Species in 'the Crux Lichenologorum' Genus Usnea (Parmeliaceae, Ascomycota).

Nearly 90% of fungal diversity, one of the most speciose branches in the tree of life, remains undescribed. Lichenized fungi as symbiotic associations are still a challenge for species delimitation, and current species diversity is vastly underestimated. The ongoing democratization of Next-Generation Sequencing is turning the tables. Particularly, reference-based RADseq allows for metagenomic filtering of the symbiont sequence and yields robust phylogenomic trees of closely related species. We implemented reference-based RADseq to disentangle the evolution of neuropogonoid lichens, which inhabit harsh environments and belong to Usnea (Parmeliaceae, Ascomycota), one of the most taxonomically intriguing genera within lichenized fungi. Full taxon coverage of neuropogonoid lichens was sampled for the first time, coupled with phenotype characterizations. More than 20,000 loci of 126 specimens were analyzed through concatenated and coalescent-based methods, including time calibrations. Our analysis addressed the major taxonomic discussions over recent decades. Subsequently, two species are newly described, namely U. aymondiana and U. fibriloides, and three species names are resurrected. The late Miocene and Pliocene-Pleistocene boundary is inferred as the timeframe for neuropogonoid lichen diversification. Ultimately, this study helped fill the gap of fungal diversity by setting a solid backbone phylogeny which raises new questions about which factors may trigger complex evolutionary scenarios.

Reference-Based Restriction-Site-Associated DNA Sequencing Data Are Useful for Species Delineation in a Recently Diverged Asexually Reproducing Species Complex (Parmeliaceae, Ascomycota).

Cryptic species are common in lichen-forming fungi and have been reported from different genera in the most speciose family, Parmeliaceae. Herein, we address species delimitation in a group of mainly asexually reproducing Parmelina species. The morphologically distinct P. pastillifera was previously found nested within a morphologically circumscribed P. tiliacea based on several loci. However, these studies demonstrated a relatively high genetic diversity within P. tiliacea sensu lato. Here, we revisit the species delimitation in the group by analyzing single-nucleotide polymorphisms (SNPs) through genome-wide assessment using Restriction-Site-Associated sequencing and population genomic methods. Our data support previous studies and provide further insight into the phylogenetic relationships of the four clades found within the complex. Based on the evidence suggesting a lack of gene flow among the clades, we recognize the four clades as distinct species, P. pastillifera and P. tiliacea sensu stricto, and two new species, P. clandestina sp. nov. and P. mediterranea sp. nov.

Early stages of speciation with gene flow in the Amazilia Hummingbird (Amazilis amazilia) subspecies complex of Western South America.

Disentangling the factors underlying the diversification of geographically variable species with a wide geographical range is essential to understanding the initial stages and drivers of the speciation process. The Amazilia Hummingbird, Amazilis amazilia, is found along the Pacific coast from northern Ecuador down to the Nazca Valley of Peru, and is currently classified as six phenotypically differentiated subspecies. We aimed to resolve the evolutionary relationships of the six subspecies, to assess the geographical pattern and extent of evolutionary divergence, and to test for introgression using both a mtDNA marker and a genome-by-sequencing dataset from 86 individuals from across the species range. The consensus phylogenetic tree separated the six subspecies into three distinct clades, corresponding with the Ecuador lowlands (A. amazilia dumerilii), the Ecuador highlands (A. amazilia alticola and A. amazilia azuay), and the Peruvian coast (A. amazilia leucophoea, A. amazilia amazilia, and A. amazilia caeruleigularis). However, an unresolved mtDNA network suggests that the diversification of the subspecies was recent and rapid. We found evidence of gene flow among the subspecies A. amazilia dumerilii, A. amazilia alticola, and A. amazilia leucophoea, with strong genetic isolation of the subspecies A. amazilia azuay in the isolated Yunguilla Valley of Ecuador. Finally, environmental data from each subspecies' capture locations were concordant with the three distinct clades. Overall, our results suggest that both expansions into new habitats and geographic isolation shaped the present-day phylogeny and range of the A. amazilia subspecies, and that A. amazilia azuay may be genetically divergent enough to be considered a separate species.

Should we hail the Red King? Evolutionary consequences of a mutualistic lifestyle in genomes of lichenized ascomycetes.

The Red Queen dynamic is often brought into play for antagonistic relationships. However, the coevolutionary effects of mutualistic interactions, which predict slower evolution for interacting organisms (Red King), have been investigated to a lesser extent. Lichens are a stable, mutualistic relationship of fungi and cyanobacteria and/or algae, which originated several times independently during the evolution of fungi. Therefore, they represent a suitable system to investigate the coevolutionary effect of mutualism on the fungal genome. We measured substitution rates and selective pressure of about 2000 protein-coding genes (plus the rDNA region) in two different classes of Ascomycota, each consisting of closely related lineages of lichenized and non-lichenized fungi. Our results show that independent lichenized clades are characterized by significantly slower rates for both synonymous and non-synonymous substitutions. We hypothesize that this evolutionary pattern is connected to the lichen life cycle (longer generation time of lichenized fungi) rather than a result of different selection strengths, which is described as the main driver for the Red Kind dynamic. This first empirical evidence of slower evolution in lichens provides an important insight on how biotic cooperative interactions are able to shape the evolution of symbiotic organisms.