Identification of candidate markers associated with agronomic traits in rice using discriminant analysis

Identification of quantitative trait loci for tillering, root, and shoot biomass at the maximum tillering stage in rice

Tillering and plant biomass are key determinants of rice crop productivity. Tillering at the vegetative stage is associated with weed competition, nutrient uptake, and methane emissions. However, little information is available on quantitative trait loci (QTLs) associated with tiller number (qTN), root biomass (qRB), and shoot biomass (qSB) at the active tillering stage which occurs approximately 6 weeks after planting. Here, we mapped tiller and biomass QTLs with ~ 250 recombinant inbred lines derived from a ‘Francis’ by ‘Rondo’ cross using data collected at the maximum tillering stage from two years of greenhouse study, and further compared these QTLs with those mapped at the harvest stage from a field study. Across these three studies, we discovered six qTNs, two qRBs, and three qSBs. Multiple linear regression further indicated that qTN1-2, qTN3-3, qTN4-1, qRB3-1, and qRB5-1 were significant at the maximum tillering stage while qTN3-2 was detected only at the harvest stage. Moreover, qTN3-1 was consistently significant across different developmental stages and growing environments. The genes identified from the peak target qTN regions included a carotenoid metabolism enzyme, a MYB transcription factor, a CBS domain-containing protein, a SAC3/GANP family protein, a TIFY motif containing protein, and an ABC transporter protein. Two genes in the qRB peak target regions included an expressed protein and a WRKY gene. This knowledge of the QTLs, associated markers, candidate genes, and germplasm resources with high TN, RB and SB is of value to rice cultivar improvement programs.

Identification of QTLs for high grain yield and component traits in new plant types of rice

A panel of 60 genotypes comprising New Plant Types (NPTs) along with indica, tropical and temperate japonica genotypes was phenotypically evaluated for four seasons in irrigated situation for grain yield per se and component traits. Twenty NPT genotypes were found promising with an average grain yield varying from 5.45 to 8.8 t/ha. A total of 85 SSR markers were used in the study to identify QTLs associated with grain yield per se and related traits. Sixty-six (77.65%) markers were found to be polymorphic. The PIC values varied from 0.516 to 0.92 with an average of 0.704. A moderate level of genetic diversity (0.39) was detected among genotypes. Variation to the tune of 8% within genotypes, 68% among the genotypes within the population and 24% among the populations were observed (AMOVA). This information may help in identification of potential parents for development of transgressive segregants with very high yield. The association analysis using GLM and MLM models led to the identification of 30 and 10 SSR markers associated with 70 and 16 QTLs, respectively. Thirty novel QTLs linked with 16 SSRs were identified to be associated with eleven traits, namely tiller number (qTL-6.1, qTL-11.1, qTL-4.1), panicle length (qPL-1.1, qPL-5.1, qPL-7.1, qPL-8.1), flag leaf length (qFLL-8.1, qFLL-9.1), flag leaf width (qFLW-6.2, qFLW-5.1, qFLW-8.1, qFLW-7.1), total no. of grains (qTG-2.2, qTG-a7.1), thousand-grain weight (qTGW-a1.1, qTGW-a9.2, qTGW-5.1, qTGW-8.1), fertile grains (qFG-7.1), seed length-breadth ratio (qSlb-3.1), plant height (qPHT-6.1, qPHT-9.1), days to 50% flowering (qFD-1.1) and grain yield per se (qYLD-5.1, qYLD-6.1a, qYLD-11.1).Some of the SSRs were co-localized with more than two traits. The highest co-localization was identified with RM5709 linked to nine traits, followed by RM297 with five traits. Similarly, RM5575, RM204, RM168, RM112, RM26499 and RM22899 were also recorded to be co-localized with more than one trait and could be rated as important for marker-assisted backcross breeding programs, for pyramiding of these QTLs for important yield traits, to produce new-generation rice for prospective increment in yield potentiality and breaking yield ceiling.

Identification and validation of FaP1D7, a putative marker associated with the biosynthesis of methyl butanoate in cultivated strawberry (Fragaria x ananassa)

