In Silico Identification and Experimental Validation of Insertion–Deletion Polymorphisms in Tomato Genome

Comparative analysis of the genome sequences of Solanum lycopersicum variety Heinz 1706 and S. pimpinellifolium accession LA 1589 using MUGSY software identified 145 695 insertion–deletion (InDel) polymorphisms. A selected set of 3029 candidate InDels (≥2 bp) across the entire tomato genome were subjected to PCR validation, and 82.4% could be verified. Of 2272 polymorphic InDels between LA 1589 and Heinz 1706, 61.6, 45.2, and 31.6% were polymorphic in 8 accessions of S. pimpinellifolium, 4 accessions of S. lycopersicum var. cerasiforme, and 10 varieties of S. lycopersicum, respectively. Genetic distance was 0.216 in S. pimpinellifolium, 0.202 in S. lycopersicum var. cerasiforme, and 0.108 in S. lycopersicum. The data suggested a reduction of genetic variation from S. pimpinellifolium to S. lycopersicum var. cerasiforme and S. lycopersicum. Cluster analysis showed that the 8 accessions of S. pimpinellifolium were in one group, whereas 4 accessions of S. lycopersicum var. cerasiforme and 10 varieties of S. lycopersicum were in the same group.     

Detection of trait donors and QTL boundaries for early blight resistance using local ancestry inference in a library of genomic sequences for tomato

By conducting hierarchical clustering along a sliding window, we generated haplotypes across hundreds of re-sequenced genomes in a few hours. We leveraged our method to define cryptic introgressions underlying disease resistance in tomato (Solanum lycopersicum L.) and to discover resistant germplasm in the tomato seed bank. The genomes of 9 accessions with early blight (Alternaria linariae) disease resistance were newly sequenced and analyzed together with published sequences for 770 tomato and wild species accessions, most of which are available in germplasm collections. Identification of common ancestral haplotypes among resistant germplasm enabled rapid fine mapping of recently discovered quantitative trait loci (QTL) conferring resistance and the identification of possible causal variants. The source of the early blight QTL EB-9 was traced to a vintage tomato named ‘Devon Surprise’. Another QTL, EB-5, as well as resistance to bacterial spot disease (Xanthomonas spp.), was traced to Hawaii 7998. A genomic survey of all accessions forecasted EB-9-derived resistance in several heirloom tomatoes, accessions of S. lycopersicum var. cerasiforme, and S. pimpinellifolium PI 37009. Our haplotype-based predictions were validated by screening the accessions against the causal pathogen. There was little evidence of EB-5 prevalence in surveyed contemporary germplasm, presenting an opportunity to bolster tomato disease resistance by adding this rare locus. Our work demonstrates practical insights that can be derived from the efficient processing of large genome-scale datasets, including rapid functional prediction of disease resistance QTL in diverse genetic backgrounds. Finally, our work finds more efficient ways to leverage public genetic resources for crop improvement.     

Association Mapping for Fruit,Plant and Leaf Morphology Traits in Eggplant

An eggplant (Solanum melongena) association panel of 191 accessions, comprising a mixture of breeding lines, old varieties and landrace selections was SNP genotyped and phenotyped for key breeding fruit and plant traits at two locations over two seasons. A genome-wide association (GWA) analysis was performed using the mixed linear model, which takes into account both a kinship matrix and the sub-population membership of the accessions. Overall, 194 phenotype/genotype associations were uncovered, relating to 30 of the 33 measured traits. These associations involved 79 SNP loci mapping to 39 distinct chromosomal regions distributed over all 12 eggplant chromosomes. A comparison of the map positions of these SNPs with those of loci derived from conventional linkage mapping showed that GWA analysis both validated many of the known controlling loci and detected a large number of new marker/trait associations. Exploiting established syntenic relationships between eggplant chromosomes and those of tomato and pepper recognized orthologous regions in ten eggplant chromosomes harbouring genes influencing breeders’ traits.     

Genomic variation in tomato,from wild ancestors to contemporary breeding accessions

Background

Domestication modifies the genomic variation of species. Quantifying this variation provides insights into the domestication process, facilitates the management of resources used by breeders and germplasm centers, and enables the design of experiments to associate traits with genes. We described and analyzed the genetic diversity of 1,008 tomato accessions including Solanum lycopersicum var. lycopersicum (SLL), S. lycopersicum var. cerasiforme (SLC), and S. pimpinellifolium (SP) that were genotyped using 7,720 SNPs. Additionally, we explored the allelic frequency of six loci affecting fruit weight and shape to infer patterns of selection.

