Pooled Sequencing Research Articles

Pervasive Linked Selection and Intermediate-Frequency Alleles Are Implicated in an Evolve-and-Resequencing Experiment of Drosophila simulans.

We develop analytical and simulation tools for evolve-and-resequencing experiments and apply them to a new study of rapid evolution in Drosophila simulans Likelihood test statistics applied to pooled population sequencing data suggest parallel evolution of 138 SNPs across the genome. This number is reduced by orders of magnitude from previous studies (thousands or tens of thousands), owing to differences in both experimental design and statistical analysis. Whole genome simulations calibrated from Drosophila genetic data sets indicate that major features of the genome-wide response could be explained by as few as 30 loci under strong directional selection with a corresponding hitchhiking effect. Smaller effect loci are likely also responding, but are below the detection limit of the experiment. Finally, SNPs showing strong parallel evolution in the experiment are intermediate in frequency in the natural population (usually 30-70%) indicative of balancing selection in nature. These loci also exhibit elevated differentiation among natural populations of D. simulans, suggesting environmental heterogeneity as a potential balancing mechanism.

Extreme copy number variation at a tRNA ligase gene affecting phenology and fitness in yellow monkeyflowers.

Copy number variation (CNV) is a major part of the genetic diversity segregating within populations, but remains poorly understood relative to single nucleotide variation. Here, we report on a tRNA ligase gene (Migut.N02091; RLG1a) exhibiting unprecedented, and fitness-relevant, CNV within an annual population of the yellow monkeyflower Mimulus guttatus. RLG1a variation was associated with multiple traits in pooled population sequencing (PoolSeq) scans of phenotypic and phenological cohorts. Resequencing of inbred lines revealed intermediate-frequency three-copy variants of RLG1a (trip+; 5/35=14%), and trip+ lines exhibited elevated RLG1a expression under multiple conditions. trip+ carriers, in addition to being over-represented in late-flowering and large-flowered PoolSeq populations, flowered later under stressful conditions in a greenhouse experiment (p<0.05). In wild population samples, we discovered an additional rare RLG1a variant (high+) that carries 250-300 copies of RLG1a totalling ~5.7Mb (20-40% of a chromosome). In the progeny of a high+ carrier, Mendelian segregation of diagnostic alleles and qPCR-based copy counts indicate that high+ is a single tandem array unlinked to the single-copy RLG1a locus. In the wild, high+ carriers had highest fitness in two particularly dry and/or hot years (2015 and 2017; both p<0.01), while single-copy individuals were twice as fecund as either CNV type in a lush year (2016: p<0.005). Our results demonstrate fluctuating selection on CNVs affecting phenological traits in a wild population, suggest that plant tRNA ligases mediate stress-responsive life-history traits, and introduce a novel system for investigating the molecular mechanisms of gene amplification.

Evaluation of genome scaffolding tools using pooled clone sequencing.

DNA sequencing technologies hold great promise in generating information that will guide scientists to understand how the genome effects human health and organismal evolution. The process of generating raw genome sequence data becomes cheaper and faster, but more error-prone. Assembly of such data into high-quality finished genome sequences remains challenging. Many genome assembly tools are available, but they differ in terms of their performance and their final output. More importantly, it remains largely unclear how to best assess the quality of assembled genome sequences. Here we evaluate the accuracies of several genome scaffolding algorithms using two different types of data generated from the genome of the same human individual: whole genome shotgun sequencing (WGS) and pooled clone sequencing (PCS). We observe that it is possible to obtain better assemblies if PCS data are used, compared to using only WGS data. However, the current scaffolding algorithms are developed only for WGS, and PCS-aware scaffolding algorithms remain an open problem.

Small and genetically highly structured populations in a long-legged bee, Rediviva longimanus, as inferred by pooled RAD-seq

