Integrated Analysis of Small RNA, Transcriptome and Degradome Sequencing Provides New Insights into Floral Development and Abscission in Yellow Lupine (Lupinus luteus L.)

A tissue-specific profile of miRNAs and their targets related to paeoniaflorin and monoterpenoids biosynthesis in Paeonia lactiflora Pall. by transcriptome, small RNAs and degradome sequencing

Paeonia lactiflora Pall. (Paeonia) has aroused many concerns due to its extensive medicinal value, in which monoterpene glucoside paeoniflorin and its derivatives are the active chemical components. However, little is known in the molecular mechanism of monoterpenoids biosynthesis, and the regulation network between small RNAs and mRNAs in monoterpenoids biosynthesis has not been investigated yet. Herein, we attempted to reveal the tissue-specific regulation network of miRNAs and their targets related to paeoniaflorin and monoterpenoids biosynthesis in Paeonia by combining mRNA and miRNA expression data with degradome analysis. In all, 289 miRNAs and 30177 unigenes were identified, of which nine miRNAs from seven miRNA families including miR396, miR393, miR835, miR1144, miR3638, miR5794 and miR9555 were verified as monoterpenoids biosynthesis-related miRNAs by degradome sequencing. Moreover, the co-expression network analysis showed that four monoterpenoids-regulating TFs, namely AP2, MYBC1, SPL12 and TCP2, were putatively regulated by five miRNAs including miR172, miR828, miR858, miR156 and miR319, respectively. The present study will improve our knowledge of the molecular mechanisms of the paeoniaflorin and monoterpenoids biosynthesis mediated by miRNA to a new level, and provide a valuable resource for further study on Paeonia.

Transcriptome and targeted hormone metabolome reveal the molecular mechanisms of flower abscission in camellia

INTRODUCTION

Camellia is among the most ornamentally valuable flowers and plants worldwide. Flower abscission typically causes significant financial losses by the horticultural landscape. Previous research has revealed that phytohormones, transcription factors, and other genes involved in floral development regulate the maintenance and mortality of flowers.

METHODS

In this study, for the first time, the transcriptomes and targeted hormone metabolomics of three developmental stages of the receptacles of two distinct camellia strains (CF: abscission strain, CHF: nonabscission strain) were analyzed to determine their roles in regulating blossom abscission in camellia.

RESULTS

ABA content was shown to be considerably upregulated throughout all phases of CF development, as were the genes implicated in the ABA production pathway and their downstream counterparts. Highly expressed genes in CF were involved in galactose metabolism, phenylpropanoid biosynthesis, amino and nucleotide sugar metabolism, pentose and glucuronate interconversions, and MAPK. Among others, highly expressed genes in CHF are associated with fructose and mannose metabolism, alpha-linolenic acid metabolism, biosynthesis of secondary metabolites, starch and sucrose metabolism, and cutin, suberin, and wax biosynthesis. A vast variety of stress response-related pathways and redox-related activities were also shown to be active in CHF. In contrast, CF dramatically activated pathways associated with lignin production, keratinogenesis, cell wall biogenesis, and ABA response. A comparative transcriptomic study of the CF and CHF pathways revealed that the downstream response pathways of hormones, including CTK, BR, IAA, ethylene, and GA, were very active in CF, indicating a significant amount of signal transduction and transcriptional regulation by CF. In addition, members of the transcription factor family, such as MYB, bHLH, MADS, and WD40, may regulate flower abscission.

DISCUSSION

A comparative transcriptome analysis of two distinct strains of camellia receptacles elucidates the molecular processes and regulatory characteristics of flower abscission and provides direction for the targeted improvement and breeding of camellia.

Integrated analysis of miRNAome transcriptome and degradome reveals miRNA-target modules governing floral florescence development and senescence across early- and late-flowering genotypes in tree peony

