QTL scanning for rice yield using a whole genome SNP array

Powerful QTL mapping and favorable allele mining in an all-in-one population: a case study of heading date

The multiparent advanced generation intercross (MAGIC) population is characterized with great potentials in power and resolution of quantitative trait locus (QTL) mapping, but single nucleotide polymorphism (SNP)-based GWAS does not fully reach its potential. In this study, a MAGIC population of 1021 lines was developed from four Xian and four Geng varieties from five subgroups of rice. A total of 44 000 genes showed functional polymorphisms among eight parents, including frameshift variations or premature stop codon variations, which provides the potential to map almost all genes of the MAGIC population. Principal component analysis results showed that the MAGIC population had a weak population structure. A high-density bin map of 24 414 bins was constructed. Segregation distortion occurred in the regions possessing the genes underlying genetic incompatibility and gamete development. SNP-based association analysis and bin-based linkage analysis identified 25 significant loci and 47 QTLs for heading date, including 14 known heading date genes. The mapping resolution of genes is dependent on genetic effects with offset distances of <55 kb for major effect genes and <123 kb for moderate effect genes. Four causal variants and noncoding structure variants were identified to be associated with heading date. Three to four types of alleles with strong, intermediate, weak, and no genetic effects were identified from eight parents, providing flexibility for the improvement of rice heading date. In most cases, japonica rice carries weak alleles, and indica rice carries strong alleles and nonfunctional alleles. These results confirm that the MAGIC population provides the exceptional opportunity to detect QTLs, and its use is encouraged for mapping genes and mining favorable alleles for breeding.

Fine mapping of the panicle length QTL qPL5 in rice

Panicle length is a crucial trait tightly associated with spikelets per panicle and grain yield in rice. To dissect the genetic basis of panicle length, a population of 161 recombinant inbred lines (RILs) was developed from the cross between an aus variety Chuan 7 (C7) and a tropical Geng variety Haoboka (HBK). C7 has a panicle length of 30 cm, 7 cm longer than that of HBK, and the panicle length was normally distributed in the RIL population. A total of six quantitative trait loci (QTLs) for panicle length were identified, and single QTLs explained the phenotypic variance from 4.9 to 18.1%. Among them, three QTLs were mapped to the regions harbored sd1, DLT, and Ehd1, respectively. To validate the genetic effect of a minor QTL qPL5, a near-isogenic F2 (NIF2) population segregated at qPL5 was developed. Interestingly, panicle length displayed bimodal distribution, and heading date also exhibited significant variation in the NIF2 population. qPL5 accounted for 66.5% of the panicle length variance. The C7 allele at qPL5 increased panicle length by 2.4 cm and promoted heading date by 5 days. Finally, qPL5 was narrowed down to an 80-kb region flanked by markers M2197 and M2205 using a large NIF2 population of 7600 plants. LOC_Os05g37540, encoding a phytochrome signal protein whose homolog in Arabidopsis enlarges panicle length, is regarded as the candidate gene because a single-nucleotide mutation (C1099T) caused a premature stop codon in HBK. The characterization of qPL5 with enlarging panicle length but promoting heading date makes its great value in breeding early mature varieties without yield penalty in rice.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01443-2.

Yield components affected by rice black-streaked dwarf virus disease in rice cultivars with different resistance levels

Introduction

Rice black-streaked dwarf virus disease (RBSDVD) is one of the most destructive rice viral diseases, leading to severe yield losses in rice production. However, little is known about the yield-related components associated with the disease and no resistance cultivars have been successfully used in rice breeding.

Methods

Seven rice cultivars were analyzed in this study, including six commercial rice varieties and a new line Zhongjian No. 201 (ZJ201) containing the resistance gene OsAP47. Resistance levels of these cultivars were evaluated by artificial inoculation and yield components were collected, including panicle length (PL), spikelets per panicle (SPP), ripened grains per panicle (RGPP), as well as panicles per square meter (PPSM) and 1000-grain weight (TGW). Seed setting rate (SSR) were calculated with the data of SPP and RGPP.

