Molecular dissection of developmental behavior of plant height in rice (Oryza sativa L.)

Meta-analysis of QTLs and candidate genes associated with seed germination in rice (Oryza sativa L.)

Seed germination is one of the critical stages of plant life, and many quantitative trait loci (QTLs) control this complex trait. Meta-analysis of QTLs is a powerful computational technique for estimating the most stable QTLs regardless of the population's genetic background. Besides, this analysis effectively narrows down the confidence interval (CI) to identify candidate genes (CGs) and marker development. In the current study, a comprehensive genome-wide meta-analysis was performed on QTLs associated with germination in rice. This analysis was conducted based on the data reported over the last two decades. In this case, various analyses were performed, including seed germination rate, plumule length, radicle length, germination percentage, coleoptile length, coleorhiza length, radicle fresh weight, germination potential, and germination index. A total of 67 QTLs were projected onto a reference map for these traits and then integrated into 32 meta-QTLs (MQTLs) to provide a genetic framework for seed germination. The average CI of MQTLs was considerably reduced from 15.125 to 8.73 cM compared to the initial QTLs. This situation identified 728 well-known functionally characterized genes and novel putative CGs for investigated traits. The fold change calculation demonstrated that 155 CGs had significant changes in expression analysis. In this case, 112 and 43 CGs were up-regulated and down-regulated during germination, respectively. This study provides an overview and compares genetic loci controlling traits related to seed germination in rice. The findings can bridge the gap between QTLs and CGs for seed germination.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01232-1.

Dynamic analysis of QTLs on plant height with single segment substitution lines in rice

Dynamic regulation of QTLs remains mysterious. Single segment substitution lines (SSSLs) and conditional QTL mapping and functional QTL mappings are ideal materials and methods to explore dynamics of QTLs for complex traits. This paper analyzed the dynamics of QTLs on plant height with SSSLs in rice. Five SSSLs were verified with plant height QTLs first. All five QTLs had significant positive effects at one or more developmental stages except QTL₁. They interacted each other, with negative effects before 49 d after transplanting and positive effects since then. The five QTLs selectively expressed in specific periods, mainly in the periods from 35 to 42 d and from 49 to 56 d after transplanting. Expressions of epistasis were dispersedly in various periods, negative effects appearing mainly before 35 d. The five QTLs brought the inflexion point ahead of schedule, accelerated growth and degradation, and changed the peak plant height, while their interactions had the opposite effects. The information will be helpful to understand the genetic mechanism for developmental traits.

Statistical Methods for the Quantitative Genetic Analysis of High-Throughput Phenotyping Data

The advent of plant phenomics, coupled with the wealth of genotypic data generated by next-generation sequencing technologies, provides exciting new resources for investigations into and improvement of complex traits. However, these new technologies also bring new challenges in quantitative genetics, namely, a need for the development of robust frameworks that can accommodate these high-dimensional data. In this chapter, we describe methods for the statistical analysis of high-throughput phenotyping (HTP) data with the goal of enhancing the prediction accuracy of genomic selection (GS). Following the Introduction in Sec. 1, Sec. 2 discusses field-based HTP, including the use of unoccupied aerial vehicles and light detection and ranging, as well as how we can achieve increased genetic gain by utilizing image data derived from HTP. Section 3 considers extending commonly used GS models to integrate HTP data as covariates associated with the principal trait response, such as yield. Particular focus is placed on single-trait, multi-trait, and genotype by environment interaction models. One unique aspect of HTP data is that phenomics platforms often produce large-scale data with high spatial and temporal resolution for capturing dynamic growth, development, and stress responses. Section 4 discusses the utility of a random regression model for performing longitudinal modeling. The chapter concludes with a discussion of some standing issues.

Unconditional and conditional analysis of epistasis between tillering QTLs based on single segment substitution lines in rice

Epistasis plays an important role in manipulating rice tiller number, but epistatic mechanism still remains a challenge. Here we showed the process of epistatic analysis between tillering QTLs. A half diallel mating scheme was conducted based on 6 single segment substitution lines and 9 dual segment pyramiding lines to allow the analysis of 4 epistatic components. Additive-additive, additive-dominance, dominance-additive, and dominance-dominance epistatic effects were estimated at 9 stages of development via unconditional QTL analysis simultaneously. Unconditional QTL effect (QTL cumulative effect before a certain stage) was then divided into several conditional QTL components (QTL net effect in a certain time interval). The results indicated that epistatic interaction was prevalent, all QTL pairs harboring epistasis and one QTL always interacting with other QTLs in various component ways. Epistatic effects were dynamic, occurring mostly within 14d and 21-35d after transplant and exhibited mainly negative effects. The genetic and developmental mechanism on several tillering QTLs was further realized and perhaps was useful for molecular pyramiding breeding and heterosis utilization for improving plant architecture.

Cruz PC, Baltazar AM, Diaz MGQ, Pang YL, Acero BL, Li Z. Identification of quantitative trait loci governing early germination and seedling vigor traits related to weed competitive ability in rice

Weed competitive ability (WCA) is vital for the improvement of grain yield under direct-seeded and aerobic rice ecosystems where weeds are a major limiting factor. Early seed germination (ESG) and early seedling vigor (ESV) are the crucial traits for WCA. This study attempted to map the quantitative trait loci (QTLs) and hotspot regions governing ESG and ESV traits. A total of 167 BC1F5 selective introgression lines developed from an early backcross population involving Weed Tolerant Rice 1 (WTR-1) as the recipient parent and Y-134 as the donor parent were phenotyped for ESG and ESV traits. Analysis of variance revealed significant differences in ESG-related traits except for root length and in ESV-related traits except for plant height at 7 days after sowing. A total of 677-high quality single nucleotide polymorphism (SNP) markers were used to analyze the marker-trait association from a 6 K SNP genotyping array. Forty-three QTLs were identified on all chromosomes, except on chromosomes 4 and 8. Thirty QTLs were contributed by a desirable allele from Y-134, whereas 13 QTLs were from WTR-1. Twenty-eight of the identified genetic loci associated with ESG and ESV traits were novel. Two QTL hotspot regions were mapped on chromosomes 11 and 12. The genomic regions of QTL hotspots were fine-tuned and a total of 13 putative candidate genes were discovered on chromosomes 11 and 12 collectively. The mapped QTLs will be useful in advancing the marker aided-selection schemes and breeding programs for the development of rice cultivars with WCA traits.

Increased Power and Accuracy of Causal Locus Identification in Time Series Genome-wide Association in Sorghum

The phenotypes of plants develop over time and change in response to the environment. New engineering and computer vision technologies track these phenotypic changes. Identifying the genetic loci regulating differences in the pattern of phenotypic change remains challenging. This study used functional principal component analysis (FPCA) to achieve this aim. Time series phenotype data were collected from a sorghum (Sorghum bicolor) diversity panel using a number of technologies including conventional color photography and hyperspectral imaging. This imaging lasted for 37 d and centered on reproductive transition. A new higher density marker set was generated for the same population. Several genes known to control trait variation in sorghum have been previously cloned and characterized. These genes were not confidently identified in genome-wide association analyses at single time points. However, FPCA successfully identified the same known and characterized genes. FPCA analyses partitioned the role these genes play in controlling phenotypes. Partitioning was consistent with the known molecular function of the individual cloned genes. These data demonstrate that FPCA-based genome-wide association studies can enable robust time series mapping analyses in a wide range of contexts. Moreover, time series analysis can increase the accuracy and power of quantitative genetic analyses.

