Rapid growth of Moso bamboo (Phyllostachys edulis): Cellular roadmaps, transcriptome dynamics, and environmental factors

Moso bamboo (Phyllostachys edulis) shows remarkably rapid growth (114.5 cm/day), but the underlying biological mechanisms remain unclear. After examining more than 12,750 internodes from more than 510 culms from 17 Moso populations, we identified internode 18 as a representative internode for rapid growth. This internode includes a 2-cm cell division zone (DZ), a cell elongation zone up to 12 cm, and a secondary cell wall (SCW) thickening zone. These zones elongated 11.8 cm, produced approximately 570,000,000 cells, and deposited ∼28 mg g-1 dry weight (DW) lignin and ∼44 mg g-1 DW cellulose daily, far exceeding vegetative growth observed in other plants. We used anatomical, mathematical, physiological, and genomic data to characterize development and transcriptional networks during rapid growth in internode 18. Our results suggest that (1) gibberellin may directly trigger the rapid growth of Moso shoots, (2) decreased cytokinin and increased auxin accumulation may trigger cell DZ elongation, and (3) abscisic acid and mechanical pressure may stimulate rapid SCW thickening via MYB83L. We conclude that internode length involves a possible tradeoff mediated by mechanical pressure caused by rapid growth, possibly influenced by environmental temperature and regulated by genes related to cell division and elongation. Our results provide insight into the rapid growth of Moso bamboo.

Molecular mapping of QTLs for grain dimension traits in Basmati rice

Basmati rice is known for its extra-long slender grains, exceptional kernel dimensions after cooking, high volume expansion, and strong aroma. Developing high yielding Basmati rice varieties with good cooking quality is a gigantic task. Therefore, identifying the genomic regions governing the grain and cooked kernel dimension traits is of utmost importance for its use in marker-assisted breeding. Although several QTLs governing grain dimension traits have been reported, limited attempts have been made to map QTLs for grain and cooked kernel dimension traits of Basmati rice. In the current study, a population of recombinant inbred lines (RIL) was generated from a cross of Sonasal and Pusa Basmati 1121 (PB1121). In the RIL population, there was a significant positive correlation among the length (RRL: rough rice length, MRL: milled rice length, CKL: cooked kernel length) and breadth (RRB: rough rice breadth, MRB: milled rice breadth and CKB: cooked kernel breadth) of the related traits, while there was significant negative correlation between them. QTL mapping has led to the identification of four major genomic regions governing MRL and CKL. Two QTLs co-localize with the earlier reported major gene GS3 and a QTL qGRL7.1, while the remaining two QTLs viz., qCKL3.2 (qMRL3.2) and qCKL4.1 (qMRL4.1) were novel. The QTL qCKL3.2 has been bracketed to a genomic region of 0.78 Mb between the markers RM15247 and RM15281. Annotation of this region identified 18 gene models, of which the genes predicted to encode pentatricopeptides and brassinosteroid insensitive 1-associated receptor kinase 1 precursor may be the putative candidate genes. Furthermore, we identified a novel QTL qKER2.1 governing kernel elongation ratio (KER) in Basmati rice.

Genome-Wide Association Mapping Reveals Novel Putative Gene Candidates Governing Reproductive Stage Heat Stress Tolerance in Rice

Temperature rise predicted for the future will severely affect rice productivity because the crop is highly sensitive to heat stress at the reproductive stage. Breeding tolerant varieties is an economically viable option to combat heat stress, for which the knowledge of target genomic regions associated with the reproductive stage heat stress tolerance (RSHT) is essential. A set of 192 rice genotypes of diverse origins were evaluated under natural field conditions through staggered sowings for RSHT using two surrogate traits, spikelet fertility and grain yield, which showed significant reduction under heat stress. These genotypes were genotyped using a 50 k SNP array, and the association analysis identified 10 quantitative trait nucleotides (QTNs) for grain yield, of which one QTN (qHTGY8.1) was consistent across the different models used. Only two out of 10 MTAs coincided with the previously reported QTLs, making the remaing eight novel. A total of 22 QTNs were observed for spikelet fertility, among which qHTSF5.1 was consistently found across three models. Of the QTNs identified, seven coincided with previous reports, while the remaining QTNs were new. The genes near the QTNs were found associated with the protein–protein interaction, protein ubiquitination, stress signal transduction, and so forth, qualifying them to be putative for RSHT. An in silico expression analysis revealed the predominant expression of genes identified for spikelet fertility in reproductive organs. Further validation of the biological relevance of QTNs in conferring heat stress tolerance will enable their utilization in improving the reproductive stage heat stress tolerance in rice.

