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(2) 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.     

Fine-scale mapping in Neurospora crassa by using genome-wide knockout strains

Fine-scale genetic mapping is often hindered by the lack of adequate markers surrounding the locus of interest. In the filamentous ascomycete Neurospora crassa, the genome has been sequenced and an effort has been made to generate genome-wide deletion strains for the entire gene set. Accordingly, the hygromycin-resistant marker in each deletion strain can be used as a mapping locus in a classical three-point cross, along with the mapping target and a standard marker. We have demonstrated the feasibility of this fine-scale mapping approach in N. crassa by refining the location of r(Sk-2).     

Positional cloning without a genome map: using 'Targeted RFLP Subtraction' to isolate dense markers tightly linked to the regA locus of Volvox carteri.

The ability to isolate genes defined by mutant phenotypes has fueled the rapid progress in understanding basic biological mechanisms and the causes of inherited diseases. Positional cloning, a commonly used method for isolating genes corresponding to mutations, is most efficiently applied to the small number of model organisms for which high resolution genetic maps exist. We demonstrate a new and generally applicable positional cloning method that obviates the need for a genetic map. The technique is based on Restriction Fragment Length Polymorphism (RFLP) Subtraction, a method that isolates RFLP markers spanning an entire genome. The new method, Targeted RFLP Subtraction (TRS), isolates markers from a specific region by combining RFLP Subtraction with a phenotypic pooling strategy. We used TRS to directly isolate dense markers tightly linked to the regA gene of the eukaryotic green alga Volvox. As a generally applicable method for saturating a small targeted region with DNA markers, TRS should facilitate gene isolation from diverse organisms and accelerate the process of physically mapping specific regions in preparation for sequence analysis.     

General method for cloning Neurospora crassa nuclear genes by complementation of mutants.

We have developed a sib selection procedure for cloning Neurospora crassa nuclear genes by complementation of mutants. This procedure takes advantage of a modified N. crassa transformation procedure that gives as many as 10,000 to 50,000 stable transformants per microgram of DNA with recombinant plasmids containing the N. crassa qa-2+ gene. Here, we describe the use of the sib selection procedure to clone genes corresponding to auxotrophic mutants, nic-1 and inl. The identities of the putative clones were confirmed by mapping their chromosomal locations in standard genetic crosses and using restriction site polymorphisms as genetic markers. Because we can obtain very high N. crassa transformation frequencies, cloning can be accomplished with as few as five subdivisions of an N. crassa genomic library. The sib selection procedure should, for the first time, permit the cloning of any gene corresponding to an N. crassa mutant for which an appropriate selection can be devised. Analogous procedures may be applicable to other filamentous fungi before the development of operational shuttle vectors.     

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

The wealth of sequence information available for Neurospora crassa and other fungi has greatly facilitated evolutionary and molecular analyses of this group. Although "reverse" genetics, in which genes are first identified by their sequence rather than by their mutant phenotypes, serves as a valuable new approach for elucidating biological processes, classical "forward" genetic analysis is still extremely useful. Unfortunately, mapping mutations and identifying the corresponding genes has typically been slow and laborious. To facilitate forward genetics in Neurospora, we have adapted microarray-based restriction-site-associated DNA (RAD) mapping for use with N. crassa oligonucleotide microarrays. This technique was used to simultaneously detect an unprecedented number of genomewide restriction site polymorphisms from two N. crassa strains: Mauriceville and Oak Ridge. Furthermore, RAD mapping was used to quickly map a previously unknown gene, defective in methylation-7 (dim-7).     

Centromere-linked microsatellite markers for linkage groups 3, 4, 6, 7, 13, and 20 of zebrafish (Danio rerio)

