Genome-wide association study of leaf photosynthesis using a high-throughput gas exchange system in rice.

Enhancing leaf photosynthetic capacity is essential for improving the yield of rice (Oryza sativa L.). Although the exploitation of natural genetic resources is considered a promising approach to enhance photosynthetic capacity, genomic factors related to the genetic diversity of leaf photosynthetic capacity have yet to be fully elucidated due to the limitation of measurement efficiency. In this study, we aimed to identify novel genomic regions for the net CO2 assimilation rate (A) by combining genome-wide association study (GWAS) and the newly developed rapid closed gas exchange system MIC-100. Using three MIC-100 systems in the field at the vegetative stage, we measured A of 168 temperate japonica rice varieties with six replicates for three years. We found that the modern varieties exhibited higher A than the landraces, while there was no significant relationship between the release year and A among the modern varieties. Our GWAS scan revealed two major peaks located on chromosomes 4 and 8, which were repeatedly detected in the different experiments and in the generalized linear modelling approach. We suggest that high-throughput gas exchange measurements combined with GWAS is a reliable approach for understanding the genetic mechanisms underlying photosynthetic diversities in crop species.

Identification of Genomic Regions for Deep-Water Resistance in Rice for Efficient Weed Control with Reduced Herbicide Use.

Deep-water (DW) management in rice fields is a promising technique for efficient control of paddy weeds with reduced herbicide use. Maintaining a water depth of 10-20 cm for several weeks can largely suppress the weed growth, though it also inhibits rice growth because the DW management is usually initiated immediately after transplanting. Improving the DW resistance of rice during the initial growth stage is essential to avoid suppressing growth. In this study, we demonstrate a large genetic variation in the above-ground biomass (AGB) after the end of DW management among 165 temperate japonica varieties developed in Japan. Because the AGB closely correlated with plant length (PL) and tiller number (TN) at the early growth stage, we analyzed genomic regions associated with PL and TN by conducting a genome-wide association study. For PL, a major peak was detected on chromosome 3 (qPL3), which includes a gene encoding gibberellin biosynthesis, OsGA20ox1. The rice varieties with increased PL had a higher expression level of OsGA20ox1 as reported previously. For TN, a major peak was detected on chromosome 4 (qTN4), which includes NAL1 gene associated with leaf morphological development and panicle number. Although there was less difference in the expression level of NAL1 between genotypes, our findings suggest that an amino acid substitution in the exon region is responsible for the phenotypic changes. We also found that the rice varieties having alternative alleles of qPL3 and qTN4 showed significantly higher AGB than the varieties with the reference alleles. Our results suggest that OsGA20ox1 and NAL1 are promising genes for improving DW resistance in rice.

Identification of Novel Quantitative Trait Loci for Culm Thickness of Rice Derived from Strong-Culm Landrace in Japan, Omachi.

Increasing the lodging resistance of rice through genetic improvement has been an important target in breeding. To further enhance the lodging resistance of high-yielding rice varieties amidst climate change, it is necessary to not only shorten culms but strengthen them as well. A landrace rice variety, Omachi, which was established more than 100 years ago, has the largest culm diameter and bending moment at breaking in the basal internodes among 135 temperate japonica accessions. Using unused alleles in such a landrace is an effective way to strengthen the culm. In this study, we performed quantitative trait locus (QTL) analysis to identify the genetic factors of culm strength of Omachi using recombinant inbred lines (RILs) derived from a cross between Omachi and Koshihikari, a standard variety in Japan. We identified three QTLs for the culm diameter of the 5th internode on chromosomes 3 (qCD3) and 7 (qCD7-1, qCD7-2). Among them, qCD7-2 was verified by QTL analysis using the F2 population derived from a cross between one of the RILs and Koshihikari. RNA-seq analysis of shoot apex raised 10 candidate genes underlying the region of qCD7-2. The increase in culm strength by accumulating Omachi alleles of qCD3, qCD7-1 and qCD7-2 was 25.0% in 2020. These QTLs for culm diameter pleiotropically increased spikelet number per panicle but did not affect days to heading or culm length. These results suggest that the Omachi alleles of qCD3, qCD7-1 and qCD7-2 are useful for breeding to increase lodging resistance and yield.

Metabolomic analysis of rice brittle culm mutants reveals each mutant- specific metabolic pattern in each organ.