Transcriptome analysis of a near-isogenic line and its recurrent parent reveals the role of Pup1 QTL in phosphorus deficiency tolerance of rice at tillering stage

Phosphorus (P) is essential for cellular processes like respiration, photosynthesis, biosynthesis of membrane phospholipids, etc. To cope with P deficiency stress, plants adopt reprograming of the expression of genes involved in different metabolic/signaling pathways for survival, growth, and development. Plants use transcriptional, post-transcriptional, and/or post-translational machinery to achieve P homeostasis. Several transcription factors (TFs), miRNAs, and P transporters play important roles in P deficiency tolerance; however, the underlying mechanisms responsible for P deficiency tolerance remain poorly understood. Studies on P starvation/deficiency responses in plants at early (seedling) stage of growth have been reported but only a few of them focused on molecular responses of the plant at advanced (tillering or reproductive) stage of growth. To decipher the strategies adopted by rice at tillering stage under P deficiency stress, a pair of contrasting genotypes [Pusa-44 (a high-yielding, P deficiency sensitive cultivar) and its near-isogenic line (NIL-23, P deficiency tolerant) for Pup1 QTL] was used for morphophysiological, biochemical, and molecular analyses. Comparative analyses of shoot and root tissues from 45-day-old plants grown hydroponically under P sufficient (16 ppm) or P deficient (4 ppm) medium confirmed some of the known morphophysiological responses. Moreover, RNA-seq analysis revealed the important roles of phosphate transporters, TFs, auxin-responsive proteins, modulation in the cell wall, fatty acid metabolism, and chromatin architecture/epigenetic modifications in providing P deficiency tolerance to NIL-23, which were brought in due to the introgression of the Pup1 QTL in Pusa-44. This study provides insights into the molecular functions of Pup1 for P deficiency tolerance, which might be utilized to improve P-use efficiency of rice for better productivity in P deficient soils. KEY MESSAGE: Introgression of Pup1 QTL in high-yielding rice cultivar modulates mainly phosphate transporters, TFs, auxin-responsive proteins, cell wall structure, fatty acid metabolism, and chromatin architecture/epigenetic modifications at tillering stage of growth under phosphorus deficiency stress.

Cellular and molecular characterizations of the irregular internode division zone formation of a slow-growing bamboo variant

The key molecular mechanisms underlying the sectionalized growth within bamboo or other grass internodes remain largely unknown. Here, we genetically and morphologically compared the culm and rhizome internode division zones (DZs) of a slow-growing bamboo variant (sgv) having dwarf internodes, with those of the corresponding wild type (WT). Histological analysis discovers that the sgv has an irregular internode DZ. However, the shoot apical meristems in height, width, outside shape, cell number and cell width of the sgv and the WT were all similar. The DZ irregularities first appeared post apical meristem development, in 1-mm sgv rhizome internodes. Thus, the sgv is a DZ irregularity bamboo variant, which has been first reported in bamboo according to our investigation. Transcriptome sequencing analysis finds that a number of cell wall biogenesis and cell division-related genes are dramatically downregulated in the sgv DZ. Interestingly, both transcriptomic and brassinosteroid (BR) contents detection, as well as quantitative real-time PCR analyses show that these irregularities have resulted from the BR signaling pathway defects. Brassinosteroid defect might also cause the erect leaves and branches as well as the irregular epidermis of the sgv. These results suggest that BR signaling pathway plays critical roles in bamboo internode DZ and leaf development from a mutant perspective and also explain the upstream mechanisms causing the dwarf internode of the sgv bamboo.

Identification of potential MTAs and candidate genes for juice quality- and yield-related traits in Saccharum clones: a genome-wide association and comparative genomic study

