Genotyping 1000 yeast strains by next-generation sequencing
- Stefan Wilkening†1,
- Manu M Tekkedil†1,
- Gen Lin†1,
- Emilie S Fritsch1,
- Wu Wei1,
- Julien Gagneur2,
- David W Lazinski3,
- Andrew Camilli3 and
- Lars M Steinmetz1Email author
© Wilkening et al.; licensee BioMed Central Ltd. 2013
Received: 7 December 2012
Accepted: 6 February 2013
Published: 9 February 2013
The throughput of next-generation sequencing machines has increased dramatically over the last few years; yet the cost and time for library preparation have not changed proportionally, thus representing the main bottleneck for sequencing large numbers of samples. Here we present an economical, high-throughput library preparation method for the Illumina platform, comprising a 96-well based method for DNA isolation for yeast cells, a low-cost DNA shearing alternative, and adapter ligation using heat inactivation of enzymes instead of bead cleanups.
Up to 384 whole-genome libraries can be prepared from yeast cells in one week using this method, for less than 15 euros per sample. We demonstrate the robustness of this protocol by sequencing over 1000 yeast genomes at ~30x coverage. The sequence information from 768 yeast segregants derived from two divergent S. cerevisiae strains was used to generate a meiotic recombination map at unprecedented resolution. Comparisons to other datasets indicate a high conservation of recombination at a chromosome-wide scale, but differences at the local scale. Additionally, we detected a high degree of aneuploidy (3.6%) by examining the sequencing coverage in these segregants. Differences in allele frequency allowed us to attribute instances of aneuploidy to gains of chromosomes during meiosis or mitosis, both of which showed a strong tendency to missegregate specific chromosomes.
Here we present a high throughput workflow to sequence genomes of large number of yeast strains at a low price. We have used this workflow to obtain recombination and aneuploidy data from hundreds of segregants, which can serve as a foundation for future studies of linkage, recombination, and chromosomal aberrations in yeast and higher eukaryotes.
KeywordsNext-generation sequencing High throughput DNA isolation Yeast DNA fragmentation Heat inactivation Recombination Aneuploidy
Results and discussion
It is critical to disrupt the cell wall of S. cerevisiae before extracting genomic DNA (gDNA) from the cells. One can use either physical disruption (strong vortexing with glass beads ), or enzymes like zymolyase and lyticase . Since strong vortexing in 96-well plates with phenol and glass beads can disrupt adhesiveness of the plate seals, this approach risks leakage and cross-contamination. Hence, we used enzymatic cell wall disruption for gDNA isolation. We isolated gDNA from up to 384 samples per day in 96-well format using a Biomek FX liquid handling robot. Combining cell pellets from 4 ml of overnight culture, the method yielded ~5.6 μg of DNA (average CV 2.6). This yield was slightly higher than from protocols that use glass beads [21, 22] and a commercial column based method (“DNeasy 96 Blood & Tissue Kit”, Qiagen) in our hands, and was highly cost-effective (0.8 €/sample). For all tested protocols, the DNA was of high quality as determined by gel imaging and absorbance ratios (260/280 and 260/230 ratios 1.8 - 2.2). Furthermore, the isolated DNA contained enough mitochondrial DNA to genotype the mitochondrial genome in most of the segregants.
Most standard library preparation methods perform purification after blunt-ending, A-addition, and ligation steps to avoid carryover of enzymes. Many recent high-throughput protocols [4, 6, 11] have replaced column-based purifications and gel size selection steps with magnetic SPRI bead cleanups . The reuse of the beads  and the use of a homemade bead mix  have also been applied to further reduce the cost of bead cleanups. Here we use heat inactivation, thereby circumventing purifications after blunt-ending, A-addition, and adapter ligation. This also reduces the risk of cross-contamination and sample loss during cleanups. The yield of the heat-inactivation protocol is comparable to the standard protocol (50- to 100-fold increase after PCR). <1% of read pairs have different barcodes on their forward and reverse sequences, indicating that the libraries obtained from this protocol have proper adapter ligation. The libraries are high-quality, with 87% mappability and 2.2% PCR duplicates, (detailed comparisons in Additional file 1: Table S2). In addition, the coverage of the S. cerevisiae genome yielded by our heat-inactivation protocol was highly uniform and comparable to libraries prepared with the standard Illumina protocol (Figure 3). A decrease in coverage was especially observed in regions with low GC content (<25%) when the heat-inactivation protocol was applied (examples are displayed in Additional file 2: Figure S1). This bias is slightly higher compared to the standard protocol (using SPRI cleanups), but was negligible for genotyping S. cerevisiae, as less than 0.5% of the 200 bp bins fall in this range. Genomes with a 30x coverage had 99% of the genome covered at > =1x and 97% at > =10x coverage. For the data shown here, we used 250 ng of fragmented DNA for the library preparation. We have also prepared libraries from starting amounts as low as 20 ng without a major loss in quality (see Additional file 1: Table S2 and Additional file 2: Figure S1). In principle, this would make the protocol compatible for RNA-Seq library preparations as well.
In this study, we used a set of 48 sequencing adapters containing 6 bp barcodes for ligation to the insert as reported by other groups [4, 6–9, 11]. After ligation, equimolar amounts of the barcoded libraries were pooled, size-selected, and amplified. The pooling of the samples before PCR resulted in moderately uneven barcode representation (Additional file 2: Figure S2), similar to previous reports [26, 27]; this, however, did not adversely affect our genotyping quality. Seven barcodes that displayed extremely poor performance in the pool were excluded in our subsequent studies (Additional file 1: Table S3). We did not observe any particular pattern among the poorly performing barcodes, except that three of them had an “AA” before the T-overhang. No significant barcode bias was observed when the samples were amplified individually and pooled at equimolar concentration before size selection and sequencing (data not shown). For sample sets with limited DNA amounts, we would therefore recommend performing the PCRs individually and to pool equimolar amounts of samples prior to size selection.
Recombination distributions between the S96xSK1 set and the S96xYJM789 set displayed a high correlation on a chromosome wide scale (0.944, P = 4.2e-08). To investigate possible differences on a local scale, we identified SNPs that are common between YJM789 and SK1 and then partitioned the S288c genome into non-overlapping bins (min 2 kb, max 3 kb) based on these SNPs. For this window size a lower correlation (0.616, P < 2.2e-16) was observed. A list of 20 regions with the largest differences in normalized recombination rates is provided in Additional file 1: Table S5.
Using the same partitioned bins as described above, we also compared the recombination rate with the genomic double-strand break (DSB) map generated by Pan et al.  (using immunoprecipitation of Spo11-bound oligos in meiotic SK1 cells). Similar to the comparison of our dataset with the S96xYJM789 dataset, we observed a good correlation on the chromosome-scale (0.726, P = 1.44e-03), but a lower one on the finer scale (0.375, P < 2.2e-16). These differences in hotspot intensities could be due to S96-specific hotspots or the possibility that not all DSBs lead to a detectable recombination event. Plotting the distance from the center of Spo11 oligo hotspots to the center of the S96xSK1 recombination events revealed a significant drop in recombination frequency in the vicinity (400-500 nucleotides) of the Spo11 hotspot (Additional file 2: Figure S3). This drop could be explained by the 5′ to 3′ resection of the resulting DNA ends, required for the repair of DSB by homologous recombination .
In this study, we present various optimization steps for whole DNA-Seq library preparation, considerably reducing the time and cost for library preparation compared to standard procedures. These include efficient high-throughput DNA isolation from yeast cells, a cost-effective alternative to standard Covaris fragmentation, and a library preparation that avoids most cleanup steps. The protocol was developed for the Illumina platform, but most of the steps are adaptable to other sequencing platforms with minor modifications. The quality of the DNA and final library was similar to that obtained by standard techniques. Although our heat inactivation step resulted in a slightly reduced coverage of regions with extreme GC content, this did not interfere with genotype calling. The genotype data was also used to map quantitative traits (Wilkening et al., in revision), for which sample size and marker resolution are critical to maximize mapping resolution and statistical power. Furthermore, we created a map of meiotic recombination points in yeast with a yet unprecedented resolution as well as a catalog of chromosomal aberrations. Despite a high conservation of recombination at a chromosome-wide scale, our results indicate differences at the local scale. We also found an unexpectedly high degree of chromosomal aberrations in this genetic background. In conclusion, our method is a rapid, high-throughput approach for genotyping many small genomes or target-enriched DNA, and our results provide a unique basis for future and current studies of aneuploidy and recombination.
DNA isolation from yeast cells
A modified version of the PrepEase Genomic DNA Isolation Kit (Affymetrix, 78855 1 KT) based on enzymatic cell wall digestion was used for the DNA isolation from yeast cells. This protocol can easily be applied to blood, bacteria or homogenized tissue or plant material by substituting buffers according the manufacturer's instructions. All of the mixing steps were performed by pipetting using a Biomek FX pipetting robot (Beckman Coulter) in 96-well plates. Cell pellets from 4 deep-well plates, each containing 1 ml overnight culture, were combined for the DNA isolation. The “Spheroplast” and “Enzyme Solution” from the kit was replaced by Qiagen’s Y1 lysis buffer (1 M sorbitol, 100 mM EDTA, pH 8.0, 14 mM β-mercaptoethanol) freshly supplemented with 2.5 μl/ml of Zymolyase (Seikagaku Inc.) and 2.5 μl/ml RNase A (10 mg/ml, Qiagen). 200 μl of this buffer was added to each pellet, mixed and incubated at 37°C for 90 min with gentle shaking every 30 min. 200 μl of water was added to each well and the plate was centrifuged at 6000 x g for 4 min (for centrifuges with maximum 3,000 x g, centrifugation times can be tripled) and the supernatant was decanted. 120 μl Homogenization Buffer was added and mixed to resuspend the pellet completely. 100 μl of chloroform and 400 μl of Protein Precipitation Buffer were added to the lysate and mixed. Plates were centrifuged at 6,000 × g for 15 min. 450 μl of the upper aqueous phase was transferred with the robot (pipetting height was optimized in advance) to a 1 ml deep-well plate containing 340 μl of isopropanol per well. The solution was mixed, left for 15 min at room temperature, and centrifuged at 6,000 × g for 15 min. After decanting the supernatant, 1 ml of cold 70% ethanol was added to the pellet, mixed and centrifuged at 6,000 × g for 10 min. The supernatant was decanted, and the tube was placed upside down on a paper towel and dried for 5 min at 37°C. The DNA pellet was resuspended in 300 μl of DNA Resuspension Buffer or Elution Buffer (EB, 10 mM Tris HCl) by shaking plates for 30 min at 37°C and later by mixing. A detailed Biomek protocol including .bmf files is provided in Additional file 3.
Additional file 4: Video S1: Showing fragmentation by Bandelin sonication. (MOV 3 MB)
End repair, dA-tailing, and ligation using heat-inactivation
Instead of the standard column or bead-based cleanup steps, we heat-inactivated the enzymes used for end repair, dA-tailing, and ligation, then added the respective enzyme (+ buffer). For this, 250 ng of fragmented gDNA were used for the library preparation in 96-well PCR plates in a volume of 17 μl. End repair for the fragments was performed by adding 3 μl of End repair master mix composed of 2 μl of End repair buffer and 1 μl End repair enzyme (NEBNext End Repair Module, NEB #E6050L) using a 8-channel pipette. The contents were mixed by vortexing, shortly spun down, and incubated in a thermocycler at 20°C for 45 min. The enzymes were then heat-inactivated at 75°C for 15 min. The contents of the plate were quickly spun down, and 2 μl of A tailing master mix containing 1 μl of Klenow Fragment (3′→5′ exo–) (NEB #M0212L), 0.5 μl of nuclease free water, and 0.5 μl of 100 mM dATP (NEB #N0440S) were added to the 20 μl reaction. The contents were mixed by vortexing, spun down, and incubated in a thermocycler at 37°C for 45 minutes. The enzymes were then heat-inactivated at 75°C for 15 minutes. 5 μl of ligation master mix containing 3 μl 10X T4 DNA ligase buffer and 2 μl T4 DNA ligase (NEB #M0202L) were added to the reaction followed by 3 μl of 7 μM multiplex barcode adapters (aliquoted into 8-strip PCR tubes or 96-well plates for convenient pipetting, see Additional file 1: Table S2 for sequences). The concentration of adapters was optimized to reduce the formation of adapter dimers. The reaction contents were mixed well, spun down, and incubated on a thermocycler at 16°C for 1 h followed by heat inactivation at 75°C for 15 min.
Pooling and size selection
After barcode ligation and heat inactivation, 48 samples were pooled together by combining 5 μl of each sample in a 1.5 ml reaction tube. The samples were cleaned up and concentrated to 40 μl using 1x Ampure XP cleanup. This pooling step reduces the sample size from 96 to two for the subsequent size selection and PCR. For size selection, 25 μl (roughly 1.25 μg) of ligated DNA was loaded on a 2% E-Gel SizeSelect (Invitrogen) and DNA fragments were collected at 350 bp and 400 bp. The DNA concentrations were then determined by Qubit HS DNA reagent (Invitrogen).
In a 50 μl reaction, 5–10 ng of the pooled libraries were amplified. We have observed that performing PCR with an excess of template DNA (>20 ng) significantly reduces the efficiency of the PCR. The PCR was performed on a thermocycler (MJ Research tetrad) containing 1x Phusion Master Mix with HF Buffer (Thermo Scientific) and 0.2 μM Illumina PE 1.0 and 2.0 primers. The low primer concentration reduced the formation of primer dimers often observed at standard primer concentrations (1.25 μM). PCR conditions were 98°C for 45 s, 10x [98°C for 15 s, 65°C for 30 s, 72°C for 30 s], 72°C for 5 min, 4°C hold.
DNA purification with SPRI beads
The amplified DNA was purified in 0.2 ml PCR strips by mixing the DNA with 1x volume of Agencourt AMPure XP beads (Beckman Coulter) and select the magnetic beads with a homemade magnetic stand (Additional file 2: Figure S6). This stand consists of neodymium magnets (Webcraft GmbH, Gottmadingen, Germany) mounted on trimmed 96-well plates and can be used in combination with an 8-channel pipet. For a high-throughput SPRI clean-up we further provide a detailed protocol of the pipetting steps for 96-well plates and Biomek robot in Additional file 3. DNA concentrations were quantified for the subsequent library preparation (see DNA quantification and quality control).
DNA quantification and quality control
The quality of individual samples of isolated DNA was determined by a photospectrometric measurement using a NanoDrop 1000 (Thermo Fisher). For quantification of genomic DNA and pre-PCR libraries in 96-well plates, we used Quant-iT PicoGreen dsDNA Reagent (Invitrogen) in optical plates (Greiner). The fluorescence was measured at 485 nm excitation and 535 nm emission in a Genios microplate reader (Tecan) according to the manufacturer’s instructions. The pooled libraries were quantified before and after PCR with a Qubit spectrofluorometer (Invitrogen) according to the manufacturer’s instructions. All pre- and post-PCR libraries were run on a High Sensitivity Bioanalyzer chip (Agilent) to determine the size distribution. After a 10-cycle PCR, we typically observed a 10-fold increase in DNA amount and a 24–30 bp increase in library size due to adapter elongation. Depending on the amplification efficiency, either the low (350 bp) or the high molecular weight (400 bp) library was selected for sequencing. The samples were then diluted to 10 nM and clustered on the Illumina cBot clustering station for paired-end sequencing on an Illumina HiSeq 2000.
Design of 48 multiplexing barcodes
We designed a set of 48 adapters, each with a different hexamer sequence just before the T-overhang, similar to Lefrancois et al . We selected 64 of 96 Illumina barcodes, which had at least a 3 bp difference compared to any of the other 63 barcodes. From this set, 48 barcodes that had an equilibrated base composition at the first two bases (for better cluster calling) were manually chosen. Following quality control analysis, we replaced the seven poorest performing barcodes with new ones (Additional file 1: Table S2).
To demultiplex, we extracted the first six bases of each read and compared it to all possible barcodes. The perfect match or best hit to one barcode with the least number of mismatches was assigned to the read. For genotyping, reads from the segregants along with both SK1 and S96 (a haploid strain isogenic to S288c) parental strains were aligned to the S288c reference genome (build R63) using Novoalign (v2.07.06; http://www.novocraft.com/), allowing for unique alignments. Thereafter GATK was used for realignment and recalibration of the bam files , and subsequent SNP calling was performed using SAMtools . The vcf file produced by SAMtools contains a list of variant positions and the individual genotype calls across all samples at each variant position. The formula that SAMtools applies for calling the genotype is dependent on allele frequency, which is not directly applicable to our study, because the allele frequency at true SNP positions is expected to be 0.5 in crosses generated from 2 parents. Instead, we used the genotype likelihood (PL stats generated by GATK) to infer the genotype. SNP positions, which correspond to a homozygous reference call in the S96 parent and a homozygous variant in the SK1 parent, are chosen first. From this set of SNPs, we excluded calls whose allele frequency is not between 0.3-0.7. These SNP calls are unreliable and often not in linkage with their surrounding SNPs, and could either be SNPs within regions that are repetitive in one but not in the other parent, or result from misaligned reads.
GC bias and coverage plots
We calculated the genome-wide, per-base coverage of the S288c genome using SAMtools. Positions where all samples had at least 1 read were considered. The density was plotted with a bandwidth of 0.1. For plotting GC bias, the genome was divided into non-overlapping 200 bp bins, and the depth was estimated by the mean values of per-base depth in these bins. Bins with less than 50% covered by at least 1 read were excluded. All analyses were run in the software R (v. 2.12.0; http://cran.r-project.org). For analyzing chromosomal abnormalities, an identical method for binning and GC correction was applied, except that a 10 kb bin size was used. GC bias correction was applied using a LOESS method, as described previously .
Recombination map analysis
For both genotype datasets (S96xSK1 and S96xYJM789) the rqtl package (with the function, est.map (maxit = 1000,error.prob = 0.01) was used to construct the genetic map for both crosses. After obtaining the genetic map, the genotypes were filtered for errors and crossovers counted for each segregant, using functions in rqtl (cleanGeno(maxdist = 2.5, maxmark = 2) followed by countXO). For 2-3 kb bins (partitioned by common SNPs), the recombination rate was calculated as genetic distance between 2 SNPs/physical distance between 2 SNPs. For identifying regions with difference in recombination rates, we normalized the rate in both, S96xSK1 and S96xYJM789 by setting the mean of each set to 1. Raw sequences for Spo11 oligo maps were download from SRA (GSE26452) and aligned to the S288c genome build R63 using bowtie2 allowing for only unique alignments.
NGS of genomic DNA
Solid phase reversible immobilization
Adaptive Focused Acoustics
Genome Analysis Toolkit
Variant call format
Double strand break
Coefficient of variation.
We thank Adam Deutschbauer (Lawrence Berkeley National Laboratory), Michelle Nguyen, and Raquel Kuehn (Stanford Genome Technology Center) for the construction of the 768 S96xSK1 segregants. We also thank Raeka Aiyar (EMBL) for assistance in writing the manuscript and Eugenio Mancera (UCSF) and Vicent Pelechano (EMBL) for fruitful discussions. This study was technically supported by the EMBL Genomics Core facility, where the libraries were sequenced. This work was supported by grants from the National Institutes of Health and the Deutsche Forschungsgemeinschaft to LMS.
- Lander ES: Initial impact of the sequencing of the human genome. Nature. 2011, 470 (7333): 187-197. 10.1038/nature09792.View ArticlePubMedGoogle Scholar
- Quail MA, Kozarewa I, Smith F, Scally A, Stephens PJ, Durbin R, Swerdlow H, Turner DJ: A large genome center's improvements to the Illumina sequencing system. Nat Methods. 2008, 5 (12): 1005-1010. 10.1038/nmeth.1270.PubMed CentralView ArticlePubMedGoogle Scholar
- Adey A, Morrison HG, Asan , Xun X, Kitzman JO, Turner EH, Stackhouse B, MacKenzie AP, Caruccio NC, Zhang X, et al: Rapid, low-input, low-bias construction of shotgun fragment libraries by high-density in vitro transposition. Genome Biol. 2010, 11 (12): R119-10.1186/gb-2010-11-12-r119.PubMed CentralView ArticlePubMedGoogle Scholar
- Andolfatto P, Davison D, Erezyilmaz D, Hu TT, Mast J, Sunayama-Morita T, Stern DL: Multiplexed shotgun genotyping for rapid and efficient genetic mapping. Genome Res. 2011, 21 (4): 610-617. 10.1101/gr.115402.110.PubMed CentralView ArticlePubMedGoogle Scholar
- Fisher S, Barry A, Abreu J, Minie B, Nolan J, Delorey TM, Young G, Fennell TJ, Allen A, Ambrogio L, et al: A scalable, fully automated process for construction of sequence-ready human exome targeted capture libraries. Genome Biol. 2011, 12 (1): R1-10.1186/gb-2011-12-1-r1.PubMed CentralView ArticlePubMedGoogle Scholar
- Lennon NJ, Lintner RE, Anderson S, Alvarez P, Barry A, Brockman W, Daza R, Erlich RL, Giannoukos G, Green L, et al: A scalable, fully automated process for construction of sequence-ready barcoded libraries for 454. Genome Biol. 2010, 11 (2): R15-10.1186/gb-2010-11-2-r15.PubMed CentralView ArticlePubMedGoogle Scholar
- Rohland N, Reich D: Cost-effective, high-throughput DNA sequencing libraries for multiplexed target capture. Genome Res. 2012, 22 (5): 939-946. 10.1101/gr.128124.111.PubMed CentralView ArticlePubMedGoogle Scholar
- Farias-Hesson E, Erikson J, Atkins A, Shen P, Davis RW, Scharfe C, Pourmand N: Semi-automated library preparation for high-throughput DNA sequencing platforms. J Biomed Biotechnol. 2010, 2010: 617469-PubMed CentralView ArticlePubMedGoogle Scholar
- Lefrancois P, Euskirchen GM, Auerbach RK, Rozowsky J, Gibson T, Yellman CM, Gerstein M, Snyder M: Efficient yeast ChIP-Seq using multiplex short-read DNA sequencing. BMC Genomics. 2009, 10: 37-10.1186/1471-2164-10-37.PubMed CentralView ArticlePubMedGoogle Scholar
- Meyer M, Kircher M: Illumina sequencing library preparation for highly multiplexed target capture and sequencing. Cold Spring Harb Protoc. 2010, 2010 (6): pdb prot5448-View ArticlePubMedGoogle Scholar
- Borgstrom E, Lundin S, Lundeberg J: Large scale library generation for high throughput sequencing. PLoS One. 2011, 6 (4): e19119-10.1371/journal.pone.0019119.PubMed CentralView ArticlePubMedGoogle Scholar
- Xu Z, Wei W, Gagneur J, Perocchi F, Clauder-Munster S, Camblong J, Guffanti E, Stutz F, Huber W, Steinmetz LM: Bidirectional promoters generate pervasive transcription in yeast. Nature. 2009, 457 (7232): 1033-1037. 10.1038/nature07728.PubMed CentralView ArticlePubMedGoogle Scholar
- Nagalakshmi U, Wang Z, Waern K, Shou C, Raha D, Gerstein M, Snyder M: The transcriptional landscape of the yeast genome defined by RNA sequencing. Science. 2008, 320 (5881): 1344-1349. 10.1126/science.1158441.PubMed CentralView ArticlePubMedGoogle Scholar
- Swinnen S, Thevelein JM, Nevoigt E: Genetic mapping of quantitative phenotypic traits in Saccharomyces cerevisiae. FEMS Yeast Res. 2011, 12 (2): 215-227.View ArticleGoogle Scholar
- Steinmetz LM, Davis RW: Maximizing the potential of functional genomics. Nat Rev Genet. 2004, 5 (3): 190-201. 10.1038/nrg1293.View ArticlePubMedGoogle Scholar
- Liti G, Carter DM, Moses AM, Warringer J, Parts L, James SA, Davey RP, Roberts IN, Burt A, Koufopanou V, et al: Population genomics of domestic and wild yeasts. Nature. 2009, 458 (7236): 337-341. 10.1038/nature07743.PubMed CentralView ArticlePubMedGoogle Scholar
- Mancera E, Bourgon R, Brozzi A, Huber W, Steinmetz LM: High-resolution mapping of meiotic crossovers and non-crossovers in yeast. Nature. 2008, 454 (7203): 479-485. 10.1038/nature07135.PubMed CentralView ArticlePubMedGoogle Scholar
- Pan J, Sasaki M, Kniewel R, Murakami H, Blitzblau HG, Tischfield SE, Zhu X, Neale MJ, Jasin M, Socci ND, et al: A hierarchical combination of factors shapes the genome-wide topography of yeast meiotic recombination initiation. Cell. 2011, 144 (5): 719-731. 10.1016/j.cell.2011.02.009.PubMed CentralView ArticlePubMedGoogle Scholar
- Hoffman CS, Winston F: A ten-minute DNA preparation from yeast efficiently releases autonomous plasmids for transformation of Escherichia coli. Gene. 1987, 57 (2–3): 267-272.View ArticlePubMedGoogle Scholar
- van Burik JA, Schreckhise RW, White TC, Bowden RA, Myerson D: Comparison of six extraction techniques for isolation of DNA from filamentous fungi. Med Mycol. 1998, 36 (5): 299-303.View ArticlePubMedGoogle Scholar
- Harju S, Fedosyuk H, Peterson KR: Rapid isolation of yeast genomic DNA: Bust n' Grab. BMC Biotechnol. 2004, 4: 8-10.1186/1472-6750-4-8.PubMed CentralView ArticlePubMedGoogle Scholar
- Loeffler J, Schmidt K, Hebart H, Schumacher U, Einsele H: Automated extraction of genomic DNA from medically important yeast species and filamentous fungi by using the MagNA pure LC system. J Clin Microbiol. 2002, 40 (6): 2240-2243. 10.1128/JCM.40.6.2240-2243.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Syed F, Grunenwald H, Caruccio N: Next-generation sequencing library preparation: simultaneous fragmentation and tagging using in vitro transposition. Nat Methods. 2009, 6: i-ii.Google Scholar
- Sexton T, Yaffe E, Kenigsberg E, Bantignies F, Leblanc B, Hoichman M, Parrinello H, Tanay A, Cavalli G: Three-dimensional folding and functional organization principles of the Drosophila genome. Cell. 2012, 148 (3): 458-472. 10.1016/j.cell.2012.01.010.View ArticlePubMedGoogle Scholar
- DeAngelis MM, Wang DG, Hawkins TL: Solid-phase reversible immobilization for the isolation of PCR products. Nucleic Acids Res. 1995, 23 (22): 4742-4743. 10.1093/nar/23.22.4742.PubMed CentralView ArticlePubMedGoogle Scholar
- Alon S, Vigneault F, Eminaga S, Christodoulou DC, Seidman JG, Church GM, Eisenberg E: Barcoding bias in high-throughput multiplex sequencing of miRNA. Genome Res. 2011, 21 (9): 1506-1511. 10.1101/gr.121715.111.PubMed CentralView ArticlePubMedGoogle Scholar
- Van Nieuwerburgh F, Soetaert S, Podshivalova K, Ay-Lin Wang E, Schaffer L, Deforce D, Salomon DR, Head SR, Ordoukhanian P: Quantitative bias in Illumina TruSeq and a novel post amplification barcoding strategy for multiplexed DNA and small RNA deep sequencing. PLoS One. 2011, 6 (10): e26969-10.1371/journal.pone.0026969.PubMed CentralView ArticlePubMedGoogle Scholar
- Browning SR, Browning BL: Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering. Am J Hum Genet. 2007, 81 (5): 1084-1097. 10.1086/521987.PubMed CentralView ArticlePubMedGoogle Scholar
- Martini E, Borde V, Legendre M, Audic S, Regnault B, Soubigou G, Dujon B, Llorente B: Genome-wide analysis of heteroduplex DNA in mismatch repair-deficient yeast cells reveals novel properties of meiotic recombination pathways. PLoS Genet. 2011, 7 (9): e1002305-10.1371/journal.pgen.1002305.PubMed CentralView ArticlePubMedGoogle Scholar
- Paull TT: Making the best of the loose ends: Mre11/Rad50 complexes and Sae2 promote DNA double-strand break resection. DNA Repair (Amst). 2010, 9 (12): 1283-1291. 10.1016/j.dnarep.2010.09.015.View ArticleGoogle Scholar
- Schacherer J, Shapiro JA, Ruderfer DM, Kruglyak L: Comprehensive polymorphism survey elucidates population structure of Saccharomyces cerevisiae. Nature. 2009, 458 (7236): 342-345. 10.1038/nature07670.PubMed CentralView ArticlePubMedGoogle Scholar
- Hunter N, Chambers SR, Louis EJ, Borts RH: The mismatch repair system contributes to meiotic sterility in an interspecific yeast hybrid. EMBO J. 1996, 15 (7): 1726-1733.PubMed CentralPubMedGoogle Scholar
- D'Amours D, Stegmeier F, Amon A: Cdc14 and condensin control the dissolution of cohesin-independent chromosome linkages at repeated DNA. Cell. 2004, 117 (4): 455-469. 10.1016/S0092-8674(04)00413-1.View ArticlePubMedGoogle Scholar
- Sullivan M, Higuchi T, Katis VL, Uhlmann F: Cdc14 phosphatase induces rDNA condensation and resolves cohesin-independent cohesion during budding yeast anaphase. Cell. 2004, 117 (4): 471-482. 10.1016/S0092-8674(04)00415-5.View ArticlePubMedGoogle Scholar
- Torres EM, Sokolsky T, Tucker CM, Chan LY, Boselli M, Dunham MJ, Amon A: Effects of aneuploidy on cellular physiology and cell division in haploid yeast. Science. 2007, 317 (5840): 916-924. 10.1126/science.1142210.View ArticlePubMedGoogle Scholar
- Anders KR, Kudrna JR, Keller KE, Kinghorn B, Miller EM, Pauw D, Peck AT, Shellooe CE, Strong IJ: A strategy for constructing aneuploid yeast strains by transient nondisjunction of a target chromosome. BMC Genet. 2009, 10: 36-PubMed CentralView ArticlePubMedGoogle Scholar
- McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, et al: The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010, 20 (9): 1297-1303. 10.1101/gr.107524.110.PubMed CentralView ArticlePubMedGoogle Scholar
- Li H: A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics. 2011, 27 (21): 2987-2993. 10.1093/bioinformatics/btr509.PubMed CentralView ArticlePubMedGoogle Scholar
- Alkan C, Kidd JM, Marques-Bonet T, Aksay G, Antonacci F, Hormozdiari F, Kitzman JO, Baker C, Malig M, Mutlu O, et al: Personalized copy number and segmental duplication maps using next-generation sequencing. Nat Genet. 2009, 41 (10): 1061-1067. 10.1038/ng.437.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.