Identification of SNPs tightly linked to the QTL for pod shattering in soybean

The pod shattering or dehiscence is essential for the propagation of pod-bearing plant species in the wild, but it causes significant yield losses during harvest of domesticated crop plants. Identifying novel molecular makers, which are linked to seed-shattering genes, is needed to employ the molecular marker-assisted selection for efficiently developing shattering-resistant soybean varieties. In this study, a genetic linkage map was constructed using 115 recombinant inbred lines (RILs) developed from crosses between the pod shattering susceptible variety, Keunol, and resistant variety, Sinpaldal. A 180 K Axiom® SoyaSNPs data and pod shattering data from two environments in 2001 and 2015 were used to identify quantitative trait loci (QTL) for pod shattering. A major QTL was identified between two flanking single nucleotide polymorphism (SNP) markers, AX-90320801 and AX-90306327 on chromosome 16 with 1.3 cM interval, 857 kb of physical range. In sequence, genotype distribution analysis was conducted using extreme phenotype RILs. This could narrow down the QTL down to 153 kb on the physical map and was designated as qPDH1-KS with 6 annotated gene models. All exons within qPDH1-KS were sequenced and the 6 polymorphic SNPs affecting the amino acid sequence were identified. To develop universally available molecular markers, 38 Korean soybean cultivars were investigated by the association study using the 6 identified SNPs. Only two SNPs were strongly associated with the pod shattering. These two identified SNPs will help to identify the pod shattering responsible gene and to develop pod shattering-resistant soybean plants using marker-assisted selection.     

Impacts of soybean-induced defenses on <Emphasis Type="Italic">Spodoptera frugiperda</Emphasis> (Lepidoptera: Noctuidae) development

As soybean (Glycine max [L.] Merrill) is an important crop throughout the world, the action of herbivorous insects responsible for economic and productivity losses is the subject of constant research. Spodoptera frugiperda (J.E. Smith) caterpillars can cause extensive damage to soybean culture; this work investigated possible harmful effects on these caterpillars associated with the possible induced defenses of soybean plants. For this purpose, we assessed the biology of the insect (leaf consumption and performance traits) and chemical composition of the soybean leaves by ultraviolet–visible spectroscopy and chemometrics of three treatments: control, mechanically wounded soybean leaves, and S. frugiperda-damaged soybean leaves. The results reveal that both types of injuries induce changes in soybean metabolism regarding the production of phenolic substances, although only the herbivore-damaged plants provoke negative effects on insect biology. Variations in carotenoid production during the circadian cycle were also found in the control group. These results confirm that the soybean plants could endure and activate chemical defense mechanisms that impair the developmental lifecycle of the insect, suggesting possibilities for sustainable control strategies.     

Mapping of genetic locus for leaf trichome in <Emphasis Type="Italic">Brassica oleracea</Emphasis>

Key message

The genetic locus for leaf trichome was identified via marker-based mapping and SNP microarray assay, and a functional marker was developed to facilitate the breeding for hairiness in Brassica oleracea.

Abstract

Plant trichomes are involved in various functions particularly in protecting plants against some biotic and abiotic damages. In the present study, an F2 segregating population was developed from the cross between a glabrous cultivated B. oleracea (CC, 2n = 18) and a hairy wild relative, B. incana (CC, 2n = 18). A 1:3 segregation pattern between glabrous and hairy plants was detected among 1063 F2 genotypes, and the locus for hairiness was mapped in a 4.3-cM genetic region using 267 SSR markers among 149 F2 genotypes, corresponding to a 17.6-Mb genomic region on chromosome C01. To narrow the genetic region for hairiness, the Brassica 60 K SNP Bead Chip Arrays were applied to genotype 64 glabrous and 30 hairy F2 plants, resulting in a 1.04-Mb single peak region located in the 17.6-Mb interval. A candidate gene, BoTRY, was identified by qRT-PCR which revealed significant higher expression in glabrous F2 genotypes as compared with that in hairy plants. A cleaved amplified polymorphic site marker was successfully developed to distinguish the sequence variations of BoTRY between hairy and glabrous plants. Our study will be helpful for molecular breeding for hairiness in B. oleracea.

     

Effects of enhanced UV-B radiation on plant physiology and growth on the Tibetan Plateau: a meta-analysis

