Evolution and association analysis of GmCYP78A10 gene with seed size/weight and pod number in soybean

UDP-glucosyltransferase 71C4 controls the flux of phenylpropanoid metabolism to shape cotton seed development

Seeds play a crucial role in plant reproduction, making it essential to identify genes that affect seed development. In this study, we focused on UDP-glucosyltransferase 71C4 (UGT71C4) in cotton, a member of the glycosyltransferase family that shapes seed width and length, thereby influencing seed index and seed cotton yield. Overexpression of UGT71C4 results in seed enlargement owing to its glycosyltransferase activity on flavonoids, which redirects metabolic flux from lignin to flavonoid metabolism. This shift promotes cell proliferation in the ovule via accumulation of flavonoid glycosides, significantly enhancing seed cotton yield and increasing the seed index from 10.66 g to 11.91 g. By contrast, knockout of UGT71C4 leads to smaller seeds through activation of the lignin metabolism pathway and redirection of metabolic flux back to lignin synthesis. This redirection leads to increased ectopic lignin deposition in the ovule, inhibiting ovule growth and development, and alters yield components, increasing the lint percentage from 41.42% to 43.40% and reducing the seed index from 10.66 g to 8.60 g. Our research sheds new light on seed size development and reveals potential pathways for enhancing seed yield.

Genome-Wide Association Analysis of Yield-Related Traits and Candidate Genes in Vegetable Soybean

Owing to the rising demand for vegetable soybean products, there is an increasing need for high-yield soybean varieties. However, the complex correlation patterns among quantitative traits with genetic architecture pose a challenge for improving vegetable soybean through breeding. Herein, a genome-wide association study (GWAS) was applied to 6 yield-related traits in 188 vegetable soybean accessions. Using a BLINK model, a total of 116 single nucleotide polymorphisms (SNPs) were identified for plant height, pod length, pod number, pod thickness, pod width, and fresh pod weight. Furthermore, a total of 220 genes were found in the 200 kb upstream and downstream regions of significant SNPs, including 11 genes encoding functional proteins. Among them, four candidate genes, Glyma.13G109100, Glyma.03G183200, Glyma.09G102200, and Glyma.09G102300 were analyzed for significant haplotype variations and to be in LD block, which encode MYB-related transcription factor, auxin-responsive protein, F-box protein, and CYP450, respectively. The relative expression of candidate genes in V030 and V071 vegetable soybean (for the plant height, pod number, and fresh pod weight of V030 were lower than those of the V071 strains) was significantly different, and these genes could be involved in plant growth and development via various pathways. Altogether, we identified four candidate genes for pod yield and plant height from vegetable soybean germplasm. This study provides insights into the genomic basis for improving soybean and crucial genomic resources that can facilitate genome-assisted high-yielding vegetable soybean breeding.

Dissection of quantitative trait nucleotides and candidate genes associated with agronomic and yield-related traits under drought stress in rapeseed varieties: integration of genome-wide association study and transcriptomic analysis

Introduction

An important strategy to combat yield loss challenge is the development of varieties with increased tolerance to drought to maintain production. Improvement of crop yield under drought stress is critical to global food security.

Methods

In this study, we performed multiomics analysis in a collection of 119 diverse rapeseed (Brassica napus L.) varieties to dissect the genetic control of agronomic traits in two watering regimes [well-watered (WW) and drought stress (DS)] for 3 years. In the DS treatment, irrigation continued till the 50% pod development stage, whereas in the WW condition, it was performed throughout the whole growing season.

Results

The results of the genome-wide association study (GWAS) using 52,157 single-nucleotide polymorphisms (SNPs) revealed 1,281 SNPs associated with traits. Six stable SNPs showed sequence variation for flowering time between the two irrigation conditions across years. Three novel SNPs on chromosome C04 for plant weight were located within drought tolerance-related gene ABCG16, and their pleiotropically effects on seed weight per plant and seed yield were characterized. We identified the C02 peak as a novel signal for flowering time, harboring 52.77% of the associated SNPs. The 288-kbps LD decay distance analysis revealed 2,232 candidate genes (CGs) associated with traits. The CGs BIG1-D, CAND1, DRG3, PUP10, and PUP21 were involved in phytohormone signaling and pollen development with significant effects on seed number, seed weight, and grain yield in drought conditions. By integrating GWAS and RNA-seq, 215 promising CGs were associated with developmental process, reproductive processes, cell wall organization, and response to stress. GWAS and differentially expressed genes (DEGs) of leaf and seed in the yield contrasting accessions identified BIG1-D, CAND1, and DRG3 genes for yield variation.

Discussion

The results of our study provide insights into the genetic control of drought tolerance and the improvement of marker-assisted selection (MAS) for breeding high-yield and drought-tolerant varieties.

Identification of the domestication gene GmCYP82C4 underlying the major quantitative trait locus for the seed weight in soybean

KEY MESSAGE

A major quantitative trait locus (QTL) for the hundred-seed weight (HSW) was identified and confirmed in the two distinct soybean populations, and the target gene GmCYP82C4 underlying this locus was identified that significantly associated with soybean seed weight, and it was selected during the soybean domestication and improvement process. Soybean is a major oil crop for human beings and the seed weight is a crucial goal of soybean breeding. However, only a limited number of target genes underlying the quantitative trait loci (QTLs) controlling seed weight in soybean are known so far. In the present study, six loci associated with hundred-seed weight (HSW) were detected in the first population of 573 soybean breeding lines by genome-wide association study (GWAS), and 64 gene models were predicted in these candidate QTL regions. The QTL qHSW_1 exhibits continuous association signals on chromosome four and was also validated by region association study (RAS) in the second soybean population (409 accessions) with wild, landrace, and cultivar soybean accessions. There were seven genes in qHSW_1 candidate region by linkage disequilibrium (LD) block analysis, and only Glyma.04G035500 (GmCYP82C4) showed specifically higher expression in flowers, pods, and seeds, indicating its crucial role in the soybean seed development. Significant differences in HSW trait were detected when the association panels are genotyped by single-nucleotide polymorphisms (SNPs) in putative GmCYP82C4 promoter region. Eight haplotypes were generated by six SNPs in GmCYP82C4 in the second soybean population, and two superior haplotypes (Hap2 and Hap4) of GmCYP82C4 were detected with average HSW of 18.27 g and 18.38 g, respectively. The genetic diversity of GmCYP82C4 was analyzed in the second soybean population, and GmCYP82C4 was most likely selected during the soybean domestication and improvement process, leading to the highest proportion of Hap2 of GmCYP82C4 both in landrace and cultivar subpopulations. The QTLs and GmCYP82C4 identified in this study provide novel genetic resources for soybean seed weight trait, and the GmCYP82C4 could be used for soybean molecular breeding to develop desirable seed weight in the future.

Genome-Wide Association Studies Prioritize Genes Controlling Seed Size and Reproductive Period Length in Soybean

Hundred-seed weight (HSW) and reproductive period length (RPL) are two major agronomic traits critical for soybean production and adaptation. However, both traits are quantitatively controlled by multiple genes that have yet to be comprehensively elucidated due to the lack of major genes; thereby, the genetic basis is largely unknown. In the present study, we conducted comprehensive genome-wide association analyses (GWAS) of HSW and RPL with multiple sets of accessions that were phenotyped across different environments. The large-scale analysis led to the identification of sixty-one and seventy-four significant QTLs for HSW and RPL, respectively. An ortholog-based search analysis prioritized the most promising candidate genes for the QTLs, including nine genes (TTG2, BZR1, BRI1, ANT, KLU, EOD1/BB, GPA1, ABA2, and ABI5) for HSW QTLs and nine genes (such as AGL8, AGL9, TOC1, and COL4) and six known soybean flowering time genes (E2, E3, E4, Tof11, Tof12, and FT2b) for RPL QTLs. We also demonstrated that some QTLs were targeted during domestication to drive the artificial selection of both traits towards human-favored traits. Local adaptation likely contributes to the increased genomic diversity of the QTLs underlying RPL. The results provide additional insight into the genetic basis of HSW and RPL and prioritize a valuable resource of candidate genes that merits further investigation to reveal the complex molecular mechanism and facilitate soybean improvement.

Genome-Wide Association Study and Identification of Candidate Genes Associated with Seed Number per Pod in Soybean