The synthesis of strigolactone is affected by endogenous ascorbic acid in transgenic rice for l-galactono-1, 4-lactone dehydrogenase suppressed or overexpressing

Rice tillering, which determines the panicle number per plant, is an important agronomic trait for grain production. In higher plants, ascorbic acid (Asc) plays a major role in ROS-scavenging activity. l-Galactono-1, 4-lactone dehydrogenase (GalLDH, EC1.3.2.3) is an enzyme that catalyzes the last step of Asc biosynthesis in plants. Previously, we have reported that homozygous L-GalLDH-suppressed transgenic rice plants (GI) display a reduced tiller number and a lower level of foliar carotenoids (Car) compared with wild type. Strigolactones (SL), which play an important role in the suppression of shoot branching, are synthesized in the roots of rice plant using Car as substrates. In this paper, the relationship between Asc, SL, the accumulation of H₂O₂, changes in antioxidant capacity, enzyme activities, and gene transcriptions related to the synthesis of SL were analyzed in transgenic rice plants for L-GalLDH suppressed (GI-1 and GI-2) and overexpressing (GO-2). The results showed that the altered level of Asc in the L-GalLDH transgenic rice plants leads to a change in redox homeostasis, resulting in a marked accumulation of H₂O₂ and decreased antioxidant capacity in GI-1 and GI-2, but lower H₂O₂ content and increased antioxidant capacity in GO-2. Meanwhile, the altered level of Asc also leads to altered enzyme activities and gene transcript abundances related to SL synthesis in L-GalLDH transgenics. These observations support the conclusion that Asc influences tiller number in the L-GalLDH transgenics by affecting H₂O₂ accumulation and antioxidant capacity, and altering those enzyme activities and gene transcript abundances related to SL synthesis.

Mapping developmental QTL for plant height in soybean [Glycine max (L.) Merr.] using a four-way recombinant inbred line population

Plant height (PH) is an important trait in soybean, as taller plants may have higher yields but may also be at risk for lodging. Many genes act jointly to influence PH throughout development. To map the quantitative trait loci (QTL) controlling PH, we used the unconditional variable method (UVM) and conditional variable method (CVM) to analyze PH data for a four-way recombinant inbred line (FW-RIL) population derived from the cross of (Kenfeng14 × Kenfeng15) × (Heinong48 × Kenfeng19). We identified 7, 8, 16, 19, 15, 27, 17, 27, 22, and 24 QTL associated with PH at 10 developmental stages, respectively. These QTL mapped to 95 genomic regions. Among these QTL, 9 were detected using UVM and CVM, and 89 and 66 were only detected by UVM or CVM, respectively. In total, 36 QTL controlling PH were detected at multiple developmental stages and these made unequal contributions to genetic variation throughout development. Among 19 novel regions discovered in our study, 7 could explain over 10% of the phenotypic variation and contained only one single QTL. The unconditional and conditional QTL detected here could be used in molecular design breeding across the whole developmental procedure.

Non-dormant Axillary Bud 1 regulates axillary bud outgrowth in sorghum

Tillering contributes to grain yield and plant architecture and therefore is an agronomically important trait in sorghum (Sorghum bicolor). Here, we identified and functionally characterized a mutant of the Non-dormant Axillary Bud 1 (NAB1) gene from an ethyl methanesulfonate-mutagenized sorghum population. The nab1 mutants have increased tillering and reduced plant height. Map-based cloning revealed that NAB1 encodes a carotenoid-cleavage dioxygenase 7 (CCD7) orthologous to rice (Oryza sativa) HIGH-TILLERING DWARF1/DWARF17 and Arabidopsis thaliana MORE AXILLARY BRANCHING 3. NAB1 is primarily expressed in axillary nodes and tiller bases and NAB1 localizes to chloroplasts. The nab1 mutation causes outgrowth of basal axillary buds; removing these non-dormant basal axillary buds restored the wild-type phenotype. The tillering of nab1 plants was completely suppressed by exogenous application of the synthetic strigolactone analog GR24. Moreover, the nab1 plants had no detectable strigolactones and displayed stronger polar auxin transport than wild-type plants. Finally, RNA-seq showed that the expression of genes involved in multiple processes, including auxin-related genes, was significantly altered in nab1. These results suggest that NAB1 functions in strigolactone biosynthesis and the regulation of shoot branching via an interaction with auxin transport.

Shortened Basal Internodes Encodes a Gibberellin 2-Oxidase and Contributes to Lodging Resistance in Rice

Breeding semi-dwarf varieties to improve lodging resistance has been proven to be enormously successful in increasing grain yield since the advent of the "green revolution." However, the breeding of the majority of semi-dwarf rice varieties in Asia has been dependent mainly on genetic introduction of the mutant alleles of SD1, which encodes a gibberellin (GA) 20-oxidase, OsGA20ox2, for catalyzing GA biosynthesis. Here, we report a new rice lodging-resistance gene, Shortened Basal Internodes (SBI), which encodes a gibberellin 2-oxidase and specifically controls the elongation of culm basal internodes through deactivating GA activity. SBI is predominantly expressed in culm basal internodes. Genetic analyses indicate that SBI is a semi-dominant gene affecting rice height and lodging resistance. SBI allelic variants display different activities and are associated with the height of rice varieties. Breeding with higher activity of the SBI allele generates new rice varieties with improved lodging resistance and increased yield. The discovery of the SBI provides a desirable gene resource for producing semi-dwarf rice phenotypes and offers an effective strategy for breeding rice varieties with enhanced lodging resistance and high yield.

Time dependent genetic analysis links field and controlled environment phenotypes in the model C4 grass Setaria

Vertical growth of plants is a dynamic process that is influenced by genetic and environmental factors and has a pronounced effect on overall plant architecture and biomass composition. We have performed six controlled growth trials of an interspecific Setaria italica x Setaria viridis recombinant inbred line population to assess how the genetic architecture of plant height is influenced by developmental queues, water availability and planting density. The non-destructive nature of plant height measurements has enabled us to monitor height throughout the plant life cycle in both field and controlled environments. We find that plant height is reduced under water limitation and high density planting and affected by growth environment (field vs. growth chamber). The results support a model where plant height is a heritable, polygenic trait and that the major genetic loci that influence plant height function independent of growth environment. The identity and contribution of loci that influence height changes dynamically throughout development and the reduction of growth observed in water limited environments is a consequence of delayed progression through the genetic program which establishes plant height in Setaria. In this population, alleles inherited from the weedy S. viridis parent act to increase plant height early, whereas a larger number of small effect alleles inherited from the domesticated S. italica parent collectively act to increase plant height later in development.