As a candidate national flower of China, tree peony has extremely high ornamental, medicinal and oil value. However, the short florescence and rarity of early-flowering and late-flowering varieties restrict further improvement of the economic value of tree peony. Specific miRNAs and their target genes engaged in tree peony floral florescence, development and senescence remain unknown. This report presents the integrated analysis of the miRNAome, transcriptome and degradome of tree peony petals collected from blooming, initial flowering, full blooming and decay stages in early-flowering variety Paeonia ostii 'Fengdan', an early-flowering mutant line of Paeonia ostii 'Fengdan' and late-flowering variety Paeonia suffruticosa 'Lianhe'. Transcriptome analysis revealed a transcript ('psu.G.00014095') which was annotated as a xyloglucan endotransglycosylase/hydrolase precursor XTH-25 and found to be differentially expressed across flower developmental stages in Paeonia ostii 'Fengdan' and Paeonia suffruticosa 'Lianhe'. The miRNA-mRNA modules were presented significant enrichment in various pathways such as plant hormone signal transduction, indole alkaloid biosynthesis, arachidonic acid metabolism, folate biosynthesis, fatty acid elongation, and the MAPK signaling pathway. Multiple miRNA-mRNA-TF modules demonstrated the potential functions of MYB-related, bHLH, Trihelix, NAC, GRAS and HD-ZIP TF families in floral florescence, development, and senescence of tree peony. Comparative spatio-temporal expression investigation of eight floral-favored miRNA-target modules suggested that transcript 'psu.T.00024044' and microRNA mtr-miR166g-5p are involved in the floral florescence, development and senescence associated agronomic traits of tree peony. The results might accelerate the understanding of the potential regulation mechanism in regards to floral florescence, development and abscission, and supply guidance for tree peony breeding of varieties with later and longer florescence characteristics.

Involvement of miRNAs in Metabolic Herbicide Resistance to Bispyribac-Sodium in Echinochloa crus-galli (L.) P. Beauv

Several mechanisms involved in weed herbicide resistance are unknown, particularly those acting at the epigenetic level, such as the capacity of small-non-coding RNAs (sncRNAs) to target messenger RNAs of genes involved in herbicide detoxification. The transcription of these sncRNAs is stimulated by epigenetic factors, thereby affecting gene expression. This study was carried out in order to evaluate, for the first time in Echinochloa crus-galli (L.) P. Beauv. (barnyardgrass), the capacity of miRNAs to regulate the expression of genes associated with bispyribac-sodium detoxification. The expression profiles of eight miRNAs with a high degree of complementarity (≥80%) with mRNAs of genes involved in herbicide detoxification (CYP450, GST and eIF4B) were determined by qRT-PCR before and after herbicide spraying. Five of the miRNAs studied (gra-miR7487c, gma-miR396f, gra-miR8759, osa-miR395f, ath-miR847) showed an increased expression after herbicide application in both susceptible and resistant biotypes. All the miRNAs, except gra-miR8759, were more highly expressed in the herbicide-resistant biotypes. In specimens with increased expression of miRNAs, we observed reduced expression of the target genes. The remaining three miRNAs (ata-miR166c-5p, ath-miR396b-5p and osa-miR5538) showed no over-expression after herbicide treatment, and no difference in expression was recorded between susceptible and resistant biotypes. Our results represent a first overview of the capacity of miRNAs to regulate the expression of genes involved in bispyribac-sodium detoxification in the genus Echinochloa. Further research is required to identify novel miRNAs and target genes to develop more focused and sustainable strategies of weed control.

Genome-wide association study of traits in sacred lotus uncovers MITE-associated variants underlying stamen petaloid and petal number variations

Understanding the genetic variants responsible for floral trait diversity is important for the molecular breeding of ornamental flowers. Widely used in water gardening for thousands of years, the sacred lotus exhibits a wide range of diversity in floral organs. Nevertheless, the genetic variations underlying various morphological characteristics in lotus remain largely unclear. Here, we performed a genome-wide association study of sacred lotus for 12 well-recorded ornamental traits. Given a moderate linkage disequilibrium level of 32.9 kb, we successfully identified 149 candidate genes responsible for seven flower traits and plant size variations, including many pleiotropic genes affecting multiple floral-organ-related traits, such as NnKUP2. Notably, we found a 2.75-kb presence-and-absence genomic fragment significantly associated with stamen petaloid and petal number variations, which was further confirmed by re-examining another independent population dataset with petal number records. Intriguingly, this fragment carries MITE transposons bound by siRNAs and is related to the expression differentiation of a nearby candidate gene between few-petalled and double-petalled lotuses. Overall, these genetic variations and candidate genes responsible for diverse lotus traits revealed by our GWAS highlight the role of transposon variations, particularly MITEs, in shaping floral trait diversity.

Integrated mRNA and Small RNA Sequencing Reveals microRNAs Associated With Xylem Development in Dalbergia odorifera

