Tracking Five Millennia of Horse Management with Extensive Ancient Genome Time Series.

Horse domestication revolutionized warfare and accelerated travel, trade, and the geographic expansion of languages. Here, we present the largest DNA time series for a non-human organism to date, including genome-scale data from 149 ancient animals and 129 ancient genomes (≥1-fold coverage), 87 of which are new. This extensive dataset allows us to assess the modern legacy of past equestrian civilizations. We find that two extinct horse lineages existed during early domestication, one at the far western (Iberia) and the other at the far eastern range (Siberia) of Eurasia. None of these contributed significantly to modern diversity. We show that the influence of Persian-related horse lineages increased following the Islamic conquests in Europe and Asia. Multiple alleles associated with elite-racing, including at the MSTN "speed gene," only rose in popularity within the last millennium. Finally, the development of modern breeding impacted genetic diversity more dramatically than the previous millennia of human management.

The origins and spread of domestic horses from the Western Eurasian steppes.

Domestication of horses fundamentally transformed long-range mobility and warfare1. However, modern domesticated breeds do not descend from the earliest domestic horse lineage associated with archaeological evidence of bridling, milking and corralling2-4 at Botai, Central Asia around 3500 BC3. Other longstanding candidate regions for horse domestication, such as Iberia5 and Anatolia6, have also recently been challenged. Thus, the genetic, geographic and temporal origins of modern domestic horses have remained unknown. Here we pinpoint the Western Eurasian steppes, especially the lower Volga-Don region, as the homeland of modern domestic horses. Furthermore, we map the population changes accompanying domestication from 273 ancient horse genomes. This reveals that modern domestic horses ultimately replaced almost all other local populations as they expanded rapidly across Eurasia from about 2000 BC, synchronously with equestrian material culture, including Sintashta spoke-wheeled chariots. We find that equestrianism involved strong selection for critical locomotor and behavioural adaptations at the GSDMC and ZFPM1 genes. Our results reject the commonly held association7 between horseback riding and the massive expansion of Yamnaya steppe pastoralists into Europe around 3000 BC8,9 driving the spread of Indo-European languages10. This contrasts with the scenario in Asia where Indo-Iranian languages, chariots and horses spread together, following the early second millennium BC Sintashta culture11,12.

Decryption of the survival "black box": gene family expansion promotes the encystment in ciliated protists.

BACKGROUND: Encystment is an important survival strategy extensively employed by microbial organisms to survive unfavorable conditions. Single-celled ciliated protists (ciliates) are popular model eukaryotes for studying encystment, whereby these cells degenerate their ciliary structures and develop cyst walls, then reverse the process under more favorable conditions. However, to date, the evolutionary basis and mechanism for encystment in ciliates is largely unknown. With the rapid development of high-throughput sequencing technologies, genome sequencing and comparative genomics of ciliates have become effective methods to provide insights into above questions. RESULTS: Here, we profiled the MAC genome of Pseudourostyla cristata, a model hypotrich ciliate for encystment studies. Like other hypotrich MAC genomes, the P. cristata MAC genome is extremely fragmented with a single gene on most chromosomes, and encodes introns that are generally small and lack a conserved branch point for pre-mRNA splicing. Gene family expansion analyses indicate that multiple gene families involved in the encystment are expanded during the evolution of P. cristata. Furthermore, genomic comparisons with other five representative hypotrichs indicate that gene families of phosphorelay sensor kinase, which play a role in the two-component signal transduction system that is related to encystment, show significant expansion among all six hypotrichs. Additionally, cyst wall-related chitin synthase genes have experienced structural changes that increase them from single-exon to multi-exon genes during evolution. These genomic features potentially promote the encystment in hypotrichs and enhance their ability to survive in adverse environments during evolution. CONCLUSIONS: We systematically investigated the genomic structure of hypotrichs and key evolutionary phenomenon, gene family expansion, for encystment promotion in ciliates. In summary, our results provided insights into the evolutionary mechanism of encystment in ciliates.

Tobacco leaf tissue rapidly detoxifies direct salt loads without activation of calcium and SOS signaling.