Uncertainties about the response of plant physiology and growth to enhanced UV-B radiation cause uncertainty to predict how plant production will vary under future radiation change on the Tibetan Plateau. Here, we used a meta-analysis approach to test the influence of UV-B radiation on plant physiology and growth. This hypothesis was tested by investigating the response of plants, which was expressed by some measurable variables. Enhanced UV-B radiation decreased plant biomass, plant height, basal diameter, leaf area index, maximal PSII efficiency, and Chl a+b, but increased intercellular CO₂ concentration, malondialdehyde (MDA), hydrogen peroxide, superoxide anion radical, peroxidase, ascorbate peroxidase, proline and UV-B absorbing compounds. The effect of enhanced UV-B radiation on net photosynthesis rate (P _n ) increased with mean annual precipitation and experimental duration. The effect of enhanced UV-B radiation on MDA decreased with experimental duration. The effect of enhanced UV-B radiation on superoxide dismutase (SOD) increased with the magnitude of enhanced UV-B radiation. Forests rather than grasslands exhibited a positive response of SOD and a negative response of P _n to enhanced UV-B radiation. Therefore, the effect of enhanced UV-B radiation on alpine plants varied with ecosystem types. Local climate conditions may regulate effects of enhanced UV-B radiation on alpine plants.     

Mapping QTL for leaf area in oil palm using genotyping by sequencing

Leaf area is an important parameter in oil palm breeding as it is positively correlated with oil yield. However, measurement of leaf area is tedious and also destructive. In the present study, a breeding population with 145 palms derived from a cross between Deli Dura and Avros Pisifera was used to construct a high-density linkage map and identify quantitative trait loci (QTL) for leaf area in oil palm. Using genotyping by sequencing, a linkage map containing 2413 SNPs was constructed. The total length of the linkage map was 1161.89 cM, with an average marker spacing of 0.48 cM. Based on the continuous phenotyping of leaf area from 2010 to 2015, two suggestive QTL for leaf area were mapped on chromosomes (Chr) 3 and 9. The allelic effects of the QTL associated with leaf area in the mapping population were estimated by linear regression using ordinary least squares method. The QTL on Chr 9 explained 13.3% of phenotypic variation for leaf area. A candidate gene, ARC5, within the narrow interval of QTL on Chr 9 was identified. The gene was significantly higher expressed in leaf than root and fruit of oil palm. This high-quality and SNP-based map supplies a base to fine map QTL for agronomic traits in oil palm, and the markers closely linked to the stable QTL may be used in marker-assisted selection in oil palm breeding.     

Identifying <Emphasis Type="Italic">FATB1a</Emphasis> deletion that causes reduced palmitic acid content in soybean N87-2122-4 to develop a functional marker for marker-assisted selection

Palmitic acid is a major saturated fatty acid in soybean oil, and consumption of saturated fat is linked to a risk of coronary diseases. Development of soybean (Glycine max) cultivars with reduced palmitic acid content is an important goal of soybean breeding. The FATB1a gene was previously found to be responsible for reduced palmitic acid in the soybean line N87-2122-4. The objective of this research was to characterize the FATB1a gene identified in N87-2122-4 and develop a breeder-friendly, functional marker to facilitate marker-assisted selection and improve breeding efficiency for reduced palmitic soybeans. With the availability of soybean genetic maps, reference genome, and gene annotations, an approximate 254 kb deleted genomic region, including the FATB1a gene, was identified. Based on the gene deletion information, we developed a TaqMan marker and tested it with a segregating F ₂ population that consisted of 140 individual plants derived from ‘Cook’ × N87-2122-4. The marker performed well and accounted for 57 % of the phenotypic variation. The marker was also validated using a panel of 121 diverse soybean lines with known fatty acid profiles. The result indicated that the marker can be used effectively in marker-assisted breeding for reduced palmitic acid in soybean.     

Responses of He-Ne laser on agronomic traits and the crosstalk between UVR8 signaling and phytochrome B signaling pathway in <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis> subjected to supplementary ultraviolet-B (UV-B) stress