Adaptation to local host plants may impact a pollinator’s population genetic structure by reducing gene flow and driving population genetic differentiation, representing an early stage of ecological speciation. South African Rediviva longimanus bees exhibit elongated forelegs, a bizarre adaptation for collecting oil from floral spurs of their Diascia hosts. Furthermore, R. longimanus foreleg length (FLL) differs significantly among populations, which has been hypothesised to result from selection imposed by inter-population variation in Diascia floral spur length. Here, we used a pooled restriction site-associated DNA sequencing (pooled RAD-seq) approach to investigate the population genetic structure of R. longimanus and to test if phenotypic differences in FLL translate into increased genetic differentiation (i) between R. longimanus populations and (ii) between phenotypes across populations. We also inferred the effects of demographic processes on population genetic structure and tested for genetic markers underpinning local adaptation.ResultsPopulations showed marked genetic differentiation (average FST = 0.165), though differentiation was not statistically associated with differences between populations in FLL. All populations exhibited very low genetic diversity and were inferred to have gone through recent bottleneck events, suggesting extremely low effective population sizes. Genetic differentiation between samples pooled by leg length (short versus long) rather than by population of origin was even higher (FST = 0.260) than between populations, suggesting reduced interbreeding between long and short-legged individuals. Signatures of selection were detected in 1119 (3.8%) of a total of 29,721 SNP markers,ConclusionsPopulations of R. longimanus appear to be small, bottlenecked and isolated. Though we could not detect the effect of local adaptation (FLL in response to floral spurs of host plants) on population genetic differentiation, short and long legged bees appeared to be partially differentiated, suggesting incipient ecological speciation. To test this hypothesis, greater resolution through the use of individual-based whole-genome analyses is now needed to quantify the degree of reproductive isolation between long and short legged bees between and even within populations.

Pooled DNA sequencing to identify SNPs associated with a major QTL for bacterial wilt resistance in Italian ryegrass (Lolium multiflorum Lam.)

Key messageSNPs and candidate genes associated with bacterial wilt resistance in Italian ryegrass were identified by sequencing the parental plants and pooled F1 progeny of a segregating population.Italian ryegrass (Lolium multiflorum Lam.) is one of the most important forage grass species in temperate regions. Its yield, quality and persistency can significantly be reduced by bacterial wilt, a serious disease caused by Xanthomonas translucens pv. graminis. Although a major QTL for bacterial wilt resistance has previously been reported, detailed knowledge on underlying genes and DNA markers to allow for efficient resistance breeding strategies is currently not available. We used pooled DNA sequencing to characterize a major QTL for bacterial wilt resistance of Italian ryegrass and to develop inexpensive sequence-based markers to efficiently target resistance alleles for marker-assisted recurrent selection. From the mapping population segregating for the QTL, DNA of 44 of the most resistant and 44 of the most susceptible F1 individuals was pooled and sequenced using the Illumina HiSeq 2000 platform. Allele frequencies of 18 × 106 single nucleotide polymorphisms (SNP) were determined in the resistant and susceptible pool. A total of 271 SNPs on 140 scaffold sequences of the reference parental genome showed significantly different allele frequencies in both pools. We converted 44 selected SNPs to KASP™ markers, genetically mapped these proximal to the major QTL and thus validated their association with bacterial wilt resistance. This study highlights the power of pooled DNA sequencing to efficiently target binary traits in biparental mapping populations. It delivers genome sequence data, SNP markers and potential candidate genes which will allow to implement marker-assisted strategies to fix bacterial wilt resistance in outcrossing breeding populations of Italian ryegrass.

Evolution of drug resistance in an antifungal-naive chronic Candida lusitaniae infection

