Bulked Segregant Analysis Coupled with Whole-Genome Sequencing (BSA-Seq) Mapping Identifies a Novel pi21 Haplotype Conferring Basal Resistance to Rice Blast Disease

Haplotypes at the sorghum ARG4 and ARG5 NLR loci confer resistance to anthracnose

Sorghum anthracnose caused by the fungus Colletotrichum sublineola (Cs) is a damaging disease of the crop. Here, we describe the identification of ANTHRACNOSE RESISTANCE GENES (ARG4 and ARG5) encoding canonical nucleotide-binding leucine-rich repeat (NLR) receptors. ARG4 and ARG5 are dominant resistance genes identified in the sorghum lines SAP135 and P9830, respectively, that show broad-spectrum resistance to Cs. Independent genetic studies using populations generated by crossing SAP135 and P9830 with TAM428, fine mapping using molecular markers, comparative genomics and gene expression studies determined that ARG4 and ARG5 are resistance genes against Cs strains. Interestingly, ARG4 and ARG5 are both located within clusters of duplicate NLR genes at linked loci separated by ~1 Mb genomic region. SAP135 and P9830 each carry only one of the ARG genes while having the recessive allele at the second locus. Only two copies of the ARG5 candidate genes were present in the resistant P9830 line while five non-functional copies were identified in the susceptible line. The resistant parents and their recombinant inbred lines carrying either ARG4 or ARG5 are resistant to strains Csgl1 and Csgrg suggesting that these genes have overlapping specificities. The role of ARG4 and ARG5 in resistance was validated through sorghum lines carrying independent recessive alleles that show increased susceptibility. ARG4 and ARG5 are located within complex loci displaying interesting haplotype structures and copy number variation that may have resulted from duplication. Overall, the identification of anthracnose resistance genes with unique haplotype stucture provides a foundation for genetic studies and resistance breeding.

Identification of Clubroot (Plasmodiophora brassicae) Resistance Loci in Chinese Cabbage (Brassica rapa ssp. pekinensis) with Recessive Character

The soil-borne pathogen Plasmodiophora brassicae is the causal agent of clubroot, a major disease in Chinese cabbage (Brassica rapa ssp. pekinensis). The host's resistance genes often confer immunity to only specific pathotypes and may be rapidly overcome. Identification of novel clubroot resistance (CR) from germplasm sources is necessary. In this study, Bap246 was tested by being crossed with different highly susceptible B. rapa materials and showed recessive resistance to clubroot. An F2 population derived from Bap246 × Bac1344 was used to locate the resistance Quantitative Trait Loci (QTL) by Bulk Segregant Analysis Sequencing (BSA-Seq) and QTL mapping methods. Two QTL on chromosomes A01 (4.67-6.06 Mb) and A08 (10.42-11.43 Mb) were found and named Cr4Ba1.1 and Cr4Ba8.1, respectively. Fifteen and eleven SNP/InDel markers were used to narrow the target regions in the larger F2 population to 4.67-5.17 Mb (A01) and 10.70-10.84 Mb (A08), with 85 and 19 candidate genes, respectively. The phenotypic variation explained (PVE) of the two QTL were 30.97% and 8.65%, respectively. Combined with gene annotation, mutation site analysis, and real-time quantitative polymerase chain reaction (qRT-PCR) analysis, one candidate gene in A08 was identified, namely Bra020861. And an insertion and deletion (InDel) marker (co-segregated) named Crr1-196 was developed based on the gene sequence. Bra013275, Bra013299, Bra013336, Bra013339, Bra013341, and Bra013357 in A01 were the candidate genes that may confer clubroot resistance in Chinese cabbage. The resistance resource and the developed marker will be helpful in Brassica breeding programs.

BSA-Seq for the Identification of Major Genes for EPN in Rice

