Genome-wide Association Analysis Tracks Bacterial Leaf Blight Resistance Loci In Rice Diverse Germplasm

Development of Novel KASP Markers for Improved Germination in Deep-Sown Direct Seeded Rice

BACKGROUND

The lack of stable-high yielding and direct-seeded adapted varieties with better germination ability from deeper soil depth and availability of molecular markers are major limitation in achieving the maximum yield potential of rice under water and resource limited conditions. Development of high-throughput and trait-linked markers are of great interest in genomics-assisted breeding. The aim of present study was to develop and validate novel KASP (Kompetitive Allele-Specific PCR) markers associated with traits improving germination and seedling vigor of deep sown direct seeded rice (DSR).

RESULTS

Out of 58 designed KASP assays, four KASP assays did not show any polymorphism in any of the eleven genetic backgrounds considered in the present study. The 54 polymorphic KASP assays were then validated for their robustness and reliability on the F1s plants developed from eight different crosses considered in the present study. The third next validation was carried out on 256 F3:F4 and 713 BC3F2:3 progenies. Finally, the reliability of the KASP assays was accessed on a set of random 50 samples from F3:F4 and 80-100 samples from BC3F2:3 progenies using the 10 random markers. From the 54 polymorphic KASP, based on the false positive rate, false negative rate, KASP utility in different genetic backgrounds and significant differences in the phenotypic values of the positive (desirable) and negative (undesirable) traits, a total of 12 KASP assays have been selected. These 12 KASP include 5 KASP on chromosome 3, 1 on chromosome 4, 3 on chromosome 7 and 3 on chromosome 8. The two SNPs lying in the exon regions of LOC_Os04g34290 and LOC_Os08g32100 led to non-synonymous mutations indicating a possible deleterious effect of the SNP variants on the protein structure.

CONCLUSION

The present research work will provide trait-linked KASP assays, improved breeding material possessing favourable alleles and breeding material in form of expected pre-direct-seeded adapted rice varieties. The marker can be utilized in introgression program during pyramiding of valuable QTLs/genes providing adaptation to rice under DSR. The functional studies of the genes LOC_Os04g34290 and LOC_Os08g32100 possessing two validated SNPs may provide valuable information about these genes.

Genome-Wide Association Study Identifies Resistance Loci for Bacterial Blight in a Collection of Asian Temperate Japonica Rice Germplasm

Growing resistant rice cultivars is the most effective strategy to control bacterial blight (BB), a devastating disease caused by Xanthomonas oryzae pv. oryzae (Xoo). Screening resistant germplasm and identifying resistance (R) genes are prerequisites for breeding resistant rice cultivars. We conducted a genome-wide association study (GWAS) to detect quantitative trait loci (QTL) associated with BB resistance using 359 East Asian temperate Japonica accessions inoculated with two Chinese Xoo strains (KS6-6 and GV) and one Philippine Xoo strain (PXO99^A). Based on the 55K SNPs Array dataset of the 359 Japonica accessions, eight QTL were identified on rice chromosomes 1, 2, 4, 10, and 11. Four of the QTL coincided with previously reported QTL, and four were novel loci. Six R genes were localized in the qBBV-11.1, qBBV-11.2, and qBBV-11.3 loci on chromosome 11 in this Japonica collection. Haplotype analysis revealed candidate genes associated with BB resistance in each QTL. Notably, LOC_Os11g47290 in qBBV-11.3, encoding a leucine-rich repeat receptor-like kinase, was a candidate gene associated with resistance to the virulent strain GV. Knockout mutants of Nipponbare with the susceptible haplotype of LOC_Os11g47290 exhibited significantly improved BB resistance. These results will be useful for cloning BB resistance genes and breeding resistant rice cultivars.

Evaluation of nine statistics to identify QTLs in bulk segregant analysis using next generation sequencing approaches

Background

Bulk segregant analysis (BSA) combined with next generation sequencing is a powerful tool to identify quantitative trait loci (QTL). The impact of the size of the study population and the percentage of extreme genotypes analysed have already been assessed. But a good comparison of statistical approaches designed to identify QTL regions using next generation sequencing (NGS) technologies for BSA is still lacking.

Results

We developed an R code to simulate QTLs in bulks of F2 contrasted lines. We simulated a range of recombination rates based on estimations using different crop species. The simulations were used to benchmark the ability of statistical methods identify the exact location of true QTLs. A single QTL led to a shift in allele frequency across a large fraction of the chromosome for plant species with low recombination rate. The smoothed version of all statistics performed best notably the smoothed Euclidean distance-based statistics was always found to be more accurate in identifying the location of QTLs. We propose a simulation approach to build confidence interval statistics for the detection of QTLs.

Conclusion

We highlight the statistical methods best suited for BSA studies using NGS technologies in crops even when recombination rate is low. We also provide simulation codes to build confidence intervals and to assess the impact of recombination for application to other studies. This computational study will help select NGS-based BSA statistics that are useful to the broad scientific community.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08718-y.

Long Intergenic Noncoding RNA (lincRNA) Discovery from Non-Strand-Specific RNA-Seq Data

Long noncoding RNAs (lncRNAs) are transcripts over 200 base pairs in length without discernible protein coding potential. Long intergenic noncoding RNAs (lincRNAs) constitute a subset of lncRNAs, which do not overlap protein coding genes. Here we describe a detailed pipeline for lincRNA discovery from publicly available non-stranded RNA-Seq datasets. The pipeline presented can be applied to any plant species for which RNA-Seq data and a reference genome sequence are available.

Two decades of association mapping: Insights on disease resistance in major crops

