Multi-enzymatic systems synergize new RCA technique amplified super-long dsDNA from DNA circle.

BACKGROUND: In the field of DNA amplification, there are great challenges in the effectively amplify of long-chain amplification, especially amplification up to several hundred kb level. RESULTS: A novel technique for the unbiased whole genome amplification from a thimbleful of DNA circles, such as low as 10 ng/ 10 µL of the circular cpDNA or low as 5 ng/ 10 µL of the plasmid, is developed, which can amplify an abundance of the whole genome sequences. Specifically, the new technique that combines rolling-amplification and triple-enzyme system presents a tightly controlled process of a series of buffers/reactions and optimized procedures, that applies from the primer-template duplexes to the Elution step. The result of this technique provides a new approach for extending RCA capacity, where it can reach 200 kb from the circular cpDNA amplification and 150 kb from the plasmid DNA amplification, that demonstrates superior breadth and evenness of genome coverage, high reproducibility, small amplification bias with the amplification efficiency. SIGNIFICANCE AND NOVELTY: This new technique will develop into one of the powerful tools for isothermal DNA amplification in vitro, genome sequencing/analysis, phylogenetic analysis, physical mapping, and other molecular biology applications.

A novel in-situ-process technique constructs whole circular cpDNA library.

BACKGROUND: The chloroplast genome (cp genome) is directly related to the study and analysis of molecular phylogeny and evolution of plants in the phylogenomics era. The cp genome, whereas, is highly plastic and exists as a heterogeneous mixture of sizes and physical conformations. It is advantageous to purify/enrich the circular chloroplast DNA (cpDNA) to reduce sequence complexity in cp genome research. Large-insert, ordered DNA libraries are more practical for genomics research than conventional, unordered ones. From this, a technique of constructing the ordered BAC library with the goal-insert cpDNA fragment is developed in this paper. RESULTS: This novel in-situ-process technique will efficiently extract circular cpDNA from crops and construct a high-quality cpDNA library. The protocol combines the in-situ chloroplast lysis for the high-purity circular cpDNA with the in-situ substitute/ligation for the high-quality cpDNA library. Individually, a series of original buffers/solutions and optimized procedures for chloroplast lysis in-situ is different than bacterial lysis in-situ; the in-situ substitute/ligation that reacts on the MCE membrane is suitable for constructing the goal-insert, ordered cpDNA library while preventing the large-insert cpDNA fragment breakage. The goal-insert, ordered cpDNA library is arrayed on the microtiter plate by three colonies with the definite cpDNA fragment that are homologous-corresponds to the whole circular cpDNA of the chloroplast. CONCLUSION: The novel in-situ-process technique amply furtherance of research in genome-wide functional analysis and characterization of chloroplasts, such as genome sequencing, bioinformatics analysis, cloning, physical mapping, molecular phylogeny and evolution.

A small heat shock protein GmHSP18.5a improves the male fertility restorability of cytoplasmic male sterility-based restorer line under high temperature stress in soybean.

Small heat shock protein (sHSP) is involved in high temperature (HT) stress response. However, the function of sHSPs in regulating male fertility of soybean under HT stress remains largely unknown. Here, we identified a sHSP gene, GmHSP18.5a, which was responded to HT stress during flowering in cytoplasmic male sterility (CMS)-based restorer line of soybean. Moreover, GmHSFA6b turned out to directly activated the expression of GmHSP18.5a by binding to the heat shock cis-element in its promoter. Overexpression of GmHSP18.5a increased male fertility in transgenic Arabidopsis, soybean CMS-based restorer line and its hybrid F1 with CMS line under HT stress. Reactive oxygen species (ROS) content detection revealed that GmHSP18.5a promoted the ROS scavenging ability of Arabidopsis inflorescence and soybean flower bud under HT stress. Enzyme activity assay and gene expression analysis indicated that GmHS18.5a mainly increased the activity of antioxidant enzymes and the expression level of ROS metabolism-related genes under HT stress. Our results indicated that GmHSP18.5a improved the male fertility restorability of CMS-based restorer line in soybean by regulating ROS metabolic pathway and reducing ROS accumulation. Our findings not only revealed the molecular mechanism of sHSP regulating the male fertility of soybean under HT stress, but also provided a theoretical basis for creating strong restorer line with thermotolerance.

