Next-generation sequencing to identify candidate genes and develop diagnostic markers for a novel Phytophthora resistance gene, RpsHC18, in soybean

Fine mapping and identification of a Fusarium wilt resistance gene FwS1 in pea

KEY MESSAGE

A Fusarium wilt resistance gene FwS1 on pea chromosome 6 was identified and mapped to a 91.4 kb region by a comprehensive genomic-based approach, and the gene Psat6g003960 harboring NB-ARC domain was identified as the putative candidate gene. Pea Fusarium wilt, incited by Fusarium oxysporum f. sp. pisi (Fop), has always been a devastating disease that causes severe yield losses and economic damage in pea-growing regions worldwide. The utilization of pea cultivars carrying resistance gene is the most efficient approach for managing this disease. In order to finely map resistance gene, F2 populations were established through the cross between Shijiadacaiwan 1 (resistant) and Y4 (susceptible). The resistance genetic analysis indicated that the Fop resistance in Shijiadacaiwan 1 was governed by a single dominant gene, named FwS1. Based on the bulked segregant analysis sequencing analyses, the gene FwS1 was initially detected on chromosome 6 (i.e., linking group II, chr6LG2), and subsequent linkage mapping with 589 F2 individuals fine-mapped the gene FwS1 into a 91.4 kb region. The further functional annotation and haplotype analysis confirmed that the gene Psat6g003960, characterized by a NB-ARC (nucleotide-binding adaptor shared by APAF-1, R proteins, and CED-4) domain, was considered as the most promising candidate gene. The encoding amino acids were altered by a "T/C" single-nucleotide polymorphism (SNP) in the first exon of the Psat6g003960, and based on this SNP locus, the molecular marker A016180 was determined to be a diagnostic marker for FwS1 by validating its specificity in both pea accessions and genetic populations with different genetic backgrounds. The FwS1 with diagnostic KASP marker A016180 could facilitate marker-assisted selection in resistance pea breeding in pea. In addition, a comparison of the candidate gene Psat6g003960 in 74SN3B and SJ1 revealed the same sequences. This finding indicated that 74SN3B carried the candidate gene for FwS1, suggesting that FwS1 and Fwf may be closely linked or an identical resistant gene against Fusarium wilt.

Identification of major QTLs for drought tolerance in soybean, together with a novel candidate gene, GmUAA6

Drought tolerance is a complex trait in soybean that is controlled by polygenetic quantitative trait loci (QTLs). In this study, wilting score, days-to-wilting, leaf relative water content, and leaf relative conductivity were used to identify QTLs associated with drought tolerance in recombinant inbred lines derived from a cross between a drought-sensitive variety, Lin, and a drought-tolerant variety, Meng. A total of 33 drought-tolerance QTLs were detected. Of these 17 were major QTLs. In addition, 15 were novel drought-tolerance QTLs. The most predominant QTL was on chromosome 11. This was detected in at least three environments. The overlapped mapping interval of the four measured traits was 0.2 cM in genetic distance (about 220 kb in physical length). Glyma.11g143500 (designated as GmUAA6), which encodes a UDP-N-acetylglucosamine transporter, was identified as the most likely candidate gene. The allele of GmUAA6 from Lin (GmUAA6Lin) was associated with improved soybean drought tolerance. Overexpression of GmUAA6Lin in Arabidopsis and soybean hairy roots enhanced drought tolerance. Furthermore, a 3-bp insertion/deletion (InDel) in the coding sequence of GmUAA6 explained up to 49.9% of the phenotypic variation in drought tolerance-related traits, suggesting that this InDel might be used in future marker-assisted selection of drought-tolerant lines in soybean breeding programs.

Soybean ZINC FINGER PROTEIN03 targets two SUPEROXIDE DISMUTASE1s and confers resistance to Phytophthora sojae

Phytophthora sojae causes Phytophthora root and stem rot disease of soybean (Glycine max), leading to huge annual yield loss worldwide, but resistance to Phytophthora sojae (Rps) genes remains elusive. Soybean cultivar "Yudou 29" is resistant to P. sojae strain PsMC1, and this study aimed to clone, identify, and characterize the Rps gene in Yudou 29 (RpsYD29) and clarify its functional mechanism. We map-based cloned RpsYD29 (ZINC FINGER PROTEIN03, GmZFP03) using the families of a cross between Yudou 29 and a P. sojae-susceptible soybean cultivar "Jikedou 2". P. sojae resistance of GmZFP03 was functionally validated by stable soybean genetic transformation and allele-phenotype association analysis. GmZFP03 was identified as a C2H2-type zinc finger protein transcription factor, showing 4 amino acid residue polymorphisms (V79F, G122-, G123-, and D125V) and remarkably different expression patterns between resistant and susceptible soybeans. Notably boosted activity and gene expression of superoxide dismutase (SOD) in resistant-type GmZFP03-expressed transgenic soybean, substantial enhancement of P. sojae resistance of wild-type soybean by exogenous SOD treatment, and GmZFP03 binding to and activation of 2 SOD1 (Glyma.03g242900 and Glyma.19g240400) promoters demonstrated the involvement of SOD1s in GmZFP03-mediated resistance to P. sojae strain PsMC1. Thus, this study cloned the soybean P. sojae-resistant GmZFP03, the product of which specifically targets 2 SOD1 promoters. GmZFP03 can be directly used for precise P. sojae-resistance soybean breeding.

BSR and Full-Length Transcriptome Approaches Identified Candidate Genes for High Seed Ratio in Camellia vietnamensis

