Seed weight differences between wild and domesticated soybeans are associated with specific changes in gene expression

Exploration of GmDof11-lncRNA13082 Module Regulating Oil Synthesis in Plants

Soybean (Glycine max [Linn.] Merr.) is an important oilseed crop. Although transcription factors (TFs) can coordinate the expression of mRNA and lncRNA, their coordination in the soybean oil synthesis pathway remains unclear. This study examined the interaction between the TF GmDof11 and lncRNA13082 and found that overexpression of GmDof11 led to an increase in the number of Arabidopsis seeds, thousand seed weight, crude protein, hydrolysis amino acid, and soluble sugar. Additionally, it reduced the triglyceride and starch contents and affected the proportion of fatty acids, increasing the contents of palmitic acid, stearic acid, and linolenic acid. The yeast two-hybrid experiments revealed that GmDof11 interacts with GmBCCP1, GmLEC1b, and GmFAB2 proteins. In the RT-qPCR analysis of transgenic soybean roots, it was found that GmDof11 can activate the production of lncRNA13082 and work in conjunction with lncRNA13082 to oversee oil synthesis and nutrient storage. Our research provides robust theoretical evidence for a comprehensive resolution of TF-lncRNA regulation in the soybean oil synthesis network.

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.

Emerging roles of long non-coding RNAs in regulating agriculturally important seed traits

Seeds have enormous economic importance as source of calories, nutrition, edible oil, and biofuels. Therefore, seed traits like seed size and shape, weight, micronutrient content, oil content, quality, post-harvest sprouting, etc., are some of the main targets in crop improvement. Designing the strategies for their improvement benefits heavily from understanding the regulatory aspects of seed development. Recent studies indicate that long non-coding RNAs (lncRNAs) are one of the important regulators of seed development. They played a significant role in crop domestication by influencing seed traits. LncRNAs are conventionally defined as non-coding RNAs greater than 200 bp in length but lacking protein coding potential. Here we highlight the emerging pieces of evidence of lncRNA-mediated regulation of seed development through diverse mechanisms, for instance, by acting as target mimics or precursors of regulatory small RNAs or through chromatin remodeling and post-transcriptional repression. We also enumerate the insights from high-throughput transcriptomic studies from developing seeds of cereal, oilseed, biofuel, and pulse crops. We highlight the lncRNA candidates and lncRNA-mediated regulatory networks regulating seed development and related agronomic traits. Further, we discuss the potential of lncRNAs for improvement of agriculturally important seed traits through marker-assisted breeding and/or transgenic approaches.

Genome-wide analysis suggests the potential role of lncRNAs during seed development and seed size/weight determination in chickpea

The integrated transcriptome data analyses suggested the plausible roles of lncRNAs during seed development in chickpea. The candidate lncRNAs associated with QTLs and those involved in miRNA-mediated seed size/weight determination in chickpea have been identified. Long non-coding RNAs (lncRNAs) are important regulators of various biological processes. Here, we identified lncRNAs at seven successive stages of seed development in small-seeded and large-seeded chickpea cultivars. In total, 4751 lncRNAs implicated in diverse biological processes were identified. Most of lncRNAs were conserved between the two cultivars, whereas only a few of them were conserved in other plants, suggesting their species-specificity. A large number of lncRNAs differentially expressed between the two chickpea cultivars associated with seed development-related processes were identified. The lncRNAs acting as precursors of miRNAs and those mimicking target protein-coding genes of miRNAs involved in seed size/weight determination, including HAIKU1, BIG SEEDS1, and SHB1, were also revealed. Further, lncRNAs located within seed size/weight associated quantitative trait loci were also detected. Overall, we present a comprehensive resource and identified candidate lncRNAs that may play important roles during seed development and seed size/weight determination in chickpea.

Genome-wide identification and analysis of long non-coding RNAs involved in fatty acid biosynthesis in young soybean pods

Long non-coding RNAs (lncRNAs) are non-coding RNAs of more than 200 nucleotides. To date, the roles of lncRNAs in soybean fatty acid synthesis have not been fully studied. Here, the low-linolenic acid mutant 'MT72' and the wild-type control 'JN18' were used as materials. The lncRNAs in young pods at 30 and 40 days (d) after flowering were systematically identified and analyzed using transcriptome sequencing technology combined with bioinformatics tools. A total of 39,324 lncRNAs and 561 differentially expressed lncRNAs were identified. A lncRNAs-miRNAs-protein-coding genes (mRNAs) network was constructed, and 46 lncRNAs, 46 miRNAs and 137 mRNAs were found to be correlated. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of 12 targeted mRNAs in the competing endogenous RNA network showed that these lncRNAs may be involved in the biological processes of fatty acid transport, lipid synthesis and cell division. Finally, the expression levels of differentially expressed lncRNAs, miRNAs and mRNAs were verified using qRT-PCR. The expression patterns of most genes were consistent with the sequencing results. In conclusion, new information was provided for the study of fatty acid synthesis by lncRNAs in young soybean pods.

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.

Pod-shattering characteristics differences between two groups of soybeans are associated with specific changes in gene expression