INTRODUCTION: Plant cell walls play an important role in providing physical strength and defence against abiotic stress. Rice brittle culm (bc) mutants are a strength-decreased mutant because of abnormal cell walls, and it has been reported that the causative genes of bc mutants affect cell wall composition. However, the metabolic alterations in each organ of bc mutants have remained unknown. OBJECTIVES: To evaluate the metabolic changes in rice bc mutants, comparative analysis of the primary metabolites was conducted. METHODS: The primary metabolites in leaves, internodes, and nodes of rice bc mutants and wild-type control were measured using CE- and LC-MS/MS. Multivariate analyses using metabolomic data was performed. RESULTS: We found that mutations in each bc mutant had different effects on metabolism. For example, higher oxalate content was observed in bc3 and bc1 bc3 mutants, suggesting that surplus carbon that was not used for cell wall components might be used for oxalate synthesis. In addition, common metabolic alterations such as a decrease of sugar nucleotides in nodes were found in bc1 and Bc6, in which the causative genes are involved in cellulose accumulation. CONCLUSION: These results suggest that metabolic analysis of the bc mutants could elucidate the functions of causative gene and improve the cell wall components for livestock feed or bioethanol production.

Pyramiding of multiple strong-culm genes originating from indica and tropical japonica to the temperate japonica rice.

Severe lodging has recurrently occurred at strong typhoon's hitting in recent climate change. The identification of quantitative trait loci and their responsible genes associated with a strong culm and their pyramiding are important for developing high-yielding varieties with a superior lodging resistance. To evaluate the effects of four strong-culm genes on lodging resistance, the temperate japonica near isogenic line (NIL) with the introgressed SCM1 or SCM2 locus of the indica variety, Habataki and the other NIL with the introgeressed SCM3 or SCM4 locus of the tropical japonica variety, Chugoku 117 were developed. Then, we developed the pyramiding lines with double,triple and quadruple combinations derived from step-by-step crosses among NIL-SCM1-NIL-SCM4. Quadruple pyramiding line (NIL-SCM1 + 2 + 3 + 4) showed the largest culm diameter and the highest culm strength among the combinations and increased spikelet number due to the pleiotropic effects of these genes. Pyramiding of strong culm genes resulted in much increased culm thickness, culm strength and spikelet number due to their additive effect. SCM1 mainly contributed to enhance their pyramiding effect. These results in this study suggest the importance of identifying the combinations of superior alleles of strong culm genes among natural variation and pyramiding these genes for improving high-yielding varieties with a superior lodging resistance.

Potential of rice landraces with strong culms as genetic resources for improving lodging resistance against super typhoons.

It is generally believed that rice landraces with long culms are susceptible to lodging, and have not been utilized for breeding to improve lodging resistance. However, little is known about the structural culm strength of landraces and their beneficial genetic loci. Therefore, in this study, genome-wide association studies (GWAS) were performed using a rice population panel including Japanese rice landraces to identify beneficial loci associated with strong culms. As a result, the landraces were found to have higher structural culm strength and greater diversity than the breeding varieties. Genetic loci associated with strong culms were identified, and it was demonstrated that haplotypes with positive effects of those loci were present in a high proportion of these landraces. These results indicated that the utilization of the strong culm-associated loci present in Japanese rice landraces may further improve the lodging resistance of modern breeding varieties that have relied on semi-dwarfism.

Maintaining higher leaf photosynthesis after heading stage could promote biomass accumulation in rice.

Leaf photosynthetic rate changes across the growing season as crop plants age. Most studies of leaf photosynthesis focus on a specific growth stage, leaving the question of which pattern of photosynthetic dynamics maximizes crop productivity unanswered. Here we obtained high-frequency data of canopy leaf CO2 assimilation rate (A) of two elite rice (Oryza sativa) cultivars and 76 inbred lines across the whole growing season. The integrated A value after heading was positively associated with crop growth rate (CGR) from heading to harvest, but that before heading was not. A curve-smoothing analysis of A after heading showed that accumulated A at > 80% of its maximum (A80) was positively correlated with CGR in analyses of all lines mixed and of lines grouped by genetic background, while the maximum A and accumulated A at ≤ 80% were less strongly correlated with CGR. We also found a genomic region (~ 12.2 Mb) that may enhance both A80 and aboveground biomass at harvest. We propose that maintaining a high A after heading, rather than having high maximum A, is a potential target for enhancing rice biomass accumulation.

