Deep Sequencing of Small RNA Reveals the Molecular Regulatory Network of AtENO2 Regulating Seed Germination

Seed germination is a key step in the new life cycle of plants. In agriculture, we regard the rapid and consistent process of seed germination as one of the necessary conditions to measure the high quality and yield of crops. ENO2 is a key enzyme in glycolysis, which also plays an important role in plant growth and abiotic stress responses. In our study, we found that the time of seed germination in AtENO2 mutation (eno2⁻) was earlier than that of wild type (WT) in Arabidopsis thaliana. Previous studies have shown that microRNAs (miRNAs) were vital in seed germination. After deep sequencing of small RNA, we found 590 differentially expressed miRNAs in total, of which 87 were significantly differentially expressed miRNAs. By predicting the target genes of miRNAs and analyzing the GO annotation, we have counted 18 genes related to seed germination, including ARF family, TIR1, INVC, RR19, TUDOR2, GA3OX2, PXMT1, and TGA1. MiR9736-z, miR5059-z, ath-miR167a-5p, ath-miR167b, ath-miR5665, ath-miR866-3p, miR10186-z, miR8165-z, ath-miR857, ath-miR399b, ath-miR399c-3p, miR399-y, miR163-z, ath-miR393a-5p, and ath-miR393b-5p are the key miRNAs regulating seed germination-related genes. Through KEGG enrichment analysis, we found that phytohormone signal transduction pathways were significantly enriched, and these miRNAs mentioned above also participate in the regulation of the genes in plant hormone signal transduction pathways, thus affecting the synthesis of plant hormones and further affecting the process of seed germination. This study laid the foundation for further exploration of the AtENO2 regulation for seed germination.     

Identification and Comparative Analysis of CBS Domain-Containing Proteins in Soybean (Glycine max) and the Primary Function of GmCBS21 in Enhanced Tolerance to Low Nitrogen Stress

Nitrogen is an important macronutrient required for plant growth, and is a limiting factor for crop productivity. Improving the nitrogen use efficiency (NUE) is therefore crucial. At present, the NUE mechanism is unclear and information on the genes associated with NUE in soybeans is lacking. cystathionine beta synthase (CBS) domain-containing proteins (CDCPs) may be implicated in abiotic stress tolerance in plants. We identified and classified a CBS domain–containing protein superfamily in soybean. A candidate gene for NUE, GmCBS21, was identified. GmCBS21 gene characteristics, the temporal expression pattern of the GmCBS21 gene, and the phenotype of GmCBS21 overexpression in transgenic Arabidopsis thaliana under low nitrogen stress were analyzed. The phenotypes suggested that the transgenic Arabidopsis thaliana seedlings performed better under the nitrogen-deficient condition. GmCBS21-overexpressing transgenic plants exhibit higher low nitrogen stress tolerance than WT plants, and this suggests its role in low nitrogen stress tolerance in plants. We conclude that GmCBS21 may serve as an excellent candidate for breeding crops with enhanced NUE and better yield.     

Combined BSA-Seq Based Mapping and RNA-Seq Profiling Reveal Candidate Genes Associated with Plant Architecture in Brassica napus

Plant architecture involves important agronomic traits affecting crop yield, resistance to lodging, and fitness for mechanical harvesting in Brassica napus. Breeding high-yield varieties with plant architecture suitable for mechanical harvesting is the main goal of rapeseed breeders. Here, we report an accession of B. napus (4942C-5), which has a dwarf and compact plant architecture in contrast to cultivated varieties. A BC₈ population was constructed by crossing a normal plant architecture line, 8008, with the recurrent parent 4942C-5. To investigate the molecular mechanisms underlying plant architecture, we performed phytohormone profiling, bulk segregant analysis sequencing (BSA-Seq), and RNA sequencing (RNA-Seq) in BC₈ plants with contrasting plant architecture. Genetic analysis indicated the plant architecture traits of 4942C-5 were recessive traits controlled by multiple genes. The content of auxin (IAA), gibberellin (GA), and abscisic acid (ABA) differed significantly between plants with contrasting plant architecture in the BC₈ population. Based on BSA-Seq analysis, we identified five candidate intervals on chromosome A01, namely those of 0 to 6.33 Mb, 6.45 to 6.48 Mb, 6.51 to 6.53 Mb, 6.77 to 6.79 Mb, and 7 to 7.01 Mb regions. The RNA-Seq analysis revealed a total of 4378 differentially expressed genes (DEGs), of which 2801 were up-regulated and 1577 were down-regulated. There, further analysis showed that genes involved in plant hormone biosynthesis and signal transduction, cell structure, and the phenylpropanoid pathway might play a pivotal role in the morphogenesis of plant architecture. Association analysis of BSA-Seq and RNA-Seq suggested that seven DEGs involved in plant hormone signal transduction and a WUSCHEL-related homeobox (WOX) gene (BnaA01g01910D) might be candidate genes responsible for the dwarf and compact phenotype in 4942C-5. These findings provide a foundation for elucidating the mechanisms underlying rapeseed plant architecture and should contribute to breed new varieties suitable for mechanization.     

