High-resolution genotyping via whole genome hybridizations to microarrays containing long oligonucleotide probes

To date, microarray-based genotyping of large, complex plant genomes has been complicated by the need to perform genome complexity reduction to obtain sufficiently strong hybridization signals. Genome complexity reduction techniques are, however, tedious and can introduce unwanted variables into genotyping assays. Here, we report a microarray-based genotyping technology for complex genomes (such as the 2.3 GB maize genome) that does not require genome complexity reduction prior to hybridization. Approximately 200,000 long oligonucleotide probes were identified as being polymorphic between the inbred parents of a mapping population and used to genotype two recombinant inbred lines. While multiple hybridization replicates provided ～97% accuracy, even a single replicate provided ～95% accuracy. Genotyping accuracy was further increased to >99% by utilizing information from adjacent probes. This microarray-based method provides a simple, high-density genotyping approach for large, complex genomes.     

Genotyping‐by‐sequencing approaches to characterize crop genomes: choosing the right tool for the right application

In the last decade, the revolution in sequencing technologies has deeply impacted crop genotyping practice. New methods allowing rapid, high‐throughput genotyping of entire crop populations have proliferated and opened the door to wider use of molecular tools in plant breeding. These new genotyping‐by‐sequencing (GBS) methods include over a dozen reduced‐representation sequencing (RRS) approaches and at least four whole‐genome resequencing (WGR) approaches. The diversity of methods available, each often producing different types of data at different cost, can make selection of the best‐suited method seem a daunting task. We review the most common genotyping methods used today and compare their suitability for linkage mapping, genomewide association studies (GWAS), marker‐assisted and genomic selection and genome assembly and improvement in crops with various genome sizes and complexity. Furthermore, we give an outline of bioinformatics tools for analysis of genotyping data. WGR is well suited to genotyping biparental cross populations with complex, small‐ to moderate‐sized genomes and provides the lowest cost per marker data point. RRS approaches differ in their suitability for various tasks, but demonstrate similar costs per marker data point. These approaches are generally better suited for de novo applications and more cost‐effective when genotyping populations with large genomes or high heterozygosity. We expect that although RRS approaches will remain the most cost‐effective for some time, WGR will become more widespread for crop genotyping as sequencing costs continue to decrease.     

Development,validation and genetic analysis of a large soybean SNP genotyping array

Cultivated soybean (Glycine max) suffers from a narrow germplasm relative to other crop species, probably because of under‐use of wild soybean (Glycine soja) as a breeding resource. Use of a single nucleotide polymorphism (SNP) genotyping array is a promising method for dissecting cultivated and wild germplasms to identify important adaptive genes through high‐density genetic mapping and genome‐wide association studies. Here we describe a large soybean SNP array for use in diversity analyses, linkage mapping and genome‐wide association analyses. More than four million high‐quality SNPs identified from high‐depth genome re‐sequencing of 16 soybean accessions and low‐depth genome re‐sequencing of 31 soybean accessions were used to select 180 961 SNPs for creation of the Axiom^® SoyaSNP array. Validation analysis for a set of 222 diverse soybean lines showed that 170 223 markers were of good quality for genotyping. Phylogenetic and allele frequency analyses of the validation set data indicated that accessions showing an intermediate morphology between cultivated and wild soybeans collected in Korea were natural hybrids. More than 90 unanchored scaffolds in the current soybean reference sequence were assigned to chromosomes using this array. Finally, dense average spacing and preferential distribution of the SNPs in gene‐rich chromosomal regions suggest that this array may be suitable for genome‐wide association studies of soybean germplasm. Taken together, these results suggest that use of this array may be a powerful method for soybean genetic analyses relating to many aspects of soybean breeding.     

Comparative physical mapping of targeted regions of the rat genome

