Million-year-old DNA sheds light on the genomic history of mammoths.

Temporal genomic data hold great potential for studying evolutionary processes such as speciation. However, sampling across speciation events would, in many cases, require genomic time series that stretch well back into the Early Pleistocene subepoch. Although theoretical models suggest that DNA should survive on this timescale1, the oldest genomic data recovered so far are from a horse specimen dated to 780-560 thousand years ago2. Here we report the recovery of genome-wide data from three mammoth specimens dating to the Early and Middle Pleistocene subepochs, two of which are more than one million years old. We find that two distinct mammoth lineages were present in eastern Siberia during the Early Pleistocene. One of these lineages gave rise to the woolly mammoth and the other represents a previously unrecognized lineage that was ancestral to the first mammoths to colonize North America. Our analyses reveal that the Columbian mammoth of North America traces its ancestry to a Middle Pleistocene hybridization between these two lineages, with roughly equal admixture proportions. Finally, we show that the majority of protein-coding changes associated with cold adaptation in woolly mammoths were already present one million years ago. These findings highlight the potential of deep-time palaeogenomics to expand our understanding of speciation and long-term adaptive evolution.

Complex patterns of gene flow and convergence in the evolutionary history of the spiral-horned antelopes (Tragelaphini).

The Tragelaphini, also known as spiral-horned antelope, is a phenotypically diverse mammalian tribe comprising a single genus, Tragelaphus. The evolutionary history of this tribe has attracted the attention of taxonomists and molecular geneticists for decades because its diversity is characterised by conflicts between morphological and molecular data as well as between mitochondrial, nuclear and chromosomal DNA. These inconsistencies point to a complex history of ecological diversification, coupled by either phenotypic convergence or introgression. Therefore, to unravel the phylogenetic relationships among spiral-horned antelopes, and to further investigate the role of divergence and gene flow in trait evolution, we sequenced genomes for all nine accepted species of the genus Tragelaphus, including a genome each for the highly divergent bushbuck lineages (T. s. scriptus and T. s. sylvaticus). We successfully reconstructed the Tragelaphus species tree, providing genome-level support for the early Pliocene divergence and monophyly of the nyala (T. angasii) and lesser kudu (T. imberbis), the monophyly of the two eland species (T. oryx and T. derbianus) and, importantly, the monophyly of kéwel (T. s. scriptus) and imbabala (T. s. sylvaticus) bushbuck. We found strong evidence for gene flow in at least four of eight nodes on the species tree. Among the six phenotypic traits assessed here, only habitat type mapped onto the species tree without homoplasy, showing that trait evolution was the result of complex patterns of divergence, introgression and convergent evolution.

A genomic assessment of the marine-speciation paradox within the toothed whale superfamily Delphinoidea.

The impact of post-divergence gene flow in speciation has been documented across a range of taxa in recent years, and may have been especially widespread in highly mobile, wide-ranging marine species, such as cetaceans. Here, we studied individual genomes from nine species across the three families of the toothed whale superfamily Delphinoidea (Delphinidae, Phocoenidae and Monodontidae). To investigate the role of post-divergence gene flow in the speciation process, we used a multifaceted approach, including (i) phylogenomics, (ii) the distribution of shared derived alleles and (iii) demographic inference. We found the divergence of lineages within Delphinoidea did not follow a process of pure bifurcation, but was much more complex. Sliding-window phylogenomics reveal a high prevalence of discordant topologies within the superfamily, with further analyses indicating these discordances arose due to both incomplete lineage sorting and gene flow. D-statistics and f-branch analyses supported gene flow between members of Delphinoidea, with the vast majority of gene flow occurring as ancient interfamilial events. Demographic analyses provided evidence that introgressive gene flow has likely ceased between all species pairs tested, despite reports of contemporary interspecific hybrids. Our study provides the first steps towards resolving the large complexity of speciation within Delphinoidea; we reveal the prevalence of ancient interfamilial gene flow events prior to the diversification of each family, and suggest that contemporary hybridisation events may be disadvantageous, as hybrid individuals do not appear to contribute to the parental species' gene pools.

Functional insights from the GC-poor genomes of two aphid parasitoids, Aphidius ervi and Lysiphlebus fabarum.