Improving rice yield is one of the most important food issues internationally. It is an undeniable goal of rice breeding, and the effective panicle number (EPN) is a key factor determining rice yield. Increasing the EPN in rice is a major way to increase rice yield. Currently, the main quantitative trait locus (QTL) for EPN in rice is limited, and there is also limited research on the gene for EPN in rice. Therefore, the excavation and analysis of major genes related to EPN in rice is of great significance for molecular breeding and yield improvement. This study used japonica rice varieties Dongfu 114 and Longyang 11 to construct an F5 population consisting of 309 individual plants. Two extreme phenotypic pools were constructed by identifying the EPN of the population, and QTL-seq analysis was performed to obtain three main effective QTL intervals for EPN. This analysis also helped to screen out 34 candidate genes. Then, EPN time expression pattern analysis was performed on these 34 genes to screen out six candidate genes with higher expression levels. Using a 3K database to perform haplotype analysis on these six genes, we selected haplotypes with significant differences in EPN. Finally, five candidate genes related to EPN were obtained.

The E3 ubiquitin ligase OsRGLG5 targeted by the Magnaporthe oryzae effector AvrPi9 confers basal resistance against rice blast

Rice blast, caused by Magnaporthe oryzae, is one of the most devastating diseases of rice. During infection, M. oryzae secretes effectors to facilitate blast development. Among these effectors, the avirulence factor AvrPi9 is recognized by Pi9, a broad-spectrum blast resistance protein that triggers Pi9-mediated resistance in rice. However, little is known about the interaction between AvrPi9 and Pi9 and how AvrPi9 exerts virulence to promote infection. In this study, we found that ectopic expression of AvrPi9 in the Pi9-lacking cultivar TP309 suppressed basal resistance against M. oryzae. Furthermore, we identified an AvrPi9-interacting protein in rice, which we named OsRGLG5, encoding a functional RING-type E3 ubiquitin ligase. During infection, AvrPi9 was ubiquitinated and degraded by OsRGLG5. Meanwhile, AvrPi9 affected the stability of OsRGLG5. Infection assays revealed that OsRGLG5 is a positive regulator of basal resistance against M. oryzae, but it is not essential for Pi9-mediated blast resistance in rice. In conclusion, our results revealed that OsRGLG5 is targeted by the M. oryzae effector AvrPi9 and positively regulates basal resistance against rice blast.

Phenotypic analysis of Longya-10 × pale flax hybrid progeny and identification of candidate genes regulating prostrate/erect growth in flax plants

Flax is a dual-purpose crop that is important for oil and fiber production. The growth habit is one of the crucial targets of selection during flax domestication. Wild hybridization between cultivated flax and wild flax can produce superior germplasms for flax breeding and facilitate the study of the genetic mechanism underlying agronomically important traits. In this study, we used pale flax, Linum grandiflorum, and L. perenne to pollinate Longya-10. Only pale flax interspecific hybrids were obtained, and the trait analysis of the F₁ and F₂ generations showed that the traits analyzed in this study exhibited disparate genetic characteristics. In the F₁ generation, only one trait, i.e., the number of capsules per plant (140) showed significant heterosis, while the characteristics of other traits were closely associated with those of the parents or a decline in hybrid phenotypes. The traits of the F₂ generation were widely separated, and the variation coefficient ranged from 9.96% to 146.15%. The quantitative trait locus underlying growth habit was preliminarily found to be situated on chromosome 2 through Bulked-segregant analysis sequencing. Then linkage mapping analysis was performed to fine-map GH2.1 to a 23.5-kb interval containing 4 genes. Among them, L.us.o.m.scaffold22.109 and L.us.o.m.scaffold22.112 contained nonsynonymous SNPs with Δindex=1. Combined with the qRT-PCR results, the two genes might be possible candidate genes for GH2.1. This study will contribute to the development of important germplasms for flax breeding, which would facilitate the elucidation of the genetic mechanisms regulating the growth habit and development of an ideal architecture for the flax plant.

