Molecular regulation of seed development and strategies for engineering seed size in crop plants

Genome-Wide Analyses of MADS-Box Genes Reveal Their Involvement in Seed Development and Oil Accumulation of Tea-Oil Tree (Camellia oleifera)

The seeds of Camellia oleifera produce high amount of oil, which can be broadly used in the fields of food, industry, and medicine. However, the molecular regulation mechanisms of seed development and oil accumulation in C. oleifera are unclear. In this study, evolutionary and expression analyses of the MADS-box gene family were performed across the C. oleifera genome for the first time. A total of 86 MADS-box genes (ColMADS) were identified, including 60 M-type and 26 MIKC members. More gene duplication events occurred in M-type subfamily (6) than that in MIKC subfamily (2), and SEP-like genes were lost from the MIKCC clade. Furthermore, 8, 15, and 17 differentially expressed ColMADS genes (DEGs) were detected between three developmental stages of seed (S1/S2, S2/S3, and S1/S3), respectively. Among these DEGs, the STK-like ColMADS12 and TT16-like ColMADS17 were highly expressed during the seed formation (S1 and S2), agreeing with their predicted functions to positively regulate the seed organogenesis and oil accumulation. While ColMADS57 and ColMADS07 showed increasing expression level with the seed maturation (S2 and S3), conforming to their potential roles in promoting the seed ripening. In all, these results revealed a critical role of MADS-box genes in the C. oleifera seed development and oil accumulation, which will contribute to the future molecular breeding of C. oleifera.

Plant Hormone Pathway Is Involved in Regulating the Embryo Development Mechanism of the Hydrangea macrophylla Hybrid

The research is aimed to elucidate the role of plant hormones in regulating the development of hybrid embryos in Hydrangea macrophylla. Fruits from the intraspecific cross of H. macrophylla 'Otaksa' × 'Coerulea' were selected at the globular, heart, and torpedo stages of embryo development. Transcriptome sequencing and differential gene expression analysis were conducted. The results showed that fruit growth followed a single "S-shaped growth curve, with globular, heart, and torpedo embryos appearing at 30, 40, and 50 d post-pollination, respectively, and the embryo maintaining the torpedo shape from 60 to 90 d. A total of 12,933 genes was quantified across the three developmental stages, with 3359, 3803, and 3106 DEGs in the S1_vs_S2, S1_vs_S3, and S2_vs_S3 comparisons, respectively. Among these, 133 genes related to plant hormone biosynthesis and metabolism were differentially expressed, regulating the synthesis and metabolism of eight types of plant hormones, including cytokinin, auxin, gibberellin, abscisic acid, and jasmonic acid. The pathways with the most differentially expressed genes were cytokinin, auxin, and gibberellin, suggesting these hormones may play crucial roles in embryo development. In the cytokinin pathway, CKX (Hma1.2p1_0579F.1_g182670.gene, Hma1.2p1_1194F.1_g265700.gene, and NewGene_12164) genes were highly expressed during the globular embryo stage, promoting rapid cell division in the embryo. In the auxin pathway, YUC (Hma1.2p1_0271F.1_g109005.gene and Hma1.2p1_0271F.1_g109020.gene) genes were progressively up-regulated during embryo growth; the early response factor AUX/IAA (Hma1.2p1_0760F.1_g214260.gene) was down-regulated, while the later transcriptional activator ARF (NewGene_21460, NewGene_21461, and Hma1.2p1_0209F.1_g089090.gene) was up-regulated, sustaining auxin synthesis and possibly preventing the embryo from transitioning to maturity. In the gibberellin pathway, GA3ox (Hma1.2p1_0129F.1_g060100.gene) expression peaked during the heart embryo stage and then declined, while the negative regulator GA2ox (Hma1.2p1_0020F.1_g013915.gene) showed the opposite trend; and the gibberellin signaling repressor DELLA (Hma1.2p1_1054F.1_g252590.gene) increased over time, potentially inhibiting embryo development and maintaining the torpedo shape until fruit maturity. These findings preliminarily uncover the factors affecting the development of hybrid H. macrophylla embryos, laying a foundation for further research into the regulatory mechanisms of H. macrophylla hybrid embryo development.

MtPIN4 plays critical roles in amino acid biosynthesis and metabolism of seed in Medicago truncatula