Breeding strawberry (Fragaria x ananassa) with enhanced fruit flavour is one of the top breeding goals of many strawberry-producing countries. Although several genes involved in the biosynthetic pathways of key aroma compounds have been identified, the development and application of molecular markers associated with fruit flavour remain limited. This study aims to identify molecular markers closely linked to genes controlling strawberry aroma. A purpose-built Subtracted Diversity Array (SDA) known as Fragaria Discovery Panel (FDP) was used for marker screening. Polymorphic sequences associated with key aroma compounds were identified from two DNA bulks with extreme phenotypes, established using 50 F₁ progeny plants derived from Juliette X 07-102-41 cross, two strawberry genotypes differing in aroma profile. A total of 49 polymorphic markers for eight key aroma compounds were detected using genotypic data of the extreme DNA bulks and phenotypic data obtained from gas chromatography-mass spectrometry (GC-MS). A similarity search against the physical maps of Fragaria vesca revealed that FaP1D7 is linked to genes potentially involved in the synthesis of methyl butanoate. A C/T SNP was detected within the feature, which could possibly be converted to a molecular tool for rapid screening of the strawberry accessions for their methyl butanoate production capacity.

Dissecting Genetic Network of Fruit Branch Traits in Upland Cotton by Association Mapping Using SSR Markers

Genetic architecture of branch traits has large influences on the morphological structure, photosynthetic capacity, planting density, and yield of Upland cotton (Gossypium hirsutum L.). This research aims to reveal the genetic effects of six branch traits, including bottom fruit branch node number (BFBNN), bottom fruit branch length (BFBL), middle fruit branch node number (MFBNN), middle fruit branch length (MFBL), upper fruit branch node number (UFBNN), and upper fruit branch length (UFBL). Association mapping was conducted for these traits of 39 lines and their 178 F₁ hybrids in three environments. There were 20 highly significant Quantitative Trait SSRs (QTSs) detected by mixed linear model approach analyzing a full genetic model with genetic effects of additive, dominance, epistasis and their environment interaction. The phenotypic variation explained by genetic effects ranged from 32.64 ~ 91.61%, suggesting these branch traits largely influenced by genetic factors.

Application of subtracted gDNA microarray-assisted Bulked Segregant Analysis for rapid discovery of molecular markers associated with day-neutrality in strawberry (Fragaria x ananassa)

A Fragaria Discovery Panel (FDP; strawberry-specific SDA) containing 287 features was constructed by subtracting the pooled gDNA of nine non-angiosperm species from the pooled gDNA of five strawberry genotypes. This FDP was used for Bulk Segregant Analysis (BSA) to enable identification of molecular markers associated with day-neutrality. Analysis of hybridisation patterns of a short day (SD) DNA bulk and three day-neutral (DN) DNA bulks varying in flowering strength allowed identification of a novel feature, FaP2E11, closely linked to CYTOKININ OXIDASE 1 (CKX1) gene possibly involved in promoting flowering under non-inductive condition. The signal intensities of FaP2E11 feature obtained from the strong DN bulk (DN1) is three fold higher than the short day bulk (SD), indicating that the putative marker may linked to a CKX1 variant allele with lower enzyme activity. We propose a model for flowering regulation based on the hypothesis that flowering strength may be regulated by the copy number of FaP2E11-linked CKX1 alleles. This study demonstrates the feasibility of the SDA-based BSA approach for the identification of molecular markers associated with day-neutrality in strawberry. This innovative strategy is an efficient and cost-effective approach for molecular marker discovery.

Association Analysis in Rice: From Application to Utilization

Association analysis based on linkage disequilibrium (LD) is an efficient way to dissect complex traits and to identify gene functions in rice. Although association analysis is an effective way to construct fine maps for quantitative traits, there are a few issues which need to be addressed. In this review, we will first summarize type, structure, and LD level of populations used for association analysis of rice, and then discuss the genotyping methods and statistical approaches used for association analysis in rice. Moreover, we will review current shortcomings and benefits of association analysis as well as specific types of future research to overcome these shortcomings. Furthermore, we will analyze the reasons for the underutilization of the results within association analysis in rice breeding.

Association mapping for important agronomic traits in core collection of rice (Oryza sativa L.) with SSR markers

Mining elite genes within rice landraces is of importance for the improvement of cultivated rice. An association mapping for 12 agronomic traits was carried out using a core collection of rice consisting of 150 landraces (Panel 1) with 274 simple sequence repeat (SSR) markers, and the mapping results were further verified using a Chinese national rice micro-core collection (Panel 2) and a collection from a global molecular breeding program (Panel 3). Our results showed that (1) 76 significant (P<0.05) trait-marker associations were detected using mixed linear model (MLM) within Panel 1 in two years, among which 32% were identical with previously mapped QTLs, and 11 significant associations had >10% explained ratio of genetic variation; (2) A total of seven aforementioned trait-marker associations were verified within Panel 2 and 3 when using a general linear model (GLM) and 55 SSR markers of the 76 significant trait-marker associations. However, no significant trait-marker association was found to be identical within three panels when using the MLM model; (3) several desirable alleles of the loci which showed significant trait-marker associations were identified. The research provided important information for further mining these elite genes within rice landraces and using them for rice breeding.

