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.

Transposon Mutagenesis Paired with Deep Sequencing of Caulobacter crescentus under Uranium Stress Reveals Genes Essential for Detoxification and Stress Tolerance.

UNLABELLED: The ubiquitous aquatic bacterium Caulobacter crescentus is highly resistant to uranium (U) and facilitates U biomineralization and thus holds promise as an agent of U bioremediation. To gain an understanding of how C. crescentus tolerates U, we employed transposon (Tn) mutagenesis paired with deep sequencing (Tn-seq) in a global screen for genomic elements required for U resistance. Of the 3,879 annotated genes in the C. crescentus genome, 37 were found to be specifically associated with fitness under U stress, 15 of which were subsequently tested through mutational analysis. Systematic deletion analysis revealed that mutants lacking outer membrane transporters (rsaFa and rsaFb), a stress-responsive transcription factor (cztR), or a ppGpp synthetase/hydrolase (spoT) exhibited a significantly lower survival rate under U stress. RsaFa and RsaFb, which are homologues of TolC in Escherichia coli, have previously been shown to mediate S-layer export. Transcriptional analysis revealed upregulation of rsaFa and rsaFb by 4- and 10-fold, respectively, in the presence of U. We additionally show that rsaFa mutants accumulated higher levels of U than the wild type, with no significant increase in oxidative stress levels. Our results suggest a function for RsaFa and RsaFb in U efflux and/or maintenance of membrane integrity during U stress. In addition, we present data implicating CztR and SpoT in resistance to U stress. Together, our findings reveal novel gene targets that are key to understanding the molecular mechanisms of U resistance in C. crescentus. IMPORTANCE: Caulobacter crescentus is an aerobic bacterium that is highly resistant to uranium (U) and has great potential to be used in U bioremediation, but its mechanisms of U resistance are poorly understood. We conducted a Tn-seq screen to identify genes specifically required for U resistance in C. crescentus. The genes that we identified have previously remained elusive using other omics approaches and thus provide significant insight into the mechanisms of U resistance by C. crescentus. In particular, we show that outer membrane transporters RsaFa and RsaFb, previously known as part of the S-layer export machinery, may confer U resistance by U efflux and/or by maintaining membrane integrity during U stress.

Rapid quantification of mutant fitness in diverse bacteria by sequencing randomly bar-coded transposons.

UNLABELLED: Transposon mutagenesis with next-generation sequencing (TnSeq) is a powerful approach to annotate gene function in bacteria, but existing protocols for TnSeq require laborious preparation of every sample before sequencing. Thus, the existing protocols are not amenable to the throughput necessary to identify phenotypes and functions for the majority of genes in diverse bacteria. Here, we present a method, random bar code transposon-site sequencing (RB-TnSeq), which increases the throughput of mutant fitness profiling by incorporating random DNA bar codes into Tn5 and mariner transposons and by using bar code sequencing (BarSeq) to assay mutant fitness. RB-TnSeq can be used with any transposon, and TnSeq is performed once per organism instead of once per sample. Each BarSeq assay requires only a simple PCR, and 48 to 96 samples can be sequenced on one lane of an Illumina HiSeq system. We demonstrate the reproducibility and biological significance of RB-TnSeq with Escherichia coli, Phaeobacter inhibens, Pseudomonas stutzeri, Shewanella amazonensis, and Shewanella oneidensis. To demonstrate the increased throughput of RB-TnSeq, we performed 387 successful genome-wide mutant fitness assays representing 130 different bacterium-carbon source combinations and identified 5,196 genes with significant phenotypes across the five bacteria. In P. inhibens, we used our mutant fitness data to identify genes important for the utilization of diverse carbon substrates, including a putative d-mannose isomerase that is required for mannitol catabolism. RB-TnSeq will enable the cost-effective functional annotation of diverse bacteria using mutant fitness profiling. IMPORTANCE: A large challenge in microbiology is the functional assessment of the millions of uncharacterized genes identified by genome sequencing. Transposon mutagenesis coupled to next-generation sequencing (TnSeq) is a powerful approach to assign phenotypes and functions to genes. However, the current strategies for TnSeq are too laborious to be applied to hundreds of experimental conditions across multiple bacteria. Here, we describe an approach, random bar code transposon-site sequencing (RB-TnSeq), which greatly simplifies the measurement of gene fitness by using bar code sequencing (BarSeq) to monitor the abundance of mutants. We performed 387 genome-wide fitness assays across five bacteria and identified phenotypes for over 5,000 genes. RB-TnSeq can be applied to diverse bacteria and is a powerful tool to annotate uncharacterized genes using phenotype data.

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.

Genome-wide nucleosome map and cytosine methylation levels of an ancient human genome.

Epigenetic information is available from contemporary organisms, but is difficult to track back in evolutionary time. Here, we show that genome-wide epigenetic information can be gathered directly from next-generation sequence reads of DNA isolated from ancient remains. Using the genome sequence data generated from hair shafts of a 4000-yr-old Paleo-Eskimo belonging to the Saqqaq culture, we generate the first ancient nucleosome map coupled with a genome-wide survey of cytosine methylation levels. The validity of both nucleosome map and methylation levels were confirmed by the recovery of the expected signals at promoter regions, exon/intron boundaries, and CTCF sites. The top-scoring nucleosome calls revealed distinct DNA positioning biases, attesting to nucleotide-level accuracy. The ancient methylation levels exhibited high conservation over time, clustering closely with modern hair tissues. Using ancient methylation information, we estimated the age at death of the Saqqaq individual and illustrate how epigenetic information can be used to infer ancient gene expression. Similar epigenetic signatures were found in other fossil material, such as 110,000- to 130,000-yr-old bones, supporting the contention that ancient epigenomic information can be reconstructed from a deep past. Our findings lay the foundation for extracting epigenomic information from ancient samples, allowing shifts in epialleles to be tracked through evolutionary time, as well as providing an original window into modern epigenomics.