(1) Background: C. vietnamensis is very suitable for growth in the low hilly areas of southern subtropical regions. Under appropriate conditions, the oil yield of C. vietnamensis can reach 1125 kg/ha (the existing varieties can reach 750 kg/ha). Moreover, the fruit of C. vietnamensis is large and the pericarp is thick (>5 cm). Therefore, a high seed ratio has become the main target economic trait for the breeding of C. vietnamensis. (2) Methods: A half-sibling population of C. vietnamensis plants with a combination of high and low seed ratios was constructed by crossing a C. vietnamensis female parent. Bulked segregant RNA analysis and full-length transcriptome sequencing were performed to determine the molecular mechanisms underlying a high seed ratio. (3) Results: Seed ratio is a complex quantitative trait with a normal distribution, which is significantly associated with four other traits of fruit (seed weight, seed number, fruit diameter, and pericarp thickness). Two candidate regions related to high seed ratio (HSR) were predicted. One spanned 140.8−148.4 Mb of chromosome 2 and was associated with 97 seed-yield-related candidate genes ranging in length from 278 to 16,628 bp. The other spanned 35.3−37.3 Mb on chromosome 15 and was associated with 38 genes ranging in length from 221 to 16,928 bp. Using the full-length transcript as a template, a total of 115 candidate transcripts were obtained, and 78 transcripts were predicted to be functionally annotated. The DEGs from two set pairs of cDNA sequencing bulks were enriched to cytochrome p450 CYP76F14 (KOG0156; GO:0055114, HSR4, HSR7), the gibberellin phytohormone pathway (GO:0016787, HSR5), the calcium signaling pathway (GO:0005509, HSR6), the polyubiquitin-PPAR signaling pathway (GO:0005515, HSR2, HSR3), and several main transcription factors (bZIP transcription factor, HSR1) in C. vietnamensis.

Progress and prospectus in genetics and genomics of Phytophthora root and stem rot resistance in soybean (Glycine max L.)

Soybean is one of the largest sources of protein and oil in the world and is also considered a "super crop" due to several industrial advantages. However, enhanced acreage and adoption of monoculture practices rendered the crop vulnerable to several diseases. Phytophthora root and stem rot (PRSR) caused by Phytophthora sojae is one of the most prevalent diseases adversely affecting soybean production globally. Deployment of genetic resistance is the most sustainable approach for avoiding yield losses due to this disease. PRSR resistance is complex in nature and difficult to address by conventional breeding alone. Genetic mapping through a cost-effective sequencing platform facilitates identification of candidate genes and associated molecular markers for genetic improvement against PRSR. Furthermore, with the help of novel genomic approaches, identification and functional characterization of Rps (resistance to Phytophthora sojae) have also progressed in the recent past, and more than 30 Rps genes imparting complete resistance to different PRSR pathotypes have been reported. In addition, many genomic regions imparting partial resistance have also been identified. Furthermore, the adoption of emerging approaches like genome editing, genomic-assisted breeding, and genomic selection can assist in the functional characterization of novel genes and their rapid introgression for PRSR resistance. Hence, in the near future, soybean growers will likely witness an increase in production by adopting PRSR-resistant cultivars. This review highlights the progress made in deciphering the genetic architecture of PRSR resistance, genomic advances, and future perspectives for the deployment of PRSR resistance in soybean for the sustainable management of PRSR disease.

Molecular Breeding and Drought Tolerance in Chickpea

Cicer arietinum L. is the third greatest widely planted imperative pulse crop worldwide, and it belongs to the Leguminosae family. Drought is the utmost common abiotic factor on plants, distressing their water status and limiting their growth and development. Chickpea genotypes have the natural ability to fight drought stress using certain strategies viz., escape, avoidance and tolerance. Assorted breeding methods, including hybridization, mutation, and marker-aided breeding, genome sequencing along with omics approaches, could be used to improve the chickpea germplasm lines(s) against drought stress. Root features, for instance depth and root biomass, have been recognized as the greatest beneficial morphological factors for managing terminal drought tolerance in the chickpea. Marker-aided selection, for example, is a genomics-assisted breeding (GAB) strategy that can considerably increase crop breeding accuracy and competence. These breeding technologies, notably marker-assisted breeding, omics, and plant physiology knowledge, underlined the importance of chickpea breeding and can be used in future crop improvement programmes to generate drought-tolerant cultivars(s).

Breeding for disease resistance in soybean: a global perspective

KEY MESSAGE

This review provides a comprehensive atlas of QTLs, genes, and alleles conferring resistance to 28 important diseases in all major soybean production regions in the world. Breeding disease-resistant soybean [Glycine max (L.) Merr.] varieties is a common goal for soybean breeding programs to ensure the sustainability and growth of soybean production worldwide. However, due to global climate change, soybean breeders are facing strong challenges to defeat diseases. Marker-assisted selection and genomic selection have been demonstrated to be successful methods in quickly integrating vertical resistance or horizontal resistance into improved soybean varieties, where vertical resistance refers to R genes and major effect QTLs, and horizontal resistance is a combination of major and minor effect genes or QTLs. This review summarized more than 800 resistant loci/alleles and their tightly linked markers for 28 soybean diseases worldwide, caused by nematodes, oomycetes, fungi, bacteria, and viruses. The major breakthroughs in the discovery of disease resistance gene atlas of soybean were also emphasized which include: (1) identification and characterization of vertical resistance genes reside rhg1 and Rhg4 for soybean cyst nematode, and exploration of the underlying regulation mechanisms through copy number variation and (2) map-based cloning and characterization of Rps11 conferring resistance to 80% isolates of Phytophthora sojae across the USA. In this review, we also highlight the validated QTLs in overlapping genomic regions from at least two studies and applied a consistent naming nomenclature for these QTLs. Our review provides a comprehensive summary of important resistant genes/QTLs and can be used as a toolbox for soybean improvement. Finally, the summarized genetic knowledge sheds light on future directions of accelerated soybean breeding and translational genomics studies.

Mapping Locus RSC11K and predicting candidate gene resistant to Soybean mosaic virus strain SC11 through linkage analysis combined with genome resequencing of the parents in soybean