Genetic diversity and association mapping of seed vigor in rice (Oryza sativa L.)

Seed vigor is closely related to direct seeding in rice (Oryza sativa L.). Previous quantitative trait locus (QTL) studies for seed vigor were mainly derived from bi-parental segregating populations and no report from natural populations. In this study, association mapping for seed vigor was performed on a selected sample of 540 rice cultivars (419 from China and 121 from Vietnam). Population structure was estimated on the basis of 262 simple sequence repeat (SSR) markers. Seed vigor was evaluated by root length (RL), shoot length (SL) and shoot dry weight in 2011 and 2012. Abundant phenotypic and genetic diversities were found in the studied population. The population was divided into seven subpopulations, and the levels of linkage disequilibrium (LD) ranged from 10 to 80 cM. We identified 27 marker-trait associations involving 18 SSR markers for three traits. According to phenotypic effects for alleles of the detected QTLs, elite alleles were mined. These elite alleles could be used to design parental combinations and the expected results would be obtained by pyramiding or substituting the elite alleles per QTL (apart from possible epistatic effects). Our results demonstrate that association mapping can complement and enhance previous QTL information for marker-assisted selection and breeding by design.

Association of molecular markers with cold tolerance and green period in zoysiagrass (Zoysia Willd.)

Cold tolerance and the green period are key traits in the breeding of zoysiagrass (Zoysia Willd.). Identification of molecular markers associated with cold tolerance and the green period of zoysiagrass will contribute to efficient selection of elite cultivars. These two traits were measured in 96 zoysiagrass accessions in 2004 and 2005-2006, respectively. The mapping population was screened with 29 pairs of simple sequence repeat (SSR) primers and 54 pairs of sequence-related amplified polymorphism (SRAP) primers. A multi-loci in silico mapping approach implemented with an empirical Bayes method was applied for association mapping of cold tolerance and green period. We detected 254 SSR polymorphic loci and 338 SRAP polymorphic loci, among which three SSR loci (Xgwm131-3B-187, Xgwm469-6D-194 and Xgwm234-5B-244) and one SRAP locus (Me11Em7-406) were significantly associated with cold tolerance with effect values of 57.83%, 38.05%, 36.92% and 37%, respectively. Three SSR loci (Xgwm132-6B-225, Xgwm111-7D-34 and Xgwm102-2D-97) and two SRAP loci (Me19Em5-359 and Me16Em8-483) were significantly associated with the green period with effect values of 79.54%, 62.59%, 99.04%, 49.01% and 82.57%. These markers will be useful for genetic improvement of the cold tolerance and green period of zoysiagrass by marker-assisted breeding.

Identifying loci influencing grain number by microsatellite screening in bread wheat (Triticum aestivum L.)

Grain number (GN) is one of three major yield-related components in wheat. We used the Chinese wheat mini core collection to undertake a genome-wide association analysis of grain number using 531 SSR markers randomly located on all 21 chromosomes. Grain numbers of all accessions were measured in four trials, i.e. two environments in four growing seasons. Association analysis based on a mixed linear model (MLM) revealed that 27 SSR loci were significantly associated with mean GN (MGN) estimated by the best linear unbiased predictor (BLUP) method. These included numerous breeder favorable alleles with strong positive effects at 23 loci. Significant or extremely significant differences were detected on MGN between varieties conveying favored allele and varieties with other alleles. Moreover, statistical simulation showed that the favored alleles have additive genetic effects. Although modern varieties combined larger numbers of favored alleles, the numbers of favored alleles were not significantly different from those in landraces, especially those alleles contributing mostly to the phenotypic variation. These results indicate that there is still considerable genetic potential for use of markers for genome selection of GN for high yield in wheat.

Profiling and association mapping of grain metabolites in a subset of the core collection of Chinese rice germplasm (Oryza sativa L.)

In this study, metabolic profiles of a set of 48 rice germplasms from the Chinese core collection were obtained by gas chromatography and time-of-flight mass spectrometry (GC-TOF-MS). Forty-one metabolites were identified and relatively quantified according to the internal standard (IS). Wide ranges of variations for all metabolites were observed among rice accessions. The maximum/minimum ratios varied from 4.73 to 211.36. The metabolites were categorized into seven groups based on their chemical characteristics. Clustering analysis and a correlation network showed that most of the metabolites had variations among rice accessions in the same direction. Using 218 molecular markers, association mapping was conducted to identify the chromosomal loci influencing the concentrations of identified metabolites. Twenty markers were identified associating with the concentrations of 29 metabolites [-lg(P) > 3]. Allelic effects were investigated in detail in two markers (RM315 and RM541) as examples.