UV-B acclimation effects and UV-B damage repair induced by a 632.8-nm He-Ne laser were investigated in Arabidopsis thaliana plants in response to supplementary UV-B stress. There was an increasing trend in growth parameters in the combination-treated plants with He-Ne laser and UV-B light compared to those stressed with enhanced UV-B light alone during different developmental stages of plants. The photosynthetic efficiency (Pn) and survival rates of seedlings were significantly higher in the combination treatments than UV-B stress alone. The expression of UVR8, phytochrome B (PhyB), and their mediated signal responsive genes such as COP1, HY5, and CHS were also significantly upregulated in plants with the laser irradiation compared with other groups without the laser. Levels of flavonol accumulation in leaves and capsule yield of He-Ne laser-treated plants were increased. The phyB-9 mutants were more sensitive to enhanced UV-B stress and had no obvious improvements in plant phenotypic development and physiological damage caused by enhanced UV-B stress after He-Ne laser irradiation. Our results suggested that UVR8 and its mediated signaling pathway via interaction with COP1 can be induced by He-Ne laser, and these processes were dependent on cytoplasmic PhyB levels in plant cells, which might be one of the most important mechanisms of He-Ne laser on UV-B protection and UV-B damage repair. These current data have also elucidated that the biostimulatory effects of He-Ne laser on Arabidopsis thaliana plants would happen not only during the early growth stage but also during the entire late developmental stage.     

Genetic mapping and QTL analysis of Botrytis resistance in <Emphasis Type="Italic">Gerbera hybrida</Emphasis>

Gerbera hybrida is an economically important cut flower. In the production and transportation of gerbera with unavoidable periods of high relative humidity, grey mould occurs and results in losses in quality and quantity of flowers. Considering the limitations of chemical use in greenhouses and the impossibility to use these chemicals in auction or after sale, breeding for resistant gerbera cultivars is considered as the best practical approach. In this study, we developed two segregating F1 populations (called S and F). Four parental linkage maps were constructed using common and parental specific SNP markers developed from expressed sequence tag sequencing. Parental genetic maps, containing 30, 29, 27 and 28 linkage groups and a consensus map covering 24 of the 25 expected chromosomes, could be constructed. After evaluation of Botrytis disease severity using three different tests, whole inflorescence, bottom (of disc florets) and ray floret, quantitative trait locus (QTL) mapping was performed using the four individual parental maps. A total of 20 QTLs (including one identical QTL for whole inflorescence and bottom tests) were identified in the parental maps of the two populations. The number of QTLs found and the explained variance of most QTLs detected reflect the complex mechanism of Botrytis disease response.     

Exogenous 5-aminolevulinic acid alleviates the detrimental effects of UV-B stress on lettuce (<Emphasis Type="Italic">Lactuca sativa</Emphasis> L) seedlings

One of the abiotic stress factors affecting plant metabolism is ultraviolet-B (UV-B) radiation. 5-Aminolevulinic acid (ALA), a key precursor of porphyrin biosynthesis, promotes plant growth and crop yields. To investigate the alleviating effects of exogenous ALA on the damages caused by UV-B exposure, two different concentrations [10 ppm (ALA1) and 25 ppm (ALA2)] of ALA were applied to lettuce seedlings for 24 h and then they were exposed to 3.3 W m^?2 UV-B. Results showed that UV-B treatment significantly decreased chlorophyll a and b (Chl a and b) concentration, enhanced the activity of antioxidant enzymes, total phenolic concentration, soluble sugar contents, expression of phenylalanine ammonia lyase (PAL) and γ-tocopherol methyltransferase (γ-TMT) genes, the concentration of malondialdehyde (MDA), hydrogen peroxide (H₂O₂), and the rate of superoxide radical (\({\text{O}}_{2}^{ - }\)) generation in the lettuce seedlings when compared to the control. Pre-treatment with exogenous ALA significantly enhanced UV-B stress tolerance in lettuce seedlings by decreasing the reactive oxygen species. On the other hand, ALA application caused more increases in the PAL and γ-TMT gene expression, antioxidant enzymes activities, Chl a and b concentration, total phenolic content, antioxidant capacity and the concentrations of soluble sugars. Obtained results indicated that UV-B radiation exerts an adverse effect on lettuce seedlings, and some of the negative effects of UV-B radiation can be alleviated by exogenous ALA.     

UV-B susceptibility and photoreactivation in embryonic development of the two-spotted spider mite, <Emphasis Type="Italic">Tetranychus urticae</Emphasis>