Soybean mosaic virus (SMV) strain SC11 was prevalent in middle China. Its resistance was controlled by a Mendelian single dominant gene RSC11K in soybean Kefeng-1. This study aimed at mapping RSC11K and identifying its candidate gene. RSC11K locus was mapped ~217 kb interval between two SNP-linkage-disequilibrium-blocks (Gm02_BLOCK_11273955_11464884 and Gm02_BLOCK_11486875_11491354) in W82.a1.v1 genome using recombinant inbred lines population derived from Kefeng-1 (Resistant) × NN1138-2 (Susceptible), but inserted with a ~245 kb segment in W82.a2.v1 genome. In the entire 462 kb RSC11K region, 429 SNPs, 142 InDels and 34 putative genes were identified with more SNPs/InDels distributed in non-functional regions. Thereinto, ten genes contained SNP/InDel variants with high and moderate functional impacts on proteins, among which Glyma.02G119700 encoded a typical innate immune receptor-like kinase involving in virus disease process and responded to SMV inoculation, therefore was recognized as RSC11K's candidate gene. The novel RSC11K locus and candidate genes may help developing SMV resistance germplasm.

Identification and molecular mapping of Rps14, a gene conferring broad-spectrum resistance to Phytophthora sojae in soybean

KEY MESSAGE

A soybean landrace carries broad-spectrum resistance to Phytophthora sojae, which is conferred by a single gene, designated Rps14, on the short arm of chromosome 3. Phytophthora sojae is the causative agent for Phytophthora root and stem rot in soybean [Glycine max (L.) Merr.] and can be managed by deployment of resistance to P. sojae (Rps) genes. PI 340,029 is a soybean landrace carrying broad-spectrum resistance to the pathogen. Analysis of an F2 population derived from a cross between PI 340,029 and a susceptible cultivar 'Williams' reveals that the resistance to P. sojae race 1 is conferred by a single gene, designated Rps14, which was initially mapped to a 4.5-cM region on the short arm of chromosome 3 by bulked segregant analysis (BSA), and subsequently narrowed to a 1.48 cM region corresponding to 229-kb in the Williams 82 reference genome (Wm82 v2.a1), using F3:4 families derived from the F2 population. Further analysis indicates that the broad-spectrum resistance carried by PI 340,029 is fully attributable to Rps14. The genomic sequences corresponding to the defined Rps14 region from a set of diverse soybean varieties exhibit drastic NBS-LRR gene copy number variation, ranging from 3 to 17 copies. Ultimate isolation of Rps14 would be critical for precise selection and deployment of the gene for soybean protection.

Genetic mapping high protein content QTL from soybean 'Nanxiadou 25' and candidate gene analysis

BACKGROUND

Soybean is a globally important legume crop that provides a primary source of high-quality vegetable protein and oil. Seed protein content (SPC) is a valuable quality trait controlled by multiple genes in soybean.

RESULTS

In this study, we performed quantitative trait loci (QTL) mapping, QTL-seq, and RNA sequencing (RNA-seq) to reveal the genes controlling protein content in the soybean by using the high protein content variety Nanxiadou 25. A total of 50 QTL for SPC distributed on 14 chromosomes except chromosomes 4, 12, 14, 17, 18, and 19 were identified by QTL mapping using 178 recombinant inbred lines (RILs). Among these QTL, the major QTL qSPC_20-1 and qSPC_20-2 on chromosome 20 were repeatedly detected across six tested environments, corresponding to the location of the major QTL detected using whole-genome sequencing-based QTL-seq. 329 candidate DEGs were obtained within the QTL region of qSPC_20-1 and qSPC_20-2 via gene expression profile analysis. Nine of which were associated with SPC, potentially representing candidate genes. Clone sequencing results showed that different single nucleotide polymorphisms (SNPs) and indels between high and low protein genotypes in Glyma.20G088000 and Glyma.16G066600 may be the cause of changes in this trait.

CONCLUSIONS

These results provide the basis for research on candidate genes and marker-assisted selection (MAS) in soybean breeding for seed protein content.

Identification of three closely linked loci conferring broad-spectrum Phytophthora sojae resistance in soybean variety Tosan-231

KEY MESSAGE

A variable genomic region containing two Harosoy-derived loci related to Rps7 and one Nemashirazu-derived locus confers broad-spectrum Phytophthora sojae resistance in Tosan-231 and is useful for developing resistant cultivars. We investigated resistance to pathotypically variable Phytophthora sojae isolates in the soybean variety Tosan-231, which has broad-spectrum resistance. Mapping analysis using descendent lines from a cross between Shuurei and Tosan-231 demonstrated that a genomic region between SSR markers BARCSOYSSR_03_0209 and BARCSOYSSR_03_0385 (termed "Region T"), confers broad-spectrum resistance in Tosan-231 and contains three closely linked resistance loci. Inoculation tests with 20 P. sojae isolates of different pathotypes and simple sequence repeat (SSR) analysis of progenitors of Tosan-231 facilitated identification and characterization of Rps genes at the three resistance loci. Two resistance genes, RpsT1 and RpsT2, were found to be derived from Harosoy carrying Rps7. This result suggested two mutually exclusive possibilities: (1) either RpsT1 or RpsT2 is Rps7, and the other is a locally functional novel gene; (2) Rps7 is not a single gene but in fact comprises RpsT1 and RpsT2. The resistance locus containing RpsT3 is derived from Nemashirazu, in which Rps genes have remained poorly defined. Moreover, we identified two genomic regions with relatively high recombination frequencies on the basis of mapping information and proposed a strategy to readily assemble useful resistance genes in or around Region T.

Identification of alkali-tolerant candidate genes using the NGS-assisted BSA strategy in rice