The regulation of seed development is critical for determining crop yield. Auxins are vital phytohormones that play roles in various aspects of plant growth and development. However, its role in amino acid biosynthesis and metabolism in seeds is not fully understood. In this study, we identified a mutant with small seeds through forward genetic screening in Medicago truncatula. The mutated gene encodes MtPIN4, an ortholog of PIN1. Using molecular approaches and integrative omics analyses, we discovered that auxin and amino acid content significantly decreased in mtpin4 seeds, highlighting the role of MtPIN4-mediated auxin distribution in amino acid biosynthesis and metabolism. Furthermore, genetic analysis revealed that the three orthologs of PIN1 have specific and overlapping functions in various developmental processes in M. truncatula. Our findings emphasize the significance of MtPIN4 in seed development and offer insights into the molecular mechanisms governing the regulation of seed size in crops. This knowledge could be applied to enhance crop quality by targeted manipulation of seed protein regulatory pathways.

Identification of Candidate Expansin Genes Associated with Seed Weight in Pomegranate (Punica granatum L.)

Seed weight is an important target trait in pomegranate breeding and culture. Expansins act by loosening plant cell walls and cellulosic materials, permitting turgor-driven cell enlargement. However, the role of expansin genes (EXPs) in pomegranate seed weight remains elusive. A total of 29 PgrEXPs were identified in the 'Dabenzi' genome. These genes were classified into four subfamilies and 14 subgroups, including 22 PgrEXPAs, 5 PgrEXPBs, 1 PgrEXPLA, and 1 PgrEXPLB. Transcriptome analysis of PgrEXPs in different tissues (root, leaf, flower, peel, and seed testa) in 'Dabenzi', and the seed testa of the hard-seeded pomegranate cultivar 'Dabenzi' and soft-seeded cultivar 'Tunisia' at three development stages showed that three PgrEXPs (PgrEXPA11, PgrEXPA22, PgrEXPA6) were highly expressed throughout seed development, especially in the sarcotesta. SNP/Indel markers of these PgrEXPs were developed and used to genotype 101 pomegranate accessions. The association of polymorphic PgrEXPs with seed weight-related traits (100-seed weight, 100-kernel weight, 100-sarcotesta weight, and the percentage of 100-sarcotesta to 100-seed weight) were analyzed. PgrEXP22 was significantly associated with 100-seed weight and 100-sarcotesta weight and is a likely candidate for regulating seed weight and sarcotesta development in particular. This study provides an effective tool for the genetic improvement of seed weight in pomegranate breeding programs.

Heat shock factor binding protein BrHSBP1 regulates seed and pod development in Brassica rapa

Plant heat shock factor binding proteins (HSBPs) are well known for their implication in the negative regulation of heat stress response (HSR) pathways. Herein, we report on the hitherto unknown functions of HSBP1 in Brassica rapa (BrHSBP1). BrHBSP1 was found to be predominant in flower buds and young leaves, while its segmental duplicate, BrHSBP1-like, was abundant in green siliques. Exposure to abiotic stress conditions, such as heat, drought, cold, and H2O2, and to phytohormones was found to differentially regulate BrHSBP1. The activity of BrHSBP1-GFP fusion proteins revealed their cellular localization in nuclei and cytosols. Transgenic overexpression of BrHSBP1 (BrHSBP1OX) improved pod and seed sizes, while CRISPR-Cas BrHSBP1 knock-out mutants (Brhsbp1_KO) were associated with aborted seed and pod development. The transcriptomic signatures of BrHSBP1OX and Brhsbp1_KO lines revealed that 360 and 2381 genes, respectively, were differentially expressed (Log2FC≥2, padj<0.05) expressed relative to control lines. In particular, developmental processes, including plant reproductive structure development (RSD)-related genes, were relatively downregulated in Brhsbp1_KO. Furthermore, yeast two-hybrid assays confirmed that BrHSBP1 can physically bind to RSD and other genes. Taking the findings together, it is clear that BrHSBP1 is involved in seed development via the modulation of RSD genes. Our findings represent the addition of a new regulatory player in seed and pod development in B. rapa.

DELLA proteins positively regulate seed size in Arabidopsis