Sugarcane is an economically important commercial crop which provides raw material for the production of sugar, jaggery, bioethanol, biomass and other by-products. Sugarcane breeding till today heavily relies on conventional breeding approaches which is time consuming, laborious and costly. Integration of marker-assisted selection (MAS) in sugarcane genetic improvement programs for difficult to select traits like sucrose content, resistance to pests and diseases and tolerance to abiotic stresses will accelerate varietal development. In the present study, association mapping approach was used to identify QTLs and genes associated with sucrose and other important yield-contributing traits. A mapping panel of 110 diverse sugarcane genotypes and 148 microsatellite primers were used for structured association mapping study. An optimal subpopulation number (ΔK) of 5 was identified by structure analysis. GWAS analysis using TASSEL identified a total of 110 MTAs which were localized into 27 QTLs by GLM and MLM (Q + K, PC + K) approaches. Among the 24 QTLs sequenced, 12 were able to identify potential candidate genes, viz., starch branching enzyme, starch synthase 4, sugar transporters and G3P-DH related to carbohydrate metabolism and hormone pathway-related genes ethylene insensitive 3-like 1, reversion to ethylene sensitive1-like, and auxin response factor associated to juice quality- and yield-related traits. Six markers, NKS 5_185, SCB 270_144, SCB 370_256, NKS 46_176 and UGSM 648_245, associated with juice quality traits and marker SMC31CUQ_304 associated with NMC were validated and identified as significantly associated to the traits by one-way ANOVA analysis. In conclusion, 24 potential QTLs identified in the present study could be used in sugarcane breeding programs after further validation in larger population. The candidate genes from carbohydrate and hormone response pathway presented in this study could be manipulated with genome editing approaches to further improve sugarcane crop.

A novel LOX3-null allele (lox3-b) originated in the aromatic Basmati rice cultivars imparts storage stability to rice bran

Rapid deterioration of rice bran due to the LOX3 enzyme catalysed oxidation of PUFA is the major bottleneck for its utilization in various downstream applications. In the present study, we have identified a set of nine novel LOX3-null rice accessions carrying a deletion of C residue in the exon2 causing a frameshift mutation resulting in a truncated non-functional LOX3 protein. Our study, further manifested the predominance of C deletion based LOX3-null allele, named lox3-b, in the aromatic rice germplasm particularly in the Indian Basmati rice group. The LOX3-null genotypes exhibited significantly reduced rancidity, after six months of storage. They also showed significantly lower percentage reduction of linoleic acid (LA), higher γ-oryzanol content and lower hexanal content. A functional dCAPS marker designed based on the deletion polymorphism clearly differentiated LOX3 and lox3-b alleles, and has the potential application in marker assisted rice breeding programmes to develop cultivars with better bran storability.

Characterization of contrasting rice (Oryza sativa L.) genotypes reveals the Pi-efficient schema for phosphate starvation tolerance

BACKGROUND

Phosphorus (P), being one of the essential components of nucleic acids, cell membranes and enzymes, indispensable for diverse cellular processes like photosynthesis/carbohydrate metabolism, energy production, redox homeostasis and signaling. Crop yield is severely affected due to Phosphate (Pi) deficiency; and to cope with Pi-deficiency, plants have evolved several strategies. Some rice genotypes are compatible with low Pi availability, whereas others are sensitive to Pi deficiency. However, the underlying molecular mechanism for low Pi tolerance remains largely unexplored.

RESULT

Several studies were carried out to understand Pi-deficiency responses in rice at seedling stage, but few of them targeted molecular aspects/responses of Pi-starvation at the advanced stage of growth. To delineate the molecular mechanisms for low Pi tolerance, a pair of contrasting rice (Oryza sativa L.) genotypes [viz. Pusa-44 (Pi-deficiency sensitive) and its near isogenic line (NIL-23, Pi-deficiency tolerant) harboring Phosphorus uptake 1 (Pup1) QTL from an aus landrace Kasalath] were used. Comparative morphological, physiological, and biochemical analyses confirmed some of the well-known findings. Transcriptome analysis of shoot and root tissues from 45-day-old rice plants grown hydroponically under P-sufficient (16 ppm Pi) or P-starved (0 ppm Pi) medium revealed that Pi-starvation stress causes global transcriptional reprogramming affecting several transcription factors, signaling pathways and other regulatory genes. We could identify several significantly up-regulated genes in roots of NIL-23 under Pi-starvation which might be responsible for the Pi starvation tolerance. Pathway enrichment analysis indicated significant role of certain phosphatases, transporters, transcription factors, carbohydrate metabolism, hormone-signaling, and epigenetic processes in improving P-starvation stress tolerance in NIL-23.

