Identification of Potential Gene Regulatory Pathways Affecting the Ratio of Four-Seed Pod in Soybean

Application of bulk segregant RNA-Seq (BSR-Seq) and allele-specific primers to study soybean powdery mildew resistance

BACKGROUND

Powdery mildew (PM) is one of the important soybean diseases, and host resistance could practically contribute to soybean PM management. To date, only the Rmd locus on chromosome (Chr) 16 was identified through traditional QTL mapping and GWAS, and it remains unclear if the bulk segregant RNA-Seq (BSR-Seq) methodology is feasible to explore additional PM resistance that might exist in other varieties.

RESULTS

BSR-Seq was applied to contrast genotypes and gene expressions between the resistant bulk (R bulk) and the susceptible bulk (S bulk), as well as the parents. The ∆(SNP-index) and G' value identified several QTL and significant SNPs/Indels on Chr06, Chr15, and Chr16. Differentially expressed genes (DEGs) located within these QTL were identified using HISAT2 and Kallisto, and allele-specific primers (AS-primers) were designed to validate the accuracy of phenotypic prediction. While the AS-primers on Chr06 or Chr15 cannot distinguish the resistant and susceptible phenotypes, AS-primers on Chr16 exhibited 82% accuracy prediction with an additive effect, similar to the SSR marker Satt431.

CONCLUSIONS

Evaluation of additional AS-primers in the linkage disequilibrium (LD) block on Chr16 further confirmed the resistant locus, derived from the resistant parental variety 'Kaohsiung 11' ('KS11'), not only overlaps with the Rmd locus with unique up-regulated LRR genes (Glyma.16G213700 and Glyma.16G215100), but also harbors a down-regulated MLO gene (Glyma.16G145600). Accordingly, this study exemplified the feasibility of BSR-Seq in studying biotrophic disease resistance in soybean, and showed the genetic makeup of soybean variety 'KS11' comprising the Rmd locus and one MLO gene.

Analysis of Genetic Diversity in Adzuki Beans (Vigna angularis): Insights into Environmental Adaptation and Early Breeding Strategies for Yield Improvement

Adzuki beans are widely cultivated in East Asia and are one of the earliest domesticated crops. In order to gain a deeper understanding of the genetic diversity and domestication history of adzuki beans, we conducted Genotyping by Sequencing (GBS) analysis on 366 landraces originating from Korea, China, and Japan, resulting in 6586 single-nucleotide polymorphisms (SNPs). Population structure analysis divided these 366 landraces into three subpopulations. These three subpopulations exhibited distinctive distributions, suggesting that they underwent extended domestication processes in their respective regions of origin. Phenotypic variance analysis of the three subpopulations indicated that the Korean-domesticated subpopulation exhibited significantly higher 100-seed weights, the Japanese-domesticated subpopulation showed significantly higher numbers of grains per pod, and the Chinese-domesticated subpopulation displayed significantly higher numbers of pods per plant. We speculate that these differences in yield-related traits may be attributed to varying emphases placed by early breeders in these regions on the selection of traits related to yield. A large number of genes related to biotic/abiotic stress resistance and defense were found in most quantitative trait locus (QTL) for yield-related traits using genome-wide association studies (GWAS). Genomic sliding window analysis of Tajima's D and a genetic differentiation coefficient (Fst) revealed distinct domestication selection signatures and genotype variations on these QTLs within each subpopulation. These findings indicate that each subpopulation would have been subjected to varied biotic/abiotic stress events in different origins, of which these stress events have caused balancing selection differences in the QTL of each subpopulation. In these balancing selections, plants tend to select genotypes with strong resistance under biotic/abiotic stress, but reduce the frequency of high-yield genotypes to varying degrees. These biotic/abiotic stressors impact crop yield and may even lead to selection purging, resulting in the loss of several high-yielding genotypes among landraces. However, this also fuels the flow of crop germplasms. Overall, balancing selection appears to have a more significant impact on the three yield-related traits compared to breeder-driven domestication selection. These findings are crucial for understanding the impact of domestication selection history on landraces and yield-related traits, aiding in the improvement of adzuki bean varieties.

Integrated Bioinformatics Analyses of PIN1, CKX, and Yield-Related Genes Reveals the Molecular Mechanisms for the Difference of Seed Number Per Pod Between Soybean and Cowpea

There is limited advancement on seed number per pod (SNPP) in soybean breeding, resulting in low yield in China. To address this issue, we identified PIN1 and CKX gene families that regulate SNPP in Arabidopsis, analyzed the differences of auxin and cytokinin pathways, and constructed interaction networks on PIN1, CKX, and yield-related genes in soybean and cowpea. First, the relative expression level (REL) of PIN1 and the plasma membrane localization and phosphorylation levels of PIN1 protein were less in soybean than in cowpea, which make auxin transport efficiency lower in soybean, and its two interacted proteins might be involved in serine hydrolysis, so soybean has lower SNPP than cowpea. Then, the CKX gene family, along with its positive regulatory factor ROCK1, had higher REL and less miRNA regulation in soybean flowers than in cowpea ones. These lead to higher cytokinin degradation level, which further reduces the REL of PIN1 and decreases soybean SNPP. We found that VuACX4 had much higher REL than GmACX4, although the two genes essential in embryo development interact with the CKX gene family. Next, a tandem duplication experienced by legumes led to the differentiation of CKX3 into CKX3a and CKX3b, in which CKX3a is a key gene affecting ovule number. Finally, in the yield-related gene networks, three cowpea CBP genes had higher RELs than two soybean CBP genes, low RELs of three soybean-specific IPT genes might lead to a decrease in cytokinin synthesis, and some negative and positive SNPP regulation were found, respectively, in soybean and cowpea. These networks may explain the SNPP difference in the two crops. We deduced that ckx3a or ckx3a ckx6 ckx7 mutants, interfering CYP88A, and over-expressed DELLA increase SNPP in soybean. This study reveals the molecular mechanism for the SNPP difference in the two crops, and provides an important idea for increasing soybean yield.