Genomic basis for skin phenotype and cold adaptation in the extinct Steller's sea cow.

Steller's sea cow, an extinct sirenian and one of the largest Quaternary mammals, was described by Georg Steller in 1741 and eradicated by humans within 27 years. Here, we complement Steller's descriptions with paleogenomic data from 12 individuals. We identified convergent evolution between Steller's sea cow and cetaceans but not extant sirenians, suggesting a role of several genes in adaptation to cold aquatic (or marine) environments. Among these are inactivations of lipoxygenase genes, which in humans and mouse models cause ichthyosis, a skin disease characterized by a thick, hyperkeratotic epidermis that recapitulates Steller's sea cows' reportedly bark-like skin. We also found that Steller's sea cows' abundance was continuously declining for tens of thousands of years before their description, implying that environmental changes also contributed to their extinction.

Three mitochondrial genomes of early-winged insects (Ephemeroptera: Baetidae and Leptophlebiidae).

Mayflies (Ephemeroptera) are a semi-aquatic insect order with comparatively few genomic data available despite their phylogenetic position at the root of the winged-insects and possession of ancestral traits. Here, we provide three mitochondrial genomes (mtgenomes) from representatives of the two most species-rich families, Baetis rutilocylindratus and Cloeon dipterum (Baetidae), and Habrophlebiodes zijinensis (Leptophlebiidae). All mtgenomes had a complete set of 13 protein-coding genes and a conserved orientation except for two inverted tRNAs in H. zijinensis. Phylogenetic reconstructions using 21 mayfly mtgenomes and representatives of seven additional orders recovered both Baetidae and Leptophlebiidae as well supported monophyletic clades, with Ephemeroptera as the sister-taxon to all other winged insects (i.e. Odonata and Neoptera).

Ecological Specialization and Evolutionary Reticulation in Extant Hyaenidae.

During the Miocene, Hyaenidae was a highly diverse family of Carnivora that has since been severely reduced to four species: the bone-cracking spotted, striped, and brown hyenas, and the specialized insectivorous aardwolf. Previous studies investigated the evolutionary histories of the spotted and brown hyenas, but little is known about the remaining two species. Moreover, the genomic underpinnings of scavenging and insectivory, defining traits of the extant species, remain elusive. Here, we generated an aardwolf genome and analyzed it together with the remaining three species to reveal their evolutionary relationships, genomic underpinnings of their scavenging and insectivorous lifestyles, and their respective genetic diversities and demographic histories. High levels of phylogenetic discordance suggest gene flow between the aardwolf lineage and the ancestral brown/striped hyena lineage. Genes related to immunity and digestion in the bone-cracking hyenas and craniofacial development in the aardwolf showed the strongest signals of selection, suggesting putative key adaptations to carrion and termite feeding, respectively. A family-wide expansion in olfactory receptor genes suggests that an acute sense of smell was a key early adaptation. Finally, we report very low levels of genetic diversity within the brown and striped hyenas despite no signs of inbreeding, putatively linked to their similarly slow decline in effective population size over the last ∼2 million years. High levels of genetic diversity and more stable population sizes through time are seen in the spotted hyena and aardwolf. Taken together, our findings highlight how ecological specialization can impact the evolutionary history, demographics, and adaptive genetic changes of an evolutionary lineage.

Elevation in tropical sky islands as the common driver in structuring genes and communities of freshwater organisms.

Tropical mountains are usually characterized by a vertically-arranged sequence of ecological belts, which, in contrast to temperate habitats, have remained relatively stable in space across the Quaternary. Such long-lasting patterning of habitats makes them ideal to test the role of environmental pressure in driving ecological and evolutionary processes. Using Sumatran freshwater mayfly communities, we test whether elevation, rather than other spatial factors (i.e. volcanoes, watersheds) structures both species within communities and genes within species. Based on the analysis of 31 mayfly (Ephemeroptera) communities and restriction-site-associated-DNA sequencing in the four most ubiquitous species, we found elevation as the major spatial component structuring both species and genes in the landscape. In other words, similar elevations across different mountains or watersheds harbor more similar species and genes than different elevations within the same mountain or watershed. Tropical elevation gradients characterized by environmental conditions that are both steep and relatively stable seasonally and over geological time scales, are thus responsible for both ecological and genetic differentiation. Our results demonstrate how in situ ecological diversification at the micro-evolutionary level might fuel alpha- and beta- components of diversity in tropical sky islands.