Results

Our results revealed a pattern of variation that strongly supported a two-step domestication process, occasional hybridization in the wild, and differentiation through human selection. These interpretations were consistent with the observed allele frequencies for the six loci affecting fruit weight and shape. Fruit weight was strongly selected in SLC in the Andean region of Ecuador and Northern Peru prior to the domestication of tomato in Mesoamerica. Alleles affecting fruit shape were differentially selected among SLL genetic subgroups. Our results also clarified the biological status of SLC. True SLC was phylogenetically positioned between SP and SLL and its fruit morphology was diverse. SLC and “cherry tomato” are not synonymous terms. The morphologically-based term “cherry tomato” included some SLC, contemporary varieties, as well as many admixtures between SP and SLL. Contemporary SLL showed a moderate increase in nucleotide diversity, when compared with vintage groups.

Conclusions

This study presents a broad and detailed representation of the genomic variation in tomato. Tomato domestication seems to have followed a two step-process; a first domestication in South America and a second step in Mesoamerica. The distribution of fruit weight and shape alleles supports that domestication of SLC occurred in the Andean region. Our results also clarify the biological status of SLC as true phylogenetic group within tomato. We detect Ecuadorian and Peruvian accessions that may represent a pool of unexplored variation that could be of interest for crop improvement.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1444-1) contains supplementary material, which is available to authorized users.     

Genome-wide association study reveals the genetic architecture of 27 agronomic traits in tomato

Tomato (Solanum lycopersicum) is a highly valuable fruit crop, and yield is one of the most important agronomic traits. However, the genetic architecture underlying tomato yield-related traits has not been fully addressed. Based on ∼4.4 million single nucleotide polymorphisms obtained from 605 diverse accessions, we performed a comprehensive genome-wide association study for 27 agronomic traits in tomato. A total of 239 significant associations corresponding to 129 loci, harboring many previously reported and additional genes related to vegetative and reproductive development, were identified, and these loci explained an average of ∼8.8% of the phenotypic variance. A total of 51 loci associated with 25 traits have been under selection during tomato domestication and improvement. Furthermore, a candidate gene, Sl-ACTIVATED MALATE TRANSPORTER15, that encodes an aluminum-activated malate transporter was functionally characterized and shown to act as a pivotal regulator of leaf stomata formation, thereby affecting photosynthesis and drought resistance. This study provides valuable information for tomato genetic research and breeding.

A large-scale genome-wide association study sheds light on genetic and genomic bases underlying important yield-related traits in tomato.     

Genetic diversity, linkage disequilibrium, and association mapping analyses of peach (Prunus persica) landraces in China

The genetic diversity, population structure, and linkage disequilibrium (LD) of peaches are greatly important in genome-wide association mapping. In the current study, 104 peach landrace accessions from six Chinese geographical regions were evaluated for fruit and phenological period. The accessions were genotyped with 53 genome-wide simple sequence repeat (SSR) markers. All SSR markers were highly polymorphic across the accessions, and a total of 340 alleles were detected, including 59 private alleles. Of the six regions studied, the northern part of China as well as the middle and lower reaches of the Changjiang River were found to be the most highly diverse genetically. Based on population structure analysis, the peaches were divided into five groups, which well agreed with the geographical distribution. Of the SSR pairs in these accessions, 18.07% (P?<?0.05) were in LD. The mean r ² value for all intrachromosomal loci pairs was 0.0149, and LD decayed at 6.01?cM. The general linear model was used to calculate the genome-wide marker-trait associations of 10 complex traits. The traits include flesh color around the stone, red pigment in the flesh, flesh texture, flesh adhesion, flesh firmness, fruit weight, chilling requirement, flowering time, ripening time, and fruit development period. These traits were estimated by analyzing the 104 landraces. Many of the associated markers were located in regions where quantitative trait loci (QTLs) were previously identified. Peach association mapping is an effective approach for identifying QTLs and may be an alternative to QTL mapping based on crosses between different lines.     

Cuticular wax variation in the tomato (Solanum lycopersicum L.), related wild species and their interspecific hybrids