The miR156b-GmSPL2b module mediates male fertility regulation of cytoplasmic male sterility-based restorer line under high-temperature stress in soybean.

High-temperature (HT) stress at flowering stage causes significant damage to soybean, including pollen abortion and fertilization failure, but few genes involved in male fertility regulation under HT stress in soybean have been characterized. Here, we demonstrated that miR156b-GmSPL2b module involved in male fertility regulation of soybean cytoplasmic male sterility (CMS)-based restorer line under HT stress. Overexpression of miR156b decreased male fertility in soybean CMS-based restorer line and its hybrid F1 with CMS line under HT stress. RNA-seq analysis found that miR156b mediated male fertility regulation in soybean under HT stress by regulating the expression of pollen development and HT response related genes. Metabolomic analysis of miR156bOE revealed reduction in flavonoid content under HT stress. Integrated transcriptomic and metabolomic analysis showed that the overexpression of miR156b caused flavonoid metabolism disorder in soybean flower bud under HT stress. Knockout of GmSPL2b also decreased the thermotolerance of soybean CMS-based restorer line during flowering. Moreover, GmSPL2b turned out to be directly bounded to the promoter of GmHSFA6b. Further verification indicated that GmHSFA6b overexpression enhanced HT tolerance in Arabidopsis during flowering. Substance content and gene expression analysis revealed that miR156b-GmSPL2b may mediate reactive oxygen species clearance by regulating flavonoid metabolism, thus participating in the regulation of male fertility in soybean under HT stress. This study not only provided important progress for understanding the molecular mechanism of miR156b-GmSPL2b regulating the male fertility of soybean CMS-based restorer line under HT stress, but also provided genetic resources and theoretical basis for creating HT-tolerant strong restorer lines.

Transcriptome Analysis Reveals the Genes Related to Pollen Abortion in a Cytoplasmic Male-Sterile Soybean (Glycine max (L.) Merr.).

Cytoplasmic male sterility (CMS) lays a foundation for the utilization of heterosis in soybean. The soybean CMS line SXCMS5A is an excellent CMS line exhibiting 100% male sterility. Cytological analysis revealed that in SXCMS5A compared to its maintainer SXCMS5B, its tapetum was vacuolated and abnormally developed. To identify the genes and metabolic pathways involving in pollen abortion of SXCMS5A, a comparative transcriptome analysis was conducted between SXCMS5A and SXCMS5B using flower buds. A total of 372,973,796 high quality clean reads were obtained from 6 samples (3 replicates for each material), and 840 differentially expressed genes (DEGs) were identified, including 658 downregulated and 182 upregulated ones in SXCMS5A compared to SXCMS5B. Among them, 13 DEGs, i.e., 12 open reading frames (ORFs) and 1 COX2, were mitochondrial genome genes in which ORF178 and ORF103c were upregulated in CMS lines and had transmembrane domain(s), therefore, identified as CMS candidate mitochondrial genes of SXCMS5A. Furthermore, numerous DEGs were associated with pollen wall development, carbohydrate metabolism, sugar transport, reactive oxygen species (ROS) metabolism and transcription factor. Some of them were further confirmed by quantitative real time PCR analysis between CMS lines with the same cytoplasmic source as SXCMS5A and their respective maintainer lines. The amount of soluble sugar and adenosine triphosphate and the activity of catalase and ascorbic acid oxidase showed that energy supply and ROS scavenging decreased in SXCMS5A compared to SXCMS5B. These findings provide valuable information for further understanding the molecular mechanism regulating the pollen abortion of soybean CMS.