BnaCRCs with domestication preference positively correlate with the seed-setting rate of canola

Canola (Brassica napus) is an important oil crop worldwide. The seed-setting rate (SS) is a critical factor in determining its yield, and the development of pistils affects pollination and seed sets. However, research on seed-setting defects has been limited owing to difficulties in the identification of phenotypes, mutations, and complex genetic mechanisms. In this study, we found a stigma defect (sd) mutant in B. napus, which had no nectary. The SS of sd mutants in the field was approximately 93.4% lower than that of the wild type. Scanning and transmission electron microscopy imaging of sd mutants showed a low density of stigma papillary cells and stigma papillary cell vacuoles that disappeared 16 h after flowering. Genetic analysis of segregated populations showed that two recessive nuclear genes are responsible for the mutant phenotype of sd. Based on re-sequencing and map-based cloning, we reduced the candidate sites on ChrA07 (BnaSSA07) and ChrC06 (BnaSSC06) to 30 and 67 kb, including six and eight predicted genes, respectively. Gene analyses showed that a pair of CRABS CLAW (CRC) homeologous genes at BnaSSA07 and BnaSSC06 were associated with the development of carpel and nectary. BnaSSA07.CRC and BnaSSC06.CRC candidate genes were found to be expressed in flower organs only, with significant differences in their expression in the pistils of the near-isogenic lines. DNA sequencing showed transposon insertions in the upstream region and intron of the candidate gene BnaSSA07.crc. We also found that BnaSSC06.crc exists widely in the natural population and we give possible reasons for its widespread existence.

Harnessing the potential of bulk segregant analysis sequencing and its related approaches in crop breeding

Most plant traits are governed by polygenes including both major and minor genes. Linkage mapping and positional cloning have contributed greatly to mapping genomic loci controlling important traits in crop species. However, they are low-throughput, time-consuming, and have low resolution due to which their efficiency in crop breeding is reduced. In this regard, the bulk segregant analysis sequencing (BSA-seq) and its related approaches, viz., quantitative trait locus (QTL)-seq, bulk segregant RNA-Seq (BSR)-seq, and MutMap, have emerged as efficient methods to identify the genomic loci/QTLs controlling specific traits at high resolution, accuracy, reduced time span, and in a high-throughput manner. These approaches combine BSA with next-generation sequencing (NGS) and enable the rapid identification of genetic loci for qualitative and quantitative assessments. Many previous studies have shown the successful identification of the genetic loci for different plant traits using BSA-seq and its related approaches, as discussed in the text with details. However, the efficiency and accuracy of the BSA-seq depend upon factors like sequencing depth and coverage, which enhance the sequencing cost. Recently, the rapid reduction in the cost of NGS together with the expected cost reduction of third-generation sequencing in the future has further increased the accuracy and commercial applicability of these approaches in crop improvement programs. This review article provides an overview of BSA-seq and its related approaches in crop breeding together with their merits and challenges in trait mapping.

GmIAA27 Encodes an AUX/IAA Protein Involved in Dwarfing and Multi-Branching in Soybean

Soybean plant height and branching affect plant architecture and yield potential in soybean. In this study, the mutant dmbn was obtained by treating the cultivar Zhongpin 661 with ethylmethane sulfonate. The dmbn mutant plants were shorter and more branched than the wild type. The genetic analysis showed that the mutant trait was controlled by a semi-dominant gene. The candidate gene was fine-mapped to a 91 kb interval on Chromosome 9 by combining BSA-seq and linkage analysis. Sequence analysis revealed that Glyma.09g193000 encoding an Aux/IAA protein (GmIAA27) was mutated from C to T in the second exon of the coding region, resulting to amino acid substitution of proline to leucine. Overexpression of the mutant type of this gene in Arabidopsis thaliana inhibited apical dominance and promoted lateral branch development. Expression analysis of GmIAA27 and auxin response genes revealed that some GH3 genes were induced. GmIAA27 relies on auxin to interact with TIR1, whereas Gmiaa27 cannot interact with TIR1 owing to the mutation in the degron motif. Identification of this unique gene that controls soybean plant height and branch development provides a basis for investigating the mechanisms regulating soybean plant architecture development.