Management of the limited number of antimicrobials currently available requires the identification of infections that contain drug-resistant isolates and the discovery of factors that promote the evolution of drug resistance. Here, we report a single fungal infection in which we have identified numerous subpopulations that differ in their alleles of a single gene that impacts drug resistance. The diversity at this locus was markedly greater than the reported heterogeneity of alleles conferring antibiotic resistance in bacterial infections. Analysis of genomes from hundreds of Clavispora (Candida) lusitaniae isolates, through individual and pooled isolate sequencing, from a single individual with cystic fibrosis revealed at least 25 nonsynonymous mutations in MRR1, which encodes a transcription factor capable of inducing fluconazole (FLZ) resistance in Candida species. Isolates with high-activity Mrr1 variants were resistant to FLZ due to elevated expression of the MDR1-encoded efflux pump. We found that high Mrr1-regulated Mdr1 activity protected against host and bacterial factors, suggesting drug resistance can be selected for indirectly and perhaps explaining the Mrr1 heterogeneity in this individual who had no prior azole exposure. Regional analysis of C. lusitaniae populations from the upper and lower lobes of the right lung suggested intermingling of subpopulations throughout. Our retrospective characterization of sputum and lung populations by pooled sequencing found that alleles that confer FLZ resistance were a minority in each pool, possibly explaining why they were undetected before unsuccessful FLZ therapy. New susceptibility testing regimes may detect problematical drug-resistant subpopulations in heterogeneous single-species infections.

A new approach based on targeted pooled DNA sequencing identifies novel mutations in patients with Inherited Retinal Dystrophies

Inherited retinal diseases (IRD) are a heterogeneous group of diseases that mainly affect the retina; more than 250 genes have been linked to the disease and more than 20 different clinical phenotypes have been described. This heterogeneity both at the clinical and genetic levels complicates the identification of causative mutations. Therefore, a detailed genetic characterization is important for genetic counselling and decisions regarding treatment. In this study, we developed a method consisting on pooled targeted next generation sequencing (NGS) that we applied to 316 eye disease related genes, followed by High Resolution Melting and copy number variation analysis. DNA from 115 unrelated test samples was pooled and samples with known mutations were used as positive controls to assess the sensitivity of our approach. Causal mutations for IRDs were found in 36 patients achieving a detection rate of 31.3%. Overall, 49 likely causative mutations were identified in characterized patients, 14 of which were first described in this study (28.6%). Our study shows that this new approach is a cost-effective tool for detection of causative mutations in patients with inherited retinopathies.

Effects of hybridization on sea turtle fitness

Sea turtle hybridization is a common phenomenon in Brazil between loggerheads (Caretta caretta) and hawksbills (Eretmochelys imbricata) as well as between loggerheads and olive ridleys (Lepidochelys olivacea). In a previous study we showed that the reproductive output of loggerhead/hawksbill hybrids is similar to that of parental species, suggesting no negative effect of hybridization at this life stage. In this study, we used pooled amplicon sequencing to assign species identity to dams and their progeny, and to investigate the fitness consequences of hybridization, using hatchling viability as a proxy for fitness. We genotyped 4829 hatchlings from egg clutches laid by 78 loggerheads, 13 hawksbills, seven loggerhead/hawksbill hybrids, and three loggerhead/olive ridley hybrids. The proportion of viable hybrid (heterozygous) hatchlings was similar to that of homozygous hatchlings (based on data at two loci), independent of the dam’s genotype. Multiple species paternity was observed in 35.7% of the nests. Both hybrid males and females were fertile and produced viable offspring, and we found no evidence for hybrid breakdown. We suggest a genome-wide study of the hybrids and parental species to better characterize hybrids, as well as studies on additional demographic and ecological parameters to further assess the effects of hybridization and its consequences for sea turtles and their environment.

Avian Influenza Virus PB1 Gene in H3N2 Viruses Evolved in Humans To Reduce Interferon Inhibition by Skewing Codon Usage toward Interferon-Altered tRNA Pools.