A large number of interesting mutations affecting development and organogenesis have been identified through genetic screens in zebrafish. Mapping of these mutations to a chromosomal region can be rapidly accomplished using half-tetrad analysis. However, knowledge of centromere-linked markers on every chromosome is essential to this mapping method. Centromeres on all 25 linkage groups have been mapped on the RAPD zebrafish genetic map. However, species specificity and the lack of codominance make RAPD markers less practical for mapping than microsatellite-based markers. On the microsatellite-based genetic map, centromere-linked markers have been identified for 19 linkage groups. No direct evidence has been published linking microsatellite markers to the centromeres of linkage groups 3, 4, 6, 7, 13, and 20. Therefore, we compared the microsatellite-based genetic map with the RAPD map to identify markers most likely linked to the centromeres of these 6 linkage groups. These candidate markers were tested for potential centromere linkage using four panels of half-tetrad embryos derived by early-pressure treatment of eggs from four different female zebrafish. We have identified microsatellite markers for linkage groups 3, 4, 6, 7, 13, and 20 to within 1.7 cM of their centromeres. These markers will greatly facilitate the rapid mapping of mutations in zebrafish by half-tetrad analysis.     

玉米相对饱和遗传连锁图谱构建与一种新的AFLPs共显性分析方法探讨

利用一个F2作图群体(X178×B73),首先构建了一个含有130个SSRs的玉米连锁框架图,然后用119个AFLPs位点增加图谱密度,得到一个全长1659·3cM,标记间平均间距6·66cM的玉米相对饱和连锁图。同时,对SSRs和AFLPs的一些遗传特性进行了分析,探讨了AFLP标记进行共显性分析的一种新方法。分析表明SSRs和AFLPs分子标记具有多态性和可靠性高等特点,是构建高密度分子标记遗传连锁图的有效技术。加密的玉米遗传连锁图谱为比较基因组研究、数量性状位点(quantitativetraitloci,QTLs)克隆、杂种优势机理研究以及标记辅助选择等提供了技术基础。     

Development and Application of a Whole-Genome Simple Sequence Repeat Panel for High-Throughput Genotyping in Soybean

Among commonly applied molecular markers, simple sequence repeats (SSRs, or microsatellites) possess advantages such as a high level of polymorphism and codominant pattern of inheritance at individual loci. To facilitate systematic and rapid genetic mapping in soybean, we designed a genotyping panel comprised 304 SSR markers selected for allelic diversity and chromosomal location so as to provide wide coverage. Most primer pairs for the markers in the panel were redesigned to yield amplicons of 80–600 bp in multiplex polymerase chain reaction (PCR) and fluorescence-based sequencer analysis, and they were labelled with one of four different fluorescent dyes. Multiplex PCR with sets of six to eight primer pairs per reaction generated allelic data for 283 of the 304 SSR loci in three different mapping populations, with the loci mapping to the same positions as previously determined. Four SSRs on each chromosome were analysed for allelic diversity in 87 diverse soybean germplasms with four-plex PCR. These 80 loci showed an average allele number and polymorphic information content value of 14.8 and 0.78, respectively. The high level of polymorphism, ease of analysis, and high accuracy of the SSR genotyping panel should render it widely applicable to soybean genetics and breeding.     

Common sites of retroviral integration in mouse hematopoietic tumors identified by high-throughput,single nucleotide polymorphism-based mapping and bacterial artificial chromosome hybridization

Retroviral insertional mutagenesis in mouse hematopoietic tumors provides a powerful cancer gene discovery tool. Here, we describe a high-throughput, single nucleotide polymorphism (SNP)-based method, for mapping retroviral integration sites cloned from mouse tumors, and a bacterial artificial chromosome (BAC) hybridization method, for localizing these retroviral integration sites to common sites of retroviral integration (CISs). Several new CISs were identified, including one CIS that mapped near Notch1, a gene that has been causally associated with human T-cell tumors. This mapping method is applicable to many different species, including ones where few genetic markers and little genomic sequence information are available. It can also be used to map endogenous proviruses.     

Plant microsatellite genotyping with 4-color fluorescent detection using multiple-tailed primers

