Frank Guzman

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.

De novo assembly of Eugenia uniflora L. transcriptome and identification of genes from the terpenoid biosynthesis pathway

Pitanga (Eugenia uniflora L.) is a member of the Myrtaceae family and is of particular interest due to its medicinal properties that are attributed to specialized metabolites with known biological activities. Among these molecules, terpenoids are the most abundant in essential oils that are found in the leaves and represent compounds with potential pharmacological benefits. The terpene diversity observed in Myrtaceae is determined by the activity of different members of the terpene synthase and oxidosqualene cyclase families. Therefore, the aim of this study was to perform a de novo assembly of transcripts from E. uniflora leaves and to annotation to identify the genes potentially involved in the terpenoid biosynthesis pathway and terpene diversity. In total, 72,742 unigenes with a mean length of 1048bp were identified. Of these, 43,631 and 36,289 were annotated with the NCBI non-redundant protein and Swiss-Prot databases, respectively. The gene ontology categorized the sequences into 53 functional groups. A metabolic pathway analysis with KEGG revealed 8,625 unigenes assigned to 141 metabolic pathways and 40 unigenes predicted to be associated with the biosynthesis of terpenoids. Furthermore, we identified four putative full-length terpene synthase genes involved in sesquiterpenes and monoterpenes biosynthesis, and three putative full-length oxidosqualene cyclase genes involved in the triterpenes biosynthesis. The expression of these genes was validated in different E. uniflora tissues.

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.

Identification and evaluation of reference genes for expression studies by RT-qPCR during embryonic development of the emerging model organism, Macrobrachium olfersii.

RT-qPCR is a sensitive and highly efficient technique that is widely used in gene expression analysis and to provide insight into the molecular mechanisms underlying embryonic development. The freshwater prawn, Macrobrachium olfersii is an emerging model organism, but, the stable reference genes of this species need to be identified and validated for RT-qPCR analysis. Thus, the aim of this study was to evaluate the expression stability of six genes (β-act, GAPDH, EF-1α, RpL8, RpS6, AK) in embryos and in adult tissues (cerebral ganglia, muscle and hepatopancreas) of M. olfersii. The expression stabilities of these genes were evaluated using geNorm, NormFinder, BestKeeper, ΔCt method and integrated tool RefFinder. In the general ranking, RpL8 and RpS6 were the most stable genes in embryos, while RpS6 and RpL8 were the most stable in a combined adult tissue analysis. Analysis of the adult tissues revealed that β-act and AK were the most stable genes in cerebral ganglia, RpL8 and AK in muscle, and RpS6 and β-act in hepatopancreas. EF-1α and GAPDH were the least stable genes and as normalizer genes in RT-qPCR affected expression of the Distal-less gene during M. olfersii development. This study provides suitable reference genes for RT-qPCR analysis and allows future studies of the gene expression in M. olfersii for understanding the molecular mechanisms of their development. To our knowledge, this is the first published study that identifies and evaluates reference genes for RT-qPCR analysis in M. olfersii and could be useful as basis for evaluations of reference genes in other prawns.

Exploring developmental gene toolkit and associated pathways in a potential new model crustacean using transcriptomic analysis

The crustaceans are one of the largest, most diverse, and most successful groups of invertebrates. The diversity among the crustaceans is also reflected in embryonic development models. However, the molecular genetics that regulates embryonic development is not known in those crustaceans that have a short germ-band development with superficial cleavage, such as Macrobrachium olfersi. This species is a freshwater decapod and has great potential to become a model for developmental biology, as well as for evolutionary and environmental studies. To obtain sequence data of M. olfersi from an embryonic developmental perspective, we performed de novo assembly and annotation of the embryonic transcriptome. Using a pooling strategy of total RNA, paired-end Illumina sequencing, and assembly with multiple k-mers, a total of 25,636,097 pair reads were generated. In total, 99,751 unigenes were identified, and 20,893 of these returned a Blastx hit. KEGG pathway analysis mapped a total of 6866 unigenes related to 129 metabolic pathways. In general, 21,845 unigenes were assigned to gene ontology (GO) categories: molecular function (19,604), cellular components (10,254), and biological processes (13,841). Of these, 2142 unigenes were assigned to the developmental process category. More specifically, a total of 35 homologs of embryonic development toolkit genes were identified, which included maternal effect (one gene), gap (six), pair-rule (six), segment polarity (seven), Hox (four), Wnt (eight), and dorsoventral patterning genes (three). In addition, genes of developmental pathways were found, including TGF-β, Wnt, Notch, MAPK, Hedgehog, Jak-STAT, VEGF, and ecdysteroid-inducible nuclear receptors. RT-PCR analysis of eight genes related to embryonic development from gastrulation to late morphogenesis/organogenesis confirmed the applicability of the transcriptome analysis.

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.