Combining DNA pooling with selective recombinant genotyping for increased efficiency in fine mapping

One of the key steps in positional cloning and marker-aided selection is to identify marker(s) tightly linked to the target gene (i.e., fine mapping). Selective genotyping such as selective recombinant genotyping (SRG) is commonly used in fine mapping for cost-saving. To further decrease genotyping effort and rapidly screen for tightly linked markers, we propose here a combined DNA pooling and SRG strategy. A two-stage pooled genotyping can be used for identifying recombinants between a pair of flanking markers more efficiently, and a joint use of bulked DNA analysis and two-stage pooling can also save cost for genotyping recombinants. The combined DNA pooling and SRG strategy can further be extended to fine mapping for polygenic traits. The numerical results based on hypothetical scenarios and an illustrative application to fine mapping of a mutant gene, called xl(t), in rice suggest that the proposed strategy can remarkably reduce genotyping amount compared with the conventional SRG.     

Fine-scale comparative genetic and physical mapping supports map-based cloning strategies for the self-incompatibility loci of perennial ryegrass (<Emphasis Type="Italic">Lolium perenne</Emphasis> L.)

Perennial ryegrass is an obligate outbreeding pasture grass of the Poaceae family, with a two-locus (S and Z) gametophytic self-incompatibility (SI) mechanism. This system has provided a major obstacle to targeted varietal development, and enhanced knowledge is expected to support more efficient breeding strategies. Comparative genetics and physical mapping approaches have been developed to permit molecular cloning of the SI genes. SI gene-linked genetic markers based on heterologous cDNA restriction fragment length polymorphisms (RFLPs) and homologous genomic DNA-derived simple sequence repeats (SSRs) were converted to single nucleotide polymorphism (SNP) format for efficient genotyping. Genetic mapping identified the location of SI loci and demonstrated macrosynteny between related grass species. S- and Z-linked bacterial artificial chromosome (BAC) clones were sequenced using massively parallel pyrosequencing technology to provide the first physical mapping data for Poaceae SI loci. The sequence assembly process suggested a lower prevalence of middle repetitive sequences in the Z locus region and hence precedence for positional cloning strategy. In silico mapping using data from rice, Brachypodium distachyon and Sorghum revealed high sequence conservation in the vicinity of the Z locus region between SI and self-compatible (SC) grass species. Physical mapping identified a total of nine genes encoded in the Z locus region. Expression profiling and nucleotide diversity assessment identified two Z-linked genes, LpTC116908 and LpDUF247, as plausible candidates for the male and female determinants of the S-Z SI system.     

An AFLP-based genome-wide mapping strategy

To efficiently determine the chromosomal location of phenotypic mutants, we designed a genome-wide mapping strategy that can be used in any crop for which a dense AFLP (Amplified Fragment Length Polymorphism) map is available or can be made. The AFLP technique is particularly suitable to initiate map-based cloning projects because it detects many markers per reaction. First a standard set of AFLP primer combinations that results in a framework of AFLP markers well dispersed over the genome is selected. These primer combinations are applied to a limited number of mutant individuals from a segregating population to register linkage and non-linkage of the AFLP markers to the gene-of-interest. Further delineation of the area of interest is accomplished by analyzing the remaining recombinants and additional mutant individuals with AFLP markers that lie within the identified region. We illustrate the usefulness of the method by mapping three rotunda (ron) leaf-form mutant loci of Arabidopsis thaliana and show that in the initial phase of map-based cloning projects a 400–600 kb interval can be identified for the average mutant locus within a few weeks. Once such an area is identified and before initiating the more time-consuming fine-mapping procedure, it is essential to examine publicly available databases for candidate genes and known mutants in the identified region. The 390-kb interval on chromosome 4 that harbors the ron2 mutation, also carries a known flower mutant, leunig (lug); upon crossing, the two mutants appeared to be allelic. When no such candidates are found, the mapping procedure should be continued. We present a strategy to efficiently select recombinants that can be used for fine mapping.Electronic Supplementary Material Supplementary material is available in the online version of this article at .Communicated by R. Hagemann     

Comparison of the power and accuracy of biallelic and microsatellite markers in population-based gene-mapping methods.