The tomato (Solanum lycopersicum L.) is one of the world's most important vegetable crop species. Among the many tomato accessions available, only a few are tolerant to abiotic stresses, which are responsible for the majority of the crop losses worldwide. Wild tomato species are then secondary gene pool in the breeding of more resistant tomato cultivars. In the current study, the composition of leaf cuticular waxes from fourteen tomato accessions, including S. lycopersicum, Solanum pennellii, Solanum pimpinellifolium, and their interspecific hybrids was studied in order to select the most adequate chemotaxonomic markers. Total cuticular wax load of S. pennellii plants was much higher than in the other plant species. Hydrocarbons were usually the most abundant wax components, followed by minor quantities of triterpenes and other compounds. Interspecific hybrids showed intermediate wax characteristics. The amount and composition of surface waxes were not correlated with the abiotic stress tolerance in S. lycopersicum. The composition of the hydrocarbon fraction was the least variable both within a single accession and between all the plants studied. Based on the results, cuticular hydrocarbons are proposed as potential chemotaxonomic markers in the classification of tomato and related species.     

Vegetative performance,leaf water potential,and partitioning of minerals and soluble sugars: Traits for ranking the NaCl-tolerance of tomato genotypes?

This study addresses the physiological response of four tomato genotypes with distinct sensitivity to high NaCl concentrations, with the aim of identifying physiological traits to rank the genotypes’ sensitivity to salt stress. The central hypothesis was that tomato genotypes grown in saline environments show a characteristic salinity-triggered absorption and translocation of ions, leading to a distinct distribution pattern of Na, K, and soluble sugars. Experiments were conducted on two commercial cultivars: one assumed to be sensitive (Solanum lycopersicum L. F1 hybrid Harzfeuer) and one known to be tolerant (S. lycopersicum L. var. edkawi) to high salt concentrations. Furthermore, two wild salt-tolerant relatives (S. pennellii and S. lycopersicum var. cerasiforme) were selected. Based on our results regarding vegetative performance and partitioning of Na, K, glucose, fructose, and sucrose, it is possible to classify the genotype S. lycopersicum F1 hybrid Harzfeuer as moderately sensitive to salt stress and the genotypes S. lycopersicum var. edkawi, S. lycopersicum var. cerasiforme, and S. pennellii as moderately resistant to 210 mM NaCl. Calculations of the percentage of modification revealed non-specific genotype responses for the amount of sodium in roots and leaves, as well as the sucrose concentration and the osmotic potential of leaves. As shown, the salt-induced changes in potassium levels in leaves, and glucose concentration in roots might be used as additional traits to discriminate genotypes regarding their salt-sensitivity. These parameters might be useful when comparisons of fruit production or vegetative performance provide no conclusive indication e.g. due to the distinct growth habits of commercial cultivars selected for fruit productivity and native genotypes having a stronger vegetative development. However, additional studies should be conducted to evaluate a large number of genotypes differing in their NaCl tolerance. Furthermore, plant responses considering the dynamic source-sink relations due to fruit load needs to be considered.     

Genetic diversity,population structure,and linkage disequilibrium of elite and local apple accessions from Belgium using the IRSC array

The identification of molecular markers associated with economic and quality traits will help improve breeding for new apple (Malus × domestica Borkh.) cultivars. Tools such as the 8K apple SNP array developed by the RosBREED consortium allow for high-throughput genotyping of SNP polymorphisms within collections. However, genetic characterization and the identification of population stratification and kinship within germplasm collections is a fundamental prerequisite for identifying robust marker–trait associations. In this study, a collection of apple germplasm originally developed for plant architectural studies and consisting of both non-commercial/local and elite accessions was genotyped using the 8K apple SNP array to identify cryptic relationships between accessions, to analyze population structure and to calculate the linkage disequilibrium (LD). A total of nine pairs of synonyms and several triploids accessions were identified within the 130 accessions genotyped. In addition, most of the known parent-child relations were confirmed, and several putative, previously unknown parent-child relations were identified among the local accessions. No clear subgroups could be identified although some separation between local and elite accessions was evident. The study of LD showed a rapid decay in our collection, indicating that a larger number of SNPs is necessary to perform whole genome association mapping. Finally, an association mapping effort for architectural traits was carried out on a small number of accessions to estimate the feasibility of this approach.     