Identification of qGL3.5, a Novel Locus Controlling Grain Length in Rice Through Bulked Segregant Analysis and Fine Mapping

Grain length (GL) directly affects the yield and quality of rice. Very few cloned GL-related genes are applied in production because their yield-increasing effects are not obvious, and the overall regulatory networks underlying the associated processes remain poorly understood. DNA samples from two bulk DNA pools (L-pool and S-pool) and their parents (KJ01 and Huaye 3) were subjected to high-throughput sequencing. Using bulked segregant analysis (BSA), qGL3.5 was mapped to a 0.34-Mb "hotspot" region on chromosome 3 that contains 37 genes related to various traits. Then, qGL3.5 was mapped to the genomic interval between the flanking markers M2 and M3 using 2786 BC₄F₂ individuals. Because the region from B5 to B6 was not the associated region under BSA-seq analysis, qGL3.5 was narrowed down to the interval between B6 and M3, which spanned 24.0-kb. Of all 37 genes with non-synonymous single-nucleotide polymorphisms (SNPs) between KJ01 and Huaye 3 based on BSA-seq analysis, only one complete annotated gene, ORF18 (Gene ID: LOC_Os03g42790.1) was found. ORF18 encodes an IBR-RING zinc-finger-related protein, with one really interesting new gene (RING) and two in between ring finger (IBR) domains. The knockout of ORF18 derived from Huaye 3 using clustered, regularly interspaced, short palindromic repeat (CRISPR)/CRISPR-associated 9 (Cas9) editing technology increased the GL of the mutant by approximately 2.2 mm. The novel locus qGL3.5 negatively regulated GL by promoting elongation of the longitudinal cell of the grain outer glume. These results provide a new genetic resource for rice grain shape breeding and a starting point for the functional characterization of the wild rice GL gene.

Bulk segregation analysis in the NGS era: a review of its teenage years

Bulk segregation analysis (BSA) utilizes a strategy of pooling individuals with extreme phenotypes to conduct economical and rapidly linked marker screening or quantitative trait locus (QTL) mapping. With the development of next-generation sequencing (NGS) technology in the past 10 years, BSA methods and technical systems have been gradually developed and improved. At the same time, the ever-decreasing costs of sequencing accelerate NGS-based BSA application in different species, including eukaryotic yeast, grain crops, economic crops, horticultural crops, trees, aquatic animals, and insects. This paper provides a landscape of BSA methods and reviews the BSA development process in the past decade, including the sequencing method for BSA, different populations, different mapping algorithms, associated region threshold determination, and factors affecting BSA mapping. Finally, we summarize related strategies in QTL fine mapping combining BSA.

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.

Identification, Analysis, and Confirmation of Seed Storability-Related Loci in Dongxiang Wild Rice (Oryza rufipogon Griff.)

Dongxiang wild rice (Oryza rufipogon Griff.) (DXWR) has strong seed storability and identifying its elite gene resources may facilitate genetic improvements in rice seed storability. In this study, we developed two backcross inbred lines (BILs) populations, with DXWR as a common donor parent and two rice varieties (F6 and R974) as recipient parents. Bulked segregant analysis via whole genome sequencing (BSA-seq) was used to identify seed storability-related loci in the DXWR and F6 population. Two main genomic regions containing 18,550,000–20,870,000 bp on chromosome 4 and 7,860,000–9,780,000 bp on chromosome 9 were identified as candidate loci of DXWR seed storability; these overlapped partially with seed storability-related quantitative trait loci (QTLs) discovered in previous studies, suggesting that these loci may provide important regions for isolating the responsible genes. In total, 448 annotated genes were predicted within the identified regions, of which 274 and 82 had nonsynonymous and frameshift mutations, respectively. We detected extensive metabolic activities and cellular processes during seed storability and confirmed the effects of the seed storability-related candidate loci using four BILs from DXWR and R974. These results may facilitate the cloning of DXWR seed storability-related genes, thereby elucidating rice seed storability and its improvement potential.