Soybean is an economically important leguminous crop, and pod dehiscence of soybean could cause huge yield loss. In this study, we measured fruit-cracking forces and percentages of dehisced pods for ten soybean accessions, then separated them into two groups as shattering-sensitive (SS) and shattering-resistant (SR) soybeans. Pod transcriptomes from these two groups were analyzed, and 225 differentially expressed genes (DEGs) were identified between SS and SR soybeans. Some of these DEGs have been previously reported to be associated with pod dehiscence in soybean. The expression patterns of selected DEGs were validated by real-time quantitative reverse transcription PCR, which confirmed the expression changes found in RNA-seq analysis. We also de novo identified 246 soybean pod-long intergenic ncRNAs (lincRNAs), 401 intronic lncRNAs, and 23 antisense lncRNAs from these transcriptomes. Furthermore, genes and lincRNAs co-expression network analysis showed that there are distinct expression patterns between SS and SR soybeans in some co-expression modules. In conclusion, we systematically investigated potential genes and molecular pathways as candidates for differences in soybean pod dehiscence and will provide a useful resource for molecular breeding of soybeans.

Utilization of Interspecific High-Density Genetic Map of RIL Population for the QTL Detection and Candidate Gene Mining for 100-Seed Weight in Soybean

Seed-weight is one of the most important traits determining soybean yield. Hence, it is prerequisite to have detailed understanding of the genetic basis regulating seed-weight for the development of improved cultivars. In this regard, the present study used high-density interspecific linkage map of NJIR4P recombinant inbred population evaluated in four different environments to detect stable Quantitative trait loci (QTLs) as well as mine candidate genes for 100-seed weight. In total, 19 QTLs distributed on 12 chromosomes were identified in all individual environments plus combined environment, out of which seven were novel and eight are stable identified in more than one environment. However, all the novel QTLs were minor (R 2 < 10%). The remaining 12 QTLs detected in this study were co-localized with the earlier reported QTLs with narrow genomic regions, and out of these only 2 QTLs were major (R 2 > 10%) viz., qSW-17-1 and qSW-17-4. Beneficial alleles of all identified QTLs were derived from cultivated soybean parent (Nannong493-1). Based on Protein ANalysis THrough Evolutionary Relationships, gene annotation information, and literature search, 29 genes within 5 stable QTLs were predicted to be possible candidate genes that might regulate seed-weight/size in soybean. However, it needs further validation to confirm their role in seed development. In conclusion, the present study provides better understanding of trait genetics and candidate gene information through the use high-density inter-specific bin map, and also revealed considerable scope for genetic improvement of 100-seed weight in soybean using marker-assisted breeding.

Diversity of a wall-associated kinase gene in wild and cultivated barley

Domestication of barley and other cereals was accompanied by an increase in seed size which has been ascribed to human selection, large seeds being preferred by early farmers or favoured by cultivation practices such as deep sowing. An alternative suggestion is that the increase in seed size was an indirect consequence of selection for plants with more vigorous growth. To begin to address the latter hypothesis we studied the diversity of HvWAK1, a wall-associated kinase gene involved in root proliferation, in 220 wild barley accessions and 200 domesticated landraces. A 3655-bp sequence comprising the gene and upstream region contained 69 single nucleotide polymorphisms (SNPs), one indel and four short tandem repeats. A network of 50 haplotypes revealed a complex evolutionary relationship, but with landraces largely restricted to two parts of the topology. SNPs in the HvWAK1 coding region resulted in nonsynonymous substitutions at nine positions in the translation product, but none of these changes were predicted to have a significant effect on the protein structure. In contrast, the region upstream of the coding sequence contained five SNPs that were invariant in the domesticated population, fixation of these SNPs decreasing the likelihood that the upstream of a pair of TATA boxes and transcription start sites would be used to promote transcription of HvWAK1. The sequence diversity therefore suggests that the cis-regulatory region of HvWAK1 might have been subject to selection during barley domestication. The extent of root proliferation has been linked with traits such as above-ground biomass, so selection for particular cis-regulatory variants of HvWAK1 would be consistent with the hypothesis that seed size increases during domestication were the indirect consequence of selection for plants with increased growth vigour.

Pol III-Dependent Cabbage BoNR8 Long ncRNA Affects Seed Germination and Growth in Arabidopsis

Long non-coding RNAs (lncRNAs) are non-protein-coding transcripts longer than 200 nt that are distributed widely in organisms and play many physiological roles. The BoNR8 lncRNA is a 272 nt long transcript yielded by RNA polymerase III in cabbage that was identified as the closest homolog of the AtR8 lncRNA in Arabidopsis. The BoNR8 lncRNA was expressed extensively in the epidermal tissue in the root elongation zone of germinated seeds, and its accumulation was induced by abiotic stresses, auxins and ABA. To investigate the correlation between the BoNR8 lncRNA and germination, BoNR8-overexpressing Arabidopsis plants (BoNR8-AtOX) were prepared. Three independent BoNR8-AtOX lines showed less primary root elongation, incomplete silique development and decreased germination rates. The germination efficiencies were affected strongly by ABA and slightly by salt stress, and ABA-related gene expression was changed in the BoNR8-AtOX lines.