Rice (Oryza sativa L.) is a saline-alkali-sensitive crop. Saline-alkali environments can seriously affect the growth, development, and yield of rice. The mechanisms of salt tolerance and alkali tolerance in rice are different; thus, it is very important to study and explore the alkali-tolerant gene loci to improve the saline-alkali tolerance of rice varieties. In this study, the japonica rice varieties Dongnong 425 (DN425) and Changbai 10 (CB10) and a hybridized recombinant inbred line (RIL) population were used as materials to be irrigated with Na2CO3 solution under field test conditions. A resistant pool (R-pool) and a sensitive pool (S-pool) were constructed by selecting the lines with extremely high and extremely low 1000-grain weight (TGW), respectively, from the RIL population under alkali treatment. Four candidate TGW regions on chromosomes (Chr.) 2 and 3 were associated using the bulked segregant analysis (BSA) strategy assisted by next-generation sequencing (NGS) technology (NGS-assisted BSA). Using the linkage analysis, QTL-qATGW2-2 in the candidate region was mapped within a range of 116 Kb between the SSR marker RM13592 and the Indel marker Indel3 of Chr. 2, which contained 18 predictive genes. The BSA sequencing results showed that Os02g39884 contained a nonsynonymous substitution mutation SNP (nsSNP), leading to the transformation of a residue from arginine (cGg) to glutamine (cAg); thus, Os02g39884 was inferred to be the candidate gene of qATGW2-2. The results of the qRT-PCR analysis also confirmed this. This paper provides important information for the rapid and accurate identification of the alkali-tolerant gene loci in rice.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01228-x.

Mining germplasm panels and phenotypic datasets to identify loci for resistance to Phytophthora sojae in soybean

Phytophthora sojae causes Phytophthora root and stem rot of soybean and has been primarily managed through deployment of qualitative Resistance to P. sojae genes (Rps genes). The effectiveness of each individual or combination of Rps gene(s) depends on the diversity and pathotypes of the P. sojae populations present. Due to the complex nature of P. sojae populations, identification of more novel Rps genes is needed. In this study, phenotypic data from previous studies of 16 panels of plant introductions (PIs) were analyzed. Panels 1 and 2 consisted of 448 Glycine max and 520 G. soja, which had been evaluated for Rps gene response with a combination of P. sojae isolates. Panels 3 and 4 consisted of 429 and 460 G. max PIs, respectively, which had been evaluated using individual P. sojae isolates with complex virulence pathotypes. Finally, Panels 5-16 (376 G. max PIs) consisted of data deposited in the USDA Soybean Germplasm Collection from evaluations with 12 races of P. sojae. Using these panels, genome-wide association (GWA) analyses were carried out by combining phenotypic and SoySNP50K genotypic data. GWA models identified two, two, six, and seven novel Rps loci with Panels 1, 2, 3, and 4, respectively. A total of 58 novel Rps loci were identified using Panels 5-16. Genetic and phenotypic dissection of these loci may lead to the characterization of novel Rps genes that can be effectively deployed in new soybean cultivars against diverse P. sojae populations.

Cytological Observations and Bulked-Segregant Analysis Coupled Global Genome Sequencing Reveal Two Genes Associated with Pollen Fertility in Tetraploid Rice

Neo-tetraploid rice with high fertility is a useful germplasm for polyploid rice breeding, which was developed from the crossing of different autotetraploid rice lines. However, little information is available on the molecular mechanism underlying the fertility of neo-tetraploid rice. Here, two contrasting populations of tetraploid rice, including one with high fertility (hereafter referred to as JG) and another with low fertility (hereafter referred to as JD), were generated by crossing Huaduo 3 (H3), a high fertility neo-tetraploid rice that was developed by crossing Jackson-4x with 96025-4x, and Huajingxian74-4x (T452), a low fertility autotetraploid rice parent. Cytological, global genome sequencing-based bulked-segregant (BSA-seq) and CRISPR/Cas9 technology were employed to study the genes associated with pollen fertility in neo-tetraploid rice. The embryo sacs of JG and JD lines were normal; however, pollen fertility was low in JD, which led to scarce fertilization and low seed setting. Cytological observations displayed low pollen fertility (25.1%) and approximately 31.3 and 27.2% chromosome lagging at metaphase I and II, and 28.8 and 24.8% chromosome straggling at anaphase I and II in JD, respectively. BSA-seq of F_2–3 generations and RNA-seq of F₄ generation detected a common fragment, i.e., 18,915,234–19,500,000, at chromosome 7, which was comprised of 78 genes associated with fertility. Among 78 genes, 9 genes had been known to be involved in meiosis and pollen development. Two mutants ny1 (LOC_Os07g32406) and ny2 (LOC_Os07g32040) were generated by CRISPR/Cas9 knockout in neo-tetraploid rice, and which exhibited low pollen fertility and abnormal chromosome behavior. Our study revealed that two unknown genes, LOC_Os07g32406 (NY1) and LOC_Os07g32040 (NY2) play an important role in pollen development of neo-tetraploid rice and provides a new perspective about the genetic mechanisms of fertility in polyploid rice.

Genetic Mapping of a Resistance Locus to Phytophthora sojae in the Korean Soybean Cultivar Daewon

Phytophthora root and stem rot reduce soybean yields worldwide. The use of R-gene type resistance is currently crucial for protecting soybean production. The present study aimed to identify the genomic location of a gene conferring resistance to Phytophthora sojae isolate 2457 in the recombinant inbred line population developed by a cross of Daepung × Daewon. Single-marker analysis identified 20 single nucleotide polymorphisms associated with resistance to the P. sojae isolate 2457, which explained ~67% of phenotypic variance. Daewon contributed a resistance allele for the locus. This region is a well-known location for Rps1 and Rps7. The present study is the first, however, to identify an Rps gene locus from a major soybean variety cultivated in South Korea. Linkage analysis also identified a 573 kb region on chromosome 3 with high significance (logarithm of odds = 13.7). This genomic region was not further narrowed down due to lack of recombinants within the interval. Based on the latest soybean genome, ten leucine-rich repeat coding genes and four serine/threonine protein kinase-coding genes are annotated in this region, which all are well-known types of genes for conferring disease resistance in crops. These genes would be candidates for molecular characterization of the resistance in further studies. The identified R-gene locus would be useful in developing P. sojae resistant varieties in the future. The results of the present study provide foundational knowledge for researchers who are interested in soybean-P. sojae interaction.