Because of their great abundance and amenability to fully automated genotyping, single-nucleotide polymorphisms (SNPs) and simple insertion/deletion are emerging as a new generation of markers for positional cloning. Although the efficiency and cost associated with the markers are important in the mapping of human disease genes, the power to detect the linkage between the marker and the disease locus, as well as the accuracy of the estimation of the map location of the disease gene, dictate the selection of the markers. Both the power and the accuracy depend not only on the type of the markers but also on other factors, such as the age of the disease mutation, the magnitude of the genetic effect, the marker-allele distribution in the population, mutation rates of marker loci, the frequency of the disease allele, the recombination fraction, and the methods for mapping the human disease genes. In this article, we develop a mathematical framework and the analytical formulas for calculation of the power and the accuracy and investigate the impact that the aforementioned factors have on the power and the accuracy, by using two population-based gene-mapping methods-likelihood-based linkage-disequilibrium mapping and the transmission/disequilibrium test, for both biallelic SNPs and microsatellites. These studies provide not only guidance in selection of the markers and in the design of the sample scheme for positional cloning but also insight into the biological bases of the mapping of human disease genes.     

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.     

水稻QTL图位克隆的特征分析

图位克隆作为基因克隆的最有效手段, 在水稻QTL克隆方面取得了很大的进展。文章通过分析近年来关于水稻QTL图位克隆的15个成功案例, 总结了水稻QTL图位克隆的几个重要特征: (1) 亲本杂交类型为种间或亚种间杂交, 双亲的目标性状差异显著; (2) 目标QTL为主效, 一般能解释大部分的表型变异; (3) 物理图距一般小于40 kb; (4)初级定位结果准确, 精细定位群体大于6 000(隐性群体不小于1 500)单株。文章还对QTL图位克隆的难点及解决方法进行了有益的探讨。     

A comprehensive SNP-based genetic analysis of inbred mouse strains

Dense genetic maps of mammalian genomes facilitate a variety of biological studies including the mapping of polygenic traits, positional cloning of monogenic traits, mapping of quantitative or qualitative trait loci, marker association, allelic imbalance, speed congenic construction, and evolutionary or phylogenetic comparison. In particular, single nucleotide polymorphisms (SNPs) have proved useful because of their abundance and compatibility with multiple high-throughput technology platforms. SNP genotyping is especially suited for the genetic analysis of model organisms such as the mouse because biallelic markers remain fully informative when used to characterize crosses between inbred strains. Here we report the mapping and genotyping of 673 SNPs (including 519 novel SNPs) in 55 of the most commonly used mouse strains. These data have allowed us to construct a phylogenetic tree that correlates and expands known genealogical relationships and clarifies the origin of strains previously having an uncertain ancestry. All 55 inbred strains are distinguishable genetically using this SNP panel. Our data reveal an uneven SNP distribution consistent with a mosaic pattern of inheritance and provide some insight into the changing dynamics of the physical architecture of the genome. Furthermore, these data represent a valuable resource for the selection of markers and the design of experiments that require the genetic distinction of any pair of mouse inbred strains such as the generation of congenic mice, positional cloning, and the mapping of quantitative or qualitative trait loci.The content of this publication does not necessarily reflect the view or policies of the Department of Health and Human Services, nor does mention of trade names, commercial products, or organizations imply endorsement by the U.S. Government.     

A PCR-Based Genetic Map for Human Chromosome 3

Oligonucleotide primers for 125 simple sequence repeat microsatellite-based genetic markers have been assayed by polymerase chain reaction (PCR) in the CEPH reference family panel. These microsatellites include 101 dinucleotide repeats as well as 24 new tetranucleotide repeats. The average heterozygosity of this marker set was 72.4%. Genetic data were analyzed with the genetic mapping package LINKAGE. A subset of these microsatellite markers define a set of 56 uniquely ordered loci (>1000:1 against local inversion) that span 271 cM. Sixty-seven additional loci were tightly linked to markers on the uniquely ordered map, but could not be ordered with such high precision. These markers were positioned by CMAP into confidence intervals. One hundred thirteen of the microsatellite markers were also tested on a chromosome 3 framework somatic cell hybrid panel that divides this chromosome into 23 cytogenetically defined regions, integrating the genetic and physical maps of this chromosome. The high density, high heterozygosity, and PCR format of this genetically and physically mapped set of markers will accelerate the mapping and positional cloning of new chromosome 3 genes.     