Molecular Analysis Uncovers the Mechanism of Fertility Restoration in Temperature-Sensitive Polima Cytoplasmic Male-Sterile Brassica napus

Temperature-sensitive male sterility is a heritable agronomic trait affected by genotype-environment interactions. In rapeseed (Brassica napus), Polima (pol) temperature-sensitive cytoplasmic male sterility (TCMS) is commonly used for two-line breeding, as the fertility of pol TCMS lines can be partially restored at certain temperatures. However, little is known about the underlying molecular mechanism that controls fertility restoration. Therefore, we aimed to investigate the fertility conversion mechanism of the pol TCMS line at two different ambient temperatures (16 °C and 25 °C). Our results showed that the anthers developed and produced vigorous pollen at 16 °C but not at 25 °C. In addition, we identified a novel co-transcript of orf224-atp6 in the mitochondria that might lead to fertility conversion of the pol TCMS line. RNA-seq analysis showed that 1637 genes were significantly differentially expressed in the fertile flowers of 596-L when compared to the sterile flower of 1318 and 596-H. Detailed analysis revealed that differentially expressed genes were involved in temperature response, ROS accumulation, anther development, and mitochondrial function. Single-molecule long-read isoform sequencing combined with RNA sequencing revealed numerous genes produce alternative splicing transcripts at high temperatures. Here, we also found that alternative oxidase, type II NAD(P)H dehydrogenases, and transcription factor Hsfs might play a crucial role in male fertility under the low-temperature condition. RNA sequencing and bulked segregant analysis coupled with whole-genome sequencing identified the candidate genes involved in the post-transcriptional modification of orf224. Overall, our study described a putative mechanism of fertility restoration in a pol TCMS line controlled by ambient temperature that might help utilise TCMS in the two-line breeding of Brassica crops.

Recent Progress in Rice Broad-Spectrum Disease Resistance

Rice is one of the most important food crops in the world. However, stable rice production is constrained by various diseases, in particular rice blast, sheath blight, bacterial blight, and virus diseases. Breeding and cultivation of resistant rice varieties is the most effective method to control the infection of pathogens. Exploitation and utilization of the genetic determinants of broad-spectrum resistance represent a desired way to improve the resistance of susceptible rice varieties. Recently, researchers have focused on the identification of rice broad-spectrum disease resistance genes, which include R genes, defense-regulator genes, and quantitative trait loci (QTL) against two or more pathogen species or many isolates of the same pathogen species. The cloning of broad-spectrum disease resistance genes and understanding their underlying mechanisms not only provide new genetic resources for breeding broad-spectrum rice varieties, but also promote the development of new disease resistance breeding strategies, such as editing susceptibility and executor R genes. In this review, the most recent advances in the identification of broad-spectrum disease resistance genes in rice and their application in crop improvement through biotechnology approaches during the past 10 years are summarized.

Rapid identification of mutations caused by fast neutron bombardment in Medicago truncatula

BACKGROUND

Fast neutron bombardment (FNB) is a very effective approach for mutagenesis and has been widely used in generating mutant libraries in many plant species. The main type of mutations of FNB mutants are deletions of DNA fragments ranging from few base pairs to several hundred kilobases, thus usually leading to the null mutation of genes. Despite its efficiency in mutagenesis, identification of the mutation sites is still challenging in many species. The traditional strategy of positional cloning is very effective in identifying the mutation but time-consuming. With the availability of genome sequences, the array-based comparative genomic hybridization (CGH) method has been developed to detect the mutation sites by comparing the signal intensities of probes between wild-type and mutant plants. Though CGH method is effective in detecting copy number variations (CNVs), the resolution and coverage of CGH probes are not adequate to identify mutations other than CNVs.