Human and animal nutrition is mainly based on seeds. Seed size is a key factor affecting seed yield and has thus been one of the primary objectives of plant breeders since the domestication of crop plants. Seed size is coordinately regulated by signals of maternal and zygotic tissues that control the growth of the seed coat, endosperm and embryo. Here, we provide previously unreported evidence for the role of DELLA proteins, key repressors of gibberellin responses, in the maternal control of seed size. The gain-of-function della mutant gai-1 produces larger seeds as a result of an increase in the cell number in ovule integuments. This leads to an increase in ovule size and, in turn, to an increase in seed size. Moreover, DELLA activity promotes increased seed size by inducing the transcriptional activation of AINTEGUMENTA, a genetic factor that controls cell proliferation and organ growth, in the ovule integuments of gai-1. Overall, our results indicate that DELLA proteins are involved in the control of seed size and suggest that modulation of the DELLA-dependent pathway could be used to improve crop yield.

Transcriptome Landscape Analyses of the Regulatory Network for Zygotic Embryo Development in Paeonia ostii

Paeonia ostii is a worldwide ornamental flower and an emerging oil crop. Zyotic embryogenesis is a critical process during seed development, and it can provide a basis for improving the efficiency of somatic embryogenesis (SE). In this study, transcriptome sequencing of embryo development was performed to investigate gene expression profiling in P. ostii and identified Differentially expressed genes (DEGs) related to transcription factors, plant hormones, and antioxidant enzymes. The results indicated that IAA (Indole-3-acetic acid), GA (Gibberellin), BR (Brassinosteroid) and ETH (Ethylene) were beneficial to early embryonic morphogenesis, while CTK (Cytokinin) and ABA (Abscisic Acid) promoted embryo morphogenesis and maturation. The antioxidant enzymes' activity was the highest in early embryos and an important participant in embryo formation. The high expression of the genes encoding fatty acid desaturase was beneficial to fast oil accumulation. Representative DEGs were selected and validated using qRT-PCR. Protein-protein interaction network (PPI) was predicted, and six central node proteins, including AUX1, PIN1, ARF6, LAX3, ABCB19, PIF3, and PIF4, were screened. Our results provided new insights into the formation of embryo development and even somatic embryo development in tree peonies.

Construction of a high-density genetic map for faba bean (Vicia faba L.) and quantitative trait loci mapping of seed-related traits

Faba bean (Vicia faba L.) is a valuable legume crop and data on its seed-related traits is required for yield and quality improvements. However, basic research on faba bean is lagging compared to that of other major crops. In this study, an F₂ faba bean population, including 121 plants derived from the cross WY7×TCX7, was genotyped using the Faba_bean_130 K targeted next-generation sequencing genotyping platform. The data were used to construct the first ultra-dense faba bean genetic map consisting of 12,023 single nucleotide polymorphisms markers covering 1,182.65 cM with an average distance of 0.098 cM. The map consisted of 6 linkage groups, which is consistent with the 6 faba bean chromosome pairs. A total of 65 quantitative trait loci (QTL) for seed-related traits were identified (3 for 100-seed weight, 28 for seed shape, 12 for seed coat color, and 22 for nutritional quality). Furthermore, 333 candidate genes that are likely to participate in the regulation of seed-related traits were also identified. Our research findings can provide a basis for future faba bean marker-assisted breeding and be helpful to further modify and improve the reference genome.

Transfer of the high-temperature adult-plant stripe rust resistance gene Yr62 in four Chinese wheat cultivars

Wheat stripe rust is one of the diseases that seriously affect wheat production worldwide. Breeding resistant cultivars is an effective way to control this disease. The wheat stripe rust resistance gene Yr62 has high-temperature adult-plant resistance (HTAP). In this study, PI 660,060, a single Yr62 gene line, was crossed with four Chinese wheat cultivars, LunXuan987 (LX987), Bainongaikang58 (AK58), ZhengMai9023 (ZM9023), and HanMai6172 (H6172). F1 seeds of four cross combinations were planted and self-crossed to develop the advance generations in the field. The seeds of each cross were mixed harvested and about 2400 to 3000 seeds were sown in each generation for F1 to F4 to maintain the maximum possible genotypes. Forty-five lines were selected and evaluated for resistance to stripe rust and agronomic traits, including plant height, number of grains per spike, and tiller number, in F5 and F6. Then, 33 lines with good agronomic traits and high disease resistance were developed to F9 generation. SSR markers Xgwm251 and Xgwm192 flank linked with the Yr62 were used to detect the presence of Yr62 in these 33 F9 lines. Of these, 22 lines were confirmed with the resistance gene Yr62. Finally, nine lines with good agronomic traits and disease resistance were successfully selected. The selected wheat lines in this study provide material support for the future breeding of wheat for stripe rust resistance.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01393-1.

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.