BACKGROUND: Parasitoid wasps have fascinating life cycles and play an important role in trophic networks, yet little is known about their genome content and function. Parasitoids that infect aphids are an important group with the potential for biological control. Their success depends on adapting to develop inside aphids and overcoming both host aphid defenses and their protective endosymbionts. RESULTS: We present the de novo genome assemblies, detailed annotation, and comparative analysis of two closely related parasitoid wasps that target pest aphids: Aphidius ervi and Lysiphlebus fabarum (Hymenoptera: Braconidae: Aphidiinae). The genomes are small (139 and 141 Mbp) and the most AT-rich reported thus far for any arthropod (GC content: 25.8 and 23.8%). This nucleotide bias is accompanied by skewed codon usage and is stronger in genes with adult-biased expression. AT-richness may be the consequence of reduced genome size, a near absence of DNA methylation, and energy efficiency. We identify missing desaturase genes, whose absence may underlie mimicry in the cuticular hydrocarbon profile of L. fabarum. We highlight key gene groups including those underlying venom composition, chemosensory perception, and sex determination, as well as potential losses in immune pathway genes. CONCLUSIONS: These findings are of fundamental interest for insect evolution and biological control applications. They provide a strong foundation for further functional studies into coevolution between parasitoids and their hosts. Both genomes are available at https://bipaa.genouest.org.

Novel Genes, Ancient Genes, and Gene Co-Option Contributed to the Genetic Basis of the Radula, a Molluscan Innovation.

The radula is the central foraging organ and apomorphy of the Mollusca. However, in contrast to other innovations, including the mollusk shell, genetic underpinnings of radula formation remain virtually unknown. Here, we present the first radula formative tissue transcriptome using the viviparous freshwater snail Tylomelania sarasinorum and compare it to foot tissue and the shell-building mantle of the same species. We combine differential expression, functional enrichment, and phylostratigraphic analyses to identify both specific and shared genetic underpinnings of the three tissues as well as their dominant functions and evolutionary origins. Gene expression of radula formative tissue is very distinct, but nevertheless more similar to mantle than to foot. Generally, the genetic bases of both radula and shell formation were shaped by novel orchestration of preexisting genes and continuous evolution of novel genes. A significantly increased proportion of radula-specific genes originated since the origin of stem-mollusks, indicating that novel genes were especially important for radula evolution. Genes with radula-specific expression in our study are frequently also expressed during the formation of other lophotrochozoan hard structures, like chaetae (hes1, arx), spicules (gbx), and shells of mollusks (gbx, heph) and brachiopods (heph), suggesting gene co-option for hard structure formation. Finally, a Lophotrochozoa-specific chitin synthase with a myosin motor domain (CS-MD), which is expressed during mollusk and brachiopod shell formation, had radula-specific expression in our study. CS-MD potentially facilitated the construction of complex chitinous structures and points at the potential of molecular novelties to promote the evolution of different morphological innovations.

Extended and Continuous Decline in Effective Population Size Results in Low Genomic Diversity in the World's Rarest Hyena Species, the Brown Hyena.

Hyenas (family Hyaenidae), as the sister group to cats (family Felidae), represent a deeply diverging branch within the cat-like carnivores (Feliformia). With an estimated population size of <10,000 individuals worldwide, the brown hyena (Parahyaena brunnea) represents the rarest of the four extant hyena species and has been listed as Near Threatened by the IUCN. Here, we report a high-coverage genome from a captive bred brown hyena and both mitochondrial and low-coverage nuclear genomes of 14 wild-caught brown hyena individuals from across southern Africa. We find that brown hyena harbor extremely low genetic diversity on both the mitochondrial and nuclear level, most likely resulting from a continuous and ongoing decline in effective population size that started ∼1 Ma and dramatically accelerated towards the end of the Pleistocene. Despite the strikingly low genetic diversity, we find no evidence of inbreeding within the captive bred individual and reveal phylogeographic structure, suggesting the existence of several potential subpopulations within the species.

Complex Admixture Preceded and Followed the Extinction of Wisent in the Wild.

