Genome-wide insertion–deletion (InDel) marker discovery and genotyping for genomics-assisted breeding applications in chickpea

We developed 21,499 genome-wide insertion–deletion (InDel) markers (2- to 54-bp in silico fragment length polymorphism) by comparing the genomic sequences of four (desi, kabuli and wild C. reticulatum) chickpea [Cicer arietinum (L.)] accessions. InDel markers showing 2- to 6-bp fragment length polymorphism among accessions were abundant (76.8%) in the chickpea genome. The physically mapped 7,643 and 13,856 markers on eight chromosomes and unanchored scaffolds, respectively, were structurally and functionally annotated. The 4,506 coding (23% large-effect frameshift mutations) and regulatory InDel markers were identified from 3,228 genes (representing 11.7% of total 27,571 desi genes), suggesting their functional relevance for trait association/genetic mapping. High amplification (97%) and intra-specific polymorphic (60–83%) potential and wider genetic diversity (15–89%) were detected by genome-wide 6,254 InDel markers among desi, kabuli and wild accessions using even a simpler cost-effective agarose gel-based assay. This signifies added advantages of this user-friendly genetic marker system for manifold large-scale genotyping applications in laboratories with limited infrastructure and resources. Utilizing 6,254 InDel markers-based high-density (inter-marker distance: 0.212 cM) inter-specific genetic linkage map (ICC 4958 × ICC 17160) of chickpea as a reference, three major genomic regions harboring six flowering and maturity time robust QTLs (16.4–27.5% phenotypic variation explained, 8.1–11.5 logarithm of odds) were identified. Integration of genetic and physical maps at these target QTL intervals mapped on three chromosomes delineated five InDel markers-containing candidate genes tightly linked to the QTLs governing flowering and maturity time in chickpea. Taken together, our study demonstrated the practical utility of developing and high-throughput genotyping of such beneficial InDel markers at a genome-wide scale to expedite genomics-assisted breeding applications in chickpea.     

Development of a large number of SSR and InDel markers and construction of a high-density genetic map based on a RIL population of pepper (<Emphasis Type="Italic">Capsicum annuum</Emphasis> L.)

Capsicum annuum, the most widely cultivated species of pepper, is used worldwide for its important nutritional and medicinal values. The construction of an intraspecific high-density genetic linkage map would be of practical value for pepper breeding. However, the numbers of PCR-based simple sequence repeat (SSR) and insertion/deletion (InDel) markers that are available are limited, and there is a need to develop a saturated, intraspecific linkage map. The non-redundant Capsicum species’ expressed sequence tag (EST) database from the National Center for Biotechnology Information was used in this study to develop a total of 902 usable EST-SSR markers. Additionally, 177,587 SSR loci were identified based on the pepper genomic information, including 9182 SSR loci 500 bp both upstream and downstream of coding regions. Another 4497 stable and reliable InDel loci were also developed. From 9182 SSR and 4497 InDel loci, 3356 pairs of genomic SSR primers and 1400 pairs of InDel primers that were evenly distributed in 12 chromosomes were selected. A high-density intraspecific genetic map of C. annuum was constructed using the F₁₀-generation recombinant inbred line of parents PM702 and FS871 as the mapping population, screening the selected 3356 pairs of genomic SSR primers and 1400 pairs of InDel primers and the 902 EST-SSR markers developed earlier, and 524 published SSR markers and 299 orthologous markers (including 263 COSII markers and 36 tomato-derived markers) used previously to develop an interspecific genetic map (C. annuum × C. frutescens). Eventually, a high-density complete genetic intraspecific linkage map of C. annuum containing 12 linkage groups and 708 molecular markers with a length of 1260.00 cM and an average map distance of 1.78 cM was produced. This intraspecific, high-density, complete genetic linkage map of C. annuum contains the largest number of SSR and InDel markers and the highest amount of saturation so far, and it will be of considerable significance for the breeding of improved cultivars of this important field crop in the future.     

Genotyping-by-sequencing map permits identification of clubroot resistance QTLs and revision of the reference genome assembly in cabbage (Brassica oleracea L.)

