Rapid identification of heterozygous mutations in Drosophila melanogaster using genomic capture sequencing

The Barley HvWRKY6 Transcription Factor Is Required for Resistance Against Pyrenophora teres f. teres

Barley is an important cereal crop worldwide because of its use in the brewing and distilling industry. However, adequate supplies of quality malting barley are threatened by global climate change due to drought in some regions and excess precipitation in others, which facilitates epidemics caused by fungal pathogens. The disease net form net blotch caused by the necrotrophic fungal pathogen Pyrenophora teres f. teres (Ptt) has emerged as a global threat to barley production and diverse populations of Ptt have shown a capacity to overcome deployed genetic resistances. The barley line CI5791 exhibits remarkably effective resistance to diverse Ptt isolates from around the world that maps to two major QTL on chromosomes 3H and 6H. To identify genes involved in this effective resistance, CI5791 seed were γ-irradiated and two mutants, designated CI5791-γ3 and CI5791-γ8, with compromised Ptt resistance were identified from an M₂ population. Phenotyping of CI5791-γ3 and -γ8 × Heartland F₂ populations showed three resistant to one susceptible segregation ratios and CI5791-γ3 × -γ8 F₁ individuals were susceptible, thus these independent mutants are in a single allelic gene. Thirty-four homozygous mutant (susceptible) CI5791-γ3 × Heartland F₂ individuals, representing 68 recombinant gametes, were genotyped via PCR genotype by sequencing. The data were used for single marker regression mapping placing the mutation on chromosome 3H within an approximate 75 cM interval encompassing the 3H CI5791 resistance QTL. Sequencing of the mutants and wild-type (WT) CI5791 genomic DNA following exome capture identified independent mutations of the HvWRKY6 transcription factor located on chromosome 3H at ∼50.7 cM, within the genetically delimited region. Post transcriptional gene silencing of HvWRKY6 in barley line CI5791 resulted in Ptt susceptibility, confirming that it functions in NFNB resistance, validating it as the gene underlying the mutant phenotypes. Allele analysis and transcript regulation of HvWRKY6 from resistant and susceptible lines revealed sequence identity and upregulation upon pathogen challenge in all genotypes analyzed, suggesting a conserved transcription factor is involved in the defense against the necrotrophic pathogen. We hypothesize that HvWRKY6 functions as a conserved signaling component of defense mechanisms that restricts Ptt growth in barley.

Mutagenesis of Puccinia graminis f. sp. tritici and Selection of Gain-of-Virulence Mutants

Wheat stem rust caused by the fungus Puccinia graminis f. sp. tritici (Pgt), is regaining prominence due to the recent emergence of virulent isolates and epidemics in Africa, Europe and Central Asia. The development and deployment of wheat cultivars with multiple stem rust resistance (Sr) genes stacked together will provide durable resistance. However, certain disease resistance genes can suppress each other or fail in particular genetic backgrounds. Therefore, the function of each Sr gene must be confirmed after incorporation into an Sr-gene stack. This is difficult when using pathogen disease assays due to epistasis from recognition of multiple avirulence (Avr) effectors. Heterologous delivery of single Avr effectors can circumvent this limitation, but this strategy is currently limited by the paucity of cloned Pgt Avrs. To accelerate Avr gene cloning, we outline a procedure to develop a mutant population of Pgt spores and select for gain-of-virulence mutants. We used ethyl methanesulphonate (EMS) to mutagenize urediniospores and create a library of > 10,000 independent mutant isolates that were combined into 16 bulks of ~658 pustules each. We sequenced random mutants and determined the average mutation density to be 1 single nucleotide variant (SNV) per 258 kb. From this, we calculated that a minimum of three independently derived gain-of-virulence mutants is required to identify a given Avr gene. We inoculated the mutant library onto plants containing Sr43, Sr44, or Sr45 and obtained 9, 4, and 14 mutants with virulence toward Sr43, Sr44, or Sr45, respectively. However, only mutants identified on Sr43 and Sr45 maintained their virulence when reinolculated onto the lines from which they were identified. We further characterized 8 mutants with virulence toward Sr43. These also maintained their virulence profile on the stem rust international differential set containing 20 Sr genes, indicating that they were most likely not accidental contaminants. In conclusion, our method allows selecting for virulent mutants toward targeted resistance (R) genes. The development of a mutant library from as little as 320 mg spores creates a resource that enables screening against several R genes without the need for multiple rounds of spore multiplication and mutagenesis.

Analysis of transient hypermorphic activity of E(spl)D during R8 specification