Application of association mapping to understanding the genetic diversity of plant germplasm resources

Compared to the conventional linkage mapping, linkage disequilibrium (LD)-mapping, using the nonrandom associations of loci in haplotypes, is a powerful high-resolution mapping tool for complex quantitative traits. The recent advances in the development of unbiased association mapping approaches for plant population with their successful applications in dissecting a number of simple to complex traits in many crop species demonstrate a flourish of the approach as a “powerful gene tagging” tool for crops in the plant genomics era of 21st century. The goal of this review is to provide nonexpert readers of crop breeding community with (1) the basic concept, merits, and simple description of existing methodologies for an association mapping with the recent improvements for plant populations, and (2) the details of some of pioneer and recent studies on association mapping in various crop species to demonstrate the feasibility, success, problems, and future perspectives of the efforts in plants. This should be helpful for interested readers of international plant research community as a guideline for the basic understanding, choosing the appropriate methods, and its application.

Maintaining the yield of edible rice in a warming world

High temperature increases the amount of chalk in rice (Oryza sativa L.) grains, which causes grains to break during polishing, lowering the amount of rice for consumption. Here, we examined the effect of elevated temperature on substrate supply to the panicle, the capacity of the panicle to produce edible grains, and underlying factors affecting yield of edible grain in two varieties. During grain-filling, substrate supply followed a bell shaped curve, and high temperature significantly shortened supply time. The rate of grain-filling did not change and paddy yield fell in both varieties. In high temperature, yield loss in IR8 was due to lighter grains relative to those grown in cool temperature, but in IR60, it was due to the early sacrifice of 30% of the spikelets. The yield of edible rice was zero for IR8 and ~60% for IR60 for the high temperature treatments, and 100% for IR60 and 70% for IR8 in the cool temperature. IR60 differs from IR8 in regulation of substrate supply, architecture of the panicles and the capacity of the panicles to alter sink size in response to the stress and these factors may be responsible for the difference in edible rice in the two varieties.

Population structure and association mapping on chromosome 7 using a diverse panel of Chinese germplasm of rice (Oryza sativa L.)

The majority of 170 rice accessions used in this study were diverse landraces or varieties from a putative mini-core collection of Chinese germplasm along with some widely used parental lines in genetic analysis or breeding (a few from abroad). The population was genotyped using 84 SSR or InDel markers on chromosome 7 and 48 markers on other chromosomes. The phenotyping of heading date, plant height and panicle length were carried out in different locations for 2 years. Based on morphological characterization, distance-based clustering and model-based estimation of marker data, the population showed a predominant structure with two subpopulations in correspondence with indica and japonica subspecies. The estimation of linkage disequilibrium in 2 Mb windows varied along chromosome 7 and showed parallel changes with inter-subspecies differentiation of marker loci (Fst). Based on the mixed linear model considering population structure and family relatedness [i.e. the (Q + K) model], one to three associated markers (P < or = 0.0001) per trait per experiment were scanned out on rice chromosome 7. Most significant loci were repeated for the data from both field experiments while two loci were associated with two or three traits. Marker-based allelic effects were shown in a couple of associated markers as examples. The application of association results in breeding program was also discussed.

Genetic variation in rice

Completion of the genomic sequencing of rice has enhanced the discovery of new genes. Wild rice relatives are good sources for extending the genetic variation of cultivated rice. Reproductive barriers are commonly found in distant crosses of rice and are attracting attention. The combination of genetic analyses and molecular tools has greatly facilitated the molecular cloning of rice genes based on the classical approach and enabled the tracking of dissemination of the alleles for domestication. Basic information for population genetics study in rice is still being collected and is expected to provide an alternative approach for finding new genes. The wide genetic variation available in wild rice relatives and the combination of various genetic approaches will allow the analysis and understanding of genetic variation at the nucleotide sequence level, as well as the discovery of novel alleles by sequence-based approaches.

Genetic diversity, structure and marker-trait associations in a collection of Italian tomato (Solanum lycopersicum L.) landraces