Dalbergia odorifera is a rare and precious rosewood specie, whose wood is a very high-quality material for valuable furniture and carving crafts. However, limited information is available about the process of wood formation in D. odorifera. To determine genes that might be closely associated with the xylem differentiation process, we analyzed the differentially expressed genes (DEGs) and microRNAs (miRNAs) from specific xylem tissues of D. odorifera by RNA sequencing (RNA-seq) and small RNA sequencing (small RNA-seq). In total, we obtained 134,221,955 clean reads from RNA-seq and 90,940,761 clean reads from small RNA-seq. By comparing the transition zone (Dotz) and sapwood (Dosw) samples, a total of 395 DEGs were identified. Further analysis revealed that DEGs encoded for WRKY transcription factors (eight genes), lignin synthesis (PER47, COMT, CCR2), cell wall composition (UXS2), gibberellin synthesis (KAO2, GA20OX1), jasmonic acid synthesis (OPR2, CYP74A), and synthesis of flavonoids (PAL2) and terpenoids (CYP71A1). Subsequently, a preliminary analysis by small RNA-seq showed that the expressions of 14 miRNAs (such as miR168a-5p, miR167f-5p, miR167h-5p, miR167e, miR390a, miR156g, novel_52, and novel_9) were significantly different between Dotz and Dosw. Further analysis revealed that the target genes of these differentially expressed miRNAs were enriched in the GO terms "amino acid binding," "cellulase activity," and "DNA beta-glucosyltransferase activity". Further, KEGG pathway annotation showed significant enrichment in "fatty acid elongation" and "biosynthesis of unsaturated fatty acids". These processes might be participating in the xylem differentiation of D. odorifera. Next, expression correlation analysis showed that nine differentially expressed miRNAs were significantly negatively associated with 21 target genes, which encoded for proteins such as pyrH, SPL6, SPL12, GCS1, and ARF8. Overall, this is the first study on miRNAs and their potential functions in the xylem development of D. odorifera, which provides a stepping stone for a detailed functional investigation of D. odorifera miRNAs.

Progress of Genomics-Driven Approaches for Sustaining Underutilized Legume Crops in the Post-Genomic Era

Legume crops, belonging to the Fabaceae family, are of immense importance for sustaining global food security. Many legumes are profitable crops for smallholder farmers due to their unique ability to fix atmospheric nitrogen and their intrinsic ability to thrive on marginal land with minimum inputs and low cultivation costs. Recent progress in genomics shows promise for future genetic gains in major grain legumes. Still it remains limited in minor legumes/underutilized legumes, including adzuki bean, cluster bean, horse gram, lathyrus, red clover, urd bean, and winged bean. In the last decade, unprecedented progress in completing genome assemblies of various legume crops and resequencing efforts of large germplasm collections has helped to identify the underlying gene(s) for various traits of breeding importance for enhancing genetic gain and contributing to developing climate-resilient cultivars. This review discusses the progress of genomic resource development, including genome-wide molecular markers, key breakthroughs in genome sequencing, genetic linkage maps, and trait mapping for facilitating yield improvement in underutilized legumes. We focus on 1) the progress in genomic-assisted breeding, 2) the role of whole-genome resequencing, pangenomes for underpinning the novel genomic variants underlying trait gene(s), 3) how adaptive traits of wild underutilized legumes could be harnessed to develop climate-resilient cultivars, 4) the progress and status of functional genomics resources, deciphering the underlying trait candidate genes with putative function in underutilized legumes 5) and prospects of novel breeding technologies, such as speed breeding, genomic selection, and genome editing. We conclude the review by discussing the scope for genomic resources developed in underutilized legumes to enhance their production and play a critical role in achieving the "zero hunger" sustainable development goal by 2030 set by the United Nations.

Integrated analysis of mRNA-seq and miRNA-seq reveals the advantage of polyploid Solidago canadensis in sexual reproduction

BACKGROUND

The invasion of Solidago canadensis probably related to polyploidy, which may promotes its potential of sexual reproductive. S. canadensis as an invasive species which rapidly widespread through yield huge numbers of seed, but the mechanism remains unknown. To better understand the advantages of sexual reproduction in hexaploid S. canadensis, transcriptome and small RNA sequencing of diploid and hexaploid cytotypes in flower bud and fruit development stages were performed in this study.

RESULTS

The transcriptome analysis showed that in the flower bud stage, 29 DEGs were MADS-box related genes with 14 up-regulated and 15 down-regulated in hexaploid S. canadensis; 12 SPL genes were detected differentially expressed with 5 up-regulated and 7 down-regulated. In the fruit development stage, 26 MADS-box related genes with 20 up-regulated and 6 down-regulated in hexaploid S. canadensis; 5 SPL genes were all up-regulated; 28 seed storage protein related genes with 18 were up-regulated and 10 down-regulated. The weighted gene co-expression network analysis (WGCNA) identified 19 modules which consisted of co-expressed DEGs with functions such as sexual reproduction, secondary metabolism and transcription factors. Furthermore, we discovered 326 miRNAs with 67 known miRNAs and 259 novel miRNAs. Some of miRNAs, such as miR156, miR156a and miR156f, which target the sexual reproduction related genes.