CONCLUSION

We report the important candidate mechanisms for Pi acquisition/solubilization, recycling, remobilization/transport, sensing/signalling, genetic/epigenetic regulation, and cell wall structural changes to be responsible for P-starvation tolerance in NIL-23. The study provides some of the novel information useful for improving phosphorus-use efficiency in rice cultivars.

Marker Assisted Development and Characterization of Herbicide Tolerant Near Isogenic Lines of a Mega Basmati Rice Variety, "Pusa Basmati 1121"

BACKGROUND

Direct-seeded rice (DSR) is a potential technology for sustainable rice farming as it saves water and labor. However, higher incidence of weed under DSR limits productivity. Therefore, there is a need to develop herbicide tolerant (HT) rice varieties.

RESULTS

We used marker assisted backcross breeding (MABB) to transfer a mutant allele of Acetohydroxy acid synthase (AHAS) gene, which confers tolerance to imidazolinone group of herbicides from the donor parent (DP) "Robin" into the genetic background of an elite popular Basmati rice variety, Pusa Basmati 1121 (PB 1121). Foreground selection was done using the AHAS gene linked Simple Sequence Repeat (SSR) marker RM6844 and background selection was performed using 112 genome-wide SSR markers polymorphic between PB 1121 and Robin. Phenotypic selection for agronomic, Basmati grain and cooking quality traits in each generation was carried out to improve the recovery of recurrent parent phenome (RPP). Finally, a set of 12 BC₄F₄ near isogenic lines (NILs), with recurrent parent genome (RPG) recovery ranging from 98.66 to 99.55% were developed and evaluated. PB 1121-HT NILs namely 1979-14-7-33-99-10, 1979-14-7-33-99-15 and 1979-14-7-33-99-66 were found superior to PB 1121 in yield with comparable grain and cooking quality traits and herbicide tolerance similar to Robin.

CONCLUSION

Overall, the present study reports successful development of HT NILs in the genetic background of popular Basmati rice variety, PB 1121 by introgression of mutated AHAS allele. This is the first report on the development of HT Basmati rice. Superior NILs are being evaluated in the national Basmati trials, the release of which is likely to provide a viable option for the adoption of DSR technology in Basmati rice cultivation.

Marker aided introgression of 'Saltol', a major QTL for seedling stage salinity tolerance into an elite Basmati rice variety 'Pusa Basmati 1509'

Marker assisted backcross breeding was used to transfer Saltol, a major QTL for seedling stage salinity tolerance from the donor FL478 to Pusa Basmati 1509 (PB 1509), a high yielding and early maturing Basmati rice variety. Foreground selection was carried out using three markers namely, AP3206f, RM3412b and RM10793, linked to Saltol. In addition, 105 genome-wide SSR markers polymorphic between FL478 and PB 1509 were used in background selection. Among the BC₃F₄ near isogenic lines (NILs) developed, recurrent parent genome recovery ranged from 96.67 to 98.57%. Multi-season evaluation identified some of the NILs showing significantly higher yield with grain and cooking quality comparable to PB 1509. All the NILs exhibited tolerance to salinity with significantly higher relative water content, membrane stability index and proline content as compared to PB 1509. The root and shoot concentration of Na⁺, K⁺ and Na⁺/K⁺ in NILs was at par with FL478 under stress conditions. The gene OsHKT1;5 located in the Saltol region showed higher expression levels under stress indicating its role in conferring salinity tolerance. Salt tolerant NILs of PB 1509 will be useful in stabilizing production in salt affected areas.

Genome-Wide Association Study Reveals Novel Marker-Trait Associations (MTAs) Governing the Localization of Fe and Zn in the Rice Grain