Limits of resolution of genetic linkage studies: implications for the positional cloning of human disease genes.

Positional cloning studies to identify disease genes are being carried out for many human genetic diseases. Such studies often include a genome-scan linkage analysis to identify the rough chromosomal location of a disease gene, fine structure genetic mapping to define and narrow the chromosomal interval in which the disease gene may be located, and physical mapping and gene identification in the genetically defined interval to clone the disease gene. During the planning of a positional cloning study, it is important to know that, if linkage is found, the genetic interval identified is likely to be sufficiently narrow to be dissected efficiently by methods of physical mapping and gene identification. Thus, we wish to know the limits of resolution of a genetic linkage study. In this paper, I determine for Mendelian diseases the distributions and moments of three measures of linkage resolution: (1) in a set of N chromosomes, the distance between the nearest crossovers that flank a disease locus, (2) the distance between the nearest genetic markers that flank the pair of flanking crossovers after a genome scan, and (3) the distance between the nearest flanking markers after additional randomly placed markers are generated and typed in an identified interval. These results provide explicit sample-size guidelines for future positional cloning studies of Mendelian diseases and make possible a more objective evaluation of whether a proposed positional cloning study is likely to be successful. I also briefly discuss the more difficult problem of linkage resolution for complex genetic diseases.     

Power calculations for selective genotyping in QTL mapping in backcross mice

Selective genotyping concerns the genotyping of a portion of individuals chosen on the basis of their phenotypic values. Often individuals are selected for genotyping from the high and low extremes of the phenotypic distribution. This procedure yields savings in cost and time by decreasing the total number of individuals genotyped. Previous work by Darvasi et al. (1993) has shown that the power to detect a QTL by genotyping 40-50 % of a population is roughly equivalent to genotyping the entire sample. However, these power studies have not accounted for different strategies of analysing the data when phenotypes of individuals in the middle are excluded, nor have they investigated the genome-wide type I error rate under these different strategies or different selection percentages. Further, these simulation studies have not considered markers over the entire genome. In this paper, we present simulation studies of power for the maximum likelihood approach to QTL mapping by Lander & Botstein (1989) in the context of selective genotyping. We calculate the power of selectively genotyping the individuals from the middle of the phenotypic distribution when performing QTL mapping over the whole mouse genome.     

Molecular mapping of Or, a gene inducing beta-carotene accumulation in cauliflower (Brassica oleracea L. var. botrytis).

The cauliflower (Brassica oleracea L. var. botrytis) Or gene is a semi-dominant, single-locus mutation that induces the accumulation of high levels of beta-carotene in various tissues of the plant, turning them orange. As part of a map-based cloning strategy, molecular mapping of the Or gene in the cauliflower genome was undertaken in a mapping population consisting of 195 F2 individuals. By using amplified fragment length polymorphism (AFLP) in conjunction with bulked segregant analysis, we identified 10 AFLP markers closely linked to the Or gene. Four of the most closely linked flanking markers were converted into restriction fragment length polymorphism (RFLP) markers. Mapping of these markers in the mapping population placed two of them at 0.5 cM from the Or locus on one side, while another marker flanked the Or gene at 1.6 cM on the other side. Three of these markers were also successfully converted into sequence-characterized amplified region (SCAR) markers. These PCR-based markers will be useful for a large-scale application in facilitating the positional cloning of the Or gene.     

Genotype‐by‐sequencing of the plant‐pathogenic fungi Pyrenophora teres and Sphaerulina musiva utilizing Ion Torrent sequence technology

Genetic and genomics tools to characterize host–pathogen interactions are disproportionately directed to the host because of the focus on resistance. However, understanding the genetics of pathogen virulence is equally important and has been limited by the high cost of de novo genotyping of species with limited marker data. Non‐resource‐prohibitive methods that overcome the limitation of genotyping are now available through genotype‐by‐sequencing (GBS). The use of a two‐enzyme restriction‐associated DNA (RAD)‐GBS method adapted for Ion Torrent sequencing technology provided robust and reproducible high‐density genotyping of several fungal species. A total of 5783 and 2373 unique loci, ‘sequence tags’, containing 16 441 and 9992 single nucleotide polymorphisms (SNPs) were identified and characterized from natural populations of Pyrenophora teres f. maculata and Sphaerulina musiva, respectively. The data generated from the P. teres f. maculata natural population were used in association mapping analysis to map the mating‐type gene to high resolution. To further validate the methodology, a biparental population of P. teres f. teres, previously used to develop a genetic map utilizing simple sequence repeat (SSR) and amplified fragment length polymorphism (AFLP) markers, was re‐analysed using the SNP markers generated from this protocol. A robust genetic map containing 1393 SNPs on 997 sequence tags spread across 15 linkage groups with anchored reference markers was generated from the P. teres f. teres biparental population. The robust high‐density markers generated using this protocol will allow positional cloning in biparental fungal populations, association mapping of natural fungal populations and population genetics studies.     