Climate change across the globe has an impact on the occurrence, prevalence, and severity of plant diseases. About 30% of yield losses in major crops are due to plant diseases; emerging diseases are likely to worsen the sustainable production in the coming years. Plant diseases have led to increased hunger and mass migration of human populations in the past, thus a serious threat to global food security. Equipping the modern varieties/hybrids with enhanced genetic resistance is the most economic, sustainable and environmentally friendly solution. Plant geneticists have done tremendous work in identifying stable resistance in primary genepools and many times other than primary genepools to breed resistant varieties in different major crops. Over the last two decades, the availability of crop and pathogen genomes due to advances in next generation sequencing technologies improved our understanding of trait genetics using different approaches. Genome-wide association studies have been effectively used to identify candidate genes and map loci associated with different diseases in crop plants. In this review, we highlight successful examples for the discovery of resistance genes to many important diseases. In addition, major developments in association studies, statistical models and bioinformatic tools that improve the power, resolution and the efficiency of identifying marker-trait associations. Overall this review provides comprehensive insights into the two decades of advances in GWAS studies and discusses the challenges and opportunities this research area provides for breeding resistant varieties.

Inducible Enrichment of Osa-miR1432 Confers Rice Bacterial Blight Resistance through Suppressing OsCaML2

MicroRNAs (miRNAs) handle immune response to pathogens by adjusting the function of target genes in plants. However, the experimentally documented miRNA/target modules implicated in the interplay between rice and Xanthomonas oryzae pv. oryzae (Xoo) are still in the early stages. Herein, the expression of osa-miR1432 was induced in resistant genotype IRBB5, but not susceptible genotype IR24, under Xoo strain PXO86 attack. Overexpressed osa-miR1432 heightened rice disease resistance to Xoo, indicated by enhancive enrichment of defense marker genes, raised reactive oxygen species (ROS) levels, repressed bacterial growth and shortened leaf lesion length, whilst the disruptive accumulation of osa-miR1432 accelerated rice susceptibility to Xoo infection. Noticeably, OsCaML2 (LOC_Os03g59770) was experimentally confirmed as a target gene of osa-miR1432, and the overexpressing OsCaML2 transgenic plants exhibited compromised resistance to Xoo infestation. Our results indicate that osa-miR1432 and OsCaML2 were differently responsive to Xoo invasion at the transcriptional level and fine-tune rice resistance to Xoo infection, which may be referable in resistance gene discovery and valuable in the pursuit of improving Xoo resistance in rice breeding.

Applying HPLC to Screening QTLs for BLB Resistance in Rice

Bacterial leaf blight (BLB) is caused by Xanthomonas oryzae pv. oryzae and is a major cause of rice yield reductions around the world. When diseased, plants produce a variety of metabolites to resist pathogens. In this study, the various defense metabolites were quantified using high-performance liquid chromatography (HPLC) after Xoo inoculation in a 120 Cheongcheong/Nagdong double haploid (CNDH) population. Quantitative trait locus (QTL) mapping was conducted using the concentration of the plant defense metabolites. HPLC analyzes the concentration of substances according to the severity of disease symptoms. Searching for BLB resistance candidate genes by applying this analysis method is very effective when mapping related genes. These resistance genes can be mapped directly to the causative pathogens. A total of 17 metabolites were detected by means of HPLC analysis after Xoo inoculation in the 120 CNDH population. QTL mapping of the metabolite concentrations resulted in the detection of the BLB resistance candidate gene, OsWRKYq6, in RM3343 of chromosome 6. OsWRKYq6 has a very high homology sequence with WRKY transcription factor 39, and when inoculated with Xoo, the relative expression level of the resistant population was higher than that of the susceptible population. Resistance genes have previously been detected using only phenotypic change data. In this study, resistance candidate genes were detected using the concentration of metabolites produced in plants after inoculation with pathogens. This newly developed analysis method can be used to effectively detect and identify genes directly involved in disease resistance for future studies.

Genome-Wide Association Study of QTLs Conferring Resistance to Bacterial Leaf Streak in Rice

Bacterial leaf streak (BLS) is a devastating rice disease caused by the bacterial pathogen, Xanthomonas oryzae pv. oryzicola (Xoc), which can result in severe damage to rice production worldwide. Based on a total of 510 rice accessions, trialed in two seasons and using six different multi-locus GWAS methods (mrMLM, ISIS EM-BLASSO, pLARmEB, FASTmrMLM, FASTmrEMMA and pKWmEB), 79 quantitative trait nucleotides (QTNs) reflecting 69 QTLs for BLS resistance were identified (LOD > 3). The QTNs were distributed on all chromosomes, with the most distributed on chromosome 11, followed by chromosomes 1 and 5. Each QTN had an additive effect of 0.20 (cm) and explained, on average, 2.44% of the phenotypic variance, varying from 0.00–0.92 (cm) and from 0.00–9.86%, respectively. Twenty-five QTNs were detected by at least two methods. Among them, qnBLS11.17 was detected by as many as five methods. Most of the QTNs showed a significant interaction with their environment, but no QTNs were detected in both seasons. By defining the QTL range for each QTN according to the LD half-decay distance, a total of 848 candidate genes were found for nine top QTNs. Among them, more than 10% were annotated to be related to biotic stress resistance, and five showed a significant response to Xoc infection. Our results could facilitate the in-depth study and marker-assisted improvement of rice resistance to BLS.

Spatial kernel models capturing field heterogeneity for accurate estimation of genetic potential