Salt stress is a major abiotic stress, responsible for declining agricultural productivity. Roots are regarded as hubs for salt detoxification, however, leaf salt concentrations may exceed those of roots. How mature leaves manage acute sodium chloride (NaCl) stress is mostly unknown. To analyze the mechanisms for NaCl redistribution in leaves, salt was infiltrated into intact tobacco leaves. It initiated pronounced osmotically-driven leaf movements. Leaf downward movement caused by hydro-passive turgor loss reached a maximum within 2 h. Salt-driven cellular water release was accompanied by a transient change in membrane depolarization but not an increase in cytosolic calcium ion (Ca2+ ) level. Nonetheless, only half an hour later, the leaves had completely regained turgor. This recovery phase was characterized by an increase in mesophyll cell plasma membrane hydrogen ion (H+ ) pumping, a salt uptake-dependent cytosolic alkalization, and a return of the apoplast osmolality to pre-stress levels. Although, transcript numbers of abscisic acid- and Salt Overly Sensitive pathway elements remained unchanged, salt adaptation depended on the vacuolar H+ /Na+ -exchanger NHX1. Altogether, tobacco leaves can detoxify sodium ions (Na+ ) rapidly even under massive salt loads, based on pre-established posttranslational settings and NHX1 cation/H+ antiport activity. Unlike roots, signaling and processing of salt stress in tobacco leaves does not depend on Ca2+ signaling.

Three closely-related subclasses Phacodiniidia Small & Lynn, 1985, Protohypotrichia Shi et al., 1999, and Euplotia Jankowski, 1979 (Protista, Ciliophora): A new contribution to their phylogeny with reconsiderations on the evolutionary hypotheses.

The huge variety of species and worldwide distribution of ciliated protists in class Spirotrichea continue to make it one of the most complicated and confused groups in Ciliophora, despite significant research interest in the unique molecular genetics of these organisms. In this study, the morphological and molecular information were integrated, and it is inferred from a new perspective for the evolutionary relationship among Phacodiniidia, Protohypotrichia, Hypotrichia and Euplotia. Our results indicate that Kiitricha and Caryotricha, two members in Protohypotrichia, may represent two parallel branches of evolution; Euplotidae and Aspidiscidae represent the most recently diverged taxa within Euplotida, followed by Certesiidae, Gastrocirrhidae, and Uronychidae. Further, representative morphological characters (e.g. fronto-ventral-transverse cirral anlagen, undulating membranes, marginal cirri and caudal cirri) were stochastically mapped on phylogenies to speculate evolutionary path and morphological characters of the evolutionary transition node groups were assumed.

The Venus flytrap trigger hair-specific potassium channel KDM1 can reestablish the K+ gradient required for hapto-electric signaling.

The carnivorous plant Dionaea muscipula harbors multicellular trigger hairs designed to sense mechanical stimuli upon contact with animal prey. At the base of the trigger hair, mechanosensation is transduced into an all-or-nothing action potential (AP) that spreads all over the trap, ultimately leading to trap closure and prey capture. To reveal the molecular basis for the unique functional repertoire of this mechanoresponsive plant structure, we determined the transcriptome of D. muscipula's trigger hair. Among the genes that were found to be highly specific to the trigger hair, the Shaker-type channel KDM1 was electrophysiologically characterized as a hyperpolarization- and acid-activated K+-selective channel, thus allowing the reuptake of K+ ions into the trigger hair's sensory cells during the hyperpolarization phase of the AP. During trap development, the increased electrical excitability of the trigger hair is associated with the transcriptional induction of KDM1. Conversely, when KDM1 is blocked by Cs+ in adult traps, the initiation of APs in response to trigger hair deflection is reduced, and trap closure is suppressed. KDM1 thus plays a dominant role in K+ homeostasis in the context of AP and turgor formation underlying the mechanosensation of trigger hair cells and thus D. muscipula's hapto-electric signaling.

Stalk cell polar ion transport provide for bladder-based salinity tolerance in Chenopodium quinoa.

Chenopodium quinoa uses epidermal bladder cells (EBCs) to sequester excess salt. Each EBC complex consists of a leaf epidermal cell, a stalk cell, and the bladder. Under salt stress, sodium (Na+ ), chloride (Cl- ), potassium (K+ ) and various metabolites are shuttled from the leaf lamina to the bladders. Stalk cells operate as both a selectivity filter and a flux controller. In line with the nature of a transfer cell, advanced transmission electron tomography, electrophysiology, and fluorescent tracer flux studies revealed the stalk cell's polar organization and bladder-directed solute flow. RNA sequencing and cluster analysis revealed the gene expression profiles of the stalk cells. Among the stalk cell enriched genes, ion channels and carriers as well as sugar transporters were most pronounced. Based on their electrophysiological fingerprint and thermodynamic considerations, a model for stalk cell transcellular transport was derived.