Retracing complex population processes that precede extreme bottlenecks may be impossible using data from living individuals. The wisent (Bison bonasus), Europe's largest terrestrial mammal, exemplifies such a population history, having gone extinct in the wild but subsequently restored by captive breeding efforts. Using low coverage genomic data from modern and historical individuals, we investigate population processes occurring before and after this extinction. Analysis of aligned genomes supports the division of wisent into two previously recognized subspecies, but almost half of the genomic alignment contradicts this population history as a result of incomplete lineage sorting and admixture. Admixture between subspecies populations occurred prior to extinction and subsequently during the captive breeding program. Admixture with the Bos cattle lineage is also widespread but results from ancient events rather than recent hybridization with domestics. Our study demonstrates the huge potential of historical genomes for both studying evolutionary histories and for guiding conservation strategies.

Thermal and electrical transport across a magnetic quantum critical point.

A quantum critical point (QCP) arises when a continuous transition between competing phases occurs at zero temperature. Collective excitations at magnetic QCPs give rise to metallic properties that strongly deviate from the expectations of Landau's Fermi-liquid description, which is the standard theory of electron correlations in metals. Central to this theory is the notion of quasiparticles, electronic excitations that possess the quantum numbers of the non-interacting electrons. Here we report measurements of thermal and electrical transport across the field-induced magnetic QCP in the heavy-fermion compound YbRh(2)Si(2) (refs 2, 3). We show that the ratio of the thermal to electrical conductivities at the zero-temperature limit obeys the Wiedemann-Franz law for magnetic fields above the critical field at which the QCP is attained. This is also expected for magnetic fields below the critical field, where weak antiferromagnetic order and a Fermi-liquid phase form below 0.07 K (at zero field). At the critical field, however, the low-temperature electrical conductivity exceeds the thermal conductivity by about 10 per cent, suggestive of a non-Fermi-liquid ground state. This apparent violation of the Wiedemann-Franz law provides evidence for an unconventional type of QCP at which the fundamental concept of Landau quasiparticles no longer holds. These results imply that Landau quasiparticles break up, and that the origin of this disintegration is inelastic scattering associated with electronic quantum critical fluctuations--these insights could be relevant to understanding other deviations from Fermi-liquid behaviour frequently observed in various classes of correlated materials.

Decline in territory size and fecundity as a response to carrying capacity in an endangered songbird.

Density-dependent processes are fundamental mechanisms for the regulation of populations. Ecological theories differ in their predictions on whether increasing population density leads to individual adjustments of survival and reproductive output or to dominance and monopolization of resources. Here, we use a natural experiment to examine which factors limit population growth in the only remaining population of the endangered pale-headed brush finch (Atlapetes pallidiceps). For three distinct phases (a phase of population suppression, 2001-2002; expansion due to conservation management, 2003-2008; and equilibrium phase, 2009-2014), we estimated demographic parameters with an integrated population model using population size, the proportion of successfully breeding pairs and their productivity, territory size, and mark-recapture data of adult birds. A low proportion of successful breeders due to brood parasitism (0.42, 95% credible interval 0.26-0.59) limited population growth before 2003; subsequent culling of the brood parasite resulted in a two-fold increase of the proportion of successful breeders during the 'expansion phase'. When the population approached the carrying capacity of its habitat, territory size declined by more than 50% and fecundity declined from 1.9 (1.54-2.27) to 1.3 (1.12-1.53) chicks per breeding pair, but the proportion of successful breeders remained constant (expansion phase: 0.85; 0.76-0.93; equilibrium phase: 0.86; 0.79-0.92). This study demonstrates that limiting resources can lead to individual adjustments instead of despotic behavior, and the individual reduction of reproductive output at high population densities is consistent with the slow life-history of many tropical species.

Phylogenomic analyses unravel annelid evolution.