According to Fisher's principles, an experimental field is typically divided into multiple blocks for local control. Although homogeneity is supposed within a block, this assumption may not be practical for large blocks, such as those including hundreds of plots. In line evaluation trials, which are essential in plant breeding, field heterogeneity must be carefully treated, because it can cause bias in the estimation of genetic potential. To more accurately estimate genotypic values in a large field trial, we developed spatial kernel models incorporating genome-wide markers, which consider continuous heterogeneity within a block and over the field. In the simulation study, the spatial kernel models were robust under various conditions. Although heritability, spatial autocorrelation range, replication number, and missing plots directly affected the estimation accuracy of genotypic values, the spatial kernel models always showed superior performance over the classical block model. We also employed these spatial kernel models for quantitative trait locus mapping. Finally, using field experimental data of bioenergy sorghum lines, we validated the performance of the spatial kernel models. The results suggested that a spatial kernel model is effective for evaluating the genetic potential of lines in a heterogeneous field.

Identification of rice (Oryza sativa L.) genes involved in sheath blight resistance via a genome-wide association study

Rice sheath blight (RSB) is an economically significant disease affecting rice yield worldwide. Genetic resistance to RSB is associated with multiple minor genes, with each providing a minor phenotypic effect, but the underlying dominant resistance genes remain unknown. A genome-wide association study (GWAS) of 259 diverse rice varieties, with genotypes based on a single nucleotide polymorphism (SNP) and haplotype, was conducted to assess their sheath blight reactions at three developmental stages (seedlings, tillering and booting). A total of 653 genes were correlated with sheath blight resistance, of which the disease resistance protein RPM1 (OsRSR1) and protein kinase domain-containing protein (OsRLCK5) were validated by overexpression and knockdown assays. We further found that the coiled-coil (CC) domain of OsRSR1 (OsRSR1-CC) and full-length OsRLCK5 interacted with serine hydroxymethyltransferase 1 (OsSHM1) and glutaredoxin (OsGRX20), respectively. It was found that OsSHM1, which has a role in the reactive oxygen species (ROS) burst, and OsGRX20 enhanced the antioxidation ability of plants. A regulation model of the new RSB resistance though the glutathione (GSH)-ascorbic acid (AsA) antioxidant system was therefore revealed. These results enhance our understanding of RSB resistance mechanisms and provide better gene resources for the breeding of disease resistance in rice.

Genome-Wide Association Mapping and Gene Expression Analysis Reveal the Negative Role of OsMYB21 in Regulating Bacterial Blight Resistance in Rice

BACKGROUND

Bacterial blight (BB), caused by Xanthomonas oryzae pv. oryzae (Xoo), is one of the most devastating diseases in rice all over the world. Due to the diversity and rapid evolution of Xoo, identification and use of the non-race specific quantitative resistance QTLs has been considered the preferred strategy for effective control of this disease. Although numerous QTLs for BB resistance have been identified, they haven't been effectively used for improvement of BB resistance in rice due to their small effects and lack of knowledge on the function of genes underlying the QTLs.

RESULTS

In the present study, a genome-wide association study of BB resistance was performed in a rice core collection from South China. A total of 17 QTLs were identified to be associated with BB resistance. Among them, 13 QTLs were newly identified in the present study and the other 4 QTLs were co-localized with the previously reported QTLs or Xa genes that confer qualitative resistance to Xoo strains. Particularly, the qBBR11-4 on chromosome 11 explained the largest phenotypic variation in this study and was co-localized with the previously identified QTLs for BB and bacterial leaf streak (BLS) resistance against diverse strains in three studies, suggesting its broad-spectrum resistance and potential value in rice breeding. Through combined analysis of differential expression and annotations of the predicted genes within qBBR11-4 between two sets of rice accessions selected based on haplotypes and disease phenotypes, we identified the transcription factor OsMYB21 as the candidate gene for qBBR11-4. The OsMYB21 overexpressing plants exhibited decreased resistance to bacterial blight, accompanied with down-regulation of several defense-related genes compared with the wild-type plants.

CONCLUSION

The results suggest that OsMYB21 negatively regulates bacterial blight resistance in rice, and this gene can be a promising target in rice breeding by using the gene editing method. In addition, the potential candidate genes for the 13 novel QTLs for BB resistance were also analyzed in this study, providing a new source for cloning of genes associated with BB resistance and molecular breeding in rice.

Single nucleotide polymorphisms reveal genetic diversity in New Mexican chile peppers (Capsicum spp.)

BACKGROUND

Chile peppers (Capsicum spp.) are among the most important horticultural crops in the world due to their number of uses. They are considered a major cultural and economic crop in the state of New Mexico in the United States. Evaluating genetic diversity in current New Mexican germplasm would facilitate genetic improvement for different traits. This study assessed genetic diversity, population structure, and linkage disequilibrium (LD) among 165 chile pepper genotypes using single nucleotide polymorphism (SNP) markers derived from genotyping-by-sequencing (GBS).

RESULTS

A GBS approach identified 66,750 high-quality SNP markers with known map positions distributed across the 12 chromosomes of Capsicum. Principal components analysis revealed four distinct clusters based on species. Neighbor-joining phylogenetic analysis among New Mexico State University (NMSU) chile pepper cultivars showed two main clusters, where the C. annuum genotypes grouped together based on fruit or pod type. A Bayesian clustering approach for the Capsicum population inferred K = 2 as the optimal number of clusters, where the C. chinense and C. frutescens grouped in a single cluster. Analysis of molecular variance revealed majority of variation to be between the Capsicum species (76.08 %). Extensive LD decay (~ 5.59 Mb) across the whole Capsicum population was observed, demonstrating that a lower number of markers would be required for implementing genome wide association studies for different traits in New Mexican type chile peppers. Tajima's D values demonstrated positive selection, population bottleneck, and balancing selection for the New Mexico Capsicum population. Genetic diversity for the New Mexican chile peppers was relatively low, indicating the need to introduce new alleles in the breeding program to broaden the genetic base of current germplasm.