RESULTS

We report a new strategy and pipeline to sensitively identify the mutation sites of FNB mutants by combining deep-coverage whole-genome sequencing (WGS), polymorphism calling, and customized filtering in Medicago truncatula. Initially, we performed a bulked sequencing for a FNB white nodule (wn) mutant and its wild-type like plants derived from a backcross population. Following polymorphism calling and filtering, validation by manual check and Sanger sequencing, we identified that SymCRK is the causative gene of white nodule mutant. We also sequenced an individual FNB mutant yellow leaves 1 (yl1) and wild-type plant. We identified that ETHYLENE-DEPENDENT GRAVITROPISM-DEFICIENT AND YELLOW-GREEN 1 (EGY1) is the candidate gene for M. truncatula yl1 mutant.

CONCLUSION

Our results demonstrated that the method reported here is rather robust in identifying the mutation sites for FNB mutants.

Fine Mapping and Candidate Gene Identification of a White Flower Gene BrWF3 in Chinese Cabbage (Brassica rapa L. ssp. pekinensis)

Flower color is an important trait in plants. However, genes responsible for the white flower trait in Chinese cabbage are rarely reported. In this study, we constructed an F₂ population derived from the Y640-288 (white flower) and Y641-87 (yellow flower) lines for the fine mapping of the white flower gene BrWF3 in Chinese cabbage. Genetic analysis indicated that BrWF3 was controlled by a single recessive gene. Using BSA-seq and KASP assays, BrWF3 was fine-mapped to an interval of 105.6 kb. Functional annotation, expression profiling, and sequence variation analyses confirmed that the AtPES2 homolog, Bra032957, was the most likely candidate gene for BrWF3. Carotenoid profiles and transmission electron microscopy analysis suggested that BrWF3 might participate in the production of xanthophyll esters (particularly violaxanthin esters), which in turn disrupt chromoplast development and the formation of plastoglobules (PGs). A SNP deletion in the third exon of BrWF3 caused the loss of protein function, and interfered with the normal assembly of PGs, which was associated with reduced expression levels of genes involved in carotenoid metabolism. Furthermore, we developed and validated the functional marker TXBH83 for BrWF3. Our results provide insight into the molecular mechanism underlying flower color pigmentation and reveal valuable information for marker-assisted selection (MAS) breeding in Chinese cabbage.

Rapid Mining of Candidate Genes for Verticillium Wilt Resistance in Cotton Based on BSA-Seq Analysis

Cotton is a globally important cash crop. Verticillium wilt (VW) is commonly known as "cancer" of cotton and causes serious loss of yield and fiber quality in cotton production around the world. Here, we performed a BSA-seq analysis using an F_2:3 segregation population to identify the candidate loci involved in VW resistance. Two QTLs (qvw-D05-1 and qvw-D05-2) related to VW resistance in cotton were identified using two resistant/susceptible bulks from the F₂ segregation population constructed by crossing the resistant cultivar ZZM2 with the susceptible cultivar J11. A total of 30stop-lost SNPs and 42 stop-gained SNPs, which included 17 genes, were screened in the qvw-D05-2 region by SnpEff analysis. Further analysis of the transcriptome data and qRT-PCR revealed that the expression level of Ghir_D05G037630 (designated as GhDRP) varied significantly at certain time points after infection with V. dahliae. The virus-induced gene silencing of GhDRP resulted in higher susceptibility of the plants to V. dahliae than the control, suggesting that GhDRP is involved in the resistance to V. dahlia infection. This study provides a method for rapid mining of quantitative trait loci and screening of candidate genes, as well as enriches the genomic information and gene resources for the molecular breeding of disease resistance in cotton.