Soybean (Glycine max [L.] Merr.) is one of the primary sources of plant protein and oil for human foods, animal feed, and industrial processing. The seed number per pod generally varies from one to four and is an important component of seed number per unit area and seed yield. We used natural variation in 264 landraces and improved cultivars or lines to identify candidate genes involved in the regulation of seed number per pod in soybean. Genome-wide association tests revealed 65 loci that are associated with seed number per pod trait. Among them, 11 could be detected in multiple environments. Candidate genes were identified for seed number per pod phenotype from the most significantly associated loci, including a gene encoding protein argonaute 4, a gene encoding histone acetyltransferase of the MYST family 1, a gene encoding chromosome segregation protein SMC-1 and a gene encoding exocyst complex component EXO84A. In addition, plant hormones were found to be involved in ovule and seed development and the regulation of seed number per pod in soybean. This study facilitates the dissection of genetic networks underlying seed number per pod in soybean, which will be useful for the genetic improvement of seed yield in soybean.

Understanding the Molecular Regulatory Networks of Seed Size in Soybean

Soybean being a major cash crop provides half of the vegetable oil and a quarter of the plant proteins to the global population. Seed size traits are the most important agronomic traits determining the soybean yield. These are complex traits governed by polygenes with low heritability as well as are highly influenced by the environment as well as by genotype x environment interactions. Although, extensive efforts have been made to unravel the genetic basis and molecular mechanism of seed size in soybean. But most of these efforts were majorly limited to QTL identification, and only a few genes for seed size were isolated and their molecular mechanism was elucidated. Hence, elucidating the detailed molecular regulatory networks controlling seed size in soybeans has been an important area of research in soybeans from the past decades. This paper describes the current progress of genetic architecture, molecular mechanisms, and regulatory networks for seed sizes of soybeans. Additionally, the main problems and bottlenecks/challenges soybean researchers currently face in seed size research are also discussed. This review summarizes the comprehensive and systematic information to the soybean researchers regarding the molecular understanding of seed size in soybeans and will help future research work on seed size in soybeans.

Global analysis of seed transcriptomes reveals a novel PLATZ regulator for seed size and weight control in soybean

Seed size and weight are important factors that influence soybean yield. Combining the weighted gene co-expression network analysis (WGCNA) of 45 soybean accessions and gene dynamic changes in seeds at seven developmental stages, we identified candidate genes that may control the seed size/weight. Among these, a PLATZ-type regulator overlapping with 10 seed weight QTLs was further investigated. This zinc-finger transcriptional regulator, named as GmPLATZ, is required for the promotion of seed size and weight in soybean. The GmPLATZ may exert its functions through direct binding to the promoters and activation of the expression of cyclin genes and GmGA20OX for cell proliferation. Overexpression of the GmGA20OX enhanced seed size/weight in soybean. We further found that the GmPLATZ binds to a 32-bp sequence containing a core palindromic element AATGCGCATT. Spacing of the flanking sequences beyond the core element facilitated GmPLATZ binding. An elite haplotype Hap3 was also identified to have higher promoter activity and correlated with higher gene expression and higher seed weight. Orthologues of the GmPLATZ from rice and Arabidopsis play similar roles in seeds. Our study reveals a novel module of GmPLATZ-GmGA20OX/cyclins in regulating seed size and weight and provides valuable targets for breeding of crops with desirable agronomic traits.

A CYP78As-small grain4-coat protein complex Ⅱ pathway promotes grain size in rice

CYP78A, a cytochrome P450 subfamily that includes rice (Oryza sativa L.) BIG GRAIN2 (BG2, CYP78A13) and Arabidopsis thaliana KLUH (KLU, CYP78A5), generate an unknown mobile growth signal (referred to as a CYP78A-derived signal) that increases grain (seed) size. However, the mechanism by which the CYP78A pathway increases grain size remains elusive. Here, we characterized a rice small grain mutant, small grain4 (smg4), with smaller grains than its wild type due to restricted cell expansion and cell proliferation in spikelet hulls. SMG4 encodes a multidrug and toxic compound extrusion (MATE) transporter. Loss of function of SMG4 causes smaller grains while overexpressing SMG4 results in larger grains. SMG4 is mainly localized to endoplasmic reticulum (ER) exit sites (ERESs) and partially localized to the ER and Golgi. Biochemically, SMG4 interacts with coat protein complex Ⅱ (COPⅡ) components (Sar1, Sec23, and Sec24) and CYP78As (BG2, GRAIN LENGTH 3.2 [GL3.2], and BG2-LIKE 1 [BG2L1]). Genetically, SMG4 acts, at least in part, in a common pathway with Sar1 and CYP78As to regulate grain size. In summary, our findings reveal a CYP78As-SMG4-COPⅡ regulatory pathway for grain size in rice, thus providing new insights into the molecular and genetic regulatory mechanism of grain size.

Molecular, genetic, and genomic basis of seed size and yield characteristics in soybean

Soybean (Glycine max L. Merr.) is a crucial oilseed cash crop grown worldwide and consumed as oil, protein, and food by humans and feed by animals. Comparatively, soybean seed yield is lower than cereal crops, such as maize, rice, and wheat, and the demand for soybean production does not keep up with the increasing consumption level. Therefore, increasing soybean yield per unit area is the most crucial breeding objective and is challenging for the scientific community. Moreover, yield and associated traits are extensively researched in cereal crops, but little is known about soybeans' genetics, genomics, and molecular regulation of yield traits. Soybean seed yield is a complex quantitative trait governed by multiple genes. Understanding the genetic and molecular processes governing closely related attributes to seed yield is crucial to increasing soybean yield. Advances in sequencing technologies have made it possible to conduct functional genomic research to understand yield traits' genetic and molecular underpinnings. Here, we provide an overview of recent progress in the genetic regulation of seed size in soybean, molecular, genetics, and genomic bases of yield, and related key seed yield traits. In addition, phytohormones, such as auxin, gibberellins, cytokinins, and abscisic acid, regulate seed size and yield. Hence, we also highlight the implications of these factors, challenges in soybean yield, and seed trait improvement. The information reviewed in this study will help expand the knowledge base and may provide the way forward for developing high-yielding soybean cultivars for future food demands.

WGCNA and transcriptome profiling reveal hub genes for key development stage seed size/oil content between wild and cultivated soybean

BACKGROUND

Soybean is one of the most important oil crops in the world. The domestication of wild soybean has resulted in significant changes in the seed oil content and seed size of cultivated soybeans. To better understand the molecular mechanisms of seed formation and oil content accumulation, WDD01514 (E1), ZYD00463 (E2), and two extreme progenies (E23 and E171) derived from RILs were used for weighted gene coexpression network analysis (WGCNA) combined with transcriptome analysis.

RESULTS

In this study, both seed weight and oil content in E1 and E171 were significantly higher than those in E2 and E23, and 20 DAF and 30 DAF may be key stages of soybean seed oil content accumulation and weight increase. Pathways such as "Photosynthesis", "Carbon metabolism", and "Fatty acid metabolism", were involved in oil content accumulation and grain formation between wild and cultivated soybeans at 20 and 30 DAF according to RNA-seq analysis. A total of 121 oil content accumulation and 189 seed formation candidate genes were screened from differentially expressed genes. WGCNA identified six modules related to seed oil content and seed weight, and 76 candidate genes were screened from modules and network. Among them, 16 genes were used for qRT-PCR and tissue specific expression pattern analysis, and their expression-levels in 33-wild and 23-cultivated soybean varieties were subjected to correlation analysis; some key genes were verified as likely to be involved in oil content accumulation and grain formation.

CONCLUSIONS

Overall, these results contribute to an understanding of seed lipid metabolism and seed size during seed development, and identify potential functional genes for improving soybean yield and seed oil quantity.

Genetic regulatory networks of soybean seed size, oil and protein contents

As a leading oilseed crop that supplies plant oil and protein for daily human life, increasing yield and improving nutritional quality (high oil or protein) are the top two fundamental goals of soybean breeding. Seed size is one of the most critical factors determining soybean yield. Seed size, oil and protein contents are complex quantitative traits governed by genetic and environmental factors during seed development. The composition and quantity of seed storage reserves directly affect seed size. In general, oil and protein make up almost 60% of the total storage of soybean seed. Therefore, soybean's seed size, oil, or protein content are highly correlated agronomical traits. Increasing seed size helps increase soybean yield and probably improves seed quality. Similarly, rising oil and protein contents improves the soybean's nutritional quality and will likely increase soybean yield. Due to the importance of these three seed traits in soybean breeding, extensive studies have been conducted on their underlying quantitative trait locus (QTLs) or genes and the dissection of their molecular regulatory pathways. This review summarized the progress in functional genome controlling soybean seed size, oil and protein contents in recent decades, and presented the challenges and prospects for developing high-yield soybean cultivars with high oil or protein content. In the end, we hope this review will be helpful to the improvement of soybean yield and quality in the future breeding process.

Novel Seed Size: A Novel Seed-Developing Gene in Glycine max