New contribution to the peritrichous genus Ophrydium (Protista, Ciliophora) with notes on the morphology, taxonomy, and phylogeny of a well-known species Ophrydium crassicaule Penard, 1922.

The peritrichous genus Ophrydium is relatively distinctive since its colonies have a gelatinous lorica. In the present work, one morphologically well-described species, Ophrydium crassicaule, was collected from a freshwater habitat in China. Both morphological characters and sequence data for SSU rDNA, ITS1-5.8S-ITS2 region and LSU rDNA were obtained. Morphologically, O. crassicaule is characterized by its elongated vase-shaped zooid, single-layered peristomial lip, dichotomously branched stalk, colony with gelatinous lorica, and infundibular polykinety 3 (P3) containing one short inner row and two long outer rows. Phylogenetic analyses revealed that the genus Ophrydium is monophyletic and nests within the family Vorticellidae. In addition, we briefly revise the nominal congeners for aiding the species identification in future studies.

Venus flytrap carnivorous lifestyle builds on herbivore defense strategies.

Although the concept of botanical carnivory has been known since Darwin's time, the molecular mechanisms that allow animal feeding remain unknown, primarily due to a complete lack of genomic information. Here, we show that the transcriptomic landscape of the Dionaea trap is dramatically shifted toward signal transduction and nutrient transport upon insect feeding, with touch hormone signaling and protein secretion prevailing. At the same time, a massive induction of general defense responses is accompanied by the repression of cell death-related genes/processes. We hypothesize that the carnivory syndrome of Dionaea evolved by exaptation of ancient defense pathways, replacing cell death with nutrient acquisition.

Redescription of a Hymenostome Ciliate, Tetrahymena setosa (Protozoa, Ciliophora) Notes on its Molecular Phylogeny.

In recent years, Tetrahymena species have been used as model organisms for research in a wide range of fields, highlighting the need for a fuller understanding of the taxonomy of this group. It is in this context that this paper uses living observation and silver staining methods to investigate the morphology and infraciliature of one Tetrahymena species, T. setosa (Schewiakoff 1892 Verh. Naturh. Med. Ver. Heidelb., 4:544) McCoy (1975) Acta Protozool., 14:253; the senior subjective synonym of T. setifera Holz and Corliss (1956) J. Protozool., 3:112; isolated from a freshwater pond in Harbin, north-eastern China. This organism can be distinguished from other described Tetrahymena species mainly by its single caudal cilium, which is about twice the length of the somatic ciliature. While the Harbin isolate appears similar to the population described by Holz and Corliss (1956) J. Protozool., 3:112, an improved diagnosis for T. setosa is given based on the previous descriptions and the Harbin population. In summary, this species can be recognized mainly by the combination of the following characters: body in vivo approximately 40 µm × 25 µm, 21-26 somatic kineties, one to four contractile vacuole pores associated with meridians 6-11 and a single caudal cilium. The small subunit ribosomal (SSU) rRNA gene and the cox1 gene sequences of Harbin population are also characterized in order to corroborate that the isolated species branches in phylogenetic trees as a T. setosa species. The phylogenetic analysis also indicated that sequences of populations of Tetrahymena species should be published with detailed morphological identifications.

Recalibrating Equus evolution using the genome sequence of an early Middle Pleistocene horse.

