High-density detection of restriction-site-associated DNA markers for rapid mapping of mutated loci in Neurospora

Key triggers of adaptive genetic variability of sessile oak [Q. petraea (Matt.) Liebl.] from the Balkan refugia: outlier detection and association of SNP loci from ddRAD-seq data

Knowledge on the genetic composition of Quercus petraea in south-eastern Europe is limited despite the species' significant role in the re-colonisation of Europe during the Holocene, and the diverse climate and physical geography of the region. Therefore, it is imperative to conduct research on adaptation in sessile oak to better understand its ecological significance in the region. While large sets of SNPs have been developed for the species, there is a continued need for smaller sets of SNPs that are highly informative about the possible adaptation to this varied landscape. By using double digest restriction site associated DNA sequencing data from our previous study, we mapped RAD-seq loci to the Quercus robur reference genome and identified a set of SNPs putatively related to drought stress-response. A total of 179 individuals from eighteen natural populations at sites covering heterogeneous climatic conditions in the southeastern natural distribution range of Q. petraea were genotyped. The detected highly polymorphic variant sites revealed three genetic clusters with a generally low level of genetic differentiation and balanced diversity among them but showed a north-southeast gradient. Selection tests showed nine outlier SNPs positioned in different functional regions. Genotype-environment association analysis of these markers yielded a total of 53 significant associations, explaining 2.4-16.6% of the total genetic variation. Our work exemplifies that adaptation to drought may be under natural selection in the examined Q. petraea populations.

Genomic landscape of the global oak phylogeny

The tree of life is highly reticulate, with the history of population divergence emerging from populations of gene phylogenies that reflect histories of introgression, lineage sorting and divergence. In this study, we investigate global patterns of oak diversity and test the hypothesis that there are regions of the oak genome that are broadly informative about phylogeny. We utilize fossil data and restriction-site associated DNA sequencing (RAD-seq) for 632 individuals representing nearly 250 Quercus species to infer a time-calibrated phylogeny of the world's oaks. We use a reversible-jump Markov chain Monte Carlo method to reconstruct shifts in lineage diversification rates, accounting for among-clade sampling biases. We then map the > 20 000 RAD-seq loci back to an annotated oak genome and investigate genomic distribution of introgression and phylogenetic support across the phylogeny. Oak lineages have diversified among geographic regions, followed by ecological divergence within regions, in the Americas and Eurasia. Roughly 60% of oak diversity traces back to four clades that experienced increases in net diversification, probably in response to climatic transitions or ecological opportunity. The strong support for the phylogeny contrasts with high genomic heterogeneity in phylogenetic signal and introgression. Oaks are phylogenomic mosaics, and their diversity may in fact depend on the gene flow that shapes the oak genome.

A RAD-sequencing approach to genome-wide marker discovery, genotyping, and phylogenetic inference in a diverse radiation of primates

Until recently, most phylogenetic and population genetics studies of nonhuman primates have relied on mitochondrial DNA and/or a small number of nuclear DNA markers, which can limit our understanding of primate evolutionary and population history. Here, we describe a cost-effective reduced representation method (ddRAD-seq) for identifying and genotyping large numbers of SNP loci for taxa from across the New World monkeys, a diverse radiation of primates that shared a common ancestor ~20-26 mya. We also estimate, for the first time, the phylogenetic relationships among 15 of the 22 currently-recognized genera of New World monkeys using ddRAD-seq SNP data using both maximum likelihood and quartet-based coalescent methods. Our phylogenetic analyses robustly reconstructed three monophyletic clades corresponding to the three families of extant platyrrhines (Atelidae, Pitheciidae and Cebidae), with Pitheciidae as basal within the radiation. At the genus level, our results conformed well with previous phylogenetic studies and provide additional information relevant to the problematic position of the owl monkey (Aotus) within the family Cebidae, suggesting a need for further exploration of incomplete lineage sorting and other explanations for phylogenetic discordance, including introgression. Our study additionally provides one of the first applications of next-generation sequencing methods to the inference of phylogenetic history across an old, diverse radiation of mammals and highlights the broad promise and utility of ddRAD-seq data for molecular primatology.

Next-Generation Sequencing-Based Approaches for Mutation Mapping and Identification in Caenorhabditis elegans

The use of next-generation sequencing (NGS) has revolutionized the way phenotypic traits are assigned to genes. In this review, we describe NGS-based methods for mapping a mutation and identifying its molecular identity, with an emphasis on applications in Caenorhabditis elegans In addition to an overview of the general principles and concepts, we discuss the main methods, provide practical and conceptual pointers, and guide the reader in the types of bioinformatics analyses that are required. Owing to the speed and the plummeting costs of NGS-based methods, mapping and cloning a mutation of interest has become straightforward, quick, and relatively easy. Removing this bottleneck previously associated with forward genetic screens has significantly advanced the use of genetics to probe fundamental biological processes in an unbiased manner.