Results and discussion

The results showed that ZJ201 displayed the highest resistance level and most of the commercial rice cultivars exhibited susceptible to RBSDVD. Yields of all the rice cultivars were significantly declined except ZJ201 and yield losses produced by RBSDVD were mainly due to the reduction of PL, SPP, RGPP, and TGW, suggesting that developments of these traits are associated with RBSDV infection. Resistant rice cultivar could reduce yield losses by maintaining normal development of these traits. Significant correlations were identified between resistance levels and the yield components except SSR and PPSM. The results provided useful clues for understanding the mechanisms of RBSDV invasion and its effect on rice production. ZJ201 was demonstrated as a resistance material that could be used in rice breeding.

Genetic Mapping to Detect Stringent QTLs Using 1k-RiCA SNP Genotyping Platform from the New Landrace Associated with Salt Tolerance at the Seedling Stage in Rice

Rice is the world’s most important food crop, providing the daily calorie intake for more than half of the world’s population. Rice breeding has always been preoccupied with maximizing yield potential. However, numerous abiotic factors, such as salt, cold, drought, and heat, significantly reduce rice productivity. Salinity, one of the major abiotic stresses, reduces rice yield worldwide. This study was conducted to determine new quantitative trait loci (QTLs) that regulate salt tolerance in rice seedlings. One F_2:3 mapping population was derived from a cross between BRRI dhan49 (a popular but sensitive rainfed rice variety) and Akundi (a salt-tolerant rice landrace in Bangladesh used as a donor parent). The 1k-Rice Custom Amplicon (1k-RiCA) single-nucleotide polymorphism (SNP) markers were used to genotype this mapping population. After removing segregation distortion and monomorphic markers, 884 SNPs generated a 1526.8 cM-long genetic linkage map with a mean marker density of 1.7 cM for the 12 linkage groups. By exploiting QGene and ICIM-ADD, a sum of 15 QTLs for nine traits was identified in salt stress on seven chromosomes. Four important genomic loci were identified (qSES1, qSL1, qSUR1 and qRL1) on chromosome 1. Out of these 15 QTLs, 14 QTLs are unique, as no other study has mapped in the same chromosomal location. We also detected 15 putative candidate genes and their functions. The ICIM-EPI approach identified 43 significant pairwise epistasis interactions between regions associated with and unassociated with QTLs. Apart from more well-known donors, Akundi serves as an important new donor source for global salt tolerance breeding initiatives, including Bangladesh. The introgression of the novel QTLs identified in this study will accelerate the development of new salt-tolerant varieties that are highly resistant to salt stress using marker-enabled breeding.

Genetic Analysis of a Collection of Rice Germplasm (Oryza sativa L.) through High-Density SNP Array Provides Useful Information for Further Breeding Practices

Traditional breeding strategies mainly focus on the evaluation of trait performance, but pay less attention to the changing genetic background. A comprehensive understanding of the genetic diversity of germplasms is crucial for the deliberate improvement of specific traits. A collection of 154 highland rice varieties were collected as the initial genetic resource in our breeding program to improve the pathogen resistance and eating and cooking qualities. These varieties were analyzed using a whole-genome SNP array and were clustered into three groups. Further analysis revealed that the favorable alleles of pathogen resistance genes are mostly absent in our collected varieties. However, it showed that most varieties possess favorable alleles of Waxy (Wx) and ALKALI DEGENERATION (ALK), which are able to enhance the eating and cooking qualities. Moreover, only about one fifth of all varieties harbors favorable the allele of fragrance gene Betainealdehyde dehydrogenase (BADH2). Together, these results give an overall view of the genetic constitution of the target traits, which provide useful information for future genetic improvement in breeding practices.

Five plants per RIL for phenotyping traits of high or moderate heritability ensure the power of QTL mapping in a rice MAGIC population

