Ana Paula Korbes

Identifying conserved and novel microRNAs in developing seeds of Brassica napus using deep sequencing.

MicroRNAs (miRNAs) are important post-transcriptional regulators of plant development and seed formation. In Brassica napus, an important edible oil crop, valuable lipids are synthesized and stored in specific seed tissues during embryogenesis. The miRNA transcriptome of B. napus is currently poorly characterized, especially at different seed developmental stages. This work aims to describe the miRNAome of developing seeds of B. napus by identifying plant-conserved and novel miRNAs and comparing miRNA abundance in mature versus developing seeds. Members of 59 miRNA families were detected through a computational analysis of a large number of reads obtained from deep sequencing two small RNA and two RNA-seq libraries of (i) pooled immature developing stages and (ii) mature B. napus seeds. Among these miRNA families, 17 families are currently known to exist in B. napus; additionally 29 families not reported in B. napus but conserved in other plant species were identified by alignment with known plant mature miRNAs. Assembled mRNA-seq contigs allowed for a search of putative new precursors and led to the identification of 13 novel miRNA families. Analysis of miRNA population between libraries reveals that several miRNAs and isomiRNAs have different abundance in developing stages compared to mature seeds. The predicted miRNA target genes encode a broad range of proteins related to seed development and energy storage. This work presents a comparative study of the miRNA transcriptome of mature and developing B. napus seeds and provides a basis for future research on individual miRNAs and their functions in embryogenesis, seed maturation and lipid accumulation in B. napus.

Identifying MicroRNAs and Transcript Targets in Jatropha Seeds

MicroRNAs, or miRNAs, are endogenously encoded small RNAs that play a key role in diverse plant biological processes. Jatropha curcas L. has received significant attention as a potential oilseed crop for the production of renewable oil. Here, a sRNA library of mature seeds and three mRNA libraries from three different seed development stages were generated by deep sequencing to identify and characterize the miRNAs and pre-miRNAs of J. curcas. Computational analysis was used for the identification of 180 conserved miRNAs and 41 precursors (pre-miRNAs) as well as 16 novel pre-miRNAs. The predicted miRNA target genes are involved in a broad range of physiological functions, including cellular structure, nuclear function, translation, transport, hormone synthesis, defense, and lipid metabolism. Some pre-miRNA and miRNA targets vary in abundance between the three stages of seed development. A search for sequences that produce siRNA was performed, and the results indicated that J. curcas siRNAs play a role in nuclear functions, transport, catalytic processes and disease resistance. This study presents the first large scale identification of J. curcas miRNAs and their targets in mature seeds based on deep sequencing, and it contributes to a functional understanding of these miRNAs.

Identification of microRNAs from Eugenia uniflora by high-throughput sequencing and bioinformatics analysis.

Background microRNAs or miRNAs are small non-coding regulatory RNAs that play important functions in the regulation of gene expression at the post-transcriptional level by targeting mRNAs for degradation or inhibiting protein translation. Eugenia uniflora is a plant native to tropical America with pharmacological and ecological importance, and there have been no previous studies concerning its gene expression and regulation. To date, no miRNAs have been reported in Myrtaceae species. Results Small RNA and RNA-seq libraries were constructed to identify miRNAs and pre-miRNAs in Eugenia uniflora. Solexa technology was used to perform high throughput sequencing of the library, and the data obtained were analyzed using bioinformatics tools. From 14,489,131 small RNA clean reads, we obtained 1,852,722 mature miRNA sequences representing 45 conserved families that have been identified in other plant species. Further analysis using contigs assembled from RNA-seq allowed the prediction of secondary structures of 25 known and 17 novel pre-miRNAs. The expression of twenty-seven identified miRNAs was also validated using RT-PCR assays. Potential targets were predicted for the most abundant mature miRNAs in the identified pre-miRNAs based on sequence homology. Conclusions This study is the first large scale identification of miRNAs and their potential targets from a species of the Myrtaceae family without genomic sequence resources. Our study provides more information about the evolutionary conservation of the regulatory network of miRNAs in plants and highlights species-specific miRNAs.