The micronutrient content in underutilized crops: the Lupinus mutabilis sweet case

Adequate intake of micronutrients is necessary to reduce widespread health issues linked to low intake of iron (Fe), zinc (Zn), boron (B), copper (Cu), and manganese (Mn). Because more than two billion people suffer from micronutrient deficiency globally, to address this problem, highly-nutritious ancestral Peruvian crops like tarwi can be an important component of food security. Thus, our work explores the tarwi micronutrient variability to select biofortified genotypes without affecting seed size and weight. Tarwi is a biofortified food because of its seeds' Fe, Zn, and B content. Furthermore, Boron showed a positive correlation between seed size and weight. At the same time, copper showed a negative correlation. Finally, six accessions (P14, P16, P21, T05, T08, and T25) that are biofortified for Fe, Zn, and B with excellent seed size and weight and with adequate levels of Cu and Mn; adding value to Peruvian biodiversity at a low cost is a starting point for a breeding program to prevent micronutrient disorders.

Dynamic transcriptome analysis suggests the key genes regulating seed development and filling in Tartary buckwheat (Fagopyrum tataricum Garetn.)

Tartary buckwheat is highly attractive for the richness of nutrients and quality, yet post-embryonic seed abortion greatly halts the yield. Seed development is crucial for determining grain yield, whereas the molecular basis and regulatory network of Tartary buckwheat seed development and filling is not well understood at present. Here, we assessed the transcriptional dynamics of filling stage Tartary buckwheat seeds at three developmental stages by RNA sequencing. Among the 4249 differentially expressed genes (DEGs), genes related to seed development were identified. Specifically, 88 phytohormone biosynthesis signaling genes, 309 TFs, and 16 expansin genes participating in cell enlargement, 37 structural genes involved in starch biosynthesis represented significant variation and were candidate key seed development genes. Cis-element enrichment analysis indicated that the promoters of differentially expressed expansin genes and starch biosynthesis genes are rich of hormone-responsive (ABA-, AUX-, ET-, and JA-), and seed growth-related (MYB, MYC and WRKY) binding sites. The expansin DEGs showed strong correlations with DEGs in phytohormone pathways and transcription factors (TFs). In total, phytohormone ABA, AUX, ET, BR and CTK, and related TFs could substantially regulate seed development in Tartary buckwheat through targeting downstream expansin genes and structural starch biosynthetic genes. This transcriptome data could provide a theoretical basis for improving yield of Tartary buckwheat.

Epitranscriptomics: An Additional Regulatory Layer in Plants' Development and Stress Response

Epitranscriptomics has added a new layer of regulatory machinery to eukaryotes, and the advancement of sequencing technology has revealed more than 170 post-transcriptional modifications in various types of RNAs, including messenger RNA (mRNA), transfer RNA (tRNA), ribosomal RNA (rRNA), and long non-coding RNA (lncRNA). Among these, N6-methyladenosine (m6A) and N5-methylcytidine (m5C) are the most prevalent internal mRNA modifications. These regulate various aspects of RNA metabolism, mainly mRNA degradation and translation. Recent advances have shown that regulation of RNA fate mediated by these epitranscriptomic marks has pervasive effects on a plant's development and responses to various biotic and abiotic stresses. Recently, it was demonstrated that the removal of human-FTO-mediated m6A from transcripts in transgenic rice and potatoes caused a dramatic increase in their yield, and that the m6A reader protein mediates stress responses in wheat and apple, indicating that regulation of m6A levels could be an efficient strategy for crop improvement. However, changing the overall m6A levels might have unpredictable effects; therefore, the identification of precise m6A levels at a single-base resolution is essential. In this review, we emphasize the roles of epitranscriptomic modifications in modulating molecular, physiological, and stress responses in plants, and provide an outlook on epitranscriptome engineering as a promising tool to ensure food security by editing specific m6A and m5C sites through robust genome-editing technology.

Editing of the starch branching enzyme gene SBE2 generates high-amylose storage roots in cassava