Growth is a dynamic process that responds to a changing environment. Most of the methods that we have for measuring are static and collecting information throughout an organisms lifecycle is labor and cost prohibitive. Advances in imaging and robotics technology have enabled novel approaches to understanding how plants adapt to the environment. Using the model grass Setaria and new methods for measuring parameters from images, we investigate the genetic architecture of plant height in response to water availability and planting density. Height is one of the most influential components of plant architecture, determining tradeoffs between competition and resource allocation and is an important trait for boosting yields. The non-destructive nature of plant height measurements has enabled us to monitor growth throughout the plant life cycle in both field and controlled environments. We identified several loci controlling height in a population derived from a wild strain of Setaria viridis and its domesticated relative Setaria italica, as well as the developmental time in which these loci act. In this population, alleles inherited from the wild parent act to increase plant height early, whereas a larger number of small effect alleles inherited from the domesticated parent collectively act to increase plant height later in development.

Characterization of the temporal and spatial expression of wheat (Triticum aestivum L.) plant height at the QTL level and their influence on yield-related traits

The temporal and spatial expression patterns of stable QTL for plant height and their influences on yield were characterized. Plant height (PH) is a complex trait in wheat (Triticum aestivum L.) that includes the spike length (SL) and the internode lengths from the first to the fifth internode, which are counted from the top and abbreviated as FIRITL, SECITL, THIITL, FOUITL, and FIFITL, respectively. This study identified eight putative additive quantitative trait loci (QTL) for PH. In addition, unconditional and conditional QTL mapping were used to analyze the temporal and spatial expression patterns of five stable QTL for PH. qPh-3A mainly regulated SL, FIRITL, and FIFITL to affect PH during the booting-heading stage (BS-HS); qPh-3D regulated all internode lengths to affect PH, especially during the BS-HS; before HS, qPh-4B mainly affected FIRITL, SECITL, THIITL, and FOUITL and qPh-5A.1 mainly affected SECITL, THIITL, and FOUITL to regulate PH; and qPh-6B mainly regulated FIRITL to affect the PH after the booting stage (BS). qPhdv-4B, a QTL for the response of PH to nitrogen stress, was stable and co-localized with qPh-4B. All five stable QTL, except for qPh-3A, were related to the 1000 kernel weight and yield per plant. Regions of qPh-3A, qPh-3D, qPh-4B, qPh-5A.1, and qPh-6B showed synteny to parts of rice chromosomes 1, 1, 3, 9, and 2, respectively. Based on comparative genomics analysis, Rht-B1b was cloned and mapped in the CI of qPh-4B. This report provides useful information for fine mapping of the stable QTL for PH and the genetic improvement of wheat plant type.

Genetic variation of growth dynamics in maize (Zea mays L.) revealed through automated non-invasive phenotyping

Hitherto, most quantitative trait loci of maize growth and biomass yield have been identified for a single time point, usually the final harvest stage. Through this approach cumulative effects are detected, without considering genetic factors causing phase-specific differences in growth rates. To assess the genetics of growth dynamics, we employed automated non-invasive phenotyping to monitor the plant sizes of 252 diverse maize inbred lines at 11 different developmental time points; 50 k SNP array genotype data were used for genome-wide association mapping and genomic selection. The heritability of biomass was estimated to be over 71%, and the average prediction accuracy amounted to 0.39. Using the individual time point data, 12 main effect marker-trait associations (MTAs) and six pairs of epistatic interactions were detected that displayed different patterns of expression at various developmental time points. A subset of them also showed significant effects on relative growth rates in different intervals. The detected MTAs jointly explained up to 12% of the total phenotypic variation, decreasing with developmental progression. Using non-parametric functional mapping and multivariate mapping approaches, four additional marker loci affecting growth dynamics were detected. Our results demonstrate that plant biomass accumulation is a complex trait governed by many small effect loci, most of which act at certain restricted developmental phases. This highlights the need for investigation of stage-specific growth affecting genes to elucidate important processes operating at different developmental phases.

Molecular Dissection of Seedling Salinity Tolerance in Rice (Oryza sativa L.) Using a High-Density GBS-Based SNP Linkage Map

BACKGROUND

Salinity is one of the many abiotic stresses limiting rice production worldwide. Several studies were conducted to identify quantitative trait loci (QTLs) for traits associated to salinity tolerance. However, due to large confidence interval for the position of QTLs, utility of reported QTLs and the associated markers has been limited in rice breeding programs. The main objective of this study is to construct a high-density rice genetic map for identification QTLs and candidate genes for salinity tolerance at seedling stage.

RESULTS

We evaluated a population of 187 recombinant inbred lines (RILs) developed from a cross between Bengal and Pokkali for nine traits related to salinity tolerance. A total of 9303 SNP markers generated by genotyping-by-sequencing (GBS) were mapped to 2817 recombination points. The genetic map had a total map length of 1650 cM with an average resolution of 0.59 cM between markers. For nine traits, a total of 85 additive QTLs were identified, of which, 16 were large-effect QTLs and the rest were small-effect QTLs. The average interval size of QTL was about 132 kilo base pairs (Kb). Eleven of the 85 additive QTLs validated 14 reported QTLs for shoot potassium concentration, sodium-potassium ratio, salt injury score, plant height, and shoot dry weight. Epistatic QTL mapping identified several pairs of QTLs that significantly contributed to the variation of traits. The QTL for high shoot K⁺ concentration was mapped near the qSKC1 region. However, candidate genes within the QTL interval were a CC-NBS-LRR protein, three uncharacterized genes, and transposable elements. Additionally, many QTLs flanked small chromosomal intervals containing few candidate genes. Annotation of the genes located within QTL intervals indicated that ion transporters, osmotic regulators, transcription factors, and protein kinases may play essential role in various salt tolerance mechanisms.

CONCLUSION

The saturation of SNP markers in our linkage map increased the resolution of QTL mapping. Our study offers new insights on salinity tolerance and presents useful candidate genes that will help in marker-assisted gene pyramiding to develop salt tolerant rice varieties.

Dissimilar Manifestation of Heterosis in Superhybrid Rice at Early-Tillering Stage under Nutrient-Deficient and Nutrient-Sufficient Condition

Heterosis has long been exploited for crop breeding; however, the genetic mechanisms, particularly the initial establishment of heterosis during the early vegetative growth phase, remain elusive. The biggest challenge for that is to exclude noise genes from the identified heterosis-related candidates. Herein, we use nutrient-deficient hybrid with no measurable growth heterosis as control. After filtering these irrelevant genes, only 336 differentially expressed genes (DEGs), which is significantly lower than in previous reports, were identified as heterosis-related genes in a superhybrid rice of Liang-You-Pei-Jiu (LYP9) at early-tillering stage. Among the DEGs that could be mapped to quantitative trait loci (QTL), approximately 72.8% could be covered by yield or growth vigor-related QTL, thereby suggesting that our DEGs were reliable and may have potential value to rice breeding. Among the 336 DEGs identified, a majority showed intermediate expression relative to that of its parental lines (i.e. additive effects), particularly, expression was frequently more similar to the paternal line rather than the maternal line (44.1% versus 32.7%); the remaining 27.1% were exclusively up- or down-regulated between the hybrid and either parent. Interestingly, up-regulated genes encoded various enzymes, whereas down-regulated genes were enriched in responses to stress, indicating that hybrids may benefit from both activating metabolism-related pathways and alleviating fitness cost through allelic interactions. Furthermore, a significantly larger proportion of divergent genes and higher nonsynonymous substitutions rates were detected in these DEGs, suggesting a potential contribution to the establishment of heterosis in superhybrid rice LYP9.