Linkage Disequilibrium and Genome-Wide Association Mapping in Tetraploid Wheat (Triticum turgidum L.)

Association mapping is a powerful tool for the identification of quantitative trait loci through the exploitation of the differential decay of linkage disequilibrium (LD) between marker loci and genes of interest in natural and domesticated populations. Using a sample of 230 tetraploid wheat lines (Triticum turgidum ssp), which included naked and hulled accessions, we analysed the pattern of LD considering 26 simple sequence repeats and 970 mostly mapped diversity array technology loci. In addition, to validate the potential for association mapping in durum wheat, we evaluated the same genotypes for plant height, heading date, protein content, and thousand-kernel weight. Molecular and phenotypic data were used to: (i) investigate the genetic and phenotypic diversity; (ii) study the dynamics of LD across the durum wheat genome, by investigating the patterns of LD decay; and (iii) test the potential of our panel to identify marker–trait associations through the analysis of four quantitative traits of major agronomic importance. Moreover, we compared and validated the association mapping results with outlier detection analysis based on population divergence. Overall, in tetraploid wheat, the pattern of LD is extremely population dependent and is related to the domestication and breeding history of durum wheat. Comparing our data with several other studies in wheat, we confirm the position of many major genes and quantitative trait loci for the traits considered. Finally, the analysis of the selection signature represents a very useful complement to validate marker–trait associations.     

Identification of Genetic Loci Associated with Quality Traits in Almond via Association Mapping

To design an appropriate association study, we need to understand population structure and the structure of linkage disequilibrium within and among populations as well as in different regions of the genome in an organism. In this study, we have used a total of 98 almond accessions, from five continents located and maintained at the Centro de Investigación y Tecnología Agroalimentaria de Aragón (CITA; Spain), and 40 microsatellite markers. Population structure analysis performed in ‘Structure’ grouped the accessions into two principal groups; the Mediterranean (Western-Europe) and the non-Mediterranean, with K = 3, being the best fit for our data. There was a strong subpopulation structure with linkage disequilibrium decaying with increasing genetic distance resulting in lower levels of linkage disequilibrium between more distant markers. A significant impact of population structure on linkage disequilibrium in the almond cultivar groups was observed. The mean r² value for all intra-chromosomal loci pairs was 0.040, whereas, the r² for the inter-chromosomal loci pairs was 0.036. For analysis of association between the markers and phenotypic traits, five models comprising both general linear models and mixed linear models were selected to test the marker trait associations. The mixed linear model (MLM) approach using co-ancestry values from population structure and kinship estimates (K model) as covariates identified a maximum of 16 significant associations for chemical traits and 12 for physical traits. This study reports for the first time the use of association mapping for determining marker-locus trait associations in a world-wide almond germplasm collection. It is likely that association mapping will have the most immediate and largest impact on the tier of crops such as almond with the greatest economic value.     

Genetic analysis of physiological components of salt tolerance conferred by Solanum rootstocks. What is the rootstock doing for the scion?

Grafting desirable crop varieties on stress-tolerant rootstocks provides an opportunity to increase crop salt tolerance. Here, a commercial hybrid tomato variety was grafted on two populations of recombinant inbred lines developed from a salt-sensitive genotype of Solanum lycopersicum var. cerasiforme, as female parent, and two salt-tolerant lines, as male parents, from S. pimpinellifolium, the P population, and S. cheesmaniae, the C population, to identify an easy screening method for identifying rootstocks conferring salt tolerance in terms of fruit yield. Potential physiological components of salt tolerance were assessed in the scion: leaf biomass, [Na⁺], nutrition, water relations and xylem ABA concentration. A significant correlation between scion fruit yield and scion leaf fresh weight, water potential or the ABA concentration was found in the C population under salinity, but the only detected QTL did not support this relationship. The rootstocks of the P population clearly affected seven traits related to the sodium, phosphorous and copper concentrations and water content of the scion leaf, showing heritability estimates around 0.4 or higher. According to heritability estimates in the P population, up to five QTLs were detected per trait. QTLs contributing over 15% to the total variance were found for P and Cu concentrations and water content of the scion leaf, and the proportion of fresh root weight. Correlation and QTL analysis suggests that rootstock-mediated improvement of fruit yield in the P population under salinity is mainly explained by the rootstock’s ability to minimise perturbations in scion water status.     