Assessing the potential of RAD-sequencing to resolve phylogenetic relationships within species radiations: The fly genus Chiastocheta (Diptera: Anthomyiidae) as a case study.

Determining phylogenetic relationships among recently diverged species has long been a challenge in evolutionary biology. Cytoplasmic DNA markers, which have been widely used, notably in the context of molecular barcoding, have not always proved successful in resolving such phylogenies. However, with the advent of next-generation-sequencing technologies and associated techniques of reduced genome representation, phylogenies of closely related species have been resolved at a much higher detail in the last couple of years. Here we examine the potential and limitations of one of such techniques-Restriction-site Associated DNA (RAD) sequencing, a method that produces thousands of (mostly) anonymous nuclear markers, in disentangling the phylogeny of the fly genus Chiastocheta (Diptera: Anthomyiidae). In Europe, this genus encompasses seven species of seed predators, which have been widely studied in the context of their ecological and evolutionary interactions with the plant Trollius europaeus (Ranunculaceae). So far, phylogenetic analyses using mitochondrial markers failed to resolve monophyly of most of the species from this recently diversified genus, suggesting that their taxonomy may need a revision. However, relying on a single, non-recombining marker and ignoring potential incongruences between mitochondrial and nuclear loci may provide an incomplete account of the lineage history. In this study, we applied both classical Sanger sequencing of three mtDNA regions and RAD-sequencing, for reconstructing the phylogeny of the genus. Contrasting with results based on mitochondrial markers, RAD-sequencing analyses retrieved the monophyly of all seven species, in agreement with the morphological species assignment. We found robust nuclear-based species assignment of individual samples, and low levels of estimated contemporary gene flow among them. However, despite recovering species' monophyly, interspecific relationships varied depending on the set of RAD loci considered, producing contradictory topologies. Moreover, coalescence-based phylogenetic analyses revealed low supports for most of the interspecific relationships. Our results indicate that despite the higher performance of RAD-sequencing in terms of species trees resolution compared to cytoplasmic markers, reconstructing inter-specific relationships among recently-diverged lineages may lie beyond the possibilities offered by large sets of RAD-sequencing markers in cases of strong gene tree incongruence.

Colonization and diversification of aquatic insects on three Macaronesian archipelagos using 59 nuclear loci derived from a draft genome.

The study of processes driving diversification requires a fully sampled and well resolved phylogeny, although a lack of phylogenetic markers remains a limitation for many non-model groups. Multilocus approaches to the study of recent diversification provide a powerful means to study the evolutionary process, but their application remains restricted because multiple unlinked loci with suitable variation for phylogenetic or coalescent analysis are not available for most non-model taxa. Here we identify novel, putative single-copy nuclear DNA (nDNA) phylogenetic markers to study the colonization and diversification of an aquatic insect species complex, Cloeon dipterum L. 1761 (Ephemeroptera: Baetidae), in Macaronesia. Whole-genome sequencing data from one member of the species complex were used to identify 59 nDNA loci (32,213 base pairs), followed by Sanger sequencing of 29 individuals sampled from 13 islands of three Macaronesian archipelagos. Multispecies coalescent analyses established six putative species. Three island species formed a monophyletic clade, with one species occurring on the Azores, Europe and North America. Ancestral state reconstruction indicated at least two colonization events from the mainland (to the Canaries, respectively Azores) and one within the archipelago (between Madeira and the Canaries). Random subsets of the 59 loci showed a positive linear relationship between number of loci and node support. In contrast, node support in the multispecies coalescent tree was negatively correlated with mean number of phylogenetically informative sites per locus, suggesting a complex relationship between tree resolution and marker variability. Our approach highlights the value of combining genomics, coalescent-based phylogeography, species delimitation, and phylogenetic reconstruction to resolve recent diversification events in an archipelago species complex.

discomark: nuclear marker discovery from orthologous sequences using draft genome data.