Molecular evolution of the lysophosphatidic acid acyltransferase (LPAAT) gene family.

Lysophosphatidic acid acyltransferases (LPAATs) perform an essential cellular function by controlling the production of phosphatidic acid (PA), a key intermediate in the synthesis of membrane, signaling and storage lipids. Although LPAATs have been extensively explored by functional and biotechnological studies, little is known about their molecular evolution and diversification. We performed a genome-wide analysis using data from several plants and animals, as well as other eukaryotic and prokaryotic species, to identify LPAAT genes and analyze their evolutionary history. We used phylogenetic and molecular evolution analysis to test the hypothesis of distinct origins for these genes. The reconstructed phylogeny supported the ancient origin of some isoforms (plant LPAAT1 and LPAATB; animal AGPAAT1/2), while others emerged more recently (plant LPAAT2/3/4/5; AGPAAT3/4/5/8). Additionally, the hypothesis of endosymbiotic origin of the plastidic isoform LPAAT1 was confirmed. LPAAT genes from plants and animals mainly experienced strong purifying selection pressures with limited functional divergence after the species-specific duplications. Gene expression analyses of LPAAT isoforms in model plants demonstrated distinct LPAAT expression patterns in these organisms. The results showed that distinct origins followed by diversification of the LPAAT genes shaped the evolution of TAG biosynthesis. The expression pattern of individual genes may be responsible for adaptation into multiple ecological niches.

Comprehensive selection of reference genes for quantitative gene expression analysis during seed development in Brassica napus

Key messageMicroRNAs have higher expression stability than protein-coding genes inB. napusseeds and are therefore good reference genes for miRNA and mRNA RT-qPCR analysis.AbstractReverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) has become the “gold standard” to gain insight into function of genes. However, the accuracy of the technique depends on appropriate reference genes for quantification analysis in different experimental conditions. Accumulation of microRNAs (miRNAs) has also been studied by RT-qPCR, but there are no reference genes currently validated for normalization of Brassica napus miRNA expression data. In this study, we selected 43 B. napus miRNAs and 18 previously validated mRNA reference genes. The expression stability of the candidate reference genes was evaluated in different tissue samples (stages of seed development, flowers, and leaves) using geNorm, NormFinder, and RefFinder analysis. The best-ranked reference genes for expression studies during seed development (miR167-1_2, miR11-1, miR159-1 and miR168-1) were used to asses the expression of miR03-1. Since candidate miRNAs showed higher expression stability than protein-coding genes in most of the tested conditions, the expression profile of DGAT1 gene was compared when normalized by the four most stable miRNAs reference genes and by the four most stable mRNA reference genes. The expected expression pattern of DGAT1 during seed development was achieved with the use of miRNA as reference genes. In conclusion, the most stable miRNA reference genes can be employed in the normalization of RT-qPCR quantification of miRNAs and protein-coding genes. This work is the first to perform a comprehensive survey of the stability of miRNA reference genes in B. napus and provides guidelines to obtain more accurate RT-qPCR results in B. napus seeds studies.

Analysis of castor bean ribosome-inactivating proteins and their gene expression during seed development

Ribosome-inactivating proteins (RIPs) are enzymes that inhibit protein synthesis after depurination of a specific adenine in rRNA. The RIP family members are classified as type I RIPs that contain an RNA-N-glycosidase domain and type II RIPs that contain a lectin domain (B chain) in addition to the glycosidase domain (A chain). In this work, we identified 30 new plant RIPs and characterized 18 Ricinus communis RIPs. Phylogenetic and functional divergence analyses indicated that the emergence of type I and II RIPs probably occurred before the monocot/eudicot split. We also report the expression profiles of 18 castor bean genes, including those for ricin and agglutinin, in five seed stages as assessed by quantitative PCR. Ricin and agglutinin were the most expressed RIPs in developing seeds although eight other RIPs were also expressed. All of the RIP genes were most highly expressed in the stages in which the endosperm was fully expanded. Although the reason for the large expansion of RIP genes in castor beans remains to be established, the differential expression patterns of the type I and type II members reinforce the existence of biological functions other than defense against predators and herbivory.