Micronutrient malnutrition due to Fe and Zn, affects around two billion people globally particularly in the developing countries. More than 90% of the Asian population is dependent on rice-based diets, which is low in these micronutrients. In the present study, a set of 192 Indian rice germplasm accessions, grown at two locations, were evaluated for Fe and Zn in brown rice (BR) and milled rice (MR). A significant variation was observed in the rice germplasm for these micronutrients. The grain Fe concentration was in the range of 6.2–23.1 ppm in BR and 0.8–12.3 ppm in MR, while grain Zn concentration was found to be in the range of 11.0–47.0 ppm and 8.2–40.8 ppm in the BR and MR, respectively. Grain Fe exhibited maximum loss upon milling with a mean retention of 24.9% in MR, while Zn showed a greater mean retention of 74.2% in MR. A genome-wide association study (GWAS) was carried out implementing the FarmCPU model to control the population structure and kinship, and resulted in the identification of 29 marker-trait associations (MTAs) with significant associations for traits viz. FeBR (6 MTAs), FeMR (7 MTAs), ZnBR (11 MTAs), and ZnMR (5 MTAs), which could explain the phenotypic variance from 2.1 to as high as 53.3%. The MTAs governing the correlated traits showed co-localization, signifying the possibility of their simultaneous improvement. The robust MTAs identified in the study could be valuable resource for enhancing Fe and Zn concentration in the rice grain and addressing the problem of Fe and Zn malnutrition among rice consumers.

WA-CMS-based iso-cytoplasmic restorers derived from commercial rice hybrids reveal distinct population structure and genetic divergence towards restorer diversification

One hundred diverse iso-cytoplasmic restorer (ICR) lines carrying WA cytoplasm indicated significant but moderate variability for agro-morphological traits as well as for the microsatellite-based allele patterns. There were two major groups of ICRs based on agro-morphological clustering. Simple sequence repeat (SSR) markers identified allelic variants with an average of 2.48 alleles per locus and the gene diversity (GD) ranged from 0.02 to 0.62 at different loci. ICR lines showed a genetic structure involving two sub-populations, POP1 and POP2. Both the subpopulations had the presence of admixture lines. Nearest ancestry-based grouping of ICRs by neighbour-joining (NJ) method showed near similar grouping as that of sub-population division. The POP2 was the largest group but with fewer admixed lines. POP1 was more distinct than POP2. Since the hybrid parents of the ICRs had limited diversity on maternal lineage, paternal lineage was concluded as the major contributor to the observed divergence and population differentiation. ICRs developed from certain hybrids were more genetically distinct than other hybrids. Even with the moderate variability, ICRs could be considered as a potential source of fertility restoration in hybrid development because of their distinct population structure and the full complement of restorer genes they contained. ICR lines with high per se performance can be utilized in hybrid rice development by estimating their combining ability.

Marker-aided Incorporation of Xa38, a Novel Bacterial Blight Resistance Gene, in PB1121 and Comparison of its Resistance Spectrum with xa13 + Xa21

Basmati rice is preferred internationally because of its appealing taste, mouth feel and aroma. Pusa Basmati 1121 (PB1121) is a widely grown variety known for its excellent grain and cooking quality in the international and domestic market. It contributes approximately USD 3 billion to India's forex earning annually by being the most traded variety. However, PB1121 is highly susceptible to bacterial blight (BB) disease. A novel BB resistance gene Xa38 was incorporated in PB1121 from donor parent PR114-Xa38 using a modified marker-assisted backcross breeding (MABB) scheme. Phenotypic selection prior to background selection was instrumental in identifying the novel recombinants with maximum recovery of recurrent parent phenome. The strategy was effective in delimiting the linkage drag to <0.5 mb upstream and <1.9 mb downstream of Xa38 with recurrent parent genome recovery upto 96.9% in the developed NILs. The NILs of PB1121 carrying Xa38 were compared with PB1121 NILs carrying xa13 + Xa21 (developed earlier in our lab) for their resistance to BB. Both NILs showed resistance against the Xoo races 1, 2, 3 and 6. Additionally, Xa38 also resisted Xoo race 5 to which xa13 + Xa21 was susceptible. The PB1121 NILs carrying Xa38 gene will provide effective control of BB in the Basmati growing region.

Improvement of Basmati rice varieties for resistance to blast and bacterial blight diseases using marker assisted backcross breeding