Genomic interventions for sustainable agriculture

Agricultural production faces a Herculean challenge to feed the increasing global population. Food production systems need to deliver more with finite land and water resources while exerting the least negative influence on the ecosystem. The unpredictability of climate change and consequent changes in pests/pathogens dynamics aggravate the enormity of the challenge. Crop improvement has made significant contributions towards food security, and breeding climate-smart cultivars are considered the most sustainable way to accelerate food production. However, a fundamental change is needed in the conventional breeding framework in order to respond adequately to the growing food demands. Progress in genomics has provided new concepts and tools that hold promise to make plant breeding procedures more precise and efficient. For instance, reference genome assemblies in combination with germplasm sequencing delineate breeding targets that could contribute to securing future food supply. In this review, we highlight key breakthroughs in plant genome sequencing and explain how the presence of these genome resources in combination with gene editing techniques has revolutionized the procedures of trait discovery and manipulation. Adoption of new approaches such as speed breeding, genomic selection and haplotype-based breeding could overcome several limitations of conventional breeding. We advocate that strengthening varietal release and seed distribution systems will play a more determining role in delivering genetic gains at farmer's field. A holistic approach outlined here would be crucial to deliver steady stream of climate-smart crop cultivars for sustainable agriculture.

Identification of a genomic region controlling thermotolerance at flowering in maize using a combination of whole genomic re-sequencing and bulked segregant analysis

A novel genomic region controlling thermotolerance at flowering was identified by the combination of whole genomic re-sequencing and bulked segregant analysis in maize. The increasing frequency of extreme high temperature has brought a great threat to the development of maize throughout its life cycle, especially during the flowering phase. However, the genetic basis of thermotolerance at flowering in maize remains poorly understood. Here, we characterized a thermotolerant maize ecotype Abe2 and dissected its genetic basis using a F_2:8 recombinant inbred line (RIL) population generated from a cross between Abe2 and B73. After continuous high temperature stress above 35 °C for 17 days, Abe2 and B73 show distinct leaf scorching phenotype under field conditions. To identify the genomic regions associated with the phenotypic variation, we applied a combination of whole genomic re-sequencing and bulked segregant analysis, and revealed 10,316,744 SNPs and 1,488,302 InDels between the two parental lines, and 2,693,054 SNPs and 313,757 InDels between the two DNA pools generated from the thermos-tolerant and the sensitive individuals of the RIL, of which, 108,655 and 17,853 SNPs may cause nonsynonymous variations. Finally, a 7.41 Mb genomic region on chromosome 1 was identified, and 7 candidate genes were annotated to participate in high temperature-related stress response. A candidate gene Zm00001d033339 encoding a serine/threonine protein kinase was proposed to be the most likely causative gene contributing to the thermotolerance at flowering by involving in stomatal movement (GO: 0010119) via Abscisic acid (ABA) pathway (KO04075). This work could provide an opportunity for gene cloning and pyramiding breeding to improve thermotolerance at flowering in maize.

Fine Mapping, Candidate Gene Identification and Co-segregating Marker Development for the Phytophthora Root Rot Resistance Gene RpsYD25

Phytophthora root rot (PRR) caused by Phytophthora sojae is a serious disease of soybean. The most effective disease-control strategy is to deploy resistant cultivars carrying Rps genes. Soybean cultivar Yudou25 can effectively resist pathotypes of P. sojae in China. Previous studies have mapped the Rps gene in Yudou25, RpsYD25, on chromosome 3. In this study, at first RpsYD25 was located between SSR markers Satt1k3 (2.2 cM) and BARCSOYSSR_03_0253 (4.5 cM) by using an F_2:3 population containing 165 families derived from Zaoshu18 and Yudou25. Then the recombination sites were identified in 1127 F_3:4 families derived from Zaoshu18 and Yudou25 using the developed PCR-based SNP, InDel and SSR markers, and RpsYD25 was finely mapped in the a 101.3 kb genomic region. In this region, a zinc ion binding and nucleic acid binding gene Glyma.03g034700 and two NBS-LRR genes Glyma.03g034800 and Glyma.03g034900 were predicted as candidate genes of RpsYD25, and five co-segregated SSR markers with RpsYD25 were identified and validated to be diagnostic markers. Combined with the resistance reaction to multiple P. sojae isolates, seven of 178 soybean genotypes were detected to contain RpsYD25 using the five co-segregated SSR markers. The soybean genotypes carrying RpsYD25 and the developed co-segregated markers can be effectively applied in the breeding for P. sojae resistance in China.

Fine mapping and gene cloning in the post-NGS era: advances and prospects

Improvement in traits of agronomic importance is the top breeding priority of crop improvement programs. Majority of these agronomic traits show complex quantitative inheritance. Identification of quantitative trait loci (QTLs) followed by fine mapping QTLs and cloning of candidate genes/QTLs is central to trait analysis. Advances in genomic technologies revolutionized our understanding of genetics of complex traits, and genomic regions associated with traits were employed in marker-assisted breeding or cloning of QTLs/genes. Next-generation sequencing (NGS) technologies have enabled genome-wide methodologies for the development of ultra-high-density genetic linkage maps in different crops, thus allowing placement of candidate loci within few kbs in genomes. In this review, we compare the marker systems used for fine mapping and QTL cloning in the pre- and post-NGS era. We then discuss how different NGS platforms in combination with advanced experimental designs have improved trait analysis and fine mapping. We opine that efficient genotyping/sequencing assays may circumvent the need for cumbersome procedures that were earlier used for fine mapping. A deeper understanding of the trait architectures of agricultural significance will be crucial to accelerate crop improvement.

Impacts of genomic research on soybean improvement in East Asia