Clubroot is a devastating disease caused by Plasmodiophora brassicae and results in severe losses of yield and quality in Brassica crops. Many clubroot resistance genes and markers are available in Brassica rapa but less is known in Brassica oleracea. Here, we applied the genotyping-by-sequencing (GBS) technique to construct a high-resolution genetic map and identify clubroot resistance (CR) genes. A total of 43,821 SNPs were identified using GBS data for two parental lines, one resistant and one susceptible lines to clubroot, and 18,187 of them showed >5× coverage in the GBS data. Among those, 4,103 were credibly genotyped for all 78 F₂ individual plants. These markers were clustered into nine linkage groups spanning 879.9 cM with an average interval of 1.15 cM. Quantitative trait loci (QTLs) survey based on three rounds of clubroot resistance tests using F_{2 : 3} progenies revealed two and single major QTLs for Race 2 and Race 9 of P. brassicae, respectively. The QTLs show similar locations to the previously reported CR loci for Race 4 in B. oleracea but are in different positions from any of the CR loci found in B. rapa. We utilized two reference genome sequences in this study. The high-resolution genetic map developed herein allowed us to reposition 37 and 2 misanchored scaffolds in the 02–12 and TO1000DH genome sequences, respectively. Our data also support additional positioning of two unanchored 3.3 Mb scaffolds into the 02–12 genome sequence.     

Genetic mapping and QTL analysis for yield and agronomic traits with an F2:3 population derived from a waxy corn × sweet corn cross

We constructed a framework map using SSR markers in the F₂ population derived from a cross between a waxy corn inbred line and a sweet corn inbred line. We constructed a genetic linkage map of the F_2:3 population employing 295 SSR markers on 158 F₂ individuals produced from the cross. The map comprised a total genomic length of 2,626.5 cM in 10 linkage groups and an average distance between markers of 8.9 cM. The number of loci per linkage group ranged from 27 (chr. 5) to 34 (chr. 7). The genetic distance per linkage group ranged from 213.6 cM (chr. 10) to 360.6 cM (chr. 2). Χ ² tests revealed that 254 markers (86.1 %) distributed over all 10 chromosomes exhibited a Mendelian segregation ratio of 1:2:1. A total of 14 quantitative trait loci (QTLs) for days to silking (DTS), plant height (PH), ear height (EH), ear height ratio (ER), ear length (L-ear), and setted ear length (L-sear) were found in the 158 F₂ progeny. They were mapped to chromosomes 1, 2, 3, 7, 8, and 10. Among them, one QTL was associated with DTS, three with PH, six with EH, one with ER, two with L-ear, and one QTL was related to L-sear. In our study, we found that four QTLs: qDTS1, qEH1a, qEH1b, and qPH1, were clustered between umc2390 and umc1603 on chromosome 1. These new QTLs identified by the present study could serve as useful molecular markers in selecting for yield and agronomic traits in maize. The results of this study may improve the identification and characterization of genes responsible for yield and agronomic traits in waxy corn and sweet corn.     

Construction of the First High-Density Genetic Linkage Map and Analysis of Quantitative Trait Loci for Growth-Related Traits in <Emphasis Type="Italic">Sinonovacula constricta</Emphasis>

The razor clam (Sinonovacula constricta) is an important aquaculture species, for which a high-density genetic linkage map would play an important role in marker-assisted selection (MAS). In this study, we constructed a high-density genetic map and detected quantitative trait loci (QTLs) for Sinonovacula constricta with an F₁ cross population by using the specific locus amplified fragment sequencing (SLAF-seq) method. A total of 315,553 SLAF markers out of 467.71 Mreads were developed. The final linkage map was composed of 7516 SLAFs (156.60-fold in the parents and 20.80-fold in each F₁ population on average). The total distance of the linkage map was 2383.85 cM, covering 19 linkage groups with an average inter-marker distance of 0.32 cM. The proportion of gaps less than 5.0 cM was on average 96.90%. A total of 16 suggestive QTLs for five growth-related traits (five QTLs for shell height, six QTLs for shell length, three QTLs for shell width, one QTL for total body weight, and one QTL for soft body weight) were identified. These QTLs were distributed on five linkage groups, and the regions showed overlapping on LG9 and LG13. In conclusion, the high-density genetic map and QTLs for S. constricta provide a valuable genetic resource and a basis for MAS.     