Genome-Wide Association of Rice Blast Disease Resistance and Yield-Related Components of Rice

Robust disease resistance may require an expenditure of energy that may limit crop yield potential. In the present study, a subset of a United States Department of Agriculture rice core collection consisting of 151 accessions was selected using a major blast resistance (R) gene, Pi-ta, marker and was genotyped with 156 simple sequence repeat (SSR) markers. Disease reactions to Magnaporthe oryzae, the causal agent of rice blast disease, were evaluated under greenhouse and field conditions, and heading date, plant height, paddy and brown seed weight in two field environments were analyzed, using an association mapping approach. A total of 21 SSR markers distributed among rice chromosomes 2 to 12 were associated with blast resistance, and 16 SSR markers were associated with seed weight, heading date, and plant height. Most noticeably, shorter plants were significantly correlated with resistance to blast, rice genomes with Pi-ta were associated with lighter seed weights, and the susceptible alleles of RM171 and RM6544 were associated with heavier seed weight. These findings unraveled a complex relationship between disease resistance and yield-related components.

Dynamic and comparative QTL analysis for plant height in different developmental stages of Brassica napus L

This report describes a dynamic QTL analysis for plant height at various stages using a large doubled haploid population and performs a QTL comparison between different populations in Brassica napus. Plant height (PH) not only plays an important role in determining plant architecture, but is also an important character related to yield. The process of determining PH occurs through a series of steps; however, no studies have focused on developmental behavior factors affecting PH in Brassica napus. In the present study, KN DH, a large doubled haploid population containing 348 lines was used for a dynamic quantitative trait locus (QTL) analysis for PH in six experiments. In all, 20 QTLs were identified at maturity, whereas 50 QTLs were detected by conditional m apping method and the same number was identified by unconditional mapping strategies. Interestingly, five unconditional QTLs ucPH.A2-2, ucPH.A3-2, ucPH.C5-1, ucPH.C6-2 and ucPH.C6-3 were identified that were consistent over the all growth stages of one or two particular experiments, and one conditional QTL cPH.A2-3 was expressed throughout the entire growth process in one experiment. A total of 70 QTLs were obtained after combining QTLs identified at maturity, by conditional and unconditional mapping strategies, in which 25 showed opposite genetic effects in different periods/stages and experiments. A consensus map containing 1357 markers was constructed to compare QTLs identified in the KN population with five previously mapped populations. Alignment of the QTLs detected in different populations onto the consensus map showed that 27 were repeatedly detected in different genetic backgrounds. These findings will enhance our understanding of the genetic control of PH regulation in B. napus, and will be useful for rapeseed genetic manipulation through molecular marker-assisted selection.

Early seedling vigour, an imperative trait for direct-seeded rice: an overview on physio-morphological parameters and molecular markers

Rapid uniform germination and accumulation of biomass during initial phase of seedling establishment is an essential phenotypic trait considered as early seedling vigour for direct seeded situation in rice irrespective of environment. Enhanced role of carbohydrate, amylase, growth hormones, antioxidant enzymes and ascorbic acid brings changes in vigour and phenotype of seedling. Early establishment and demanding life form dominate the surroundings. Crop plant that has better growth overdrives the weed plant and suppresses its growth. Seedling early vigour is the characteristic of seed quality and describes the rapid, uniform germination and the establishment of strong seedlings in any environmental condition. The phenotype of modern rice varieties has been changed into adaptable for transplanted rice with thirst toward water and selection pressure for semi-dwarf architecture resulting in reduced early vigour. Decreasing freshwater availability and rising labour cost drives the search for a suitable alternative management system to enhance grain yield productivity for the burgeoning world population. In view of these issues, much attention has been focused on dry direct-seeded rice, because it demands low input. A rice cultivar with a strong seedling vigour trait is desirable in case of direct seeding. However, seedling vigour has not been selected in crop improvement programmes in conventional breeding due to its complex nature and quantitative inheritance. Molecular markers have been proven effective in increasing selection efficiency, particularly for quantitative traits that are simply inherited. Marker-assisted selection approach has facilitated efficient and precise transfer of genes/QTL(s) into many crop species and suggests a speedy and efficient technique over conventional breeding and selection methods. In this review, we present the findings and investigations in the field of seedling vigour in rice that includes the nature of inheritance of physio-morphological and biochemical traits and QTLs to assist plant breeders who work for direct-seeded rice.

Dynamic QTL and epistasis analysis on seedling root traits in upland cotton

Roots are involved in acquisition of water and nutrients, as well as in providing structural support to plant. The root system provides a dynamic model for developmental analysis. Here, we investigated quantitative trait loci (QTL), dynamic conditional QTL and epistatic interactions for seedling root traits using an upland cotton F2 population and a constructed genetic map. Totally, 37 QTLs for root traits, 35 dynamic conditional QTLs based on the net increased amount of root traits (root tips, forks, length, surface area and volume) (i) after transplanting 10 days compared to 5 days, and (ii) after transplanting 15 days to 10 days were detected. Obvious dynamic characteristic of QTL and dynamic conditional QTL existed at different developmental stages of root because QTL and dynamic conditional QTL had not been detected simultaneously. We further confirmed that additive and dominance effects of QTL qRSA-chr1-1 in interval time 5 to 10 DAT (days after transplant) offset the effects in 10 to 15 DAT. Lots of two-locus interactions for root traits were identified unconditionally or dynamically, and a few epistatic interactions were only detected simultaneously in interval time of 5-10 DAT and 10-15 DAT, suggesting different interactive genetic mechanisms on root development at different stages. Dynamic conditional QTL and epistasis effects provide new attempts to understand the dynamics of roots and provide clues for root architecture selection in upland cotton.

Adult plant development in triticale (× triticosecale wittmack) is controlled by dynamic genetic patterns of regulation

Many biologically and agronomically important traits are dynamic and show temporal variation. In this study, we used triticale (× Triticosecale Wittmack) as a model crop to assess the genetic dynamics underlying phenotypic plasticity of adult plant development. To this end, a large mapping population with 647 doubled haploid lines derived from four partially connected families from crosses among six parents was scored for developmental stage at three different time points. Using genome-wide association mapping, we identified main effect and epistatic quantitative trait loci (QTL) at all three time points. Interestingly, some of these QTL were identified at all time points, whereas others appear to only contribute to the genetic architecture at certain developmental stages. Our results illustrate the temporal contribution of QTL to the genetic control of adult plant development and more generally, the temporal genetic patterns of regulation that underlie dynamic traits.

Generalized linear model for mapping discrete trait loci implemented with LASSO algorithm

