Comparative Genetic Mapping Research Articles

Postzygotic barriers persist despite ongoing introgression in hybridizing Mimulus species.

The evolution of postzygotic isolation is thought to be a key step in maintaining species boundaries upon secondary contact, yet the dynamics and persistence of hybrid incompatibilities in naturally hybridizing species are not well understood. Here, we explore these issues using genetic mapping in three independent populations of recombinant inbred lines between naturally hybridizing monkeyflowers, Mimulus guttatus and Mimulus nasutus, from the sympatric Catherine Creek population. We discover that the three M. guttatus founders differ dramatically in admixture history, with nearly a quarter of one founder's genome introgressed from M. nasutus. Comparative genetic mapping in the three RIL populations reveals three new putative inversions, each one segregating among the M. guttatus founders, two due to admixture. We find strong, genome-wide transmission ratio distortion in all RILs, but patterns are highly variable among the three populations. At least some of this distortion appears to be explained by epistatic selection favouring parental genotypes, but tests of inter-chromosomal linkage disequilibrium also reveal multiple candidate Dobzhansky-Muller incompatibilities. We also map several genetic loci for hybrid pollen viability, including two interacting pairs that coincide with peaks of distortion. Remarkably, even with this limited sample of three M. guttatus lines, we discover abundant segregating variation for hybrid incompatibilities with M. nasutus, suggesting this population harbours diverse contributors to postzygotic isolation. Moreover, even with substantial admixture, hybrid incompatibilities between Mimulus species persist, suggesting postzygotic isolation might be a potent force in maintaining species barriers in this system.

One hundred years of comparative genetic and physical mapping in cultivated oat (Avena sativa)