The rich fossil record of equids has made them a model for evolutionary processes. Here we present a 1.12-times coverage draft genome from a horse bone recovered from permafrost dated to approximately 560-780 thousand years before present (kyr BP). Our data represent the oldest full genome sequence determined so far by almost an order of magnitude. For comparison, we sequenced the genome of a Late Pleistocene horse (43 kyr BP), and modern genomes of five domestic horse breeds (Equus ferus caballus), a Przewalski's horse (E. f. przewalskii) and a donkey (E. asinus). Our analyses suggest that the Equus lineage giving rise to all contemporary horses, zebras and donkeys originated 4.0-4.5 million years before present (Myr BP), twice the conventionally accepted time to the most recent common ancestor of the genus Equus. We also find that horse population size fluctuated multiple times over the past 2 Myr, particularly during periods of severe climatic changes. We estimate that the Przewalski's and domestic horse populations diverged 38-72 kyr BP, and find no evidence of recent admixture between the domestic horse breeds and the Przewalski's horse investigated. This supports the contention that Przewalski's horses represent the last surviving wild horse population. We find similar levels of genetic variation among Przewalski's and domestic populations, indicating that the former are genetically viable and worthy of conservation efforts. We also find evidence for continuous selection on the immune system and olfaction throughout horse evolution. Finally, we identify 29 genomic regions among horse breeds that deviate from neutrality and show low levels of genetic variation compared to the Przewalski's horse. Such regions could correspond to loci selected early during domestication.

Evolutionary Patterns and Processes: Lessons from Ancient DNA.

Ever since its emergence in 1984, the field of ancient DNA has struggled to overcome the challenges related to the decay of DNA molecules in the fossil record. With the recent development of high-throughput DNA sequencing technologies and molecular techniques tailored to ultra-damaged templates, it has now come of age, merging together approaches in phylogenomics, population genomics, epigenomics, and metagenomics. Leveraging on complete temporal sample series, ancient DNA provides direct access to the most important dimension in evolution­time, allowing a wealth of fundamental evolutionary processes to be addressed at unprecedented resolution. This review taps into the most recent findings in ancient DNA research to present analyses of ancient genomic and metagenomic data.

Calcium sensor kinase activates potassium uptake systems in gland cells of Venus flytraps.

The Darwin plant Dionaea muscipula is able to grow on mineral-poor soil, because it gains essential nutrients from captured animal prey. Given that no nutrients remain in the trap when it opens after the consumption of an animal meal, we here asked the question of how Dionaea sequesters prey-derived potassium. We show that prey capture triggers expression of a K(+) uptake system in the Venus flytrap. In search of K(+) transporters endowed with adequate properties for this role, we screened a Dionaea expressed sequence tag (EST) database and identified DmKT1 and DmHAK5 as candidates. On insect and touch hormone stimulation, the number of transcripts of these transporters increased in flytraps. After cRNA injection of K(+)-transporter genes into Xenopus oocytes, however, both putative K(+) transporters remained silent. Assuming that calcium sensor kinases are regulating Arabidopsis K(+) transporter 1 (AKT1), we coexpressed the putative K(+) transporters with a large set of kinases and identified the CBL9-CIPK23 pair as the major activating complex for both transporters in Dionaea K(+) uptake. DmKT1 was found to be a K(+)-selective channel of voltage-dependent high capacity and low affinity, whereas DmHAK5 was identified as the first, to our knowledge, proton-driven, high-affinity potassium transporter with weak selectivity. When the Venus flytrap is processing its prey, the gland cell membrane potential is maintained around -120 mV, and the apoplast is acidified to pH 3. These conditions in the green stomach formed by the closed flytrap allow DmKT1 and DmHAK5 to acquire prey-derived K(+), reducing its concentration from millimolar levels down to trace levels.

Tracking the origins of Yakutian horses and the genetic basis for their fast adaptation to subarctic environments.

Yakutia, Sakha Republic, in the Siberian Far East, represents one of the coldest places on Earth, with winter record temperatures dropping below -70 °C. Nevertheless, Yakutian horses survive all year round in the open air due to striking phenotypic adaptations, including compact body conformations, extremely hairy winter coats, and acute seasonal differences in metabolic activities. The evolutionary origins of Yakutian horses and the genetic basis of their adaptations remain, however, contentious. Here, we present the complete genomes of nine present-day Yakutian horses and two ancient specimens dating from the early 19th century and ∼5,200 y ago. By comparing these genomes with the genomes of two Late Pleistocene, 27 domesticated, and three wild Przewalski's horses, we find that contemporary Yakutian horses do not descend from the native horses that populated the region until the mid-Holocene, but were most likely introduced following the migration of the Yakut people a few centuries ago. Thus, they represent one of the fastest cases of adaptation to the extreme temperatures of the Arctic. We find cis-regulatory mutations to have contributed more than nonsynonymous changes to their adaptation, likely due to the comparatively limited standing variation within gene bodies at the time the population was founded. Genes involved in hair development, body size, and metabolic and hormone signaling pathways represent an essential part of the Yakutian horse adaptive genetic toolkit. Finally, we find evidence for convergent evolution with native human populations and woolly mammoths, suggesting that only a few evolutionary strategies are compatible with survival in extremely cold environments.