Developmental errors are often induced in the embryos of many organisms by environmental stress. Ultraviolet-B radiation (UV-B) is one of the most serious environmental stressors in embryonic development. Here, we investigated susceptibility to UV-B (0.5 kJ m^?2) in embryos of the two-spotted spider mite, Tetranychus urticae, to examine the potential use of UV-B in control of this important agricultural pest worldwide. Peak susceptibility to UV-B (0% hatchability) was found in T. urticae eggs 36–48 h after oviposition at 25 °C, which coincides with the stages of morphogenesis forming the germ band and initial limb primordia. However, hatchability recovered to?~?80% when eggs irradiated with UV-B were subsequently exposed to visible radiation (VIS) at 10.2 kJ m^?2, driving photoreactivation (the photoenzymatic repair of DNA damage). The recovery effect decreased to 40–70% hatchability, depending on the embryonic developmental stage, when VIS irradiation was delayed for 4 h after the end of exposure to UV-B. Thus UV-B damage to T. urticae embryos is critical, particularly in the early stages of morphogenesis, and photoreactivation functions to mitigate UV-B damage, even in the susceptible stages, but immediate VIS irradiation is needed after exposure to UV-B. These findings suggest that nighttime irradiation with UV-B can effectively kill T. urticae eggs without subsequent photoreactivation and may be useful in the physical control of this species.     

Early characterization of somatic hybrids from symmetric protoplast electrofusion of <Emphasis Type="Italic">Solanum pinnatisectum</Emphasis> Dun. and <Emphasis Type="Italic">Solanum tuberosum</Emphasis> L.

Variability of 31 somatic hybrids of Solanum pinnatisectum Dun. with Solanum tuberosum L. for leaf morphology, plant vigor, resistance to Phytophthora infestans, ploidy level, and cytoplasm type was evaluated in vitro. The composition of these somatic hybrids was as follows: [S. pinnatisectum Dun. (2n = 2x = 24; cytoplasmic type Wγ) + S. tuberosum L. (2n = 4x = 48; cytoplasmic type Tß)]. Based on leaf morphology and plant growth vigor, plants were divided into three groups, including plants close to tbr parent with unlobed leaves, small plants with scarcely dissected leaves, and vigorous plants with asymmetrically and pinnately lobed leaves. Nine of the somatic hybrids were found to be highly resistant to P. infestans. Somatic hybrids were either tetraploid or hexaploid, with hexaploids being predominant. The cytoplasm of somatic hybrids was either Tßγ or Wßγ, with Tßγ being more common. Overall, in contrast to leaf morphology and growth vigor, level of resistance to P. infestans was not related to either ploidy level or type of cytoplasm. These findings demonstrate that early in vitro selection of promising hybrids can be useful in breeding programs.     

Salt stress induced soybean <Emphasis Type="Italic">GmIFS1</Emphasis> expression and isoflavone accumulation and salt tolerance in transgenic soybean cotyledon hairy roots and tobacco

Effects of isoflavones on plant salt tolerance were investigated in soybean (Glycine max L. Merr. cultivar N23674) and tobacco (Nicotiana tabacum L.). Leaf area, fresh weight, net photosynthetic rate (Pn), and transpiration rate (Tr) of soybean N23674 plants treated with 80 mM NaCl were significantly reduced, while a gene (GmIFS1) encoding for 2-hydroxyisoflavone synthase was highly induced, and isoflavone contents significantly increased in leaves and seeds. To test the impact of isoflavones to salt tolerance, transgenic soybean cotyledon hairy roots expressing GmIFS1 (hrGmIFS1) were produced. Salt stress slightly increased isoflavone content in hairy roots of the transgenic control harboring the empty vector but substantially reduced the maximum root length, root fresh weight, and relative water content (RWC). The isoflavone content in hrGmIFS1 roots, however, was significantly higher, and the above-mentioned root growth parameters decreased much less. The GmIFS1 gene was also transformed into tobacco plants; plant height and leaf fresh weight of transgenic GmIFS1 tobacco plants were much greater than control plants after being treated with 85 mM NaCl. Leaf antioxidant capacity of transgenic tobacco was significantly higher than the control plants. Our results suggest that salt stress-induced GmIFS1 expression increased isoflavone accumulation in soybean and improved salt tolerance in transgenic soybean hairy roots and tobacco plants.     

Changes in growth and photosynthetic capacity of rice with increased UV-B radiation