Annelida, the ringed worms, is a highly diverse animal phylum that includes more than 15,000 described species and constitutes the dominant benthic macrofauna from the intertidal zone down to the deep sea. A robust annelid phylogeny would shape our understanding of animal body-plan evolution and shed light on the bilaterian ground pattern. Traditionally, Annelida has been split into two major groups: Clitellata (earthworms and leeches) and polychaetes (bristle worms), but recent evidence suggests that other taxa that were once considered to be separate phyla (Sipuncula, Echiura and Siboglinidae (also known as Pogonophora)) should be included in Annelida. However, the deep-level evolutionary relationships of Annelida are still poorly understood, and a robust reconstruction of annelid evolutionary history is needed. Here we show that phylogenomic analyses of 34 annelid taxa, using 47,953 amino acid positions, recovered a well-supported phylogeny with strong support for major splits. Our results recover chaetopterids, myzostomids and sipunculids in the basal part of the tree, although the position of Myzostomida remains uncertain owing to its long branch. The remaining taxa are split into two clades: Errantia (which includes the model annelid Platynereis), and Sedentaria (which includes Clitellata). Ancestral character trait reconstructions indicate that these clades show adaptation to either an errant or a sedentary lifestyle, with alteration of accompanying morphological traits such as peristaltic movement, parapodia and sensory perception. Finally, life history characters in Annelida seem to be phylogenetically informative.

Endogenous murine leukemia retroviral variation across wild European and inbred strains of house mouse.

BACKGROUND: Endogenous murine leukemia retroviruses (MLVs) are high copy number proviral elements difficult to comprehensively characterize using standard low throughput sequencing approaches. However, high throughput approaches generate data that is challenging to process, interpret and present. RESULTS: Next generation sequencing (NGS) data was generated for MLVs from two wild caught Mus musculus domesticus (from mainland France and Corsica) and for inbred laboratory mouse strains C3H, LP/J and SJL. Sequence reads were grouped using a novel sequence clustering approach as applied to retroviral sequences. A Markov cluster algorithm was employed, and the sequence reads were queried for matches to specific xenotropic (Xmv), polytropic (Pmv) and modified polytropic (Mpmv) viral reference sequences. CONCLUSIONS: Various MLV subtypes were more widespread than expected among the mice, which may be due to the higher coverage of NGS, or to the presence of similar sequence across many different proviral loci. The results did not correlate with variation in the major MLV receptor Xpr1, which can restrict exogenous MLVs, suggesting that endogenous MLV distribution may reflect gene flow more than past resistance to infection.

The Arabidopsis Kinome: phylogeny and evolutionary insights into functional diversification.

BACKGROUND: Protein kinases constitute a particularly large protein family in Arabidopsis with important functions in cellular signal transduction networks. At the same time Arabidopsis is a model plant with high frequencies of gene duplications. Here, we have conducted a systematic analysis of the Arabidopsis kinase complement, the kinome, with particular focus on gene duplication events. We matched Arabidopsis proteins to a Hidden-Markov Model of eukaryotic kinases and computed a phylogeny of 942 Arabidopsis protein kinase domains and mapped their origin by gene duplication. RESULTS: The phylogeny showed two major clades of receptor kinases and soluble kinases, each of which was divided into functional subclades. Based on this phylogeny, association of yet uncharacterized kinases to families was possible which extended functional annotation of unknowns. Classification of gene duplications within these protein kinases revealed that representatives of cytosolic subfamilies showed a tendency to maintain segmentally duplicated genes, while some subfamilies of the receptor kinases were enriched for tandem duplicates. Although functional diversification is observed throughout most subfamilies, some instances of functional conservation among genes transposed from the same ancestor were observed. In general, a significant enrichment of essential genes was found among genes encoding for protein kinases. CONCLUSIONS: The inferred phylogeny allowed classification and annotation of yet uncharacterized kinases. The prediction and analysis of syntenic blocks and duplication events within gene families of interest can be used to link functional biology to insights from an evolutionary viewpoint. The approach undertaken here can be applied to any gene family in any organism with an annotated genome.

Genetic depletion at adaptive but not neutral loci in an endangered bird species.