The production of high-amylose cassava through CRISPR/Cas9-mediated mutagenesis of the starch branching enzyme gene SBE2 was firstly achieved. High-amylose cassava (Manihot esculenta Crantz) is desirable for starch industrial applications and production of healthier processed food for human consumption. In this study, we report the production of high-amylose cassava through CRISPR/Cas9-mediated mutagenesis of the starch branching enzyme 2 (SBE2). Mutations in two targeted exons of SBE2 were identified in all regenerated plants; these mutations, which included nucleotide insertions, and short or long deletions in the SBE2 gene, were classified into eight mutant lines. Three mutants, M6, M7 and M8, with long fragment deletions in the second exon of SBE2 showed no accumulation of SBE2 protein. After harvest from the field, significantly higher amylose (up to 56% in apparent amylose content) and resistant starch (up to 35%) was observed in these mutants compared with the wild type, leading to darker blue coloration of starch granules after quick iodine staining and altered starch viscosity with a higher pasting temperature and peak time. Further ¹H-NMR analysis revealed a significant reduction in the degree of starch branching, together with fewer short chains (degree of polymerization [DP] 15-25) and more long chains (DP>25 and especially DP>40) of amylopectin, which indicates that cassava SBE2 catalyzes short chain formation during amylopectin biosynthesis. Transition from A- to B-type crystallinity was also detected in the starches. Our study showed that CRISPR/Cas9-mediated mutagenesis of starch biosynthetic genes in cassava is an effective approach for generating novel varieties with valuable starch properties for food and industrial applications.

Comparative Analysis of Seed Transcriptome and Coexpression Analysis Reveal Candidate Genes for Enhancing Seed Size/Weight in Brassica juncea

Seed size/weight is a multigenic trait that is governed by complex transcriptional regulatory pathways. An understanding of the genetic basis of seed size is of great interest in the improvement of seed yield and quality in oilseed crops. A global transcriptome analysis was performed at the initial stages of seed development in two lines of Brassica juncea, small-seeded EH-2 and large-seeded PJ. The anatomical analyses revealed significant differences in cell number and cell size in the outer layer of the seed coat between EH-2 and PJ. Pairwise comparisons at each developmental stage identified 5,974 differentially expressed genes (DEGs) between the two lines, of which 954 genes belong to different families of transcription factors. Two modules were found to be significantly correlated with an increased seed size using weighted gene coexpression network analysis. The DEG and coexpression datasets were integrated with the thousand seed weight (Tsw) quantitative trait loci (QTL) mapped earlier in the EPJ (EH-2 × PJ) doubled haploid (DH) population, which identified forty potential key components controlling seed size. The candidate genes included genes regulating the cell cycle, cell wall biogenesis/modification, solute/sugar transport, and hormone signaling. The results provide a valuable resource to widen the current understanding of regulatory mechanisms underlying seed size in B. juncea.

Transgenerational Genetic Effects Help Explain Latitudinal Variation in Seed Mass and Germination Timing in Plantago lanceolata

We know little about the underlying genetic control of phenotypic patterns of seed traits across large-scale geographic and environmental gradients. Such knowledge is important for understanding the evolution of populations within species and for improving species conservation. Therefore, to test for genetic variation in Plantago lanceolata, we made reciprocal crosses between northern and southern genotypes that span the species’ range in Europe. The results provide evidence of transgenerational genetic effects on seed mass and germination timing. Northern mothers produced larger seeds with delayed germination, in contrast to southern mothers, which produced smaller seeds with accelerated germination. A maternal latitude affected both the seed coat, solely maternal tissue, and embryo/endosperm tissues. Thus, latitudinal variation in seed size and germination timing can be explained, in part, by the direct influence of maternal genotype, independent of zygotic genes that parents pass directly to the embryo and endosperm. Data suggest that researchers exploring the existence and evolution of large-scale geographic variation within species test for transgenerational genetic effects. In addition, data suggest that transgenerational control of seed traits should be considered when developing procedures designed to facilitate species conservation and restoration.

Transcriptome Analysis of Seed Weight Plasticity in Brassica napus