The Evolution and Functional Roles of miR408 and Its Targets in Plants

MicroRNA408 (miR408) is an ancient and highly conserved miRNA, which is involved in the regulation of plant growth, development and stress response. However, previous research results on the evolution and functional roles of miR408 and its targets are relatively scattered, and there is a lack of a systematic comparison and comprehensive summary of the detailed evolutionary pathways and regulatory mechanisms of miR408 and its targets in plants. Here, we analyzed the evolutionary pathway of miR408 in plants, and summarized the functions of miR408 and its targets in regulating plant growth and development and plant responses to various abiotic and biotic stresses. The evolutionary analysis shows that miR408 is an ancient and highly conserved microRNA, which is widely distributed in different plants. miR408 regulates the growth and development of different plants by down-regulating its targets, encoding blue copper (Cu) proteins, and by transporting Cu to plastocyanin (PC), which affects photosynthesis and ultimately promotes grain yield. In addition, miR408 improves tolerance to stress by down-regulating target genes and enhancing cellular antioxidants, thereby increasing the antioxidant capacity of plants. This review expands and promotes an in-depth understanding of the evolutionary and regulatory roles of miR408 and its targets in plants.     

Exploring the Effect of Methyl Jasmonate on the Expression of microRNAs Involved in Biosynthesis of Active Compounds of Rosemary Cell Suspension Cultures through RNA-Sequencing

Our aim in the experiment was to study the effects of methyl jasmonates (MeJA) on the active compounds of rosemary suspension cells, the metabolites’ change of contents under different concentrations of MeJA, including 0 (CK), 10 (M10), 50 (M50) and 100 μM MeJA (M100). The results demonstrated that MeJA treatments promoted the accumulation of rosmarinic acid (RA), carnosic acid (CA), flavonoids, jasmonate (JA), gibberellin (GA), and auxin (IAA); but reduced the accumulations of abscisic acid (ABA), salicylic acid (SA), and aspartate (Asp). In addition, 50 and 100 μM MeJA promoted the accumulation of alanine (Ala) and glutamate (Glu), and 50 μM MeJA promoted the accumulation of linoleic acid and alpha-linolenic acid in rosemary suspension cells. Comparative RNA-sequencing analysis of different concentrations of MeJA showed that a total of 30, 61, and 39 miRNAs were differentially expressed in the comparisons of CKvsM10, CKvsM50, CKvsM100, respectively. The analysis of the target genes of the differentially expressed miRNAs showed that plant hormone signal transduction, linoleic acid, and alpha-linolenic acid metabolism-related genes were significantly enriched. In addition, we found that miR160a-5p target ARF, miR171d_1 and miR171f_3 target DELLA, miR171b-3p target ETR, and miR156a target BRI1, which played a key role in rosemary suspension cells under MeJA treatments. qRT-PCR of 12 differentially expressed miRNAs and their target genes showed a high correlation between the RNA-seq and the qRT-PCR result. Amplification culture of rosemary suspension cells in a 5 L stirred bioreactor showed that cell biomass accumulation in the bioreactor was less than that in the shake flask under the same conditions, and the whole cultivation period was extended to 14 d. Taken together, MeJA promoted the synthesis of the active compounds in rosemary suspension cells in a wide concentration range via concentration-dependent differential expression patterns. This study provided an overall view of the miRNAs responding to MeJA in rosemary.     

Distinct Evolutionary Profiles and Functions of microRNA156 and microRNA529 in Land Plants