Landraces of temperate japonica rice have superior alleles for improving culm strength associated with lodging resistance.

Lodging can reduce grain yield and quality in cereal crops including rice (Oryza sativa L.). To achieve both high biomass production and lodging resistance, the breeding of new cultivars with strong culms is a promising strategy. However, little is known about the diversity of culm strength in temperate japonica rice and underlying genetic factors. Here, we report a wide variation of culm strength among 135 temperate japonica cultivars, and some landraces having the strongest culms among these cultivars. The genome-wide association study (GWAS) identified 55 quantitative trait loci for culm strength and morphological traits, and revealed several candidate genes. The superior allele of candidate gene for culm thickness, OsRLCK191, was found in many landraces but had not inherited to the modern improved cultivars. Our results suggest that landraces of temperate japonica rice have unutilized superior alleles for contributing future improvements of culm strength and lodging resistance.

Next generation long-culm rice with superior lodging resistance and high grain yield, Monster Rice 1.

During late 1960s Green Revolution, researchers utilized semidwarf 1 (sd1) to improve the yield and lodging resistance in rice (Oryza sativa L.). However, sd1 has a negative effect to culm strength and biomass production. To increase yield dramatically in 21th century, development of next generation long-culm rice for non-lodging and high grain yield independent of sd1 has been needed. The present study developed Monster Rice 1, a long-culm and heavy-panicle type of rice line and compared it with Takanari, a high-yielding semidwarf rice variety about yield and lodging resistance associated traits. Brown rice yield and bending moment at breaking of the basal elongated internode were higher in Monster Rice 1 than those in Takanari due to a large number of spikelets per panicle and thicker culm. Furthermore, to identify QTLs with superior alleles for these traits, QTL and haplotype analyses were performed using F2 population and recombinant inbred lines derived from a cross between Monster Rice 1 and Takanari. The results from this study suggest that long-culm and heavy-panicle type of rice with a superior lodging resistance by culm strength can perform its high yield potential by using these identified QTLs contributing yield and lodging resistance.

Influence of rice-husk biochar and Bacillus pumilus strain TUAT-1 on yield, biomass production, and nutrient uptake in two forage rice genotypes.

Biochar is widely used as a soil amendment to increase crop yields. However, the impact of the interaction between the biochar and microbial inoculants (e.g., biofertilizer) on plant nutrient uptake and yield in forage rice is not fully understood. A greenhouse study was conducted to evaluate the synergistic effects of rice-husk biochar and Bacillus pumilus strain TUAT-1 biofertilizer application on growth, yield, and nutrient uptake in two forage rice genotypes; Fukuhibiki and the newly bred line, LTAT-29. Positive effects of biochar and biofertilizer, alone or in a combination, on growth traits, nutrient uptake, and yield components were dependent on the rice genotypes. Biochar and TUAT-1 biofertilizer influenced the overall growth of plants positively and increased straw and above-ground biomass in both genotypes. However, although biochar application significantly increased grain yield in LTAT-29, this was not the case in Fukuhibiki. Biochar and TUAT-1 biofertilizer, either alone or combined, significantly affected plant nutrient uptake but the effect largely depended on rice genotype. Results of this study indicate that biochar amendment and TUAT-1 biofertilizer can enhance forage rice productivity depending on genotypes, and therefore, there is a need to consider plant genetic composition when evaluating the potential for crop response to these soil amendments before application on a commercial scale.

Genetic architecture of leaf photosynthesis in rice revealed by different types of reciprocal mapping populations.

The improvement of leaf net photosynthetic rate (An) is a major challenge in enhancing crop productivity. However, the genetic control of An among natural genetic accessions is still poorly understood. The high-yielding indica cultivar Takanari has the highest An of all rice cultivars, 20-30% higher than that of the high-quality japonica cultivar Koshihikari. By using reciprocal backcross inbred lines and chromosome segment substitution lines derived from a cross between Takanari and Koshihikari, we identified three quantitative trait loci (QTLs) where the Takanari alleles enhanced An in plants with a Koshihikari genetic background and five QTLs where the Koshihikari alleles enhanced An in plants with a Takanari genetic background. Two QTLs were expressed in plants with both backgrounds (type I QTL). The expression of other QTLs depended strongly on genetic background (type II QTL). These beneficial alleles increased stomatal conductance, the initial slope of An versus intercellular CO2 concentration, or An at CO2 saturation. Pyramiding of these alleles consistently increased An. Some alleles positively affected biomass production and grain yield. These alleles associated with photosynthesis and yield can be a valuable tool in rice breeding programs via DNA marker-assisted selection.