Gm6PGDH1, a Cytosolic 6-Phosphogluconate Dehydrogenase, Enhanced Tolerance to Phosphate Starvation by Improving Root System Development and Modifying the Antioxidant System in Soybean.

Phosphorus plays an important role in plant growth and development, and is an important limiting factor for crop yield. Although previous studies have shown that 6-phosphogluconate dehydrogenase (6PGDH) plays an important role in plant resistance to adversity, its response to low phosphorus (P) stress remains unknown. In this study, we reported the cloning and characterization of a cytosolic 6PGDH gene, Gm6PGDH1, which enhanced the tolerance to phosphate (Pi) starvation by improving root system development and modifying the antioxidant system in transgenic plants. Gm6PGDH1 was highly expressed in the root at full bloom stage, and strongly induced by Pi starvation. The results from intact soybean composite plant and soybean plant, both containing a Gm6PGDH1-overexpressing construct, showed that Gm6PGDH1 was involved in root system development, and subsequently affected P uptake under Pi-deficient conditions. Meanwhile, the accumulation of reactive oxygen species (ROS) in the root tip of transgenic soybean was reduced, and the activity of ROS-scavenging enzymes was enhanced compared with those of the wild type under Pi-deficient conditions. Interestingly, we found that the overexpression of Gm6PGDH1 weakened the response of several other important Pi-answer genes to Pi starvation, such as some purple acid phosphatases (PAPs) and redox-related genes. In addition, the results from a virus-induced gene silencing (VIGS) indicated that Gm6PGDH1 might have functional redundancy in soybean, and the results from a heterogeneous transformation system showed that overexpressing Gm6PGDH1 also enhanced tolerance to Pi starvation in transgenic Arabidopsis. Together, these results suggested the great potential of Gm6PGDH1 in crop breeding for low Pi tolerance.

A small RNA of miR2119b from soybean CMS line acts as a negative regulator of male fertility in transgenic Arabidopsis.

The miR2119 is involved in the growth, development and abiotic stress response of some legumes, including Medicago truncatula, Phaseolus vulgaris and soybean (Glycine max (L.) Merr.). Our previous small RNA sequencing analysis showed that miR2119b was up-regulated in the flower buds of soybean cytoplasmic male sterile (CMS) line compared with its maintainer line, but the role and mechanism of miR2119b in the regulation of soybean male fertility are still unclear. In this study, the gma-miR2119b and its target gene alcohol dehydrogenase 1.3b (ADH1.3b) were characterized and found to be highly expressed in the flowers of soybean CMS line and its maintainer. Transgenic Arabidopsis plants overexpressing gma-miR2119b exhibit male fertility abnormalities, including pollen fertility and germination rate decreased. Enzyme activity detection found the ADH and catalase (CAT) enzyme activities in inflorescence of gma-miR2119b overexpressed plants were lower than those of wild-type. Bioinformatics and gene expression analysis showed that gma-miR2119b/GmADH1.3b module was responsive to high temperature (HT) stress during flowering. After HT stress, the gma-miR2119b overexpressed plants showed male sterility, including shorter filament, sterile pollen, indehiscent anther and non seed. Moreover, some key genes involved in HT response and reactive oxygen species (ROS) signal regulation pathway, including heat shock protein70, galactinol synthase 1 and CAT, showed down-regulated expression in transgenic plants under HT stress, suggesting that gma-miR2119b regulates male fertility via HT-ROS signaling pathway under HT stress. It was speculated that the gma-miR2119b acted as a negative regulator of male fertility in plants by regulating ADH1, HT-induced and ROS scavenging genes expression.

Confirmation of GmPPR576 as a fertility restorer gene of cytoplasmic male sterility in soybean.