Soybean-seed development is controlled in multiple ways, as in many known regulating genes. Here, we identify a novel gene, Novel Seed Size (NSS), involved in seed development, by analyzing a T-DNA mutant (S006). The S006 mutant is a random mutant of the GmFTL4pro:GUS transgenic line, with phenotypes with small and brown seed coats. An analysis of the metabolomics and transcriptome combined with RT-qPCR in the S006 seeds revealed that the brown coat may result from the increased expression of chalcone synthase 7/8 genes, while the down-regulated expression of NSS leads to small seed size. The seed phenotypes and a microscopic observation of the seed-coat integument cells in a CRISPR/Cas9-edited mutant nss1 confirmed that the NSS gene conferred small phenotypes of the S006 seeds. As mentioned in an annotation on the Phytozome website, NSS encodes a potential DNA helicase RuvA subunit, and no such genes were previously reported to be involved in seed development. Therefore, we identify a novel gene in a new pathway controlling seed development in soybeans.

Fine mapping and candidate gene analysis of proportion of four-seed pods by soybean CSSLs

Soybean yield, as one of the most important and consistent breeding goals, can be greatly affected by the proportion of four-seed pods (PoFSP). In this study, QTL mapping was performed by PoFSP data and BLUE (Best Linear Unbiased Estimator) value of the chromosome segment substitution line population (CSSLs) constructed previously by the laboratory from 2016 to 2018, and phenotype-based bulked segregant analysis (BSA) was performed using the plant lines with PoFSP extreme phenotype. Totally, 5 ICIM QTLs were repeatedly detected, and 6 BSA QTLs were identified in CSSLs. For QTL (qPoFSP13-1) repeated in ICIM and BSA results, the secondary segregation populations were constructed for fine mapping and the interval was reduced to 100Kb. The mapping results showed that the QTL had an additive effect of gain from wild parents. A total of 14 genes were annotated in the delimited interval by fine mapping. Sequence analysis showed that all 14 genes had genetic variation in promoter region or CDS region. The qRT-PCR results showed that a total of 5 candidate genes were differentially expressed between the plant lines having antagonistic extreme phenotype (High PoFSP > 35.92%, low PoFSP< 17.56%). The results of haplotype analysis showed that all five genes had two or more major haplotypes in the resource population. Significant analysis of phenotypic differences between major haplotypes showed all five candidate genes had haplotype differences. And the genotypes of the major haplotypes with relatively high PoFSP of each gene were similar to those of wild soybean. The results of this study were of great significance to the study of candidate genes affecting soybean PoFSP, and provided a basis for the study of molecular marker-assisted selection (MAS) breeding and four-seed pods domestication.

Integrative genomic and transcriptomic analyses of a bud sport mutant 'Jinzao Wuhe' with the phenotype of large berries in grapevines

Background

Bud sport mutation occurs frequently in fruit plants and acts as an important approach for grapevine improvement and breeding. 'Jinzao Wuhe' is a bud sport of the elite cultivar 'Himord Seedless' with obviously enlarged organs and berries. To date, the molecular mechanisms underlying berry enlargement caused by bud sport in grapevines remain unclear.

Methods

Whole genome resequencing (WGRS) was performed for two pairs of bud sports and their maternal plants with similar phenotype to identify SNPs, InDels and structural variations (SVs) as well as related genes. Furthermore, transcriptomic sequencing at different developmental stages and weighted gene co-expression network analysis (WGCNA) for 'Jinzao Wuhe' and its maternal plant 'Himord Seedless' were carried out to identify the differentially expressed genes (DEGs), which were subsequently analyzed for Gene Ontology (GO) and function annotation.

Results

In two pairs of enlarged berry bud sports, a total of 1,334 SNPs, 272 InDels and 74 SVs, corresponding to 1,022 target genes related to symbiotic microorganisms, cell death and other processes were identified. Meanwhile, 1,149 DEGs associated with cell wall modification, stress-response and cell killing might be responsible for the phenotypic variation were also determined. As a result, 42 DEGs between 'Himord Seedless' and 'Jinzao Wuhe' harboring genetic variations were further investigated, including pectin esterase, cellulase A, cytochromes P450 (CYP), UDP-glycosyltransferase (UGT), zinc finger protein, auxin response factor (ARF), NAC transcription factor (TF), protein kinase, etc. These candidate genes offer important clues for a better understanding of developmental regulations of berry enlargement in grapevine.

Conclusion

Our results provide candidate genes and valuable information for dissecting the underlying mechanisms of berry development and contribute to future improvement of grapevine cultivars.

Systematic trait dissection in oilseed rape provides a comprehensive view, further insight, and exact roadmap for yield determination

Background

Yield is the most important and complex trait that is influenced by numerous relevant traits with very complicated interrelations. While there are a large number of studies on the phenotypic relationship and genetic basis of yield traits, systematic studies with further dissection focusing on yield are limited. Therefore, there is still lack of a comprehensive and in-depth understanding of the determination of yield.

Results

In this study, yield was systematically dissected at the phenotypic, genetic to molecular levels in oilseed rape (Brassica napus L.). The analysis of correlation, network, and principal component for 21 traits in BnaZN-RIL population showed that yield was determined by a complex trait network with key contributors. The analysis of the constructed high-density single nucleotide polymorphism (SNP) linkage map revealed the concentrated distribution of distorted and heterozygous markers, likely due to selection on genes controlling the growth period and yield heterosis. A total of 134 consensus quantitative trait loci (QTL) were identified for 21 traits, of which all were incorporated into an interconnecting QTL network with dozens of hub-QTL. Four representative hub-QTL were further dissected to the target or candidate genes that governed the causal relationships between the relevant traits.

Conclusions

The highly consistent results at the phenotypic, genetic, and molecular dissecting demonstrated that yield was determined by a multilayer composite network that involved numerous traits and genes showing complex up/down-stream and positive/negative regulation. This provides a systematic view, further insight, and exact roadmap for yield determination, which represents a significant advance toward the understanding and dissection of complex traits.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13068-022-02134-w.

An enhancing effect attributed to a nonsynonymous mutation in SOYBEAN SEED SIZE 1, a SPINDLY-like gene, is exploited in soybean domestication and improvement

Soybean (Glycine max) was domesticated from its wild relative Glycine soja. One-hundred-seed weight is one of the most important domesticated traits determining soybean yield; however, its underlying genetic basis remains elusive. We characterized a soybean seed size 1 (sss1) mutant featuring large seeds compared to its wild-type background. Positional cloning revealed that the candidate gene GmSSS1 encoded a SPINDLY homolog and was co-located in a well-identified quantitative trait locus (QTL)-rich region on chromosome 19. Knocking out GmSSS1 resulted in small seeds, while overexpressing GmSSS1/Gmsss1 induced large seeds. Modulating GmSSS1/Gmsss1 in transgenic plants can positively influence cell expansion and cell division. Relative to GmSSS1, one mutation leading to an E to Q substitution at the 182^nd residue in Gmsss1 conferred an enhancing effect on seed weight. GmSSS1 underwent diversification in wild-type and cultivated soybean, and the alleles encoding the Gmsss1-type substitution of 182^nd -Q, which originated along the central and downstream parts of the Yellow River, were selected and expanded during soybean domestication and improvement. We cloned the causative gene for the sss1 mutant, which is linked with a seed weight QTL, identified an elite allele of this gene for increasing seed weight, and provided new insights into soybean domestication and breeding.

Natural allelic variation of GmST05 controlling seed size and quality in soybean

Seed size is one of the most important agronomic traits determining the yield of crops. Cloning the key genes controlling seed size and pyramiding their elite alleles will facilitate yield improvement. To date, few genes controlling seed size have been identified in soybean, a major crop that provides half of the plant oil and one quarter of the plant protein globally. Here, through a genome-wide association study of over 1800 soybean accessions, we determined that natural allelic variation at GmST05 (Seed Thickness 05) predominantly controlled seed thickness and size in soybean germplasm. Further analyses suggested that the two major haplotypes of GmST05 differed significantly at the transcriptional level. Transgenic experiments demonstrated that GmST05 positively regulated seed size and influenced oil and protein contents, possibly by regulating the transcription of GmSWEET10a. Population genetic diversity analysis suggested that allelic variations of GmST05 were selected during geographical differentiation but have not been fixed. In summary, natural variation in GmST05 determines transcription levels and influences seed size and quality in soybean, making it an important gene resource for soybean molecular breeding.

Downregulation of a gibberellin 3β-hydroxylase enhances photosynthesis and increases seed yield in soybean