Genetic mapping and QTL analysis of agronomic traits in Indian <Emphasis Type="Italic">Mucuna pruriens</Emphasis> using an intraspecific F<Subscript>2</Subscript> population

Mucuna pruriens is a well-recognized agricultural and horticultural crop with important medicinal use. However, antinutritional factors in seed and adverse morphological characters have negatively affected its cultivation. To elucidate the genetic control of agronomic traits, an intraspecific genetic linkage map of Indian M. pruriens has been developed based on amplified fragment length polymorphism (AFLP) markers using 200 F ₂ progenies derived from a cross between wild and cultivated genotypes. The resulting linkage map comprised 129 AFLP markers dispersed over 13 linkage groups spanning a total distance of 618.88 cM with an average marker interval of 4.79 cM. For the first time, three QTLs explaining about 6.05–14.77% of the corresponding total phenotypic variation for three quantitative (seed) traits and, eight QTLs explaining about 25.96% of the corresponding total phenotypic variation for three qualitative traits have been detected on four linkage groups. The map presented here will pave a way for mapping of genes/QTLs for the important agronomic and horticultural traits contrasting between the parents used in this study.     

QTL Analysis of Head Splitting Resistance in Cabbage (Brassica oleracea L. var. capitata) Using SSR and InDel Makers Based on Whole-Genome Re-Sequencing

Head splitting resistance (HSR) in cabbage is an important trait closely related to both quality and yield of head. However, the genetic control of this trait remains unclear. In this study, a doubled haploid (DH) population derived from an intra-cross between head splitting-susceptible inbred cabbage line 79–156 and resistant line 96–100 was obtained and used to analyze inheritance and detect quantitative trait loci (QTLs) for HSR using a mixed major gene/polygene inheritance analysis and QTL mapping. HSR can be attributed to additive-epistatic effects of three major gene pairs combined with those of polygenes. Negative and significant correlations were also detected between head Hsr and head vertical diameter (Hvd), head transverse diameter (Htd) and head weight (Hw). Using the DH population, a genetic map was constructed with simple sequence repeat (SSR) and insertion–deletion (InDel) markers, with a total length of 1065.9 cM and average interval length of 4.4 cM between adjacent markers. Nine QTLs for HSR were located on chromosomes C3, C4, C7, and C9 based on 2 years of phenotypic data using both multiple-QTL mapping and inclusive composite interval mapping. The identified QTLs collectively explained 39.4 to 59.1% of phenotypic variation. Three major QTLs (Hsr 3.2, 4.2, 9.2) showing a relatively larger effect were robustly detected in different years or with different mapping methods. The HSR trait was shown to have complex genetic mechanisms. Results from QTL mapping and classical genetic analysis were consistent. The QTLs obtained in this study should be useful for molecular marker-assisted selection in cabbage breeding and provide a foundation for further research on HSR genetic regulation.     

An Improved Brassica rapa Genetic Linkage Map and Locus-specific Variations in a Doubled Haploid Population

In this study, we describe the construction of an improved Chinese cabbage genetic linkage map by integrating simple sequence repeats (SSRs) and insertion/deletion polymorphisms (InDels) into a previously published map of a doubled haploid (DH) population. The population was derived from a cross between the Chinese cabbage line BY (Brassica rapa ssp. pekinensis) and a European turnip line MM (Brassica rapa L. ssp. rapifera). A total of 629 markers were aligned to ten linkage groups, with a total map length of 1,173.8 cM, and an average distance between markers of 1.87 cm. Of the 126 SSRs and 133 InDels mapped, 46 and 34 were novel, respectively. A comparison of the linkage map with the B. rapa genome showed that more than 93 % of the markers, including 112 SSRs and 129 InDels, could be anchored unambiguously to a specific location on one of the ten chromosomes. In most cases, the order of markers on the linkage map and physical map was similar; however, the majority of linkage groups contained a number of markers whose positions were either transposed or had moved slightly forwards or backwards. During microspore culture, it was observed that 11 SSRs and one InDel showed either variation in size, or the appearance of new marker bands in the DH lines. As a first step to addressing this SSR/InDel marker instability, six SSR and one InDel loci were sequenced, which revealed that the size variation was due mainly to changes in repeat-motif number or to the insertion/deletion of new fragments of DNA.     