ABSTRACTInfluenza A viruses cause an annual contagious respiratory disease in humans and are responsible for periodic high-mortality human pandemics. Pandemic influenza A viruses usually result from the reassortment of gene segments between human and avian influenza viruses. These avian influenza virus gene segments need to adapt to humans. Here we focus on the human adaptation of the synonymous codons of the avian influenza virus PB1 gene of the 1968 H3N2 pandemic virus. We generated recombinant H3N2 viruses differing only in codon usage of PB1 mRNA and demonstrated that codon usage of the PB1 mRNA of recent H3N2 virus isolates enhances replication in interferon (IFN)-treated human cells without affecting replication in untreated cells, thereby partially alleviating the interferon-induced antiviral state. High-throughput sequencing of tRNA pools explains the reduced inhibition of replication by interferon: the levels of some tRNAs differ between interferon-treated and untreated human cells, and evolution of the codon usage of H3N2 PB1 mRNA is skewed toward interferon-altered human tRNA pools. Consequently, the avian influenza virus-derived PB1 mRNAs of modern H3N2 viruses have acquired codon usages that better reflect tRNA availabilities in IFN-treated cells. Our results indicate that the change in tRNA availabilities resulting from interferon treatment is a previously unknown aspect of the antiviral action of interferon, which has been partially overcome by human-adapted H3N2 viruses.

Transmission ratio distortion is frequent in Arabidopsis thaliana controlled crosses

The equal probability of transmission of alleles from either parent during sexual reproduction is a central tenet of genetics and evolutionary biology. Yet, there are many cases where this rule is violated. The preferential transmission of alleles or genotypes is termed transmission ratio distortion (TRD). Examples of TRD have been identified in many species, implying that they are universal, but the resolution of species-wide studies of TRD are limited. We have performed a species-wide screen for TRD in over 500 segregating F2 populations of Arabidopsis thaliana using pooled reduced-representation genome sequencing. TRD was evident in up to a quarter of surveyed populations. Most populations exhibited distortion at only one genomic region, with some regions being repeatedly affected in multiple populations. Our results begin to elucidate the species-level architecture of biased transmission of genetic material in A. thaliana, and serve as a springboard for future studies into the biological basis of TRD in this species.

Lack of efficacy of monepantel against trichostrongyle nematodes in a UK sheep flock

Monepantel resistance was diagnosed during routine monitoring of the effectiveness of a farm’s roundworm control strategy. Weaned lambs had become ill thrifty and developed diarrhoea, despite the routine use of monepantel. This clinical presentation was caused by trichostrongylosis. The faecal egg count reduction was 76.7% (95% CI: 55.1–82.2%) following treatment with 2.5 mg/kg monepantel. Predominantly Trichostrongylus vitrinus along with small proportions of Oesophagostomum venulosum and Trichostrongylus vitrinus were identified by deep amplicon sequencing of pools of larvae recovered from pre and post monepantel treatment coprocultures and on postmortem examinations. The undifferentiated FECRT showed resistance to monepantel, but not to levamisole, ivermectin, or moxidectin. Examination of farm anthelmintic treatment and animal movement records suggested that treatments before movement onto silage aftermaths, putatively with low numbers of susceptible nematodes in refugia, may have placed a high selection pressure on monepantel resistance. Effective control of parasitic gastroenteritis using anthelmintic drugs is a prerequisite for sustainable sheep production. This case reiterates the need for care when combining anthelmintic treatments with movements to safe grazing, and the value of monitoring of anthelmintic efficacy as part of iterative planned animal health management.

Genomewide signatures of selection in Epichloë reveal candidate genes for host specialization.