Drosophila atonal (ato) is required for the specification of founding R8 photoreceptors during retinal development. ato is regulated via dual eye-specific enhancers; ato-3’ is subject to initial induction whereas 5’-ato facilitates Notch-mediated autoregulation. Notch is further utilized to induce bHLH repressors of the E(spl) locus to restrict Ato from its initial broad expression to individual cells. Although Notch operates in two, distinct phases, it has remained unclear how the two phases maintain independence from one another. The difference in these two phases has attributed to the hypothesized delayed expression of E(spl). However, immunofluorescence data indicate that E(spl) are expressed during early Ato patterning, suggesting a more sophisticated underlying mechanism. To probe this mechanism, we provide evidence that although E(spl) exert no influence on ato-3’, E(spl) repress 5’-ato and deletion of the E(spl) locus elicits precocious 5’-ato activity. Thus, E(spl) imposes a delay to the timing in which Ato initiates autoregulation. We next sought to understand this finding in the context of E(spl)^D, which encodes a dysregulated variant of E(spl)M8 that perturbs R8 patterning, though, as previously reported, only in conjunction with the mutant receptor N^spl. We established a genetic interaction between E(spl)^D and roughened eye (roe), a known modulator of Notch signaling in retinogenesis. This link further suggests a dosage-dependence between E(spl) and the proneural activators Ato and Sens, as indicated via interaction assays in which E(spl)^D renders aberrant R8 patterning in conjunction with reduced proneural dosage. In total, the biphasicity of Notch signaling relies, to some degree, on the post-translational regulation of individual E(spl) members and, importantly, that post-translational regulation is likely necessary to modulate the level of E(spl) activity throughout the progression of Ato expression.

Targeted capture in evolutionary and ecological genomics

The rapid expansion of next-generation sequencing has yielded a powerful array of tools to address fundamental biological questions at a scale that was inconceivable just a few years ago. Various genome-partitioning strategies to sequence select subsets of the genome have emerged as powerful alternatives to whole-genome sequencing in ecological and evolutionary genomic studies. High-throughput targeted capture is one such strategy that involves the parallel enrichment of preselected genomic regions of interest. The growing use of targeted capture demonstrates its potential power to address a range of research questions, yet these approaches have yet to expand broadly across laboratories focused on evolutionary and ecological genomics. In part, the use of targeted capture has been hindered by the logistics of capture design and implementation in species without established reference genomes. Here we aim to (i) increase the accessibility of targeted capture to researchers working in nonmodel taxa by discussing capture methods that circumvent the need of a reference genome, (ii) highlight the evolutionary and ecological applications where this approach is emerging as a powerful sequencing strategy and (iii) discuss the future of targeted capture and other genome-partitioning approaches in the light of the increasing accessibility of whole-genome sequencing. Given the practical advantages and increasing feasibility of high-throughput targeted capture, we anticipate an ongoing expansion of capture-based approaches in evolutionary and ecological research, synergistic with an expansion of whole-genome sequencing.

Getting started in mapping-by-sequencing

Next-generation sequencing (NGS) technologies allow the cost-effective sequencing of whole genomes and have expanded the scope of genomics to novel applications, such as the genome-wide characterization of intraspecific polymorphisms and the rapid mapping and identification of point mutations. Next-generation sequencing platforms, such as the Illumina HiSeq2000 platform, are now commercially available at affordable prices and routinely produce an enormous amount of sequence data, but their wide use is often hindered by a lack of knowledge on how to manipulate and process the information produced. In this review, we focus on the strategies that are available to geneticists who wish to incorporate these novel approaches into their research but who are not familiar with the necessary bioinformatic concepts and computational tools. In particular, we comprehensively summarize case studies where the use of NGS technologies has led to the identification of point mutations, a strategy that has been dubbed "mapping-by-sequencing", and review examples from plants and other model species such as Caenorhabditis elegans, Saccharomyces cerevisiae, and Drosophila melanogaster. As these technologies are becoming cheaper and more powerful, their use is also expanding to allow mutation identification in species with larger genomes, such as many crop plants.

Using next-generation sequencing to isolate mutant genes from forward genetic screens

The long-lasting success of forward genetic screens relies on the simple molecular basis of the characterized phenotypes, which are typically caused by mutations in single genes. Mapping the location of causal mutations using genetic crosses has traditionally been a complex, multistep procedure, but next-generation sequencing now allows the rapid identification of causal mutations at single-nucleotide resolution even in complex genetic backgrounds. Recent advances of this mapping-by-sequencing approach include methods that are independent of reference genome sequences, genetic crosses and any kind of linkage information, which make forward genetics amenable for species that have not been considered for forward genetic screens so far.