It has been commonly accepted that soybean domestication originated in East Asia. Although East Asia has the historical merit in soybean production, the USA has become the top soybean producer in the world since 1950s. Following that, Brazil and Argentina have been the major soybean producers since 1970s and 1990s, respectively. China has once been the exporter of soybean to Japan before 1990s, yet she became a net soybean importer as Japan and the Republic of Korea do. Furthermore, the soybean yield per unit area in East Asia has stagnated during the past decade. To improve soybean production and enhance food security in these East Asian countries, much investment has been made, especially in the breeding of better performing soybean germplasms. As a result, China, Japan, and the Republic of Korea have become three important centers for soybean genomic research. With new technologies, the rate and precision of the identification of important genomic loci associated with desired traits from germplasm collections or mutants have increased significantly. Genome editing on soybean is also becoming more established. The year 2019 marked a new era for crop genome editing in the commercialization of the first genome-edited plant product, which is a high-oleic-acid soybean oil. In this review, we have summarized the latest developments in soybean breeding technologies and the remarkable progress in soybean breeding-related research in China, Japan, and the Republic of Korea.

Fine mapping of a Phytophthora-resistance locus RpsGZ in soybean using genotyping-by-sequencing

BACKGROUND

Phytophthora root rot (PRR) caused by Phytophthora sojae (P. sojae) is one of the most serious limitations to soybean production worldwide. The identification of resistance gene(s) and their incorporation into elite varieties is an effective approach for breeding to prevent soybean from being harmed by this disease. A valuable mapping population of 228 F8:11 recombinant inbred lines (RILs) derived from a cross of the resistant cultivar Guizao1 and the susceptible cultivar BRSMG68 and a high-density genetic linkage map with an average distance of 0.81 centimorgans (cM) between adjacent bin markers in this population were used to map and explore candidate gene(s).

RESULTS

PRR resistance in Guizao1 was found to be controlled by a single Mendelian locus and was finely mapped to a 367.371-kb genomic region on chromosome 3 harbouring 19 genes, including 7 disease resistance (R)-like genes, in the reference Willliams 82 genome. Quantitative real-time PCR assays of possible candidate genes revealed that Glyma.03 g05300 was likely involved in PRR resistance.

CONCLUSIONS

These findings from the fine mapping of a novel Rps locus will serve as a basis for the cloning and transfer of resistance genes in soybean and the breeding of P. sojae-resistant soybean cultivars through marker-assisted selection.

Discovery of two novel and adjacent QTLs on chromosome B02 controlling resistance against bacterial wilt in peanut variety Zhonghua 6

Two novel and adjacent genomics and candidate genes for bacterial wilt resistance were identified on chromosome B02 in peanut variety Zhonghua 6 using both traditional QTL mapping and QTL-seq methods. Peanut (Arachis hypogaea) is an important oilseed crop worldwide. Utilization of genetic resistance is the most economic and effective approach to control bacterial wilt, one of the most devastating plant diseases, in peanut production. To accelerate the genetic improvement of bacterial wilt resistance (BWR) in peanut breeding programs, quantitative trait locus (QTL) mapping has been conducted for two resistant varieties. In this context, we deployed linkage mapping as well as sequencing-based mapping approach, QTL-seq, to identify genomic regions and candidate genes for BWR in another highly resistant variety Zhonghua 6. The recombination inbred line population (268 progenies) from the cross Xuhua 13 × Zhonghua 6 was used in BWR evaluation across five environments. QTL mapping using both SSR- and SNP-based genetic maps identified a stable QTL (qBWRB02-1) on chromosome B02 with 37.79-78.86% phenotypic variation explained (PVE) across five environments. The QTL-seq facilitated further dissection of qBWRB02-1 into two adjacent genomic regions, qBWRB02-1-1 (2.81-4.24 Mb) and qBWRB02-1-2 (6.54-8.75 Mb). Mapping of newly developed Kompetitive allele-specific PCR (KASP) markers on the genetic map confirmed their stable expressions across five environments. The effects of qBWRB02-1-1 (49.43-68.86% PVE) were much higher than qBWRB02-1-2 (3.96-6.48% PVE) and other previously reported QTLs. Nineteen putative candidate genes affected by 49 non-synonymous SNPs were identified for qBWRB02-1-1, and ten of them were predicted to code for disease resistance proteins. The major and stable QTL qBWRB02-1-1 and validated KASP markers could be deployed in genomics-assisted breeding (GAB) to develop improved peanut varieties with enhanced BWR.

PyBSASeq: a simple and effective algorithm for bulked segregant analysis with whole-genome sequencing data

BACKGROUND

Bulked segregant analysis (BSA), coupled with next-generation sequencing, allows the rapid identification of both qualitative and quantitative trait loci (QTL), and this technique is referred to as BSA-Seq here. The current SNP index method and G-statistic method for BSA-Seq data analysis require relatively high sequencing coverage to detect significant single nucleotide polymorphism (SNP)-trait associations, which leads to high sequencing cost.

RESULTS

We developed a simple and effective algorithm for BSA-Seq data analysis and implemented it in Python; the program was named PyBSASeq. Using PyBSASeq, the significant SNPs (sSNPs), SNPs likely associated with the trait, were identified via Fisher's exact test, and then the ratio of the sSNPs to total SNPs in a chromosomal interval was used to detect the genomic regions that condition the trait of interest. The results obtained this way are similar to those generated via the current methods, but with more than five times higher sensitivity. This approach was termed the significant SNP method here.

CONCLUSIONS

The significant SNP method allows the detection of SNP-trait associations at much lower sequencing coverage than the current methods, leading to ~ 80% lower sequencing cost and making BSA-Seq more accessible to the research community and more applicable to the species with a large genome.

GmWRKY40, a member of the WRKY transcription factor genes identified from Glycine max L., enhanced the resistance to Phytophthora sojae

BACKGROUND

The WRKY proteins are a superfamily of transcription factors and members play essential roles in the modulation of diverse physiological processes, such as growth, development, senescence and response to biotic and abiotic stresses. However, the biological roles of the majority of the WRKY family members remains poorly understood in soybean relative to the research progress in model plants.

RESULTS