CONCLUSION

Our study provides a global view of the advantages of sexual reproduction in hexaploid S. canadensis based on the molecular mechanisms, which may promote hexaploid S. canadensis owing higher yield and fruit quality in the process of sexual reproduction and higher germination rate of seeds, and finally conductive to diffusion, faster propagation process and enhanced invasiveness.

Use Chou's 5-steps rule to identify protein post-translational modification and its linkage to secondary metabolism during the floral development of Lonicera japonica Thunb

Lonicera japonica Thunb. is widely used in traditional medicine systems of East Asian and attracts a large amount of studies on the biosynthesis of its active components. Currently, there is little understanding regarding the regulatory mechanisms behind the accumulation of secondary metabolites during its developmental stages. In this study, published transcriptomic and proteomic data were mined to evaluate potential linkage between protein modification and secondary metabolism during the floral development. Stronger correlations were observed between differentially expressed genes (DEGs) and their corresponding differentially abundant proteins (DAPs) in the comparison of juvenile bud stage (JBS)/third green stage (TGS) vs. silver flowering stage (SFS). Seventy-five and 76 cor-rDEGs and cor-rDAPs (CDDs) showed opposite trends at both transcriptional and translational levels when comparing their levels at JBS and TGS relative to those at SFS. CDDs were mainly involved in elements belonging to the protein metabolism and the TCA cycle. Protein-protein interaction analysis indicated that the interacting proteins in the major cluster were primarily involved in TCA cycle and protein metabolism. In the simple phenylpropanoids biosynthetic pathway of SFS, both phospho-2-dehydro-3-deoxyheptonate aldolase (PDA) and glutamate/aspartate-prephenate aminotransferase (AAT) were decreased at the protein level, but increased at the gene level. A confirmatory experiment indicated that protein ubiquitination and succinylation were more prominent during the early floral developmental stages, in correlation with simple phenylpropanoids accumulation. Taken together, those data indicates that phenylpropanoids metabolism and floral development are putatively regulated through the ubiquitination and succinylation modifications of PDA, AAT, and TCA cycle proteins in L. japonica.

Non-Coding RNAs in Legumes: Their Emerging Roles in Regulating Biotic/Abiotic Stress Responses and Plant Growth and Development

Noncoding RNAs, including microRNAs (miRNAs), small interference RNAs (siRNAs), circular RNA (circRNA), and long noncoding RNAs (lncRNAs), control gene expression at the transcription, post-transcription, and translation levels. Apart from protein-coding genes, accumulating evidence supports ncRNAs playing a critical role in shaping plant growth and development and biotic and abiotic stress responses in various species, including legume crops. Noncoding RNAs (ncRNAs) interact with DNA, RNA, and proteins, modulating their target genes. However, the regulatory mechanisms controlling these cellular processes are not well understood. Here, we discuss the features of various ncRNAs, including their emerging role in contributing to biotic/abiotic stress response and plant growth and development, in addition to the molecular mechanisms involved, focusing on legume crops. Unravelling the underlying molecular mechanisms and functional implications of ncRNAs will enhance our understanding of the coordinated regulation of plant defences against various biotic and abiotic stresses and for key growth and development processes to better design various legume crops for global food security.

Origin and Evolutionary Dynamics of the miR2119 and ADH1 Regulatory Module in Legumes

MicroRNAs are important regulators of gene expression in eukaryotes. Previously, we reported that in Phaseolus vulgaris, the precursor for miR2119 is located in the same gene as miR398a, conceiving a dicistronic MIR gene. Both miRNA precursors are transcribed and processed from a single transcript resulting in two mature microRNAs that regulate the mRNAs encoding ALCOHOL DEHYDROGENASE 1 (ADH1) and COPPER-ZINC SUPEROXIDE DISMUTASE 1 (CSD1). Genes for miR398 are distributed throughout the spermatophytes; however, miR2119 is only found in Leguminosae species, indicating its recent emergence. Here, we used public databases to explore the presence of the miR2119 sequence in several plant species. We found that miR2119 is present only in specific clades within the Papilionoideae subfamily, including important crops used for human consumption and forage. Within this subfamily, MIR2119 and MIR398a are found together as a single gene in the genomes of the Millettioids and Hologalegina. In contrast, in the Dalbergioids MIR2119 is located in a different locus from MIR398a, suggesting this as the ancestral genomic organization. To our knowledge, this is a unique example where two separate MIRNA genes have merged to generate a single polycistronic gene. Phylogenetic analysis of ADH1 gene sequences in the Papilionoideae subfamily revealed duplication events resulting in up to four ADH1 genes in certain species. Notably, the presence of MIR2119 correlates with the conservation of target sites in particular ADH1 genes in each clade. Our results suggest that post-transcriptional regulation of ADH1 genes by miR2119 has contributed to shaping the expansion and divergence of this gene family in the Papilionoideae. Future experimental work on ADH1 regulation by miR2119 in more legume species will help to further understand the evolutionary history of the ADH1 gene family and the relevance of miRNA regulation in this process.