Seed yield, determined mainly by seed numbers and seed weight, is the primary target of soybean breeding. Identifying the genes underlying yield-related traits is of great significance. Through joint linkage mapping and a genome-wide association study for 100-seed weight, we cloned GmGA3ox1, a gene encoding gibberellin 3β-hydroxylase, which is the key enzyme in the gibberellin synthesis pathway. Genome resequencing identified a beneficial GmGA3ox1 haplotype contributing to high seed weight, which was further confirmed by soybean transformants. CRISPR/Cas9-generated gmga3ox1 mutants showed lower seed weight, but promoted seed yield by increasing seed numbers. The gmga3ox1 mutants reduced gibberellin biosynthesis while enhancing photosynthesis. Knockout of GmGA3ox1 resulted in the upregulation of numerous photosynthesis-related genes, particularly the GmRCA family encoding ribulose-1,5-bispho-sphate carboxylase-oxygenase (Rubisco) activases. The basic leucine zipper transcription factors GmbZIP97 and GmbZIP159, which were both upregulated in the gmga3ox1 mutants and induced by the gibberellin synthesis inhibitor uniconazole, could bind to the promoter of GmRCAβ and activate its expression. Analysis of genomic sequences with over 2700 soybean accessions suggested that GmGA3ox1 is being gradually utilized in modern breeding. Our results elucidated the important role of GmGA3ox1 in soybean yield. These findings reveal important clues for future high-yield breeding in soybean and other crops.

A lineage-specific arginine in POS1 is required for fruit size control in Physaleae (Solanaceae) via gene co-option

Solanaceae have important economic value mainly due to their edible fruits. Physalis organ size 1/cytokinin response factor 3 (POS1/CRF3), a unique gene in Solanaceae, is involved in fruit size variation in Physalis but not in Solanum. However, the underlying mechanisms remain elusive. Here, we found that POS1/CRF3 was likely created via the fusion of CRF7 and CRF8 duplicates. Multiple genetic manipulations revealed that only POS1 and Capsicum POS1 (CaPOS1) functioned in fruit size control via the positive regulation of cell expansion. Comparative studies in a phylogenetic framework showed the directional enhancement of POS1-like expression in the flowers and fruits of Physaleae and the specific gain of certain interacting proteins associated with cell expansion by POS1 and CaPOS1. A lineage-specific single nucleotide polymorphism (SNP) caused the 68th amino acid histidine in the POS1 orthologs of non-Physaleae (Nicotiana and Solanum) to change to arginine in Physaleae (Physalis and Capsicum). Substituting the arginine in Physaleae POS1-like by histidine completely abolished their function in the fruits and the protein-protein interaction (PPI) with calreticulin-3. Transcriptomic comparison revealed the potential downstream pathways of POS1, including the brassinosteroid biosynthesis pathway. However, POS1-like may have functioned ancestrally in abiotic stress within Solanaceae. Our work demonstrated that heterometric expression and a SNP caused a single amino acid change to establish new PPIs, which contributed to the co-option of POS1 in multiple regulatory pathways to regulate cell expansion and thus fruit size in Physaleae. These results provide new insights into fruit morphological evolution and fruit yield control.

The Organ Size and Morphological Change During the Domestication Process of Soybean

Soybean is one of the most important legume crops that can provide the rich source of protein and oil for human beings and livestock. In the twenty-one century, the total production of soybean is seriously behind the needs of a growing world population. Cultivated soybean [Glycine max (L.) Merr.] was domesticated from wild soybean (G. soja Sieb. and Zucc.) with the significant morphology and organ size changes in China around 5,000 years ago, including twisted stems to erect stems, small seeds to large seeds. Then it was spread worldwide to become one of the most popular and important crops. The release of the reference soybean genome and omics data provides powerful tools for researchers and breeders to dissect the functional genes and apply the germplasm in their work. Here, we summarized the function genes related to yield traits and organ size in soybean, including stem growth habit, leaf size and shape, seed size and weight. In addition, we also summarized the selection of organ traits during soybean domestication. In the end, we also discussed the application of new technology including the gene editing on the basic research and breeding of soybean, and the challenges and research hotspots in the future.

Overexpression of Rice BSR2 Confers Disease Resistance and Induces Enlarged Flowers in Torenia fournieri Lind

Plant pathogens evade basal defense systems and attack different organs and tissues of plants. Genetic engineering of plants with genes that confer resistance against pathogens is very effective in pathogen control. Conventional breeding for disease resistance in ornamental crops is difficult and lagging relative to that in non-ornamental crops due to an inadequate number of disease-resistant genes. Therefore, genetic engineering of these plants with defense-conferring genes is a practical approach. We used rice BSR2 encoding CYP78A15 for developing transgenic Torenia fournieri Lind. lines. The overexpression of BSR2 conferred resistance against two devastating fungal pathogens, Rhizoctonia solani and Botrytis cinerea. In addition, BSR2 overexpression resulted in enlarged flowers with enlarged floral organs. Histological observation of the petal cells suggested that the enlargement in the floral organs could be due to the elongation and expansion of the cells. Therefore, the overexpression of BSR2 confers broad-spectrum disease resistance and induces the production of enlarged flowers simultaneously. Therefore, this could be an effective strategy for developing ornamental crops that are disease-resistant and economically more valuable.

Increasing Fruit Weight by Editing a Cis-Regulatory Element in Tomato KLUH Promoter Using CRISPR/Cas9

CRISPR/Cas-mediated genome editing is a powerful approach to accelerate yield enhancement to feed growing populations. Most applications focus on "negative regulators" by targeting coding regions and promoters to create nulls or weak loss-of-function alleles. However, many agriculturally important traits are conferred by gain-of-function alleles. Therefore, creating gain-of-function alleles for "positive regulators" by CRISPR will be of great value for crop improvement. CYP78A family members are the positive regulators of organ weight and size in crops. In this study, we engineered allelic variation by editing tomato KLUH promoter around a single-nucleotide polymorphism (SNP) that is highly associated with fruit weight. The SNP was located in a conserved putative cis-regulatory element (CRE) as detected by the homology-based prediction and the Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq). Twenty-one mutant alleles with various insertion and deletion sizes were generated in the LA1589 background. Five mutant alleles (m2+4bp , m3+1bp , m5-1bp , m13-8bp , and m14-9bp ) showed a consistent increase in fruit weight and a significant decrease in the proportion of small fruits in all experimental evaluations. Notably, m2+4bp and m3+1bp homozygote significantly increase fruit weight by 10.7-15.7 and 8.7-16.3%, respectively. Further analysis of fruit weight based on fruit position on the inflorescence indicated that the five beneficial alleles increase the weight of all fruits along inflorescence. We also found that allele types and transcriptional changes of SlKLUH were poor predictors of the changes in fruit weight. This study not only provides a way of identifying conserved CRE but also highlights enormous potential for CRISPR/Cas-mediated cis-engineering of CYP78A members in yield improvement.

Chromosome-Scale Assembly of the Dendrobium nobile Genome Provides Insights Into the Molecular Mechanism of the Biosynthesis of the Medicinal Active Ingredient of Dendrobium

Orchids constitute approximately 10% of flowering plant species. However, only about 10 orchid genomes have been published. Metabolites are the main way through which orchids respond to their environment. Dendrobium nobile, belonging to Dendrobium, the second largest genus in Orchidaceae, has high ornamental, medicinal, and ecological value. D. nobile is the source of many popular horticultural varieties. Among the Dendrobium species, D. nobile has the highest amount of dendrobine, which is regarded as one of the criteria for evaluating medicinal quality. Due to lack of data and analysis at the genomic level, the biosynthesis pathways of dendrobine and other related medicinal ingredients in D. nobile are unknown. In this paper, we report a chromosome-scale reference genome of D. nobile to facilitate the investigation of its genomic characteristics for comparison with other Dendrobium species. The assembled genome size of D. nobile was 1.19 Gb. Of the sequences, 99.45% were anchored to 19 chromosomes. Furthermore, we identified differences in gene number and gene expression patterns compared with two other Dendrobium species by integrating whole-genome sequencing and transcriptomic analysis [e.g., genes in the polysaccharide biosynthesis pathway and upstream of the alkaloid (dendrobine) biosynthesis pathway]. Differences in the TPS and CYP450 gene families were also found among orchid species. All the above differences might contribute to the species-specific medicinal ingredient biosynthesis pathways. The metabolic pathway-related analysis will provide further insight into orchid responses to the environment. Additionally, the reference genome will provide important insights for further molecular elucidation of the medicinal active ingredients of Dendrobium and enhance the understanding of orchid evolution.

Progress in soybean functional genomics over the past decade