Host specialization is a key process in ecological divergence and speciation of plant-associated fungi. The underlying determinants of host specialization are generally poorly understood, especially in endophytes, which constitute one of the most abundant components of the plant microbiome. We addressed the genetic basis of host specialization in two sympatric subspecies of grass-endophytic fungi from the Epichloë typhina complex: subsp. typhina and clarkii. The life cycle of these fungi entails unrestricted dispersal of gametes and sexual reproduction before infection of a new host, implying that the host imposes a selective barrier on viability of the progeny. We aimed to detect genes under divergent selection between subspecies, experiencing restricted gene flow due to adaptation to different hosts. Using pooled whole-genome sequencing data, we combined FST and DXY population statistics in genome scans and detected 57 outlier genes showing strong differentiation between the two subspecies. Genomewide analyses of nucleotide diversity (π), Tajima's D and dN/dS ratios indicated that these genes have evolved under positive selection. Genes encoding secreted proteins were enriched among the genes showing evidence of positive selection, suggesting that molecular plant-fungus interactions are strong drivers of endophyte divergence. We focused on five genes encoding secreted proteins, which were further sequenced in 28 additional isolates collected across Europe to assess genetic variation in a larger sample size. Signature of positive selection in these isolates and putative identification of pathogenic function supports our findings that these genes represent strong candidates for host specialization determinants in Epichloë endophytes. Our results highlight the role of secreted proteins as key determinants of host specialization.

The identification of growth, immune related genes and marker discovery through transcriptome in the yellow drum (Nibea albiflora).

Yellow drum (Nibea albiflora) is a commercially important marine fish, which is widely distributed in the coastal waters of China, Japan and Korea. Wild yellow drum resources have dramatically declined due to overfishing and ocean pollution. Genetic data can contribute to biodiversity conservation and protection. And molecular markers can play important roles in genetic breeding and aid in germplasm preservation in fish. In this study, 11 tissues (brain, heart, liver, kidney, muscle, head kidney, skin, fin, spleen, gonad and air bladder) were collected for pooled RNA sequencing. The unigenes were assembled using Trinity and EvidentialGene, and were then aligned to nr, nt, Swiss-Prot GO, KEGG, and KOG for annotation. Molecular markers (e.g. simple sequence repeat, SSR and single nucleotide polymorphism, SNP) were detected using MIcroSAtellite identification tool (MISA) and Genome Analysis Tool Kit (GATK). All clean reads were assembled into 109,209 transcripts, and 31,183 unigenes were generated after pruning and classifying, ranging from 201 to 19,857bp in length (1230bp in average), and 26,728 (85.7%) assembled unigenes had significant hits in public databases. Total of 27 and 103 unigenes were respectively identified as involved in growth- and immune-related pathways in the N. albiflora transcriptome. In addition, we identified a considerable quantity of molecular markers, including 11,484 SSRs and 56,186 SNPs. The growth- and immune-relevant genes and the molecular markers identified here provided a meaningful reference gene set and laid a foundation for future genetic selection and breeding for this species.

The Evolution of Polymorphic Hybrid Incompatibilities in House Mice.

Resolving the mechanistic and genetic bases of reproductive barriers between species is essential to understanding the evolutionary forces that shape speciation. Intrinsic hybrid incompatibilities are often treated as fixed between species, yet there can be considerable variation in the strength of reproductive isolation between populations. The extent and causes of this variation remain poorly understood in most systems. We investigated the genetic basis of variable hybrid male sterility (HMS) between two recently diverged subspecies of house mice, Mus musculus domesticus and Mus musculus musculus We found that polymorphic HMS has a surprisingly complex genetic basis, with contributions from at least five autosomal loci segregating between two closely related wild-derived strains of M. m. musculus One of the HMS-linked regions on chromosome 4 also showed extensive introgression among inbred laboratory strains and transmission ratio distortion (TRD) in hybrid crosses. Using additional crosses and whole genome sequencing of sperm pools, we showed that TRD was limited to hybrid crosses and was not due to differences in sperm motility between M. m. musculus strains. Based on these results, we argue that TRD likely reflects additional incompatibilities that reduce hybrid embryonic viability. In some common inbred strains of mice, selection against deleterious interactions appears to have unexpectedly driven introgression at loci involved in epistatic hybrid incompatibilities. The highly variable genetic basis to F1 hybrid incompatibilities between closely related mouse lineages argues that a thorough dissection of reproductive isolation will require much more extensive sampling of natural variation than has been commonly utilized in mice and other model systems.