CONCLUSIONS

Genetic diversity among New Mexican chile peppers was evaluated using GBS-derived SNP markers and genetic relatedness on the species level was observed. Introducing novel alleles from other breeding programs or from wild species could help increase diversity in current germplasm. We present valuable information for future association mapping and genomic selection for different traits for New Mexican chile peppers for genetic improvement through marker-assisted breeding.

Genome-wide association mapping for resistance to bacterial blight and bacterial leaf streak in rice

Using genome-wide SNP association mapping, a total of 77 and 7 loci were identified for rice bacterial blight and bacterial leaf streak resistance, respectively, which may facilitate rice resistance improvement. Bacterial blight (BB) and bacterial leaf streak (BLS) caused by Gram-negative bacteria Xanthomonas oryzae pv. oryzae (Xoo) and X. oryzae pv. oryzicola (Xoc), respectively, are two economically important diseases negatively affecting rice production. To mine new sources of resistance, a set of rice germplasm collection consisting of 895 re-sequenced accessions from the 3000 Rice Genomes Project (3 K RGP) were screened for BB and BLS resistance under field conditions. Higher levels of BB resistance were observed in aus/boro subgroup, whereas the japonica, temperate japonica and tropical japonica subgroups possessed comparatively high levels of resistance to BLS. A genome-wide association study (GWAS) mined 77 genomic loci significantly associated with BB and 7 with BLS resistance. The phenotypic variance (R²) explained by these loci ranged from 0.4 to 30.2%. Among the loci, 7 for BB resistance were co-localized with known BB resistance genes and one for BLS resistance overlapped with a previously reported BLS resistance QTL. A search for the candidates in other novel loci revealed several defense-related genes that may be involved in resistance to BB and BLS. High levels of phenotypic resistance to BB or BLS could be attributed to the accumulation of the resistance (R) alleles at the associated loci, indicating their potential value in rice resistance breeding via gene pyramiding. The GWAS analysis validated the known genes underlying BB and BLS resistance and identified novel loci that could enrich the current resistance gene pool. The resources with strong resistance and significant SNPs identified in this study are potentially useful in breeding for BB and BLS resistance.

TALE-triggered and iTALE-suppressed Xa1-mediated resistance to bacterial blight is independent of rice transcription factor subunits OsTFIIAγ1 or OsTFIIAγ5

Xa1-mediated resistance to rice bacterial blight, caused by Xanthomonas oryzae pv. oryzae (Xoo), is triggered by transcription activator-like effectors (TALEs) and suppressed by interfering TALEs (iTALEs). TALEs interact with the rice transcription factor OsTFIIAγ1 or OsTFIIAγ5 (Xa5) to activate expression of target resistance and/or susceptibility genes. However, it is not clear whether OsTFIIAγ is involved in TALE-triggered and iTALE-suppressed Xa1-mediated resistance. In this study, genome-edited mutations in OsTFIIAγ5 or OsTFIIAγ1 of Xa1-containing rice 'IRBB1' and Xa1-transgenic plants of xa5-containing rice 'IRBB5' did not impair the activation or suppression of Xa1-mediated resistance. Correspondingly, the expression pattern of Xa1 in mutated OsTFIIAγ5 and OsTFIIAγ1 rice lines and 'IRBB1' rice was similar. In contrast, the expression of OsSWEET11 was repressed in rice lines mutated in OsTFIIAγ5 and OsTFIIAγ1. Bimolecular fluorescence complementation (BiFC) and co-immunoprecipitation assays showed that both TALE PthXo1 and iTALE Tal3a interacted with OsTFIIAγ1 and OsTFIIAγ5 in plant nuclei. These results indicated that TALE-triggered and iTALE-suppressed Xa1-mediated resistance to bacterial blight is independent of OsTFIIAγ1 or OsTFIIAγ5 in rice, and suggest that an unknown factor is potentially involved in the interaction of Xa1, TALEs and iTALEs.

Identification of Bacterial Blight Resistance Loci in Rice (Oryza sativa L.) against Diverse Xoo Thai Strains by Genome-Wide Association Study

Bacterial leaf blight (BLB) is a serious disease affecting global rice agriculture caused by Xanthomonas oryzae pv. oryzae (Xoo). Most resistant rice lines are dependent on single genes that are vulnerable to resistance breakdown caused by pathogen mutation. Here we describe a genome-wide association study of 222 predominantly Thai rice accessions assayed by phenotypic screening against 20 Xoo isolates. Loci corresponding to BLB resistance were detected using >142,000 SNPs. We identified 147 genes according to employed significance thresholds across chromosomes 1–6, 8, 9 and 11. Moreover, 127 of identified genes are located on chromosomal regions outside estimated Linkage Disequilibrium influences of known resistance genes, potentially indicating novel BLB resistance markers. However, significantly associated SNPs only occurred across a maximum of six Xoo isolates indicating that the development of broad-spectrum Xoo strain varieties may prove challenging. Analyses indicated a range of gene functions likely underpinning BLB resistance. In accordance with previous studies of accession panels focusing on indica varieties, our germplasm displays large numbers of SNPs associated with resistance. Despite encouraging data suggesting that many loci contribute to resistance, our findings corroborate previous inferences that multi-strain resistant varieties may not be easily realised in breeding programs without resorting to multi-locus strategies.

Apache Spark based kernelized fuzzy clustering framework for single nucleotide polymorphism sequence analysis

