GenoTypeMapper: graphical genotyping on genetic and sequence-based maps
Plant Methods volume 16, Article number: 123 (2020)
The rising availability of assemblies of large genomes (e.g. bread and durum wheat, barley) and their annotations deliver the basis to graphically present genome organization of parents and progenies on a physical scale. Genetic maps are a very important tool for breeders but often represent distorted models of the actual chromosomes, e.g., in centromeric and telomeric regions. This biased picture might lead to imprecise assumptions and estimations about the size and complexity of genetic regions and the selection of suitable molecular markers for the incorporation of traits in breeding populations or near-isogenic lines (NILs). Some software packages allow the graphical illustration of genotypic data, but to the best of our knowledge, suitable software packages that allow the comparison of genotypic data on the physical and genetic scale are currently unavailable.
We developed a simple Java-based-software called GenoTypeMapper (GTM) for comparing genotypic data on genetic and physical maps and tested it for effectiveness on data of two NILs that carry QTL-regions for drought stress tolerance from wild emmer on chromosome 2BS and 7AS. Both NILs were more tolerant to drought stress than their recurrent parents but exhibited additional undesirable traits such as delayed heading time.
In this article, we illustrate that the software easily allows users to display and identify additional chromosomal introgressions in both NILs originating from the wild emmer parent. The ability to detect and diminish linkage drag can be of particular interest for pre-breeding purposes and the developed software is a well-suited tool in this respect. The software is based on a simple allele-matching algorithm between the offspring and parents of a crossing scheme. Despite this simple approach, GTM seems to be the only software that allows us to analyse, illustrate and compare genotypic data of offspring of different crossing schemes with up to four parents in two different maps. So far, up to 500 individuals with a maximum number of 50,000 markers can be examined with the software. The main limitation that hampers the performance of the software is the number of markers that are examined in parallel. Since each individual must be analysed separately, a maximum of ten individuals can currently be displayed in a single run. On a computer with an Intel five processor of the 8th generation, GTM can reliably either analyse a single individual with up to 12,000 markers or ten individuals with up to 3,600 markers in less than five seconds. Future work aims to improve the performance of the software so that more complex crossing schemes with more parents and more markers can be analysed.
High-throughput plant genotyping technologies enable the testing of a few thousands to hundreds of thousands of markers in large numbers of samples in parallel. Such technologies include next-generation sequencing (NGS) and array-based technologies, such as the Diversity Array Technology (DArT) , genotyping by sequencing (GBS) , the Illumina Bead Array/ Infinium Technology [3, 4] or the Affymetrix Gene Chip/Axiom Technology [5, 6]. Predesigned Illumina Bead Arrays and Affymetrix chips do exist for numerous important organisms and crops, such as apples, barley, brassica, cherry, chickpea, cotton, cowpea, grape, lettuce, maize, oat, peach, pear, peanut, pepper, potato, rice, rose, rye, ryegrass, soybean, strawberry, sunflower, tomato and wheat .
Notably, only a fraction of this information can be included in a genetic map due to several limitations. The resolution of a genetic map depends on the number of progenies analysed, the number of recombination events observed between polymorphic markers [8, 9] and the number of markers along the chromosomes . Markers that do not show recombination events are assigned to the same genetic position; hence, the distances of genetic and physical maps may differ significantly from each other, especially close to centromeric and telomeric regions [11,12,13]. Genetic maps often serve as templates for the localization of major genes and quantitative traits loci (QTLs) on chromosomes. The identification of such regions depends on differences in phenotypic data that can be assigned to recombination events observed between different markers . Probability of meiosis and recombination events correlate with the number of individuals in a population. Therefore, the number of markers that are used to construct a genetic map can be adjusted to the size of the mapping population and the number of expected recombination events . The biased picture of a genetic map might lead to imprecise assumptions and estimations about the size and complexity of genetic regions and the selection of suitable molecular markers for the incorporation of traits in breeding populations or near-isogenic lines (NILs). Furthermore, genotypic data that cannot be used for the construction of a genetic map, e.g., monomorphic, heterozygous or failed data in one of the parents of a biparental mapping population, are traditionally discarded at a very early stage of the mapping procedure.
The recently published physical pseudomolecules of cereals, such as einkorn (Triticum monoccum) (4.94 Gb.) , wild emmer wheat (T. turgidum ssp. dicoccoides) (12 Gb) , durum wheat (T. turgidum ssp. durum)  goat grass (Aegilops tauschii) (4.3 Gb.) , barley (Hordeum vulgare) (4.75 Gb.)  and bread wheat (T. aestivum) (15.4–15.8 Gb) [20, 21] deliver a valuable tool to solve this distorted picture and to include missing genotypic information to finally gain a better understanding of the physical scale of the regions of interest. Software packages such as Graphical Genotypes (GGT)  and Flapjack  or complex genome browsers such as the Integrative genome viewer (IGV)  allow the graphical illustration of genotypic data, but do not allow the comparison of genotypic data of offspring of different crossing schemes in different maps. The same is true for software packages that focus on the detection and illustration of specific genomic introgressions on the genetic or physical scale, like IView  or the Physical Introgression Browser . Tools such as the R-based MareyMap  allow the comparison of genetic and physical maps but do not allow analysis of genotypic data. These consequently cannot consider the distortions between physical and genetic maps and therefore miss the information that can be obtained by integrating additional genotypic information from equivalent physical regions into the genetic map, i.e., heterozygous markers. We therefore developed a combined approach that allows graphical genotyping on the physical and genetic scale and implemented it in a small platform-independent Java application, named GenoTypeMapper (GTM). The software was developed and tested under Windows 10 but should run on any operating system (OS) that has the Java virtual machine (JVM) installed.
GTM allows users to filter different allelic types (e.g., heterozygous alleles) to plot genotypic data in different ways to obtain a customized view of genomic and/ or genetic regions and their allelic composition.
Here, we present the software in detail followed by results obtained on SNP genotyping information of two NILs (NIL-U-2B-1, NIL-B-7A-2) that carry QTL-regions for improved drought stress tolerance on chromosome 2BS and 7AS from the drought tolerant wild emmer accession # G18-16 . Both NILs were shown to be more tolerant to drought stress than their recurrent parents but show phenotypes that may stem from linkage drag, e.g., an increased plant height or an increased number of days to heading (DPH) [28, 29]. We used GTM to perform genotypic analysis on the genetic and physical scales of both NILs and the respective recombinant inbred lines (RILs) (see Additional file 1: Figure S1) [30,31,32].
GTM was developed and tested under Windows 10 but should run on any OS that has the JVM installed. A simple JAR file of the program, a setup file for windows, test data (template files) and a small quick-start tutorial are provided on www.genotypemapper.org. The test data includes the original genotypic data from the F2 plants, NILs and F2 pig individuals that are discussed in the article. Additionally, extra test data sets are provided on the website to test further options for GTM, e.g., the comparison of two genetic maps. More details about the data are also summarized in a small ReadMe file on the website.
Input data can be loaded in the form of a GTM-specific tab delimited text file with at least nine columns (Table 1).
The first two columns of the GTM input file contain chromosome and marker names, and the third and fourth column contain physical and/or genetic positions of the markers.
Column five can contain logarithm of odds (LOD) information, and columns six and seven are dedicated to genotype information from the parents of a typical bi-parental population, e.g., a Double Haploid (DH) or F2 population (Fig. 1). Recurrent parents and individuals, which are used for the establishment of NILs (Fig. 1), can be specified in the eighth and ninth columns. Genotypic data of the remaining accessions can be added to the tenth and following columns. Please note that example input files are available on www.genotypemapper.org.
User interface (UI)
After loading input data via the “File” menu, a small dialogue pops up. This dialogue allows the user to select the genotypes and chromosomes of interest. In default mode, a maximum of ten individuals and an arbitrary number of chromosomes can be selected. The alternative single chromosome mode is activated by selecting only a single genotype and chromosome. Subsequently, GTM parameters can be adapted in the “Analyse” menu via the menu items “Show maps and components”, “Parameters” and “Graphical genotyping”. The first menu item allows the user to select the information that shall be displayed, such as the genetic- and or physical map. Please note that some additional components, such as marker- and loci names, can only be displayed in the chromosome mode. “Parameters” allows the adaptation of the illustrations by modifying variables such as the zooming factor or the spacing between the physical and genetic map. Illustrations can be enlarged up to ten or 30 times in the default mode or chromosome mode, respectively. In addition to this, the physical size of the same genome can be adapted relative to the genetic map of the chromosome by dividing the physical genome by a user defined factor. Finally, the last menu item, “Graphical Genotyping”, allows the analysis and interpretation of the data that are specified from the sixth to the last column of the input file to determine the origin of alleles in the respective individuals of interest (IOIs). Working scenarios in the respective modes are illustrated in Figs. 2 and 3. Basic information about the number of genotypes that will be subjected to the analysis and the number of available markers and chromosomes are summarized under the menu item “information” (Fig. 2). A small tutorial that guides the user through this process is implemented in the software in the “HELP” menu or available on the website www.genotypemapper.org.
Allelic discrimination analysis
For each marker, the genotype of the IOI and the information available for parents is compared to determine the origin of its allele. Depending on such a comparison, five, six, ten or eleven different cases might be determined by using genotypic data of either one, two, three or four different parents (Fig. 1). These cases are explored in more detail in Table 2.
If polymorphic genotype information from one or two parents of a biparental population is provided, the IOI alleles origin can either be dedicated to parent P1 (case 1_1) or parent P2 (case 1_2) (Fig. 1, Table 2). Genotype information of three or four parents, e.g., of a triple cross, allows the identification of IOI alleles that were derived from the recurrent parent P3 (case 3_7) or from parents P2 /P3 (3_9) or P1/ P3 (3_10) (Fig. 1, Table 2). In addition to this, heterozygosity can be detected in parent P1 (case_ 1_5), parent P2 (case_2_6), parent P3 (case_3_8), parent P4 (case 4_11) (Table 2). Failed or monomorphic data (case 1_3) are filtered out by default to omit mistakes in the allelic discrimination analysis (Table 2).
The colour of each case can be defined by the user to highlight either marker positions and/or marker intervals of adjacent markers that share the same allele information. No colour is added to intervals of adjacent markers with different allele information to omit imprecise representation of the genotypic data (Additional file 2: Figure S2). Markers that belong to a specific allele type of interest can be filtered out and displayed. The image or the markers can be shown on the screen and exported as a PNG image or as a tab delimited text file. The user might change the size and orientation of the image via the menu item “Layout” (Fig. 2).
Case study—data preparation
All available marker sequences of the 15K iSelect chip  were aligned to the wild emmer genome with a local version of the Basic Local Alignment Search Tool (BLAST) v2.10 . Nucleotide BLAST (BLASTn) search was used with the following default parameters:
Expectation threshold = 10
Wordsize = 28
Match/Mismatch score = 1, −2
Gapcosts = linear
Filtering of low complexity regions = Yes
Markers with BLASTn-hits on different chromosomes were subsequently compared with the durum consensus map  and discarded if they revealed multiple BLASTn hits on the consensus map chromosomes in a distance of more than one 1,000,000 base pairs (bp). A flow chart of this work scheme was added to the attachments (Additional file 3: Figure S3). The 15K iSelect data of the corresponding NILs (NIL-U-2B-1, NIL-B-7A-2) and their recurrent parents (BarNir, Uzan) and crossing parents (F7RIL12, F7RIL55 G18-16 and LDN) were used (TraitGenetics GmbH, Gatersleben, Germany). Genotypic data that did not fail or that were not monomorphic in these eight samples were regarded as putatively informative if the markers could be assigned to a physical or a genetic position (Additional file 3: Figure S3).
Results and discussion
All 12,908 15K iSelect marker sequences  were used for BLASTn search (Additional file 3: Figure S3). For 12,445 of these, BLASTn-hits were obtained. A set of 2,813 markers revealed hits to only one chromosome. Eighty-nine markers gave hits to no specific chromosome at all. Out of the remaining 9,632 markers with hits to more than one chromosome, 4,898 markers could be linked to the corresponding genetic chromosomes of the durum consensus map . Finally, sets of 2,495 and 4,898 markers were screened to exclude those with multiple hits on the target chromosomes in a range of more than 1,000,000 bp, leading to a final set of 7,245 markers (Additional file 3: Figure S3, Additional file 4: Table S1 and Additional file 5: Table S2). Two-thousand-one-hundred-twenty of the 7,245 physically anchored iSelect markers were monomorphic or harboured failed genotypic data and were therefore regarded as non-informative. Of the remaining 5,125 markers, 2,971 were polymorphic between G18-16 and LDN. Anchoring the QTL-flanking markers that were previously used to transfer the QTL-regions into the NILs  to the physical map of Zavitan  revealed that they span over the centromeric region on both chromosomes (Additional file 6: Table S3). In line with these results, about two-third of the wild emmer chromosome 2B and 7A were transferred from F7RIL55 and F7RIL12 (Fig. 4) into NIL-U-2B-1 and NIL-B-7A-2, respectively (Fig. 5). Wild emmer fragments were also detected on chromosome 2A, 3A and 5A in NIL-U-2B-1 and 4B in NIL-B-7A-2 (Fig. 5). Although a remarkable reduction in the genetic background of the wild emmer parent in both NILs was observed at BC3F5 after the MAS procedure (Additional file 1: Figure S1, Fig. 5), these regions might cause linkage drag. To pinpoint regions that are potentially involved in linkage drag, primer pairs of well-described domestication genes  that are located on the NIL-U-2B-1 chromosome arms with introgressions from the wild emmer parent (Fig. 5) were aligned to the reference genome of wild emmer. It turned out that NIL-U-2B-1 contains wild emmer alleles of the domestication genes PpdA1 (TRIDC2AG009040)  and PPdB1 (TRIDC2BG010800)  on chromosomes 2AS and 2BS, respectively. In addition to this, a heterozygous version of the so called Lux (TRIDC3AG074490) gene, which was described to additively interact with PPdB1 , was identified at the very distal end of chromosome 3AL. Notably, NIL-U-2B-1 was also shown to suffer from delayed heading time . Presence/absence analysis of the respective regions in future segmental RILs might allow the identification of the exact causative loci. The allelic origin of other popular domestication genes, such as VrnA1  or Q on chromosome 5A  or the popular dwarfing gene RhtB1 , were also tested. However, none of these was transferred from G18-16 into the NILs.
This paper describes the development of software, that easily allows the display and comparison of genotypic data on genetic and physical maps by including a maximum amount of marker information, which is usually removed at the very early stages of the mapping procedure.
We presented a case study elucidating the use of GTM by using genotypic information of NILs carrying introgressions of wild emmer wheat, aiming to improve drought resistance in wheat cultivars [28, 29]. Both NILs were shown to be more tolerant to drought than their recurrent parents but showed phenotypes that may stem from a linkage drag, e.g., an increased plant height or an increased number of days to heading (DPH) [28, 29]. GTM enabled the visualization of the genotypic and physical results of a MAS procedure, which started from wild emmer wheat into tetra- or hexaploid wheat cultivars (Additional file 1: Figure S1). It depicts the special relevance and usability of physical and genetic genotyping for pre-breeding purposes to handle linkage drag and heterozygosity transferred from crop wild relatives.
Marker information illustrated by GTM is the product of the displayed genotypes, the chromosomes and the number of markers analysed. In the current version of the program, a maximum of ten genotypes can be displayed on the screen in the chromosome mode. If more than one genotype is analysed, the allelic discrimination analysis needs to be repeated for each individual. GTM uses multithreading to handle these situations efficiently.
To doublecheck GTM computational performance, artificial datasets with 3083, 6166, 9249 and 12,332 markers were created and used to explore GTM performance.
On a computer with an Intel five processor of the 8th generation, GTM can analyse a single genotype with 3,083 polymorphic markers on 14 chromosomes in less than 752 ms (ms) (Additional file 7: Figure S4). Ten individuals with the same amount of markers can be analysed in 1,429 ms, which is approximately twice as much of the time that is needed for a single individual (Additional file 7: Figure S4). The computation time increases linearly with the amount of analysed individuals. In addition to this, the number of markers has a significant impact on GTM computation time per individual and might hamper its performance if more than a single individual is analysed in parallel (Additional file 7: Figure S4). The tests were run under the standard conditions of the 1.8.0_231 JRE with a maximum heap space of 268,435,456 Byte.
We conclude that GTM can easily analyse up to ten genotypes on datasets with up to 6,000–7,000 markers in less than six seconds (Additional file 7: Figure S4). Analysing more markers with GTM is possible but requires more computation time. In such a case, we recommend analysing fewer genotypes and/or chromosomes in parallel.
Currently, genotypic data of 500 individuals with up to 50,000 markers can be uploaded to the software.
GTM’s allelic discrimination analysis works accurate, but sometimes different allele information might appear at the same genetic position, e.g., if consensus marker positions are used. In those cases, adequate and reliable predictions about the correct corresponding allele type are theoretically not possible/difficult and require the user’s expertise (Additional file 2: Figure S2).
We illustrate how such an analysis might be performed in detail (Additional file 3: Figure S3) with two and four different parents (Figs. 4 and 5). One additional example that illustrates how genotyping analysis might be performed if genotype information of only a single parent is available (Fig. 1 and Table 2) is given in the attachment for a NIL that carries an introgression harbouring the leaf rust resistance gene Lr1  (Additional file 8: Figure S5). Markers in this example were anchored to genetic positions via the wheat consensus map . Positions of the markers in the Chinese Spring genome are publicly available [44, 45].
In the case of this study, only 7,245, i.e., merely 60% of the 12,908 iSelect markers, could be assigned to reliable physical positions. This is very likely due to the following reasons: (i) part of the iSelect markers are designed based on D genome sequences; (ii) markers might have been assigned to different chromosomes in the Zavitan genome than in the durum consensus map; (iii) iSelect markers deliver bad BLAST-hits that likely fall below filtering criteria, and (iv) markers of gene families residing on several chromosomes and/or on homologous chromosomes on the AA and BB genomes cannot be assigned to a clear physical position.
Please note that for many BeadArrays or Gene chips, physical positions on the reference genomes are provided on the Illumina website (www.illumina.com) (e.g., maize) or Affymetrix website (www.Affymterix.com) (e.g., strawberry, apple, wheat, cotton, soybean, maize). In addition to this, physical positions of diverse marker types can also be obtained from databases, such as GrainGenes . In such cases, the identification of physical positions of the markers is not required. Anchoring of markers to physical position in the genome can also easily be accomplished with a simple BLASTN- search. In some cases, this approach might be sufficient .
GTM also offers the possibility to analyse single genetic or physical maps. In addition to this, two genetic maps can be compared. Respective test datasets are available on www.genotypemapper.org. Please note that GTM can also be used for other organisms than plants. As an example, we genotyped two F2 pig individuals that were derived from a cross between a “Göttingen Minipig” boar and “Yorkshire” gilt . The results were added to the supplement of this article (Additional file 9: Figure S6). We conclude that the software can be used for human and animal genetics as well.
Availability and requirements
Project name: GenoTypeMapper (GTM).
Project home: www.genotypemapper.org
Operating system(s): Platform independent
Programming language: Java
Other requirements: Java Runtime Environment 1.8.0_231 (or higher)
License: End use license agreements (EULA)
Any restrictions to use by non-academics: no
Availability of data and materials
The program, user manual and example data sets are freely available on https://www.genotypemapper.org.
Basic locale alignment search tool
Days planting to heading
Langdon (name of the drought susceptible durum cultivar)
Logarithm of odds
Marker assisted selection
Near isogenic line
Jaccoud D, Peng K, Feinstein D, Kilian A. Diversity arrays: a solid state technology for sequence information independent genotyping. Nucleic Acids Res. 2001;29:E25. https://doi.org/10.1093/nar/29.4.e25.
Elshire RJ, Glaubitz JC, Sun Q, Poland JA, Kawamoto K, Buckler ES, Mitchell SE. A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS ONE. 2011;6:e19379. https://doi.org/10.1371/journal.pone.0019379.
Shen R, Fan J-B, Campbell D, Chang W, Chen J, Doucet D, et al. High-throughput SNP genotyping on universal bead arrays. Mutat Res. 2005;573:70–82. https://doi.org/10.1016/j.mrfmmm.2004.07.022.
Steemers FJ, Gunderson KL. Whole genome genotyping technologies on the BeadArray platform. Biotechnol J. 2007;2:41–9. https://doi.org/10.1002/biot.200600213.
Matsuzaki H, Dong S, Loi H, Di X, Liu G, Hubbell E, et al. Genotyping over 100,000 SNPs on a pair of oligonucleotide arrays. Nat Methods. 2004;1:109–11. https://doi.org/10.1038/nmeth718.
Hoffmann TJ, Kvale MN, Hesselson SE, Zhan Y, Aquino C, Cao Y, et al. Next generation genome-wide association tool: design and coverage of a high-throughput European-optimized SNP array. Genomics. 2011;98:79–89. https://doi.org/10.1016/j.ygeno.2011.04.005.
Rasheed A, Hao Y, Xia X, Khan A, Xu Y, Varshney RK, He Z. Crop breeding chips and genotyping platforms: progress, challenges, and perspectives. Mol Plant. 2017;10:1047–64. https://doi.org/10.1016/j.molp.2017.06.008.
Morgan TH. Random segregation versus coupling in Mendelian inheritance. Science. 1911;34:384. https://doi.org/10.1126/science.34.873.384.
Sturtevant AH. The linear arrangement of six sex-linked factors in Drosophila, as shown by their mode of association. J Exp Zool. 1913;14:43–59. https://doi.org/10.1002/jez.1400140104.
Smeds L, Mugal CF, Qvarnström A, Ellegren H. High-resolution mapping of crossover and non-crossover recombination events by whole-genome re-sequencing of an avian pedigree. PLoS Genet. 2016;12:e1006044. https://doi.org/10.1371/journal.pgen.1006044.
Shen C, Li X, Zhang R, Lin Z. Genome-wide recombination rate variation in a recombination map of cotton. PLoS ONE. 2017;12:e0188682. https://doi.org/10.1371/journal.pone.0188682.
Chen M, Presting G, Barbazuk WB, Goicoechea JL, Blackmon B, Fang G, et al. An integrated physical and genetic map of the rice genome. Plant Cell. 2002;14:537–45. https://doi.org/10.1105/tpc.010485.
Schnable PS, Ware D, Fulton RS, Stein JC, Wei F, Pasternak S, et al. The B73 maize genome: complexity, diversity, and dynamics. Science. 2009;326:1112–5. https://doi.org/10.1126/science.1178534.
Kumar J, Gupta DS, Gupta S, Dubey S, Gupta P, Kumar S. Quantitative trait loci from identification to exploitation for crop improvement. Plant Cell Rep. 2017;36:1187–213. https://doi.org/10.1007/s00299-017-2127-y.
Ling H-Q, Ma B, Shi X, Liu H, Dong L, Sun H, et al. Genome sequence of the progenitor of wheat A subgenome Triticum urartu. Nature. 2018;557:424–8. https://doi.org/10.1038/s41586-018-0108-0.
Avni R, Nave M, Barad O, Baruch K, Twardziok SO, Gundlach H, et al. Wild emmer genome architecture and diversity elucidate wheat evolution and domestication. Science. 2017;357:93–7. https://doi.org/10.1126/science.aan0032.
Maccaferri M, Harris NS, Twardziok SO, Pasam RK, Gundlach H, Spannagl M, et al. Durum wheat genome highlights past domestication signatures and future improvement targets. Nat Genet. 2019;51:885–95. https://doi.org/10.1038/s41588-019-0381-3.
Luo M-C, Gu YQ, Puiu D, Wang H, Twardziok SO, Deal KR, et al. Genome sequence of the progenitor of the wheat D genome Aegilops tauschii. Nature. 2017;551:498–502. https://doi.org/10.1038/nature24486.
Mascher M, Gundlach H, Himmelbach A, Beier S, Twardziok SO, Wicker T, et al. A chromosome conformation capture ordered sequence of the barley genome. Nature. 2017;544:427–33. https://doi.org/10.1038/nature22043.
Appels R, Eversole K, Feuillet C, Keller B, Rogers J, Stein N, et al. Shifting the limits in wheat research and breeding using a fully annotated reference genome. Science. 2018. https://doi.org/10.1126/science.aar7191.
Thind AK, Wicker T, Müller T, Ackermann PM, Steuernagel B, Wulff BBH, et al. Chromosome-scale comparative sequence analysis unravels molecular mechanisms of genome dynamics between two wheat cultivars. Genome Biol. 2018;19:104. https://doi.org/10.1186/s13059-018-1477-2.
van Berloo R. GGT 2.0: versatile software for visualization and analysis of genetic data. J Hered. 2008;99:232–6. https://doi.org/10.1093/jhered/esm109.
Milne I, Shaw P, Stephen G, Bayer M, Cardle L, Thomas WTB, et al. Flapjack–graphical genotype visualization. Bioinformatics. 2010;26:3133–4. https://doi.org/10.1093/bioinformatics/btq580.
Thorvaldsdóttir H, Robinson JT, Mesirov JP. Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinformatics. 2013;14:178–92. https://doi.org/10.1093/bib/bbs017.
Bottoms CA, Flint-Garcia S, McMullen MD. IView: introgression library visualization and query tool. BMC Bioinformatics. 2010;11(Suppl 6):S28. https://doi.org/10.1186/1471-2105-11-S6-S28.
Aflitos SA, Sanchez-Perez G, de Ridder D, Fransz P, Schranz ME, de Jong H, Peters SA. Introgression browser: high-throughput whole-genome SNP visualization. Plant J. 2015;82:174–82. https://doi.org/10.1111/tpj.12800.
Rezvoy C, Charif D, Guéguen L, Marais GAB. MareyMap: an R-based tool with graphical interface for estimating recombination rates. Bioinformatics. 2007;23:2188–9. https://doi.org/10.1093/bioinformatics/btm315.
Merchuk-Ovnat L, Barak V, Fahima T, Ordon F, Lidzbarsky GA, Krugman T, Saranga Y. Ancestral QTL alleles from wild emmer wheat improve drought resistance and productivity in modern wheat cultivars. Front Plant Sci. 2016;7:452. https://doi.org/10.3389/fpls.2016.00452.
Merchuk-Ovnat L, Fahima T, Ephrath JE, Krugman T, Saranga Y. Ancestral QTL alleles from wild emmer wheat enhance root development under drought in modern wheat. Front Plant Sci. 2017;8:703. https://doi.org/10.3389/fpls.2017.00703.
Peleg Z, Saranga Y, Suprunova T, Ronin Y, Röder MS, Kilian A, et al. High-density genetic map of durum wheat x wild emmer wheat based on SSR and DArT markers. Theor Appl Genet. 2008;117:103–15. https://doi.org/10.1007/s00122-008-0756-9.
Peleg Z, Fahima T, Krugman T, Abbo S, Yakir D, Korol AB, Saranga Y. Genomic dissection of drought resistance in durum wheat x wild emmer wheat recombinant inbreed line population. Plant Cell Environ. 2009;32:758–79. https://doi.org/10.1111/j.1365-3040.2009.01956.x.
Fatiukha A, Deblieck M, Klymiuk V, Merchuk-Ovnat L, Peleg Z, Ordon F, et al. Genomic architecture of phenotypic plasticity of complex traits in tetraploid wheat in response to water stress. bioRxiv. 2019. https://doi.org/10.1101/565820.
Soleimani B, Lehnert H, Keilwagen J, Plieske J, Ordon F, Naseri Rad S, et al. Comparison between core set selection methods using different Illumina marker platforms: a case study of assessment of diversity in wheat. Front Plant Sci. 2020. https://doi.org/10.3389/fpls.2020.01040.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10. https://doi.org/10.1016/S0022-2836(05)80360-2.
Maccaferri M, Cane' MA, Sanguineti MC, Salvi S, Colalongo MC, Massi A, et al. A consensus framework map of durum wheat (Triticum durum Desf.) suitable for linkage disequilibrium analysis and genome-wide association mapping. BMC Genom. 2014;15:873. https://doi.org/10.1186/1471-2164-15-873.
Gill BS, Appels R, Botha-Oberholster A-M, Buell CR, Bennetzen JL, Chalhoub B, et al. A workshop report on wheat genome sequencing: International Genome Research on Wheat Consortium. Genetics. 2004;168:1087–96. https://doi.org/10.1534/genetics.104.034769.
Muterko A, Kalendar R, Cockram J, Balashova I. Discovery, evaluation and distribution of haplotypes and new alleles of the Photoperiod-A1 gene in wheat. Plant Mol Biol. 2015;88:149–64. https://doi.org/10.1007/s11103-015-0313-2.
Díaz A, Zikhali M, Turner AS, Isaac P, Laurie DA. Copy number variation affecting the Photoperiod-B1 and Vernalization-A1 genes is associated with altered flowering time in wheat (Triticum aestivum). PLoS ONE. 2012;7:e33234. https://doi.org/10.1371/journal.pone.0033234.
Zhou W, Wu S, Ding M, Li J, Shi Z, Wei W, et al. Mapping of Ppd-B1, a major candidate gene for late heading on wild emmer chromosome Arm 2BS and assessment of its interactions with early heading QTLs on 3AL. PLoS ONE. 2016;11:e0147377. https://doi.org/10.1371/journal.pone.0147377.
Simons KJ, Fellers JP, Trick HN, Zhang Z, Tai Y-S, Gill BS, Faris JD. Molecular characterization of the major wheat domestication gene Q. Genetics. 2006;172:547–55. https://doi.org/10.1534/genetics.105.044727.
Pearce S, Saville R, Vaughan SP, Chandler PM, Wilhelm EP, Sparks CA, et al. Molecular characterization of Rht-1 dwarfing genes in hexaploid wheat. Plant Physiol. 2011;157:1820–31. https://doi.org/10.1104/pp.111.183657.
Feuillet C, Messmer M, Schachermayr G, Keller B. Genetic and physical characterization of the LR1 leaf rust resistance locus in wheat (Triticum aestivum L.). Mol Gen Genet. 1995;248:553–62. https://doi.org/10.1007/BF02423451.
Wang S, Wong D, Forrest K, Allen A, Chao S, Huang BE, et al. Characterization of polyploid wheat genomic diversity using a high-density 90,000 single nucleotide polymorphism array. Plant Biotechnol J. 2014;12:787–96. https://doi.org/10.1111/pbi.12183.
Alaux M, Rogers J, Letellier T, Flores R, Alfama F, Pommier C, et al. Linking the International Wheat Genome Sequencing Consortium bread wheat reference genome sequence to wheat genetic and phenomic data. Genome Biol. 2018;19:111. https://doi.org/10.1186/s13059-018-1491-4.
Matthews DE, Carollo VL, Lazo GR, Anderson OD. GrainGenes, the genome database for small-grain crops. Nucleic Acids Res. 2003;31:183–6. https://doi.org/10.1093/nar/gkg058.
Fazlikhani L, Keilwagen J, Kopahnke D, Deising H, Ordon F, Perovic D. High Resolution Mapping of RphMBR1012 Conferring Resistance to Puccinia hordei in Barley (Hordeum vulgare L.). Front Plant Sci. 2019;10:640. https://doi.org/10.3389/fpls.2019.00640.
Pant SD, Karlskov-Mortensen P, Jacobsen MJ, Cirera S, Kogelman LJA, Bruun CS, et al. Comparative analyses of QTLs influencing obesity and metabolic phenotypes in pigs and humans. PLoS ONE. 2015;10:e0137356. https://doi.org/10.1371/journal.pone.0137356.
This work was financially supported by the German Federal Ministry of Food and Agriculture (BMEL) based on the decision of the Parliament of the Federal Republic of Germany (FKZ: 2813IL03) and the Israel Ministry of Agriculture and Rural Development, Chief Scientist Foundation (Grants 837-0079-10, 837-0162-14). YS is the incumbent of the Haim Gvati Chair in Agriculture. Open access funding provided by Projekt DEAL.
Ethics approval and consent to participate
Consent for publication
All authors agreed to publish this manuscript.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Establishment of near isogenic lines carrying wild emmer QTLs for drought stress tolerance. (1) The wild emmer acc. G18-16 was crossed with the cv. Langdon (LDN) to obtain a segregating F6 mapping population. A map with 197 microsatellites and 493 DArT markers on 690 loci was constructed to identify QTLs for drought stress tolerance associated traits [30, 31]. (2/3) Markers that flank the QTL-regions of interest were then used to introgress the QTLs into the background of elite wheat cultivars. The plants were subsequently backcrossed three times to diminish the genomic fraction of the donor parent. Finally, BC3F3 and BC3F4 progenies were phenotyped under water-limited and well-watered conditions together with their recurrent parents . NIL-U-2B-1 and NIL-B-7A-2 not only showed better tolerance than their recurrent elite parents but also exhibited negative traits . (4) 15k iSelect genotypic data of F7 descendants of the original mapping population are used to reconstruct a new high-resolution genetic map . Green and purple boxes represent male and female plants, respectively.
Interval- and marker-based colouring in GTM. (A) In a classical genetic mapping procedure, heterozygous & monomorphic markers of parents are discarded. The remaining molecular markers form linkage groups of the same allele types so that GTM can highlight intervals that are flanked by markers with the same allele information. Intervals that are flanked by markers with different allele types are not coloured because the exact position of recombination between both markers is unknown. (B) If genetic consensus marker positions are used, markers with different allele information can be assigned to the same genetic position. In those cases, colouring of intervals is not possible, and marker-based colouring might lead to inaccurate results (see Marker M1, M2, M3). Filtering of the allele types of interest might be the method of choice to analyse that kind of genotype, as illustrated in (C) and (D).
BLAST-scheme that illustrates the assignment of physical loci of 15k-iSelect markers in the Triticum turgidum ssp. dicoccoides genome. Numbers obtained in the analysis are coloured in orange. For more details, see text.
15k iSelect -markers that were assigned to physical or genetic loci.
In total, 7,245 iSelect markers with physical and genetic positions.
Genetic and physical positions of flanking microsatellite and DArT markers of QTL regions 2BS and 7AS.
GTM’s performance relative to the number of markers and genotypes that are analysed.
Allelic discrimination analysis with one parent in GenoTypeMapper. Wheat near isogenic line Tc-Lr1 with an introgression of the leaf rust resistance gene Lr1 in the background of the parent “Thatcher” on chromosome 5DL. Only genotype data of the Thatcher line was available, so allele analysis with one parent was performed (Figure 2, Table 2). Please notice that monomorphic markers between Thatcher and the second parent were not detected due to the missing genotype information of one parent. Therefore, filtering and colouring of the introgressed markers was the method of choice to illustrate the introgressed regions. This analysis illustrates that introgressions of the unknown parent also occurred to other chromosomes. The physical genome shows that centromeric regions were introgressed from the resistant parent into the NIL. Please note: This test-dataset is quite large and comprises more than 30540 markers.
Physical genotyping of two F2 pigs: The F2 boars resulted from a cross of a Minipig (MP) (boar) and a Yorkshire (YS) (gilt). In the parental generation seven purebred Yorkshire (YY) sows were mated to 14 Goettingen Minipig (MM) boars. Among the YM F1 animals, 26 gilts and 13 boars were mated to produce 279 F2 animals. Genotypic data (60k iSelect data) for the animals were previously published and reused in this analysis . The marker data for this illustration were obtained in two steps: First, consensus marker information from the parental YS boars and MP gilts was extracted. Then, markers that failed in YS or MP were removed. Information about the colouring of respective alleles is illustrated in the figure. It should be pointed out, that the Y chromosome is not shown due to the fact that it has to be derived from the MP.
About this article
Cite this article
Deblieck, M., Fatiukha, A., Grundman, N. et al. GenoTypeMapper: graphical genotyping on genetic and sequence-based maps. Plant Methods 16, 123 (2020). https://doi.org/10.1186/s13007-020-00665-7