Genetic screens to identify new Notch pathway mutants in Drosophila

Notch signaling controls a wide range of developmental processes, including proliferation, apoptosis, and cell fate specification during both development and adult tissue homeostasis. The functional versatility of the Notch signaling pathway is tightly linked with the complexity of its regulation in different cellular contexts. To unravel the complexity of Notch signaling, it is important to identify the different components of the Notch signaling pathway. A powerful strategy to accomplish this task is based on genetic screens. Given that the developmental context of signaling is important, these screens should be customized to specific cell populations or tissues. Here, I describe how to perform F1 clonal forward genetic screens in Drosophila to identify novel components of the Notch signaling pathway. These screens combine a classical EMS (ethyl methanesulfonate) chemical mutagenesis protocol along with clonal analysis via FRT-mediated mitotic recombination. These F1 clonal screens allow rapid phenotypic screening within clones of mutant cells induced at specific developmental stages and in tissues of interest, bypassing the pleiotropic effects of isolated mutations. More importantly, since EMS mutations have been notoriously difficult to map to specific genes in the past, I briefly discuss mapping methods that allow rapid identification of the causative mutations.

Natural variations and genome-wide association studies in crop plants

Natural variants of crops are generated from wild progenitor plants under both natural and human selection. Diverse crops that are able to adapt to various environmental conditions are valuable resources for crop improvements to meet the food demands of the increasing human population. With the completion of reference genome sequences, the advent of high-throughput sequencing technology now enables rapid and accurate resequencing of a large number of crop genomes to detect the genetic basis of phenotypic variations in crops. Comprehensive maps of genome variations facilitate genome-wide association studies of complex traits and functional investigations of evolutionary changes in crops. These advances will greatly accelerate studies on crop designs via genomics-assisted breeding. Here, we first discuss crop genome studies and describe the development of sequencing-based genotyping and genome-wide association studies in crops. We then review sequencing-based crop domestication studies and offer a perspective on genomics-driven crop designs.

Forward genetics by genome sequencing reveals that rapid cyanide release deters insect herbivory of Sorghum bicolor

Whole genome sequencing has allowed rapid progress in the application of forward genetics in model species. In this study, we demonstrated an application of next-generation sequencing for forward genetics in a complex crop genome. We sequenced an ethyl methanesulfonate-induced mutant of Sorghum bicolor defective in hydrogen cyanide release and identified the causal mutation. A workflow identified the causal polymorphism relative to the reference BTx623 genome by integrating data from single nucleotide polymorphism identification, prior information about candidate gene(s) implicated in cyanogenesis, mutation spectra, and polymorphisms likely to affect phenotypic changes. A point mutation resulting in a premature stop codon in the coding sequence of dhurrinase2, which encodes a protein involved in the dhurrin catabolic pathway, was responsible for the acyanogenic phenotype. Cyanogenic glucosides are not cyanogenic compounds but their cyanohydrins derivatives do release cyanide. The mutant accumulated the glucoside, dhurrin, but failed to efficiently release cyanide upon tissue disruption. Thus, we tested the effects of cyanide release on insect herbivory in a genetic background in which accumulation of cyanogenic glucoside is unchanged. Insect preference choice experiments and herbivory measurements demonstrate a deterrent effect of cyanide release capacity, even in the presence of wild-type levels of cyanogenic glucoside accumulation. Our gene cloning method substantiates the value of (1) a sequenced genome, (2) a strongly penetrant and easily measurable phenotype, and (3) a workflow to pinpoint a causal mutation in crop genomes and accelerate in the discovery of gene function in the postgenomic era.

Organ-specific gene expression: the bHLH protein Sage provides tissue specificity to Drosophila FoxA