Sixteen rice ( Oryza sativa L.) cultivars from 7 different geographical regions were grown in greenhouses at the Univ. of Maryland with and without supplemental ultraviolet-B (UV-B) radiation to determine alterations in biomass, morphology and maximum photosynthesis that would be anticipated from potential reductions in the stratospheric ozone column. A wide range of UV-B effects were observed, with the Philippines cultivar Carreon (5993) and the Sri Lankan cultivar Kurkaruppan (15449) showing the greatest decrease and increase, respectively, in total biomass with supplemental UV-B radiation. Approximately one-third of all cultivars tested showed a statistically significant decrease in total biomass with UV-B radiation. For these sensitive cultivars, leaf area and tiller number were also significantly reduced. Photosynthetic capacity as determined by oxygen evolution declined for some cultivars, but the correlation between changes in photosynthesis and biomass with increasing UV-B was equivocal. Results from this experiment indicate that: (1) a number of rice cultivars are sensitive to potential increases in UV-B radiation: and (2) the diversity exhibited by rice in response to increased levels of UV-B suggests that selective breeding might be successfully used to develop UV-B-tolerant rice cultivars.     

Development of F1 hybrid population and the high-density linkage map for European aspen (<Emphasis Type="Italic">Populus tremula</Emphasis> L.) using RADseq technology

Background

Restriction-site associated DNA sequencing (RADseq) technology was recently employed to identify a large number of single nucleotide polymorphisms (SNP) for linkage mapping of a North American and Eastern Asian Populus species. However, there is also the need for high-density genetic linkage maps for the European aspen (P. tremula) as a tool for further mapping of quantitative trait loci (QTLs) and marker-assisted selection of the Populus species native to Europe.

Results

We established a hybrid F1 population from the cross of two aspen parental genotypes diverged in their phenological and morphological traits. We performed RADseq of 122 F1 progenies and two parents yielding 15,732 high-quality SNPs that were successfully identified using the reference genome of P. trichocarpa. 2055 SNPs were employed for the construction of maternal and paternal linkage maps. The maternal linkage map was assembled with 1000 SNPs, containing 19 linkage groups and spanning 3054.9 cM of the genome, with an average distance of 3.05 cM between adjacent markers. The paternal map consisted of 1055 SNPs and the same number of linkage groups with a total length of 3090.56 cM and average interval distance of 2.93 cM. The linkage maps were employed for QTL mapping of one-year-old seedlings height variation. The most significant QTL (LOD = 5.73) was localized to LG5 (96.94 cM) of the male linkage map, explaining 18% of the phenotypic variation.

Conclusions

The set of 15,732 SNPs polymorphic in aspen and high-density genetic linkage maps constructed for the P. tremula intra-specific cross will provide a valuable source for QTL mapping and identification of candidate genes facilitating marker-assisted selection in European aspen.

     

Construction of an ultra-high density consensus genetic map,and enhancement of the physical map from genome sequencing in <Emphasis Type="Italic">Lupinus angustifolius</Emphasis>

Key message

An ultra-high density genetic map containing 34,574 sequence-defined markers was developed in Lupinus angustifolius. Markers closely linked to nine genes of agronomic traits were identified. A physical map was improved to cover 560.5 Mb genome sequence.

Abstract

Lupin (Lupinus angustifolius L.) is a recently domesticated legume grain crop. In this study, we applied the restriction-site associated DNA sequencing (RADseq) method to genotype an F₉ recombinant inbred line population derived from a wild type × domesticated cultivar (W × D) cross. A high density linkage map was developed based on the W × D population. By integrating sequence-defined DNA markers reported in previous mapping studies, we established an ultra-high density consensus genetic map, which contains 34,574 markers consisting of 3508 loci covering 2399 cM on 20 linkage groups. The largest gap in the entire consensus map was 4.73 cM. The high density W × D map and the consensus map were used to develop an improved physical map, which covered 560.5 Mb of genome sequence data. The ultra-high density consensus linkage map, the improved physical map and the markers linked to genes of breeding interest reported in this study provide a common tool for genome sequence assembly, structural genomics, comparative genomics, functional genomics, QTL mapping, and molecular plant breeding in lupin.

     

Genetic mapping of <Emphasis Type="Italic">psl</Emphasis> locus and quantitative trait loci for angular leaf spot resistance in cucumber (<Emphasis Type="Italic">Cucumis sativus</Emphasis> L.)