Duplicate and diverge: the evolution of plant genome microstructure

The use of positional approaches for the isolation of genes from most crop species is difficult due to the large size of their genomes. If the order of genes in segments of the genomes is similar in different plants, it might be feasible to use smaller genomes as templates upon which to base strategies for the positional cloning of genes from other species. Comparative genetic mapping, using markers such as restriction-fragment length polymorphisms, has revealed extensive conservation of long-range genome organization (macrostructure) between related species. But is the organization of the tens or hundreds of genes between the genetic markers also conserved? Recent results suggest that the fine-scale structure (microstructure) of plant genomes is more dynamic than previously assumed from investigations of the macrostructure.     

An optimal DNA pooling strategy for progressive fine mapping

We present a cost-effective DNA pooling strategy for fine mapping of a single Mendelian gene in controlled crosses. The theoretical argument suggests that it is potentially possible for a single-stage pooling approach to reduce the overall experimental expense considerably by balancing costs for genotyping and sample collection. Further, the genotyping burden can be reduced through multi-stage pooling. Numerical results are provided for practical guidelines. For example, the genotyping effort can be reduced to only a small fraction of that needed for individual genotyping at a small loss of estimation accuracy or at a cost of increasing sample sizes slightly when recombination rates are 0.5% or less. An optimal two-stage pooling scheme can reduce the amount of genotyping to 19.5%, 14.5% and 6.4% of individual genotyping efforts for identifying a gene within 1, 0.5, and 0.1 cM, respectively. Finally, we use a genetic data set for mapping the rice xl(t) gene to demonstrate the feasibility and efficiency of the DNA pooling strategy. Taken together, the results demonstrate that this DNA pooling strategy can greatly reduce the genotyping burden and the overall cost in fine mapping experiments.     

Identifying the susceptibility gene(s) in a set of trait-linked genes using genotype data

There are generally three steps to isolate a disease linkage-susceptibility gene: genome-wide scan, fine mapping, and, last, positional cloning. The last step is time consuming and involves intensive laboratory work. In some cases, fine mapping cannot proceed further on a set of markers because they are tightly linked. For years, genetic statisticians have been trying different ways to narrow the fine-mapping results to provide some guidance for the next step of laboratory work. Although these methods are practical and efficient, most of them are based on IBD data, which usually can be inferred only from the genotype data with some uncertainty. The corresponding methods thus have no greater power than one using genotype data directly. Also, IBD-based methods apply only to relative pair data. Here, using genotype data, we have developed a statistical hypothesis-testing method to pinpoint a SNP, or SNPs, suspected of responsibility for a disease trait linkage among a set of SNPs tightly linked in a region. Our method uses genotype data of affected individuals or case-control studies, which are widely available in the laboratory. The testing statistic can be constructed using any genotype-based disease-marker disequilibrium measure and is asymptotically distributed as a chi-square mixture. This method can be used for singleton data, relative pair data, or general pedigree data. We have applied the method to simulated data as well as a real data set; it gives satisfactory results.     

High-throughput scanning of the rat genome using interspersed repetitive sequence-PCR markers

We report the establishment of a hybridization-based marker system for the rat genome based on the PCR amplification of interspersed repetitive sequences (IRS). Overall, 351 IRS markers were mapped within the rat genome. The IRS marker panel consists of 210 nonpolymorphic and 141 polymorphic markers that were screened for presence/absence polymorphism patterns in 38 different rat strains and substrains that are commonly used in biomedical research. The IRS marker panel was demonstrated to be useful for rapid genome screening in experimental rat crosses and high-throughput characterization of large-insert genomic library clones. Information on corresponding YAC clones is made available for this IRS marker set distributed over the whole rat genome. The two existing rat radiation hybrid maps were integrated by placing the IRS markers in both maps. The genetic and physical mapping data presented provide substantial information for ongoing positional cloning projects in the rat.     