We extended the concept of fluorescent microsatellite genotyping with a single-universal tailed primer to the simultaneous use of three different tailed primers to allow multiplexed 4-color detection for medium throughput genotyping of plant species. The method was tested on Eucalyptus DNA samples using three forward primer sequences of human microsatellite markers labeled with different fluorescent dyes. The robustness of the method was tested for the simultaneous detection and genetic analysis of microsatellites in a genetic mapping experiment. This method allows reliable and cost-effective genotyping with the same level of multiplexing attained in regular microsatellite fluorescent detection assays. Besides the enhanced quality of the genotypic data provided by the fluorescent detection method when compared to colorimetric ones, the economy brought about by this method becomes greater with an increasing number of microsatellite markers. This method has been particularly useful for genotyping populations of several tropical tree species addressing community-wide population genetics and conservation questions.     

基于转录组数据高通量发掘扶桑绵粉蚧微卫星引物

【目的】扶桑绵粉蚧Phenacoccus solenopsis Tinsley是我国重要的检疫性害虫。简单重复序列(simple sequence repeat,SSR)研究在遗传图谱和物理图谱的构建、分子标记辅助育种、品种鉴定、基因定位、遗传多样性、动植物分类和进化等方面具有重要意义。筛选的SSR引物将为扶桑绵粉蚧遗传多样性分析、进化分析及入侵生物学等奠定基础。【方法】利用高通量搜索的方法对扶桑绵粉蚧转录组中28 120条unigenes的数据进行搜索。【结果】共找到1 781个SSR位点。扶桑绵粉蚧转录组中SSRs的主要重复类型是单核苷酸重复,占SSR总数的89.44%;其次是三核苷酸重复,占SSR总数的7.52%。单核苷酸重复里主要是A/T基序,占了总量的87.42%。基于筛选的SSRs,运用Primer 3软件进行引物的批量设计,共有481个unigenes成功设计引物,共设计出1 228对引物。【结论】研究表明利用扶桑绵粉蚧转录组数据开发SSR标记是可行的,本研究开发的引物将为扶桑绵粉蚧遗传多样性分析、进化分析及入侵生物学等奠定基础。     

A genetic linkage map and comparative genome analysis of common carp (Cyprinus carpio L.) using microsatellites and SNPs

A genetic linkage map is a powerful research tool for mapping traits of interest and is essential to understanding genome evolution. The aim of this study is to provide an expanded genetic linkage map of common carp to effectively carry out quantitative trait loci analysis and conduct comparative mapping analysis between lineages. Here, we constructed a genetic linkage map of common carp (Cyprinus carpio L.) using microsatellite and single-nucleotide polymorphism (SNP) markers in a 159 sibling family. A total of 246 microsatellites and 306 SNP polymorphic markers were genotyped in this family. Linkage analysis using JoinMap 4.0 organized 427 markers (186 microsatellites and 241 SNPs) to 50 linkage groups, ranging in size from 1.4 to 130.1 cM. Each group contained 2-30 markers. The linkage map covered a genetic distance of 2,039.2 cM and the average interval for markers within the linkage groups was approximately 6.4 cM. In addition, comparative genome analysis within five model teleost fish revealed a high percentage (74.7%) of conserved loci corresponding to zebrafish chromosomes. In most cases, each zebrafish chromosome comprised two common carp linkage groups. The comparative analysis also revealed independent chromosome rearrangements in common carp and zebrafish. The linkage map will be of great assistance in mapping genes of interest and serve as a reference to approach comparative mapping and enable further insights into the comprehensive investigations of genome evolution of common carp.     

A genetic linkage map of the model legume Lotus japonicus and strategies for fast mapping of new loci

