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.     

Identification and Functional Analysis of MicroRNAs and Their Targets in Platanus acerifolia under Lead (Pb) Stress

MicroRNAs (miRNAs) play important regulatory roles in development and stress responses in plants. Lead (Pb) is a non-essential element that is highly toxic to living organisms. Platanus acerifolia is grown as a street tree in cities throughout temperate regions for its importance in improving the urban ecological environment. MiRNAs that respond to abiotic stresses have been identified in plants; however, until now, the influence of Pb stress on P. acerifolia miRNAs has not been reported. To identify miRNAs and predict their target genes under Pb stress, two small RNA and two degradome libraries were constructed from Pb-treated and Pb-free leaves of P. acerifolia seedlings. After sequencing, 55 known miRNAs and 129 novel miRNAs were obtained, and 104 target genes for the miRNAs were identified by degradome sequencing. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses were performed to predict the functions of the targets. The expressions of eight differentially expressed miRNAs were validated by quantitative real-time polymerase chain reaction (qRT-PCR). This is the first report about P. acerifolia miRNAs and their target genes under Pb stress. This study has provided data for further research into molecular mechanisms involved in resistance of P. acerifolia to Pb stress.     

Exogenous Brassinosteroid Facilitates Xylem Development in Pinus massoniana Seedlings

Brassinosteroids (BRs) are known to be essential regulators for wood formation in herbaceous plants and poplar, but their roles in secondary growth and xylem development are still not well-defined, especially in pines. Here, we treated Pinus massoniana seedlings with different concentrations of exogenous BRs, and assayed the effects on plant growth, xylem development, endogenous phytohormone contents and gene expression within stems. Application of exogenous BR resulted in improving development of xylem more than phloem, and promoting xylem development in a dosage-dependent manner in a certain concentration rage. Endogenous hormone determination showed that BR may interact with other phytohormones in regulating xylem development. RNA-seq analysis revealed that some conventional phenylpropanoid biosynthesis- or lignin synthesis-related genes were downregulated, but the lignin content was elevated, suggesting that new lignin synthesis pathways or other cell wall components should be activated by BR treatment in P. massoniana. The results presented here reveal the foundational role of BRs in regulating plant secondary growth, and provide the basis for understanding molecular mechanisms of xylem development in P. massoniana.     

Conserved miRNAs and Their Response to Salt Stress in Wild Eggplant Solanum linnaeanum Roots

The Solanaceae family includes some important vegetable crops, and they often suffer from salinity stress. Some miRNAs have been identified to regulate gene expression in plant response to salt stress; however, little is known about the involvement of miRNAs in Solanaceae species. To identify salt-responsive miRNAs, high-throughput sequencing was used to sequence libraries constructed from roots of the salt tolerant species, Solanum linnaeanum, treated with and without NaCl. The sequencing identified 98 conserved miRNAs corresponding to 37 families, and some of these miRNAs and their expression were verified by quantitative real-time PCR. Under the salt stress, 11 of the miRNAs were down-regulated, and 3 of the miRNAs were up-regulated. Potential targets of the salt-responsive miRNAs were predicted to be involved in diverse cellular processes in plants. This investigation provides valuable information for functional characterization of miRNAs in S. linnaeanum, and would be useful for developing strategies for the genetic improvement of the Solanaceae crops.     

Transcriptome Analyses Identify Potential Key microRNAs and Their Target Genes Contributing to Ovarian Reserve