A critical barrier to improving crop yield is the trade-off between seed weight (SW) and seed number (SN), which has been commonly reported in several crops, including Brassica napus. Despite the agronomic relevance of this issue, the molecular factors involved in the interaction between SW and SN are largely unknown in crops. In this work, we performed a detailed transcriptomic analysis of 48 seed samples obtained from two rapeseed spring genotypes subjected to different source-sink (S-S) ratios in order to examine the relationship between SW and SN under different field conditions. A multifactorial analysis of the RNA-seq data was used to identify a group of 1014 genes exclusively regulated by the S-S ratio. We found that a reduction in the S-S ratio during seed filling induces the expression of genes involved in sucrose transport, seed weight, and stress responses. Moreover, we identified five co-expression modules that are positively correlated with SW and negatively correlated with SN. Interestingly, one of these modules was significantly enriched in transcription factors (TFs). Furthermore, our network analysis predicted several NAC TFs as major hubs underlying SW and SN compensation. Taken together, our study provides novel insights into the molecular factors associated with the SW-SN relationship in rapeseed and identifies TFs as potential targets when improving crop yield.

Genome-wide identification and functional characterization of LEA genes during seed development process in linseed flax (Linum usitatissimum L.)

BACKGROUND

LEA proteins are widely distributed in the plant and animal kingdoms, as well as in micro-organisms. LEA genes make up a large family and function in plant protection against a variety of adverse conditions.

RESULTS

Bioinformatics approaches were adopted to identify LEA genes in the flax genome. In total, we found 50 LEA genes in the genome. We also conducted analyses of the physicochemical parameters and subcellular location of the genes and generated a phylogenetic tree. LuLEA genes were unevenly mapped among 15 flax chromosomes and 90% of the genes had less than two introns. Expression profiles of LuLEA showed that most LuLEA genes were expressed at a late stage of seed development. Functionally, the LuLEA1 gene reduced seed size and fatty acid contents in LuLEA1-overexpressed transgenic Arabidopsis lines.

CONCLUSION

Our study adds valuable knowledge about LEA genes in flax which can be used to improve related genes of seed development.

Integrated microRNA and transcriptome profiling reveal key miRNA-mRNA interaction pairs associated with seed development in Tartary buckwheat (Fagopyrum tataricum)

BACKGROUND

Tartary buckwheat seed development is an extremely complex process involving many gene regulatory pathways. MicroRNAs (miRNAs) have been identified as the important negative regulators of gene expression and performed crucial regulatory roles in various plant biological processes. However, whether miRNAs participate in Tartary buckwheat seed development remains unexplored.

RESULTS

In this study, we first identified 26 miRNA biosynthesis genes in the Tartary buckwheat genome and described their phylogeny and expression profiling. Then we performed small RNA (sRNA) sequencing for Tartary buckwheat seeds at three developmental stages to identify the miRNAs associated with seed development. In total, 230 miRNAs, including 101 conserved and 129 novel miRNAs, were first identified in Tartary buckwheat, and 3268 target genes were successfully predicted. Among these miRNAs, 76 exhibited differential expression during seed development, and 1534 target genes which correspond to 74 differentially expressed miRNAs (DEMs) were identified. Based on integrated analysis of DEMs and their targets expression, 65 miRNA-mRNA interaction pairs (25 DEMs corresponding to 65 target genes) were identified that exhibited significantly opposite expression during Tartary buckwheat seed development, and 6 of the miRNA-mRNA pairs were further verified by quantitative real-time polymerase chain reaction (qRT-PCR) and ligase-mediated rapid amplification of 5' cDNA ends (5'-RLM-RACE). Functional annotation of the 65 target mRNAs showed that 56 miRNA-mRNA interaction pairs major involved in cell differentiation and proliferation, cell elongation, hormones response, organogenesis, embryo and endosperm development, seed size, mineral elements transport, and flavonoid biosynthesis, which indicated that they are the key miRNA-mRNA pairs for Tartary buckwheat seed development.

CONCLUSIONS

Our findings provided insights for the first time into miRNA-mediated regulatory pathways in Tartary buckwheat seed development and suggested that miRNAs play important role in Tartary buckwheat seed development. These findings will be help to study the roles and regulatory mechanism of miRNAs in Tartary buckwheat seed development.

A proteomic analysis of peanut seed at different stages of underground development to understand the changes of seed proteins