Association of Agronomic Traits with SNP Markers in Durum Wheat (Triticum turgidum L. durum (Desf.))

Association mapping is a powerful approach to detect associations between traits of interest and genetic markers based on linkage disequilibrium (LD) in molecular plant breeding. In this study, 150 accessions of worldwide originated durum wheat germplasm (Triticum turgidum spp. durum) were genotyped using 1,366 SNP markers. The extent of LD on each chromosome was evaluated. Association of single nucleotide polymorphisms (SNP) markers with ten agronomic traits measured in four consecutive years was analyzed under a mix linear model (MLM). Two hundred and one significant association pairs were detected in the four years. Several markers were associated with one trait, and also some markers were associated with multiple traits. Some of the associated markers were in agreement with previous quantitative trait loci (QTL) analyses. The function and homology analyses of the corresponding ESTs of some SNP markers could explain many of the associations for plant height, length of main spike, number of spikelets on main spike, grain number per plant, and 1000-grain weight, etc. The SNP associations for the observed traits are generally clustered in specific chromosome regions of the wheat genome, mainly in 2A, 5A, 6A, 7A, 1B, and 6B chromosomes. This study demonstrates that association mapping can complement and enhance previous QTL analyses and provide additional information for marker-assisted selection.     

Linkage and Association Mapping for Two Major Traits Used in the Maritime Pine Breeding Program: Height Growth and Stem Straightness

BackgroundIncreasing our understanding of the genetic architecture of complex traits, through analyses of genotype-phenotype associations and of the genes/polymorphisms accounting for trait variation, is crucial, to improve the integration of molecular markers into forest tree breeding. In this study, two full-sib families and one breeding population of maritime pine were used to identify quantitative trait loci (QTLs) for height growth and stem straightness, through linkage analysis (LA) and linkage disequilibrium (LD) mapping approaches.ResultsThe populations used for LA consisted of two unrelated three-generation full-sib families (n = 197 and n = 477). These populations were assessed for height growth or stem straightness and genotyped for 248 and 217 markers, respectively. The population used for LD mapping consisted of 661 founders of the first and second generations of the breeding program. This population was phenotyped for the same traits and genotyped for 2,498 single-nucleotide polymorphism (SNP) markers corresponding to 1,652 gene loci. The gene-based reference genetic map of maritime pine was used to localize and compare the QTLs detected by the two approaches, for both traits. LA identified three QTLs for stem straightness and two QTLs for height growth. The LD study yielded seven significant associations (P ≤ 0.001): four for stem straightness and three for height growth. No colocalisation was found between QTLs identified by LA and SNPs detected by LD mapping for the same trait.ConclusionsThis study provides the first comparison of LA and LD mapping approaches in maritime pine, highlighting the complementary nature of these two approaches for deciphering the genetic architecture of two mandatory traits of the breeding program.     

Association mapping of agronomic traits on chromosome 2A of wheat

Association mapping is a method to test the association between molecular markers and quantitative trait loci (QTL) based on linkage disequilibrium (LD). In this study, the collection of 108 wheat germplasm accessions form China were evaluated for their plant heights, spike length, spikelets per spike, grains per spike, thousand kernel weight and spikelets density in 3 years at three locations. And they were genotyped with 85 SSR markers and 40 EST-SSR markers. The population structure was inferred on the basis of unlinked 48 SSR markers and 40 EST-SSR markers. The extent of LD on chromosome 2A was 2.3 cM. Association of 37 SSR loci on chromosomes 2A with six agronomic traits was analysed with a mixed linear model. A total of 14 SSR loci were significantly associated with agronomic traits. Some of the associated markers were located in the QTL region detected in previous linkage mapping analysis. Our results demonstrated that association mapping can enhance QTL information and achieves higher resolution with short LD extent.     

Cell wall glycosidase activities and protein content variations during fruit development and ripening in three texture contrasted tomato cultivars

