Identification of potential miRNAs and their targets in Vriesea carinata (Poales, Bromeliaceae)

Beyond transcription factors: more regulatory layers affecting soybean gene expression under abiotic stress

Abiotic stresses such as nutritional imbalance, salt, light intensity, and high and low temperatures negatively affect plant growth and development. Through the course of evolution, plants developed multiple mechanisms to cope with environmental variations, such as physiological, morphological, and molecular adaptations. Epigenetic regulation, transcription factor activity, and post-transcriptional regulation operated by RNA molecules are mechanisms associated with gene expression regulation under stress. Epigenetic regulation, including histone and DNA covalent modifications, triggers chromatin remodeling and changes the accessibility of transcription machinery leading to alterations in gene activity and plant homeostasis responses. Soybean is a legume widely produced and whose productivity is deeply affected by abiotic stresses. Many studies explored how soybean faces stress to identify key elements and improve productivity through breeding and genetic engineering. This review summarizes recent progress in soybean gene expression regulation through epigenetic modifications and circRNAs pathways, and points out the knowledge gaps that are important to study by the scientific community. It focuses on epigenetic factors participating in soybean abiotic stress responses, and chromatin modifications in response to stressful environments and draws attention to the regulatory potential of circular RNA in post-transcriptional processing.

Novel approaches on identification of conserved miRNAs for broad-spectrum Potyvirus control measures

Potyviridae comprises more than 200 ssRNA viruses, many of which have a broad host range and geographical distributions. Potyvirids (members of Potyviridae) infect several economically important plants such as saffron, cardamom, cucumber, pepper, potato, tomato, yam, etc. Cumulatively, potyvirids cause a substantial economic loss. The major bottleneck in developing an efficient antiviral strategy is that viruses quickly evade host immunity owing to their higher mutation and recombination rates. Due to this reason, the emergence of newer and improved broad-spectrum approaches to combat viral infections is essential. The use of microRNA's (miRNA) to circumvent viral infection against animal viruses has been successfully employed. Fewer studies reported the development of efficient miRNA-based antivirus resistant strategies against plant viruses and none focused on multiple virus resistance. We focused on potyviruses since studies are limited and identification of conserved miRNAs among various host plants would be an initiative to design broad-spectrum antivirus strategies. In this study, we predicted evolutionarily conserved miRNAs by BLAST searching of reported miRNAs from 15 plants against the GSS and EST sequences of banana. A total of nine miRNAs were predicted and screened in nine diverse potyvirids' hosts (Banana, Tomato, Green gram, Jasmine, Chilli, Coriander, Onion, Rose and Colocasia) belonging to eight different orders (Zingiberales, Solanales, Fabales, Lamiales, Apiales, Asperagales, Rosales and Alismatales). Results suggested that miR168 and miR162 are conserved among all the selected plants. This comprehensive study laid the foundations to design broad-spectrum antivirus resistance using miRNAs. To conclude miR168 and miR162 are conserved among many plants and play a crucial role in evading virus infection which could be used as a potential candidate for developing antiviral strategies against potyvirid infections.

Small RNA sequencing revealed various microRNAs involved in ethylene-triggered flowering process in Aechmea fasciata

Ethylene-triggered flowering is a common phenomenon in plants of the family Bromeliaceae, but its molecular mechanism remains unclear. As a classical group of small RNAs, microRNAs play an essential role in the regulation of flowering. In this study, we found that various miRNAs participate in the ethylene-triggered flowering process in Aechmea fasciata via small RNA sequencing using juvenile and adult plants treated with ethylene for 24 hours. Finally, 63 known miRNAs, 52 novel miRNAs and 1721 target genes were identified or predicted. Expression changes of specific miRNAs were validated by qRT-PCR and northern blotting. Some predicted targets, including SPL, GAMYB and ARF, were verified in RLM-RACE experiments. Gene Ontology (GO) and KEGG analysis showed that numerous developmental and RNA-related processes were enriched. Integrated analysis of the transcriptomic data with small RNA sequencing revealed that numerous miRNAs and targets involved in ethylene-triggered flowering in A. fasciata. Our study is helpful for illuminating the molecular basis of the ethylene-triggered flowering phenomenon in Bromeliaceae.

De novo assembly of Vriesea carinata leaf transcriptome to identify candidate cysteine-proteases

Vriesea carinata is an endemic bromeliad from the Brazilian Atlantic Forest. It has trichome and tank system in their leaves which allows to absorb water and nutrients. It belongs to Bromeliaceae family, which includes several species highly enriched of cysteine-proteases (CysPs). These proteolytic enzymes regulate processes as senescence, cell differentiation, pathogen-linked programmed cell death and mobilization of proteins. Although, their biological importance, there are not genomic resources in V. carinata that can help to identify and understand their molecular mechanisms involved in different biological processes. Thus high-throughput transcriptome sequencing of V. carinata is necessary to generate sequences for the purpose of gene discovery and functional genomic studies. In the present study, we sequenced and assembled the V. carinata transcriptome to the identification of CysPs. A total of 43,232 contigs were assembled for the leaf tissue. BLAST analysis indicated that 23,803 contigs exhibited similarity to non-redundant Viridiplantae proteins. 28.24% of the contigs were classified into the COG database, and gene ontology categorized them into 61 functional groups. A metabolic pathway analysis with KEGG revealed 9679 contigs assigned to 31 metabolic pathways. Among 16 full-length CysPs identified, 11 were evaluated in respect to their expression patterns in the leaf apex, base and inflorescence tissues. The results showed differential expression levels of legumain, metacaspase, pyroglutamyl and papain-like CysPs depending of the leaf region. These results provide a global overview of V. carinata gene functions and expression activities of CysPs in those tissues.