Context Researchers have been accumulating information concerning the locations of genes and quantitative trait loci (QTLs) in cultivated oat (Avena sativa L.) for more than 100 years. Aims The aim of this work was to create an inventory of genes and QTLs found in cultivated hexaploid oat and produce tools to make this resource more useful. Methods By using the positions of perfectly matched, single nucleotide polymorphism markers, each centimorgan (cM) location along the consensus map was assigned to a location on the OT3098 v2 physical map found on the GrainGenes database website (https://wheat.pw.usda.gov/jb/?data=/ggds/oat-ot3098v2-pepsico). This information was then used to assign physical locations to the genes and QTLs in the inventory, where possible. Key results A table comparing the major genetic maps of hexaploid oats to each other, to the 2018 oat consensus map, and to physical chromosomes was produced. Genome browser tracks aligning the consensus map regions and the locations of the genes and QTLs to OT3098 v2 were added to GrainGenes. Conclusions Many oat genes and QTLs identified using genetic mapping could be assigned positions on physical oat chromosomes. However, many of these assigned regions are quite long, owing to the presence of large areas of reduced recombination. Specific examples of identified patterns of recombination between the genetic and physical maps and validated gene and QTL locations are discussed. Implications These resources will assist researchers performing comparative genetic and physical mapping in oat.

Heterochiasmy and Sex Chromosome Evolution in Silene

The evolution of a non-recombining sex-specific region is a key step in sex chromosome evolution. Suppression of recombination between the (proto-) X- and Y-chromosomes in male meiosis creates a non-recombining Y-linked region (NRY), while the X-chromosome continues to recombine in females. Lack of recombination in the NRY defines its main properties-genetic degeneration and accumulation of repetitive DNA, making X and Y chromosomes very different from each other. How and why recombination suppression on sex chromosomes evolves remains controversial. A strong difference in recombination rates between the sexes (heterochiasmy) can facilitate or even cause recombination suppression. In the extreme case-complete lack of recombination in the heterogametic sex (achiasmy)-the entire sex-specific chromosome is automatically non-recombining. In this study, I analyse sex-specific recombination rates in a dioecious plant Silene latifolia (Caryophyllaceae), which evolved separate sexes and sex chromosomes ~11 million years ago. I reconstruct high-density RNAseq-based genetic maps including over five thousand genic markers for the two sexes separately. The comparison of the male and female maps reveals only modest heterochiasmy across the genome, with the exception of the sex chromosomes, where recombination is suppressed in males. This indicates that heterochiasmy likely played only a minor, if any, role in NRY evolution in S. latifolia, as recombination suppression is specific to NRY rather than to the entire genome in males. Other mechanisms such as structural rearrangements and/or epigenetic modifications were likely involved, and comparative genome analysis and genetic mapping in multiple Silene species will help to shed light on the mechanism(s) of recombination suppression that led to the evolution of sex chromosomes.

Two Genomic Loci Control Three Eye Colors in the Domestic Pigeon (Columba livia).

The iris of the eye shows striking color variation across vertebrate species, and may play important roles in crypsis and communication. The domestic pigeon (Columba livia) has three common iris colors, orange, pearl (white), and bull (dark brown), segregating in a single species, thereby providing a unique opportunity to identify the genetic basis of iris coloration. We used comparative genomics and genetic mapping in laboratory crosses to identify two candidate genes that control variation in iris color in domestic pigeons. We identified a nonsense mutation in the solute carrier SLC2A11B that is shared among all pigeons with pearl eye color, and a locus associated with bull eye color that includes EDNRB2, a gene involved in neural crest migration and pigment development. However, bull eye is likely controlled by a heterogeneous collection of alleles across pigeon breeds. We also found that the EDNRB2 region is associated with regionalized plumage depigmentation (piebalding). Our study identifies two candidate genes for eye colors variation, and establishes a genetic link between iris and plumage color, two traits that vary widely in the evolution of birds and other vertebrates.

High macro-collinearity between Crassostrea angulata and C. gigas genomes was revealed by comparative genetic mapping with transferable EST-SNP markers

Comparative genomics has become an important strategy for transferring genetic information from a well-characterized model species to one that is less described and can facilitate in functional gene mining and genomic evolution. Crassostrea angulata and Crassostrea gigas are closed related species that they can hybridize to produce fertile offspring, and C. gigas has accumulated a considerable number of genetic and genomic resources while C. angulata lags far behind. In the current study, 684 single nucleotide polymorphism (SNP) markers developed from C. gigas expressed sequence tags (ESTs) were used for transferability test in C. angulata, and 635 loci were successfully amplified. Then a genetic linkage map of C. angulata was for the first time constructed using these transferable SNP markers. A total of 988 transferable SNP markers were used for polymorphism screening in a mapping family, and 273 loci were distributed on the genetic map which was composed of 10 linkage groups and estimated 1040.29 cM in total length. The average spacing between markers was 3.96 cM. By comparing the shared EST-SNP marker locations on homologous linkage groups, comparative mapping analysis between the C. angulata genetic map and an integrated map of C. gigas revealed putative macro-collinearity in the two oyster genomes with little difference in markers order. These transferable SNP markers could be applied in genetic diversity and population structure studies, and the genetic map would be useful for quantitative trait loci mapping and genetic improvement through marker-assisted selection for C. angulata.

Widespread Recombination Suppression Facilitates Plant Sex Chromosome Evolution.

Classical models suggest that recombination rates on sex chromosomes evolve in a stepwise manner to localize sexually antagonistic variants in the sex in which they are beneficial, thereby lowering rates of recombination between X and Y chromosomes. However, it is also possible that sex chromosome formation occurs in regions with preexisting recombination suppression. To evaluate these possibilities, we constructed linkage maps and a chromosome-scale genome assembly for the dioecious plant Rumex hastatulus. This species has a polymorphic karyotype with a young neo-sex chromosome, resulting from a Robertsonian fusion between the X chromosome and an autosome, in part of its geographic range. We identified the shared and neo-sex chromosomes using comparative genetic maps of the two cytotypes. We found that sex-linked regions of both the ancestral and the neo-sex chromosomes are embedded in large regions of low recombination. Furthermore, our comparison of the recombination landscape of the neo-sex chromosome to its autosomal homolog indicates that low recombination rates mainly preceded sex linkage. These patterns are not unique to the sex chromosomes; all chromosomes were characterized by massive regions of suppressed recombination spanning most of each chromosome. This represents an extreme case of the periphery-biased recombination seen in other systems with large chromosomes. Across all chromosomes, gene and repetitive sequence density correlated with recombination rate, with patterns of variation differing by repetitive element type. Our findings suggest that ancestrally low rates of recombination may facilitate the formation and subsequent evolution of heteromorphic sex chromosomes.

Multiple chromosomal inversions contribute to adaptive divergence of a dune sunflower ecotype.

Both models and case studies suggest that chromosomal inversions can facilitate adaptation and speciation in the presence of gene flow by suppressing recombination between locally adapted alleles. Until recently, however, it has been laborious and time-consuming to identify and genotype inversions in natural populations. Here we apply RAD sequencing data and newly developed population genomic approaches to identify putative inversions that differentiate a sand dune ecotype of the prairie sunflower (Helianthus petiolaris) from populations found on the adjacent sand sheet. We detected seven large genomic regions that exhibit a different population structure than the rest of the genome and that vary in frequency between dune and nondune populations. These regions also show high linkage disequilibrium and high heterozygosity between, but not within, arrangements, consistent with the behaviour of large inversions, an inference subsequently validated in part by comparative genetic mapping. Genome-environment association analyses show that key environmental variables, including vegetation cover and soil nitrogen, are significantly associated with inversions. The inversions colocate with previously described "islands of differentiation," and appear to play an important role in adaptive divergence and incipient speciation within H. petiolaris.

A Similar Genetic Architecture Underlies the Convergent Evolution of the Selfing Syndrome in Capsella.

Whether, and to what extent, phenotypic evolution follows predictable genetic paths remains an important question in evolutionary biology. Convergent evolution of similar characters provides a unique opportunity to address this question. The transition to selfing and the associated changes in flower morphology are among the most prominent examples of repeated evolution in plants. In this study, we take advantage of the independent transitions to self-fertilization in the genus Capsella to compare the similarities between parallel modifications of floral traits and test for genetic and developmental constraints imposed on flower evolution in the context of the selfing syndrome. Capsella rubella and Capsella orientalis emerged independently but evolved almost identical flower characters. Not only is the evolutionary outcome identical but the same developmental strategies underlie the convergent reduction of flower size. This has been associated with convergent evolution of gene expression changes. The transcriptomic changes common to both selfing lineages are enriched in genes with low network connectivity and with organ-specific expression patterns. Comparative genetic mapping also suggests that, at least in the case of petal size evolution, these similarities have a similar genetic basis. Based on these results, we hypothesize that the limited availability of low-pleiotropy paths predetermines closely related species to similar evolutionary outcomes.

Genome Improvement and Genetic Map Construction for Aethionema arabicum, the First Divergent Branch in the Brassicaceae Family.

The genus Aethionema is a sister-group to the core-group of the Brassicaceae family that includes Arabidopsis thaliana and the Brassica crops. Thus, Aethionema is phylogenetically well-placed for the investigation and understanding of genome and trait evolution across the family. We aimed to improve the quality of the reference genome draft version of the annual species Aethionema arabicum. Second, we constructed the first Ae. arabicum genetic map. The improved reference genome and genetic map enabled the development of each other. We started with the initially published genome (version 2.5). PacBio and MinION sequencing together with genetic map v2.5 were incorporated to produce the new reference genome v3.0. The improved genome contains 203 MB of sequence, with approximately 94% of the assembly made up of called (non-gap) bases, assembled into 2,883 scaffolds (with only 6% of the genome made up of non-called bases (Ns)). The N50 (10.3 MB) represents an 80-fold increase over the initial genome release. We generated a Recombinant Inbred Line (RIL) population that was derived from two ecotypes: Cyprus and Turkey (the reference genotype. Using a Genotyping by Sequencing (GBS) approach, we generated a high-density genetic map with 749 (v2.5) and then 632 SNPs (v3.0) was generated. The genetic map and reference genome were integrated, thus greatly improving the scaffolding of the reference genome into 11 linkage groups. We show that long-read sequencing data and genetics are complementary, resulting in an improved genome assembly in Ae. arabicum. They will facilitate comparative genetic mapping work for the Brassicaceae family and are also valuable resources to investigate wide range of life history traits in Aethionema.

GOOGA: A platform to synthesize mapping experiments and identify genomic structural diversity.

Understanding genomic structural variation such as inversions and translocations is a key challenge in evolutionary genetics. We develop a novel statistical approach to comparative genetic mapping to detect large-scale structural mutations from low-level sequencing data. The procedure, called Genome Order Optimization by Genetic Algorithm (GOOGA), couples a Hidden Markov Model with a Genetic Algorithm to analyze data from genetic mapping populations. We demonstrate the method using both simulated data (calibrated from experiments on Drosophila melanogaster) and real data from five distinct crosses within the flowering plant genus Mimulus. Application of GOOGA to the Mimulus data corrects numerous errors (misplaced sequences) in the M. guttatus reference genome and confirms or detects eight large inversions polymorphic within the species complex. Finally, we show how this method can be applied in genomic scans to improve the accuracy and resolution of Quantitative Trait Locus (QTL) mapping.

Comparative genetic mapping reveals synteny and collinearity between the American cranberry and diploid blueberry genomes

Cranberry and blueberry are closely related and recently domesticated fruit crops in the genus Vaccinium. Both have a presumed American origin and likely evolved from a common ancestor; however, details of their adaptive radiation and the extent of their genomic divergence remains little understood. To better understand their evolutionary and genomic relationships, a set of 323 cross-transferable simple sequence repeat (SSR) markers were identified, added to existing marker datasets, and used to construct linkage maps for cranberry (582 SSRs) and an interspecific diploid blueberry population (V. darrowii x V. corymbosum) x V. corymbosum (409 markers, densest blueberry SSR map currently available). The maps allowed for the first comparative genetic mapping study in Vaccinium, and revealed a surprisingly high degree of macro-synteny and collinearity between the cranberry and blueberry genomes. Approximately 93% of the blueberry linkage map was collinear with cranberry, while the remaining 7% (66.3 cM) was spread across 15 non-collinear regions detected in eight of the 12 linkage groups. These observations suggest that large-scale genome differentiation between the cranberry and blueberry genomes has not occurred during their evolution, and that sequence information will be highly transferable between the species in future genetic research and breeding. Finally, the set of 323 cross-transferable SSRs and linkage maps they were used to construct can serve as a shared resource for the Vaccinium research community, enabling additional comparative mapping studies, the identification and transfer of quantitative trait loci and candidate genes between species, and future exploration of evolutionary relationships in Vaccinium.

High-Throughput RNA-Seq Data Analysis of the Single Nucleotide Polymorphisms (SNPs) and Zygomorphic Flower Development in Pea (Pisum sativum L.).

Pea (Pisum sativum L.) is a model plant that has been used in classical genetics and organ development studies. However, its large and complex genome has hindered research investigations in pea. Here, we generated transcriptomes from different tissues or organs of three pea accessions using next-generation sequencing to assess single nucleotide polymorphisms (SNPs), and further investigated petal differentially expressed genes to elucidate the mechanisms regulating floral zygomorphy. Eighteen samples were sequenced, which yielded a total of 617 million clean reads, and de novo assembly resulted in 87,137 unigenes. A total of 9044 high-quality SNPs were obtained among the three accessions, and a consensus map was constructed. We further discovered several dorsoventral asymmetrically expressed genes that were confirmed by qRT-PCR among different petals, including previously reported three CYC-like proliferating cell factor (TCP) genes. One MADS-box gene was highly expressed in dorsal petals, and several MYB factors were predominantly expressed among dorsal, lateral, and/or ventral petals, together with a ventrally expressed TCP gene. In sum, our comprehensive database complements the existing resources for comparative genetic mapping and facilitates future investigations in legume zygomorphic flower development.

Genetic Mapping and Phylogenetic Analysis Reveal Intraspecific Variation in Sex Chromosomes of the Virginian Strawberry.

With their extraordinary diversity in sexual systems, flowering plants offer unparalleled opportunities to understand sex determination and to reveal generalities in the evolution of sex chromosomes. Comparative genetic mapping of related taxa with good phylogenetic resolution can delineate the extent of sex chromosome diversity within plant groups, and lead the way to understanding the evolutionary drivers of such diversity. The North American octoploid wild strawberries provide such an opportunity. We performed linkage mapping using targeted sequence capture for the subdioecious western Fragaria virginiana ssp. platypetala and compared the location of its sex-determining region (SDR) to those of 2 other (sub)dioecious species, the eastern subspecies, F. virginiana ssp. virginiana (whose SDR is at 0-5.5 Mb on chromosome VI of the B2 subgenome), and the sister species F. chiloensis (whose SDR is at 37 Mb on chromosome VI of the Av subgenome). Male sterility was dominant in F. virginiana ssp. platypetala and mapped to a chromosome also in homeologous group VI. Likewise, one major quantitative trait locus (QTL) for female fertility overlapped the male sterility region. However, the SDR mapped to yet another subgenome (B1), and to a different location (13 Mb), but similar to the location inferred in one population of the naturally occurring hybrid between F. chiloensis and F. virginiana (F. ×ananassa ssp. cuneifolia). Phylogenetic analysis of chromosomes across the octoploid taxa showed consistent subgenomic composition reflecting shared evolutionary history but also reinforced within-species variation in the SDR-carrying chromosome, suggesting either repeated evolution, or recent turnovers in SDR.

Major Contribution of Flowering Time and Vegetative Growth to Plant Production in Common Bean As Deduced from a Comparative Genetic Mapping.

Determinacy growth habit and accelerated flowering traits were selected during or after domestication in common bean. Both processes affect several presumed adaptive traits such as the rate of plant production. There is a close association between flowering initiation and vegetative growth; however, interactions among these two crucial developmental processes and their genetic bases remain unexplored. In this study, with the aim to establish the genetic relationships between these complex processes, a multi-environment quantitative trait locus (QTL) mapping approach was performed in two recombinant inbred line populations derived from inter-gene pool crosses between determinate and indeterminate genotypes. Additive and epistatic QTLs were found to regulate flowering time, vegetative growth, and rate of plant production. Moreover, the pleiotropic patterns of the identified QTLs evidenced that regions controlling time to flowering traits, directly or indirectly, are also involved in the regulation of plant production traits. Further QTL analysis highlighted one QTL, on the lower arm of the linkage group Pv01, harboring the Phvul.001G189200 gene, homologous to the Arabidopsis thaliana TERMINAL FLOWER1 (TFL1) gene, which explained up to 32% of phenotypic variation for time to flowering, 66% for vegetative growth, and 19% for rate of plant production. This finding was consistent with previous results, which have also suggested Phvul.001G189200 (PvTFL1y) as a candidate gene for determinacy locus. The information here reported can also be applied in breeding programs seeking to optimize key agronomic traits, such as time to flowering, plant height and an improved reproductive biomass, pods, and seed size, as well as yield.

The locus for resistance to Asian soybean rust in PI 587855

AbstractAsian soybean rust (ASR) caused by Phakopsora pachyrhizi is one of the most serious soybean (Glycine max) diseases in tropical and subtropical areas. A soybean line, PI 587855, showed a resistance phenotype against ASR pathogens in Japan and South America at high frequency; however, little is known of the genetic control of this resistance and chromosomal location of the corresponding locus. Therefore, the aim of this study was to study the inheritance of PI 587855 resistance and map the corresponding locus with SSR markers aiming to use the linked markers in marker‐assisted selection. In the segregating population, resistance to ASR appeared to be controlled by a single dominant gene. The ASR resistance locus was mapped near to the chromosomal region where the resistant loci, Rpp1 and Rpp1‐b, were previously mapped. Comparative genetic mapping and disease reaction profiles of other seven lines carrying Rpp1 or Rpp1‐b to four Brazilian ASR isolates revealed that the resistance reaction exhibited by PI 587855 was similar to that of Rpp1‐b‐carrying varieties which have useful resistance to South American ASR strains.

FISH and AgNor mapping of the 45S and 5S rRNA genes in wild and cultivated species of Capsicum (Solananceae).

Chromosome number and position of rDNA were studied in 12 wild and cultivated species of the genus Capsicum with chromosome numbers x = 12 and x = 13 (22 samples). For the first time in these species, the 5S and 45S rRNA loci were localized and physically mapped using two-color fluorescence in situ hybridization and AgNOR banding. We focused on the comparison of the results obtained with both methods with the aim of accurately revealing the real functional rRNA genes. The analyzes were based on a previous work that reported that the 18S-5.8S-25S loci mostly coincide with GC-rich heterochromatic regions and likely have given rise to satellite DNAs, which are not active genes. These data show the variability of rDNA within karyotypes of the genus Capsicum, providing anchor points for (comparative) genetic maps. In addition, the obtained information might be useful for studies on evolution of repetitive DNA.

The Pace of Hybrid Incompatibility Evolution in House Mice.

Hybrids between species are often sterile or inviable. This form of reproductive isolation is thought to evolve via the accumulation of mutations that interact to reduce fitness when combined in hybrids. Mathematical formulations of this "Dobzhansky-Muller model" predict an accelerating buildup of hybrid incompatibilities with divergence time (the "snowball effect"). Although the Dobzhansky-Muller model is widely accepted, the snowball effect has only been tested in two species groups. We evaluated evidence for the snowball effect in the evolution of hybrid male sterility among subspecies of house mice, a recently diverged group that shows partial reproductive isolation. We compared the history of subspecies divergence with patterns of quantitative trait loci (QTL) detected in F2 intercrosses between two pairs of subspecies (Mus musculus domesticus with M. m. musculus and M. m. domesticus with M. m. castaneus). We used a recently developed phylogenetic comparative method to statistically measure the fit of these data to the snowball prediction. To apply this method, QTL were partitioned as either shared or unshared in the two crosses. A heuristic partitioning based on the overlap of QTL confidence intervals produced unambiguous support for the snowball effect. An alternative approach combining data among crosses favored the snowball effect for the autosomes, but a linear accumulation of incompatibilities for the X chromosome. Reasoning that the X chromosome analyses are complicated by low mapping resolution, we conclude that hybrid male sterility loci have snowballed in house mice. Our study illustrates the power of comparative genetic mapping for understanding mechanisms of speciation.

Plant contributions to our understanding of sex chromosome evolution.

A minority of angiosperms have male and female flowers separated in distinct individuals (dioecy), and most dioecious plants do not have cytologically different (heteromorphic) sex chromosomes. Plants nevertheless have several advantages for the study of sex chromosome evolution, as genetic sex determination has evolved repeatedly and is often absent in close relatives. I review sex-determining regions in non-model plant species, which may help us to understand when and how (and, potentially, test hypotheses about why) recombination suppression evolves within young sex chromosomes. I emphasize high-throughput sequencing approaches that are increasingly being applied to plants to test for non-recombining regions. These data are particularly illuminating when combined with sequence data that allow phylogenetic analyses, and estimates of when these regions evolved. Together with comparative genetic mapping, this has revealed that sex-determining loci and sex-linked regions evolved independently in many plant lineages, sometimes in closely related dioecious species, and often within the past few million years. In reviewing recent progress, I suggest areas for future work, such as the use of phylogenies to allow the informed choice of outgroup species suitable for inferring the directions of changes, including testing whether Y chromosome-like regions are undergoing genetic degeneration, a predicted consequence of losing recombination.

The functional basis of wing patterning in Heliconius butterflies: the molecules behind mimicry.

Wing-pattern mimicry in butterflies has provided an important example of adaptation since Charles Darwin and Alfred Russell Wallace proposed evolution by natural selection >150 years ago. The neotropical butterfly genus Heliconius played a central role in the development of mimicry theory and has since been studied extensively in the context of ecology and population biology, behavior, and mimicry genetics. Heliconius species are notable for their diverse color patterns, and previous crossing experiments revealed that much of this variation is controlled by a small number of large-effect, Mendelian switch loci. Recent comparative analyses have shown that the same switch loci control wing-pattern diversity throughout the genus, and a number of these have now been positionally cloned. Using a combination of comparative genetic mapping, association tests, and gene expression analyses, variation in red wing patterning throughout Heliconius has been traced back to the action of the transcription factor optix. Similarly, the signaling ligand WntA has been shown to control variation in melanin patterning across Heliconius and other butterflies. Our understanding of the molecular basis of Heliconius mimicry is now providing important insights into a variety of additional evolutionary phenomena, including the origin of supergenes, the interplay between constraint and evolvability, the genetic basis of convergence, the potential for introgression to facilitate adaptation, the mechanisms of hybrid speciation in animals, and the process of ecological speciation.

Comparative genetic mapping revealed powdery mildew resistance gene MlWE4 derived from wild emmer is located in same genomic region of Pm36 and Ml3D232 on chromosome 5BL

Powdery mildew, caused by Blumeria graminis f. sp. tritici, is one of the most devastating wheat diseases. Wild emmer wheat (Triticum turgidum ssp. dicoccoides) is a promising source of disease resistance for wheat. A powdery mildew resistance gene conferring resistance to B. graminis f. sp. tritici isolate E09, originating from wild emmer wheat, has been transferred into the hexaploid wheat line WE4 through crossing and backcrossing. Genetic analyses indicated that the powdery mildew resistance was controlled by a single dominant gene, temporarily designated MlWE4. By mean of comparative genomics and bulked segregant analysis, a genetic linkage map of MlWE4 was constructed, and MlWE4 was mapped on the distal region of chromosome arm 5BL. Comparative genetic linkage maps showed that genes MlWE4, Pm36 and Ml3D232 were co-segregated with markers XBD37670 and XBD37680, indicating they are likely the same gene or alleles in the same locus. The co-segregated markers provide a starting point for chromosome landing and map-based cloning of MlWE4, Pm36 and Ml3D232.