In order to obtain more valuable insights into the protein dynamics and accumulation of allergens in seeds during underground development, we performed a proteomic study on developing peanut seeds at seven different stages. A total of 264 proteins with altered abundance and contained at least one unique peptide was detected by matrix-assisted laser desorption ionization time-of-flight/time-of-flight mass spectrometry (MALDI-TOF/TOF MS). All identified proteins were classified into five functional categories as level 1 and 20 secondary functional categories as level 2. Among them, 88 identified proteins (IPs) were related to carbohydrate/ amino acid/ lipid transport and metabolism, indicating that carbohydrate/amino acid/ lipid metabolism played a key role in the underground development of peanut seeds. Hierarchical cluster analysis showed that all IPs could be classified into eight cluster groups according to the abundance profiles, suggesting that the modulatory patterns of these identified proteins were complicated during seed development. The largest group contained 41 IPs, the expression of which decreased at R 2 and reached a maximum at R3 but gradually decreased from R4. A total of 14 IPs were identified as allergen-like proteins by BLAST with A genome (Arachis duranensis) or B genome (Arachis ipaensis) translated allergen sequences. Abundance profile analysis of 14 identified allergens showed that the expression of all allergen proteins was low or undetectable by 2-DE at the early stages (R1 to R4), and began to accumulate from the R5 stage and gradually increased. Network analysis showed that most of the significant proteins were involved in active metabolic pathways in early development. Real time RT-PCR analysis revealed that transcriptional regulation was approximately consistent with expression at the protein level for 8 selected identified proteins. In addition, some amino acid sequences that may be associated with new allergens were also discussed.

Heterotrimeric Gα subunit regulates plant architecture, organ size and seed weight in the oilseed Brassica juncea

Two BjuGα proteins exhibit conserved GTP-binding and GTP-hydrolysis activities, and function in maintaining overall plant architecture and controlling multiple yield-related traits in the oilseed Brassica juncea. Heterotrimeric G-protein (Gα, Gβ and Gγ) are key signal transducers, well characterized in model plants Arabidopsis and rice. However, our knowledge about the roles played by G-proteins in regulating various growth and developmental traits in polyploid crops, having a complex G-protein signalling network, is quite sparse. In the present study, two Gα encoding genes (BjuA.Gα1 and BjuB.Gα1) were isolated from the allotetraploid Brassica juncea, a globally cultivated oilseed crop of the Brassicaceae family. BjuGα1 genes share a close evolutionary relationship, and the encoded proteins exhibit highly conserved G-protein activities while showing expression differentiation, wherein BjuA.Gα1 was the highly abundant transcript during plant growth and developmental stages. RNAi based suppression of BjuGα1 displayed compromised effects on most of the tested vegetative and reproductive parameters, particularly plant height (32-58%), flower and siliques dimensions, and seed weight (11-13%). Further, over-expression of a constitutively active Gα, lacking the GTPase activity, produced plants with increased height, organ size and seed weight (7-25%), without altering seed quality traits like fatty acid composition, glucosinolates, oil and protein contents. Our study demonstrates that BjuGα1 proteins control overall plant architecture and multiple yield-related traits in the oilseed B. juncea, suggesting that BjuGα1 could be a promising target for crop improvement.

Quantitative Trait Loci for Seed Size Variation in Cucurbits - A Review

Cucurbits (Cucurbitaceae family) include many economically important fruit vegetable crops such as watermelon, pumpkin/squash, cucumber, and melon. Seed size (SS) is an important trait in cucurbits breeding, which is controlled by quantitative trait loci (QTL). Recent advances have deciphered several signaling pathways underlying seed size variation in model plants such as Arabidopsis and rice, but little is known on the genetic basis of SS variation in cucurbits. Here we conducted literature review on seed size QTL identified in watermelon, pumpkin/squash, cucumber and melon, and inferred 14, 9 and 13 consensus SS QTL based on their physical positions in respective draft genomes. Among them, four from watermelon (ClSS2.2, ClSS6.1, ClSS6.2, and ClSS8.2), two from cucumber (CsSS4.1 and CsSS5.1), and one from melon (CmSS11.1) were major-effect, stable QTL for seed size and weight. Whole genome sequence alignment revealed that these major-effect QTL were located in syntenic regions across different genomes suggesting possible structural and functional conservation of some important genes for seed size control in cucurbit crops. Annotation of genes in the four watermelon consensus SS QTL regions identified genes that are known to play important roles in seed size control including members of the zinc finger protein and the E3 ubiquitin-protein ligase families. The present work highlights the utility of comparative analysis in understanding the genetic basis of seed size variation, which may help future mapping and cloning of seed size QTL in cucurbits.