New insights on the evolution of Leafy cotyledon1 (LEC1) type genes in vascular plants

NF-Y is a conserved oligomeric transcription factor found in all eukaryotes. In plants, this regulator evolved with a broad diversification of the genes coding for its three subunits (NF-YA, NF-YB and NF-YC). The NF-YB members can be divided into Leafy Cotyledon1 (LEC1) and non-LEC1 types. Here we presented a comparative genomic study using phylogenetic analyses to validate an evolutionary model for the origin of LEC-type genes in plants and their emergence from non-LEC1-type genes. We identified LEC1-type members in all vascular plant genomes, but not in amoebozoa, algae, fungi, metazoa and non-vascular plant representatives, which present exclusively non-LEC1-type genes as constituents of their NF-YB subunits. The non-synonymous to synonymous nucleotide substitution rates (Ka/Ks) between LEC1 and non-LEC1-type genes indicate the presence of positive selection acting on LEC1-type members to the fixation of LEC1-specific amino acid residues. The phylogenetic analyses demonstrated that plant LEC1-type genes are evolutionary divergent from the non-LEC1-type genes of plants, fungi, amoebozoa, algae and animals. Our results point to a scenario in which LEC1-type genes have originated in vascular plants after gene expansion in plants. We suggest that processes of neofunctionalization and/or subfunctionalization were responsible for the emergence of a versatile role for LEC1-type genes in vascular plants, especially in seed plants. LEC1-type genes besides being phylogenetic divergent also present different expression profile when compared with non-LEC1-type genes. Altogether, our data provide new insights about the LEC1 and non-LEC1 evolutionary relationship during the vascular plant evolution.

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

The miRNAs play important roles in regulation of gene expression at the post-transcriptional level. A small RNA and RNA-seq of libraries were constructed to identify miRNAs in Vriesea carinata, a native bromeliad species from Brazilian Atlantic Rainforest. Illumina technology was used to perform high throughput sequencing and data was analyzed using bioinformatics tools. We obtained 2,191,509 mature miRNAs sequences representing 54 conserved families in plant species. Further analysis allowed the prediction of secondary structures for 19 conserved and 16 novel miRNAs. Potential targets were predicted from pre-miRNAs by sequence homology and validated using RTqPCR approach. This study provides the first identification of miRNAs and their potential targets of a bromeliad species.

Rice peroxisomal ascorbate peroxidase knockdown affects ROS signaling and triggers early leaf senescence

H2O2, which is continually produced by aerobic metabolism, is a cytotoxic molecule when in high levels. However, low levels can act as a signaling molecule able to regulate the expression of stress responses, senescence, programmed cell death, plant growth, and development. Ascorbate peroxidase (APX) enzyme plays an essential role in the control of intracellular H2O2 levels. Here, the function of a gene encoding a peroxisomal APX (OsAPX4) from rice (Oryza sativa L.) was studied. OsAPX4 gene expression can be detected in roots and panicles, but the highest expression level occurs in leaves. Silencing of OsAPX4 and OsAPX3 expression in RNAiOsAPX4 did not affect the growth of plants under growth chamber conditions, but aging transgenic plants interestingly displayed an early senescence phenotype. Leaf fragments from silenced plants were also more sensitive to induced senescence conditions. RNAiOsAPX4 plants did not present detectable changes in intracellular H2O2 levels, but biochemical analyses showed that transgenic plants displayed some decreased APX activity in the chloroplastic fraction. Also, the peroxisomal enzyme glycolate oxidase exhibited lower activity, whereas catalase activity was similar to non-transformed rice. The results imply that OsAPX4 gene has an important role in leaf senescence pathway mediated by ROS signaling.