Theobroma cacao L.: a genetic linkage map and quantitative trait loci analysis

A genetic linkage map of Theobroma cacao (cocoa) has been constructed from 131 backcross trees derived from a cross between a single tree of the variety Catongo and an F1 tree from the cross of Catongo by Pound 12. The map comprises 138 markers: 104 RAPD loci, 32 RFLP loci and two morphologic loci. Ten linkage groups were found which cover 1068 centimorgans (cM). Only six (4%) molecular-marker loci show a significant deviation from the expected 11 segregation ratio.The average distance between two adjacent markers is 8.3 cM. The final genome-size estimates based on two-point linkage data ranged from 1078 to 1112 cM for the cocoa genome. This backcross progeny segregates for two apparently single gene loci controlling (1) anthocyanidin synthesis (Anth) in seeds, leaves and flowers and (2) self-compatibility (Autoc). The Anth locus was found to be 25 cM from Autoc and two molecular markers co-segregate with Anth. The genetic linkage map was used to localize QTLs for early flowering, trunk diameter, jorquette height and ovule number in the BC₁ generation using both single-point ANOVA and interval mapping. A minimum number of 2–4 QTLs (P<0.01) involved in the genetic expression of the traits studied was detected. Coincident map locations of a QTL for jorquette height and trunk diameter suggests the possibility of pleiotropic effects in cocoa for these traits. The combined estimated effects of the different mapped QTLs explained between 11.2% and 25.8% of the phenotypic variance observed in the BC₁ population.     

Construction of a genetic linkage map and identification of molecular markers in peach rootstocks for response to peach tree short life syndrome

Peach tree short life (PTSL) is a devastating disease syndrome of peach [Prunus persica (L.) Batsch] caused by multiple factors; the molecular biology of its tolerance/susceptibility is unknown. The difficulty of studying PTSL is that tree survival or death is not obvious until 3 to 5 years after planting when the symptoms of PTSL first appear. Tolerance to PTSL was unknown in Prunus until the rootstock Guardian® ‘BY520-9’ was introduced into commercial orchards in 1994. To study the genetics of the response to PTSL, a controlled F₂ cross was made between Guardian® ‘BY520-9’ selection 3-17-7 (PTSL-tolerant) and Nemaguard (PTSL-susceptible). An F₁ hybrid was then selfed to generate an F₂ population expected to segregate for PTSL response. One hundred fifty-one AFLPs and 21 SSRs, including anchor loci from the Prunus reference genetic map, were used to construct a molecular genetic map based on 100 F₂ seedlings. This map covers a genetic distance of 737 cM with an average marker spacing of 4.7 cM and will be used as a framework to construct a highly saturated molecular genetic map. Of the 140 mapped AFLP markers, 38 were associated with PTSL response, as identified previously by bulked segregant analysis. The distribution of the markers associated with PTSL response on the newly constructed genetic map was compared with the recently published Prunus resistance map. This comparison revealed that some resistance gene analogs and several PTSL-associated AFLP markers were located in the same regions in several Prunus linkage groups: G1, G2, G4, G5, and G6. This peach rootstock map can also be viewed and compared with other Prunus maps in comparative map viewer CMap in the Genome Database for Rosaceae (GDR) at http://www.rosaceae.org     

The first genetic and comparative map of white lupin (Lupinus albus L.): identification of QTLs for anthracnose resistance and flowering time, and a locus for alkaloid content.