Transcriptome Analysis Reveals Key Seed-Development Genes in Common Buckwheat (Fagopyrum esculentum)

Seed development is an essential and complex process, which is involved in seed size change and various nutrients accumulation, and determines crop yield and quality. Common buckwheat (Fagopyrum esculentum Moench) is a widely cultivated minor crop with excellent economic and nutritional value in temperate zones. However, little is known about the molecular mechanisms of seed development in common buckwheat (Fagopyrum esculentum). In this study, we performed RNA-Seq to investigate the transcriptional dynamics and identify the key genes involved in common buckwheat seed development at three different developmental stages. A total of 4619 differentially expressed genes (DEGs) were identified. Based on the results of Gene Ontology (GO) and KEGG analysis of DEGs, many key genes involved in the seed development, including the Ca²⁺ signal transduction pathway, the hormone signal transduction pathways, transcription factors (TFs), and starch biosynthesis-related genes, were identified. More importantly, 18 DEGs were identified as the key candidate genes for seed size through homologous query using the known seed size-related genes from different seed plants. Furthermore, 15 DEGs from these identified as the key genes of seed development were selected to confirm the validity of the data by using quantitative real-time PCR (qRT-PCR), and the results show high consistency with the RNA-Seq results. Taken together, our results revealed the underlying molecular mechanisms of common buckwheat seed development and could provide valuable information for further studies, especially for common buckwheat seed improvement.

De novo transcriptome assembly and identification of genes related to seed size in common buckwheat (Fagopyrum esculentum M.)

Common buckwheat (Fagopyrum esculentum M.) belongs to the eudicot family Polygonaceae, Fagopyrum Mill, and its seeds have high nutritional value. The mechanism of seed development of common buckwheat remains unclear at the molecular level and no genes related to seed size have been identified. In this study, we performed genome-wide transcriptome sequencing and analysis using common buckwheat seeds at 5 days post anthesis (DPA) and 10 DPA from two cultivars (large-seeded and small-seeded). A total of 259,895 transcripts were assembled, resulting in 187,034 unigenes with average length of 1097 bp and N50 of 1538 bp. Based on gene expression profiles, 9127 differentially expressed genes (DEGs) were identified and analyzed in GO enrichment and KEGG analysis. In addition, genes related to seed size in the IKU pathway, ubiquitin-proteasome pathway, MAPK signaling pathway, TFs and phytohormones were identified and analyzed. AP2 and bZIP transcription factors, BR-signal and ABA were considered to be important regulators of seed size. This study provides a valuable genetic resource for future identification and functional analysis of candidate genes regulating seed size in common buckwheat and will be useful for improving seed yield in common buckwheat through molecular breeding in the future.

The genome assembly and annotation of yellowhorn (Xanthoceras sorbifolium Bunge)

BACKGROUND

Yellowhorn (Xanthoceras sorbifolium Bunge), a deciduous shrub or small tree native to north China, is of great economic value. Seeds of yellowhorn are rich in oil containing unsaturated long-chain fatty acids that have been used for producing edible oil and nervonic acid capsules. However, the lack of a high-quality genome sequence hampers the understanding of its evolution and gene functions.

FINDINGS

In this study, a whole genome of yellowhorn was sequenced and assembled by integration of Illumina sequencing, Pacific Biosciences single-molecule real-time sequencing, 10X Genomics linked reads, Bionano optical maps, and Hi-C. The yellowhorn genome assembly was 439.97 Mb, which comprised 15 pseudo-chromosomes covering 95.42% (419.84 Mb) of the assembled genome. The repetitive fractions accounted for 56.39% of the yellowhorn genome. The genome contained 21,059 protein-coding genes. Of them, 18,503 (87.86%) genes were found to be functionally annotated with ≥1 "annotation" term by searching against other databases. Transcriptomic analysis showed that 341, 135, 125, 113, and 100 genes were specifically expressed in hermaphrodite flower, staminate flower, young fruit, leaf, and shoot, respectively. Phylogenetic analysis suggested that yellowhorn and Dimocarpus longan diverged from their most recent common ancestor ∼46 million years ago.

CONCLUSIONS

The availability and subsequent annotation of the yellowhorn genome, as well as the identification of tissue-specific functional genes, provides a valuable reference for plant comparative genomics, evolutionary studies, and molecular design breeding.