FoxA transcription factors play major roles in organ-specific gene expression, regulating, for example, glucagon expression in the pancreas, GLUT2 expression in the liver, and tyrosine hydroxylase expression in dopaminergic neurons. Organ-specific gene regulation by FoxA proteins is achieved through cooperative regulation with a broad array of transcription factors with more limited expression domains. Fork head (Fkh), the sole Drosophila FoxA family member, is required for the development of multiple distinct organs, yet little is known regarding how Fkh regulates tissue-specific gene expression. Here, we characterize Sage, a bHLH transcription factor expressed exclusively in the Drosophila salivary gland (SG). We show that Sage is required for late SG survival and normal tube morphology. We find that many Sage targets, identified by microarray analysis, encode SG-specific secreted cargo, transmembrane proteins, and the enzymes that modify these proteins. We show that both Sage and Fkh are required for the expression of Sage target genes, and that co-expression of Sage and Fkh is sufficient to drive target gene expression in multiple cell types. Sage and Fkh drive expression of the bZip transcription factor Senseless (Sens), which boosts expression of Sage-Fkh targets, and Sage, Fkh and Sens colocalize on SG chromosomes. Importantly, expression of Sage-Fkh target genes appears to simply add to the tissue-specific gene expression programs already established in other cell types, and Sage and Fkh cannot alter the fate of most embryonic cell types even when expressed early and continuously.

Genetically engineered animal models for in vivo target identification and validation in oncology

In vitro approaches using human cancer cell lines aimed to identify and validate oncology targets, have pinpointed a number of key targets and signalling pathways which control cell growth and cell death. However, tumors are more than insular masses of proliferating cancer cells. Instead they are complex tissues composed of multiple distinct cell types that participate in homotypic and heterotypic interactions and depend upon each other for their growth. Therefore, many targets in oncology need to be validated in the context of the whole animal. This review provides an overview on how animal models can be generated and used for target identification and validation in vivo.

Whole Genome Sequencing and a New Bioinformatics Platform Allow for Rapid Gene Identification in D. melanogaster EMS Screens