The desert plant Phoenix dactylifera closes stomata via nitrate-regulated SLAC1 anion channel.

Date palm Phoenix dactylifera is a desert crop well adapted to survive and produce fruits under extreme drought and heat. How are palms under such harsh environmental conditions able to limit transpirational water loss? Here, we analysed the cuticular waxes, stomata structure and function, and molecular biology of guard cells from P. dactylifera. To understand the stomatal response to the water stress phytohormone of the desert plant, we cloned the major elements necessary for guard cell fast abscisic acid (ABA) signalling and reconstituted this ABA signalosome in Xenopus oocytes. The PhoenixSLAC1-type anion channel is regulated by ABA kinase PdOST1. Energy-dispersive X-ray analysis (EDXA) demonstrated that date palm guard cells release chloride during stomatal closure. However, in Cl- medium, PdOST1 did not activate the desert plant anion channel PdSLAC1 per se. Only when nitrate was present at the extracellular face of the anion channel did the OST1-gated PdSLAC1 open, thus enabling chloride release. In the presence of nitrate, ABA enhanced and accelerated stomatal closure. Our findings indicate that, in date palm, the guard cell osmotic motor driving stomatal closure uses nitrate as the signal to open the major anion channel SLAC1. This initiates guard cell depolarization and the release of anions together with potassium.

A Single-Pore Residue Renders the Arabidopsis Root Anion Channel SLAH2 Highly Nitrate Selective.

In contrast to animal cells, plants use nitrate as a major source of nitrogen. Following the uptake of nitrate, this major macronutrient is fed into the vasculature for long-distance transport. The Arabidopsis thaliana shoot expresses the anion channel SLOW ANION CHANNEL1 (SLAC1) and its homolog SLAC1 HOMOLOGOUS3 (SLAH3), which prefer nitrate as substrate but cannot exclude chloride ions. By contrast, we identified SLAH2 as a nitrate-specific channel that is impermeable for chloride. To understand the molecular basis for nitrate selection in the SLAH2 channel, SLAC1 and SLAH2 were modeled to the structure of HiTehA, a distantly related bacterial member. Structure-guided site-directed mutations converted SLAC1 into a SLAH2-like nitrate-specific anion channel and vice versa. Our findings indicate that two pore-occluding phenylalanines constrict the pore. The selectivity filter of SLAC/SLAH anion channels is determined by the polarity of pore-lining residues located on alpha helix 3. Changing the polar character of a single amino acid side chain (Ser-228) to a nonpolar residue turned the nitrate-selective SLAH2 into a chloride/nitrate-permeable anion channel. Thus, the molecular basis of the anion specificity of SLAC/SLAH anion channels seems to be determined by the presence and constellation of polar side chains that act in concert with the two pore-occluding phenylalanines.

Morphology and Phylogenetic Placement of Three New Zoothamnium species (Ciliophora: Peritrichia) from Coastal Waters of Southern China.

The morphology, infraciliature, and silverline system of three peritrichous ciliates, Zoothamnium bucciniiformum sp. n., Zoothamnium florens sp. n., and Zoothamnium zhanjiangense sp. n., were investigated based on both living and silver-stained specimens. Zoothamnium bucciniiformum sp. n., collected from coastal waters (salinity 30‰) off Zhanjiang, southern China, can be distinguished by the following characters: dichotomously branched stalk, peristomial lip with medial circumferential infolding, contractile vacuole apically positioned, 32-49 silverlines between the anterior end and the aboral trochal band, 15-26 between the aboral trochal band and the scopula; two kineties in peniculus 3, not parallel to each other. Zoothamnium florens sp. n., collected from a mangrove wetland (salinity 13‰) off Zhanjiang, is characterized by its large conical zooid, tuberculate peristomial lip, asymmetrical dichotomously branched colony, 59-81 silverlines between the anterior end and the aboral trochal band and 29-36 between the aboral trochal band and the scopula. Zoothamnium zhanjiangense, collected from a mangrove wetland (salinity about 9.5‰) off Zhanjiang, differs from its congeners by the alternately branched stalk, peristomial lip with medial circumferential infolding, 40-63 silverlines from the peristomial area to the aboral trochal band and 13-24 from the aboral trochal band to the scopula. The comparison and analysis of SSU rDNA sequences also support present identifications.