In this study, we identified and characterized GmWRKY40, which is a group IIc WRKY gene. Transient expression analysis revealed that the GmWRKY40 protein is located in the nucleus of plant cells. Expression of GmWRKY40 was strongly induced in soybean following infection with Phytophthora sojae, or treatment with methyl jasmonate, ethylene, salicylic acid, and abscisic acid. Furthermore, soybean hairy roots silencing GmWRKY40 enhanced susceptibility to P. sojae infection compared with empty vector transgenic roots. Moreover, suppression of GmWRKY40 decreased the accumulation of reactive oxygen species (ROS) and modified the expression of several oxidation-related genes. Yeast two-hybrid experiment combined with RNA-seq analysis showed that GmWRKY40 interacted with 8 JAZ proteins with or without the WRKY domain or zinc-finger domain of GmWRKY40, suggesting there were different interaction patterns among these interacted proteins.

CONCLUSIONS

Collectively, these results suggests that GmWRKY40 functions as a positive regulator in soybean plants response to P. sojae through modulating hydrogen peroxide accumulation and JA signaling pathway.

Dissecting Resistance to Phytophthora cinnamomi in Interspecific Hybrid Chestnut Crosses Using Sequence-Based Genotyping and QTL Mapping

The soilborne oomycete Phytophthora cinnamomi-which causes root rot, trunk cankers, and stem lesions on an estimated 5,000 plant species worldwide-is a lethal pathogen of American chestnut (Castanea dentata) as well as many other woody plant species. P. cinnamomi is particularly damaging to chestnut and chinquapin trees (Castanea spp.) in the southern portion of its native range in the United States due to relatively mild climatic conditions that are conductive to disease development. Introduction of resistant genotypes is the most practical solution for disease management in forests because treatment with fungicides and eradication of the pathogen are neither practical nor economically feasible in natural ecosystems. Using backcross families derived from crosses of American chestnuts with two resistant Chinese chestnut cultivars Mahogany and Nanking, we constructed linkage maps and identified quantitative trait loci (QTLs) for resistance to P. cinnamomi that had been introgressed from these Chinese chestnut cultivars. In total, 957 plants representing five cohorts of three hybrid crosses were genotyped by sequencing and phenotyped by standardized inoculation and visual examination over a 6-year period from 2011 to 2016. Eight parental linkage maps comprising 7,715 markers were constructed, and 17 QTLs were identified on four linkage groups (LGs): LG_A, LG_C, LG_E, and LG_K. The most consistent QTLs were detected on LG_E in seedlings from crosses with both 'Mahogany' and 'Nanking' and LG_K in seedlings from 'Mahogany' crosses. Two consistent large and medium effect QTLs located ∼10 cM apart were present in the middle and at the lower end of LG_E; other QTLs were considered to have small effects. These results imply that the genetic architecture of resistance to P. cinnamomi in Chinese chestnut × American chestnut hybrid progeny may resemble the P. sojae-soybean pathosystem, with a few dominant QTLs along with quantitatively inherited partial resistance conferred by multiple small-effect QTLs.

Discovery of genomic regions and candidate genes controlling shelling percentage using QTL-seq approach in cultivated peanut (Arachis hypogaea L.)

Cultivated peanut (Arachis hypogaea L.) is an important grain legume providing high-quality cooking oil, rich proteins and other nutrients. Shelling percentage (SP) is the 2nd most important agronomic trait after pod yield and this trait significantly affects the economic value of peanut in the market. Deployment of diagnostic markers through genomics-assisted breeding (GAB) can accelerate the process of developing improved varieties with enhanced SP. In this context, we deployed the QTL-seq approach to identify genomic regions and candidate genes controlling SP in a recombinant inbred line population (Yuanza 9102 × Xuzhou 68-4). Four libraries (two parents and two extreme bulks) were constructed and sequenced, generating 456.89-790.32 million reads and achieving 91.85%-93.18% genome coverage and 14.04-21.37 mean read depth. Comprehensive analysis of two sets of data (Yuanza 9102/two bulks and Xuzhou 68-4/two bulks) using the QTL-seq pipeline resulted in discovery of two overlapped genomic regions (2.75 Mb on A09 and 1.1 Mb on B02). Nine candidate genes affected by 10 SNPs with non-synonymous effects or in UTRs were identified in these regions for SP. Cost-effective KASP (Kompetitive Allele-Specific PCR) markers were developed for one SNP from A09 and three SNPs from B02 chromosome. Genotyping of the mapping population with these newly developed KASP markers confirmed the major control and stable expressions of these genomic regions across five environments. The identified candidate genomic regions and genes for SP further provide opportunity for gene cloning and deployment of diagnostic markers in molecular breeding for achieving high SP in improved varieties.

Characterization and Rapid Gene-Mapping of Leaf Lesion Mimic Phenotype of spl-1 Mutant in Soybean (Glycine max (L.) Merr.)

In plants, lesion mimic mutants (LMMs) reveal spontaneous disease-like lesions in the absence of pathogen that constitutes powerful genetic material to unravel genes underlying programmed cell death (PCD), particularly the hypersensitive response (HR). However, only a few LMMs are reported in soybean, and no related gene has been cloned until now. In the present study, we isolated a new LMM named spotted leaf-1 (spl-1) from NN1138-2 cultivar through ethyl methanesulfonate (EMS) treatment. The present study revealed that lesion formation might result from PCD and excessive reactive oxygen species (ROS) accumulation. The chlorophyll content was significantly reduced but antioxidant activities, viz., superoxide dismutase (SOD), peroxidase (POD) and catalase (CAT), as well as the malondialdehyde (MDA) contents, were detected higher in spl-1 than in the wild-type. According to segregation analysis of mutant phenotype in two genetic populations, viz., W82×spl-1 and PI378692×spl-1, the spotted leaf phenotype of spl-1 is controlled by a single recessive gene named lm1. The lm1 locus governing mutant phenotype of spl-1 was first identified in 3.15 Mb genomic region on chromosome 04 through MutMap analysis, which was further verified and fine mapped by simple sequence repeat (SSR) marker-based genetic mapping. Genetic linkage analysis narrowed the genomic region (lm1 locus) for mutant phenotype to a physical distance of ~76.23 kb. By searching against the Phytozome database, eight annotated candidate genes were found within the lm1 region. qRT-PCR expression analysis revealed that, among these eight genes, only Glyma.04g242300 showed highly significant expression levels in wild-type relative to the spl-1 mutant. However, sequencing data of the CDS region showed no nucleotide difference between spl-1 and its wild type within the coding regions of these genes but might be in the non-coding regions such as 5' or 3' UTR. Hence, the data of the present study are in favor of Glyma.04g242300 being the possible candidate genes regulating the mutant phenotype of spl-1. However, further validation is needed to prove this function of the gene as well as its role in PCD, which in turn would be helpful to understand the mechanism and pathways involved in HR disease resistance of soybean.

