Screening and identification of miRNAs related to sexual differentiation of strobili in Ginkgo biloba by integration analysis of small RNA, RNA, and degradome sequencing

MicroRNAs constitute an additional layer in plant response to simultaneous bio- and abiotic stresses as exemplified by UV-B radiation and flg22-treatment on Arabidopsis thaliana

Plants are confronted with various environmental stresses and develop sophisticated adaptive mechanisms. Our previous work demonstrated that the crosstalk of flg22 and ultraviolet (UV)-B-induced signalling cascades reprograms the expression of flavonol pathway genes (FPGs), benefiting plant defence responses. Although several transcription factors have been identified to be involved in this crosstalk, the underlying mechanism is largely unclear. Here, we analyzed microRNAs (miRNAs) and identified 126, 129 and 113 miRNAs with altered abundances compared to untreated control in flg22-, UV-B- and flg22/UV-B-treated seedlings, respectively. Two distinct modules were identified: The first consists of 10 miRNAs repressed by UV-B but up-regulated by flg22, and the second with five miRNAs repressed by flg22 but up-regulated by UV-B. In Arabidopsis, the knockdown of miR858a, a representative of module I, increased the abundance of CHS (a marker gene for FPGs), whereas its overexpression reduced CHS. Conversely, knockout of miR164b from module II decreased CHS and its overexpression increased CHS transcript levels. These data suggest a decisive role of miRNAs in the crosstalk. In the next, we described the interaction between miR858a and its target MYB111 (a positive regulator of FPGs) from module I in detail. We showed that MYB111 was profoundly post-transcriptionally regulated by miR858a during the crosstalk, whose expression was specifically but antagonistically controlled by UVR8- and FLS2-mediated signallings. Moreover, transcriptional monitoring using the GUS reporter gene demonstrates that miRNA-mediated posttranscriptional regulation is the main driving force in reprogramming the expression of FPGs and regulates plant adaptation to multiple concurrent environmental stresses.

Genome-wide identification and expression pattern analysis of the TCP transcription factor family in Ginkgo biloba

Plant-specific TCP transcription factors play an essential role in plant growth and development. They can regulate leaf curvature, flower symmetry and the synthesis of secondary metabolites. The flavonoids in Ginkgo biloba leaf are one of the main medicinally bioactivate compounds, which have pharmacological and beneficial health effects for humans. In this study, a total of 13 TCP genes were identified in G. biloba, and 5 of them belonged to PCF subclades (GbTCP03, GbTCP07, GbTCP05, GbTCP13, GbTCP02) while others belonged to CIN (GbTCP01, GbTCP04, GbTCP06, GbTCP08, GbTCP09, GbTCP10, GbTCP11, GbTCP12) subclades according to phylogenetic analysis. Numerous cis-acting elements related to various biotic and abiotic signals were predicted on the promoters by cis-element analysis, suggesting that the expression of GbTCPs might be co-regulated by multiple signals. Transcript abundance analysis exhibited that most of GbTCPs responded to multiple phytohormones. Among them, the relative expression levels of GbTCP06, GbTCP11, and GbTCP13 were found to be significantly influenced by exogenous ABA, SA and MeJA application. In addition, a total of 126 miRNAs were predicted to target 9 TCPs (including GbTCP01, GbTCP02, GbTCP04, GbTCP05, GbTCP06, GbTCP08, GbTCP11, GbTCP12, GbTCP13). The correlation analysis between the expression level of GbTCPs and the flavonoid contents showed that GbTCP03, GbTCP04, GbTCP07 might involve in flavonoid biosynthesis in G. biloba. In short, this study mainly provided a theoretical foundation for better understanding the potential function of TCPs in G. biloba.

Identification of miRNAs and Target Genes at Key Stages of Sexual Differentiation in Androdioecious Osmanthus fragrans

Androdioecy is the crucial transition state in the evolutionary direction of hermaphroditism to dioecy, however, the molecular mechanisms underlying the formation of this sex system remain unclear. While popular in China for its ornamental and cultural value, Osmanthus fragrans has an extremely rare androdioecy breeding system, meaning that there are both male and hermaphroditic plants in a population. To unravel the mechanisms underlying the formation of androdioecy, we performed small RNA sequencing studies on male and hermaphroditic O. fragrans. A total of 334 miRNAs were identified, of which 59 were differentially expressed. Functional categorization revealed that the target genes of differentially expressed miRNAs were mainly involved in the biological processes of reproductive development and the hormone signal transduction pathway. We speculated that the miRNA160, miRNA167, miRNA393 and miRNA396 families may influence the sex differentiation in O. fragrans. Overall, our study is the first exploration of miRNAs in the growth and development process of O. fragrans, and is also the first study of androdioecious plants from the miRNA sequencing perspective. The analysis of miRNAs and target genes that may be involved in the sex differentiation process lay a foundation for the ultimate discovery of the androdioecious molecular mechanism in O. fragrans.

Transcriptome mining of genes in Zanthoxylum armatum revealed ZaMYB86 as a negative regulator of prickly development