In soybean, heterosis achieved through the three-line system has been gradually applied in breeding to increase yield, but the underlying molecular mechanism remains unknown. We conducted a genetic analysis using the pollen fertility of offspring of the cross NJCMS1A×NJCMS1C. All the pollen of F1 plants was semi-sterile; in F2, the ratio of pollen-fertile plants to pollen-semi-sterile plants was 208:189. This result indicates that NJCMS1A is gametophyte sterile, and the fertility restoration of NJCMS1C to NJCMS1A is a quality trait controlled by a single gene locus. Using bulked segregant analysis, the fertility restorer gene Rf in NJCMS1C was located on chromosome 16 between the markers BARCSOYSSR_16_1067 and BARCSOYSSR_16_1078. Sequence analysis of genes in that region showed that GmPPR576 was non-functional in rf cultivars. GmPPR576 has one functional allele in Rf cultivars but three non-functional alleles in rf cultivars. Phylogenetic analysis showed that the GmPPR576 locus evolved rapidly with the presence of male-sterile cytoplasm. GmPPR576 belongs to the RFL fertility restorer gene family and is targeted to the mitochondria. GmPPR576 was knocked out in soybean N8855 using CRISPR/Cas9. The T1 plants showed sterile pollen, and T2 plants produced few pods at maturity. The results indicate that GmPPR576 is the fertility restorer gene of NJCMS1A.

Heat-Responsive miRNAs Participate in the Regulation of Male Fertility Stability in Soybean CMS-Based F1 under High Temperature Stress.

MicroRNAs (miRNAs), a class of noncoding small RNAs (sRNAs), are widely involved in the response to high temperature (HT) stress at both the seedling and flowering stages. To dissect the roles of miRNAs in regulating male fertility in soybean cytoplasmic male sterility (CMS)-based F1 under HT, sRNA sequencing was performed using flower buds from HT-tolerant and HT-sensitive CMS-based F1 combinations (NF1 and YF1, respectively). A total of 554 known miRNAs, 59 new members of known miRNAs, 712 novel miRNAs, and 1145 target genes of 580 differentially expressed miRNAs (DEMs) were identified under normal temperature and HT conditions. Further integrated analysis of sRNA and transcriptome sequencing found that 21 DEMs and 15 differentially expressed target genes, such as gma-miR397a/Laccase 2, gma-miR399a/Inorganic phosphate transporter 1-4, and gma-miR4413a/PPR proteins, mitochondrial-like, were negatively regulated under HT stress. Furthermore, all members of the gma-miR156 family were suppressed by HT stress in both NF1 and YF1, but were highly expressed in YF1 under HT condition. The negative correlation between gma-miR156b and its target gene squamosa promoter-binding protein-like 2b was confirmed by expression analysis, and overexpression of gma-miR156b in Arabidopsis led to male sterility under HT stress. With these results, we proposed that miRNAs play an important role in the regulation of male fertility stability in soybean CMS-based F1 under HT stress.

Comparative analysis of mitochondrial genomes of soybean cytoplasmic male-sterile lines and their maintainer lines.

In soybean, only one mitochondrial genome of cultispecies has been completely obtained. To explore the effect of mitochondrial genome on soybean cytoplasmic male sterility (CMS), two CMS lines and three maintainer lines were used for sequencing. Comparative analysis showed that mitochondrial genome of the CMS line was more compact than that of its maintainer line, but genes were highly conserved. Conserved and unique sequence coexisted in the genomes. Mitochondrial genomes contained different sequence lengths and copy numbers of repeats between CMS line and maintainer line. Large and short repeats mediated intramolecular and intermolecular recombination in mitochondria. Unique sequences and genes were also involved in recombination process and constituted a complex network. orf178 and orf261 were identified as CMS-associated candidate genes. They had sequence characteristics of reported CMS genes in other crops and could be transcribed in CMS lines but not in maintainer lines. This report reveals mitochondrial genome of soybean CMS lines and compares complete mitochondrial sequence between CMS lines and their maintainer lines. The information will be helpful in further understanding the characteristics of soybean mitochondrial genome and the mechanism underlying CMS.

Genome-wide identification and characterization of GRAS genes in soybean (Glycine max).