Soybean is one of the most important oilseed and fodder crops. Benefiting from the efforts of soybean breeders and the development of breeding technology, large number of germplasm has been generated over the last 100 years. Nevertheless, soybean breeding needs to be accelerated to meet the needs of a growing world population, to promote sustainable agriculture and to address future environmental changes. The acceleration is highly reliant on the discoveries in gene functional studies. The release of the reference soybean genome in 2010 has significantly facilitated the advance in soybean functional genomics. Here, we review the research progress in soybean omics (genomics, transcriptomics, epigenomics and proteomics), germplasm development (germplasm resources and databases), gene discovery (genes that are responsible for important soybean traits including yield, flowering and maturity, seed quality, stress resistance, nodulation and domestication) and transformation technology during the past decade. At the end, we also briefly discuss current challenges and future directions.

Genetic and Molecular Regulation Mechanisms in the Formation and Development of Vegetable Fruit Shape

Vegetable crops have a long history of cultivation worldwide and rich germplasm resources. With its continuous development and progress, molecular biology technology has been applied to various fields of vegetable crop research. Fruit is an important organ in vegetable crops, and fruit shape can affect the yield and commercialization of vegetables. In nature, fruits show differences in size and shape. Based on fruit shape diversity, the growth direction and coordination mechanism of fruits remain unclear. In this review, we discuss the latest research on fruit shape. In addition, we compare the current theories on the molecular mechanisms that regulate fruit growth, size, and shape in different vegetable families.

Fine mapping of a QTL locus (QNFSP07-1) and analysis of candidate genes for four-seeded pods in soybean

Soybean [Glycine max (L.) Merr.] is an important grain and oil crop in the world, and it is the main source of high-quality protein. The number of four-seeded pods is a quantitative trait in soybean and is closely related to yield in terms of breeding. Therefore, it is of great significance to study the inheritance of four-seed pods and to excavate related genes for improving soybean yield. In this study, individuals with high ratio of four-seed pods which from chromosome segment substitution lines (CSSLs) that can be stably inherited were selected as the parent, and Suinong 14 (SN14) was used as recurrent parent to construct secondary mapping population via marker-assisted selection. From 2006 to 2017, QTL analysis was performed using secondary mapping populations, and the initial QTL mapping interval was 0.67 Mb and was located on Gm07. Based on the initial QTL mapping results, individuals that were heterozygous at the interval (36,116,118-37,399,738 bp) were screened in 2018, and the heterozygous individuals were subjected to inbreeding to obtain 13 F3 populations, with a target interval of 321 kb. Gene annotation was performed on the fine mapping interval, and 27 genes were obtained. Among 27 genes, Glyma.07G200900 and Glyma.07G201200 were identified as candidate genes. qRT-PCR was used to measure the expression of the 2 candidate genes at different developmental stages of soybean, and the expression levels of the 2 candidate genes in terms of cell division (axillary buds, COTs, EMs) were higher than those in terms of cell expansion (MM, LM), and these genes play a positive regulatory role in the formation of four-seeded pods. Haplotype analysis of 2 candidate genes which shows that Glyma.07G201200 has two excellent haplotypes, and the significance level between the two excellent haplotypes at p < 0.05. Those results provide the information for gene map-based cloning and molecular marker-assisted breeding of the number of four-seeded pod in soybean.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01265-6.

GmFULa improves soybean yield by enhancing carbon assimilation without altering flowering time or maturity

KEY MESSAGE

GmFULa improved soybean yield by enhancing carbon assimilation. Meanwhile, different from known yield-related genes, it did not alter flowering time or maturity. Soybean (Glycine max (L.) Merr.) is highly demanded by a continuously growing human population. However, increasing soybean yield is a major challenge. FRUITFULL (FUL), a MADS-box transcription factor, plays important roles in multiple developmental processes, especially fruit and pod development, which are crucial for soybean yield formation. However, the functions of its homologs in soybean are not clear. Here, through haplotype analysis, we found that one haplotype of the soybean homolog GmFULa (GmFULa-H02) is dominant in cultivated soybeans, suggesting that GmFULa-H02 was highly selected during domestication and varietal improvement of soybean. Interestingly, transgenic overexpression of GmFULa enhanced vegetative growth with more biomass accumulated and ultimately increased the yield but without affecting the plant height or changing the flowering time and maturity, indicating that it enhances the efficiency of dry matter accumulation. It also promoted the yield factors like branch number, pod number and 100-seed weight, which ultimately increased the yield. It increased the palisade tissue cell number and the chlorophyll content to promote photosynthesis and increase the soluble sugar content in leaves and fresh seeds. Furthermore, GmFULa were found to be sublocalized in the nucleus and positively regulate sucrose synthases (SUSs) and sucrose transporters (SUTs) by binding with the conserved CArG boxes in their promoters. Overall, these results showed GmFULa promotes the capacity of assimilation and the transport of the resultant assimilates to increase yield, and provided insights into the link between GmFULa and sucrose synthesis with transport-related molecular pathways that control seed yield.

A Major and Stable Quantitative Trait Locus qSS2 for Seed Size and Shape Traits in a Soybean RIL Population

Seed size and shape traits are important determinants of seed yield and appearance quality in soybean [Glycine max (L.) Merr.]. Understanding the genetic architecture of these traits is important to enable their genetic improvement through efficient and targeted selection in soybean breeding, and for the identification of underlying causal genes. To map seed size and shape traits in soybean, a recombinant inbred line (RIL) population developed from K099 (small seed size) × Fendou 16 (large seed size), was phenotyped in three growing seasons. A genetic map of the RIL population was developed using 1,485 genotyping by random amplicon sequencing-direct (GRAS-Di) and 177 SSR markers. Quantitative trait locus (QTL) mapping was conducted by inclusive composite interval mapping. As a result, 53 significant QTLs for seed size traits and 27 significant QTLs for seed shape traits were identified. Six of these QTLs (qSW8.1, qSW16.1, qSLW2.1, qSLT2.1, qSWT1.2, and qSWT4.3) were identified with LOD scores of 3.80-14.0 and R ² of 2.36%-39.49% in at least two growing seasons. Among the above significant QTLs, 24 QTLs were grouped into 11 QTL clusters, such as, three major QTLs (qSL2.3, qSLW2.1, and qSLT2.1) were clustered into a major QTL on Chr.02, named as qSS2. The effect of qSS2 was validated in a pair of near isogenic lines, and its candidate genes (Glyma.02G269400, Glyma.02G272100, Glyma.02G274900, Glyma.02G277200, and Glyma.02G277600) were mined. The results of this study will assist in the breeding programs aiming at improvement of seed size and shape traits in soybean.

Differential expression of SlKLUH controlling fruit and seed weight is associated with changes in lipid metabolism and photosynthesis-related genes

The sizes of plant organs such as fruit and seed are crucial yield components. Tomato KLUH underlies the locus fw3.2, an important regulator of fruit and seed weight. However, the mechanism by which the expression levels of KLUH affect organ size is poorly understood. We found that higher expression of SlKLUH increased cell proliferation in the pericarp within 5 d post-anthesis in tomato near-isogenic lines. Differential gene expression analyses showed that lower expression of SlKLUH was associated with increased expression of genes involved in lipid metabolism. Lipidomic analysis revealed that repression of SlKLUH mainly increased the contents of certain non-phosphorus glycerolipids and phospholipids and decreased the contents of four unknown lipids. Co-expression network analyses revealed that lipid metabolism was possibly associated with but not directly controlled by SlKLUH, and that this gene instead controls photosynthesis-related processes. In addition, many transcription factors putatively involved in the KLUH pathway were identified. Collectively, we show that SlKLUH regulates fruit and seed weight which is associated with altered lipid metabolism. The results expand our understanding of fruit and seed weight regulation and offer a valuable resource for functional studies of candidate genes putatively involved in regulation of organ size in tomato and other crops.

An evolutionary population structure model reveals pleiotropic effects of GmPDAT for traits related to seed size and oil content in soybean

Seed oil traits in soybean that are of benefit to human nutrition and health have been selected for during crop domestication. However, these domesticated traits have significant differences across various evolutionary types. In this study, we found that the integration of evolutionary population structure (evolutionary types) with genome-wide association studies increased the power of gene detection, and it identified one locus for traits related to seed size and oil content on chromosome 13. This domestication locus, together with another one in a 200-kb region, was confirmed by the GEMMA and EMMAX software. The candidate gene, GmPDAT, had higher expressional levels in high-oil and large-seed accessions than in low-oil and small-seed accessions. Overexpression lines had increased seed size and oil content, whereas RNAi lines had decreased seed size and oil content. The molecular mechanism of GmPDAT was deduced based on results from linkage analysis for triacylglycerols and on histocytological comparisons of transgenic soybean seeds. Our results illustrate a new approach for identifying domestication genes with pleiotropic effects.