Female endocrinological symptoms, such as premature ovarian inefficiency (POI) are caused by diminished ovarian reserve and chemotherapy. The etiology of POI remains unknown, but this can lead to infertility. This has accelerated the search for master regulator genes or other molecules that contribute as enhancers or silencers. The impact of regulatory microRNAs (miRNAs) on POI has gained attention; however, their regulatory function in this condition is not well known. RNA sequencing was performed at four stages, 2-(2 W), 6-(6 W), 15-(15 W), and 20-(20 W) weeks, on ovarian tissue samples and 5058 differentially expressed genes (DEGs) were identified. Gene expression and enrichment were analyzed based on the gene ontology and KEGG databases, and their association with other proteins was assessed using the STRING database. Gene set enrichment analysis was performed to identify the key target genes. The DEGs were most highly enriched in 6 W and 15 W groups. Figla, GDF9, Nobox, and Pou51 were significantly in-creased at 2 W compared with levels at 6 W and 20 W, whereas the expression of Foxo1, Inha, and Taf4b was significantly de-creased at 20 W. Ccnd2 and Igf1 expression was maintained at similar levels in each stage. In total, 27 genes were upregulated and 26 genes interacted with miRNAs; moreover, stage-specific upregulated and downregulated interactions were demonstrated. Increased and decreased miRNAs were identified at each stage in the ovaries. The constitutively expressed genes, Ccnd2 and Igf1, were identified as the major targets of many miRNAs (p < 0.05), and Fshr and Foxo3 interacted with miRNAs, namely mmu-miR-670-3p and mmu-miR-153-3p. miR-26a-5p interacted with Piwil2, and its target genes were downregulated in the 20 W mouse ovary. In this study, we aimed to identify key miRNAs and their target genes encompassing the reproductive span of mouse ovaries using mRNA and miRNA sequencing. These results indicated that gene sets are regulated in the reproductive stage-specific manner via interaction with miRNAs. Furthermore, consistent expression of Ccnd2 and Igf1 is considered crucial for the ovarian reserve and is regulated by many interactive miRNAs.     

Deep Sequencing Analysis of miRNA Expression in Breast Muscle of Fast-Growing and Slow-Growing Broilers

Growth performance is an important economic trait in chicken. MicroRNAs (miRNAs) have been shown to play important roles in various biological processes, but their functions in chicken growth are not yet clear. To investigate the function of miRNAs in chicken growth, breast muscle tissues of the two-tail samples (highest and lowest body weight) from Recessive White Rock (WRR) and Xinghua Chickens (XH) were performed on high throughput small RNA deep sequencing. In this study, a total of 921 miRNAs were identified, including 733 known mature miRNAs and 188 novel miRNAs. There were 200, 279, 257 and 297 differentially expressed miRNAs in the comparisons of WRRh vs. WRRl, WRRh vs. XHh, WRRl vs. XHl, and XHh vs. XHl group, respectively. A total of 22 highly differentially expressed miRNAs (fold change > 2 or < 0.5; p-value < 0.05; q-value < 0.01), which also have abundant expression (read counts > 1000) were found in our comparisons. As far as two analyses (WRRh vs. WRRl, and XHh vs. XHl) are concerned, we found 80 common differentially expressed miRNAs, while 110 miRNAs were found in WRRh vs. XHh and WRRl vs. XHl. Furthermore, 26 common miRNAs were identified among all four comparisons. Four differentially expressed miRNAs (miR-223, miR-16, miR-205a and miR-222b-5p) were validated by quantitative real-time RT-PCR (qRT-PCR). Regulatory networks of interactions among miRNAs and their targets were constructed using integrative miRNA target-prediction and network-analysis. Growth hormone receptor (GHR) was confirmed as a target of miR-146b-3p by dual-luciferase assay and qPCR, indicating that miR-34c, miR-223, miR-146b-3p, miR-21 and miR-205a are key growth-related target genes in the network. These miRNAs are proposed as candidate miRNAs for future studies concerning miRNA-target function on regulation of chicken growth.     

Identification of microRNAs in the Lyme Disease Vector Ixodes scapularis