BACKGROUND: GRAS proteins are crucial transcription factors, which are plant-specific and participate in various plant biological processes. Thanks to the rapid progress of the whole genome sequencing technologies, the GRAS gene families in different plants have been broadly explored and studied. However, comprehensive research on the soybean (Glycine max) GRAS gene family is relatively lagging. RESULTS: In this study, 117 Glycine max GRAS genes (GmGRAS) were identified. Further phylogenetic analyses showed that the GmGRAS genes could be categorized into nine gene subfamilies: DELLA, HAM, LAS, LISCL, PAT1, SCL3, SCL4/7, SCR and SHR. Gene structure analyses turned out that the GmGRAS genes lacked introns and were relatively conserved. Conserved domains and motif patterns of the GmGRAS members in the same subfamily or clade exhibited similarities. Notably, the expansion of the GmGRAS gene family was driven both by gene tandem and segmental duplication events. Whereas, segmental duplications took the major role in generating new GmGRAS genes. Moreover, the synteny and evolutionary constraints analyses of the GRAS proteins among soybean and distinct species (two monocots and four dicots) provided more detailed evidence for GmGRAS gene evolution. Cis-element analyses indicated that the GmGRAS genes may be responsive to diverse environmental stresses and regulate distinct biological processes. Besides, the expression patterns of the GmGRAS genes were varied in various tissues, during saline and dehydration stresses and during seed germination processes. CONCLUSIONS: We conducted a systematic investigation of the GRAS genes in soybean, which may be valuable in paving the way for future GmGRAS gene studies and soybean breeding.

Comparative Transcriptomics Analysis and Functional Study Reveal Important Role of High-Temperature Stress Response Gene GmHSFA2 During Flower Bud Development of CMS-Based F1 in Soybean.

High-temperature (HT) is one of the most important environmental factors that negatively impact the yield of some soybean cytoplasmic male sterility (CMS)-based hybrid (F1) combinations. The response of soybean to HT, especially at the male organ development stage, is poorly understood. To investigate the molecular mechanisms of the response from soybean CMS-based F1 male organ to HT, a detailed transcriptomics analysis was performed during flower bud development of soybean HT-tolerant and HT-sensitive CMS-based F1 combinations (NF1 and YF1) under normal-temperature and HT conditions. Obvious HT damage was observed by subjecting YF1 with HT, such as indehiscent anthers and decreased pollen fertility, whereas the male fertility of NF1 was normal. In total, 8,784 differentially expressed genes (DEGs) were found to respond to HT stress, which were mainly associated with anther/pollen wall development, carbohydrate metabolism and sugar transport, and auxin signaling. The quantitative real-time PCR (qRT-PCR) analysis and substance content detection also revealed that HT caused male fertility defects in YF1 by altering pectin metabolism, auxin, and sugar signaling pathways. Most importantly, the sugar signaling-PIF-auxin signaling pathway may underlie the instability of male fertility in YF1 under HT. Furthermore, HT induced the expression of heat shock factor (HSF) and heat shock protein (HSP) gene families. Overexpression of GmHSFA2 in Arabidopsis can promote the expression of HT protective genes (such as HSP20) by binding to the HSE motifs in their promoters, so as to improve the HT tolerance during flowering. Our results indicated that GmHSFA2 acted as a positive regulator, conferring HT tolerance improvement in soybean CMS-based F1. GmHSFA2 may be directly involved in the activation of male fertility protection mechanism in the soybean CMS-based F1 under HT stress.

Metabolomics Studies on Cytoplasmic Male Sterility during Flower Bud Development in Soybean.