Integrated Analysis of Small RNA, Transcriptome, and Degradome Sequencing Reveals the Water-Deficit and Heat Stress Response Network in Durum Wheat

Water-deficit and heat stress negatively impact crop production. Mechanisms underlying the response of durum wheat to such stresses are not well understood. With the new durum wheat genome assembly, we conducted the first multi-omics analysis with next-generation sequencing, providing a comprehensive description of the durum wheat small RNAome (sRNAome), mRNA transcriptome, and degradome. Single and combined water-deficit and heat stress were applied to stress-tolerant and -sensitive Australian genotypes to study their response at multiple time-points during reproduction. Analysis of 120 sRNA libraries identified 523 microRNAs (miRNAs), of which 55 were novel. Differentially expressed miRNAs (DEMs) were identified that had significantly altered expression subject to stress type, genotype, and time-point. Transcriptome sequencing identified 49,436 genes, with differentially expressed genes (DEGs) linked to processes associated with hormone homeostasis, photosynthesis, and signaling. With the first durum wheat degradome report, over 100,000 transcript target sites were characterized, and new miRNA-mRNA regulatory pairs were discovered. Integrated omics analysis identified key miRNA-mRNA modules (particularly, novel pairs of miRNAs and transcription factors) with antagonistic regulatory patterns subject to different stresses. GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis revealed significant roles in plant growth and stress adaptation. Our research provides novel and fundamental knowledge, at the whole-genome level, for transcriptional and post-transcriptional stress regulation in durum wheat.

LuluDB-The Database Created Based on Small RNA, Transcriptome, and Degradome Sequencing Shows the Wide Landscape of Non-coding and Coding RNA in Yellow Lupine (Lupinus luteus L.) Flowers and Pods

Yellow lupine (Lupinus luteus L.) belongs to a legume family that benefits from symbiosis with nitrogen-fixing bacteria. Its seeds are rich in protein, which makes it a valuable food source for animals and humans. Yellow lupine is also the model plant for basic research on nodulation or abscission of organs. Nevertheless, the knowledge about the molecular regulatory mechanisms of its generative development is still incomplete. The RNA-Seq technique is becoming more prominent in high-throughput identification and expression profiling of both coding and non-coding RNA sequences. However, the huge amount of data generated with this method may discourage other scientific groups from making full use of them. To overcome this inconvenience, we have created a database containing analysis-ready information about non-coding and coding L. luteus RNA sequences (LuluDB). LuluDB was created on the basis of RNA-Seq analysis of small RNA, transcriptome, and degradome libraries obtained from yellow lupine cv. Taper flowers, pod walls, and seeds in various stages of development, flower pedicels, and pods undergoing abscission or maintained on the plant. It contains sequences of miRNAs and phased siRNAs identified in L. luteus, information about their expression in individual samples, and their target sequences. LuluDB also contains identified lncRNAs and protein-coding RNA sequences with their organ expression and annotations to widely used databases like GO, KEGG, NCBI, Rfam, Pfam, etc. The database also provides sequence homology search by BLAST using, e.g., an unknown sequence as a query. To present the full capabilities offered by our database, we performed a case study concerning transcripts annotated as DCL 1–4 (DICER LIKE 1–4) homologs involved in small non-coding RNA biogenesis and identified miRNAs that most likely regulate DCL1 and DCL2 expression in yellow lupine. LuluDB is available at http://luluseqdb.umk.pl/basic/web/index.php.

Legume Genetics and Biology: From Mendel's Pea to Legume Genomics

Legumes have played an important part in cropping systems since the dawn of agriculture, both as human food and as animal feed. The legume family is arguably one of the most abundantly domesticated crop plant families. Their ability to symbiotically fix nitrogen and improve soil fertility has been rewarded since antiquity and makes them a key protein source. The pea was the original model organism used in Mendel’s discovery of the laws of inheritance, making it the foundation of modern plant genetics. This Special Issue provides up-to-date information on legume biology, genetic advances, and the legacy of Mendel.