The study of phenotypic and genetic diversity in landrace collections is important for germplasm conservation. In addition, the characterisation of very diversified materials with molecular markers offers a unique opportunity to define significant marker-trait associations of biological and agronomic interest. Here, 50 tomato landraces (mainly collected in central Italy), nine vintage and modern cultivars, and two wild outgroups were grown at two locations in central Italy and characterised for 15 morpho-physiological traits and 29 simple sequence repeat (SSR) loci. The markers were selected to include a group of loci in regions harbouring reported quantitative trait loci (QTLs) that affect fruit size and/or shape (Q-SSRs) and a group of markers that have not been mapped or shown to have a priori known linkage (NQ-SSRs). As revealed by univariate and multivariate analyses of morphological data, the landraces grouped according to vegetative and reproductive traits, with emphasis on fruit size, shape and final destination of the product. Compared to the low molecular polymorphism reported in tomato modern cultivars, our data reveal a high level of molecular diversity in landraces. Such diversity has allowed the inference of the existence of a genetic structure that was factored into the association analysis. As the proportion of significant associations is higher between the Q-SSR subset of markers and the subset of traits related to fruit size and shape than for all of the other combinations, we conclude that this approach is valid for establishing true-positive marker-trait relationships in tomato.

Rice molecular breeding laboratories in the genomics era: Current status and future considerations

Using DNA markers in plant breeding with marker-assisted selection (MAS) could greatly improve the precision and efficiency of selection, leading to the accelerated development of new crop varieties. The numerous examples of MAS in rice have prompted many breeding institutes to establish molecular breeding labs. The last decade has produced an enormous amount of genomics research in rice, including the identification of thousands of QTLs for agronomically important traits, the generation of large amounts of gene expression data, and cloning and characterization of new genes, including the detection of single nucleotide polymorphisms. The pinnacle of genomics research has been the completion and annotation of genome sequences for indica and japonica rice. This information-coupled with the development of new genotyping methodologies and platforms, and the development of bioinformatics databases and software tools-provides even more exciting opportunities for rice molecular breeding in the 21st century. However, the great challenge for molecular breeders is to apply genomics data in actual breeding programs. Here, we review the current status of MAS in rice, current genomics projects and promising new genotyping methodologies, and evaluate the probable impact of genomics research. We also identify critical research areas to "bridge the application gap" between QTL identification and applied breeding that need to be addressed to realize the full potential of MAS, and propose ideas and guidelines for establishing rice molecular breeding labs in the postgenome sequence era to integrate molecular breeding within the context of overall rice breeding and research programs.

Bayesian association mapping of multiple quantitative trait loci and its application to the analysis of genetic variation among Oryza sativa L. germplasms

One way to use a crop germplasm collection directly to map QTLs without using line-crossing experiments is the whole genome association mapping. A major problem with association mapping is the presence of population structure, which can lead to both false positives and failure to detect genuine associations (i.e., false negatives). Particularly in highly selfing species such as Asian cultivated rice, high levels of population structure are expected and therefore the efficiency of association mapping remains almost unknown. Here, we propose an approach that combines a Bayesian method for mapping multiple QTLs with a regression method that directly incorporates estimates of population structure. That is, the effects due to both multiple QTLs and population structure were included in our statistical model. We evaluated the efficiency of our approach in simulated- and real-trait analyses of a rice germplasm collection. Simulation analyses based on real marker data showed that our model could suppress both false-positive and false-negative rates and the error of estimation of genetic effects over single QTL models, indicating that our model has statistically desirable attributes over single QTL models. As real traits, we analyzed the size and shape of milled rice grains and found significant markers that may be linked to QTLs reported previously. Association mapping should have good prospects in highly selfing species such as rice if proper methods are adopted. Our approach will be useful for the whole genome association mapping of various selfing crop species.

Genetic and molecular basis of grain size and grain number and its relevance to grain productivity in higher plants

Grain size and grain number constitute 2 important components of grain yield. In particular, the grain size also influences the end-use quality (e.g., flour yield and protein content) and attracts consumer preference. These 2 traits are also the components of the domestication syndrome of crop plants. A number of important studies have recently been conducted to understand the genetic and molecular basis of these 2 important yield-contributing traits. Information generated from these studies was collected and synthesized for the benefit of plant biologists, particularly plant breeders. In the present article, this information is briefly reviewed and the prospects of using this information for improvement of grain productivity in crop plants are discussed.

Haplotype diversity: the link between statistical and biological association

In the rapidly growing field of association mapping in plants, the use of (marker) haplotypes rather than single markers can be an effective way of improving detection power. Here, we highlight the information that can be obtained from deducing the historical relationships between haplotypes. The ordering of haplotype classes according to deduced historical relationships should further enhance association detection power, but can also be used to predict the genotypic and phenotypic values of unobserved germplasm.