Transcriptomic and photosynthetic responses to grafting of the Nod1 gene in nodulated and non-nodulated soybeans

Chromosome-scale assembly and analysis of Melilotus officinalis genome for SSR development and nodulation genes analysis

Melilotus officinalis is an important legume crop with forage and Chinese medicinal value. The unknown genome of M. officinalis restricted the domestication and utilization of the species and its germplasm resource diversity. A chromosome-scale assembly of the M. officinalis genome was assembled and analysed. The 976.27 Mb of genome was divided into eight chromosomes covering 99.16% of the whole genome. A total of 50022 genes were predicted in the genome. M. officinalis and Melilotus albus shared a common ancestor 0.5-5.65 million years ago (MYA). A genome-wide doubling event occurred 68.93 MYA according to the synonymous nucleotide-substitution values. A total of 552102 tandem repeats were predicted, and 46004 SSR primers of TRs with 10 or more base pairs were developed and designed. The elucidation of the M. officinalis genome provides a compelling model system for studying the genetic, evolutionary and biosynthesis of this legume.

Combined Analysis of BSA-Seq Based Mapping, RNA-Seq, and Metabolomic Unraveled Candidate Genes Associated with Panicle Grain Number in Rice (Oryza sativa L.)

Rice grain yield is a complex and highly variable quantitative trait consisting of several key components, including the grain weight, the effective panicles per unit area, and the grain number per panicle (GNPP). The GNPP is a significant contributor to grain yield controlled by multiple genes (QTL) and is crucial for improvement. Attempts have been made to find genes for this trait, which has always been a challenging and arduous task through conventional methods. We combined a BSA analysis, RNA profiling, and a metabolome analysis in the present study to identify new candidate genes involved in the GNPP. The F₂ population from crossing R4233 (high GNPP) and Ce679 (low GNPP) revealed a frequency distribution fitting two segregated genes. Three pools, including low, middle, and high GNPP, were constructed and a BSA analysis revealed six candidate regions spanning 5.38 Mb, containing 739 annotated genes. Further, a conjunctive analysis of BSA-Seq and RNA-Seq showed 31 differentially expressed genes (DEGs) in the candidate intervals. Subsequently, a metabolome analysis showed 1024 metabolites, with 71 significantly enriched, including 44 up and 27 downregulated in Ce679 vs. R4233. A KEGG enrichment analysis of these 31 DEGs and 71 differentially enriched metabolites (DEMs) showed two genes, Os12g0102100 and Os01g0580500, significantly enriched in the metabolic pathways’ biosynthesis of secondary metabolites, cysteine and methionine metabolism, and fatty acid biosynthesis. Os12g0102100, which encodes for the alcohol dehydrogenase superfamily and a zinc-containing protein, is a novel gene whose contribution to the GNPP is not yet elucidated. This gene coding for mitochondrial trans-2-enoyl-CoA reductase is involved in the biosynthesis of myristic acid, also known as tetradecanoic acid. The Os01g0580500 coding for the enzyme 1-aminoclopropane-1-carboxylate oxidase (OsACO7) is responsible for the final step of the ethylene biosynthesis pathway through the conversion of 1-aminocyclopropane-1-carboxylic acid (ACC) into ethylene. Unlike Os12g0102100, this gene was significantly upregulated in R4233, downregulated in Ce679, and significantly enriched in two of the three metabolite pathways. This result pointed out that these two genes are responsible for the difference in the GNPP in the two cultivars, which has never been identified. Further validation studies may disclose the physiological mechanisms through which they regulate the GNPP in rice.