We report the first genetic linkage map of white lupin (Lupinus albus L.). An F8 recombinant inbred line population developed from Kiev mutant x P27174 was mapped with 220 amplified fragment length polymorphism and 105 gene-based markers. The genetic map consists of 28 main linkage groups (LGs) that varied in length from 22.7 cM to 246.5 cM and spanned a total length of 2951 cM. There were seven additional pairs and 15 unlinked markers, and 12.8% of markers showed segregation distortion at P < 0.05. Syntenic relationships between Medicago truncatula and L. albus were complex. Forty-five orthologous markers that mapped between M. truncatula and L. albus identified 17 small syntenic blocks, and each M. truncatula chromosome aligned to between one and six syntenic blocks in L. albus. Genetic mapping of three important traits: anthracnose resistance, flowering time, and alkaloid content allowed loci governing these traits to be defined. Two quantitative trait loci (QTLs) with significant effects were identified for anthracnose resistance on LG4 and LG17, and two QTLs were detected for flowering time on the top of LG1 and LG3. Alkaloid content was mapped as a Mendelian trait to LG11.     

A haplotype inference algorithm for trios based on deterministic sampling

Background

The three-dimensional shape of grain, measured as grain length, width, and thickness (GL, GW, and GT), is one of the most important components of grain appearance in rice. Determining the genetic basis of variations in grain shape could facilitate efficient improvements in grain appearance. In this study, an F_7:8 recombinant inbred line population (RIL) derived from a cross between indica and japonica cultivars (Nanyangzhan and Chuan7) contrasting in grain size was used for quantitative trait locus (QTL) mapping. A genetic linkage map was constructed with 164 simple sequence repeat (SSR) markers. The major aim of this study was to detect a QTL for grain shape and to fine map a minor QTL, qGL7.

Results

Four QTLs for GL were detected on chromosomes 3 and 7, and 10 QTLs for GW and 9 QTLs for GT were identified on chromosomes 2, 3, 5, 7, 9 and 10, respectively. A total of 28 QTLs were identified, of which several are reported for the first time; four major QTLs and six minor QTLs for grain shape were also commonly detected in both years. The minor QTL, qGL7, exhibited pleiotropic effects on GL, GW, GT, 1000-grain weight (TGW), and spikelets per panicle (SPP) and was further validated in a near isogenic F₂ population (NIL-F₂). Finally, qGL7 was narrowed down to an interval between InDel marker RID711 and SSR marker RM6389, covering a 258-kb region in the Nipponbare genome, and cosegregated with InDel markers RID710 and RID76.

Conclusion

Materials with very different phenotypes were used to develop mapping populations to detect QTLs because of their complex genetic background. Progeny tests proved that the minor QTL, qGL7, could display a single mendelian characteristic. Therefore, we suggested that minor QTLs for traits with high heritability could be isolated using a map-based cloning strategy in a large NIL-F₂ population. In addition, combinations of different QTLs produced diverse grain shapes, which provide the ability to breed more varieties of rice to satisfy consumer preferences.     

Analysis of genetic mapping in a waxy/dent maize RIL population using SSR and SNP markers

A maize genetic linkage map was generated using SSR and SNP markers in a F_7:8 recombinant inbred line (RIL) population derived from a cross of waxy corn (KW7) and dent corn (Mo17). A total of 465 markers, including 459 SSR and 6 SNP markers, were assigned to 10 linkage groups which spanned 2,656.5 cM with an average genetic distance between markers of 5.7 cM, and the number of loci per linkage group ranged from 39 to 55. The SSR (85.4%) and SNP (83.3%) markers showed Mendelian segregation ratios in the RIL population at a 5% significance threshold. In linkage analysis of six SNP loci associated with kernel starch synthesis genes (ae1, bt2, sh1, sh2, su1, and wx1), all six loci were successfully mapped and are closely linked with SSR markers in chromosomes 3 (sh2), 4 (su1 and bt2), 5 (ae1), and 9 (sh1 and wx1). The SSR markers linked with genes in starch synthesis may be utilized in marker assisted breeding programs. The resulting genetic map will be useful in dissection of quantitative traits and the identification of superior QTLs from the waxy hybrid corn. Additionally, these data support further genetic analysis and development of maize breeding programs.     