In vivo insertion pool sequencing identifies virulence factors in a complex fungal-host interaction.

Large-scale insertional mutagenesis screens can be powerful genome-wide tools if they are streamlined with efficient downstream analysis, which is a serious bottleneck in complex biological systems. A major impediment to the success of next-generation sequencing (NGS)-based screens for virulence factors is that the genetic material of pathogens is often underrepresented within the eukaryotic host, making detection extremely challenging. We therefore established insertion Pool-Sequencing (iPool-Seq) on maize infected with the biotrophic fungus U. maydis. iPool-Seq features tagmentation, unique molecular barcodes, and affinity purification of pathogen insertion mutant DNA from in vivo-infected tissues. In a proof of concept using iPool-Seq, we identified 28 virulence factors, including 23 that were previously uncharacterized, from an initial pool of 195 candidate effector mutants. Because of its sensitivity and quantitative nature, iPool-Seq can be applied to any insertional mutagenesis library and is especially suitable for genetically complex setups like pooled infections of eukaryotic hosts.

Effective discovery of rare variants by pooled target capture sequencing: A comparative analysis with individually indexed target capture sequencing

Identification of all genetic variants associated with complex traits is one of the most important goals in modern human genetics. Genome-wide association studies (GWAS) have been successfully applied to identify common variants, which thus far explain only small portion of heritability. Interests in rare variants have been increasingly growing as an answer for this missing heritability. While next-generation sequencing allows detection of rare variants, its cost is still prohibitively high to sequence a large number of human DNA samples required for rare variant association studies. In this study, we evaluated the sensitivity and specificity of sequencing for pooled DNA samples of multiple individuals (Pool-seq) as a cost-effective and robust approach for rare variant discovery. We comparatively analyzed Pool-seq vs. individual-seq of indexed target capture of up to 960 genes in ∼1000 individuals, followed by independent genotyping validation studies. We found that Pool-seq was as effective and accurate as individual-seq in detecting rare variants and accurately estimating their minor allele frequencies (MAFs). Our results suggest that Pool-seq can be used as an efficient and cost-effective method for discovery of rare variants for population-based sequencing studies in individual laboratories.

Insight into higher-level phylogeny of Neuropterida: Evidence from secondary structures of mitochondrial rRNA genes and mitogenomic data.

It is well known that the rRNA structure information is important to assist phylogenetic analysis through identifying homologous positions to improve alignment accuracy. In addition, the secondary structure of some conserved motifs is highly stable among distantly related taxa, which can provide potentially informative characters for estimating phylogeny. In this paper, we applied the high-throughput pooled sequencing approach to the determination of neuropteran mitogenomes. Four complete mitogenome sequences were obtained: Micromus angulatus (Hemerobiidae), Chrysoperla nipponensis (Chrysopidae), Rapisma sp. (Ithonidae), and Thaumatosmylus sp. (Osmylidae). This allowed us to sample more complete mitochondrial RNA gene sequences. Secondary structure diagrams for the complete mitochondrial small and large ribosomal subunit RNA genes of eleven neuropterid species were predicted. Comparative analysis of the secondary structures indicated a closer relationship of Megaloptera and Neuroptera. This result was congruent with the resulting phylogeny inferred from sequence alignments of all 37 mitochondrial genes, namely the hypothesis of (Raphidioptera + (Megaloptera + Neuroptera)).

Exploring DNA variant segregation types in pooled genome sequencing enables effective mapping of weeping trait in Malus.