One of the most important cucumber diseases is bacterial angular leaf spot (ALS), whose increased occurrence in open-field production has been observed over the last years. To map ALS resistance genes, a recombinant inbred line (RIL) mapping population was developed from a narrow cross of cucumber line Gy14 carrying psl resistance gene and susceptible B10 line. Parental lines and RILs were tested under growth chamber conditions as well as in the field for angular leaf spot symptoms. Based on simple sequence repeat and DArTseq, genotyping a genetic map was constructed, which contained 717 loci in seven linkage groups, spanning 599.7 cM with 0.84 cM on average between markers. Monogenic inheritance of the lack of chlorotic halo around the lesions, which is typical for ALS resistance and related with the presence of recessive psl resistance gene, was confirmed. The psl locus was mapped on cucumber chromosome 5. Two major quantitative trait loci (QTL) psl5.1 and psl5.2 related to disease severity were found and located next to each other on chromosome 5; moreover, psl5.1 was co-located with psl locus. Identified QTL were validated in the field experiment. Constructed genetic map and markers linked to ALS resistance loci are novel resources that can contribute to cucumber breeding programs.     

Development of an HRM-based,safe and high-throughput genotyping system for two <Emphasis Type="Italic">low phytic acid</Emphasis> mutations in soybean

Two low phytic acid (lpa) mutants, Gm-lpa-ZC-2 (ZC-lpa) and Gm-lpa-TW-1 (TW-lpa), resulting from a G → A mutation in GmIPK1 and a 2-bp deletion in GmMIPS1, respectively, were previously developed to increase the nutritional value and environmental friendliness of soybean meal. Two functional CAPS markers were subsequently developed for genotyping plants carrying the two mutant genes; however, both are costly and time consuming and hence unsuitable for large-scale breeding use. In the present work, by integrating a quick DNA extraction protocol with an optimized high-resolution melting curve (HRM) analysis, we developed a fast and high-throughput genotyping system for the two mutations. In this system, (1) DNAs are extracted within half an hour using a protocol that only requires freezing and heating of leaf disks in two non-toxic solutions and can be directly used for PCR; (2) for genotyping, asymmetric PCRs with competitive primers are performed, and the samples are then discriminated and grouped through HRM analysis; and (3) all steps are performed in a 96-well plate, and hence adaptable to high-throughput genotyping. Although the system was developed for two lpa mutations, the general principle should be applicable to any other genes in soybean.     

Fine-mapping of QTLs for individual and total isoflavone content in soybean (<Emphasis Type="Italic">Glycine max L.</Emphasis>) using a high-density genetic map

Key message

Fifteen stable QTLs were identified using a high-density soybean genetic map across multiple environments. One major QTL, qIF5-1, contributing to total isoflavone content explained phenotypic variance 49.38, 43.27, 46.59, 45.15 and 52.50%, respectively.

Abstract

Soybeans (Glycine max L.) are a major source of dietary isoflavones. To identify novel quantitative trait loci (QTL) underlying isoflavone content, and to improve the accuracy of marker-assisted breeding in soybean, a valuable mapping population comprised of 196 F_7:8–10 recombinant inbred lines (RILs, Huachun 2 × Wayao) was utilized to evaluate individual and total isoflavone content in plants grown in four different environments in Guangdong. A high-density genetic linkage map containing 3469 recombination bin markers based on 0.2 × restriction site-associated DNA tag sequencing (RAD-seq) technology was used to finely map QTLs for both individual and total isoflavone contents. Correlation analyses showed that total isoflavone content, and that of five individual isoflavone, was significantly correlated across the four environments. Based on the high-density genetic linkage map, a total of 15 stable quantitative trait loci (QTLs) associated with isoflavone content across multiple environments were mapped onto chromosomes 02, 05, 07, 09, 10, 11, 13, 16, 17, and 19. Further, one of them, qIF5-1, localized to chromosomes 05 (38,434,171–39,045,620 bp) contributed to almost all isoflavone components across all environments, and explained 6.37–59.95% of the phenotypic variance, especially explained 49.38, 43.27, 46.59, 45.15 and 52.50% for total isoflavone. The results obtained in the present study will pave the way for a better understanding of the genetics of isoflavone accumulation and reveals the scope available for improvement of isoflavone content through marker-assisted selection.