Many endangered species suffer from the loss of genetic diversity, but some populations may be able to thrive even if genetically depleted. To investigate the underlying genetic processes of population bottlenecks, we apply an innovative approach for assessing genetic diversity in the last known population of the endangered Pale-headed Brushfinch (Atlapetes pallidiceps) in Ecuador. First, we measure genetic diversity at eleven neutral microsatellite loci and adaptive SNP variation in five Toll-like receptor (TLR) immune system genes. Bottleneck tests confirm genetic drift as the main force shaping genetic diversity in this species and indicate a 99 % reduction in population size dating back several hundred years. Second, we compare contemporary microsatellite diversity with historic museum samples of A. pallidiceps, finding no change in genetic diversity. Third, we compare genetic diversity in the Pale-headed Brushfinch with two co-occurring-related brushfinch species (Atlapetes latinuchus, Buarremon torquatus), finding a reduction of up to 91% diversity in the immune system genes but not in microsatellites. High TLR diversity is linked to decreased survival probabilities in A. pallidiceps. Low TLR diversity is thus probably an adaptation to the specific selection regime within its currently very restricted distribution (approximately 200 ha), but could severely restrict the adaptive potential of the species in the long run. Our study illustrates the importance of investigating both neutral and adaptive markers to assess the effect of population bottlenecks and for recommending specific management plans in endangered species.

Colonial-driven extinction of the blue antelope despite genomic adaptation to low population size.

Low genomic diversity is generally indicative of small population size and is considered detrimental by decreasing long-term adaptability.1,2,3,4,5,6 Moreover, small population size may promote gene flow with congeners and outbreeding depression.7,8,9,10,11,12,13 Here, we examine the connection between habitat availability, effective population size (Ne), and extinction by generating a 40× nuclear genome from the extinct blue antelope (Hippotragus leucophaeus). Historically endemic to the relatively small Cape Floristic Region in southernmost Africa,14,15 populations were thought to have expanded and contracted across glacial-interglacial cycles, tracking suitable habitat.16,17,18 However, we found long-term low Ne, unaffected by glacial cycles, suggesting persistence with low genomic diversity for many millennia prior to extinction in ∼AD 1800. A lack of inbreeding, alongside high levels of genetic purging, suggests adaptation to this long-term low Ne and that human impacts during the colonial era (e.g., hunting and landscape transformation), rather than longer-term ecological processes, were central to its extinction. Phylogenomic analyses uncovered gene flow between roan (H. equinus) and blue antelope, as well as between roan and sable antelope (H. niger), approximately at the time of divergence of blue and sable antelope (∼1.9 Ma). Finally, we identified the LYST and ASIP genes as candidates for the eponymous bluish pelt color of the blue antelope. Our results revise numerous aspects of our understanding of the interplay between genomic diversity and evolutionary history and provide the resources for uncovering the genetic basis of this extinct species' unique traits.

Single-Ion Kondo scaling of the coherent Fermi liquid regime in Ce(1-x)La(x)Ni2Ge2.

Thermodynamic and transport properties of the La-diluted Kondo lattice CeNi(2)Ge(2) were studied in a wide temperature range. The Ce-rich alloys Ce(1-x)La(x)Ni(2)Ge(2) were found to exhibit distinct features of the coherent heavy Fermi liquid. At intermediate compositions (0.7≤x≤0.9), non-Fermi liquid properties have been observed, followed by the local Fermi liquid behavior in the dilute limit. The 4f-electron contribution to the specific heat was found to follow the predictions of the Kondo-impurity model in both the local as well as the coherent regimes, with the characteristic Kondo temperature decreasing rapidly from about 30 K for the parent compound CeNi(2)Ge(2) to about 1 K in the most dilute samples. The specific heat does not show any evidence for the emergence of a new characteristic energy scale related to the formation of the coherent Kondo lattice.

Genome-scale phylogenetics: inferring the plant tree of life from 18,896 gene trees.

Phylogenetic analyses using genome-scale data sets must confront incongruence among gene trees, which in plants is exacerbated by frequent gene duplications and losses. Gene tree parsimony (GTP) is a phylogenetic optimization criterion in which a species tree that minimizes the number of gene duplications induced among a set of gene trees is selected. The run time performance of previous implementations has limited its use on large-scale data sets. We used new software that incorporates recent algorithmic advances to examine the performance of GTP on a plant data set consisting of 18,896 gene trees containing 510,922 protein sequences from 136 plant taxa (giving a combined alignment length of >2.9 million characters). The relationships inferred from the GTP analysis were largely consistent with previous large-scale studies of backbone plant phylogeny and resolved some controversial nodes. The placement of taxa that were present in few gene trees generally varied the most among GTP bootstrap replicates. Excluding these taxa either before or after the GTP analysis revealed high levels of phylogenetic support across plants. The analyses supported magnoliids sister to a eudicot + monocot clade and did not support the eurosid I and II clades. This study presents a nuclear genomic perspective on the broad-scale phylogenic relationships among plants, and it demonstrates that nuclear genes with a history of duplication and loss can be phylogenetically informative for resolving the plant tree of life.