Bioinformatic Exploration of the Targets of Xylem Sap miRNAs in Maize under Cadmium Stress

Cadmium (Cd) has the potential to be chronically toxic to humans through contaminated crop products. MicroRNAs (miRNAs) can move systemically in plants. To investigate the roles of long-distance moving xylem miRNAs in regulating maize response to Cd stress, three xylem sap small RNA (sRNA) libraries were constructed for high-throughput sequencing to identify potential mobile miRNAs in Cd-stressed maize seedlings and their putative targets in maize transcriptomes. In total, about 199 miRNAs (20–22 nucleotides) were identified in xylem sap from maize seedlings, including 97 newly discovered miRNAs and 102 known miRNAs. Among them, 10 miRNAs showed differential expression in xylem sap after 1 h of Cd treatment. Two miRNAs target prediction tools, psRNAtarget (reporting the inhibition pattern of cleavage) and DPMIND (discovering Plant MiRNA-Target Interaction with degradome evidence), were used in combination to identify, via bioinformatics, the targets of 199 significantly expressed miRNAs in maize xylem sap. The integrative results of these two bioinformatic tools suggested that 27 xylem sap miRNAs inhibit 34 genes through cleavage with degradome evidence. Moreover, nearly 300 other genes were also the potential miRNAs cleavable targets without available degradome data support, and the majority of them were enriched in abiotic stress response, cell signaling, transcription regulation, as well as metal handling. These approaches and results not only enhanced our understanding of the Cd-responsive long-distance transported miRNAs from the view of xylem sap, but also provided novel insights for predicting the molecular genetic mechanisms mediated by miRNAs.

Novel and Conserved miRNAs Among Brazilian Pine and Other Gymnosperms

The knowledge about plant miRNAs has increased exponentially, with thousands of miRNAs been reported in different plant taxa using high throughput sequencing technologies and bioinformatic tools. Nevertheless, several groups of plants remain unexplored, and the gap of knowledge about conifer miRNAs is considerable. There is no sequence or functional information available on miRNAs in Araucariaceae. This group is represented in Brazil by only one species, Araucaria angustifolia, an endangered species known as Brazilian pine. In the present study, Brazilian pine has its transcriptome explored with respect to small RNAs, representing the first description in a member of the Araucariaceae family. The screening for conserved miRNAs in Brazilian pine revealed 115 sequences of 30 miRNA families. A total of 106 precursors sequences were predicted. Forty one comprised conserved miRNAs from 16 families, whereas 65 were annotated as novel miRNAs. The comparison of Brazilian pine precursors with sRNA libraries of other five conifer species indicates that 9 out 65 novel miRNAs are conserved among gymnosperms, while 56 seems to be specific for Brazilian pine or restricted to Araucariaceae family. Analysis comparing novel Brazilian pine miRNAs precursors and Araucaria cunninghamii RNA-seq data identified seven orthologs between both species. Mature miRNA identified by bioinformatics predictions were validated using stem-loop RT-qPCR assays. The expression pattern of conserved and novel miRNAs was analyzed in five different tissues of 3-month-old Araucaria seedlings. The present study provides insights about the nature and composition of miRNAs in an Araucariaceae species, with valuable information on miRNAs diversity and conservation in this taxon.

Genome-wide association analysis of the lipid and fatty acid metabolism regulatory network in the mesocarp of oil palm (Elaeis guineensis Jacq.) based on small noncoding RNA sequencing