Currently, the power of QTL mapping is mainly dependent on the quality of phenotypic data in a given population, regardless of the statistical method, as the quality of genotypic data is easily guaranteed in the laboratory. Increasing the sample size per line used for phenotyping is a good way to improve the quality of phenotypic data. However, accommodating a large-scale mapping population takes a large area of rice field, which frequently results in high costs and extra environmental noises. To acquire a reasonable small sample size without a penalty in mapping power, we conducted three experiments with a 4-way MAGIC population and measured phenotypes of 5, 10, and 20 plants per RIL. Three traits including heading date, plant height, and tillers per plant were focused. With SNP- and bin-based QTL mapping, 3 major and 3 minor QTLs for heading date with high heritability and 2 major QTLs for plant height with moderate heritability were commonly detected across the three experiments, but no QTL for tillers per plant with low heritability were commonly identified. In addition, bin-based QTL mapping was more powerful than SNP-based mapping and able to rank the genetic effects of parental alleles. Thus, 5 plants per RIL for phenotyping ensure the power of QTL mapping for traits of high or moderate heritability, and bin-based QTL mapping is recommended for multiparent populations.

Identification and verification of grain shape QTLs by SNP array in rice

Grain shape strongly influences the economic value and grain yield of rice. Thus, identifying quantitative trait loci (QTLs) for grain shape has been a longstanding goal in rice genetic research and breeding programs. Single nucleotide polymorphism (SNP) markers are ubiquitous in the rice genome and are more abundant and evenly distributed on the 12 rice chromosomes than traditional markers. An F₂ population was genotyped using the RICE6K SNP array to elucidate the mechanisms governing grain shape. Thirty-five QTLs for grain shape were detected on 11 of 12 chromosomes over 2 years. The major QTL cluster qGS7 was detected in both years and displayed strong genetic effects on grain length and width, showing consistency with GL7/GW7. Some minor QTLs were also detected, and the effects of four QTLs on seed size were then validated using BC₁F₆ populations with residual heterozygous lines in each QTL region. Our findings provide insights into the molecular basis of grain shape as well as additional resources and approaches for producing hybrid high-yield rice varieties.

Dominant complementary interaction between OsC1 and two tightly linked genes, Rb1 and Rb2, controls the purple leaf sheath in rice

Two tightly linked genes for rice purple leaf sheath were identified via map-based cloning. Further analysis indicated that these two genes together with OsC1 co-regulating the purple leaf sheath. The purple color of the leaf sheath in rice is dependent on the accumulation of anthocyanins such as cyanidin 3-O-glucoside (C3G) and peonidin 3-O-glucoside (P3G). Although many genes related to leaf sheath color have been mapped, the genetic basis for leaf sheath color is not yet clear. Here, PSH1 (purple leaf sheath 1) was mapped to chromosome 1 using an F₂ and a RIL population. Map-based cloning and transformation assays further divided PSH1 as two tightly linked bHLH genes, Rb1 and Rb2. Ectopic expression of these two genes resulted in substantial accumulation of C3G and P3G in the leaf blade, leaf sheath and pericarp. Single gene mutants displayed a faded purple leaf sheath or green leaf sheath in the top half of the leaf sheath, but double mutants displayed a green leaf sheath, indicating that both genes have dosage effects on anthocyanin synthesis. However, overexpression of Rb1 and Rb2 sharply decreased grain filling. A segregation ratio of green to purple was 15:1 observed in the F₂ population from parents Minghui 63 and Xizang 2, which both had green leaf sheaths; these results demonstrate that dominant complementary interaction between OsC1 and Rb (Rb1 and Rb2) controls the purple leaf sheath. These findings systematically uncovered the genetic basis of leaf sheath color and provided alternative genes for breeding anthocyanin-rich rice.

Bin-based genome-wide association analyses improve power and resolution in QTL mapping and identify favorable alleles from multiple parents in a four-way MAGIC rice population