Finding the superior allele of japonica-type for increasing stem lodging resistance in indica rice varieties using chromosome segment substitution lines.

BACKGROUND: In cereal crops, stem lodging can be classified into two types: stem-breaking type and stem-bending type. To improve stem-lodging resistance, the strong culm traits of superior lodging-resistant varieties must be characterized. The identification of quantitative trait loci (QTLs) and the corresponding genes associated with the parameters for bending moment at breaking (M) and flexural rigidity (FR) is expected to enable the efficient development of lodging-resistant varieties. A set of Chromosome Segment Substitution Lines (CSSLs) derived from the cross between Takanari and Koshihikari were used in this study to identify QTLs associated with lodging resistance. RESULTS: The indica variety Takanari possesses large M due to its large section modulus (SM) despite its small bending stress (BS), whereas Takanari also has large FR due to its large secondary moment of inertia (SMI) and Young's modulus (YM). The QTLs for BS were assigned to chromosomes 3, 5, 6, 8, 9, 10, 11, and 12. Koshihikari alleles increased BS in these QTLs. The YM was increased by substitution of the Koshihikari chromosomal segments on chromosomes 2, 10, and 11. Other QTLs mapped to chromosomes 7 and 12, such that the Koshihikari alleles contributed to the decrease of YM. QTLs for cellulose density were assigned to chromosomes 1, 3, and 5, which were replaced by substitutions of Koshihikari segments. The QTLs for hemicellulose, cellulose, and holocellulose densities identified on chromosome 5 overlapped with those for BS, indicating the positive effect of the Koshihikari segment on increasing BS. CONCLUSIONS: These results suggested that the QTLs for the densities of cell wall materials in japonica varieties contributed to increased BS and might be utilized for improving lodging resistance in indica varieties of rice.

QTLs maintaining grain fertility under salt stress detected by exome QTL-seq and interval mapping in barley.

Enhancing salt stress tolerance is a key strategy for increasing global food production. We previously found that long-term salinity stress significantly reduced grain fertility in the salt-sensitive barley (Hordeum vulgare) accession, 'OUC613', but not in the salt-tolerant accession, 'OUE812', resulting in large differences in grain yield. Here, we examined the underlying causes of the difference in grain fertility between these accessions under long-term treatment with 150 or 200 mM NaCl from the seedling stage to harvest and identified quantitative trait loci (QTLs) for maintaining grain fertility. In an artificial pollination experiment of the two accessions, grain fertility was significantly reduced only in OUC613 plants produced using pollen from plants grown under NaCl stress, suggesting that the low grain fertility of OUC613 was mainly due to reduced pollen fertility. Using QTL-seq combined with exome-capture sequencing and composite interval mapping of recombinant inbred lines derived from a cross between OUE812 and OUC613, we identified a QTL (qRP-2Hb) for grain fertility on chromosome 2H. The QTL region includes two genes encoding an F-box protein and a TIFY protein that are associated with male sterility, highlighting the importance of this region for maintaining grain fertility under salt stress.

Leaf photosynthetic rate and mesophyll cell anatomy changes during ontogenesis in backcrossed indica × japonica rice inbred lines.

The high-yielding indica rice variety, 'Takanari', has the high rate of leaf photosynthesis compared with the commercial japonica varieties. Among backcrossed inbred lines from a cross between 'Takanari' and a japonica variety, 'Koshihikari', two lines, BTK-a and BTK-b, showed approximately 20% higher photosynthetic rate than that of 'Takanari' for a flag leaf at full heading. This is a highest recorded rate of rice leaf photosynthesis. Here, the timing and cause of the increased leaf photosynthesis in the BTK lines were investigated by examining the photosynthesis and related parameters, as well as mesophyll cell anatomy during ontogenesis. Their photosynthetic rate was greater than that of 'Takanari' in the 13th leaf, as well as the flag leaf, but there were no differences in the 7th and 10th leaves. There were no consistent differences in the stomatal conductance, or the leaf nitrogen and Rubisco contents in the 13th and flag leaves. The total surface area of mesophyll cells per leaf area (TAmes) in the 13th and flag leaves increased significantly in the BTK lines due to the increased number and developed lobes of mesophyll cells compared with in 'Takanari'. The mesophyll conductance (g m) became greater in the BTK lines compared with 'Takanari' in the flag leaves but not in the 10th leaves. A close correlation was observed between TAmes and g m. We concluded that the increased mesophyll conductance through the development of mesophyll cells during the reproductive period is a probable cause of the greater photosynthetic rate in the BTK lines.