MicroRNAs (miRNAs) are a class of small non-coding RNAs involved in many biological processes, including the immune pathways that control bacterial, parasitic, and viral infections. Pathogens probably modify host miRNAs to facilitate successful infection, so they might be useful targets for vaccination strategies. There are few data on differentially expressed miRNAs in the black-legged tick Ixodes scapularis after infection with Borrelia burgdorferi, the causative agent of Lyme disease in the United States. Small RNA sequencing and qRT-PCR analysis were used to identify and validate differentially expressed I. scapularis salivary miRNAs. Small RNA-seq yielded 133,465,828 (≥18 nucleotides) and 163,852,135 (≥18 nucleotides) small RNA reads from Borrelia-infected and uninfected salivary glands for downstream analysis using the miRDeep2 algorithm. As such, 254 miRNAs were identified across all datasets, 25 of which were high confidence and 51 low confidence known miRNAs. Further, 23 miRNAs were differentially expressed in uninfected and infected salivary glands: 11 were upregulated and 12 were downregulated upon pathogen infection. Gene ontology and network analysis of target genes of differentially expressed miRNAs predicted roles in metabolic, cellular, development, cellular component biogenesis, and biological regulation processes. Several Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, including sphingolipid metabolism; valine, leucine and isoleucine degradation; lipid transport and metabolism; exosome biogenesis and secretion; and phosphate-containing compound metabolic processes, were predicted as targets of differentially expressed miRNAs. A qRT-PCR assay was utilized to validate the differential expression of miRNAs. This study provides new insights into the miRNAs expressed in I. scapularis salivary glands and paves the way for their functional manipulation to prevent or treat B. burgdorferi infection.     

FvNST1b NAC Protein Induces Secondary Cell Wall Formation in Strawberry

Secondary cell wall thickening plays a crucial role in plant growth and development. Diploid woodland strawberry (Fragaria vesca) is an excellent model for studying fruit development, but its molecular control of secondary wall thickening is largely unknown. Previous studies have shown that Arabidopsis NAC secondary wall thickening promoting factor1 (AtNST1) and related proteins are master regulators of xylem fiber cell differentiation in multiple plant species. In this study, a NST1-like gene, FvNST1b, was isolated and characterized from strawberry. Sequence alignment and phylogenetic analysis showed that the FvNST1b protein contains a highly conserved NAC domain, and it belongs to the same family as AtNST1. Overexpression of FvNST1b in wild-type Arabidopsis caused extreme dwarfism, induced ectopic thickening of secondary walls in various tissues, and upregulated the expression of genes related to secondary cell wall synthesis. In addition, transient overexpression of FvNST1b in wild-type Fragaria vesca fruit produced cells resembling tracheary elements. These results suggest that FvNST1b positively regulates secondary cell wall formation as orthologous genes from other species.     

Differential MicroRNA Expression in Porcine Endometrium Related to Spontaneous Embryo Loss during Early Pregnancy

Litter size is an important indicator to measure the production capacity of commercial pigs. Spontaneous embryo loss is an essential factor in determining sow litter size. In early pregnancy, spontaneous embryo loss in porcine is as high as 20–30% during embryo implantation. However, the specific molecular mechanism underlying spontaneous embryo loss at the end of embryo implantation remains unknown. Therefore, we comprehensively used small RNA sequencing technology, bioinformatics analysis, and molecular experiments to determine the microRNA (miRNA) expression profile in the healthy and arresting embryo implantation site of porcine endometrium on day of gestation (DG) 28. A total of 464 miRNAs were identified in arresting endometrium (AE) and healthy endometrium (HE), and 139 differentially expressed miRNAs (DEMs) were screened. We combined the mRNA sequencing dataset from the SRA database to predict the target genes of these miRNAs. A quantitative real-time PCR assay identified the expression levels of miRNAs and mRNAs. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were performed on differentially expressed target genes of DEMs, mainly enriched in epithelial development and amino acids metabolism-related pathways. We performed fluorescence in situ hybridization (FISH) and the dual-luciferase report gene assay to confirm miRNA and predicted target gene binding. miR-205 may inhibit its expression by combining 3′-untranslated regions (3′ UTR) of tubulointerstitial nephritis antigen-like 1 (TINAGL1). The resulting inhibition of angiogenesis in the maternal endometrium ultimately leads to the formation of arresting embryos during the implantation period. This study provides a reference for the effect of miRNA on the successful implantation of pig embryos in early gestation.     

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.