Generalized estimating equation (GEE) algorithm under a heterogeneous residual variance model is an extension of the iteratively reweighted least squares (IRLS) method for continuous traits to discrete traits. In contrast to mixture model-based expectation–maximization (EM) algorithm, the GEE algorithm can well detect quantitative trait locus (QTL), especially large effect QTLs located in large marker intervals in the manner of high computing speed. Based on a single QTL model, however, the GEE algorithm has very limited statistical power to detect multiple QTLs because of ignoring other linked QTLs. In this study, the fast least absolute shrinkage and selection operator (LASSO) is derived for generalized linear model (GLM) with all possible link functions. Under a heterogeneous residual variance model, the LASSO for GLM is used to iteratively estimate the non-zero genetic effects of those loci over entire genome. The iteratively reweighted LASSO is therefore extended to mapping QTL for discrete traits, such as ordinal, binary, and Poisson traits. The simulated and real data analyses are conducted to demonstrate the efficiency of the proposed method to simultaneously identify multiple QTLs for binary and Poisson traits as examples.

Identification and fine mapping of quantitative trait loci for seed vigor in germination and seedling establishment in rice

Seed vigor is an index of seed quality that is used to describe the rapid and uniform germination and the establishment of strong seedlings in any environmental conditions. Strong seed vigor in low-temperature germination conditions is particularly important in direct-sowing rice production systems. However, seed vigor has not been selected as an important breeding trait in traditional breeding programs due to its quantitative inherence. In this study, we identified and mapped eight quantitative trait loci (QTLs) for seed vigor by using a recombinant inbred population from a cross between rice (Oryza sativa L. ssp. indica) cultivars ZS97 and MH63. Conditional QTL analysis identified qSV-1, qSV-5b, qSV-6a, qSV-6b, and qSV-11 influenced seedling establishment and that qSV-5a, qSV-5c, and qSV-8 influenced only germination. Of these, qSV-1, qSV-5b, qSV-6a, qSV-6b, and qSV-8 were low-temperature-specific QTLs. Two major-effective QTLs, qSV-1, and qSV-5c were narrowed down to 1.13-Mbp and 400-kbp genomic regions, respectively. The results provide tightly linked DNA markers for the marker-assistant pyramiding of multiple positive alleles for increased seed vigor in both normal and low-temperature germination environments.

Conditional and unconditional QTL mapping of drought-tolerance-related traits of wheat seedling using two related RIL populations

For discovering the quantitative trait loci (QTLs) contributing to early seedling growth and drought tolerance during germination, conditional and unconditional analyses of 12 traits of wheat seedlings: coleoptile length, seedling height, longest root length, root number, seedling fresh weight, stem and leaves fresh weight, root fresh weight, seedling dry weight, stem and leaves dry weight, root dry weight, root to shoot fresh weight ratio, root-to-shoot dry weight ratio, were conducted under two water conditions using two F8:9 recombinant inbred line (RIL) populations. The results of unconditional analysis are as follows: 88 QTLs accounting for 3.33-77.01% of the phenotypic variations were detected on chromosomes 1A, 1B, 1D, 2A, 2B, 2D, 3A, 3B, 4A, 4B, 4D, 5A, 5B, 5D, 6A, 6B, 6D, 7A, 7B and 7D. Among these QTLs, 19 were main-effect QTLs with a contribution rate greater than 10%. The results of the conditional QTL analysis of 12 traits under osmotic stress on normal water conditions were as follows: altogether 22 QTLs concerned with drought tolerance were detected on chromosomes 1B, 2A, 2B, 3B, 4A, 5D, 6A, 6D, 7B, and 7D. Of these QTLs, six were main-effect QTLs. These 22 QTLs were all special loci directly concerned with drought tolerance and most of them could not be detected by unconditional analysis. The finding of these QTLs has an important significance for fine-mapping technique, map-based cloning, and molecular marker-assisted selection of early seedling traits, such as growth and drought tolerance.

Genetic dynamics underlying phenotypic development of biomass yield in triticale

BACKGROUND

The nature of dynamic traits with their phenotypic plasticity suggests that they are under the control of a dynamic genetic regulation. We employed a precision phenotyping platform to non-invasively assess biomass yield in a large mapping population of triticale at three developmental stages.

RESULTS

Using multiple-line cross QTL mapping we identified QTL for each of these developmental stages which explained a considerable proportion of the genotypic variance. Some QTL were identified at each developmental stage and thus contribute to biomass yield throughout the studied developmental phases. Interestingly, we also observed QTL that were only identified for one or two of the developmental stages illustrating a temporal contribution of these QTL to the trait. In addition, epistatic QTL were detected and the epistatic interaction landscape was shown to dynamically change with developmental progression.

CONCLUSIONS

In summary, our results reveal the temporal dynamics of the genetic architecture underlying biomass accumulation in triticale and emphasize the need for a temporal assessment of dynamic traits.

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.

Mapping dynamic QTL for plant height in triticale

BACKGROUND

Plant height is a prime example of a dynamic trait that changes constantly throughout adult development. In this study we utilised a large triticale mapping population, comprising 647 doubled haploid lines derived from 4 families, to phenotype for plant height by a precision phenotyping platform at multiple time points.

RESULTS

Using multiple-line cross QTL mapping we identified main effect and epistatic QTL for plant height for each of the time points. Interestingly, some QTL were detected at all time points whereas others were specific to particular developmental stages. Furthermore, the contribution of the QTL to the genotypic variance of plant height also varied with time as exemplified by a major QTL identified on chromosome 6A.

CONCLUSIONS

Taken together, our results in the small grain cereal triticale reveal the importance of considering temporal genetic patterns in the regulation of complex traits such as plant height.

A model framework for identifying genes that guide the evolution of heterochrony

Heterochrony, the phylogenic change in the time of developmental events or rate of development, has been thought to play an important role in producing phenotypic novelty during evolution. Increasing evidence suggests that specific genes are implicated in heterochrony, guiding the process of developmental divergence, but no quantitative models have been instrumented to map such heterochrony genes. Here, we present a computational framework for genetic mapping by which to characterize and locate quantitative trait loci (QTLs) that govern heterochrony described by four parameters, the timing of the inflection point, the timing of maximum acceleration of growth, the timing of maximum deceleration of growth, and the length of linear growth. The framework was developed from functional mapping, a dynamic model derived to map QTLs for the overall process and pattern of development. By integrating an optimality algorithm, the framework allows the so-called heterochrony QTLs (hQTLs) to be tested and quantified. Specific pipelines are given for testing how hQTLs control the onset and offset of developmental events, the rate of development, and duration of a particular developmental stage. Computer simulation was performed to examine the statistical properties of the model and demonstrate its utility to characterize the effect of hQTLs on population diversification due to heterochrony. By analyzing a genetic mapping data in rice, the framework identified an hQTL that controls the timing of maximum growth rate and duration of linear growth stage in plant height growth. The framework provides a tool to study how genetic variation translates into phenotypic innovation, leading a lineage to evolve, through heterochrony.

Mixed linear model approach for mapping quantitative trait loci underlying crop seed traits