Zanthoxylum armatum DC. is an important economic tree species. Prickle is a type of trichome with special morphology, and there are a lot of prickles on the leaves of Z. armatum, which seriously restricts the development of Z. armatum industry. In this study, the leaves of Z. armatum cv. Zhuye (ZY) and its budding variety 'Rongchangwuci' (WC) (A less prickly mutant variety) at different developmental stages were used as materials, and the transcriptome sequencing data were analyzed. A total of 96,931 differentially expressed genes (DEGs) were identified among the samples, among which 1560 were candidate DEGs that might be involved in hormone metabolism. The contents of JA, auxin and CK phytohormones in ZY leaves were significantly higher than those in WC leaves. Combined with weighted gene co-expression network analysis, eight genes (MYC, IAA, ARF, CRE/AHK, PP2C, ARR-A, AOS and LOX) were identified, including 25 transcripts, which might affect the metabolism of the three hormones and indirectly participate in the formation of prickles. Combining with the proteins successfully reported in other plants to regulate trichome formation, ZaMYB86, a transcription factor of R2R3 MYB family, was identified through local Blast and phylogenetic tree analysis, which might regulate prickle formation of Z. armatum. Overexpression of ZaMYB86 in mutant A. thaliana resulted in the reduction of trichomes in A. thaliana leaves, which further verified that ZaMYB86 was involved in the formation of pickles. Yeast two-hybrid results showed that ZaMYB86 interacted with ZaMYB5. Furthermore, ZaMYB5 was highly homologous to AtMYB5, a transcription factor that regulated trichomes development, in MYB DNA binding domain. Taken together, these results indicated that ZaMYB86 and ZaMYB5 act together to regulate the formation of prickles in Z. armatum. Our findings provided a new perspective for revealing the molecular mechanism of prickly formation.

Significance of miRNA in enhancement of flavonoid biosynthesis

Flavonoid metabolism shows very strong plasticity in plant development and coping with the changing environment. Flavonoid biosynthesis is regulated by many metabolic pathways, including transcriptional regulation, post-transcriptional control, post-translational regulationand epigenetic regulation. miRNA is a form of endogenous noncoding single-strand small molecule RNA that primarily regulates the expression of target genes horizontally after transcription through splicing and translational suppression. It also plays an important role in regulating plant growth and development, secondary metabolism and biotic and abiotic stress. miRNA can regulate the formation of flavonoids by acting on structural genes or indirectly by using an MBW transcription complex comprising MYB-bHLH-WD40. This study summarizes the biosynthesis and mechanisms of miRNA, and provides a summary of the mechanisms of miRNAs involved in production of flavonoids, in order to elucidate the biosynthesis pathway and complex regulatory network of plant flavonoids. We aim to provide new insights into improving the content of flavonoid active ingredients in plants.

The emerging role of small RNAs in ovule development, a kind of magic

In plants, small RNAs have been recognized as key genetic and epigenetic regulators of development. Small RNAs are usually 20 to 30 nucleotides in length and they control, in a sequence specific manner, the transcriptional or post-transcriptional expression of genes. In this review, we present a comprehensive overview of the most recent findings about the function of small RNAs in ovule development, including megasporogenesis and megagametogenesis, both in sexual and apomictic plants. We discuss recent studies on the role of miRNAs, siRNAs and trans-acting RNAs (ta-siRNAs) in early female germline differentiation. The mechanistic complexity and unique regulatory features are reviewed, and possible directions for future research are provided.

MicroRNAs May Play an Important Role in Sexual Reversal Process of Chinese Soft-Shelled Turtle, Pelodiscus sinensis

The Chinese soft-shelled (Pelodiscus sinensis) turtle exhibits obvious sex dimorphism, which leads to the higher economic and nutritional value of male individuals. Exogenous hormones can cause the transformation from male to female phenotype during gonadal differentiation. However, the molecular mechanism related to the sexual reversal process is unclear. In this study, we compared the difference between the small RNAs of male, female, and pseudo-female turtles by small RNA-seq to understand the sexual reversal process of Chinese soft-shelled turtles. A certain dose of estrogen can cause the transformation of Chinese soft-shelled turtles from male to female, which are called pseudo-female individuals. The result of small RNA-seq has revealed that the characteristics of pseudo-females are very similar to females, but are strikingly different from males. The number of the microRNAs (miRNAs) of male individuals was significantly less than the number of female individuals or pseudo-female individuals, while the expression level of miRNAs of male individuals were significantly higher than the other two types. Furthermore, we found 533 differentially expressed miRNAs, including 173 up-regulated miRNAs and 360 down-regulated miRNAs, in the process of transformation from male to female phenotype. Cluster analysis of the total 602 differential miRNAs among females, males, and pseudo-females showed that miRNAs played a crucial role during the sexual differentiation. Among these differential miRNAs, we found 12 miRNAs related to gonadal development and verified their expression by qPCR. The TR-qPCR results confirmed the differential expression of 6 of the 12 miRNAs: miR-26a-5p, miR-212-5p, miR-202-5p, miR-301a, miR-181b-3p and miR-96-5p were involved in sexual reversal process, which was consistent with the results of omics. Using these six miRNAs and some of their target genes, we constructed a network diagram related to gonadal development. We suggest that these miRNAs may play an important role in the process of effective sex reversal, which would contribute to the breeding of all male strains of Chinese soft-shelled turtles.