A genetic map for the model legume Lotus japonicus has been developed. The F(2) mapping population was established from an interspecific cross between L. japonicus and L. filicaulis. A high level of DNA polymorphism between these parents was the source of markers for linkage analysis and the map is based on a framework of amplified fragment length polymorphism (AFLP) markers. Additional markers were generated by restriction fragment length polymorphism (RFLP) and sequence-specific PCR. A total of 524 AFLP markers, 3 RAPD markers, 39 gene-specific markers, 33 microsatellite markers, and six recessive symbiotic mutant loci were mapped. This genetic map consists of six linkage groups corresponding to the six chromosomes in L. japonicus. Fluorescent in situ hybridization (FISH) with selected markers aligned the linkage groups to chromosomes as described in the accompanying article by Pedrosa et al. 2002(this issue). The length of the linkage map is 367 cM and the average marker distance is 0.6 cM. Distorted segregation of markers was found in certain sections of the map and linkage group I could be assembled only by combining colormapping and cytogenetics (FISH). A fast method to position genetic loci employing three AFLP primer combinations yielding 89 markers was developed and evaluated by mapping three symbiotic loci, Ljsym1, Ljsym5, and Ljhar1-3.     

Fast and Accurate Construction of Ultra-Dense Consensus Genetic Maps Using Evolution Strategy Optimization

Our aim was to develop a fast and accurate algorithm for constructing consensus genetic maps for chip-based SNP genotyping data with a high proportion of shared markers between mapping populations. Chip-based genotyping of SNP markers allows producing high-density genetic maps with a relatively standardized set of marker loci for different mapping populations. The availability of a standard high-throughput mapping platform simplifies consensus analysis by ignoring unique markers at the stage of consensus mapping thereby reducing mathematical complicity of the problem and in turn analyzing bigger size mapping data using global optimization criteria instead of local ones. Our three-phase analytical scheme includes automatic selection of ~100-300 of the most informative (resolvable by recombination) markers per linkage group, building a stable skeletal marker order for each data set and its verification using jackknife re-sampling, and consensus mapping analysis based on global optimization criterion. A novel Evolution Strategy optimization algorithm with a global optimization criterion presented in this paper is able to generate high quality, ultra-dense consensus maps, with many thousands of markers per genome. This algorithm utilizes "potentially good orders" in the initial solution and in the new mutation procedures that generate trial solutions, enabling to obtain a consensus order in reasonable time. The developed algorithm, tested on a wide range of simulated data and real world data (Arabidopsis), outperformed two tested state-of-the-art algorithms by mapping accuracy and computation time.     

Use of the Caulobacter crescentus genome sequence to develop a method for systematic genetic mapping

The functional analysis of sequenced genomes will be facilitated by the development of tools for the rapid mapping of mutations. We have developed a systematic approach to genetic mapping in Caulobacter crescentus that is based on bacteriophage-mediated transduction of strategically placed antibiotic resistance markers. The genomic DNA sequence was used to identify sites distributed evenly around the chromosome at which plasmids could be nondisruptively integrated. DNA fragments from these sites were amplified by PCR and cloned into a kanamycin-resistant (Kan(r)) suicide vector. Delivery of these plasmids into C. crescentus resulted in integration via homologous recombination. A set of 41 strains containing Kan(r) markers at 100-kb intervals was thereby generated. These strains serve as donors for generalized transduction using bacteriophage phiCr30, which can transduce at least 120 kb of DNA. Transductants are selected with kanamycin and screened for loss of the mutant phenotype to assess linkage between the marker and the site of the mutation. The dependence of cotransduction frequency on sequence distance was evaluated using several markers and mutant strains. With these data as a standard, previously unmapped mutations were readily localized to DNA sequence intervals equivalent to less than 1% of the genome. Candidate genes within the interval were then examined further by subcloning and complementation analysis. Mutations resulting in sensitivity to ampicillin, in nutritional auxotrophies, or temperature-sensitive growth were mapped. This approach to genetic mapping should be applicable to other bacteria with sequenced genomes for which generalized transducing phage are available.     

Informative SSR markers for commercial variety discrimination in watermelon (Citrullus lanatus)