The crop seed is a complex organ that may be composed of the diploid embryo, the triploid endosperm and the diploid maternal tissues. According to the genetic features of seed characters, two genetic models for mapping quantitative trait loci (QTLs) of crop seed traits are proposed, with inclusion of maternal effects, embryo or endosperm effects of QTL, environmental effects and QTL-by-environment (QE) interactions. The mapping population can be generated either from double back-cross of immortalized F2 (IF2) to the two parents, from random-cross of IF2 or from selfing of IF2 population. Candidate marker intervals potentially harboring QTLs are first selected through one-dimensional scanning across the whole genome. The selected candidate marker intervals are then included in the model as cofactors to control background genetic effects on the putative QTL(s). Finally, a QTL full model is constructed and model selection is conducted to eliminate false positive QTLs. The genetic main effects of QTLs, QE interaction effects and the corresponding P-values are computed by Markov chain Monte Carlo algorithm for Gaussian mixed linear model via Gibbs sampling. Monte Carlo simulations were performed to investigate the reliability and efficiency of the proposed method. The simulation results showed that the proposed method had higher power to accurately detect simulated QTLs and properly estimated effect of these QTLs. To demonstrate the usefulness, the proposed method was used to identify the QTLs underlying fiber percentage in an upland cotton IF2 population. A computer software, QTLNetwork-Seed, was developed for QTL analysis of seed traits.

Conditional QTL mapping of protein content in wheat with respect to grain yield and its components

Grain protein content in wheat (Triticum aestivum L.) is generally considered a highly heritable character that is negatively correlated with grain yield and yield-related traits. Quantitative trait loci (QTL) for protein content was mapped using data on protein content and protein content conditioned on the putatively interrelated traits to evaluate possible genetic interrelationships between protein content and yield, as well as yield-related traits. Phenotypic data were evaluated in a recombinant inbred line population with 302 lines derived from a cross between the Chinese cultivar Weimai 8 and Luohan 2. Inclusive composite interval mapping using IciMapping 3.0 was employed for mapping unconditional and conditional QTL with additives. A strong genetic relationship was found between protein content and grain yield, and yield-related traits. Unconditional QTL mapping analysis detected seven additive QTL for protein content, with additive effects ranging in absolute size from 0.1898% to 0.3407% protein content, jointly accounting for 43.45% of the trait variance. Conditional QTL mapping analysis indicated two QTL independent from yield, which can be used in marker-assisted selection for increasing yield without affecting grain protein content. Three additional QTL with minor effects were identified in the conditional mapping. Of the three QTLs, two were identified when protein content was conditioned on yield, which had pleiotropic effects on those two traits. Conditional QTL mapping can be used to dissect the genetic interrelationship between two traits at the individual QTL level for closely correlated traits. Further, conditional QTL mapping can reveal additional QTL with minor effects that are undetectable in unconditional mapping.

Dynamic QTL analysis and candidate gene mapping for waterlogging tolerance at maize seedling stage

Soil waterlogging is one of the major abiotic stresses adversely affecting maize growth and yield. To identify dynamic expression of genes or quantitative trait loci (QTL), QTL associated with plant height, root length, root dry weight, shoot dry weight and total dry weight were identified via conditional analysis in a mixed linear model and inclusive composite interval mapping method at three respective periods under waterlogging and control conditions. A total of 13, 19 and 23 QTL were detected at stages 3D|0D (the period during 0-3 d of waterlogging), 6D|3D and 9D|6D, respectively. The effects of each QTL were moderate and distributed over nine chromosomes, singly explaining 4.14-18.88% of the phenotypic variation. Six QTL (ph6-1, rl1-2, sdw4-1, sdw7-1, tdw4-1 and tdw7-1) were identified at two consistent stages of seedling development, which could reflect a continuous expression of genes; the remaining QTL were detected at only one stage. Thus, expression of most QTL was influenced by the developmental status. In order to provide additional evidence regarding the role of corresponding genes in waterlogging tolerance, mapping of Expressed Sequence Tags markers and microRNAs were conducted. Seven candidate genes were observed to co-localize with the identified QTL on chromosomes 1, 4, 6, 7 and 9, and may be important candidate genes for waterlogging tolerance. These results are a good starting point for understanding the genetic basis for selectively expressing of QTL in different stress periods and the common genetic control mechanism of the co-localized traits.

Conditional and unconditional QTL mapping of drought-tolerance-related traits of wheat seedling using two related RIL populations

For discovering the quantitative trait loci (QTLs) contributing to early seedling growth and drought tolerance during germination, conditional and unconditional analyses of 12 traits of wheat seedlings: coleoptile length, seedling height, longest root length, root number, seedling fresh weight, stem and leaves fresh weight, root fresh weight, seedling dry weight, stem and leaves dry weight, root dry weight, root to shoot fresh weight ratio, root-to-shoot dry weight ratio, were conducted under two water conditions using two F8:9 recombinant inbred line (RIL) populations. The results of unconditional analysis are as follows: 88 QTLs accounting for 3.33-77.01% of the phenotypic variations were detected on chromosomes 1A, 1B, 1D, 2A, 2B, 2D, 3A, 3B, 4A, 4B, 4D, 5A, 5B, 5D, 6A, 6B, 6D, 7A, 7B and 7D. Among these QTLs, 19 were main-effect QTLs with a contribution rate greater than 10%. The results of the conditional QTL analysis of 12 traits under osmotic stress on normal water conditions were as follows: altogether 22 QTLs concerned with drought tolerance were detected on chromosomes 1B, 2A, 2B, 3B, 4A, 5D, 6A, 6D, 7B, and 7D. Of these QTLs, six were main-effect QTLs. These 22 QTLs were all special loci directly concerned with drought tolerance and most of them could not be detected by unconditional analysis. The finding of these QTLs has an important significance for fine-mapping technique, map-based cloning, and molecular marker-assisted selection of early seedling traits, such as growth and drought tolerance.

Efficacy of pyramiding elite alleles for dynamic development of plant height in common wheat

Plant height is an important botanical feature closely related to yield. Two populations consisting of 118 and 262 accessions respectively were used to identify elite alleles for plant height and to validate their allelic effects. Plant height was measured from the early booting to the flowering stages. Simple sequence repeat markers for candidate quantitative trait locus (QTL) regions with large effects identified in a doubled haploid (DH) population (Hanxuan 10 × Lumai 14) were selected for further verification by association analysis. Nine loci significantly (P < 0.001) associated with plant height were detected 13 times in the population with 118 accessions. Three loci (Xgwm11-1B, Xwmc349-4B and Xcfd23-4D) were identified in three, two and two periods of plant height growth, respectively. Markers Xbarc168-2D, Xgwm249-2D, Xwmc349-4B, Xcfd23-4D and Xgwm410-5A located at or near additive QTL regions in the DH population proved to coincide with known Rht loci. The results showed a consistency between linkage analysis and association mapping, and also confirmed the value of fine mapping of QTL through combined linkage and association analyses. For final plant height, the alleles Xgwm11-1B₂₀₈ , Xwmc349-4B₁₀₃ and Xcfd23-4D₂₀₂ exhibited negative effects, i.e. reducing plant height; Xwmc349-4B₁₀₁ and Xcfd23-4D₂₀₅ showed significant positive effects. A second larger population (262 accessions) was used to validate the effects of these large-effect alleles and the efficacy of pyramiding in eight environments (year × site × water regime combinations). Strong correlations between final plant height and numbers of large-effect alleles indicated that the alleles contributed additively to plant height. The additive effects showed that pyramiding elite alleles for target traits has significant potential for wheat breeding.