Identification of Additive-Epistatic QTLs Conferring Seed Traits in Soybean Using Recombinant Inbred Lines

Seed weight and shape are important agronomic traits that affect soybean quality and yield. In the present study, we used image analysis software to evaluate 100-seed weight and seed shape traits (length, width, perimeter, projection area, length/width, and weight/projection area) of 155 novel recombinant inbred soybean lines (NJRISX) generated by crossing "Su88-M21" and "XYXHD". We examined quantitative trait loci (QTLs) associated with the six traits (except seed weight per projection area), and identified 42 additive QTLs (5-8 QTLs per trait) accounting for 24.9-37.5% of the phenotypic variation (PV). Meanwhile, 2-4 epistatic QTL pairs per trait out of a total of 18 accounted for 2.5-7.2% of the PV; and unmapped minor QTLs accounted for the remaining 35.0-56.7% of the PV. A total of 28 additive and 11 epistatic QTL pairs were concentrated in nine joint QTL segments (JQSs), indicating that QTLs associated with seed weight and shape are closely related and interacted. An interaction was also detected between additive and epistatic QTL pairs and environment, which made significant contributions of 1.4-9.5% and 0.4-0.8% to the PV, respectively. We annotated 18 candidate genes in the nine JQSs, which were important for interpreting the close relationships among the six traits. These findings indicate that examining the interactions between closely related traits rather than only analyzing individual trait provides more useful insight into the genetic system of the interrelated traits for which there has been limited QTL information.

Enhanced resistance to fungal and bacterial diseases in tomato and Arabidopsis expressing BSR2 from rice

The overexpression of rice BSR2 would offer a simple and effective strategy to protect plants from multiple devastating diseases in tomato and Arabidopsis. Many devastating plant diseases are caused by pathogens possessing a wide host range. Fungal Botrytis cinerea and Rhizoctonia solani, as well as bacterial Pseudomonas syringae and Ralstonia pseudosolanacearum are four such pathogens that infect hundreds of plant species, including agronomically important crops, and cause serious diseases, leading to severe economic losses. However, reports of genes that can confer resistance to broad host-range pathogens via traditional breeding methods are currently limited. We previously reported that Arabidopsis plants overexpressing rice BROAD-SPECTRUM RESISTANCE2 (BSR2/CYP78A15) showed tolerance not only to bacterial P. syringae pv. tomato DC3000 but also to fungal Colletotrichum higginsianum and R. solani. Rice plants overexpressing BSR2 displayed tolerance to two R. solani anastomosis groups. In the present study, first, BSR2-overexpressing (OX) Arabidopsis plants were shown to be additionally tolerant to B. cinerea, R. solani, and R. pseudosolanacearum. Next, tomato 'Micro-Tom' was used as a model to determine whether such tolerance by BSR2 can be introduced into dicot crops to prevent infection from pathogens possessing wide host range. BSR2-OX tomato displayed broad-spectrum disease tolerance to fungal B. cinerea and R. solani, as well as to bacterial P. syringae and R. pseudosolanacearum. Additionally, undesirable traits such as morphological changes were not detected. Thus, BSR2 overexpression can offer a simple and effective strategy to protect crops from multiple destructive diseases.

Overexpression of an ethylene-forming ACC oxidase (ACO) gene precedes the Minute Hilum seed coat phenotype in Glycine max

BACKGROUND

To elucidate features of seed development, we investigated the transcriptome of a soybean isoline from the germplasm collection that contained an introgressed allele known as minute hilum (mi) which confers a smaller hilum region where the seed attaches to the pod and also results in seed coat cracking surrounding the hilum region.

RESULTS

RNAs were extracted from immature seed from an extended hilum region (i.e., the hilum and a small ring of tissue surrounding the hilum in which the cracks form) at three different developmental stages:10-25, 25-50 and 50-100 mg seed fresh weight in two independent replicates for each stage. The transcriptomes of these samples from both the Clark isoline containing the mi allele (PI 547628, UC413, ii R t mi G), and its recurrent Clark 63 parent isoline (PI 548532, UC7, ii R T Mi g), which was used for six generations of backcrossing, were compared for differential expression of 88,648 Glyma models of the soybean genome Wm82.a2. The RNA sequence data obtained from the 12 cDNA libraries were subjected to padj value < 0.05 and at least two-fold expression differences to select with confidence genes differentially expressed in the hilum-containing tissue of the seed coat between the two lines. Glyma.09G008400 annotated as encoding an ethylene forming enzyme, ACC oxidase (ACO), was found to be highly overexpressed in the mi hilum region at 165 RPKMs (reads per kilobase per million mapped reads) compared to the standard line at just 0.03 RPKMs. Evidence of changes in expression of genes downstream of the ethylene pathway included those involved in auxin and gibberellin hormone action and extensive differences in expression of cell wall protein genes. These changes are postulated to determine the restricted hilum size and cracking phenotypes.

CONCLUSIONS

We present transcriptome and phenotypic evidence that substantially higher expression of an ethylene-forming ACO gene likely shifts hormone balance and sets in motion downstream changes resulting in a smaller hilum phenotype and the cracks observed in the minute hilum (mi) isoline as compared to its recurrent parent.

Three-dimensional genetic networks among seed oil-related traits, metabolites and genes reveal the genetic foundations of oil synthesis in soybean

Although the biochemical and genetic basis of lipid metabolism is clear in Arabidopsis, there is limited information concerning the relevant genes in Glycine max (soybean). To address this issue, we constructed three-dimensional genetic networks using six seed oil-related traits, 52 lipid metabolism-related metabolites and 54 294 SNPs in 286 soybean accessions in total. As a result, 284 and 279 candidate genes were found to be significantly associated with seed oil-related traits and metabolites by phenotypic and metabolic genome-wide association studies and multi-omics analyses, respectively. Using minimax concave penalty (MCP) and smoothly clipped absolute deviation (SCAD) analyses, six seed oil-related traits were found to be significantly related to 31 metabolites. Among the above candidate genes, 36 genes were found to be associated with oil synthesis (27 genes), amino acid synthesis (four genes) and the tricarboxylic acid (TCA) cycle (five genes), and four genes (GmFATB1a, GmPDAT, GmPLDα1 and GmDAGAT1) are already known to be related to oil synthesis. Using this information, 133 three-dimensional genetic networks were constructed, 24 of which are known, e.g. pyruvate-GmPDAT-GmFATA2-oil content. Using these networks, GmPDAT, GmAGT and GmACP4 reveal the genetic relationships between pyruvate and the three major nutrients, and GmPDAT, GmZF351 and GmPgs1 reveal the genetic relationships between amino acids and seed oil content. In addition, GmCds1, along with average temperature in July and the rainfall from June to September, influence seed oil content across years. This study provides a new approach for the construction of three-dimensional genetic networks and reveals new information for soybean seed oil improvement and the identification of gene function.

Marker-Trait Association Analysis of Seed Traits in Accessions of Common Bean (Phaseolus vulgaris L.) in China

Seed weight and seed size are the key agronomic traits that determine yield in common bean. To investigate the genetic architecture of four seed traits (100-seed weight, seed length, seed width, and seed height) of common bean in China, marker-trait association analysis of these seed traits was performed in a nationwide population of 395 common bean accessions using 116 polymorphic SSR markers. The four seed traits were evaluated in six trials across three environments. Seed size varied among the environments. Population structure was evaluated based on SSR markers and phaseolin, which divided the accessions into two main subpopulations representing the two known gene pools. Seed weight and seed size had a strong relationship with population clustering. In addition, in a Genome-wide association studies (GWAS), 21 significantly associated markers were identified for the four seed traits with two models, namely, general linear model (GLM) and mixed linear model (MLM). Some markers had pleiotropic effects, i.e., controlled more than one trait. The significant quantitative trait loci identified in this study could be used in marker-assisted breeding to accelerate the genetic improvement of yield in common bean.

Genome-Wide Association Study of Grain Size Traits in Indica Rice Multiparent Advanced Generation Intercross (MAGIC) Population

Rice grain size plays a crucial role in determining grain quality and yield. In this study, two multiparent advanced generation intercross (MAGIC) populations, DC1 and BIM, were evaluated for grain size across three environments and genotyped with 55K array-based SNP detection and genotype-by-sequencing (GBS), respectively, to identify QTLs and SNPs associated with grain length, grain width, grain length-width ratio, grain thickness, and thousand grain weight. A total of 18 QTLs were identified for the five grain size-related traits and explained 6.43-63.35% of the total phenotypic variance. Twelve of these QTLs colocalized with the cloned genes, GS3, GW5/qSW5, GW7/GL7/SLG7, and GW8/OsSPL16, of which the first two genes showed the strongest effect for grain length and grain width, respectively. Four potential new genes were also identified from the QTLs, which exhibited both genetic background independency and environment stability and could be validated in future studies. Moreover, the significant SNP markers identified are valuable for direct utilization in marker-assisted breeding to improve rice grain size.