Forward genetic screens in Drosophila melanogaster using ethyl methanesulfonate (EMS) mutagenesis are a powerful approach for identifying genes that modulate specific biological processes in an in vivo setting. The mapping of genes that contain randomly-induced point mutations has become more efficient in Drosophila thanks to the maturation and availability of many types of genetic tools. However, classic approaches to gene mapping are relatively slow and ultimately require extensive Sanger sequencing of candidate chromosomal loci. With the advent of new high-throughput sequencing techniques, it is increasingly efficient to directly re-sequence the whole genome of model organisms. This approach, in combination with traditional chromosomal mapping, has the potential to greatly simplify and accelerate mutation identification in mutants generated in EMS screens. Here we show that next-generation sequencing (NGS) is an accurate and efficient tool for high-throughput sequencing and mutation discovery in Drosophila melanogaster. As a test case, mutant strains of Drosophila that exhibited long-term survival of severed peripheral axons were identified in a forward EMS mutagenesis. All mutants were recessive and fell into a single lethal complementation group, which suggested that a single gene was responsible for the protective axon degenerative phenotype. Whole genome sequencing of these genomes identified the underlying gene ect4. To improve the process of genome wide mutation identification, we developed Genomes Management Application (GEM.app, https://genomics.med.miami.edu), a graphical online user interface to a custom query framework. Using a custom GEM.app query, we were able to identify that each mutant carried a unique non-sense mutation in the gene ect4 (dSarm), which was recently shown by Osterloh et al. to be essential for the activation of axonal degeneration. Our results demonstrate the current advantages and limitations of NGS in Drosophila and we introduce GEM.app as a simple yet powerful genomics analysis tool for the Drosophila community. At a current cost of <$1,000 per genome, NGS should thus become a standard gene discovery tool in EMS induced genetic forward screens.

Exon capture and bulk segregant analysis: rapid discovery of causative mutations using high-throughput sequencing

Background

Exome sequencing has transformed human genetic analysis and may do the same for other vertebrate model systems. However, a major challenge is sifting through the large number of sequence variants to identify the causative mutation for a given phenotype. In models like Xenopus tropicalis, an incomplete and occasionally incorrect genome assembly compounds this problem. To facilitate cloning of X. tropicalis mutants identified in forward genetic screens, we sought to combine bulk segregant analysis and exome sequencing into a single step.

Results

Here we report the first use of exon capture sequencing to identify mutations in a non-mammalian, vertebrate model. We demonstrate that bulk segregant analysis coupled with exon capture sequencing is not only able to identify causative mutations but can also generate linkage information, facilitate the assembly of scaffolds, identify misassembles, and discover thousands of SNPs for fine mapping.

Conclusion

Exon capture sequencing and bulk segregant analysis is a rapid, inexpensive method to clone mutants identified in forward genetic screens. With sufficient meioses, this method can be generalized to any model system with a genome assembly, polished or unpolished, and in the latter case, it also provides many critical genomic resources.

A diploid wheat TILLING resource for wheat functional genomics

BACKGROUND

Triticum monococcum L., an A genome diploid einkorn wheat, was the first domesticated crop. As a diploid, it is attractive genetic model for the study of gene structure and function of wheat-specific traits. Diploid wheat is currently not amenable to reverse genetics approaches such as insertion mutagenesis and post-transcriptional gene silencing strategies. However, TILLING offers a powerful functional genetics approach for wheat gene analysis.

RESULTS

We developed a TILLING population of 1,532 M2 families using EMS as a mutagen. A total of 67 mutants were obtained for the four genes studied. Waxy gene mutation frequencies are known to be 1/17.6 - 34.4 kb DNA in polyploid wheat TILLING populations. The T. monococcum diploid wheat TILLING population had a mutation frequency of 1/90 kb for the same gene. Lignin biosynthesis pathway genes- COMT1, HCT2, and 4CL1 had mutation frequencies of 1/86 kb, 1/92 kb and 1/100 kb, respectively. The overall mutation frequency of the diploid wheat TILLING population was 1/92 kb.

CONCLUSION

The mutation frequency of a diploid wheat TILLING population was found to be higher than that reported for other diploid grasses. The rate, however, is lower than tetraploid and hexaploid wheat TILLING populations because of the higher tolerance of polyploids to mutations. Unlike polyploid wheat, most mutants in diploid wheat have a phenotype amenable to forward and reverse genetic analysis and establish diploid wheat as an attractive model to study gene function in wheat. We estimate that a TILLING population of 5, 520 will be needed to get a non-sense mutation for every wheat gene of interest with 95% probability.

Using Drosophila melanogaster as a Model for Genotoxic Chemical Mutational Studies with a New Program, SnpSift

This paper describes a new program SnpSift for filtering differential DNA sequence variants between two or more experimental genomes after genotoxic chemical exposure. Here, we illustrate how SnpSift can be used to identify candidate phenotype-relevant variants including single nucleotide polymorphisms, multiple nucleotide polymorphisms, insertions, and deletions (InDels) in mutant strains isolated from genome-wide chemical mutagenesis of Drosophila melanogaster. First, the genomes of two independently isolated mutant fly strains that are allelic for a novel recessive male-sterile locus generated by genotoxic chemical exposure were sequenced using the Illumina next-generation DNA sequencer to obtain 20- to 29-fold coverage of the euchromatic sequences. The sequencing reads were processed and variants were called using standard bioinformatic tools. Next, SnpEff was used to annotate all sequence variants and their potential mutational effects on associated genes. Then, SnpSift was used to filter and select differential variants that potentially disrupt a common gene in the two allelic mutant strains. The potential causative DNA lesions were partially validated by capillary sequencing of polymerase chain reaction-amplified DNA in the genetic interval as defined by meiotic mapping and deletions that remove defined regions of the chromosome. Of the five candidate genes located in the genetic interval, the Pka-like gene CG12069 was found to carry a separate pre-mature stop codon mutation in each of the two allelic mutants whereas the other four candidate genes within the interval have wild-type sequences. The Pka-like gene is therefore a strong candidate gene for the male-sterile locus. These results demonstrate that combining SnpEff and SnpSift can expedite the identification of candidate phenotype-causative mutations in chemically mutagenized Drosophila strains. This technique can also be used to characterize the variety of mutations generated by genotoxic chemicals.

Library-free methylation sequencing with bisulfite padlock probes

Targeted quantification of DNA methylation allows for interrogation of the most informative loci across many samples quickly and cost-effectively. Here we report improved bisulfite padlock probes (BSPPs) with a design algorithm to generate efficient padlock probes, a library-free protocol that dramatically reduces sample-preparation cost and time and is compatible with automation, and an efficient bioinformatics pipeline to accurately obtain both methylation levels and genotypes from sequencing of bisulfite-converted DNA.

Whole-Genome Sequencing to Identify Mutants and Polymorphisms in Chlamydomonas reinhardtii

Whole-genome sequencing (WGS) provides a new platform for the identification of mutations that produce a mutant phenotype. We used Illumina sequencing to identify the mutational profile of three Chlamydomonas reinhardtii mutant strains. The three strains have more than 38,000 changes from the reference genome. NG6 is aflagellate and maps to 269 kb with only one nonsynonymous change; the V(12)E mutation falls in the FLA8 gene. Evidence that NG6 is a fla8 allele comes from swimming revertants that are either true or pseudorevertants. NG30 is aflagellate and maps to 458 kb that has six nonsynonomous changes. Evidence that NG30 has a causative nonsense allele in IFT80 comes from rescue of the nonswimming phenotype with a fragment bearing only this gene. This gene has been implicated in Jeune asphyxiating thoracic dystrophy. Electron microscopy of ift80-1 (NG30) shows a novel basal body phenotype. A bar or cap is observed over the distal end of the transition zone, which may be an intermediate in preparing the basal body for flagellar assembly. In the acetate-requiring mutant ac17, we failed to find a nonsynonymous change in the 676 kb mapped region, which is incompletely assembled. In these strains, 43% of the changes occur on two of the 17 chromosomes. The excess on chromosome 6 surrounds the mating-type locus, which has numerous rearrangements and suppressed recombination, and the changes extend beyond the mating-type locus. Unexpectedly, chromosome 16 shows an unexplained excess of single nucleotide polymorphisms and indels. Overall, WGS in combination with limited mapping allows fast and accurate identification of point mutations in Chlamydomonas.

Genetic manipulation of genes and cells in the nervous system of the fruit fly

Research in the fruit fly Drosophila melanogaster has led to insights in neural development, axon guidance, ion channel function, synaptic transmission, learning and memory, diurnal rhythmicity, and neural disease that have had broad implications for neuroscience. Drosophila is currently the eukaryotic model organism that permits the most sophisticated in vivo manipulations to address the function of neurons and neuronally expressed genes. Here, we summarize many of the techniques that help assess the role of specific neurons by labeling, removing, or altering their activity. We also survey genetic manipulations to identify and characterize neural genes by mutation, overexpression, and protein labeling. Here, we attempt to acquaint the reader with available options and contexts to apply these methods.

Exome-wide DNA capture and next generation sequencing in domestic and wild species

BACKGROUND

Gene-targeted and genome-wide markers are crucial to advance evolutionary biology, agriculture, and biodiversity conservation by improving our understanding of genetic processes underlying adaptation and speciation. Unfortunately, for eukaryotic species with large genomes it remains costly to obtain genome sequences and to develop genome resources such as genome-wide SNPs. A method is needed to allow gene-targeted, next-generation sequencing that is flexible enough to include any gene or number of genes, unlike transcriptome sequencing. Such a method would allow sequencing of many individuals, avoiding ascertainment bias in subsequent population genetic analyses.We demonstrate the usefulness of a recent technology, exon capture, for genome-wide, gene-targeted marker discovery in species with no genome resources. We use coding gene sequences from the domestic cow genome sequence (Bos taurus) to capture (enrich for), and subsequently sequence, thousands of exons of B. taurus, B. indicus, and Bison bison (wild bison). Our capture array has probes for 16,131 exons in 2,570 genes, including 203 candidate genes with known function and of interest for their association with disease and other fitness traits.

RESULTS

We successfully sequenced and mapped exon sequences from across the 29 autosomes and X chromosome in the B. taurus genome sequence. Exon capture and high-throughput sequencing identified thousands of putative SNPs spread evenly across all reference chromosomes, in all three individuals, including hundreds of SNPs in our targeted candidate genes.

CONCLUSIONS

This study shows exon capture can be customized for SNP discovery in many individuals and for non-model species without genomic resources. Our captured exome subset was small enough for affordable next-generation sequencing, and successfully captured exons from a divergent wild species using the domestic cow genome as reference.

Multiplexed shotgun genotyping for rapid and efficient genetic mapping

We present a new approach to genotyping based on multiplexed shotgun sequencing that can identify recombination breakpoints in a large number of individuals simultaneously at a resolution sufficient for most mapping purposes, such as quantitative trait locus (QTL) mapping and mapping of induced mutations. We first describe a simple library construction protocol that uses just 10 ng of genomic DNA per individual and makes the approach accessible to any laboratory with standard molecular biology equipment. Sequencing this library results in a large number of sequence reads widely distributed across the genomes of multiplexed bar-coded individuals. We develop a Hidden Markov Model to estimate ancestry at all genomic locations in all individuals using these data. We demonstrate the utility of the approach by mapping a dominant marker allele in D. simulans to within 105 kb of its true position using 96 F1-backcross individuals genotyped in a single lane on an Illumina Genome Analyzer. We further demonstrate the utility of our method by genetically mapping more than 400 previously unassembled D. simulans contigs to linkage groups and by evaluating the quality of targeted introgression lines. At this level of multiplexing and divergence between strains, our method allows estimation of recombination breakpoints to a median of 38-kb intervals. Our analysis suggests that higher levels of multiplexing and/or use of strains with lower levels of divergence are practicable.