Abnormal reactive oxygen species (ROS) may mediate cytoplasmic male sterility (CMS). To observe the effect of ROS on soybean CMS, metabolite content and antioxidant enzyme activity in the flower buds between soybean N8855-derived CMS line and its maintainer were compared. Of the 612 metabolites identified, a total of 74 metabolites were significantly differentiated in flower buds between CMS line and its maintainer. The differential metabolites involved 32 differential flavonoids, 13 differential phenolamides, and 1 differential oxidized glutathione (GSSG) belonging to a non-enzymatic ROS scavenging system. We observed lower levels of flavonoids and antioxidant enzyme activities in flower buds of the CMS line than in its maintainer. Our results suggest that deficiencies of enzymatic and non-enzymatic ROS scavenging systems in soybean CMS line cannot eliminate ROS in anthers effectively, excessive accumulation of ROS triggered programmed cell death and ultimately resulted in pollen abortion of soybean CMS line.

Exploration of miRNA-mediated fertility regulation network of cytoplasmic male sterility during flower bud development in soybean.

Cytoplasmic male sterility (CMS) plays an important role in the production of soybean hybrid seeds. MicroRNAs (miRNAs) are a class of non-coding endogenous ~ 21 nt small RNAs that play crucial roles in flower and pollen development by targeting genes in plants. To dissect the function of miRNAs in soybean CMS, a total of 558 known miRNAs, 10 novel miRNAs, and 466 target genes were identified in flower buds of the soybean CMS line NJCMS1A and its restorer line NJCMS1C through small RNA sequencing and degradome analysis. In addition, miRNA-mediated editing events were also observed, and the two most frequently observed editing types (A → G and C → U) were validated by cloning and sequencing. And as the base editing occurred, some targets were filtered, such as gma-miR2118b-P6GT with Glyma.08G122000.2. Further integrated analysis of transcriptome and small RNA found some miRNAs and their targets' expression patterns showing a negative correlation, such as gma-miR156b/GmSPL9a and gma-miR4413b/GmPPR. Furthermore, opposite expression pattern was observed between gma-miR156b and GmSPL9 during early stage of flower bud development. Taken together, the regulatory network of gma-miR156b/GmSPL9 and gma-miR4413b/GmPPR with flower bud development in soybean CMS was developed. Most importantly, previous reports showed that these targets might be related to pollen development and male sterility, suggesting that both conserved and species-specific miRNAs might act as regulators of flower bud development in soybean CMS. These findings may provide a better understanding of the miRNA-mediated regulatory networks of CMS mechanisms in soybean.

Comparative analysis of circular RNAs between soybean cytoplasmic male-sterile line NJCMS1A and its maintainer NJCMS1B by high-throughput sequencing.

BACKGROUND: Cytoplasmic male sterility (CMS) is a natural phenomenon of pollen abortion caused by the interaction between cytoplasmic genes and nuclear genes. CMS is a simple and effective pollination control system, and plays an important role in crop heterosis utilization. Circular RNAs (circRNAs) are a vital type of non-coding RNAs, which play crucial roles in microRNAs (miRNAs) function and post-transcription control. To explore the expression profile and possible functions of circRNAs in the soybean CMS line NJCMS1A and its maintainer NJCMS1B, high-throughput deep sequencing coupled with RNase R enrichment strategy was conducted. RESULTS: CircRNA libraries were constructed from flower buds of NJCMS1A and its maintainer NJCMS1B with three biological replicates. A total of 2867 circRNAs were identified, with 1009 circRNAs differentially expressed between NJCMS1A and NJCMS1B based on analysis of high-throughput sequencing. Of the 12 randomly selected circRNAs with different expression levels, 10 showed consistent expression patterns based on high-throughput sequencing and quantitative real-time PCR analyses. Tissue specific expression patterns were also verified with two circRNAs by quantitative real-time PCR. Most parental genes of differentially expressed circRNAs were mainly involved in biological processes such as metabolic process, biological regulation, and reproductive process. Moreover, 83 miRNAs were predicted from the differentially expressed circRNAs, some of which were strongly related to pollen development and male fertility; The functions of miRNA targets were analyzed using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG), and the target mRNAs were significantly enriched in signal transduction and programmed cell death. Furthermore, a total of 165 soybean circRNAs were predicted to contain at least one internal ribosome entry site (IRES) element and an open reading frame, indicating their potential to encode polypeptides or proteins. CONCLUSIONS: Our study indicated that the circRNAs might participate in the regulation of flower and pollen development, which could provide a new insight into the molecular mechanisms of CMS in soybean.