SSR markers were used for variety discrimination and genetic assessment in watermelon varieties. Genetic characterization of 49 watermelon varieties was investigated using 30 SSR markers developed from melon and watermelon. A total of 121 polymorphic amplified fragments were obtained by using 30 SSR markers. The average polymorphism information content (PIC) was 0.502 ranging from 0.223 to 0.800. One hundred twenty one SSR loci were used to calculate Jaccard’s distance coefficients for unweighted pair group method using the arithmetic averages (UPGMA) cluster analysis. A clustering group of varieties, based on the results of SSR analysis, were categorized into 5 major groups corresponding to morphological traits. Inheritance mode of 2 SSR markers was investigated to F₁ plants and F₂ populations of 2 crosses. Parental alleles were transmitted from F1 plants and F2 populations. Therefore, these marker sets may prove to be effectively applicable to genetic assessment of germplasm, genome mapping, and fingerprinting of watermelon varieties.     

Precision Mapping of Quantitative Trait Loci

Adequate separation of effects of possible multiple linked quantitative trait loci (QTLs) on mapping QTLs is the key to increasing the precision of QTL mapping. A new method of QTL mapping is proposed and analyzed in this paper by combining interval mapping with multiple regression. The basis of the proposed method is an interval test in which the test statistic on a marker interval is made to be unaffected by QTLs located outside a defined interval. This is achieved by fitting other genetic markers in the statistical model as a control when performing interval mapping. Compared with the current QTL mapping method (i.e., the interval mapping method which uses a pair or two pairs of markers for mapping QTLs), this method has several advantages. (1) By confining the test to one region at a time, it reduces a multiple dimensional search problem (for multiple QTLs) to a one dimensional search problem. (2) By conditioning linked markers in the test, the sensitivity of the test statistic to the position of individual QTLs is increased, and the precision of QTL mapping can be improved. (3) By selectively and simultaneously using other markers in the analysis, the efficiency of QTL mapping can be also improved. The behavior of the test statistic under the null hypothesis and appropriate critical value of the test statistic for an overall test in a genome are discussed and analyzed. A simulation study of QTL mapping is also presented which illustrates the utility, properties, advantages and disadvantages of the method.     

A microsatellite marker based linkage map of tobacco

We report the first linkage map of tobacco (Nicotiana tabacum L.) generated through microsatellite markers. The microsatellite markers were predominantly derived from genomic sequences of the Tobacco Genome Initiative (TGI) through bioinformatics screening for microsatellite motives. A total of 684 primer pairs were screened for functionality in a panel of 16 tobacco lines. Of those, 637 primer pairs were functional. Potential parents for mapping populations were evaluated for their polymorphism level through genetic similarity analysis. The similarity analysis revealed that the known groups of tobacco varieties (Burley, Flue-cured, Oriental and Dark) form distinct clusters. A mapping population, based on a cross between varieties Hicks Broad Leaf and Red Russian, and consisting of 186 F2 individuals, was selected for mapping. A total of 282 functional microsatellite markers were polymorphic in this population and 293 loci could be mapped together with the morphological trait flower color. Twenty-four tentative linkage groups spanning 1,920 cM could be identified. This map will provide the basis for the genetic mapping of traits in tobacco and for further analyses of the tobacco genome. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Genetic mapping of the non-nodulation phenotype of the mutant MN-1008 in tetraploid alfalfa (Medicago sativa)

Abstract. Roots of the non-nodulating Medicago sativa mutant MN-1008 neither undergo root-hair curling, cortical cell division nor any of the early molecular events that accompany nodule initiation and development following rhizobial infection or treatment with Nod factor. These observations suggested that the mutation(s) impaired a pivotal function in Nod factor perception or in the signal transduction pathway. In this paper we show that the genetic lesion conditioning the recessive non-nodulation phenotype in the tetraploid alfalfa mutant MN-1008 can be localized to a single region on LG5 of the M. sativa genetic map. This conclusion is based on genetic analyses conducted at the tetraploid level, involving both segregation analysis and genetic mapping of the trait with respect to molecular DNA markers. The genetic mapping of the Nod(-) phenotype was performed in a segregating tetraploid F2 population, taking advantage of the availability of an advanced genetic map for diploid alfalfa. Two tightly linked flanking markers have been identified which will facilitate the physical mapping and cloning of the gene(s) that underlie(s) the non-nodulation phenotype.