Fine Mapping of Carbon Assimilation Rate 8, a Quantitative Trait Locus for Flag Leaf Nitrogen Content, Stomatal Conductance and Photosynthesis in Rice.

Increasing the rate of leaf photosynthesis is one important approach for increasing grain yield in rice (Oryza sativa). Exploiting the natural variation in CO2 assimilation rate (A) between rice cultivars using quantitative genetics is one promising means to identify genes contributing to higher photosynthesis. In this study, we determined precise location of Carbon Assimilation Rate 8 (CAR8) by crossing a high-yielding indica cultivar with a Japanese commercial cultivar. Fine mapping suggested that CAR8 encodes a putative Heme Activator Protein 3 (OsHAP3) subunit of a CCAAT-box-binding transcription factor called OsHAP3H. Sequencing analysis revealed that the indica allele of CAR8 has a 1-bp deletion at 322 bp from the start codon, resulting in a truncated protein of 125 amino acids. In addition, CAR8 is identical to DTH8/Ghd8/LHD1, which was reported to control rice flowering date. The increase of A is largely due to an increase of RuBP regeneration rate via increased leaf nitrogen content, and partially explained by reduced stomatal limitation via increased stomatal conductance relative to A. This allele also increases hydraulic conductivity, which would promote higher stomatal conductance. This indicates that CAR8 affects multiple physiological aspects relating to photosynthesis. The detailed analysis of molecular functions of CAR8 would help to understand the association between photosynthesis and flowering and demonstrate specific genetic mechanisms that can be exploited to improve photosynthesis in rice and potentially other crops.

Precise estimation of genomic regions controlling lodging resistance using a set of reciprocal chromosome segment substitution lines in rice.

Severe lodging has occurred in many improved rice varieties after the recent strong typhoons in East and Southeast Asian countries. The indica variety Takanari possesses strong culm characteristics due to its large section modulus, which indicates culm thickness, whereas the japonica variety Koshihikari is subject to substantial bending stress due to its thick cortical fibre tissue. To detect quantitative trait loci (QTLs) for lodging resistance and to eliminate the effects of genetic background, we used reciprocal chromosome segment substitution lines (CSSLs) derived from a cross between Koshihikari and Takanari. The oppositional effects of QTLs for section modulus were confirmed in both genetic backgrounds on chromosomes 1, 5 and 6, suggesting that these QTLs are not affected by the genetic background and are controlled independently by a single factor. The candidate region of a QTL for section modulus included SD1. The section modulus of NIL-sd1 was lower than that of Koshihikari, whereas the section modulus of NIL-SD1 was higher than that of Takanari. This result indicated that those regions regulate the culm thickness. The reciprocal effects of the QTLs for cortical fibre tissue thickness were confirmed in both genetic backgrounds on chromosome 9 using CSSLs.

Changes in hydraulic conductance cause the difference in growth response to short-term salt stress between salt-tolerant and -sensitive black gram (Vigna mungo) varieties.

Black gram (Vigna mungo) is an important crop in Asia, However, most black gram varieties are salt-sensitive. The causes of varietal differences in salt-induced growth reduction between two black gram varieties, 'U-Taung-2' (salt-tolerant; BT) and 'Mut Pe Khaing To' (salt-sensitive; BS), were examined the potential for the first step toward the genetic improvement of salt tolerance. Seedlings grown in vermiculite irrigated with full-strength Hoagland solution were treated with 0mM NaCl (control) or 225 mM NaCl for up to 10 days. In the 225 mM NaCl treatment, plant growth rate, net assimilation rate, mean leaf area, leaf water potential, and leaf photosynthesis were reduced more in BS than in BT plants. Leaf water potential was closely related to leaf photosynthesis, net assimilation rate, and increase in leaf area. In response to salinity stress, hydraulic conductance of the root, stem, and petiole decreased more strongly in BS than in BT plants. The reduction in stem and petiole hydraulic conductance was caused by cavitation, whereas the reduction in root hydraulic conductance in BS plants was caused by a reduction in root surface area and hydraulic conductivity. We conclude that the different reduction in hydraulic conductance is a cause of the differences in the growth response between the two black gram varieties under short-term salt stress.