Construction of Commercial Sweet Cherry Linkage Maps and QTL Analysis for Trunk Diameter

A cross between the sweet cherry (Prunus avium) cultivars ‘Wanhongzhu’ and ‘Lapins’ was performed to create a mapping population suitable for the construction of a linkage map. The specific-locus amplified fragment (SLAF) sequencing technique used as a single nucleotide polymorphism (SNP) discovery platform and generated 701 informative genotypic assays; these, along with 16 microsatellites (SSRs) and the incompatibility (S) gene, were used to build a map which comprised 8 linkage groups (LGs) and covered a genetic distance of 849.0 cM. The mean inter-marker distance was 1.18 cM and there were few gaps > 5 cM in length. Marker collinearity was maintained with the established peach genomic sequence. The map was used to show that trunk diameter (TD) is under the control of 4 loci, mapping to 3 different LGs. Different locus influenced TD at a varying stage of the tree’s development. The high density ‘W×L’ genetic linkage map has the potential to enable high-resolution identification of QTLs of agronomically relevant traits, and accelerate sweet cherry breeding.

High-density interspecific genetic maps of kiwifruit and the identification of sex-specific markers

Kiwifruit (Actinidia chinensis Planchon) is an important specialty fruit crop that suffers from narrow genetic diversity stemming from recent global commercialization and limited cultivar improvement. Here, we present high-density RAD-seq-based genetic maps using an interspecific F₁ cross between Actinidia rufa ‘MT570001’ and A. chinensis ‘Guihai No4’. The A. rufa (maternal) map consists of 2,426 single-nucleotide polymorphism (SNP) markers with a total length of 2,651 cM in 29 linkage groups (LGs) corresponding to the 29 chromosomes. The A. chinensis (paternal) map consists of 4,214 SNP markers over 3,142 cM in 29 LGs. Using these maps, we were able to anchor an additional 440 scaffolds from the kiwifruit draft genome assembly. Kiwifruit is functionally dioecious, which presents unique challenges for breeding and production. Three sex-specific simple sequence repeats (SSR) markers can be used to accurately sex type male and female kiwifruit in breeding programmes. The sex-determination region (SDR) in kiwifruit was narrowed to a 1-Mb subtelomeric region on chromosome 25. Localizing the SDR will expedite the discovery of genes controlling carpel abortion in males and pollen sterility in females.

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.

A framework phylogeny of the American oak clade based on sequenced RAD data

