Ancestral sequence reconstruction (ASR) is a powerful tool to study the evolution of proteins and thus gain deep insight into the relationships among protein sequence, structure, and function. A major barrier to its broad use is the complexity of the task: it requires multiple software packages, complex file manipulations, and expert phylogenetic knowledge. Here we introduce topiary, a software pipeline that aims to overcome this barrier. To use topiary, users prepare a spreadsheet with a handful of sequences. Topiary then: (1) Infers the taxonomic scope for the ASR study and finds relevant sequences by BLAST;(2) Does taxonomically informed sequence quality control and redundancy reduction; (3) Constructs a multiple sequence alignment; (4) Generates a maximum-likelihood gene tree; (5) Reconciles the gene tree to the species tree; (6) Reconstructs ancestral amino acid sequences; and (7) Determines branch supports. The pipeline returns annotated evolutionary trees, spreadsheets with sequences, and graphical summaries of ancestor quality. This is achieved by integrating modern phylogenetics software (Muscle5, RAxML-NG, GeneRax, and PastML) with online databases (NCBI and the Open Tree of Life). In this paper, we introduce non-expert readers to the steps required for ASR, describe the specific design choices made in topiary, provide a detailed protocol for users, and then validate the pipeline using datasets from a broad collection of protein families.

Summary MNase-seq (micrococcal nuclease sequencing) is used to map nucleosome positions in eukaryotic genomes to study the relationship between chromatin structure and DNA-dependent processes. Current protocols require at least two days to isolate nucleosome-protected DNA fragments. We have developed a streamlined protocol for S. cerevisiae and other fungi which takes only three hours. Modified protocols were developed for wild fungi and mammalian cells. This method for rapidly producing sequencing-ready nucleosome footprints from several organisms makes MNase-seq faster and easier, with less chemical waste.

Eukaryotic DNA is packaged into nucleosomes, the smallest repeating unit of chromatin. The positions of nucleosomes determine the relative accessibility of genomic DNA. Several protocols exist for mapping nucleosome positions in eukaryotic genomes in order to study the relationship between chromatin structure and DNA-dependent processes. These nucleosome mapping protocols can be laborious and, at minimum, require two to three days to isolate nucleosome-protected DNA fragments. We have developed a streamlined protocol for mapping nucleosomes from S. cerevisiae liquid culture or from patches on solid agar. This method isolates nucleosome-sized footprints in three hours using 1.5 ml tubes with minimal chemical waste. We validate that these footprints match those produced by previously published methods and we demonstrate that our protocol works for N. crassa and S. pombe. A slightly modified protocol can be used for isolation of nucleosome-protected DNA fragments from a variety of wild fungal specimens thereby providing a simple, easily multiplexed and unified strategy to map nucleosome positions in model and non-model fungi. Finally, we demonstrate recovery of nucleosome footprints from the diploid myeloid leukemia cell line PLB-985 in less than three hours using an abbreviated version of the same protocol. With reduced volume and incubation times and a streamlined workflow, the described method should be compatible with high-throughput, automated creation of MNase-seq libraries. We believe this simple validated method for rapidly producing sequencing-ready nucleosome footprints from a variety of organisms will make nucleosome mapping studies widely accessible to researchers globally.

No abstract