High-throughput sequencing has laid the foundation for fast and cost-effective development of phylogenetic markers. Here we present the program discomark, which streamlines the development of nuclear DNA (nDNA) markers from whole-genome (or whole-transcriptome) sequencing data, combining local alignment, alignment trimming, reference mapping and primer design based on multiple sequence alignments to design primer pairs from input orthologous sequences. To demonstrate the suitability of discomark, we designed markers for two groups of species, one consisting of closely related species and one group of distantly related species. For the closely related members of the species complex of Cloeon dipterum s.l. (Insecta, Ephemeroptera), the program discovered a total of 78 markers. Among these, we selected eight markers for amplification and Sanger sequencing. The exon sequence alignments (2526 base pairs) were used to reconstruct a well-supported phylogeny and to infer clearly structured haplotype networks. For the distantly related species, we designed primers for the insect order Ephemeroptera, using available genomic data from four sequenced species. We developed primer pairs for 23 markers that are designed to amplify across several families. The discomark program will enhance the development of new nDNA markers by providing a streamlined, automated approach to perform genome-scale scans for phylogenetic markers. The program is written in Python, released under a public licence (GNU GPL version 2), and together with a manual and example data set available at: https://github.com/hdetering/discomark.

Molecular phylogeny and timing of diversification in Alpine Rhithrogena (Ephemeroptera: Heptageniidae).

BACKGROUND: Larvae of the Holarctic mayfly genus Rhithrogena Eaton, 1881 (Ephemeroptera, Heptageniidae) are a diverse and abundant member of stream and river communities and are routinely used as bio-indicators of water quality. Rhithrogena is well diversified in the European Alps, with a number of locally endemic species, and several cryptic species have been recently detected. While several informal species groups are morphologically well defined, a lack of reliable characters for species identification considerably hampers their study. Their relationships, origin, timing of speciation and mechanisms promoting their diversification in the Alps are unknown. RESULTS: Here we present a species-level phylogeny of Rhithrogena in Europe using two mitochondrial and three nuclear gene regions. To improve sampling in a genus with many cryptic species, individuals were selected for analysis according to a recent DNA-based taxonomy rather than traditional nomenclature. A coalescent-based species tree and a reconstruction based on a supermatrix approach supported five of the species groups as monophyletic. A molecular clock, mapped on the most resolved phylogeny and calibrated using published mitochondrial evolution rates for insects, suggested an origin of Alpine Rhithrogena in the Oligocene/Miocene boundary. A diversification analysis that included simulation of missing species indicated a constant speciation rate over time, rather than any pronounced periods of rapid speciation. Ancestral state reconstructions provided evidence for downstream diversification in at least two species groups. CONCLUSIONS: Our species-level analyses of five gene regions provide clearer definitions of species groups within European Rhithrogena. A constant speciation rate over time suggests that the paleoclimatic fluctuations, including the Pleistocene glaciations, did not significantly influence the tempo of diversification of Alpine species. A downstream diversification trend in the hybrida and alpestris species groups supports a previously proposed headwater origin hypothesis for aquatic insects.

Parallel ecological diversification in Antarctic notothenioid fishes as evidence for adaptive radiation.

Antarctic notothenioid fishes represent a rare example of a marine species flock. They evolved special adaptations to the extreme environment of the Southern Ocean including antifreeze glycoproteins. Although lacking a swim bladder, notothenioids have diversified from their benthic ancestor into a wide array of water column niches, such as epibenthic, semipelagic, cryopelagic and pelagic habitats. Applying stable carbon (C) and nitrogen (N) isotope analyses to gain information on feeding ecology and foraging habitats, we tested whether ecological diversification along the benthic-pelagic axis followed a single directional trend in notothenioids, or whether it evolved independently in several lineages. Population samples of 25 different notothenioid species were collected around the Antarctic Peninsula, the South Orkneys and the South Sandwich Islands. The C and N stable isotope signatures span a broad range (mean δ(13) C and δ(15) N values between -25.4‰ and -21.9‰ and between 8.5‰ and 13.8‰, respectively), and pairwise niche overlap between four notothenioid families was highly significant. Analysis of isotopic disparity-through-time on the basis of Bayesian inference and maximum-likelihood phylogenies, performed on a concatenated mitochondrial (cyt b) and nuclear gene (myh6, Ptr and tbr1) data set (3148 bp), showed that ecological diversification into overlapping feeding niches has occurred multiple times in parallel in different notothenioid families. This convergent diversification in habitat and trophic ecology is a sign of interspecific competition and characteristic for adaptive radiations.