This paper introduces a kernel based fuzzy clustering approach to deal with the non-linear separable problems by applying kernel Radial Basis Functions (RBF) which maps the input data space non-linearly into a high-dimensional feature space. Discovering clusters in the high-dimensional genomics data is extremely challenging for the bioinformatics researchers for genome analysis. To support the investigations in bioinformatics, explicitly on genomic clustering, we proposed high-dimensional kernelized fuzzy clustering algorithms based on Apache Spark framework for clustering of Single Nucleotide Polymorphism (SNP) sequences. The paper proposes the Kernelized Scalable Random Sampling with Iterative Optimization Fuzzy c-Means (KSRSIO-FCM) which inherently uses another proposed Kernelized Scalable Literal Fuzzy c-Means (KSLFCM) clustering algorithm. Both the approaches completely adapt the Apache Spark cluster framework by localized sub-clustering Resilient Distributed Dataset (RDD) method. Additionally, we are also proposing a preprocessing approach for generating numeric feature vectors for huge SNP sequences and making it a scalable preprocessing approach by executing it on an Apache Spark cluster, which is applied to real-world SNP datasets taken from open-internet repositories of two different plant species, i.e., soybean and rice. The comparison of the proposed scalable kernelized fuzzy clustering results with similar works shows the significant improvement of the proposed algorithm in terms of time and space complexity, Silhouette index, and Davies-Bouldin index. Exhaustive experiments are performed on various SNP datasets to show the effectiveness of proposed KSRSIO-FCM in comparison with proposed KSLFCM and other scalable clustering algorithms, i.e., SRSIO-FCM, and SLFCM.

Genomics-assisted breeding for successful development of multiple-stress-tolerant, climate-smart rice for southern and southeastern Asia

Rice (Oryza sativa L.) in rainfed marginal environments is prone to multiple abiotic and biotic stresses, which can occur in combination in a single cropping season and adversely affect rice growth and yield. The present study was undertaken to develop high-yielding, climate-resilient rice that can provide tolerance to multiple biotic and abiotic stresses. An assembled first-crossing scheme was employed to transfer 15 quantitative trait loci (QTL) and genes-qDTY_1.1 , qDTY_2.1 , qDTY_3.1 , qDTY_12.1 (drought), Sub1 (submergence), Gm4 (gall midge), Pi9, Pita2 (blast), Bph3, Bph17 (brown plant hoppers), Xa4, xa5, xa13, Xa21, and Xa23 (bacterial leaf blight)-from eight different parents using genomics-assisted breeding. A funnel mating design was employed to assemble all the targeted QTL and genes into a high-yielding breeding line IR 91648-B-1-B-3-1. Gene-based linked markers were used in each generation from intercrossing to the F₆ generation for tracking the presence of desirable alleles of targeted QTL and genes. Single-plant selections were performed from F₂ onwards to select desirable recombinants possessing alleles of interest with suitable phenotypes. Phenotyping of 95 homozygous F₆ lines carrying six to 10 QTL and genes was performed for nonstress, reproductive-stage (RS) drought, blast, bacterial leaf blight (BLB), gall midge (GM), and for grain quality parameters such as chalkiness, amylose content (AC), gelatinization temperature (GT), and head rice recovery (HRR). Finally, 56 F₇ homozygous lines were found promising for multiple-location evaluation for grain yield (GY) and other traits. These multiple-stress-tolerant lines with the desired grain quality profiling can be targeted for varietal release in southern and southeastern Asia through national release systems.

Genome-Wide Association Study Dissects Resistance Loci against Bacterial Blight in a Diverse Rice Panel from the 3000 Rice Genomes Project

BACKGROUND

Bacterial blight (BB), caused by Xanthomonas oryzae pv. oryzae (Xoo) is one of the most devastating bacterial diseases of rice in temperate and tropical regions. Breeding and deployment of resistant cultivars carrying major resistance (R) genes has been the most effective approach for BB management. However, because of specific interaction of each R gene with the product of the corresponding pathogen avirulence or effector gene, new pathogen strains that can overcome the deployed resistance often emerge rapidly. To deal with ever-evolving Xoo, it is necessary to identify novel R genes and resistance quantitative trait loci (QTL).

RESULTS

BB resistance of a diverse panel of 340 accessions from the 3000 Rice Genomes Project (3 K RGP) was evaluated by artificial inoculation with four representative Xoo strains, namely Z173 (C4), GD1358 (C5), V from China and PXO339 (P9a) from Philippines. Using the 3 K RG 4.8mio filtered SNP Dataset, a total of 11 QTL associated with BB resistance on chromosomes 4, 5, 11 and 12 were identified through a genome-wide association study (GWAS). Among them, eight resistance loci, which were narrowed down to relatively small genomic intervals, coincided with previously reported QTL or R genes, e.g. xa5, xa25, xa44(t). The other three QTL were putative novel loci associated with BB resistance. Linear regression analysis showed a dependence of BB lesion length on the number of favorable alleles, suggesting that pyramiding QTL using marker-assisted selection would be an effective approach for improving resistance. In addition, the Hap2 allele of LOC_Os11g46250 underlying qC5-11.1 was validated as positively regulating resistance against strain C5.

CONCLUSIONS

Our findings provide valuable information for the genetic improvement of BB resistance and application of germplasm resources in rice breeding programs.

Genome-Wide Association Analysis Identifies Resistance Loci for Bacterial Leaf Streak Resistance in Rice (Oryza sativa L.)