Genome-wide comparative analysis of DNA methylation between soybean cytoplasmic male-sterile line NJCMS5A and its maintainer NJCMS5B.

BACKGROUND: DNA methylation is an important epigenetic modification. It can regulate the expression of many key genes without changing the primary structure of the genomic DNA, and plays a vital role in the growth and development of the organism. The genome-wide DNA methylation profile of the cytoplasmic male sterile (CMS) line in soybean has not been reported so far. RESULTS: In this study, genome-wide comparative analysis of DNA methylation between soybean CMS line NJCMS5A and its maintainer NJCMS5B was conducted by whole-genome bisulfite sequencing. The results showed 3527 differentially methylated regions (DMRs) and 485 differentially methylated genes (DMGs), including 353 high-credible methylated genes, 56 methylated genes coding unknown protein and 76 novel methylated genes with no known function were identified. Among them, 25 DMRs were further validated that the genome-wide DNA methylation data were reliable through bisulfite treatment, and 9 DMRs were confirmed the relationship between DNA methylation and gene expression by qRT-PCR. Finally, 8 key DMGs possibly associated with soybean CMS were identified. CONCLUSIONS: Genome-wide DNA methylation profile of the soybean CMS line NJCMS5A and its maintainer NJCMS5B was obtained for the first time. Several specific DMGs which participated in pollen and flower development were further identified to be probably associated with soybean CMS. This study will contribute to further understanding of the molecular mechanism behind soybean CMS.

Differential proteomics analysis to identify proteins and pathways associated with male sterility of soybean using iTRAQ-based strategy.

To further elucidate the molecular mechanism of cytoplasmic male sterility (CMS) in soybean, a differential proteomic analysis was completed between the CMS line NJCMS1A and its maintainer NJCMS1B using iTRAQ-based strategy. As a result, 180 differential abundance proteins (DAPs) were identified, of which, 60 were down-regulated and 120 were up-regulated in NJCMS1A compared with NJCMS1B. Bioinformatic analysis showed that 167 DAPs were annotated in 41 Gene Ontology functional groups, 106 DAPs were classified into 20 clusters of orthologous groups of protein categories, and 128 DAPs were enrichment in 53 KEGG pathways. Fifteen differential level proteins/genes with the same expression pattern were identified in the further conjoint analysis of DAPs and the previously reported differential expression genes. Moreover, multiple reaction monitoring test, qRT-PCR analysis and enzyme activity assay validated that the iTRAQ results were reliable. Based on functional analysis of DAPs, we concluded that male sterility in NJCMS1A might be related to insufficiencies in energy supply, unbalance of protein synthesis and degradation, disruption of flavonoid synthesis, programmed cell death, abnormalities of substance metabolism, etc. These results might facilitate our understanding of the molecular mechanisms behind CMS in soybean. BIOLOGICAL SIGNIFICANCE: Soybean is an important global crop that provides protein and oil. Heterosis is a significantly potential approach to increase the yield of soybean. Cytoplasmic male sterility (CMS) plays a vital role in the production of hybrid seeds. However, the genetic and molecular mechanisms of male sterility in soybean still need to be further elucidated. In the present paper, a differential proteomic analysis was carried out and the results showed that several key proteins involved in key pathways were associated with male sterility in soybean. This work provides a new insight to understand the genetic and molecular mechanisms underlying CMS in soybean.

Identification of miRNAs and their targets by high-throughput sequencing and degradome analysis in cytoplasmic male-sterile line NJCMS1A and its maintainer NJCMS1B of soybean.