Beyond the "Code": A Guide to the Description and Documentation of Biodiversity in Ciliated Protists (Alveolata, Ciliophora).

Recent advances in molecular technology have revolutionized research on all aspects of the biology of organisms, including ciliates, and created unprecedented opportunities for pursuing a more integrative approach to investigations of biodiversity. However, this goal is complicated by large gaps and inconsistencies that still exist in the foundation of basic information about biodiversity of ciliates. The present paper reviews issues relating to the taxonomy of ciliates and presents specific recommendations for best practice in the observation and documentation of their biodiversity. This effort stems from a workshop that explored ways to implement six Grand Challenges proposed by the International Research Coordination Network for Biodiversity of Ciliates (IRCN-BC). As part of its commitment to strengthening the knowledge base that supports research on biodiversity of ciliates, the IRCN-BC proposes to populate The Ciliate Guide, an online database, with biodiversity-related data and metadata to create a resource that will facilitate accurate taxonomic identifications and promote sharing of data.

Speciation with gene flow in equids despite extensive chromosomal plasticity.

Horses, asses, and zebras belong to a single genus, Equus, which emerged 4.0-4.5 Mya. Although the equine fossil record represents a textbook example of evolution, the succession of events that gave rise to the diversity of species existing today remains unclear. Here we present six genomes from each living species of asses and zebras. This completes the set of genomes available for all extant species in the genus, which was hitherto represented only by the horse and the domestic donkey. In addition, we used a museum specimen to characterize the genome of the quagga zebra, which was driven to extinction in the early 1900s. We scan the genomes for lineage-specific adaptations and identify 48 genes that have evolved under positive selection and are involved in olfaction, immune response, development, locomotion, and behavior. Our extensive genome dataset reveals a highly dynamic demographic history with synchronous expansions and collapses on different continents during the last 400 ky after major climatic events. We show that the earliest speciation occurred with gene flow in Northern America, and that the ancestor of present-day asses and zebras dispersed into the Old World 2.1-3.4 Mya. Strikingly, we also find evidence for gene flow involving three contemporary equine species despite chromosomal numbers varying from 16 pairs to 31 pairs. These findings challenge the claim that the accumulation of chromosomal rearrangements drive complete reproductive isolation, and promote equids as a fundamental model for understanding the interplay between chromosomal structure, gene flow, and, ultimately, speciation.

Prehistoric genomes reveal the genetic foundation and cost of horse domestication.

The domestication of the horse ∼ 5.5 kya and the emergence of mounted riding, chariotry, and cavalry dramatically transformed human civilization. However, the genetics underlying horse domestication are difficult to reconstruct, given the near extinction of wild horses. We therefore sequenced two ancient horse genomes from Taymyr, Russia (at 7.4- and 24.3-fold coverage), both predating the earliest archeological evidence of domestication. We compared these genomes with genomes of domesticated horses and the wild Przewalski's horse and found genetic structure within Eurasia in the Late Pleistocene, with the ancient population contributing significantly to the genetic variation of domesticated breeds. We furthermore identified a conservative set of 125 potential domestication targets using four complementary scans for genes that have undergone positive selection. One group of genes is involved in muscular and limb development, articular junctions, and the cardiac system, and may represent physiological adaptations to human utilization. A second group consists of genes with cognitive functions, including social behavior, learning capabilities, fear response, and agreeableness, which may have been key for taming horses. We also found that domestication is associated with inbreeding and an excess of deleterious mutations. This genetic load is in line with the "cost of domestication" hypothesis also reported for rice, tomatoes, and dogs, and it is generally attributed to the relaxation of purifying selection resulting from the strong demographic bottlenecks accompanying domestication. Our work demonstrates the power of ancient genomes to reconstruct the complex genetic changes that transformed wild animals into their domesticated forms, and the population context in which this process took place.