Excessive softening is the main factor limiting fruit shelf life and storage. It is generally acceptable now that softening of fruit which occurs during the ripening is due to synergistic actions of several enzymes on cell wall polysaccharides. As a subject for this study, we have assayed some glycosidase activities using three tomato species (Lycopersicon esculentum) contrasted for their texture phenotypes; the cherry tomato line Cervil (Solanum lycopersicum var. cerasiforme), a common taste tomato line Levovil (S. lycopersicum Mill.) and VilB a modern line, large, firmer and with good storage capability. Four glycosidase activities namely α-galactosidase, β-galactosidase, β-mannosidase and β-glucosidase were extracted from tomato’s cell wall of the three species. Cell wall protein from fruits pericarp was extracted and compared among the three cultivars at the following stages; 14 days post anthesis (14DPA) fruit; 21 days post anthesis (21DPA), turning (breaker), red and over ripe. When glycolytic activities were also compared among these cultivars at the precited development stages, gross variations were noticed from stage to stage and also from species to species in accordance with the fruit firmness status. Interestingly, VilB cultivar, the firmer among the other two, though possessed the highest total protein content, exhibited the lowest enzymatic activities. Taken together, these results may therefore allow us to conclude that studies of glycolytic activities in a single tomato cultivar cannot be generalized to all species. On the other hand, relating fruit development to glycosidase activities should logically be coupled to these enzymes from cell wall compartment.     

Genetic diversity and population structure of a diverse set of rice germplasm for association mapping

Germplasm diversity is the mainstay for crop improvement and genetic dissection of complex traits. Understanding genetic diversity, population structure, and the level and distribution of linkage disequilibrium (LD) in target populations is of great importance and a prerequisite for association mapping. In this study, 100 genome-wide simple sequence repeat (SSR) markers were used to assess genetic diversity, population structure, and LD of 416 rice accessions including landraces, cultivars and breeding lines collected mostly in China. A model-based population structure analysis divided the rice materials into seven subpopulations. 63% of the SSR pairs in these accessions were in LD, which was mostly due to an overall population structure, since the number of locus pairs in LD was reduced sharply within each subpopulation, with the SSR pairs in LD ranging from 5.9 to 22.9%. Among those SSR pairs showing significant LD, the intrachromosomal LD had an average of 25–50 cM in different subpopulations. Analysis of the phenotypic diversity of 25 traits showed that the population structure accounted for an average of 22.4% of phenotypic variation. An example association mapping for starch quality traits using both the candidate gene mapping and genome-wide mapping strategies based on the estimated population structure was conducted. Candidate gene mapping confirmed that the Wx and starch synthase IIa (SSIIa) genes could be identified as strongly associated with apparent amylose content (AAC) and pasting temperature (PT), respectively. More importantly, we revealed that the Wx gene was also strongly associated with PT. In addition to the major genes, we found five and seven SSRs were associated with AAC and PT, respectively, some of which have not been detected in previous linkage mapping studies. The results suggested that the population may be useful for the genome-wide marker–trait association mapping. This new association population has the potential to identify quantitative trait loci (QTL) with small effects, which will aid in dissecting complex traits and in exploiting the rich diversity present in rice germplasm.     

Mapping QTLs for improving grain yield using the USDA rice mini-core collection

Yield is the most important and complex trait for genetic improvement in crops, and marker-assisted selection enhances the improvement efficiency. The USDA rice mini-core collection derived from over 18,000 accessions of global origins is an ideal panel for association mapping. We phenotyped 203 O. sativa accessions for 14 agronomic traits and identified 5 that were highly and significantly correlated with grain yield per plant: plant height, plant weight, tillers, panicle length, and kernels/branch. Genotyping with 155 genome-wide molecular markers demonstrated 5 main cluster groups. Linkage disequilibrium (LD) decayed at least 20 cM and marker pairs with significant LD ranged from 4.64 to 6.06% in four main groups. Model comparisons revealed that different dimensions of principal component analysis affected yield and its correlated traits for mapping accuracy, and kinship did not improve the mapping in this collection. Thirty marker–trait associations were highly significant, 4 for yield, 3 for plant height, 6 for plant weight, 9 for tillers, 5 for panicle length and 3 for kernels/branch. Twenty-one markers contributed to the 30 associations, because 8 markers were co-associated with 2 or more traits. Allelic analysis of OSR13, RM471 and RM7003 for their co-associations with yield traits demonstrated that allele 126 bp of RM471 and 108 bp of RM7003 should receive greater attention, because they had the greatest positive effect on yield traits. Tagging the QTLs responsible for multiple yield traits may simultaneously help dissect the complex yield traits and elevate the efficiency to improve grain yield using marker-assisted selection in rice.