Oil palm (Elaeis guineensis Jacq.) is the highest oil-yielding crop in the plant kingdom and accumulates 90% of palm oil in the mesocarp. However, the regulatory mechanisms of lipid and fatty acid (FA) metabolism in oil palm are just beginning to be understood, and more studies are needed, especially in the understanding of small noncoding RNA (ncRNA) and mRNA. Based on the deep sequencing of small noncoding RNAs and the degradome in five developmental mesocarp stages, 452 microRNAs (miRNAs), including 170 conserved known-miRNAs (kn-miRNAs) and 282 novel-miRNA (nov-miRNAs), were identified. After predicting the targets of those miRNAs to 37 FA synthesis-related genes, we found that 22 kn-miRNAs and 14 nov-miRNAs might be involved in FA metabolism pathways. Among them, eg-miR156c, eg-miR397, eg-miR444b and nov-miR129 regulated FA synthesis in plastids and the transport of FA-ACP from plastids to the endoplasmic reticulum by targeting acetyl-CoA carboxylase 1 (ACC1), long-chain acyl-CoA synthetase 9 (LACS9), LACS4 and enoyl-ACP reductase (ENR), respectively. Nov-miR138 and nov-miR59 targeted glycerol-3-phosphate acyltransferase (GPAT), and nov-miR274 targeted phosphatidate phosphatase 1 (PAP1). Both target genes are involved in triacylglycerol synthesis in the endoplasmic reticulum. Eg-miR156e and eg-miR156j played pivotal roles by targeting β-ketoacyl-CoA synthase 12 (KCS12), and nov-miR201 targets very-long-chain enoyl-CoA reductase (ECR). Several miRNAs were also predicted to indirectly regulate FA synthesis and lipid metabolism through the squamosa promoter-binding protein-like gene (SPL), NAC and MYB transcription factors. As a whole, indications of a complex and extensive miRNA-mRNA regulatory network associated with FA metabolism in the mesocarp of the oil palm is presented. The results help to broaden the knowledge of potential mechanisms that might be regulated by miRNAs through modulation of the expression of FA-related target gene metabolism in the oil palm.

Development, characterization, and transferability of SSR markers for Vriesea carinata (Bromeliaceae) based on RNA sequencing

PREMISE OF THE STUDY

Expressed sequence tag-simple sequence repeat (EST-SSR) markers were isolated for Vriesea carinata, an endemic bromeliad from the Brazilian Atlantic Forest. These SSR loci may be used to investigate the genetic diversity and population structure of this species and related bromeliads.

METHODS AND RESULTS

Based on the transcriptome data of V. carinata, 30 primer pairs were designed and selected for initial validation. Of these primer pairs, 16 generated suitable SSR loci in 69 individuals. The number of alleles per locus ranged from one to 13; the levels of observed and expected heterozygosity per locus ranged from 0.000 to 1.000 and from 0.000 to 0.935, respectively. All loci produced heterologous amplification. Transferability of the loci was tested in 15 species belonging to three Bromeliaceae subfamilies.

CONCLUSIONS

The developed EST-SSR markers revealed polymorphism in the four studied populations and could be useful to investigate the genetic diversity of V. carinata and related species. The markers may also be suitable for novel gene annotation and discovery.

Integration of the Pokeweed miRNA and mRNA Transcriptomes Reveals Targeting of Jasmonic Acid-Responsive Genes

The American pokeweed plant, Phytolacca americana, displays broad-spectrum resistance to plant viruses and is a heavy metal hyperaccumulator. However, little is known about the regulation of biotic and abiotic stress responses in this non-model plant. To investigate the control of miRNAs in gene expression, we sequenced the small RNA transcriptome of pokeweed treated with jasmonic acid (JA), a hormone that mediates pathogen defense and stress tolerance. We predicted 145 miRNAs responsive to JA, most of which were unique to pokeweed. These miRNAs were low in abundance and condition-specific, with discrete expression change. Integration of paired mRNA-Seq expression data enabled us to identify correlated, novel JA-responsive targets that mediate hormone biosynthesis, signal transduction, and pathogen defense. The expression of approximately half the pairs was positively correlated, an uncommon finding that we functionally validated by mRNA cleavage. Importantly, we report that a pokeweed-specific miRNA targets the transcript of OPR3, novel evidence that a miRNA regulates a JA biosynthesis enzyme. This first large-scale small RNA study of a Phytolaccaceae family member shows that miRNA-mediated control is a significant component of the JA response, associated with widespread changes in expression of genes required for stress adaptation.

High-throughput sequencing and degradome analysis reveal neutral evolution of Cercis gigantea microRNAs and their targets

High-throughput sequencing and degradome analysis for Cercis gigantea identified 194 known miRNAs and 23 novel miRNAs with 61 targets. The comparison results of highly conserved miRNAs and non-conserved miRNAs implied that C. gigantea miRNAs were subjected to the neutral evolution. MicroRNAs play a key role in post-transcriptionally regulating gene expression during plant growth, development and other various biological processes. Although numerous miRNAs have been identified and documented, to date, there are no reports on Cercis gigantea (C. gigantea) miRNAs. In this study, C. gigantea miRNAs and their target genes were investigated by extracting RNA from young roots, tender stems, young leaves, and flower buds of C. gigantea to establish a small RNA and a degradome library to further sequence. This study identified 194 known miRNAs belonging to 52 miRNA families and 23 novel miRNAs. Among these, 158 miRNAs from 27 miRNA families were highly conserved and existed in a plurality of plants. In addition, 60 different targets for 30 known families and one target for novel miRNA were identified by high-throughput sequencing and degradome analysis in C. gigantea. The comparison results revealed that highly conserved miRNAs have higher expression levels, more family members and more targets than non-conserved miRNAs, indicating that C. gigantea miRNAs were subjected to the neutral evolution. Meanwhile, these conserved miRNAs were also found to be involved in auxin signal transduction, regulation of transcription, and other developmental processes, which will help further understanding regulatory mechanisms of C. gigantea miRNAs.