To unlock the power of next generation sequencing-based bulked segregant analysis in allele discovery in out-crossing woody species, and to understand the genetic control of the weeping trait, an F1 population from the cross 'Cheal's Weeping' × 'Evereste' was used to create two genomic DNA pools 'weeping' (17 progeny) and 'standard' (16 progeny). Illumina pair-end (2 × 151 bp) sequencing of the pools to a 27.1× (weeping) and a 30.4× (standard) genome (742.3 Mb) coverage allowed detection of 84562 DNA variants specific to 'weeping', 92148 specific to 'standard', and 173169 common to both pools. A detailed analysis of the DNA variant genotypes in the pools predicted three informative segregation types of variants: <lm×mm> (type I) in weeping pool-specific variants, and <lm×ll> (type II) and <hk×hk> (type III) in variants common to both pools, where the first allele is assumed to be weeping linked and the allele shown in bold is a variant in relation to the reference genome. Conducting variant allele frequency and density-based mappings revealed four genomic regions with a significant association with weeping: a major locus, Weeping (W), on chromosome 13 and others on chromosomes 10 (W2), 16 (W3), and 5 (W4). The results from type I variants were noisier and less certain than those from type II and type III variants, demonstrating that although type I variants are often the first choice, type II and type III variants represent an important source of DNA variants that can be exploited for genetic mapping in out-crossing woody species. Confirmation of the mapping of W and W2, investigation into their genetic interactions, and identification of expressed genes in the W and W2 regions provided insight into the genetic control of weeping and its expressivity in Malus.

A Small RNA Isolation and Sequencing Protocol and Its Application to Assay CRISPR RNA Biogenesis in Bacteria.

Next generation high-throughput sequencing has enabled sensitive and unambiguous analysis of RNA populations in cells. Here, we describe a method for isolation and strand-specific sequencing of small RNA pools from bacteria that can be multiplexed to accommodate multiple biological samples in a single experiment. Small RNAs are isolated by polyacrylamide gel electrophoresis and treated with T4 polynucleotide kinase. This allows for 3' adapter ligation to CRISPR RNAs, which don't have pre-existing 3'-OH ends. Pre-adenylated adapters are then ligated using T4 RNA ligase 1 in the absence of ATP and with a high concentration of polyethylene glycol (PEG). The 3' capture step enables precise determination of the 3' ends of diverse RNA molecules. Additionally, a random hexamer in the ligated adapter helps control for potential downstream amplification bias. Following reverse-transcription, the cDNA product is circularized and libraries are prepared by PCR. We show that the amplified library need not be visible by gel electrophoresis for efficient sequencing of the desired product. Using this method, we routinely prepare RNA sequencing libraries from minute amounts of purified small RNA. This protocol is tailored to assay for CRISPR RNA biogenesis in bacteria through sequencing of mature CRISPR RNAs, but can be used to sequence diverse classes of small RNAs. We also provide a fully worked example of our data processing pipeline, with instructions for running the provided scripts.

How to optimize the precision of allele and haplotype frequency estimates using pooled-sequencing data.

Sequencing pools of individuals rather than individuals separately reduces the costs of estimating allele frequencies at many loci in many populations. Theoretical and empirical studies show that sequencing pools comprising a limited number of individuals (typically fewer than 50) provides reliable allele frequency estimates, provided that the DNA pooling and DNA sequencing steps are carefully controlled. Unequal contributions of different individuals to the DNA pool and the mean and variance in sequencing depth both can affect the standard error of allele frequency estimates. To our knowledge, no study separately investigated the effect of these two factors on allele frequency estimates; so that there is currently no method to a priori estimate the relative importance of unequal individual DNA contributions independently of sequencing depth. We develop a new analytical model for allele frequency estimation that explicitly distinguishes these two effects. Our model shows that the DNA pooling variance in a pooled sequencing experiment depends solely on two factors: the number of individuals within the pool and the coefficient of variation of individual DNA contributions to the pool. We present a new method to experimentally estimate this coefficient of variation when planning a pooled sequencing design where samples are either pooled before or after DNA extraction. Using this analytical and experimental framework, we provide guidelines to optimize the design of pooled sequencing experiments. Finally, we sequence replicated pools of inbred lines of the plant Medicago truncatula and show that the predictions from our model generally hold true when estimating the frequency of known multilocus haplotypes using pooled sequencing.