Andreia Carina Turchetto-Zolet

Evolutionary view of acyl-CoA diacylglycerol acyltransferase (DGAT), a key enzyme in neutral lipid biosynthesis

BackgroundTriacylglycerides (TAGs) are a class of neutral lipids that represent the most important storage form of energy for eukaryotic cells. DGAT (acyl-CoA: diacylglycerol acyltransferase; EC 2.3.1.20) is a transmembrane enzyme that acts in the final and committed step of TAG synthesis, and it has been proposed to be the rate-limiting enzyme in plant storage lipid accumulation. In fact, two different enzymes identified in several eukaryotic species, DGAT1 and DGAT2, are the main enzymes responsible for TAG synthesis. These enzymes do not share high DNA or protein sequence similarities, and it has been suggested that they play non-redundant roles in different tissues and in some species in TAG synthesis. Despite a number of previous studies on the DGAT1 and DGAT2 genes, which have emphasized their importance as potential obesity treatment targets to increase triacylglycerol accumulation, little is known about their evolutionary timeline in eukaryotes. The goal of this study was to examine the evolutionary relationship of the DGAT1 and DGAT2 genes across eukaryotic organisms in order to infer their origin.ResultsWe have conducted a broad survey of fully sequenced genomes, including representatives of Amoebozoa, yeasts, fungi, algae, musses, plants, vertebrate and invertebrate species, for the presence of DGAT1 and DGAT2 gene homologs. We found that the DGAT1 and DGAT2 genes are nearly ubiquitous in eukaryotes and are readily identifiable in all the major eukaryotic groups and genomes examined. Phylogenetic analyses of the DGAT1 and DGAT2 amino acid sequences revealed evolutionary partitioning of the DGAT protein family into two major DGAT1 and DGAT2 clades. Protein secondary structure and hydrophobic-transmembrane analysis also showed differences between these enzymes. The analysis also revealed that the MGAT2 and AWAT genes may have arisen from DGAT2 duplication events.ConclusionsIn this study, we identified several DGAT1 and DGAT2 homologs in eukaryote taxa. Overall, the data show that DGAT1 and DGAT2 are present in most eukaryotic organisms and belong to two different gene families. The phylogenetic and evolutionary analyses revealed that DGAT1 and DGAT2 evolved separately, with functional convergence, despite their wide molecular and structural divergence.

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.

The evolution of pyrroline-5-carboxylate synthase in plants: a key enzyme in proline synthesis

Many plants synthesize and accumulate proline in response to osmotic stress conditions. A central enzyme in the proline biosynthesis is the bifunctional enzyme Δ1-pyrroline-5-carboxylate synthase (P5CS) that includes two functional catalytic domains: the γ-glutamyl kinase and the glutamic-γ-semialdehyde dehydrogenase. This enzyme catalyzes the first two steps of the proline biosynthetic pathway and plays a central role in the regulation of this process in plants. To determine the evolutionary events that occurred in P5CS genes, partial sequences from four Neotropical trees were cloned and compared to those of other plant taxa. Molecular phylogenetic analysis indicated that P5CS duplication events have occurred several times following the emergence of flowering plants and at different frequencies throughout the evolution of monocots and dicots. Despite the high number of conserved residues in plant P5CS sequences, positive selection was observed at different regions of P5CS paralogous genes and also when dicots and monocots were contrasted.

Heavy metal‐associated isoprenylated plant protein (HIPP): characterization of a family of proteins exclusive to plants

Metallochaperones are key proteins for the safe transport of metallic ions inside the cell. HIPPs (heavy metal‐associated isoprenylated plant proteins) are metallochaperones that contain a metal binding domain (HMA) and a C–terminal isoprenylation motif. In this study, we provide evidence that proteins of this family are found only in vascular plants and may be separated into five distinct clusters. HIPPs may be involved in (a) heavy metal homeostasis and detoxification mechanisms, especially those involved in cadmium tolerance, (b) transcriptional responses to cold and drought, and (c) plant–pathogen interactions. In particular, our results show that the rice (Oryza sativa) HIPP OsHIPP41 gene is highly expressed in response to cold and drought stresses, and its product is localized in the cytosol and the nucleus. The results suggest that HIPPs play an important role in the development of vascular plants and in plant responses to environmental changes.

Genome-wide annotation of the soybean WRKY family and functional characterization of genes involved in response to Phakopsora pachyrhizi infection

BackgroundMany previous studies have shown that soybean WRKY transcription factors are involved in the plant response to biotic and abiotic stresses. Phakopsora pachyrhizi is the causal agent of Asian Soybean Rust, one of the most important soybean diseases. There are evidences that WRKYs are involved in the resistance of some soybean genotypes against that fungus. The number of WRKY genes already annotated in soybean genome was underrepresented. In the present study, a genome-wide annotation of the soybean WRKY family was carried out and members involved in the response to P. pachyrhizi were identified.ResultsAs a result of a soybean genomic databases search, 182 WRKY-encoding genes were annotated and 33 putative pseudogenes identified. Genes involved in the response to P. pachyrhizi infection were identified using superSAGE, RNA-Seq of microdissected lesions and microarray experiments. Seventy-five genes were differentially expressed during fungal infection. The expression of eight WRKY genes was validated by RT-qPCR. The expression of these genes in a resistant genotype was earlier and/or stronger compared with a susceptible genotype in response to P. pachyrhizi infection. Soybean somatic embryos were transformed in order to overexpress or silence WRKY genes. Embryos overexpressing a WRKY gene were obtained, but they were unable to convert into plants. When infected with P. pachyrhizi, the leaves of the silenced transgenic line showed a higher number of lesions than the wild-type plants.ConclusionsThe present study reports a genome-wide annotation of soybean WRKY family. The participation of some members in response to P. pachyrhizi infection was demonstrated. The results contribute to the elucidation of gene function and suggest the manipulation of WRKYs as a strategy to increase fungal resistance in soybean plants.

Identification and in silico characterization of soybean trihelix-GT and bHLH transcription factors involved in stress responses

Environmental stresses caused by either abiotic or biotic factors greatly affect agriculture. As for soybean [Glycine max (L.) Merril], one of the most important crop species in the world, the situation is not different. In order to deal with these stresses, plants have evolved a variety of sophisticated molecular mechanisms, to which the transcriptional regulation of target-genes by transcription factors is crucial. Even though the involvement of several transcription factor families has been widely reported in stress response, there still is a lot to be uncovered, especially in soybean. Therefore, the objective of this study was to investigate the role of bHLH and trihelix-GT transcription factors in soybean responses to environmental stresses. Gene annotation, data mining for stress response, and phylogenetic analysis of members from both families are presented herein. At least 45 bHLH (from subgroup 25) and 63 trihelix-GT putative genes reside in the soybean genome. Among them, at least 14 bHLH and 11 trihelix-GT seem to be involved in responses to abiotic/biotic stresses. Phylogenetic analysis successfully clustered these with members from other plant species. Nevertheless, bHLH and trihelix-GT genes encompass almost three times more members in soybean than in Arabidopsis or rice, with many of these grouping into new clades with no apparent near orthologs in the other analyzed species. Our results represent an important step towards unraveling the functional roles of plant bHLH and trihelix-GT transcription factors in response to environmental cues.

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.

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.

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.