A scalable, fully automated process for construction of sequence-ready barcoded libraries for 454.

We present an automated, high throughput library construction process for 454 technology. Sample handling errors and cross-contamination are minimized via end-to-end barcoding of plasticware, along with molecular DNA barcoding of constructs. Automation-friendly magnetic bead-based size selection and cleanup steps have been devised, eliminating major bottlenecks and significant sources of error. Using this methodology, one technician can create 96 sequence-ready 454 libraries in 2 days, a dramatic improvement over the standard method.

Transcriptome and metabolite responses to predation in a South pacific soft coral.

Sinularia polydactyla, a dioecious, abundant soft coral in the South Pacific, exhibits biochemical phenotypic plasticity in secondary metabolite production in relation to predation intensity. However, it is unclear to what extent changes in secondary metabolites, such as 11beta-acetoxypukalide, may result from specific, induced pathway activities at the level of gene expression. To investigate both chemical changes and differences in mRNA diversity in response to predation stress, artificial predation experiments were conducted in situ on colonies of S. polydactyla. Multivariate statistical analyses of coral biochemical metabolites and our kinetic transcriptome profiling technique indicate that that the induction of 11beta-acetoxypukalide by predation stress likely results from the upregulation of either one dominant transcript or a very small set of transcripts, indicative of a targeted upregulation rather than a generalized, genetic stress response. Overall, this work establishes a routine method for integrating high-throughput transcriptome and metabolome data sets to allow for the identification of metabolites whose intracellular concentrations can be readily linked to gene expression events in response to specific treatments in non-model organisms.

Fundulus as the premier teleost model in environmental biology: opportunities for new insights using genomics.

A strong foundation of basic and applied research documents that the estuarine fish Fundulus heteroclitus and related species are unique laboratory and field models for understanding how individuals and populations interact with their environment. In this paper we summarize an extensive body of work examining the adaptive responses of Fundulus species to environmental conditions, and describe how this research has contributed importantly to our understanding of physiology, gene regulation, toxicology, and ecological and evolutionary genetics of teleosts and other vertebrates. These explorations have reached a critical juncture at which advancement is hindered by the lack of genomic resources for these species. We suggest that a more complete genomics toolbox for F. heteroclitus and related species will permit researchers to exploit the power of this model organism to rapidly advance our understanding of fundamental biological and pathological mechanisms among vertebrates, as well as ecological strategies and evolutionary processes common to all living organisms.

A functional approach to transcriptome profiling: linking gene expression patterns to metabolites that matter.

Secondary metabolites or natural products have been isolated from many marine organisms. These metabolites often have important bioactive functions; however, very little information is available regarding the biosynthesis and regulation of many secondary metabolites. At a time when use of marine-derived metabolites is rapidly expanding in industry and pharmacological fields, a better understanding of the genetic mechanisms controlling secondary metabolite production is necessary. We review the recent development of a novel transcriptome profiling methodology that allows for rapid and high-throughput screening of changes in mRNA sequence pools. The application of genomics-based techniques and the integration of both biochemical and molecular data sets in marine organisms complement ongoing drug discovery efforts.

Profiling transcriptome complexity and secondary metabolite synthesis in a benthic soft coral, Sinularia polydactyla.

Sinularia polydactyla, an abundant Indo-Pacific soft coral species, exhibits biochemical phenotypic plasticity, prompting investigations into differences in mRNA diversity and complexity in response to predation stress. Changes in transcriptome complexity of S. polydactyla cDNA libraries were measured using reannealing rate assays that employ an informatics-based analysis of kinetic profiles. This method allows for quick, high-throughput analysis of sequence complexity and has been used to compare transcriptome-level differences in other marine invertebrates. Here, S. polydactyla colonies were transplanted between two sites exhibiting high and low predation levels. Statistically significant differences between bite scar counts found on different transplant groups suggest site-specific variation in predation. Changes in mRNA pool complexity were quantified to indicate shifts in secondary metabolite concentration between treatment groups. Examining the complexity of the mRNA pool in this soft coral is one of the first steps toward understanding the mechanisms of phenotypic plasticity at a biochemical and molecular level.

Gene expression profiling of two related soft corals, Sinularia polydactyla and S. maxima, and their putative hybrid at different life-history stages.

The soft coral species, Sinularia polydactyla and Sinularia maxima are abundant in shallow coral reef habitats in Guam. Both species are gonochoric and mass spawning events can result in an apparent S. polydactyla-S. maxima hybrid exhibiting morphological characteristics of both parents. These morphological differences prompted our investigation of potential differences in gene expression patterns among these closely related species. In concert with interspecific level differences, intraspecific level differences in transcriptome diversity and complexity were also studied among juvenile, adult male, and adult female S. polydactyla colonies. To uncover these transcriptome-level differences, RNA was extracted from RNAlater(R)-preserved samples and cDNA libraries constructed using a nano-scale synthesis strategy. Changes in transcriptome complexity (mRNA sequence composition) of these libraries were measured using reannealing rate assays that employ an informatics-based analysis of kinetic profiles. This method allows for quick, high-throughput analysis of sequence complexity and has been used to compare transcriptome-level differences in other marine invertebrates. Comparisons of transcriptome data revealed diagnostic differences in the transcriptome pools. Examination of transcriptome changes in closely related species experiencing similar environmental parameters may enable us to understand developmental, morphological, and sex-specific changes in gene expression patterns.