A whole genome bin map was developed for a MAGIC population. Association studies for heading date at bin level exhibited powerful QTL mapping and identified favorable alleles. The presumed advantages of multiparent advanced generation intercross (MAGIC) population in quantitative trait locus (QTL) mapping were not fully utilized in the previous studies in which genome-wide association studies (GWAS) were conducted at only single nucleotide polymorphism level. In this study, we genotyped a rice four-way MAGIC population of 247 F₇ lines and their parents by sequencing. A total of 5934 bins with an average length of 65 kb were constructed and covered 97% of the genome. The MAGIC population showed low population structure and balanced parental contributions. A bin-based GWAS for heading date identified 4 QTLs in three environments. Three major QTLs were mapped exactly to the bins where the major heading date genes DTH3, Ghd7.1 and Ghd8 were located. Multiple comparisons showed that different parental alleles had varied genetic effects. Like DTH3, the alleles of the Guichao 2/YJSM, IR34 and Cypress had larger, intermediate and no effects, respectively. Based on comparative sequencing of 8 known heading date genes undetected in this MAGIC population, only Ghd7 exhibited diverse function among parents. The failure in Ghd7 mapping was well explained by its interaction with Hd1 because Ghd7 had no effects on heading date when combined with the nonfunctional hd1 carried by all four parents. Overall, bin-based GWAS have more mapping power and higher resolution with a MAGIC population and provide favorable alleles to breeders. The use of more diversified parents is encouraged to develop a MAGIC population for detecting more QTLs for important agronomical traits.

QTL fine-mapping of soybean (Glycine max L.) leaf type associated traits in two RILs populations

BACKGROUND

The different leaf type associated traits of soybean (Glycine max L.) including leaf area, leaf length, leaf width, leaf shape and petiole length are considered to be associated with seed yield. In order to identify quantitative trait loci (QTLs) affecting leaf type traits, two advanced recombinant inbred line (RIL, ZH, Zhonghuang 24 × Huaxia 3; GB, Guizao 1 × Brazil 13) populations were introduced to score phenotypic values in plants across nine different environments (years, seasons, locations and soybean growth stages). Two restriction site-associated DNA sequencing (RAD-seq) based high-density genetic linkage maps with an average distance of 1.00 centimorgan (cM) between adjacent bin markers were utilized for QTL fine mapping.

RESULTS

Correlation analysis showed that most of the traits were correlated with each other and regulated both by hereditary and environmental factors. A total of 190 QTLs were identified for leaf type associated traits in the two populations, of which 14 loci were found to be environmentally stable. Moreover, these detected QTLs were categorized into 34 QTL hotspots, and four important QTL hotspots with phenotypic variance ranging from 3.89-23.13% were highlighted. Furthermore, Glyma04g05840, Glyma19g37820, Glyma14g07140 and Glyma19g39340 were predicted in the intervals of the stable loci and important QTL hotspots for leaf type traits by adopting Gene Ontology (GO) enrichment analysis.

CONCLUSIONS

Our findings of the QTLs and the putative genes will be beneficial to gain new insights into the genetic basis for soybean leaf type traits and may further accelerate the breeding process for reasonable leaf type soybean.

Intragenic recombination between two non-functional semi-dwarf 1 alleles produced a functional SD1 allele in a tall recombinant inbred line in rice

Intragenic recombination is one of the most important sources of genetic variability. In our previous study, RI92 a tall line (160 cm of plant height) was observed in the cross progeny between two semi-dwarf indica cultivars Zhenshan 97 and Minghui 63. Genome-wide genotyping and sequencing indicated that the genome constitution of RI92 was completely from both parents. Bulk segregant analysis in a BC₃F₂ population revealed that “green revolution gene” semi-dwarf 1 (sd1) was most likely the gene controlling the tall plant height in RI92. Sequencing analysis of SD1 revealed that an intragenic recombination occurred between two parental non-functional sd1 alleles and generated a functional SD1 in RI92. Four-fold high recombination rate in SD1 located bins to the genome-wide average was observed in two RIL populations, indicating recombination hotspot in the SD1 region. Intragenic recombination creates new alleles in the progeny distinct from parental alleles and diversifies natural variation.