Genetic dissection of developmental behavior of grain weight in wheat under diverse temperature and water regimes

As a quantitatively inherited trait related to high yield potential, grain weight (GW) development in wheat is constrained by abiotic stresses such as limited water supply and high temperature. Data from a doubled haploid population, derived from a cross of (Hanxuan 10 × Lumai 14), grown in four environments were used to explore the genetic basis of GW developmental behavior in unconditional and conditional quantitative trait locus (QTL) analyses using a mixed linear model. Thirty additive QTLs and 41 pairs of epistatic QTLs were detected, and were more frequently observed on chromosomes 1B, 2A, 2D, 4A, 4B and 7B. No single QTL was continually active during all stages or periods of grain growth. The QTLs with additive effects (A-QTLs) expressed in the period S1|S0 (the period from the flowering to the seventh day after) formed a foundation for GW development. GW development at these stages can be used as an index for screening superior genotypes under diverse abiotic stresses in a wheat breeding program. One QTL, i.e. Qgw.cgb-6A.2, showed high adaptability for water-limited and heat-stress environments. Many A-QTLs interacted with more than one other QTL in the two genetic models, such as Qgw.cgb-4B.2 interacted with five QTLs, showing that the genetic architecture underlying GW development involves a collective expression of genes with additive and epistatic effects.

A statistical procedure to map high-order epistasis for complex traits

Genetic interactions or epistasis have been thought to play a pivotal role in shaping the formation, development and evolution of life. Previous work focused on lower-order interactions between a pair of genes, but it is obviously inadequate to explain a complex network of genetic interactions and pathways. We review and assess a statistical model for characterizing high-order epistasis among more than two genes or quantitative trait loci (QTLs) that control a complex trait. The model includes a series of start-of-the-art standard procedures for estimating and testing the nature and magnitude of QTL interactions. Results from simulation studies and real data analysis warrant the statistical properties of the model and its usefulness in practice. High-order epistatic mapping will provide a routine procedure for charting a detailed picture of the genetic regulation mechanisms underlying the phenotypic variation of complex traits.

Conditional QTL underlying resistance to late blight in a diploid potato population

A large number of quantitative trait loci (QTL) for resistance to late blight of potato have been reported with a "conventional" method in which each phenotypic trait reflects the cumulative genetic effects for the duration of the disease process. However, as genes controlling response to disease may have unique contributions with specific temporal features, it is important to consider the phenotype as dynamic. Here, using the net genetic effects evidenced at consecutive time points during disease development, we report the first conditional mapping of QTL underlying late blight resistance in potato under five environments in Peru. Six conditional QTL were mapped, one each on chromosome 2, 7 and 12 and three on chromosome 9. These QTL represent distinct contributions to the phenotypic variation at different stages of disease development. By comparison, when conventional mapping was conducted, only one QTL was detected on chromosome 9. This QTL was the same as one of the conditional QTL. The results imply that conditional QTL reflect genes that function at particular stages during the host-pathogen interaction. The dynamics revealed by conditional QTL mapping could contribute to the understanding of the molecular mechanism of late blight resistance and these QTL could be used to target genes for marker development or manipulation to improve resistance.

Efficacy of microarray profiling data combined with QTL mapping for the identification of a QTL gene controlling the initial growth rate in rice

Seedling vigor, which is controlled by many quantitative trait loci (QTLs), is one of several important agronomic traits for direct-seedling rice systems. However, isolating these QTL genes is laborious and expensive. Here, we combined QTL mapping and microarray profiling to identify QTL genes for seedling vigor. By performing QTL mapping using 82 backcross inbred lines (BILs) of the Koshihikari (japonica) and Habataki (indica) cultivars for the rice initial growth, we identified two QTLs, early-stage plant development1/2 (qEPD1 and qEPD2), whose Koshihikari alleles promote plant height and/or leaf sheath length. Phenotypic analysis of the two substituted lines carrying the Habataki qEPD1 or qEPD2 allele revealed that qEPD2 functioned more dominantly for the initial growth of rice. From the microarray experiment, 55 and 45 candidate genes were found in the qEPD1 and qEPD2 genomic regions, which are expressed differentially between each substitution line (SL) and Koshihikari. Gibberellin 20 oxidase-2 (OsGA20ox2), which is identical to Semi Dwarf1 (SD1), was included among the 55 candidate genes of qEPD1, whereas its paralog, OsGA20ox1, was included among the 45 candidate genes of qEPD2. Consistently, introduction of the Koshihikari OsGA20ox1 allele into SL(qEPD2) increaseed its plant height and leaf sheath length significantly relative to the introduction of the Habataki OsGA20ox1 allele. Therefore, microarray profiling could be useful for rapidly screening QTL candidate genes. We concluded that OsGA20ox1 and OsGA20ox2 (SD1) function during the initial growth of rice, but OsGA20ox1 plays a dominant role in increasing plant height and leaf sheath length at the initial growth stage.

Association mapping of dynamic developmental plant height in common wheat

Drought as a major abiotic stress often occurs from stem elongation to the grain filling stage of wheat in northern China. Plant height (PH) is a suitable trait to model the dissection of drought tolerance. The purposes of the present study were to validate molecular markers for PH developmental behavior and identify elite alleles of molecular markers. After the phenotyping of 154 accessions for PH dynamic development under well-watered (WW) and drought stressed (DS) conditions, and the genotyping of 60 SSR markers from six candidate chromosome regions related to PH found in our previous linkage mapping studies, both parameters PH and drought tolerance coefficient (DTC) calculated by the conditional analysis were used for association mapping. A total of 46 significant association signals (P < 0.01) were identified in 23 markers, and phenotypic variation ranged from 7 to 50%. Among them, four markers Xgwm261-2D, Xgwm495-4B, Xbarc109-4B and Xcfd23-4D were detected under both water regimes. Furthermore, 10 markers were associated with DTC, and four with both parameters PH and DTC at the same plant development stage. The results revealed different allelic effects of associated markers; for example, the 155 bp Xgwm495-4B allele was associated with a reduced height of -11.2 cm under DS and -15.3 cm under WW, whereas the 167 bp allele exhibited increased height effects of 3.9 and 8.1 cm, respectively. This study demonstrates a strong power of joint association analysis and linkage mapping for the identification of important genes in wheat.

Dynamic QTL analysis of linolenic acid content in different developmental stages of soybean seed