The comparative mapping and sequencing of vertebrate genomes is now a key priority for the Human Genome Project. In addition to finishing the human genome sequence and generating a `working draft' of the mouse genome sequence, significant attention is rapidly turning to the analysis of other model organisms, such as the laboratory rat (Rattus norvegicus). As a complement to genome-wide mapping and sequencing efforts, it is often important to generate detailed maps and sequence data for specific regions of interest. Using an adaptation of our previously described approach for constructing mouse comparative and physical maps, we have established a general strategy for targeted mapping of the rat genome. Specifically, we constructed a framework comparative map of human Chromosome (Chr) 7 and the orthologous regions of the rat genome, as well as two large (>1-Mb) P1-derived artificial chromosome (PAC)-based physical maps. Generation of these physical maps involved the use of mouse-derived probes that cross-hybridized with rat PAC clones. The first PAC map encompasses the cystic fibrosis transmembrane conductance regulator gene (Cftr), while the second map allows a three-species comparison of a genomic region containing intra- and inter-chromosomal evolutionary rearrangements. The studies reported here further demonstrate that cross-species hybridization between related animals, such as rat and mouse, can be readily used for the targeted construction of clone-based physical maps, thereby accelerating the analysis of biologically interesting regions of vertebrate genomes. Received: 5 December 2000 / Accepted: 27 February 2001     

Inter- and intra-genomic homology of the Brassica genomes: implications for their origin and evolution

In order to determine the homologous regions shared by the cultivated Brassica genomes, linkage maps of the diploid cultivated B. rapa (A genome, n = 10), B. nigra (B genome, n = 8) and B. oleracea (C genome, n = 9), were compared. We found intergenomic conserved regions but with extensitve reordering among the genomes. Eighteen linkage groups from all three species could be associated on the basis of homologous segments based on at least three common markers. Intragenomic homologous conservation was also observed for some of the chromosomes of the A, B and C genomes. A possible chromosome phylogenetic pathway based on an ancestral genome of at least five, and no more than seven chromosomes, was drawn from the chromosomal inter-relationships observed. These results demonstrate that extensive duplication and rearrangement have been involved in the formation of the Brassica genomes from a smaller ancestral genome.     

Ultrahigh-Density Linkage Map for Cultivated Cucumber (Cucumis sativus L.) Using a Single-Nucleotide Polymorphism Genotyping Array

Genotyping arrays are tools for high-throughput genotyping, which is beneficial in constructing saturated genetic maps and therefore high-resolution mapping of complex traits. Since the report of the first cucumber genome draft, genetic maps have been constructed mainly based on simple-sequence repeats (SSRs) or on combinations of SSRs and sequence-related amplified polymorphism (SRAP). In this study, we developed the first cucumber genotyping array consisting of 32,864 single-nucleotide polymorphisms (SNPs). These markers cover the cucumber genome with a median interval of ~2 Kb and have expected genotype calls in parents/F₁ hybridizations as a training set. The training set was validated with Fluidigm technology and showed 96% concordance with the genotype calls in the parents/F₁ hybridizations. Application of the genotyping array was illustrated by constructing a 598.7 cM genetic map based on a ‘9930’ × ‘Gy14’ recombinant inbred line (RIL) population comprised of 11,156 SNPs. Marker collinearity between the genetic map and reference genomes of the two parents was estimated at R² = 0.97. We also used the array-derived genetic map to investigate chromosomal rearrangements, regional recombination rate, and specific regions with segregation distortions. Finally, 82% of the linkage-map bins were polymorphic in other cucumber variants, suggesting that the array can be applied for genotyping in other lines. The genotyping array presented here, together with the genotype calls of the parents/F₁ hybridizations as a training set, should be a powerful tool in future studies with high-throughput cucumber genotyping. An ultrahigh-density linkage map constructed by this genotyping array on RIL population may be invaluable for assembly improvement, and for mapping important cucumber QTLs.     

Integrating cytogenetics and genomics in comparative evolutionary studies of cichlid fish

ABSTRACT: BACKGROUND: The availability of a large number of recently sequenced vertebrate genomes opens new avenues to integrate cytogenetics and genomics in comparative and evolutionary studies. Cytogenetic mapping can offer alternative means to identify conserved synteny shared by distinct genomes and also to define genome regions that are still not fine characterized even after wide-ranging nucleotide sequence efforts. An efficient way to perform comparative cytogenetic mapping is based on BAC clones mapping by fluorescence in situ hybridization. In this report, to address the knowledge gap on the genome evolution in cichlid fishes, BAC clones of an Oreochromis niloticus library covering the linkage groups (LG) 1, 3, 5, and 7 were mapped onto the chromosomes of 9 African cichlid species. The cytogenetic mapping data were also integrated with BAC-end sequences information of O. niloticus and comparatively analyzed against the genome of other fish species and vertebrates. RESULTS: The location of BACs from LG1, 3, 5, and 7 revealed a strong chromosomal conservation among the analyzed cichlid species genomes, which evidenced a synteny of the markers of each LG. Comparative in silico analysis also identified large genomic blocks that were conserved in distantly related fish groups and also in other vertebrates. CONCLUSIONS: Although it has been suggested that fishes contain plastic genomes with high rates of chromosomal rearrangements and probably low rates of synteny conservation, our results evidence that large syntenic chromosome segments have been maintained conserved during evolution, at least for the considered markers. Additionally, our current cytogenetic mapping efforts integrated with genomic approaches conduct to a new perspective to address important questions involving chromosome evolution in fishes.     

Identification of repeat structure in large genomes using repeat probability clouds

The identification of repeat structure in eukaryotic genomes can be time-consuming and difficult because of the large amount of information (3 × 10⁹ bp) that needs to be processed and compared. We introduce a new approach based on exact word counts to evaluate, de novo, the repeat structure present within large eukaryotic genomes. This approach avoids sequence alignment and similarity search, two of the most time-consuming components of traditional methods for repeat identification. Algorithms were implemented to efficiently calculate exact counts for any length oligonucleotide in large genomes. Based on these oligonucleotide counts, oligonucleotide excess probability clouds, or “P-clouds,” were constructed. P-clouds are composed of clusters of related oligonucleotides that occur, as a group, more often than expected by chance. After construction, P-clouds were mapped back onto the genome, and regions of high P-cloud density were identified as repetitive regions based on a sliding window approach. This efficient method is capable of analyzing the repeat content of the entire human genome on a single desktop computer in less than half a day, at least 10-fold faster than current approaches. The predicted repetitive regions strongly overlap with known repeat elements as well as other repetitive regions such as gene families, pseudogenes, and segmental duplicons. This method should be extremely useful as a tool for use in de novo identification of repeat structure in large newly sequenced genomes.     

Mapping in the realm of polyploidy: The wheat model

Wheat is an allopolyploid containing three distinct but genetically related (homoeologous) genomes, A, B and D. Because of polyploid inheritance and large genome size (16×10¹² bp), the wheat genome is thought to be intractable to map-based cloning of agronomic and other genes of interest. We propose a targeted geneti mapping strategy that combines linkage and physical mapping and may facilitate map-based cloning. High-density linkage maps are either generated in wheat or in diploid Triticum tauschii, the donor of the D genome to wheat. Molecular marker-based chromosome maps are constructed, using an array of deletion lines in wheat. The conventional genetic linkage maps are aligned with chromosome maps to construct cytogenetic ladder maps (CLMs). The CLMs allow region-specific mapping and convert genetic distances into physical distances. The information from CLMs suggests that many genes in wheat are present in clusters that are highly recombiogenic, small, and may be amenable to cloning by chromosome walking. Therefore, the effective genome size of wheat is relatively small in comparison to the whole genome. The utility of using the smaller genome of rice for mapping and homologous gene cloning is discussed.     

Microsatellite markers for genome analysis in Brassica. II. Assignment of rapeseed microsatellites to the A and C genomes and genetic mapping in Brassica oleracea L.

Microsatellites are highly polymorphic and efficient markers for the analysis of plant genomes. Primer specificity, however, may restrict the applicability of these markers even between closely related species for comparative mapping studies. We have demonstrated that the majority of microsatellites identified in oilseed rape (Brassica napus L; AC genome) correspond to loci which can be easily assigned to the A and C progenitor genomes. A study with 63 primer pairs has shown that 54% detect two loci, one from each genome, while 25% and 21%, respectively, are either A or C genome-specific. The distribution of rapeseed microsatellites in the C genome was investigated by genetic mapping in Brassica oleracea L. Ninety two dinucleotide microsatellites were screened for polymorphism in an F₂ population derived from a cross between collard and cauliflower, for which an RFLP map has been constructed previously. Thirty three primer pairs (35.7%) have yielded either unspecific or no PCR products whereas the remaining primer pairs amplified one or more distinct loci. The level of polymorphism found in the mapping population was 49.2%. A total of 29 primer pairs disclosed 34 loci of which 31 are evenly distributed on 8 of the 9 B. oleracea linkage groups. For the remaining three markers linkage could not be established. Our results showed that microsatellite markers from the composite genome of B. napus can serve as a useful marker system in genetic studies and for plant-breeding objectives in B. oleracea. Received: 14 April 2000 / Accepted: 3 July 2000     

An extended anchored linkage map and virtual mapping for the American mink genome based on homology to human and dog

In this report we present an extended linkage map of the American mink (Neovison vison) consisting of 157 microsatellite markers and comprising at least one linkage group for each of the autosomes. Each linkage group has been assigned to a chromosome and oriented by fluorescence in situ hybridization (FISH) and/or by means of human/dog/mink comparative homology. The average interval between markers is 8.5 cM and the linkage groups collectively span 1340 cM. In addition, 217 and 275 mink microsatellites have been placed on human and dog genomes, respectively. In conjunction with the existing comparative human/dog/mink data, these assignments represent useful virtual maps for the American mink genome. Comparison of the current human/dog assembled sequential map with the existing Zoo-FISH-based human/dog/mink maps helped to refine the human/dog/mink comparative map. Furthermore, comparison of the human and dog genome assemblies revealed a number of large synteny blocks, some of which are corroborated by data from the mink linkage map.     

Reference-free SNP calling: improved accuracy by preventing incorrect calls from repetitive genomic regions

ABSTRACT: BACKGROUND: Single nucleotide polymorphisms (SNPs) are the most abundant type of genetic variation in eukaryotic genomes and have recently become the marker of choice in a wide variety of ecological and evolutionary studies. The advent of next-generation sequencing (NGS) technologies has made it possible to efficiently genotype a large number of SNPs in the non-model organisms with no or limited genomic resources. Most NGS-based genotyping methods require a reference genome to perform accurate SNP calling. Little effort, however, has yet been devoted to developing or improving algorithms for accurate SNP calling in the absence of a reference genome. RESULTS: Here we describe an improved maximum likelihood (ML) algorithm called iML, which can achieve high genotyping accuracy for SNP calling in the non-model organisms without a reference genome. The iML algorithm incorporates the mixed Poisson/normal model to detect composite read clusters and can efficiently prevent incorrect SNP calls resulting from repetitive genomic regions. Through analysis of simulation and real sequencing datasets, we demonstrate that in comparison with ML or a threshold approach, iML can remarkably improve the accuracy of de novo SNP genotyping and is especially powerful for the reference-free genotyping in diploid genomes with high repeat contents. CONCLUSIONS: The iML algorithm can efficiently prevent incorrect SNP calls resulting from repetitive genomic regions, and thus outperforms the original ML algorithm by achieving much higher genotyping accuracy. Our algorithm is therefore very useful for accurate de novo SNP genotyping in the non-model organisms without a reference genome.     

Genomic Characterization of DArT Markers Based on High-Density Linkage Analysis and Physical Mapping to the Eucalyptus Genome

Diversity Arrays Technology (DArT) provides a robust, high throughput, cost-effective method to query thousands of sequence polymorphisms in a single assay. Despite the extensive use of this genotyping platform for numerous plant species, little is known regarding the sequence attributes and genome-wide distribution of DArT markers. We investigated the genomic properties of the 7,680 DArT marker probes of a Eucalyptus array, by sequencing them, constructing a high density linkage map and carrying out detailed physical mapping analyses to the Eucalyptus grandis reference genome. A consensus linkage map with 2,274 DArT markers anchored to 210 microsatellites and a framework map, with improved support for ordering, displayed extensive collinearity with the genome sequence. Only 1.4 Mbp of the 75 Mbp of still unplaced scaffold sequence was captured by 45 linkage mapped but physically unaligned markers to the 11 main Eucalyptus pseudochromosomes, providing compelling evidence for the quality and completeness of the current Eucalyptus genome assembly. A highly significant correspondence was found between the locations of DArT markers and predicted gene models, while most of the 89 DArT probes unaligned to the genome correspond to sequences likely absent in E. grandis, consistent with the pan-genomic feature of this multi-Eucalyptus species DArT array. These comprehensive linkage-to-physical mapping analyses provide novel data regarding the genomic attributes of DArT markers in plant genomes in general and for Eucalyptus in particular. DArT markers preferentially target the gene space and display a largely homogeneous distribution across the genome, thereby providing superb coverage for mapping and genome-wide applications in breeding and diversity studies. Data reported on these ubiquitous properties of DArT markers will be particularly valuable to researchers working on less-studied crop species who already count on DArT genotyping arrays but for which no reference genome is yet available to allow such detailed characterization.     

Molecular genetics of growth and development in Populus. III. A genetic linkage map of a hybrid poplar composed of RFLP,STS, and RAPD markers

We have evaluated three DNA-based marker types for linkage map construction in Populus: RFLPs detected by Southern blot hybridization, STSs detected by a combination of PCR and RFLP analysis, and RAPDs. The mapping pedigree consists of three generations, with the F₁ produced by interspecific hybridization between a P. trichocarpa female and a P. deltoides male. The F₂ generation was made by inbreeding to the maximum degree permitted by the dioecious mating system of Populus. The applicability of STSs and RAPDs outside the mapping pedigree has been investigated, showing that these PCR-based marker systems are well-suited to breeding designs involving interspecific hybridization. A Populus genome map (343 markers) has been constructed from a combination of all three types. The length of the Populus genome is estimated to be 2400–2800 cM.Abbreviations RFLP restriction fragment length polymorphism - STS sequence-tagged site - PCR polymerase chain reaction - RAPD random amplified polymorphic DNA     

Development and initial characterization of a HAPPY panel for mapping the X. tropicalis genome

HAPPY mapping was designed to pursue the analysis of approximately random HAPloid DNA breakage samples using the PolYmerase chain reaction for mapping genomes. In the present study, we improved the method and integrated two other molecular techniques into the process: whole genome amplification and the Sequenom SNP (single nucleotide polymorphism) genotyping assay in order to facilitate whole genome mapping of X. tropicalis. The former technique amplified enough DNA materials to genotype a large number of markers, while the latter allowed for relatively high throughput marker genotyping with multiplex assays on the HAPPY lines. A total of 58 X. tropicalis genes were genotyped on an initial panel of 383 HAPPY lines, which contributed to formation of a working panel of 146 lines. Further genotyping of 29 markers on the working panel led to construction of a HAPPY map for the X. tropicalis genome. We believe that our improved HAPPY method described in the present study has paved the way for the community to map different genomes with a simple, but powerful approach.     

AFLP and SSR polymorphism in a <Emphasis Type="Italic">Coffea</Emphasis> interspecific backcross progeny [(<Emphasis Type="Italic">C. heterocalyx</Emphasis> × <Emphasis Type="Italic">C. canephora</Emphasis>) × <Emphasis Type="Italic">C. canephora</Emphasis>]

An interspecific cross (BC 1) involving a species with one of the largest genomes in the Coffea genus [Coffea heterocalyx (HET), qDNA = 1.74 pg] and a species with a medium-sized genome [Coffea canephora (CAN), qDNA = 1.43 pg] was studied using two types of molecular markers, AFLP and SSR. One hundred and eighty eight AFLP bands and 34 SSR primer pairs were suitable for mapping. The total map length was 1,360 cM with 190 loci distributed in 15 linkage groups. The results were compared to those obtained previously on an interspecific BC 1 progeny involving a species with a medium-sized genome (Coffea liberica var dewevrei, DEW) and a species with one of the smallest genomes (Coffea pseudozanguebariae, PSE). They are discussed relative to three main points: (1) the relevance of the different marker types, (2) the genomic distribution of AFLP and SSR markers, and (3) the relation between AFLP polymorphism and genome size.Communicated by H.F. Linskens     

MIG-seq is an effective method for high-throughput genotyping in wheat (Triticum spp.)

MIG-seq (Multiplexed inter-simple sequence repeats genotyping by sequencing) has been developed as a low cost genotyping technology, although the number of polymorphisms obtained is assumed to be minimal, resulting in the low application of this technique to analyses of agricultural plants. We applied MIG-seq to 12 plant species that include various crops and investigated the relationship between genome size and the number of bases that can be stably sequenced. The genome size and the number of loci, which can be sequenced by MIG-seq, are positively correlated. This is due to the linkage between genome size and the number of simple sequence repeats (SSRs) through the genome. The applicability of MIG-seq to population structure analysis, linkage mapping, and quantitative trait loci (QTL) analysis in wheat, which has a relatively large genome, was further evaluated. The results of population structure analysis for tetraploid wheat showed the differences among collection sites and subspecies, which agreed with previous findings. Additionally, in wheat biparental mapping populations, over 3,000 SNPs/indels with low deficiency were detected using MIG-seq, and the QTL analysis was able to detect recognized flowering-related genes. These results revealed the effectiveness of MIG-seq for genomic analysis of agricultural plants with large genomes, including wheat.     

An integrated molecular linkage map of diploid wheat based on a <Emphasis Type="Italic">Triticum boeoticum</Emphasis> × <Emphasis Type="Italic">T. monococcum</Emphasis> RIL population

Diploid A genome species of wheat harbour immense variability for biotic stresses and productivity traits, and these could be transferred efficiently to hexaploid wheat through marker assisted selection, provided the target genes are tagged at diploid level first. Here we report an integrated molecular linkage map of A genome diploid wheat based on 93 recombinant inbred lines (RILs) derived from Triticum boeoticum × Triticum monococcum inter sub-specific cross. The parental lines were analysed with 306 simple sequence repeat (SSR) and 194 RFLP markers, including 66 bin mapped ESTs. Out of 306 SSRs tested for polymorphism, 74 (24.2%) did not show amplification (null) in both the parents. Overall, 171 (73.7%) of the 232 remaining SSR and 98 (50.5%) of the 194 RFLP markers were polymorphic. Both A and D genome specific SSR markers showed similar transferability to A genome of diploid wheat species. The 176 polymorphic markers, that were assayed on a set of 93 RILs, yielded 188 polymorphic loci and 177 of these as well as two additional morphological traits mapped on seven linkage groups with a total map length of 1,262 cM, which is longer than most of the available A genome linkage maps in diploid and hexaploid wheat. About 58 loci showed distorted segregation with majority of these mapping on chromosome 2A^m. With a few exceptions, the position and order of the markers was similar to the ones in other maps of the wheat A genome. Chromosome 1A^m of T. monococcum and T. boeoticum showed a small paracentric inversion relative to the A genome of hexaploid wheat. The described linkage map could be useful for gene tagging, marker assisted gene introgression from diploid into hexaploid wheat as well as for map based cloning of genes from diploid A genome species and orthologous genes from hexaploid wheat.     

Mapping of genome-wide resistance gene analogs (RGAs) in maize (<Emphasis Type="Italic">Zea mays</Emphasis> L.)

Isolation and mapping of genome-wide resistance (R) gene analogs (RGAs) is of importance in identifying candidate(s) for a particular resistance gene/QTL. Here we reported our result in mapping totally 228 genome-wide RGAs in maize. By developing RGA-tagged markers and subsequent genotyping a population consisting of 294 recombinant inbred lines (RILs), 67 RGAs were genetically mapped on maize genome. Meanwhile, in silico mapping was conducted to anchor 113 RGAs by comparing all 228 RGAs to those anchored EST and BAC/BAC-end sequences via tblastx search (E-value < 10⁻²⁰). All RGAs from different mapping efforts were integrated into the existing SSR linkage map. After accounting for redundancy, the resultant RGA linkage map was composed of 153 RGAs that were mapped onto 172 loci on maize genome, and the mapped RGAs accounted for approximate three quarters of the genome-wide RGAs in maize. The extensive co-localizations were observed between mapped RGAs and resistance gene/QTL loci, implying the usefulness of this RGA linkage map in R gene cloning via candidate gene approach. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Wenkai Xiao, Jing Zhao and Shengci Fan have contributed equally to this research.