Marker assisted backcross breeding was employed to incorporate the blast resistance genes, Pi2 and Pi54 and bacterial blight (BB) resistance genes xa13 and Xa21 into the genetic background of Pusa Basmati 1121 (PB1121) and Pusa Basmati 6. Foreground selection for target gene(s) was followed by arduous phenotypic and background selection which fast-tracked the recovery of recurrent parent genome (RPG) to an extent of 95.8% in one of the near-isogenic lines (NILs) namely, Pusa 1728-23-33-31-56, which also showed high degree of resemblance to recurrent parent, PB6 in phenotype. The phenotypic selection prior to background selection provided an additional opportunity for identifying the novel recombinants viz., Pusa 1884-9-12-14 and Pusa 1884-3-9-175, superior to parental lines in terms of early maturity, higher yield and improved quality parameters. There was no significant difference between the RPG recovery estimated based on SSR or SNP markers, however, the panel of SNPs markers was considered as the better choice for background selection as it provided better genome coverage and included SNPs in the genic regions. Multi-location evaluation of NILs depicted their stable and high mean performance in comparison to the respective recurrent parents. The Pi2+Pi54 carrying NILs were effective in combating a pan-India panel of Magnaporthe oryzae isolates with high level of field resistance in northern, eastern and southern parts of India. Alongside, the PB1121-NILs and PB6-NILs carrying BB resistance genes xa13+Xa21 were resistant against Xanthomonas oryzae pv. oryzae races of north-western, southern and eastern parts of the country. Three of NILs developed in this study, have been promoted to final stage of testing during the ​Kharif 2015 in the Indian National Basmati Trial.

Development and evaluation of near-isogenic lines for major blast resistance gene(s) in Basmati rice

A set of NILs carrying major blast resistance genes in a Basmati rice variety has been developed. Also, the efficacy of pyramids over monogenic NILs against rice blast pathogen Magnaporthe oryzae has been demonstrated. Productivity and quality of Basmati rice is severely affected by rice blast disease. Major genes and QTLs conferring resistance to blast have been reported only in non-Basmati rice germplasm. Here, we report incorporation of seven blast resistance genes from the donor lines DHMASQ164-2a (Pi54, Pi1, Pita), IRBLz5-CA (Pi2), IRBLb-B (Pib), IRBL5-M (Pi5) and IRBL9-W (Pi9) into the genetic background of an elite Basmati rice variety Pusa Basmati 1 (PB1). A total of 36 near-isogenic lines (NILs) comprising of 14 monogenic, 16 two-gene pyramids and six three-gene pyramids were developed through marker-assisted backcross breeding (MABB). Foreground, recombinant and background selection was used to identify the plants with target gene(s), minimize the linkage drag and increase the recurrent parent genome (RPG) recovery (93.5-98.6 %), respectively, in the NILs. Comparative analysis performed using 50,051 SNPs and 500 SSR markers revealed that the SNPs provided better insight into the RPG recovery. Most of the monogenic NILs showed comparable performance in yield and quality, concomitantly, Pusa1637-18-7-6-20 (Pi9), was significantly superior in yield and stable across four different environments as compared to recurrent parent (RP) PB1. Further, among the pyramids, Pusa1930-12-6 (Pi2+Pi5) showed significantly higher yield and Pusa1633-7-8-53-6-8 (Pi54+Pi1+Pita) was superior in cooking quality as compared to RP PB1. The NILs carrying gene Pi9 were found to be the most effective against the concoction of virulent races predominant in the hotspot locations for blast disease. Conversely, when analyzed under artificial inoculation, three-gene pyramids expressed enhanced resistance as compared to the two-gene and monogenic NILs.

Molecular evaluation of genetic diversity and association studies in rice (Oryza sativa L.)

In the present study, we tested rice genotypes that included un(der)exploited landraces of Tamil Nadu along with indica and japonica test cultivars to ascertain their genetic diversity structure. Highly polymorphic microsatellite markers were used for generating marker segregation data. A novel measure, allele discrimination index, was used to determine subpopulation differentiation power of each marker. Phenotypic data were collected for yield and component traits. Pattern of molecular differentiation separated indica and japonica genotypes; indica genotypes had two subpopulations within. Landraces were found to have indica genome, but formed a separate subgroup with low linkage disequilibrium. The landraces further separated into distinct group in both hierarchical clustering analysis using neighbour-joining method as well as in the model based population structure analysis. Japonica and the remaining indica cultivars formed two other distinct groups. Linkage disequilibrium observed in the whole population was considerably reduced in subpopulations. Low linkage disequilibrium of landforms suggests their narrow adaptation in local geographical niche. Many population specific alleles could be identified particularly for japonica cultivars and landraces. Association analysis revealed nine marker-trait associations with three agronomic traits, of which 67% were previously reported. Although the testing landraces together with known cultivars had permitted genomewide association mapping, the experiment offers scope to study more landraces collected from the entire geographical region for drawing more reliable information.