SSR-based genetic linkage map construction in pistachio using an interspecific F<Subscript>1</Subscript> population and QTL analysis for leaf and shoot traits

Pistachio is one of the most commercially important nut trees in the world. To characterize the genetic controls of horticultural traits and facilitate marker-assisted breeding in pistachio, we constructed an SSR-based linkage map using an interspecific F₁ population derived from a cross between the cultivar “Siirt” (Pistacia vera L.) and the monoecious Pa-18 genotype of Pistacia atlantica Desf. This population was also used for the first QTL analysis in pistachio on leaf and shoot characters. In total, 1312 SSR primers were screened, and 388 loci were successfully integrated into parental linkage maps. The Siirt maternal map contained 306 markers, while the “Pa-18” paternal map included 285 markers along the 15 linkage groups. The Siirt map spanned 1410.4 cM, with an average marker distance of 4.6 cM; the Pa-18 map covered 1362.5 cM with an average marker distance of 4.8 cM. Phenotypic data were collected during the growing seasons of 2015 and 2016 for four traits: leaf length (LL), leaf width (LW), leaf length/leaf width ratio (LWR), number of leaflet pairs (NLL), and young shoot color (YSC). A total of 17 QTLs were identified in the parental maps. Four QTLs for LL and LW were located on LG2 and LG4, while four QTLs for LWR ratio on LG13 and LG14, two QTLs for NLL and two QTLs for YSC were on LG7 and LG9, respectively, with similar positions in both parental maps. The SSR markers, linkage maps, and QTLs reported here will provide a valuable resource for future molecular and genetic studies in pistachio.     

A naturally occurring insertion in the <Emphasis Type="Italic">RsFLC2</Emphasis> gene associated with late-bolting trait in radish (<Emphasis Type="Italic">Raphanus sativus</Emphasis> L.)

Premature flowering reduces the yield and quality of the harvested fleshy taproot in radish. However, there has been little molecular marker research on the radish late-bolting trait. In this study, F₂ and F_2:3 populations derived from a cross of “Ninengo” (late-bolting) and “Maer” (early-bolting) were analyzed to map late-bolting genes. Five hundred insertion and deletion (InDel) markers were designed according to the whole-genome resequencing data of the two parents. A genetic map was constructed based on the F₂ population, and a late-bolting gene was detected in a 1.1-cM region between the markers InDel520 and InDel535 on chromosome R02 that explained the highest (76.4%) phenotypic variance. RsFLC2 was identified as a candidate gene in this region. Notably, “Ninengo” contains a 1627-bp insertion near the 5′ end of the first intron of RsFLC2. Allelic variation analyses in the F₂ population further validated that RsFLC2was associated with the late-bolting trait in radish. The expression pattern of RsFLC2 was significantly different between “Ninengo” and “Maer” during vernalization. Vernalization suppressed RsFLC2 expression, and the 1627-bp insertion in the first intron weakened gene repression in “Ninengo” plants, resulting in late-bolting. This study lays a foundation for uncovering the molecular mechanism of late-bolting and marker-assisted selection for breeding late-bolting varieties of radish.     

A SNP-based genetic linkage map of <Emphasis Type="Italic">Capsicum baccatum</Emphasis> and its comparison to the <Emphasis Type="Italic">Capsicum annuum</Emphasis> reference physical map