Bacterial leaf streak (BLS) caused by Xanthomonas oryzae pv. oryzicola (Xoc) is one of the most devastating diseases in rice production areas, especially in humid tropical and subtropical zones throughout Asia and worldwide. A genome-wide association study (GWAS) analysis conducted on a collection of 236 diverse rice accessions, mainly indica varieties, identified 12 quantitative trait loci (QTLs) on chromosomes 1, 2, 3, 4, 5, 8, 9 and 11, conferring resistance to five representative isolates of Thai Xoc. Of these, five QTLs conferred resistance to more than one Xoc isolates. Two QTLs, qBLS5.1 and qBLS2.3, were considered promising QTLs for broad-spectrum resistance to BLS. The xa5 gene was proposed as a potential candidate gene for qBLS5.1 and three genes, encoding pectinesterase inhibitor (OsPEI), eukaryotic zinc-binding protein (OsRAR1), and NDP epimerase function, were proposed as candidate genes for qBLS2.3. Results from this study provide an insight into the potential QTLs and candidate genes for BLS resistance in rice. The recessive xa5 gene is suggested as a potential candidate for strong influence on broad-spectrum resistance and as a focal target in rice breeding programs for BLS resistance.

Xa1 Allelic R Genes Activate Rice Blight Resistance Suppressed by Interfering TAL Effectors

Xanthomonas oryzae pathovar oryzae (Xoo) uses transcription activator-like effectors (TALEs) to cause bacterial blight (BB) in rice. In turn, rice has evolved several mechanisms to resist BB by targeting TALEs. One mechanism involves the nucleotide-binding leucine-rich repeat (NLR) resistance gene Xa1 and TALEs. Reciprocally, Xoo has evolved TALE variants, C-terminally truncated versions (interfering TALEs or iTALEs), to overcome Xa1 resistance. However, it remains unknown to what extent the two co-adaptive mechanisms mediate Xoo-rice interactions. In this study, we cloned and characterized five additional Xa1 allelic R genes, Xa2, Xa31(t), Xa14, CGS-Xo1₁₁ , and Xa45(t) from a collection of rice accessions. Sequence analysis revealed that Xa2 and Xa31(t) from different rice cultivars are identical. These genes and their predicted proteins were found to be highly conserved, forming a group of Xa1 alleles. The XA1 alleles could be distinguished by the number of C-terminal tandem repeats consisting of 93 amino acid residues and ranged from four in XA14 to seven in XA45(t). Xa1 allelic genes were identified in the 3000 rice genomes surveyed. On the other hand, iTALEs could suppress the resistance mediated by Xa1 allelic R genes, and iTALE genes were prevalent (∼95%) in Asian, but not in African Xoo strains. Our findings demonstrate the prominence of a defense mechanism in which rice depends on Xa1 alleles and a counteracting mechanism in which Xoo relies on iTALEs for BB.

Genome-wide association study identifies an NLR gene that confers partial resistance to Magnaporthe oryzae in rice

Because of the frequent breakdown of major resistance (R) genes, identification of new partial R genes against rice blast disease is an important goal of rice breeding. In this study, we used a core collection of the Rice Diversity Panel II (C-RDP-II), which contains 584 rice accessions and are genotyped with 700 000 single-nucleotide polymorphism (SNP) markers. The C-RDP-II accessions were inoculated with three blast strains collected from different rice-growing regions in China. Genome-wide association study identified 27 loci associated with rice blast resistance (LABRs). Among them, 22 LABRs were not associated with any known blast R genes or QTLs. Interestingly, a nucleotide-binding site leucine-rich repeat (NLR) gene cluster exists in the LABR12 region on chromosome 4. One of the NLR genes is highly conserved in multiple partially resistant rice cultivars, and its expression is significantly up-regulated at the early stages of rice blast infection. Knockout of this gene via CRISPR-Cas9 in transgenic plants partially reduced blast resistance to four blast strains. The identification of this new non-strain specific partial R gene, tentatively named rice blast Partial Resistance gene 1 (PiPR1), provides genetic material that will be useful for understanding the partial resistance mechanism and for breeding durably resistant cultivars against blast disease of rice.

Marker Assisted Forward Breeding to Combine Multiple Biotic-Abiotic Stress Resistance/Tolerance in Rice

BACKGROUND

Unfavorable climatic changes have led to an increased threat of several biotic and abiotic stresses over the past few years. Looking at the massive damage caused by these stresses, we undertook a study to develop high yielding climate-resilient rice, using genes conferring resistance against blast (Pi9), bacterial leaf blight (BLB) (Xa4, xa5, xa13, Xa21), brown planthopper (BPH) (Bph3, Bph17), gall midge (GM) (Gm4, Gm8) and QTLs for drought tolerance (qDTY_1.1 and qDTY_3.1) through marker-assisted forward breeding (MAFB) approach.

RESULT

Seven introgression lines (ILs) possessing a combination of seven to ten genes/QTLs for different biotic and abiotic stresses have been developed using marker-assisted selection (MAS) breeding method in the background of Swarna with drought QTLs. These ILs were superior to the respective recurrent parent in agronomic performance and also possess preferred grain quality with intermediate to high amylose content (AC) (23-26%). Out of these, three ILs viz., IL1 (Pi9+ Xa4+ xa5+ Xa21+ Bph17+ Gm8+ qDTY_1.1+ qDTY_3.1), IL6 (Pi9+ Xa4+ xa5+ Xa21+ Bph3+ Bph17+ Gm4+ Gm8+ qDTY_1.1+ qDTY_3.1) and IL7 (Pi9+ Xa4+ xa5+ Bph3+ Gm4+ qDTY_1.1+ qDTY_3.1) had shown resistance\tolerance for multiple biotic and abiotic stresses both in the field and glasshouse conditions. Overall, the ILs were high yielding under various stresses and importantly they also performed well in non-stress conditions without any yield penalty.