Construction of high-density genetic map and QTL mapping of yield-related and two quality traits in soybean RILs population by RAD-sequencing

BACKGROUND

One of the overarching goals of soybean breeding is to develop lines that combine increased yield with improved quality characteristics. High-density-marker QTL mapping can serve as an effective strategy to identify novel genomic information to facilitate crop improvement. In this study, we genotyped a recombinant inbred line (RIL) population (Zhonghuang 24 × Huaxia 3) using a restriction-site associated DNA sequencing (RAD-seq) approach. A high-density soybean genetic map was constructed and used to identify several QTLs that were shown to influence six yield-related and two quality traits.

RESULTS

A total of 47,472 single-nucleotide polymorphisms (SNPs) were detected for the RILs that were integrated into 2639 recombination bin units, with an average distance of 1.00 cM between adjacent markers. Forty seven QTLs for yield-related traits and 13 QTLs for grain quality traits were found to be distributed on 16 chromosomes in the 2 year studies. Among them, 18 QTLs were stable, and were identified in both analyses. Twenty six QTLs were identified for the first time, with a single QTL (qNN19a) in a 56 kb region explaining 32.56% of phenotypic variation, and an additional 10 of these were novel, stable QTLs. Moreover, 8 QTL hotpots on four different chromosomes were identified for the correlated traits.

CONCLUSIONS

With RAD-sequencing, some novel QTLs and important QTL clusters for both yield-related and quality traits were identified based on a new, high-density bin linkage map. Three predicted genes were selected as candidates that likely have a direct or indirect influence on both yield and quality in soybean. Our findings will be helpful for understanding common genetic control mechanisms of co-localized traits and to select cultivars for further analysis to predictably modulate soybean yield and quality simultaneously.

QTLs for heading date and plant height under multiple environments in rice

Both heading date and plant height are important traits related to grain yield in rice. In this study, a recombinant inbred lines (RILs) population was used to map quantitative trait loci (QTLs) for both traits under 3 long-day (LD) environments and 1 short-day (SD) environment. A total of eight QTLs for heading date and three QTLs for plant height were detected by composite interval mapping under LD conditions. Additional one QTL for heading date and three QTLs for plant height were identified by Two-QTL model under LD conditions. Among them, major QTLs qHd7.1, qHd7.2 and qHd8 for heading date, and qPh1 and qPh7.1 for plant height were commonly detected. qHd7.1 and qHd7.2 were mapped to small regions of less than 1 cM. Genome position comparison of previously cloned genes with QTLs detected in this study revealed that qHd5 and qPh3.1 were two novel QTLs. The alleles of these QTLs increasing trait values were dispersed in both parents, which well explained the transgressive segregation observed in this population. In addition, the interaction between qHd7.1 and qHd8 was detected under all LD conditions. Multiple-QTL model analysis revealed that all QTLs and their interactions explained over 80% of heading date variation and 50% of plant height variation. Two heading date QTLs were detected under SD condition. Of them, qHd10 were commonly identified under LD condition. The difference in QTL detection between LD and SD conditions indicated most heading date QTLs are sensitive to photoperiod. These findings will benefit breeding design for heading date and plant height in rice.

SNP-based discovery of salinity-tolerant QTLs in a bi-parental population of rice (Oryza sativa)