Capsicum baccatum L., one of five domesticated species of Capsicum, is a valuable species in chili pepper breeding. In particular, it is a source of disease resistance against anthracnose and powdery mildew. Genetic maps and molecular markers are important to improve the efficiency of crop breeding programs. Recently, using genetic maps several researchers have identified quantitative trait loci (QTLs) for important horticultural traits and have cloned genes of interest. In this study, we constructed a genetic map of C. baccatum in an intraspecific population from a cross between ‘Golden-aji’ and ‘PI594137.’ A total of 395 high-resolution melting markers were developed based on single-nucleotide polymorphisms identified by comparing genome sequences generated through next-generation resequencing of the parents, ‘Golden-aji’ and ‘PI594137.’ The genetic linkage map contained 12 linkage groups, covered a total distance of 1056.2 cM, and had an average distance of 2.67 cM between markers. In addition, the final map was compared to the reference physical map of C. annuum ‘CM334.’ Interestingly, two major reciprocal translocations between chromosomes 3 and 5 and between chromosomes 3 and 9 were found, suggesting that these translocations might act as a genetic barrier between C. annuum and C. baccatum. Translocations between chromosomes 1 and 8 were also observed, as were previously reported in C. chinense, C. frutescens, and wild C. annuum. The synteny of other chromosomes was maintained, on the whole, except for several small inversions. The information on this genetic map will be helpful to analyze QTLs for important traits such as anthracnose resistance in C. baccatum and to study the causes of genetic barriers between C. annuum and C. baccatum.     

Mapping of QTLs for morpho-agronomic and seed quality traits in a RIL population of common bean (Phaseolus vulgaris L.)

The objective of this research was to determine the quantitative trait loci (QTLs) controlling phenological traits (days to flowering, days to end of flowering, days to harvest as green pod, and days to maturity), seed size traits (seed length, seed height, seed width, and seed weight), and seed quality traits (water absorption, and coat proportion), in common bean. A population of 104 F₇ recombinant inbred lines (RILs) derived from an inter-gene pool cross between Xana, and Cornell 49242, was used to develop a genetic linkage map including 175 AFLPs, 27 microsatellites, 30 SCARs, 33 ISSRs, 12 RAPDs, 13 loci codifying for seed proteins, and the four genes Fin,fin (growth habit); Asp,asp (seed coat shininess); P,p (seed color); and I,i (resistance to bean common mosaic virus). The map has a total length of 1,042 cM distributed across 11 linkage groups aligned to those of the core linkage map of bean using common molecular markers as anchor points. The QTL analyses were carried out over three environments using the mean environment data with composite interval mapping. Thirty-one QTLs for ten traits were found to be significant in at least one environment and in the mean environment data, the number of significant QTLs identified per trait ranging from two to five. Twenty-seven of these QTLs mapped forming clusters in eight different chromosomal regions. The rationale for this clustered mapping and the possible relationship between some QTLs for phenological traits and the genes Fin and I are discussed.     

Quantitative trait locus mapping under irrigated and drought treatments based on a novel genetic linkage map in mungbean (<Emphasis Type="Italic">Vigna radiata</Emphasis> L.)

Key message

A novel genetic linkage map was constructed using SSR markers and stable QTLs were identified for six drought tolerance related-traits using single-environment analysis under irrigation and drought treatments.

Abstract

Mungbean (Vigna radiata L.) is one of the most important leguminous food crops. However, mungbean production is seriously constrained by drought. Isolation of drought-responsive genetic elements and marker-assisted selection breeding will benefit from the detection of quantitative trait locus (QTLs) for traits related to drought tolerance. In this study, we developed a full-coverage genetic linkage map based on simple sequence repeat (SSR) markers using a recombinant inbred line (RIL) population derived from an intra-specific cross between two drought-resistant varieties. This novel map was anchored with 313 markers. The total map length was 1010.18 cM across 11 linkage groups, covering the entire genome of mungbean with a saturation of one marker every 3.23 cM. We subsequently detected 58 QTLs for plant height (PH), maximum leaf area (MLA), biomass (BM), relative water content, days to first flowering, and seed yield (Yield) and 5 for the drought tolerance index of 3 traits in irrigated and drought environments at 2 locations. Thirty-eight of these QTLs were consistently detected two or more times at similar linkage positions. Notably, qPH5A and qMLA2A were consistently identified in marker intervals from GMES5773 to MUS128 in LG05 and from Mchr11-34 to the HAAS_VR_1812 region in LG02 in four environments, contributing 6.40–20.06% and 6.97–7.94% of the observed phenotypic variation, respectively. None of these QTLs shared loci with previously identified drought-related loci from mungbean. The results of these analyses might facilitate the isolation of drought-related genes and help to clarify the mechanism of drought tolerance in mungbean.