BACKGROUND: Cytoplasmic male sterility (CMS) provides crucial breeding materials that facilitate hybrid seed production in various crops, and thus plays an important role in the study of hybrid vigor (heterosis), in plants. However, the CMS regulatory network in soybean remains unclear. MicroRNAs (miRNAs) play crucial roles in flower and pollen development by targeting genes that regulate their expression in plants. To identify the miRNAs and their targets that exist in the soybean CMS line NJCMS1A and its maintainer NJCMS1B, high-throughput sequencing and degradome analysis were conducted in this study. RESULTS: Two small RNA libraries were constructed from the flower buds of the soybean CMS line NJCMS1A and its maintainer NJCMS1B. A total of 105 new miRNAs present on the other arm of known pre-miRNAs, 23 new miRNA members, 158 novel miRNAs and 160 high-confidence soybean miRNAs were identified using high-throughput sequencing. Among the identified miRNAs, 101 differentially expressed miRNAs with greater than two-fold changes between NJCMS1A and NJCMS1B were discovered. The different expression levels of selected miRNAs were confirmed by stem-loop quantitative real-time PCR. A degradome analysis showed that 856 targets were predicted to be targeted by 296 miRNAs, including a squamosa promoter-binding protein-like transcription factor family protein, a pentatricopeptide repeat-containing protein, and an auxin response factor, which were previously shown to be involved in floral organ or anther development in plants. Additionally, some targets, including a MADS-box transcription factor, NADP-dependent isocitrate dehydrogenase and NADH-ubiquinone oxidoreductase 24 kDa subunit, were identified, and they may have some relationship with the programmed cell death, reactive oxygen species accumulation and energy deficiencies, which might lead to soybean male sterility. CONCLUSIONS: The present study is the first to use deep sequencing technology to identify miRNAs and their targets in the flower buds of the soybean CMS line NJCMS1A and its maintainer NJCMS1B. The results revealed that the miRNAs might participate in flower and pollen development, which could facilitate our understanding of the molecular mechanisms behind CMS in soybean.

Comparative Transcriptome Analysis between the Cytoplasmic Male Sterile Line NJCMS1A and Its Maintainer NJCMS1B in Soybean (Glycine max (L.) Merr.).

BACKGROUND: The utilization of soybean heterosis is probably one of the potential approaches in future yield breakthrough as was the situation in rice breeding in China. Cytoplasmic male sterility (CMS) plays an important role in the production of hybrid seeds. However, the molecular mechanism of CMS in soybean remains unclear. RESULTS: The comparative transcriptome analysis between cytoplasmic male sterile line NJCMS1A and its near-isogenic maintainer NJCMS1B in soybean was conducted using Illumina sequencing technology. A total of 88,643 transcripts were produced in Illumina sequencing. Then 56,044 genes were obtained matching soybean reference genome. Three hundred and sixty five differentially expressed genes (DEGs) between NJCMS1A and NJCMS1B were screened by threshold, among which, 339 down-regulated and 26 up-regulated in NJCMS1A compared to in NJCMS1B. Gene Ontology (GO) annotation showed that 242 DEGs were annotated to 19 functional categories. Clusters of Orthologous Groups of proteins (COG) annotation showed that 265 DEGs were classified into 19 categories. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that 46 DEGs were assigned to 33 metabolic pathways. According to functional and metabolic pathway analysis combined with reported literatures, the relations between some key DEGs and the male sterility of NJCMS1A were discussed. qRT-PCR analysis validated that the gene expression pattern in RNA-Seq was reliable. Finally, enzyme activity assay showed that energy supply was decreased in NJCMS1A compared to in NJCMS1B. CONCLUSIONS: We concluded that the male sterility of NJCMS1A might be related to the disturbed functions and metabolism pathways of some key DEGs, such as DEGs involved in carbohydrate and energy metabolism, transcription factors, regulation of pollen development, elimination of reactive oxygen species (ROS), cellular signal transduction, and programmed cell death (PCD) etc. Future research will focus on cloning and transgenic function validation of possible candidate genes associated with soybean CMS.