A genome-wide association and meta-analysis reveal regions associated with seed size in cowpea [Vigna unguiculata (L.) Walp]

This paper combined GWAS, meta-analysis and sequence homology comparison with common bean to identify regions associated with seed size variation in domesticated cowpea. Seed size is an important trait for yield and commercial value in dry-grain cowpea. Seed size varies widely among different cowpea accessions, and the genetic basis of such variation is not yet well understood. To better decipher the genetic basis of seed size, a genome-wide association study (GWAS) and meta-analysis were conducted on a panel of 368 cowpea diverse accessions from 51 countries. Four traits, including seed weight, length, width and density were evaluated across three locations. Using 51,128 single nucleotide polymorphisms covering the cowpea genome, 17 loci were identified for these traits. One locus was common to weight, width and length, suggesting pleiotropy. By integrating synteny-based analysis with common bean, six candidate genes (Vigun05g036000, Vigun05g039600, Vigun05g204200, Vigun08g217000, Vigun11g187000, and Vigun11g191300) which are implicated in multiple functional categories related to seed size such as endosperm development, embryo development, and cell elongation were identified. These results suggest that a combination of GWAS meta-analysis with synteny comparison in a related plant is an efficient approach to identify candidate gene (s) for complex traits in cowpea. The identified loci and candidate genes provide useful information for improving cowpea varieties and for molecular investigation of seed size.

Genome-wide association study of inflorescence length of cultivated soybean based on the high-throughout single-nucleotide markers

As an important and complex trait, inflorescence length (IL) of soybean [Glycine max (L.) Merr.] significantly affected seed yields. Therefore, elucidating molecular basis of inflorescence architecture, especially for IL, was important for improving soybean yield potentials. Longer IL meaned to have more pod and seed in soybean. Hence, increasing IL and improving yield are targets for soybean breeding. In this study, a association panel, comprising 283 diverse samples, was used to dissect the genetic basis of IL based on genome-wide association analysis (GWAS) and haplotype analysis. GWAS and haplotype analysis were conducted through high-throughout single-nucleotide polymorphisms (SNP) developed by SLAF-seq methodology. A total of 39, 057 SNPs (minor allele frequency ≥ 0.2 and missing data ≤ 10%) were utilized to evaluate linkage disequilibrium (LD) level in the tested association panel. A total of 30 association signals were identified to be associated with IL via GWAS. Among them, 13 SNPs were novel, and another 17 SNPs were overlapped or located near the linked regions of known quantitative trait nucleotide (QTN) with soybean seed yield or yield component. The functional genes, located in the 200-kb genomic region of each peak SNP, were considered as candidate genes, such as the cell division/ elongation, specific enzymes, and signaling or transport of specific proteins. These genes have been reported to participant in the regulation of IL. Ten typical long-IL lines and ten typical short-IL lines were re-sequencing, and then, six SNPs from five genes were obtained based on candidate gene-based association. In addition, 42 haplotypes were defined based on haplotype analysis. Of them, 11 haplotypes were found to regulate long IL (> 14 mm) in soybean. The identified 30 QTN with beneficial alleles and their candidate genes might be valuable for dissecting the molecular mechanisms of IL and further improving the yield potential of soybean.

A CACTA-like transposable element in the upstream region of BnaA9.CYP78A9 acts as an enhancer to increase silique length and seed weight in rapeseed

Rapeseed (Brassica napus L.) is a model plant for polyploid crop research and the second-leading source of vegetable oil worldwide. Silique length (SL) and seed weight are two important yield-influencing traits in rapeseed. Using map-based cloning, we isolated qSLWA9, which encodes a P450 monooxygenase (BnaA9.CYP78A9) and functions as a positive regulator of SL. The expression level of BnaA9.CYP78A9 in silique valves of the long-silique variety is much higher than that in the regular-silique variety, which results in elongated cells and a prolonged phase of silique elongation. Plants of the long-silique variety and transgenic plants with high expression of BnaA9.CYP78A9 had a higher concentration of auxin in the developing silique; this induced a number of auxin-related genes but no genes in well-known auxin biosynthesis pathways, suggesting that BnaA9.CYP78A9 may influence auxin concentration by affecting auxin metabolism or an unknown auxin biosynthesis pathway. A 3.7-kb CACTA-like transposable element (TE) inserted in the 3.9-kb upstream regulatory sequence of BnaA9.CYP78A9 elevates the expression level, suggesting that the CACTA-like TE acts as an enhancer to stimulate high gene expression and silique elongation. Marker and sequence analysis revealed that the TE in B. napus had recently been introgressed from Brassica rapa by interspecific hybridization. The insertion of the TE is consistently associated with long siliques and large seeds in both B. napus and B. rapa collections. However, the frequency of the CACTA-like TE in rapeseed varieties is still very low, suggesting that this allele has not been widely used in rapeseed breeding programs and would be invaluable for yield improvement in rapeseed breeding.

Molecular Networks of Seed Size Control in Plants

The size of seeds affects not only evolutionary fitness but also grain yield of crops. Understanding the mechanisms controlling seed size has become an important research field in plant science. Seed size is determined by the integrated signals of maternal and zygotic tissues, which control the coordinated growth of the embryo, endosperm, and seed coat. Recent advances have identified several signaling pathways that control seed size through maternal tissues, including or involving the ubiquitin-proteasome pathway, G-protein signaling, mitogen-activated protein kinase (MAPK) signaling, phytohormone perception and homeostasis, and some transcriptional regulators. Meanwhile, growth of the zygotic tissues is regulated in part by the HAIKU (IKU) pathway and phytohormones. This review provides a general overview of current findings in seed size control and discusses the emerging molecular mechanisms and regulatory networks found to be involved.

Maternal control of seed weight in rapeseed (Brassica napus L.): the causal link between the size of pod (mother, source) and seed (offspring, sink)

Seed size/weight is one of the key traits related to plant domestication and crop improvement. In rapeseed (Brassica napus L.) germplasm, seed weight shows extensive variation, but its regulatory mechanism is poorly understood. To identify the key mechanism of seed weight regulation, a systematic comparative study was performed. Genetic, morphological and cytological evidence showed that seed weight was controlled by maternal genotype, through the regulation of seed size mainly via cell number. The physiological evidence indicated that differences in the pod length might result in differences in pod wall photosynthetic area, carbohydrates and the final seed weight. We also identified two pleiotropic major quantitative trait loci that acted indirectly on seed weight via their effects on pod length. RNA-seq results showed that genes related to pod development and hormones were significantly differentially expressed in the pod wall; genes related to development, cell division, nutrient reservoir and ribosomal proteins were all up-regulated in the seeds of the large-seed pool. Finally, we proposed a potential seed weight regulatory mechanism that is specific to rapeseed and novel in plants. The results demonstrate a causal link between the size of the pod (mother, source) and the seed (offspring, sink) in rapeseed, which provides novel insight into the maternal control of seed weight and will open a new research field in plants.

The rice CYP78A gene BSR2 confers resistance to Rhizoctonia solani and affects seed size and growth in Arabidopsis and rice

The fungal pathogen Rhizoctonia solani causes devastating diseases in hundreds of plant species. Among these, R. solani causes sheath blight, one of the three major diseases in rice. To date, few genes have been reported that confer resistance to R. solani. Here, rice-FOX Arabidopsis lines identified as having resistance to a bacterial pathogen, Pseudomonas syringae pv. tomato DC3000, and a fungal pathogen, Colletotrichum higginsianum were screened for disease resistance to R. solani. BROAD-SPECTRUM RESISTANCE2 (BSR2), a gene encoding an uncharacterized cytochrome P450 protein belonging to the CYP78A family, conferred resistance to R. solani in Arabidopsis. When overexpressed in rice, BSR2 also conferred resistance to two R. solani anastomosis groups. Both Arabidopsis and rice plants overexpressing BSR2 had slower growth and produced longer seeds than wild-type control plants. In contrast, BSR2-knockdown rice plants were more susceptible to R. solani and displayed faster growth and shorter seeds in comparison with the control. These results indicate that BSR2 is associated with disease resistance, growth rate and seed size in rice and suggest that its function is evolutionarily conserved in both monocot rice and dicot Arabidopsis.

