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 MESSAGE: MicroRNAs have higher expression stability than protein-coding genes in B. napus seeds and are therefore good reference genes for miRNA and mRNA RT-qPCR analysis. Reverse 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.

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.

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.

Two GCC boxes and AP2/ERF-domain transcription factor ORA59 in jasmonate/ethylene-mediated activation of the PDF1.2 promoter in Arabidopsis.

Plant defense against microbial pathogens depends on the action of several endogenously produced hormones, including jasmonic acid (JA) and ethylene (ET). In defense against necrotrophic pathogens, the JA and ET signaling pathways synergize to activate a specific set of defense genes including PLANT DEFENSIN1.2 (PDF1.2). The APETALA2/Ethylene Response Factor (AP2/ERF)-domain transcription factor ORA59 acts as the integrator of the JA and ET signaling pathways and is the key regulator of JA- and ET-responsive PDF1.2 expression. The present study was aimed at the identification of elements in the PDF1.2 promoter conferring the synergistic response to JA/ET and interacting with ORA59. We show that the PDF1.2 promoter was activated synergistically by JA and the ET-releasing agent ethephon due to the activity of two GCC boxes. ORA59 bound in vitro to these GCC boxes and trans-activated the PDF1.2 promoter in transient assays via these two boxes. Using the chromatin immunoprecipitation technique we were able to show that ORA59 bound the PDF1.2 promoter in vivo. Finally, we show that a tetramer of a single GCC box conferred JA/ethephon-responsive expression, demonstrating that the JA and ET signaling pathways converge to a single type of GCC box. Therefore ORA59 and two functionally equivalent GCC box binding sites form the module that enables the PDF1.2 gene to respond synergistically to simultaneous activation of the JA and ET signaling pathways.

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.

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.

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.