Detection of QTL for exudation rate at ripening stage in rice and its contribution to hydraulic conductance.

Dry matter production of crops is determined by how much light they intercept and how efficiently they use it for carbon fixation; i.e., photosynthesis. The high-yielding rice cultivar, Akenohoshi, maintains a high photosynthetic rate in the middle of the day owing to its high hydraulic conductance in comparison with the elite commercial rice cultivar, Koshihikari. We developed 94 recombinant inbred lines derived from Akenohoshi and Koshihikari and measured their exudation rate to calculate hydraulic conductance to osmotic water transport in a paddy field. A quantitative trait locus (QTL) for exudation rate was detected on the long arm of chromosome 2 at the heading and ripening stages. We developed chromosome segment substitution lines which carried Akenohoshi segments in the Koshihikari genetic background, and measured hydraulic conductance to both osmotic and passive water transport. The QTL was confirmed to be located within a region of about 4.2Mbp on the distal end of long arm of chromosome 2. The Akenohoshi allele increased root surface area and hydraulic conductance, but didn't increase hydraulic conductivity of a plant.

Natural Variation in the Flag Leaf Morphology of Rice Due to a Mutation of the NARROW LEAF 1 Gene in Oryza sativa L.

We investigated the natural variations in the flag leaf morphology of rice. We conducted a principal component analysis based on nine flag leaf morphology traits using 103 accessions from the National Institute of Agrobiological Sciences Core Collection. The first component explained 39% of total variance, and the variable with highest loading was the width of the flag leaf (WFL). A genome-wide association analysis of 102 diverse Japanese accessions revealed that marker RM6992 on chromosome 4 was highly associated with WFL. In analyses of progenies derived from a cross between Takanari and Akenohoshi, the most significant quantitative trait locus (QTL) for WFL was in a 10.3-kb region containing the NARROW LEAF 1 (NAL1) gene, located 0.4 Mb downstream of RM6992. Analyses of chromosomal segment substitution lines indicated that a mutation (G1509A single-nucleotide mutation, causing an R233H amino acid substitution in NAL1) was present at the QTL. This explained 13 and 20% of total variability in WFL and the distance between small vascular bundles, respectively. The mutation apparently occurred during rice domestication and spread into japonica, tropical japonica, and indica subgroups. Notably, one accession, Phulba, had a NAL1 allele encoding only the N-terminal, or one-fourth, of the wild-type peptide. Given that the Phulba allele and the histidine-type allele showed essentially the same phenotype, the histidine-type allele was regarded as malfunctional. The phenotypes of transgenic plants varied depending on the ratio of histidine-type alleles to arginine-type alleles, raising the possibility that H(233)-type products function differently from and compete with R(233)-type products.

Isolation of a novel lodging resistance QTL gene involved in strigolactone signaling and its pyramiding with a QTL gene involved in another mechanism.

Lodging has been a major roadblock to attaining increased crop productivity. In an attempt to understand the mechanism for culm strength in rice, we isolated an effective quantitative trait locus (QTL), STRONG CULM3 (SCM3), the causal gene of which is identical to rice TEOSINTE BRANCHED1 (OsTB1), a gene previously reported to positively control strigolactone (SL) signaling. A near-isogenic line (NIL) carrying SCM3 showed enhanced culm strength and increased spikelet number despite the expected decrease in tiller number, indicating that SL also has a positive role in enhancing culm strength and spikelet number. We produced a pyramiding line carrying SCM3 and SCM2, another QTL encoding APO1 involved in panicle development. The NIL-SCM2+SCM3 showed a much stronger culm than NIL-SCM2 and NIL-SCM3 and an increased spikelet number caused by the additive effect of these QTLs. We discuss the importance of utilizing suitable alleles of these STRONG CULM QTLs without inducing detrimental traits for breeding.