Linolenic acid (LN) in soybean (Glycine max L. Merr.) seed mainly contributes to the undesirable odors and flavors commonly associated with poor oil quality. LN deposition at various stages of soybean seed development had not been reported by 2010. The objects of this study were (1) to identify and measure quantitative trait loci (QTL) underlying LN content and (2) to estimate the QTL effects expressed from earlier seed developmental stages to drying seed of soybean. One hundred and twenty-five F(5:8) and F(5:9) recombinant inbred lines derived from the cross of soybean cultivars 'Hefeng 25' and 'Dongnong L5' were used for the identification of QTL underlying LN content from the 37 day (D) to 86D stages after flowering, at Harbin in 2008 and 2009. QTL × Environment interactions (QE) effects were evaluated using a mixed genetic model (Zhu in J Zhejiang Univ (Natural Science) 33:327-335, 1999). Twelve unconditional QTL and 12 conditional QTL associated with LN content were identified at different developmental stages. Most of the QTL explained <10% of phenotypic variation of LN content. Unconditional QTL QLNF-1, QLNC2-1, QLND1b-1, QLNA2-1 and QLNH-1 influenced LN content across different development stages and environments. Conditional QTL QLNF-1, QLNC2-1 and QLNH-1 were identified in multiple developmental stages and environments. Conditional and unconditional QTL clustered in neighboring intervals on linkage groups A2, C2 and D1b. Ten QTL with conditional additive main effects (a) and/or conditional additive × environment interaction effects (ae) at specific developmental stage were identified on nine linkage groups. Of them, six QTL only possessed additive main effects and seven QTL had significant ae effects in different developmental stages. A total of 13 epistatic pairwise QTL were identified by conditional mapping in different developmental stages. Two pairs of QTL only showed aa effects and five pairs of QTL only showed aae effects at different developmental stages. QTL with aa effects, as well as their environmental interaction effects, appeared to vary at different developmental stages.

A functional-structural model of rice linking quantitative genetic information with morphological development and physiological processes

BACKGROUND AND AIMS

Although quantitative trait loci (QTL) analysis of yield-related traits for rice has developed rapidly, crop models using genotype information have been proposed only relatively recently. As a first step towards a generic genotype-phenotype model, we present here a three-dimensional functional-structural plant model (FSPM) of rice, in which some model parameters are controlled by functions describing the effect of main-effect and epistatic QTLs.

METHODS

The model simulates the growth and development of rice based on selected ecophysiological processes, such as photosynthesis (source process) and organ formation, growth and extension (sink processes). It was devised using GroIMP, an interactive modelling platform based on the Relational Growth Grammar formalism (RGG). RGG rules describe the course of organ initiation and extension resulting in final morphology. The link between the phenotype (as represented by the simulated rice plant) and the QTL genotype was implemented via a data interface between the rice FSPM and the QTLNetwork software, which computes predictions of QTLs from map data and measured trait data.

KEY RESULTS

Using plant height and grain yield, it is shown how QTL information for a given trait can be used in an FSPM, computing and visualizing the phenotypes of different lines of a mapping population. Furthermore, we demonstrate how modification of a particular trait feeds back on the entire plant phenotype via the physiological processes considered.

CONCLUSIONS

We linked a rice FSPM to a quantitative genetic model, thereby employing QTL information to refine model parameters and visualizing the dynamics of development of the entire phenotype as a result of ecophysiological processes, including the trait(s) for which genetic information is available. Possibilities for further extension of the model, for example for the purposes of ideotype breeding, are discussed.

Unconditional and conditional QTL mapping for the developmental behavior of tiller number in rice (Oryza sativa L.)

A single segment substitution population of 26 lines and their recipient parent Hua-jing-xian 74 (HJX74) were selected as experimental materials for analyzing the developmental behavior of tiller number in rice. By the unconditional QTL (quantitative trait locus) mapping method, a total number of 14 SSSLs were detected with QTLs controlling rice tiller number. The number of QTLs significantly affecting tiller number and their effect values estimated differed across measuring stages. More QTLs could be detected based on time-dependent measures of different stages. By the conditional QTL mapping method, it is possible to reveal net expression of gene in a time interval. 14 QTLs on tiller number expressed their effects in dynamic patterns of themselves during whole ontogeny. They exhibited mainly negative effects within 7 days after transplanting. During 7-21 days, QTLs were in active status and expressed larger positive effects. In the mid-period of 21-35 days, they had opposite genetic effects to wither tillers. Since then these QTLs expressed positive effects again to cause the appearance of noneffective tillers. The dynamics of QTL effects was in agreement with the actual change of tillers. Mapping QTL combining unconditional with conditional analysis for time-dependent measures is helpful to understand roundly the genetic bases for the development of quantitative traits.

Genetic dissection of the developmental behaviours of plant height in wheat under diverse water regimes

Plant height (PH), a crucial trait related to yield potential in crop plants, is known to be typically quantitatively inherited. However, its full expression can be inhibited by a limited water supply. In this study, the genetic basis of the developmental behaviour of PH was assessed in a 150-line wheat (Triticum aestivum L.) doubled haploid population (Hanxuan 10 x Lumai 14) grown in 10 environments (year x site x water regime combinations) by unconditional and conditional quantitative trait locus (QTL) analyses in a mixed linear model. Genes that were expressed selectively during ontogeny were identified. No single QTL was continually active in all periods of PH growth, and QTLs with additive effects (A-QTLs) expressed in the period S1|S0 (the period from the original point to the jointing stage) formed a foundation for PH development. Additive main effects (a effects), which were mostly expressed in S1|S0, were more important than epistatic main effects (aa effects) or QTL x environment interaction (QE) effects, suggesting that S1|S0 was the most significant development period affecting PH growth. A few QTLs, such as QPh.cgb-6B.7, showed high adaptability for water-limited environments. Many QTLs, including four A-QTLs (QPh.cgb-2D.1, QPh.cgb-4B.1, QPh.cgb-4D.1, and QPh.cgb-5A.7) coincident with previously identified reduced height (Rht) genes (Rht8, Rht1, Rht2, and Rht9), interacted with more than one other QTL, indicating that the genetic architecture underlying PH development is a network of genes with additive and epistatic effects. Therefore, based on multilocus combinations in S1|S0, superior genotypes were predicted for guiding improvements in breeding for PH.

Bayesian analysis for genetic architecture of dynamic traits

The dissection of the genetic architecture of quantitative traits, including the number and locations of quantitative trait loci (QTL) and their main and epistatic effects, has been an important topic in current QTL mapping. We extend the Bayesian model selection framework for mapping multiple epistatic QTL affecting continuous traits to dynamic traits in experimental crosses. The extension inherits the efficiency of Bayesian model selection and the flexibility of the Legendre polynomial model fitting to the change in genetic and environmental effects with time. We illustrate the proposed method by simultaneously detecting the main and epistatic QTLs for the growth of leaf age in a doubled-haploid population of rice. The behavior and performance of the method are also shown by computer simulation experiments. The results show that our method can more quickly identify interacting QTLs for dynamic traits in the models with many numbers of genetic effects, enhancing our understanding of genetic architecture for dynamic traits. Our proposed method can be treated as a general form of mapping QTL for continuous quantitative traits, being easier to extend to multiple traits and to a single trait with repeat records.