Genetic Mapping and Molecular Characterization of a Broad-spectrum Phytophthora sojae Resistance Gene in Chinese Soybean

Phytophthora root rot (PRR) causes serious annual soybean yield losses worldwide. The most effective method to prevent PRR involves growing cultivars that possess genes conferring resistance to Phytophthora sojae (Rps). In this study, QTL-sequencing combined with genetic mapping was used to identify RpsX in soybean cultivar Xiu94-11 resistance to all P. sojae isolates tested, exhibiting broad-spectrum PRR resistance. Subsequent analysis revealed RpsX was located in the 242-kb genomic region spanning the RpsQ locus. However, a phylogenetic investigation indicated Xiu94-11 carrying RpsX is distantly related to the cultivars containing RpsQ, implying RpsX and RpsQ have different origins. An examination of candidate genes revealed RpsX and RpsQ share common nonsynonymous SNP and a 144-bp insertion in the Glyma.03g027200 sequence encoding a leucine-rich repeat (LRR) region. Glyma.03g027200 was considered to be the likely candidate gene of RpsQ and RpsX. Sequence analyses confirmed that the 144-bp insertion caused by an unequal exchange resulted in two additional LRR-encoding fragments in the candidate gene. A marker developed based on the 144-bp insertion was used to analyze the genetic population and germplasm, and proved to be useful for identifying the RpsX and RpsQ alleles. This study implies that the number of LRR units in the LRR domain may be important for PRR resistance in soybean.

Complete resistance to powdery mildew and partial resistance to downy mildew in a Cucumis hystrix introgression line of cucumber were controlled by a co-localized locus

Key message A single recessive gene for complete resistance to powdery mildew and a major-effect QTL for partial resistance to downy mildew were co-localized in a Cucumis hystrix introgression line of cucumber. Downy mildew (DM) and powdery mildew (PM) are two major foliar diseases in cucumber. DM resistance (DMR) and PM resistance (PMR) may share common components; however, the genetic relationship between them remains unclear. IL52, a Cucumis hystrix introgression line of cucumber which has been reported to possess DMR, was recently identified to exhibit PMR as well. In this study, a single recessive gene pm for PMR was mapped to an approximately 468-kb region on chromosome 5 with 155 recombinant inbred lines (RILs) and 193 F₂ plants derived from the cross between a susceptible line 'changchunmici' and IL52. Interestingly, pm was co-localized with the major-effect DMR QTL dm5.2 confirmed by combining linkage analysis and BSA-seq, which was consistent with the observed linkage of DMR and PMR in IL52. Further, phenotype-genotype correlation analysis of DMR and PMR in the RILs indicated that the co-localized locus pm/dm5.2 confers complete resistance to PM and partial resistance to DM. Seven candidate genes for DMR were identified within dm5.2 by BSA-seq analysis, of which Csa5M622800.1, Csa5M622830.1 and Csa5M623490.1 were also the same candidate genes for PMR. A single nucleotide polymorphism that is present in the 3' untranslated region (3'UTR) of Csa5M622830.1 co-segregated perfectly with PMR. The GATA transcriptional factor gene Csa5M622830.1 may be a likely candidate gene for DMR and PMR. This study has provided a clear evidence for the relationship between DMR and PMR in IL52 and sheds new light on the potential value of IL52 for cucumber DMR and PMR breeding program.

Fine Mapping and Identification of a Novel Phytophthora Root Rot Resistance Locus RpsZS18 on Chromosome 2 in Soybean

Phytophthora root rot (PRR) caused by Phytophthora sojae is a major soybean disease that causes severe economic losses worldwide. Using soybean cultivars carrying a Rps resistance gene is the most effective strategy for controlling this disease. We previously detected a novel Phytophthora resistance gene, RpsZS18, on chromosome 2 of the soybean cultivar Zaoshu18. The aim of the present study was to identify and finely map RpsZS18. We used 232 F2:3 families generated from a cross between Zaoshu18 (resistant) and Williams (susceptible) as the mapping population. Simple sequence repeat (SSR) markers distributed on chromosome 2 were used to map RpsZS18. First, 12 SSR markers linked with RpsZS18 were identified by linkage analyses, including two newly developed SSR markers, ZCSSR33 and ZCSSR46, that flanked the gene at distances of 0.9 and 0.5 cM, respectively. Second, PCR-based InDel markers were developed based on sequence differences between the two parents and used to further narrow down the mapping region of RpsZS18 to 71.3 kb. Third, haplotype analyses were carried out in the RpsZS18 region using 14 soybean genotypes with whole-genome resequencing. We detected six genes with unique haplotype sequences in Zaoshu18. Finally, quantitative real-time PCR assays of the six genes revealed an EF-hand calcium-binding domain containing protein encoding gene (Glyma.02g245700), a pfkB carbohydrate kinase encoding gene (Glyma.02g245800), and a gene with no functional annotation (Glyma.02g246300), are putative candidate PRR resistance genes. This study provides useful information for breeding P. sojae-resistant soybean cultivars.