Previous phylogenetic studies in oaks (Quercus, Fagaceae) have failed to resolve the backbone topology of the genus with strong support. Here, we utilize next-generation sequencing of restriction-site associated DNA (RAD-Seq) to resolve a framework phylogeny of a predominantly American clade of oaks whose crown age is estimated at 23-33 million years old. Using a recently developed analytical pipeline for RAD-Seq phylogenetics, we created a concatenated matrix of 1.40 E06 aligned nucleotides, constituting 27,727 sequence clusters. RAD-Seq data were readily combined across runs, with no difference in phylogenetic placement between technical replicates, which overlapped by only 43-64% in locus coverage. 17% (4,715) of the loci we analyzed could be mapped with high confidence to one or more expressed sequence tags in NCBI Genbank. A concatenated matrix of the loci that BLAST to at least one EST sequence provides approximately half as many variable or parsimony-informative characters as equal-sized datasets from the non-EST loci. The EST-associated matrix is more complete (fewer missing loci) and has slightly lower homoplasy than non-EST subsampled matrices of the same size, but there is no difference in phylogenetic support or relative attribution of base substitutions to internal versus terminal branches of the phylogeny. We introduce a partitioned RAD visualization method (implemented in the R package RADami; http://cran.r-project.org/web/packages/RADami) to investigate the possibility that suboptimal topologies supported by large numbers of loci--due, for example, to reticulate evolution or lineage sorting--are masked by the globally optimal tree. We find no evidence for strongly-supported alternative topologies in our study, suggesting that the phylogeny we recover is a robust estimate of large-scale phylogenetic patterns in the American oak clade. Our study is one of the first to demonstrate the utility of RAD-Seq data for inferring phylogeny in a 23-33 million year-old clade.

Application of single nucleotide polymorphisms to non-model species: a technical review

Single nucleotide polymorphisms (SNPs) have gained wide use in humans and model species and are becoming the marker of choice for applications in other species. Technology that was developed for work in model species may provide useful tools for SNP discovery and genotyping in non-model organisms. However, SNP discovery can be expensive, labour intensive, and introduce ascertainment bias. In addition, the most efficient approaches to SNP discovery will depend on the research questions that the markers are to resolve as well as the focal species. We discuss advantages and disadvantages of several past and recent technologies for SNP discovery and genotyping and summarize a variety of SNP discovery and genotyping studies in ecology and evolution.

Deep sequencing strategies for mapping and identifying mutations from genetic screens

The development of next-generation sequencing technologies has enabled rapid and cost effective whole genome sequencing. This technology has allowed researchers to shortcut time-consuming and laborious methods used to identify nucleotide mutations in forward genetic screens in model organisms. However, causal mutations must still be mapped to a region of the genome so as to aid in their identification. This can be achieved simultaneously with deep sequencing through various methods. Here we discuss alternative deep sequencing strategies for simultaneously mapping and identifying causal mutations in Caenorhabditis elegans from mutagenesis screens. Focusing on practical considerations, such as the particular mutant phenotype obtained, this review aims to aid the reader in choosing which strategy to adopt to successfully clone their mutant.

Rapid development of molecular markers by next-generation sequencing linked to a gene conferring phomopsis stem blight disease resistance for marker-assisted selection in lupin (Lupinus angustifolius L.) breeding

Selection for phomopsis stem blight disease (PSB) resistance is one of the key objectives in lupin (Lupinus angustifolius L.) breeding programs. A cross was made between cultivar Tanjil (resistant to PSB) and Unicrop (susceptible). The progeny was advanced into F(8) recombinant inbred lines (RILs). The RIL population was phenotyped for PSB disease resistance. Twenty plants from the RIL population representing disease resistance and susceptibility was subjected to next-generation sequencing (NGS)-based restriction site-associated DNA sequencing on the NGS platform Solexa HiSeq2000, which generated 7,241 single nucleotide polymorphisms (SNPs). Thirty-three SNP markers showed the correlation between the marker genotypes and the PSB disease phenotype on the 20 representative plants, which were considered as candidate markers linked to a putative R gene for PSB resistance. Seven candidate markers were converted into sequence-specific PCR markers, which were designated as PhtjM1, PhtjM2, PhtjM3, PhtjM4, PhtjM5, PhtjM6 and PhtjM7. Linkage analysis of the disease phenotyping data and marker genotyping data on a F(8) population containing 187 RILs confirmed that all the seven converted markers were associated with the putative R gene within the genetic distance of 2.1 CentiMorgan (cM). One of the PCR markers, PhtjM3, co-segregated with the R gene. The seven established PCR markers were tested in the 26 historical and current commercial cultivars released in Australia. The numbers of "false positives" (showing the resistance marker allele band but lack of the putative R gene) for each of the seven PCR markers ranged from nil to eight. Markers PhtjM4 and PhtjM7 are recommended in marker-assisted selection for PSB resistance in the Australian national lupin breeding program due to its wide applicability on breeding germplasm and close linkage to the putative R gene. The results demonstrated that application of NGS technology is a rapid and cost-effective approach in development of markers for molecular plant breeding.

Two different high throughput sequencing approaches identify thousands of de novo genomic markers for the genetically depleted Bornean elephant

High throughput sequencing technologies are being applied to an increasing number of model species with a high-quality reference genome. The application and analyses of whole-genome sequence data in non-model species with no prior genomic information are currently under way. Recent sequencing technologies provide new opportunities for gathering genomic data in natural populations, laying the empirical foundation for future research in the field of conservation and population genomics. Here we present the case study of the Bornean elephant, which is the most endangered subspecies of Asian elephant and exhibits very low genetic diversity. We used two different sequencing platforms, the Roche 454 FLX (shotgun) and Illumina, GAIIx (Restriction site associated DNA, RAD) to evaluate the feasibility of the two methodologies for the discovery of de novo markers (single nucleotide polymorphism, SNPs and microsatellites) using low coverage data. Approximately, 6,683 (shotgun) and 14,724 (RAD) SNPs were detected within our elephant sequence dataset. Genotyping of a representative sample of 194 SNPs resulted in a SNP validation rate of ∼ 83 to 94% and 17% of the loci were polymorphic with a low diversity (H_o = 0.057). Different numbers of microsatellites were identified through shotgun (27,226) and RAD (868) techniques. Out of all di-, tri-, and tetra-microsatellite loci, 1,706 loci had sufficient flanking regions (shotgun) while only 7 were found with RAD. All microsatellites were monomorphic in the Bornean but polymorphic in another elephant subspecies. Despite using different sample sizes, and the well known differences in the two platforms used regarding sequence length and throughput, the two approaches showed high validation rate. The approaches used here for marker development in a threatened species demonstrate the utility of high throughput sequencing technologies as a starting point for the development of genomic tools in a non-model species and in particular for a species with low genetic diversity.

Application of next-generation sequencing for rapid marker development in molecular plant breeding: a case study on anthracnose disease resistance in Lupinus angustifolius L

Background

In the last 30 years, a number of DNA fingerprinting methods such as RFLP, RAPD, AFLP, SSR, DArT, have been extensively used in marker development for molecular plant breeding. However, it remains a daunting task to identify highly polymorphic and closely linked molecular markers for a target trait for molecular marker-assisted selection. The next-generation sequencing (NGS) technology is far more powerful than any existing generic DNA fingerprinting methods in generating DNA markers. In this study, we employed a grain legume crop Lupinus angustifolius (lupin) as a test case, and examined the utility of an NGS-based method of RAD (restriction-site associated DNA) sequencing as DNA fingerprinting for rapid, cost-effective marker development tagging a disease resistance gene for molecular breeding.

Results

Twenty informative plants from a cross of RxS (disease resistant x susceptible) in lupin were subjected to RAD single-end sequencing by multiplex identifiers. The entire RAD sequencing products were resolved in two lanes of the 16-lanes per run sequencing platform Solexa HiSeq2000. A total of 185 million raw reads, approximately 17 Gb of sequencing data, were collected. Sequence comparison among the 20 test plants discovered 8207 SNP markers. Filtration of DNA sequencing data with marker identification parameters resulted in the discovery of 38 molecular markers linked to the disease resistance gene Lanr1. Five randomly selected markers were converted into cost-effective, simple PCR-based markers. Linkage analysis using marker genotyping data and disease resistance phenotyping data on a F₈ population consisting of 186 individual plants confirmed that all these five markers were linked to the R gene. Two of these newly developed sequence-specific PCR markers, AnSeq3 and AnSeq4, flanked the target R gene at a genetic distance of 0.9 centiMorgan (cM), and are now replacing the markers previously developed by a traditional DNA fingerprinting method for marker-assisted selection in the Australian national lupin breeding program.

Conclusions

We demonstrated that more than 30 molecular markers linked to a target gene of agronomic trait of interest can be identified from a small portion (1/8) of one sequencing run on HiSeq2000 by applying NGS based RAD sequencing in marker development. The markers developed by the strategy described in this study are all co-dominant SNP markers, which can readily be converted into high throughput multiplex format or low-cost, simple PCR-based markers desirable for large scale marker implementation in plant breeding programs. The high density and closely linked molecular markers associated with a target trait help to overcome a major bottleneck for implementation of molecular markers on a wide range of germplasm in breeding programs. We conclude that application of NGS based RAD sequencing as DNA fingerprinting is a very rapid and cost-effective strategy for marker development in molecular plant breeding. The strategy does not require any prior genome knowledge or molecular information for the species under investigation, and it is applicable to other plant species.

Inferring phylogenies from RAD sequence data

Reduced-representation genome sequencing represents a new source of data for systematics, and its potential utility in interspecific phylogeny reconstruction has not yet been explored. One approach that seems especially promising is the use of inexpensive short-read technologies (e.g., Illumina, SOLiD) to sequence restriction-site associated DNA (RAD)--the regions of the genome that flank the recognition sites of restriction enzymes. In this study, we simulated the collection of RAD sequences from sequenced genomes of different taxa (Drosophila, mammals, and yeasts) and developed a proof-of-concept workflow to test whether informative data could be extracted and used to accurately reconstruct "known" phylogenies of species within each group. The workflow consists of three basic steps: first, sequences are clustered by similarity to estimate orthology; second, clusters are filtered by taxonomic coverage; and third, they are aligned and concatenated for "total evidence" phylogenetic analysis. We evaluated the performance of clustering and filtering parameters by comparing the resulting topologies with well-supported reference trees and we were able to identify conditions under which the reference tree was inferred with high support. For Drosophila, whole genome alignments allowed us to directly evaluate which parameters most consistently recovered orthologous sequences. For the parameter ranges explored, we recovered the best results at the low ends of sequence similarity and taxonomic representation of loci; these generated the largest supermatrices with the highest proportion of missing data. Applications of the method to mammals and yeasts were less successful, which we suggest may be due partly to their much deeper evolutionary divergence times compared to Drosophila (crown ages of approximately 100 and 300 versus 60 Mya, respectively). RAD sequences thus appear to hold promise for reconstructing phylogenetic relationships in younger clades in which sufficient numbers of orthologous restriction sites are retained across species.

Population genomic analysis of model and nonmodel organisms using sequenced RAD tags

The evolutionary processes of mutation, migration, genetic drift, and natural selection shape patterns of genetic variation among individuals, populations, and species, and they can do so differentially across genomes. The field of population genomics provides a comprehensive genome-scale view of these processes, even beyond traditional model organisms. Until recently, genome-wide studies of genetic variation have been prohibitively expensive. However, next-generation sequencing (NGS) technologies are revolutionizing the field of population genomics, allowing for genetic analysis at scales not previously possible even in organisms for which few genomic resources presently exist. To speed this revolution in evolutionary genetics, we and colleagues developed Restriction site Associated DNA (RAD) sequencing, a method that uses Illumina NGS to simultaneously type and score tens to hundreds of thousands of single nucleotide polymorphism (SNP) markers in hundreds of individuals for minimal investment of resources. The core molecular protocol is described elsewhere in this volume, which can be modified to suit a diversity of evolutionary genetic questions. In this chapter, we outline the conceptual framework of population genomics, relate genomic patterns of variation to evolutionary processes, and discuss how RAD sequencing can be used to study population genomics. In addition, we discuss bioinformatic considerations that arise from unique aspects of NGS data as compared to traditional marker based approaches, and we outline some general analytical approaches for RAD-seq and similar data, including a computational pipeline that we developed called Stacks. This software can be used for the analysis of RAD-seq data in organisms with and without a reference genome. Nonetheless, the development of analytical tools remains in its infancy, and further work is needed to fully quantify sampling variance and biases in these data types. As data-gathering technology continues to advance, our ability to understand genomic evolution in natural populations will be limited more by conceptual and analytical weaknesses than by the amount of molecular data.

RAD paired-end sequencing for local de novo assembly and SNP discovery in non-model organisms

Restriction-site Associated DNA (RAD) markers are rapidly becoming a standard for SNP discovery and genotyping studies even in organisms without a sequenced reference genome. It is difficult, however, to identify genes nearby RAD markers of interest or move from SNPs identified by RAD to a high-throughput genotyping assay. Paired-end sequencing of RAD fragments can alleviate these problems by generating a set of paired sequences that can be locally assembled into high-quality contigs up to 1 kb in length. These contigs can then be used for SNP identification, homology searching, or high-throughput assay primer design. In this chapter, we offer suggestions on how to design a RAD paired-end (RAD-PE) sequencing project and the protocol for creating paired-end RAD libraries suitable for Illumina sequencers.

Rapid mapping and identification of mutations in Caenorhabditis elegans by restriction site-associated DNA mapping and genomic interval pull-down sequencing

Forward genetic screens provide a powerful approach for inferring gene function on the basis of the phenotypes associated with mutated genes. However, determining the causal mutation by traditional mapping and candidate gene sequencing is often the rate-limiting step, especially when analyzing many mutants. We report two genomic approaches for more rapidly determining the identity of the affected genes in Caenorhabditis elegans mutants. First, we report our use of restriction site-associated DNA (RAD) polymorphism markers for rapidly mapping mutations after chemical mutagenesis and mutant isolation. Second, we describe our use of genomic interval pull-down sequencing (GIPS) to selectively capture and sequence megabase-sized portions of a mutant genome. Together, these two methods provide a rapid and cost-effective approach for positional cloning of C. elegans mutant loci, and are also applicable to other genetic model systems.

Mapping with RAD (restriction-site associated DNA) markers to rapidly identify QTL for stem rust resistance in Lolium perenne

A mapping population was created to detect quantitative trait loci (QTL) for resistance to stem rust caused by Puccinia graminis subsp. graminicola in Lolium perenne. A susceptible and a resistant plant were crossed to produce a pseudo-testcross population of 193 F(1) individuals. Markers were produced by the restriction-site associated DNA (RAD) process, which uses massively parallel and multiplexed sequencing of reduced-representation libraries. Additional simple sequence repeat (SSR) and sequence-tagged site (STS) markers were combined with the RAD markers to produce maps for the female (738 cM) and male (721 cM) parents. Stem rust phenotypes (number of pustules per plant) were determined in replicated greenhouse trials by inoculation with a field-collected, genetically heterogeneous population of urediniospores. The F(1) progeny displayed continuous distribution of phenotypes and transgressive segregation. We detected three resistance QTL. The most prominent QTL (qLpPg1) is located near 41 cM on linkage group (LG) 7 with a 2-LOD interval of 8 cM, and accounts for 30-38% of the stem rust phenotypic variance. QTL were detected also on LG1 (qLpPg2) and LG6 (qLpPg3), each accounting for approximately 10% of phenotypic variance. Alleles of loci closely linked to these QTL originated from the resistant parent for qLpPg1 and from both parents for qLpPg2 and qLpPg3. Observed quantitative nature of the resistance may be due to partial-resistance effects against all pathogen genotypes, or qualitative effects completely preventing infection by only some genotypes in the genetically mixed inoculum. RAD markers facilitated rapid construction of new genetic maps in this outcrossing species and will enable development of sequence-based markers linked to stem rust resistance in L. perenne.

Bulk segregant analysis followed by high-throughput sequencing reveals the Neurospora cell cycle gene, ndc-1, to be allelic with the gene for ornithine decarboxylase, spe-1

With the advent of high-throughput DNA sequencing, it is now straightforward and inexpensive to generate high-density small nucleotide polymorphism (SNP) maps. Here we combined high-throughput sequencing with bulk segregant analysis to expedite mutation mapping. The general map location of a mutation can be identified by a single backcross to a strain enriched in SNPs compared to a standard wild-type strain. Bulk segregant analysis simultaneously increases the likelihood of determining the precise nature of the mutation. We present here a high-density SNP map between Neurospora crassa Mauriceville-1-c (FGSC2225) and OR74A (FGSC2489), the strains most typically used by Neurospora researchers to carry out mapping crosses. We further have demonstrated the utility of the Mauriceville sequence and our approach by mapping the mutation responsible for the only existing temperature-sensitive (ts) cell cycle mutation in Neurospora, nuclear division cycle-1 (ndc-1). The single T-to-C point mutation maps to the gene encoding ornithine decarboxylase (ODC), spe-1 (NCU01271), and changes a Phe to a Ser residue within a highly conserved motif next to the catalytic site of the enzyme. By growth on spermidine and complementation with a wild-type spe-1 gene, we showed that the defect in spe-1 is responsible for the ts ndc-1 mutation. Based on our results, we propose changing ndc-1 to spe-1(ndc), which reflects that this mutation results in an ODC with a specific nuclear division defect.

Construction and application for QTL analysis of a Restriction Site Associated DNA (RAD) linkage map in barley

BACKGROUND

Linkage maps are an integral resource for dissection of complex genetic traits in plant and animal species. Canonical map construction follows a well-established workflow: an initial discovery phase where genetic markers are mined from a small pool of individuals, followed by genotyping of selected mapping populations using sets of marker panels. A newly developed sequence-based marker technology, Restriction site Associated DNA (RAD), enables synchronous single nucleotide polymorphism (SNP) marker discovery and genotyping using massively parallel sequencing. The objective of this research was to assess the utility of RAD markers for linkage map construction, employing barley as a model system. Using the published high density EST-based SNP map in the Oregon Wolfe Barley (OWB) mapping population as a reference, we created a RAD map using a limited set of prior markers to establish linakge group identity, integrated the RAD and prior data, and used both maps for detection of quantitative trait loci (QTL).

RESULTS

Using the RAD protocol in tandem with the Illumina sequence by synthesis platform, a total of 530 SNP markers were identified from initial scans of the OWB parental inbred lines--the "dominant" and "recessive" marker stocks--and scored in a 93 member doubled haploid (DH) mapping population. RAD sequence data from the structured population was converted into allele genotypes from which a genetic map was constructed. The assembled RAD-only map consists of 445 markers with an average interval length of 5 cM, while an integrated map includes 463 RAD loci and 2383 prior markers. Sequenced RAD markers are distributed across all seven chromosomes, with polymorphic loci emanating from both coding and noncoding regions in the Hordeum genome. Total map lengths are comparable and the order of common markers is identical in both maps. The same large-effect QTL for reproductive fitness traits were detected with both maps and the majority of these QTL were coincident with a dwarfing gene (ZEO) and the VRS1 gene, which determines the two-row and six-row germplasm groups of barley.

CONCLUSIONS

We demonstrate how sequenced RAD markers can be leveraged to produce high quality linkage maps for detection of single gene loci and QTLs. By combining SNP discovery and genotyping into parallel sequencing events, RAD markers should be a useful molecular breeding tool for a range of crop species. Expected improvements in cost and throughput of second and third-generation sequencing technologies will enable more powerful applications of the sequenced RAD marker system, including improvements in de novo genome assembly, development of ultra-high density genetic maps and association mapping.

Restriction Site Tiling Analysis: accurate discovery and quantitative genotyping of genome-wide polymorphisms using nucleotide arrays

A method for the simultaneous identification of polymorphic loci and the quantitative genotyping of thousands of loci in individuals is presented.

High-throughput genotype data can be used to identify genes important for local adaptation in wild populations, phenotypes in lab stocks, or disease-related traits in human medicine. Here we advance microarray-based genotyping for population genomics with Restriction Site Tiling Analysis. The approach simultaneously discovers polymorphisms and provides quantitative genotype data at 10,000s of loci. It is highly accurate and free from ascertainment bias. We apply the approach to uncover genomic differentiation in the purple sea urchin.

Identification of DIM-7, a protein required to target the DIM-5 H3 methyltransferase to chromatin

Functionally distinct chromatin domains are delineated by distinct posttranslational modifications of histones, and in some organisms by differences in DNA methylation. Proper establishment and maintenance of chromatin domains is critical but not well understood. We previously demonstrated that heterochromatin in the filamentous fungus Neurospora crassa is marked by cytosine methylation directed by trimethylated Lysine 9 on histone H3 (H3K9me3). H3K9me3 is the product of the DIM-5 Lysine methyltransferase and is recognized by a protein complex containing heterochromatin protein-1 and the DIM-2 DNA methyltransferase. To identify additional components that control the establishment and function of DNA methylation and heterochromatin, we built a strain harboring two selectable reporter genes that are silenced by DNA methylation and employed this strain to select for mutants that are defective in DNA methylation (dim). We report a previously unidentified gene (dim-7) that is essential for H3K9me3 and DNA methylation. DIM-7 homologs are found only in fungi and are highly divergent. We found that DIM-7 interacts with DIM-5 in vivo and demonstrated that a conserved domain near the N terminus of DIM-7 is required for its stability. In addition, we found that DIM-7 is essential for recruitment of DIM-5 to form heterochromatin.

Genome-wide investigation of reproductive isolation in experimental lineages and natural species of Neurospora: identifying candidate regions by microarray-based genotyping and mapping

Inherent incompatibilities between genetic components from genomes of different species may cause intrinsic reproductive isolation. In evolution experiments designed to instigate speciation in laboratory populations of the filamentous fungus Neurospora, we previously discovered a pair of incompatibility loci (dfe and dma) that interact negatively to cause severe defects in sexual reproduction. Here we show that the dfe-dma incompatibility also is a significant cause of genetic isolation between two naturally occurring species of Neurospora (N. crassa and N. intermedia). The strong incompatibility interaction has a simple genetic basis (two biallelic loci) and antagonistic epistasis occurs between heterospecific alleles only, consistent with the Dobzhansky-Muller model of genic incompatibility. We developed microarray-based, restriction-site associated DNA (RAD) markers that identified approximately 1500 polymorphisms between the genomes of the two species, and constructed the first interspecific physical map of Neurospora. With this new mapping resource, the approximate genomic locations of the incompatibility loci were determined using three different approaches: genome scanning, bulk-segregant analyses, and introgression. These population, quantitative, and classical genetics methods concordantly identified two candidate regions, narrowing the search for each incompatibility locus to only approximately 2% of the nuclear genome. This study demonstrates how advances in high-throughput, genome-wide genotyping can be applied to mapping reproductive isolation genes and speciation research.

Tools for fungal proteomics: multifunctional neurospora vectors for gene replacement, protein expression and protein purification

The completion of genome-sequencing projects for a number of fungi set the stage for detailed investigations of proteins. We report the generation of versatile expression vectors for detection and isolation of proteins and protein complexes in the filamentous fungus Neurospora crassa. The vectors, which can be adapted for other fungi, contain C- or N-terminal FLAG, HA, Myc, GFP, or HAT-FLAG epitope tags with a flexible poly-glycine linker and include sequences for targeting to the his-3 locus in Neurospora. To introduce mutations at native loci, we also made a series of knock-in vectors containing epitope tags followed by the selectable marker hph (resulting in hygromycin resistance) flanked by two loxP sites. We adapted the Cre/loxP system for Neurospora, allowing the selectable marker hph to be excised by introduction of Cre recombinase into a strain containing a knock-in cassette. Additionally, a protein purification method was developed on the basis of the HAT-FLAG tandem affinity tag system, which was used to purify HETEROCHROMATIN PROTEIN 1 (HP1) and associated proteins from Neurospora. As expected on the basis of yeast two-hybrid and co-immunoprecipitation (Co-IP) experiments, the Neurospora DNA methyltransferase DIM-2 was found in a complex with HP1. Features of the new vectors allowed for verification of an interaction between HP1 and DIM-2 in vivo by Co-IP assays on proteins expressed either from their native loci or from the his-3 locus.

A high-density single nucleotide polymorphism map for Neurospora crassa

We report the discovery and validation of a set of single nucleotide polymorphisms (SNPs) between the reference Neurospora crassa strain Oak Ridge and the Mauriceville strain (FGSC 2555), of sufficient density to allow fine mapping of most loci. Sequencing of Mauriceville cDNAs and alignment to the completed genomic sequence of the Oak Ridge strain identified 19,087 putative SNPs. Of these, a subset was validated by cleaved amplified polymorphic sequence (CAPS), a simple and robust PCR-based assay that reliably distinguishes between SNP alleles. Experimental confirmation resulted in the development of 250 CAPS markers distributed evenly over the genome. To demonstrate the applicability of this map, we used bulked segregant analysis followed by interval mapping to locate the csp-1 mutation to a narrow region on LGI. Subsequently, we refined mapping resolution to 74 kbp by developing additional markers, resequenced the candidate gene, NCU02713.3, in the mutant background, and phenocopied the mutation by gene replacement in the WT strain. Together, these techniques demonstrate a generally applicable and straightforward approach for the isolation of novel genes from existing mutants. Data on both putative and validated SNPs are deposited in a customized public database at the Broad Institute, which encourages augmentation by community users.

Genome-wide high throughput analysis of DNA methylation in eukaryotes

Cytosine methylation is the quintessential epigenetic mark. Two well-established methods, bisulfite sequencing and methyl-DNA immunoprecipitation (MeDIP) lend themselves to the genome-wide analysis of DNA methylation by high throughput sequencing. Here we provide an overview and brief review of these methods. We summarize our experience with MeDIP followed by high throughput Illumina/Solexa sequencing, exemplified by the analysis of the methylated fraction of the Neurospora crassa genome ("methylome"). We provide detailed methods for DNA isolation, processing and the generation of in vitro libraries for Illumina/Solexa sequencing. We discuss potential problems in the generation of sequencing libraries. Finally, we provide an overview of software that is appropriate for the analysis of high throughput sequencing data generated by Illumina/Solexa-type sequencing by synthesis, with a special emphasis on approaches and applications that can generate more accurate depictions of sequence reads that fall in repeated regions of a chosen reference genome.

Rapid SNP discovery and genetic mapping using sequenced RAD markers

Single nucleotide polymorphism (SNP) discovery and genotyping are essential to genetic mapping. There remains a need for a simple, inexpensive platform that allows high-density SNP discovery and genotyping in large populations. Here we describe the sequencing of restriction-site associated DNA (RAD) tags, which identified more than 13,000 SNPs, and mapped three traits in two model organisms, using less than half the capacity of one Illumina sequencing run. We demonstrated that different marker densities can be attained by choice of restriction enzyme. Furthermore, we developed a barcoding system for sample multiplexing and fine mapped the genetic basis of lateral plate armor loss in threespine stickleback by identifying recombinant breakpoints in F₂ individuals. Barcoding also facilitated mapping of a second trait, a reduction of pelvic structure, by in silico re-sorting of individuals. To further demonstrate the ease of the RAD sequencing approach we identified polymorphic markers and mapped an induced mutation in Neurospora crassa. Sequencing of RAD markers is an integrated platform for SNP discovery and genotyping. This approach should be widely applicable to genetic mapping in a variety of organisms.

Array-based high-throughput DNA markers for crop improvement

The last two decades have witnessed a remarkable activity in the development and use of molecular markers both in animal and plant systems. This activity started with low-throughput restriction fragment length polymorphisms and culminated in recent years with single nucleotide polymorphisms (SNPs), which are abundant and uniformly distributed. Although the latter became the markers of choice for many, their discovery needed previous sequence information. However, with the availability of microarrays, SNP platforms have been developed, which allow genotyping of thousands of markers in parallel. Besides SNPs, some other novel marker systems, including single feature polymorphisms, diversity array technology and restriction site-associated DNA markers, have also been developed, where array-based assays have been utilized to provide for the desired ultra-high throughput and low cost. These microarray-based markers are the markers of choice for the future and are already being used for construction of high-density maps, quantitative trait loci (QTL) mapping (including expression QTLs) and genetic diversity analysis with a limited expense in terms of time and money. In this study, we briefly describe the characteristics of these array-based marker systems and review the work that has already been done involving development and use of these markers, not only in simple eukaryotes like yeast, but also in a variety of seed plants with simple or complex genomes.

Simple sequence repeats in Neurospora crassa: distribution, polymorphism and evolutionary inference

Background

Simple sequence repeats (SSRs) have been successfully used for various genetic and evolutionary studies in eukaryotic systems. The eukaryotic model organism Neurospora crassa is an excellent system to study evolution and biological function of SSRs.

Results

We identified and characterized 2749 SSRs of 963 SSR types in the genome of N. crassa. The distribution of tri-nucleotide (nt) SSRs, the most common SSRs in N. crassa, was significantly biased in exons. We further characterized the distribution of 19 abundant SSR types (AST), which account for 71% of total SSRs in the N. crassa genome, using a Poisson log-linear model. We also characterized the size variation of SSRs among natural accessions using Polymorphic Index Content (PIC) and ANOVA analyses and found that there are genome-wide, chromosome-dependent and local-specific variations. Using polymorphic SSRs, we have built linkage maps from three line-cross populations.

Conclusion

Taking our computational, statistical and experimental data together, we conclude that 1) the distributions of the SSRs in the sequenced N. crassa genome differ systematically between chromosomes as well as between SSR types, 2) the size variation of tri-nt SSRs in exons might be an important mechanism in generating functional variation of proteins in N. crassa, 3) there are different levels of evolutionary forces in variation of amino acid repeats, and 4) SSRs are stable molecular markers for genetic studies in N. crassa.