Efficient and accurate construction of genetic linkage maps from the minimum spanning tree of a graph

Genetic linkage maps are cornerstones of a wide spectrum of biotechnology applications, including map-assisted breeding, association genetics, and map-assisted gene cloning. During the past several years, the adoption of high-throughput genotyping technologies has been paralleled by a substantial increase in the density and diversity of genetic markers. New genetic mapping algorithms are needed in order to efficiently process these large datasets and accurately construct high-density genetic maps. In this paper, we introduce a novel algorithm to order markers on a genetic linkage map. Our method is based on a simple yet fundamental mathematical property that we prove under rather general assumptions. The validity of this property allows one to determine efficiently the correct order of markers by computing the minimum spanning tree of an associated graph. Our empirical studies obtained on genotyping data for three mapping populations of barley (Hordeum vulgare), as well as extensive simulations on synthetic data, show that our algorithm consistently outperforms the best available methods in the literature, particularly when the input data are noisy or incomplete. The software implementing our algorithm is available in the public domain as a web tool under the name MSTmap.     

Optimal strategies for mapping complex diseases in the presence of multiple loci.

Recent advances in genome technology have led to mapping and subsequent isolation, by positional cloning, of a number of genes for common and/or complex human diseases. It therefore will be possible to utilize information about a known locus in the search for additional, perhaps less penetrant, genes for a particular disease. It is also unclear, under these situations, what the optimal sampling strategy should be. To address these questions, we have calculated the expected LOD score for localizing one locus in a variety of two-locus models of disease, for four different pedigree structures, and under three different scenarios regarding knowledge/testing of one of the two loci. These design considerations are evaluated by use of a cost function that incorporates the costs of ascertaining different family structures, the relative costs of genotyping and mutation testing family members, and the amount of information provided by each family structure and testing scenario. The results indicate that, in most cases, affected sib pairs are a particularly poor strategy, especially when linkage or mutation data are available at the known locus. We also demonstrate that prescreening the sample of families for mutations at known susceptibility loci is, in general, a cost-effective strategy.     

Molecular map of the Chlamydomonas reinhardtii nuclear genome

We have prepared a molecular map of the Chlamydomonas reinhardtii genome anchored to the genetic map. The map consists of 264 markers, including sequence-tagged sites (STS), scored by use of PCR and agarose gel electrophoresis, and restriction fragment length polymorphism markers, scored by use of Southern blot hybridization. All molecular markers tested map to one of the 17 known linkage groups of C. reinhardtii. The map covers approximately 1,000 centimorgans (cM). Any position on the C. reinhardtii genetic map is, on average, within 2 cM of a mapped molecular marker. This molecular map, in combination with the ongoing mapping of bacterial artificial chromosome (BAC) clones and the forthcoming sequence of the C. reinhardtii nuclear genome, should greatly facilitate isolation of genes of interest by using positional cloning methods. In addition, the presence of easily assayed STS markers on each arm of each linkage group should be very useful in mapping new mutations in preparation for positional cloning.     

Positional cloning by fast-track SNP-mapping in Drosophila melanogaster

Positional cloning of chemically induced mutations is the rate-limiting step in forward genetic screens in Drosophila. Single-nucleotide polymorphisms (SNPs) are useful markers to locate a mutated region in the genome. Here, we provide a protocol for high-throughput, high-resolution SNP mapping that enables rapid and cost-effective positional cloning in Drosophila. In stage 1 of the protocol, we use highly multiplexed tag-array mini-sequencing assays to map mutations to an interval of 1-2 Mb. In these assays, SNPs are genotyped by primer extension using fluorescently labeled dideoxy-nucleotides. Fluorescent primers are captured and detected on a microarray. In stage 2, we selectively isolate recombinants within the identified 1-2 Mb interval for fine mapping of mutations to about 50 kb. We have previously demonstrated the applicability of this protocol by mapping 14 muscle morphogenesis mutants within 4 months, which represents a significant acceleration compared with other commonly used mapping strategies that may take years.