Evolutionary Divergence and Radula Diversification in Two Ecomorphs from an Adaptive Radiation of Freshwater Snails.

(1) Background: Adaptive diversification of complex traits plays a pivotal role in the evolution of organismal diversity. In the freshwater snail genus Tylomelania, adaptive radiations were likely promoted by trophic specialization via diversification of their key foraging organ, the radula. (2) Methods: To investigate the molecular basis of radula diversification and its contribution to lineage divergence, we used tissue-specific transcriptomes of two sympatric Tylomelania sarasinorum ecomorphs. (3) Results: We show that ecomorphs are genetically divergent lineages with habitat-correlated abundances. Sequence divergence and the proportion of highly differentially expressed genes are significantly higher between radula transcriptomes compared to the mantle and foot. However, the same is not true when all differentially expressed genes or only non-synonymous SNPs are considered. Finally, putative homologs of some candidate genes for radula diversification (hh, arx, gbb) were also found to contribute to trophic specialization in cichlids and Darwin's finches. (4) Conclusions: Our results are in line with diversifying selection on the radula driving Tylomelania ecomorph divergence and indicate that some molecular pathways may be especially prone to adaptive diversification, even across phylogenetically distant animal groups.

Correction: Insight into the introduction of domestic cattle and the process of Neolithization to the Spanish region Galicia by genetic evidence.

[This corrects the article DOI: 10.1371/journal.pone.0249537.].

Collision risk of bats with small wind turbines: Worst-case scenarios near roosts, commuting and hunting structures.

Small wind turbines (SWTs) have become increasingly common within the last decade, but their impact on wildlife, especially bats, is largely unknown. We conducted an operational experiment by sequentially placing a mobile SWT with five different operational modes at six sites of high bat activity, including roosts, commuting structures, and highly frequented hunting areas. Bat flight trajectories around the SWT were documented at each site during five consecutive nights using a specifically designed high-spatial-resolution 3D camera. The recordings showed high bat activity levels close to the SWT (7,065 flight trajectories within a 10-m radius). The minimum distance to the rotor of each trajectory varied between 0 and 18 m, with a mean of 4.6 m across all sites. Linear mixed models created to account for site differences showed that, compared to a reference pole without a SWT, bats flew 0.4 m closer to the rotor (95% CI 0.3-0.6 m) if it was out of operation and 0.3 m closer (95% CI 0.1-0.4 m) if it was moving slowly. Exploratory behavior was frequently observed, with many bats deviating from their original flight trajectory to approach the rotor. Among 7,850 documented trajectories, 176 crossed the rotor, including 65 while it was in motion. The collision of one P. pygmaeus individual occurred during the experiment. These results demonstrate that, despite the generally strong ability of bats to evade moving rotor blades, bat casualties at SWTs placed at sites of high bat activity can reach or exceed the current threshold levels set for large wind turbines. As SWTs provide less energy than large turbines, their negative impact on bats should be minimized by avoidance measures such as a bat-friendly site selection or curtailment algorithms.

Insight into the introduction of domestic cattle and the process of Neolithization to the Spanish region Galicia by genetic evidence.

Domestic cattle were brought to Spain by early settlers and agricultural societies. Due to missing Neolithic sites in the Spanish region of Galicia, very little is known about this process in this region. We sampled 18 cattle subfossils from different ages and different mountain caves in Galicia, of which 11 were subject to sequencing of the mitochondrial genome and phylogenetic analysis, to provide insight into the introduction of cattle to this region. We detected high similarity between samples from different time periods and were able to compare the time frame of the first domesticated cattle in Galicia to data from the connecting region of Cantabria to show a plausible connection between the Neolithization of these two regions. Our data shows a close relationship of the early domesticated cattle of Galicia and modern cow breeds and gives a general insight into cattle phylogeny. We conclude that settlers migrated to this region of Spain from Europe and introduced common European breeds to Galicia.