CONCLUSION

The current study clearly illustrated the success of MAS in combining tolerance to multiple biotic and abiotic stresses while maintaining higher yield potential and preferred grain quality. Developed ILs with seven to ten genes in the current study showed superiority to recurrent parent Swarna+drought for multiple-biotic stresses (blast, BLB, BPH and GM) together with yield advantages of 1.0 t ha^- 1 under drought condition, without adverse effect on grain quality traits under non-stress.

A novel locus from the wild allotetraploid rice species Oryza latifolia Desv. confers bacterial blight (Xanthomonas oryzae pv. oryzae) resistance in rice (O. sativa)

Bacterial blight caused by Xanthomonas oryzae pv. oryzae (Xoo) is a major limiting factor to rice productivity worldwide. Genetic control through the identification of novel sources of bacterial blight resistance and their utilization in resistance breeding remains the most effective and economical strategy to manage the disease. Here we report the identification of a novel locus from the wild Oryza species, Oryza latifolia, conferring a race-specific resistance to Philippine Xoo race 9A (PXO339). The locus was identified from two introgression lines i.e. WH12-2252 and WH12-2256 that segregated from O. latifolia monosomic alien addition lines (MAALs). The discrete segregation ratio of susceptible and resistant phenotypes in the F2 (χ2[3:1] = 0.22 at p>0.05) and F3 (χ2[3:1] = 0.36 at p>0.05) populations indicates that PXO339 resistance in the MAAL-derived introgression lines (MDILs) is controlled by a single, recessive gene. Genotyping of a total of 216 F2, 1130 F3 and 288 F4 plants derived from crossing either of the MDILs with the recurrent parent used to generate the MAALs narrowed the candidate region to a 1,817 kb locus that extends from 10,425 to 12,266 kb in chromosome 12. Putative candidate genes that were identified by data mining and comparative sequence analysis can provide targets for further studies on mapping and cloning of the causal gene for PXO339 resistance in the MDILs. To our knowledge, this is the first report of a genetic locus from the allotetraploid wild rice, O. latifolia conferring race-specific resistance to bacterial blight.

Rice Stress-Resistant SNP Database

BACKGROUND

Rice (Oryza sativa L.) yield is limited inherently by environmental stresses, including biotic and abiotic stresses. Thus, it is of great importance to perform in-depth explorations on the genes that are closely associated with the stress-resistant traits in rice. The existing rice SNP databases have made considerable contributions to rice genomic variation information but none of them have a particular focus on integrating stress-resistant variation and related phenotype data into one web resource.

RESULTS