Breeding for salt tolerance is the most promising approach to enhance the productivity of saline prone areas. However, polygenic inheritance of salt tolerance in rice acts as a bottleneck in conventional breeding for salt tolerance. Hence, we set our goals to construct a single nucleotide polymorphism (SNP)-based molecular map employing high-throughput SNP marker technology and to investigate salinity tolerant QTLs with closest flanking markers using an elite rice background. Seedling stage salinity responses were assessed in a population of 281 recombinant inbred lines (RILs) derived from the cross between At354 (salt tolerant) and Bg352 (salt susceptible), by 11 morpho-physiological indices under a hydroponic system. Selected extreme 94 RILs were genotyped using Illumina Infinium rice 6K SNP array and densely saturated molecular map spanning 1460.81 cM of the rice genome with an average interval of 1.29 cM between marker loci was constructed using 1135 polymorphic SNP markers. The results revealed 83 significant QTLs for 11 salt responsive traits explaining 12.5-46.7 % of phenotypic variation in respective traits. Of them, 72 QTLs responsible for 10 traits were co-localized together forming 14 QTL hotspots at 14 different genomic regions. The all QTL hotspots were flanked less than 1 Mb intervals and therefore the SNP loci associated with these QTL hotspots would be important in candidate gene discovery for salt tolerance.

Identification of stable QTLs for vegetative and reproductive traits in the microvine (Vitis vinifera L.) using the 18 K Infinium chip

BACKGROUND

The increasing temperature associated with climate change impacts grapevine phenology and development with critical effects on grape yield and composition. Plant breeding has the potential to deliver new cultivars with stable yield and quality under warmer climate conditions, but this requires the identification of stable genetic determinants. This study tested the potentialities of the microvine to boost genetics in grapevine. A mapping population of 129 microvines derived from Picovine x Ugni Blanc flb, was genotyped with the Illumina® 18 K SNP (Single Nucleotide Polymorphism) chip. Forty-three vegetative and reproductive traits were phenotyped outdoors over four cropping cycles, and a subset of 22 traits over two cropping cycles in growth rooms with two contrasted temperatures, in order to map stable QTLs (Quantitative Trait Loci).

RESULTS

Ten stable QTLs for berry development and quality or leaf area were identified on the parental maps. A new major QTL explaining up to 44 % of total variance of berry weight was identified on chromosome 7 in Ugni Blanc flb, and co-localized with QTLs for seed number (up to 76 % total variance), major berry acids at green lag phase (up to 35 %), and other yield components (up to 25 %). In addition, a minor QTL for leaf area was found on chromosome 4 of the same parent. In contrast, only minor QTLs for berry acidity and leaf area could be found as moderately stable in Picovine. None of the transporters recently identified as mutated in low acidity apples or Cucurbits were included in the several hundreds of candidate genes underlying the above berry QTLs, which could be reduced to a few dozen candidate genes when a priori pertinent biological functions and organ specific expression were considered.

CONCLUSIONS

This study combining the use of microvine and a high throughput genotyping technology was innovative for grapevine genetics. It allowed the identification of 10 stable QTLs, including the first berry acidity QTLs reported so far in a Vitis vinifera intra-specific cross. Robustness of a set of QTLs was assessed with respect to temperature variation.

A whole-genome SNP array (RICE6K) for genomic breeding in rice

The advances in genotyping technology provide an opportunity to use genomic tools in crop breeding. As compared to field selections performed in conventional breeding programmes, genomics-based genotype screen can potentially reduce number of breeding cycles and more precisely integrate target genes for particular traits into an ideal genetic background. We developed a whole-genome single nucleotide polymorphism (SNP) array, RICE6K, based on Infinium technology, using representative SNPs selected from more than four million SNPs identified from resequencing data of more than 500 rice landraces. RICE6K contains 5102 SNP and insertion-deletion (InDel) markers, about 4500 of which were of high quality in the tested rice lines producing highly repeatable results. Forty-five functional markers that are located inside 28 characterized genes of important traits can be detected using RICE6K. The SNP markers are evenly distributed on the 12 chromosomes of rice with the average density of 12 SNPs per 1 Mb and can provide information for polymorphisms between indica and japonica subspecies as well as varieties within indica and japonica groups. Application tests of RICE6K showed that the array is suitable for rice germplasm fingerprinting, genotyping bulked segregating pools, seed authenticity check and genetic background selection. These results suggest that RICE6K provides an efficient and reliable genotyping tool for rice genomic breeding.