MicroRNA156 (miR156) and miR529 have high sequence similarity and recognize overlapping sites in the same target genes, SQUAMOSA promoter binding protein-like (SPL or SBP box) genes, making it difficult to accurately distinguish their roles in regulatory networks that affect numerous biological functions. Here, we collected data about miR156 and miR529 family members from representative land plants and performed sequence comparisons, phylogenetic analysis, small RNA sequencing, and parallel analysis of RNA ends (PARE) analysis to dissect their evolutionary and functional differences. Although miR156 and miR529 are highly similar, there are differences in their mismatch-sensitive regions, which are essential for target recognition. In land plants, miR156 precursors are conserved mainly within the hairpin region, whereas miR529 precursors are conserved outside the hairpin region, including both the 5’ and 3’ arms. Phylogenetic analysis showed that MIR156 and MIR529 evolved independently, through divergent evolutionary patterns. The two genes also exhibit different expression patterns, with MIR529 preferentially expressed in reproductive tissues and MIR156 in other tissues. PARE analysis revealed that miR156 and miR529 possess specific targets in addition to common targets in maize, pointing to functional differences between them. Based on our findings, we developed a method for the rapid identification of miR529 and miR156 family members and uncovered the evolutionary divergence of these families, providing insights into their different regulatory roles in plant growth and development.     

Soybean miR159-GmMYB33 Regulatory Network Involved in Gibberellin-Modulated Resistance to Heterodera glycines

Soybean cyst nematode (SCN, Heterodera glycines) is an obligate sedentary biotroph that poses major threats to soybean production globally. Recently, multiple miRNAome studies revealed that miRNAs participate in complicated soybean-SCN interactions by regulating their target genes. However, the functional roles of miRNA and target genes regulatory network are still poorly understood. In present study, we firstly investigated the expression patterns of miR159 and targeted GmMYB33 genes. The results showed miR159-3p downregulation during SCN infection; conversely, GmMYB33 genes upregulated. Furthermore, miR159 overexpressing and silencing soybean hairy roots exhibited strong resistance and susceptibility to H. glycines, respectively. In particular, miR159-GAMYB genes are reported to be involve in GA signaling and metabolism. Therefore, we then investigated the effects of GA application on the expression of miR159-GAMYB module and the development of H. glycines. We found that GA directly controls the miR159-GAMYB module, and exogenous GA application enhanced endogenous biologically active GA₁ and GA₃, the abundance of miR159, lowered the expression of GmMYB33 genes and delayed the development of H. glycines. Moreover, SCN infection also results in endogenous GA content decreased in soybean roots. In summary, the soybean miR159-GmMYB33 module was directly involved in the GA-modulated soybean resistance to H. glycines.     

Epigenetic Distribution of Recombinant Plant Chromosome Fragments in a Human–Arabidopsis Hybrid Cell Line

Methylation systems have been conserved during the divergence of plants and animals, although they are regulated by different pathways and enzymes. However, studies on the interactions of the epigenomes among evolutionarily distant organisms are lacking. To address this, we studied the epigenetic modification and gene expression of plant chromosome fragments (~30 Mb) in a human–Arabidopsis hybrid cell line. The whole-genome bisulfite sequencing results demonstrated that recombinant Arabidopsis DNA could retain its plant CG methylation levels even without functional plant methyltransferases, indicating that plant DNA methylation states can be maintained even in a different genomic background. The differential methylation analysis showed that the Arabidopsis DNA was undermethylated in the centromeric region and repetitive elements. Several Arabidopsis genes were still expressed, whereas the expression patterns were not related to the gene function. We concluded that the plant DNA did not maintain the original plant epigenomic landscapes and was under the control of the human genome. This study showed how two diverging genomes can coexist and provided insights into epigenetic modifications and their impact on the regulation of gene expressions between plant and animal genomes.     

Identification of Novel miRNAs and Their Target Genes in the Response to Abscisic Acid in Arabidopsis

miRNAs are involved in various biological processes, including adaptive responses to abiotic stress. To understand the role of miRNAs in the response to ABA, ABA-responsive miRNAs were identified by small RNA sequencing in wild-type Arabidopsis, as well as in abi1td, mkkk17, and mkkk18 mutants. We identified 10 novel miRNAs in WT after ABA treatment, while in abi1td, mkkk17, and mkkk18 mutants, three, seven, and nine known miRNAs, respectively, were differentially expressed after ABA treatment. One novel miRNA (miRn-8) was differentially expressed in the mkkk17 mutant. Potential target genes of the miRNA panel were identified using psRNATarget. Sequencing results were validated by quantitative RT-PCR of several known and novel miRNAs in all genotypes. Of the predicted targets of novel miRNAs, seven target genes of six novel miRNAs were further validated by 5′ RLM-RACE. Gene ontology analyses showed the potential target genes of ABA-responsive known and novel miRNAs to be involved in diverse cellular processes in plants, including development and stomatal movement. These outcomes suggest that a number of the identified miRNAs have crucial roles in plant responses to environmental stress, as well as in plant development, and might have common regulatory roles in the core ABA signaling pathway.