Rice Stress-Resistant SNP database (http://bioinformatics.fafu.edu.cn/RSRS) mainly focuses on SNPs specific to biotic and abiotic stress-resistant ability in rice, and presents them in a unified web resource platform. The Rice Stress-Resistant SNP (RSRS) database contains over 9.5 million stress-resistant SNPs and 797 stress-resistant candidate genes in rice, which were detected from more than 400 stress-resistant rice varieties. We incorporated the SNPs function, genome annotation and phenotype information into this database. Besides, the database has a user-friendly web interface for users to query, browse and visualize a specific SNP efficiently. RSRS database allows users to query the SNP information and their relevant annotations for individual variety or more varieties. The search results can be visualized graphically in a genome browser or displayed in formatted tables. Users can also align SNPs between two or more rice accessions.

CONCLUSION

RSRS database shows great utility for scientists to further characterize the function of variants related to environmental stress-resistant ability in rice.

1k-RiCA (1K-Rice Custom Amplicon) a novel genotyping amplicon-based SNP assay for genetics and breeding applications in rice

BACKGROUND

While a multitude of genotyping platforms have been developed for rice, the majority of them have not been optimized for breeding where cost, turnaround time, throughput and ease of use, relative to density and informativeness are critical parameters of their utility. With that in mind we report the development of the 1K-Rice Custom Amplicon, or 1k-RiCA, a robust custom sequencing-based amplicon panel of ~ 1000-SNPs that are uniformly distributed across the rice genome, designed to be highly informative within indica rice breeding pools, and tailored for genomic prediction in elite indica rice breeding programs.

RESULTS

Empirical validation tests performed on the 1k-RiCA showed average marker call rates of 95% with marker repeatability and concordance rates of 99%. These technical properties were not affected when two common DNA extraction protocols were used. The average distance between SNPs in the 1k-RiCA was 1.5 cM, similar to the theoretical distance which would be expected between 1,000 uniformly distributed markers across the rice genome. The average minor allele frequencies on a panel of indica lines was 0.36 and polymorphic SNPs estimated on pairwise comparisons between indica by indica accessions and indica by japonica accessions were on average 430 and 450 respectively. The specific design parameters of the 1k-RiCA allow for a detailed view of genetic relationships and unambiguous molecular IDs within indica accessions and good cost vs. marker-density balance for genomic prediction applications in elite indica germplasm. Predictive abilities of Genomic Selection models for flowering time, grain yield, and plant height were on average 0.71, 0.36, and 0.65 respectively based on cross-validation analysis. Furthermore the inclusion of important trait markers associated with 11 different genes and QTL adds value to parental selection in crossing schemes and marker-assisted selection in forward breeding applications.

CONCLUSIONS

This study validated the marker quality and robustness of the 1k-RiCA genotypic platform for genotyping populations derived from indica rice subpopulation for genetic and breeding purposes including MAS and genomic selection. The 1k-RiCA has proven to be an alternative cost-effective genotyping system for breeding applications.

Methodology: ssb-MASS: a single seed-based sampling strategy for marker-assisted selection in rice

BACKGROUND

Integrated breeding approaches such as combining marker-assisted selection and rapid line fixation through single-seed-descent, can effectively increase the frequency of desirable alleles in a breeding program and increase the rate of genetic gain for quantitative traits by shortening the breeding cycle. However, with most genotyping being outsourced to 3rd party service providers' nowadays, sampling has become the bottleneck for many breeding programs. While seed-chipping as prevailed as an automatable seed sampling protocol in many species, the symmetry of rice seeds makes this solution as laborious and costly as sampling leaf tissue. The aim of this study is to develop, validate and deploy a single seed sampling strategy for marker-assisted selection of fixed lines in rice that is more efficient, cost-effective and convenient compared to leaf-based sampling protocols without compromising the accuracy of the marker-assisted selection results.

RESULTS

Evaluations replicated across accessions and markers showed that a single rice seed is sufficient to generate enough DNA (7-8 ng/μL) to run at least ten PCR trait-markers suitable for marker-assisted selection strategies in rice. The DNA quantity and quality extracted from single seeds from fixed lines (F₆) with different physical and/or chemical properties were not significantly different. Nor were there significant differences between single seeds collected 15 days after panicle initiation compared to those harvested at maturity. A large-scale comparison between single seed and leaf-based methodologies showed not only high levels of genotypic concordance between both protocols (~ 99%) but also higher SNP call rates in single seed (99.24% vs. 97.5% in leaf). A cost-benefit analysis showed that this single seed sampling strategy decreased the cost of sampling fourfold. An advantage of this approach is that desirable genotypes can be selected before investing in planting activities reducing the cost associated with field operations.

CONCLUSION

This study reports the development of a cost-effective and simple single seed genotyping strategy that facilitates the adoption and deployment of marker-assisted selection strategies in rice. This will allow breeders to increase the frequency of favorable alleles and combine rapid generation advancement techniques much more cost-effectively accelerating the process and efficiency of parental selection and varietal development.

Identification of new loci for salt tolerance in soybean by high-resolution genome-wide association mapping

BACKGROUND

Salinity is an abiotic stress that negatively affects soybean [Glycine max (L.) Merr.] seed yield. Although a major gene for salt tolerance was identified and consistently mapped to chromosome (Chr.) 3 by linkage mapping studies, it does not fully explain genetic variability for tolerance in soybean germplasm. In this study, a genome-wide association study (GWAS) was performed to map genomic regions for salt tolerance in a diverse panel of 305 soybean accessions using a single nucleotide polymorphism (SNP) dataset derived from the SoySNP50K iSelect BeadChip. A second GWAS was also conducted in a subset of 234 accessions using another 3.7 M SNP dataset derived from a whole-genome resequencing (WGRS) study. In addition, three gene-based markers (GBM) of the known gene, Glyma03g32900, on Chr. 3 were also integrated into the two datasets. Salt tolerance among soybean lines was evaluated by leaf scorch score (LSS), chlorophyll content ratio (CCR), leaf sodium content (LSC), and leaf chloride content (LCC).

RESULTS

For both association studies, a major locus for salt tolerance on Chr. 3 was confirmed by a number of significant SNPs, of which three gene-based SNP markers, Salt-20, Salt14056 and Salt11655, had the highest association with all four traits studied. Also, additional genomic regions on Chrs. 1, 8, and 18 were found to be associated with various traits measured in the second GWAS using the WGRS-derived SNP dataset.

CONCLUSIONS

A region identified on Chr. 8 was identified to be associated with all four traits and predicted as a new minor locus for salt tolerance in soybean. The candidate genes harbored in this minor locus may help reveal the molecular mechanism involved in salt tolerance and to improve tolerance in soybean cultivars. The significant SNPs will be useful for marker-assisted selection for salt tolerance in soybean breeding programs.

Identification of novel recessive gene xa44(t) conferring resistance to bacterial blight races in rice by QTL linkage analysis using an SNP chip

KEY MESSAGE

Using QTL analysis and fine mapping, the novel recessive gene xa44(t) conferring resistance to BB was identified and the expression level of the gene was confirmed through qRT-PCR analysis. Bacterial blight (BB) disease caused by Xanthomonas oryzae pv. oryzae (Xoo) is a major factor causing rice yield loss in most rice-cultivating countries, especially in Asia. The deployment of cultivars with resistance to BB is the most effective method to control the disease. However, the evolution of new Xoo or pathotypes altered by single-gene-dependent mutations often results in breakdown of resistance. Thus, efforts to identify novel R-genes with sustainable BB resistance are urgently needed. In this study, we identified three quantitative trait loci (QTLs) on chromosomes 1, 4, and 11, from an F2 population of 493 individuals derived from a cross between IR73571-3B-11-3-K3 and Ilpum using a 7K SNP chip. Of these QTLs, one major QTL, qBB_11, on chromosome 11 explained 61.58% of the total phenotypic variance in the population, with an LOD value of 113.59, based on SNPs 11964077 and 11985463. The single major R-gene, with recessive gene action, was designated xa44(t) and was narrowed down to a 120-kb segment flanked within 28.00 Mbp to 28.12 Mbp. Of nine ORFs present in the target region, two ORFs revealed significantly different expression levels of the candidate genes. These candidate genes (Os11g0690066 and Os11g0690466) are described as "serine/threonine protein kinase domain containing protein" and "hypothetical protein," respectively. The results will be useful to further understand BB resistance mechanisms and provide new sources of resistance, together with DNA markers for MAS breeding to improve BB resistance in rice.