Evolution and Expression Divergence of the CYP78A Subfamily Genes in Soybean

Gene expression divergence is an important evolutionary driving force for the retention of duplicate genes. In this study, we identified three CYP78A subfamily genes in soybean, GmCYP78A70, GmCYP78A57 and GmCYP78A72, which experienced different duplication events. GmCYP78A70 was mainly expressed in leaf tissue and the vegetative phase, whereas GmCYP78A57 was mainly expressed in floral tissue and seed, i.e., the reproductive phase. Expression of GmCYP78A72 could be detected in all the tissues and phases mentioned above. The expression levels of GmCYP78A70 and GmCYP78A57 in different soybean cultivars showed positive correlations with leaf size and 100-seed weight, respectively. The population genetics analysis indicated that the three genes had experienced different selective pressures during domestication and improved breeding of soybean. Deciphering the function of this subfamily of genes may well prove useful to breeders for improving soybean’s agronomic traits.

Global identification, structural analysis and expression characterization of cytochrome P450 monooxygenase superfamily in rice

BACKGROUND

The cytochrome P450 monooxygenases (CYP450, CYP, P450) catalyze numerous monooxygenation/hydroxylation reactions in biochemical pathways. Although CYP superfamily has been systematically studied in a few species, the genome-scale research about it in rice has not been done.

RESULTS

In this study, a total of 355 CYPs encoded by 326 genes were identified in japonica genome. The OsCYP genes are classified into 10 clans including 45 families according to phylogenetic analysis. More than half of the genes are distributed in 53 tandem duplicated gene clusters. Intron-exon structure of OsCYPs exhibits highly conserved and specificity within a family, and divergences of duplicate genes in gene structure result in non-functionalization, neo-functionalization or sub-functionalization. Selection pressure analysis showed that rice CYPs are under purifying selection. The microarray data analysis shows that some genes are tissue-specific expression, such as OsCYP710A5 and OsCYP71X14 in endosperm, OsCYP99A3 and OsCYP78A16 in root and OsCYP93G2 and OsCYP97D7 in leaf. Analysis of RNA-seq data derived from rice leaf developmental gradient indicates that some OsCYPs exhibit zone-specific expression patterns. OsCYP87C2, OsCYP96B5, OsCYP96B8 and OsCYP84A5 were specifically expressed in leaf base and transitional zone. The transcripts of lineages II and IV-1 members were highly abundant in maturing zone. Eighty three OsCYPs are differentially expressed in response to drought stress, of which OsCYP51G3, OsCYP709C9, OsCYP709C5, OsCYP81A6, OsCYP72A18 and OsCYP704A5 are strongly induced and OsCYP78A16, OsCYP89C9 and OsCYP704A5 are down-regulated significantly, and some of the results were validated by qPCR. And 23 up-regulated and 17 down-regulated genes are specific to Osbhlh148 mutation under drought stress. Compared to those in wild type, the changes in transcript levels of several genes are slight in the mutant, such as OsCYP51G3, OsCYP94C2, OsCYP709C9 and OsCYP709C5.

CONCLUSION

The whole-genomic analysis of rice P450 superfamily provides a clue to understanding biological function of OsCYPs in development regulation and drought stress response, and is helpful to rice molecular breeding.

Differential expression of a WRKY gene between wild and cultivated soybeans correlates to seed size

Soybean (Glycine max) probably originated from the wild soybean (Glycine soja). Glycine max has a significantly larger seed size, but the underlying genomic changes are largely unknown. Candidate regulatory genes were preliminarily proposed by data co-localizing RNA sequencing with the quantitative loci (QTLs) for seed size. The soybean gene locus SoyWRKY15a and its orthologous genes from G. max (GmWRKY15a) and G. soja (GsWRKY15a) were analyzed in detail. The coding sequences were nearly identical between the two orthologs, but GmWRKY15a was significantly more highly expressed than GsWRKY15a. Four haplotypes (H1-H4) were found and they varied in the size of a CT-core microsatellite locus in the 5'-untranslated region of this gene. H1 (with six CT-repeats) was the only allelic version found in G. max, while H3 (with five CT-repeats) was the dominant G. soja allele. Differential expression of this gene in soybean pods was correlated with CT-repeat variation, and manipulation of the CT copy number altered the reporter gene expression, suggesting a regulatory role for the simple sequence repeats. Seed weight of wild soybeans harboring H1 was significantly greater than that of soybeans having haplotypes H2, H3, or H4, and seed weight was correlated with gene expression, suggesting the influence of GsWRKY15a in controlling seed size. However, the seed size might be refractory to increased SoyWRKY15a expression in cultivated soybeans. The evolutionary significance of SoyWRKY15a variation in soybean seed domestication is discussed.

Soybean domestication: the origin, genetic architecture and molecular bases

Domestication provides an important model for the study of evolution, and information learned from domestication research aids in the continued improvement of crop species. Recent progress in de novo assembly and whole-genome resequencing of wild and cultivated soybean genomes, in addition to new archeological discoveries, sheds light on the origin of this important crop and provides a clearer view on the modes of artificial selection that drove soybean domestication and diversification. This novel genomic information enables the search for polymorphisms that underlie variation in agronomic traits and highlights genes that exhibit a signature of selection, leading to the identification of a number of candidate genes that may have played important roles in soybean domestication, diversification and improvement. These discoveries provide a novel point of comparison on the evolutionary bases of important agronomic traits among different crop species.

Identification of regulatory networks and hub genes controlling soybean seed set and size using RNA sequencing analysis

To understand the gene expression networks controlling soybean seed set and size, transcriptome analyses were performed in three early seed developmental stages, using two genotypes with contrasting seed size. The two-dimensional data set provides a comprehensive and systems-level view on dynamic gene expression networks underpinning soybean seed set and subsequent development. Using pairwise comparisons and weighted gene coexpression network analyses, we identified modules of coexpressed genes and hub genes for each module. Of particular importance are the discoveries of specific modules for the large seed size variety and for seed developmental stages. A large number of candidate regulators for seed size, including those involved in hormonal signaling pathways and transcription factors, were transiently and specifically induced in the early developmental stages. The soybean homologs of a brassinosteroid signaling receptor kinase, a brassinosteroid-signaling kinase, were identified as hub genes operating in the seed coat network in the early seed maturation stage. Overexpression of a candidate seed size regulatory gene, GmCYP78A5, in transgenic soybean resulted in increased seed size and seed weight. Together, these analyses identified a large number of potential key regulators controlling soybean seed set, seed size, and, consequently, yield potential, thereby providing new insights into the molecular networks underlying soybean seed development.

Altered expression of maize PLASTOCHRON1 enhances biomass and seed yield by extending cell division duration

Maize is the highest yielding cereal crop grown worldwide for grain or silage. Here, we show that modulating the expression of the maize PLASTOCHRON1 (ZmPLA1) gene, encoding a cytochrome P450 (CYP78A1), results in increased organ growth, seedling vigour, stover biomass and seed yield. The engineered trait is robust as it improves yield in an inbred as well as in a panel of hybrids, at several locations and over multiple seasons in the field. Transcriptome studies, hormone measurements and the expression of the auxin responsive DR5rev:mRFPer marker suggest that PLA1 may function through an increase in auxin. Detailed analysis of growth over time demonstrates that PLA1 stimulates the duration of leaf elongation by maintaining dividing cells in a proliferative, undifferentiated state for a longer period of time. The prolonged duration of growth also compensates for growth rate reduction caused by abiotic stresses.

Arabidopsis KLU homologue GmCYP78A72 regulates seed size in soybean

Soybean (Glycine max) is one of the most important crops in the world, and its yield is largely determined by grain weight and grain size. However, the genes that regulate soybean seed size have not been identified. CYP78A, which is highly conserved within terrestrial plants, regulates organ development. In Arabidopsis, AtCYP78A5/KLU has been shown to determine seed size. In the present study, soybean CYP78A72 (GmCYP78A72), one of the orthologs of KLU, was over-expressed in both Arabidopsis and soybean to examine its function in plant development. GmCYP78A72 heterologous expression in Arabidopsis resulted in enlarged sepals, petals, seeds and carpel. Over-expression of GmCYP78A72 in soybean resulted in increased pea size, which is an extremely desirable trait for enhancing productivity. Moreover, knock-down of GmCYP78A72 does not reduce grain size. However, silencing of GmCYP78A57, GmCYP78A70 and GmCYP78A72 genes in triplet reduces the seed size significantly indicating functional redundancy of these three GmCYP78A genes. In conclusion, we investigated the role of CYP78A in soybean seed regulation, and our strategy can be effectively used to engineer large seed traits in soybean varieties as well as other crops.