Pervasive transcription of plant organelle genomes: functional noncoding transcriptomes?

Plant mitochondrial and plastid genomes typically show pervasive, genome-wide transcription. Little is known, however, about the utility of organelle noncoding RNAs, which often make up most of the transcriptome. Here, we suggest that long-read sequencing data combined with dedicated RNA databases could help identify putative functional organelle noncoding transcripts.

Functional Characterization of ent-Copalyl Diphosphate Synthase and Kaurene Synthase Genes from Coffea arabica L

The biochemical profile of coffee beans translates directly into quality traits, nutraceutical and health promoting properties of the coffee beverage. Ent-kaurene is the ubiquitous precursor for gibberellin biosynthesis in plants, but it also serves as an intermediate in specialized (i.e., secondary) diterpenoid metabolism that leads to a diversity of more than 1,000 different metabolites. Nutraceutical effects on human health attributed to diterpenes include antioxidant, anticarcinogenic, and anti-inflammatory properties. Cafestol (CAF) and kahweol (KAH) are two diterpenes found exclusively in the Coffea genus. Our objective was to identify and functionally characterize genes involved in the central step of ent-kaurene production. We identified 17 putative terpene synthase genes in the transcriptome of Coffea arabica. Two ent-copalyl diphosphate synthase (CaCPS) and three kaurene synthase (CaKS) were selected and manually annotated. Transcript expression profiles of CaCPS1 and CaKS3 best matched the CAF and KAH metabolite profiles in different tissues. CaCPS1 and CaKS3 proteins were heterologously expressed and functionally characterized. CaCPS1 catalyzes the cyclization of geranylgeranyl diphosphate (GGPP) to ent-copalyl diphosphate (ent-CPP), which is converted to ent-kaurene by CaKS3. Knowledge about the central steps of diterpene formation in coffee provides a foundation for future characterization of the subsequent enzymes involved in CAF and KAH biosynthesis.

Genome-Wide Analysis of Lipoxygenase (LOX) Genes in Angiosperms

Lipoxygenases (LOXs) are enzymes that catalyze the addition of an oxygen molecule to unsaturated fatty acids, thus forming hydroperoxides. In plants, these enzymes are encoded by a multigene family found in several organs with varying activity patterns, by which they are classified as LOX9 or LOX13. They are involved in several physiological functions, such as growth, fruit development, and plant defense. Despite several studies on genes of the LOX family in plants, most studies are restricted to a single species or a few closely related species. This study aimed to analyze the diversity, evolution, and expression of LOX genes in angiosperm species. We identified 247 LOX genes among 23 species of angiosperms and basal plants. Phylogenetic analyses identified clades supporting LOX13 and two main clades for LOX9: LOX9_A and LOX9_B. Eudicot species such as Tarenaya hassleriana, Capsella rubella, and Arabidopsis thaliana did not present LOX9_B genes; however, LOX9_B was present in all monocots used in this study. We identified that there were potential new subcellular localization patterns and conserved residues of oxidation for LOX9 and LOX13 yet unexplored. In summary, our study provides a basis for the further functional and evolutionary study of lipoxygenases in angiosperms.

Genome-wide association study for resistance to Pseudomonas syringae pv. garcae in Coffea arabica

Bacteria halo blight (BHB), a coffee plant disease caused by Pseudomonas syringae pv. garcae, has been gaining importance in producing mountain regions and mild temperatures areas as well as in coffee nurseries. Most Coffea arabica cultivars are susceptible to this disease. In contrast, a great source of genetic diversity and resistance to BHB are found in C. arabica Ethiopian accessions. Aiming to identify quantitative trait nucleotides (QTNs) associated with resistance to BHB and the influence of these genomic regions during the domestication of C. arabica, we conducted an analysis of population structure and a Genome-Wide Association Study (GWAS). For this, we used genotyping by sequencing (GBS) and phenotyping for resistance to BHB of a panel with 120 C. arabica Ethiopian accessions from a historical FAO collection, 11 C. arabica cultivars, and the BA-10 genotype. Population structure analysis based on single-nucleotide polymorphisms (SNPs) markers showed that the 132 accessions are divided into 3 clusters: most wild Ethiopian accessions, domesticated Ethiopian accessions, and cultivars. GWAS, using the single-locus model MLM and the multi-locus models mrMLM, FASTmrMLM, FASTmrEMMA, and ISIS EM-BLASSO, identified 11 QTNs associated with resistance to BHB. Among these QTNs, the four with the highest values of association for resistance to BHB are linked to g000 (Chr_0_434_435) and g010741 genes, which are predicted to encode a serine/threonine-kinase protein and a nucleotide binding site leucine-rich repeat (NBS-LRR), respectively. These genes displayed a similar transcriptional downregulation profile in a C. arabica susceptible cultivar and in a C. arabica cultivar with quantitative resistance, when infected with P. syringae pv. garcae. However, peaks of upregulation were observed in a C. arabica cultivar with qualitative resistance, for both genes. Our results provide SNPs that have potential for application in Marker Assisted Selection (MAS) and expand our understanding about the complex genetic control of the resistance to BHB in C. arabica. In addition, the findings contribute to increasing understanding of the C. arabica domestication history.

Computational Analysis of Transposable Elements and CircRNAs in Plants

This chapter provides two main contributions: (1) a description of computational tools and databases used to identify and analyze transposable elements (TEs) and circRNAs in plants; and (2) data analysis on public TE and circRNA data. Our goal is to highlight the primary information available in the literature on circular noncoding RNAs and transposable elements in plants. The exploratory analysis performed on publicly available circRNA and TEs data help discuss four sequence features. Finally, we investigate the association on circRNAs:TE in plants in the model organism Arabidopsis thaliana.

An Atlas of Plant Transposable Elements

Advances in genomic sequencing have recently offered vast opportunities for biological exploration, unraveling the evolution and improving our understanding of Earth biodiversity. Due to distinct plant species characteristics in terms of genome size, ploidy and heterozygosity, transposable elements (TEs) are common characteristics of many genomes. TEs are ubiquitous and dispersed repetitive DNA sequences that frequently impact the evolution and composition of the genome, mainly due to their redundancy and rearrangements. For this study, we provided an atlas of TE data by employing an easy-to-use portal ( APTE website ). To our knowledge, this is the most extensive and standardized analysis of TEs in plant genomes. We evaluated 67 plant genomes assembled at chromosome scale, recovering a total of 49,802,023 TE records, representing a total of 47,992,091,043 (~47,62%) base pairs (bp) of the total genomic space. We observed that new types of TEs were identified and annotated compared to other data repositories. By establishing a standardized catalog of TE annotation on 67 genomes, new hypotheses, exploration of TE data and their influences on the genomes may allow a better understanding of their function and processes. All original code and an example of how we developed the TE annotation strategy is available on GitHub ( Extended data).

Hexanoic acid: a new potential substitute for copper-based agrochemicals against citrus canker

AIMS

The aim of the study is to evaluate hexanoic acid (HA) as an alternative to manage citrus canker.

METHODS AND RESULTS

The minimal growth inhibitory concentration of HA against Xanthomonas citri subsp. citri was determined at 2·15 mmol l^-1 using a respiratory activity assay. Growth curves at different pH values showed that growth inhibition was not due to media acidification induced by HA. The germination rate and root elongation of Lactuca sativa seeds exposed to different concentrations of HA (varying from 0·86 to 5·16 mmol l^-1 ) were assessed to screen for phytotoxicity. The acid exhibited low phytotoxicity for L. sativa at 1·29 and 2·58 mmol l^-1 . To evaluate the ability of HA to protect citrus against X. citri infection, leaves of Citrus sinensis were sprayed with the acid and subsequently challenged with X. citri. HA at 3·44 mmol l^-1 was able to protect citrus against infection, showing a reduction of three orders of magnitude in the number of citrus canker lesions per cm² when compared to the untreated negative control.

CONCLUSION

HA is a potential alternative to copper for citrus canker management.

SIGNIFICANCE AND IMPACT OF THE STUDY

HA inhibits X. citri growth, exhibits low phytotoxicity and is an alternative to copper for the protection of citrus plants against bacterial infection.

The urea transporter DUR3 is differentially regulated by abiotic and biotic stresses in coffee plants

The high costs of N fertilizers in the coffee production emphasizes the need to optimize fertilization practices and improve nitrogen use efficiency. Urea is widespread in nature, characterizing itself as a significant source of nitrogen for the growth and development of several organisms. Thus, the characterization of genes involved in urea transport in coffee plants is an important research topic for the sustainable production of this valuable cash crop. In the current study, we evaluated the expression of the DUR3 gene under abiotic and biotic stresses in coffee plants. Here, we show that the expression of a high-affinity urea transporter gene (CaDUR3) was up-regulated by N starvation in leaves and roots of two out of three C. arabica cultivars examined. Moreover, the CaDUR3 gene was differentially expressed in coffee plants under different abiotic and biotic stresses. In plants of cv. IAPAR59, CaDUR3 showed an increased expression in leaves after exposure to water deficit and heat stress, while it was downregulated in plants under salinity. Upon infection with H. vastatrix (coffee rust), the CaDUR3 was markedly up-regulated at the beginning of the infection process in the disease susceptible Catuaí Vermelho 99 in comparison with the resistant cultivar. These results indicate that besides urea acquisition and N-remobilization, CaDUR3 gene may be closely involved in the response to various stresses.

TERL: classification of transposable elements by convolutional neural networks

Transposable elements (TEs) are the most represented sequences occurring in eukaryotic genomes. Few methods provide the classification of these sequences into deeper levels, such as superfamily level, which could provide useful and detailed information about these sequences. Most methods that classify TE sequences use handcrafted features such as k-mers and homology-based search, which could be inefficient for classifying non-homologous sequences. Here we propose an approach, called transposable elements pepresentation learner (TERL), that preprocesses and transforms one-dimensional sequences into two-dimensional space data (i.e., image-like data of the sequences) and apply it to deep convolutional neural networks. This classification method tries to learn the best representation of the input data to classify it correctly. We have conducted six experiments to test the performance of TERL against other methods. Our approach obtained macro mean accuracies and F1-score of 96.4% and 85.8% for superfamilies and 95.7% and 91.5% for the order sequences from RepBase, respectively. We have also obtained macro mean accuracies and F1-score of 95.0% and 70.6% for sequences from seven databases into superfamily level and 89.3% and 73.9% for the order level, respectively. We surpassed accuracy, recall and specificity obtained by other methods on the experiment with the classification of order level sequences from seven databases and surpassed by far the time elapsed of any other method for all experiments. Therefore, TERL can learn how to predict any hierarchical level of the TEs classification system and is about 20 times and three orders of magnitude faster than TEclass and PASTEC, respectively https://github.com/muriloHoracio/TERL. Contact:murilocruz@alunos.utfpr.edu.br.

mirtronDB: a mirtron knowledge base

Motivation

Mirtrons arise from short introns with atypical cleavage by using the splicing mechanism. In the current literature, there is no repository centralizing and organizing the data available to the public. To fill this gap, we developed mirtronDB, the first knowledge database dedicated to mirtron, and it is available at http://mirtrondb.cp.utfpr.edu.br/. MirtronDB currently contains a total of 1407 mirtron precursors and 2426 mirtron mature sequences in 18 species.

Results

Through a user-friendly interface, users can now browse and search mirtrons by organism, organism group, type and name. MirtronDB is a specialized resource that provides free and user-friendly access to knowledge on mirtron data.

Availability and implementation

MirtronDB is available at http://mirtrondb.cp.utfpr.edu.br/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Low-Copy Genes in Terpenoid Metabolism: The Evolution and Expression of MVK and DXR Genes in Angiosperms

Terpenoids are a diverse class of metabolites that impact plant metabolism in response to environmental cues. They are synthesized either via a predominantly cytosolic (MVA) pathway or a plastidic pathway (MEP). In Arabidopsis, several enzymes from the MVA and MEP pathways are encoded by gene families, excluding MVK and DXR, which are single-copy genes. In this study, we assess the diversity, evolution and expression of DXR and MVK genes in selected angiosperms and Coffea arabica in particular. Evolutionary analysis revealed that DXR and MVK underwent purifying selection, but the selection effect for DXR was stronger than it was for MVK. Digital gene expression (DGE) profile analysis of six species revealed that expression levels of MVK in flowers and roots were high, whereas for DXR peak values were observed in leaves. In C. arabica, both genes were highly expressed in flowers, and CaDXR was upregulated in response to methyl jasmonate. C. arabica DGE data were validated by assessing gene expression in selected organs, and by plants treated with hexanoic acid (Hx) using RT-qPCR. MVK expression was upregulated in roots treated with Hx. CaDXR was downregulated in leaves by Hx treatment in a genotype-specific manner, indicating a differential response to priming.

ceRNAs in plants: computational approaches and associated challenges for target mimic research

The competing endogenous RNA hypothesis has gained increasing attention as a potential global regulatory mechanism of microRNAs (miRNAs), and as a powerful tool to predict the function of many noncoding RNAs, including miRNAs themselves. Most studies have been focused on animals, although target mimic (TMs) discovery as well as important computational and experimental advances has been developed in plants over the past decade. Thus, our contribution summarizes recent progresses in computational approaches for research of miRNA:TM interactions. We divided this article in three main contributions. First, a general overview of research on TMs in plants is presented with practical descriptions of the available literature, tools, data, databases and computational reports. Second, we describe a common protocol for the computational and experimental analyses of TM. Third, we provide a bioinformatics approach for the prediction of TM motifs potentially cross-targeting both members within the same or from different miRNA families, based on the identification of consensus miRNA-binding sites from known TMs across sequenced genomes, transcriptomes and known miRNAs. This computational approach is promising because, in contrast to animals, miRNA families in plants are large with identical or similar members, several of which are also highly conserved. From the three consensus TM motifs found with our approach: MIM166, MIM171 and MIM159/319, the last one has found strong support on the recent experimental work by Reichel and Millar [Specificity of plant microRNA TMs: cross-targeting of mir159 and mir319. J Plant Physiol 2015;180:45-8]. Finally, we stress the discussion on the major computational and associated experimental challenges that have to be faced in future ceRNA studies.

Decarboxylation mechanisms of C4 photosynthesis in Saccharum spp.: increased PEPCK activity under water-limiting conditions

BACKGROUND

C4 plants have been classified into three subtypes based on the enzymes used to decarboxylate C4 acids in the bundle sheath cells (NADP-ME, NAD-ME and PEPCK pathways). Evidences indicate that, depending on environmental factors, C4 plants may exhibit a certain degree of flexibility in the use of the decarboxylation mechanisms. In this context, the objective was to extend the knowledge on the degree of flexibility between the pathways of decarboxylation in sugarcane, a NADP-ME species, at different levels of water deficit.

RESULTS

An experiment was carried out with two cultivars - RB92579 (tolerant to water deficit) and SP80-3280 (susceptible to water deficit) subjected to moderate level (- 1.5 to - 1.8 MPa), severe level (below - 2.0 MPa) and recovery (48 h after rehydration) and changes in the activities of the enzymes involved in the three C4 mechanisms and in gene expression were investigated. Our results showed that sugarcane uses the PEPCK pathway as a decarboxylation mechanism in addition to the NADP-ME, which was more evident under water deficit conditions for both cultivars.

CONCLUSIONS

The results obtained here, show that sugarcane increases the use of the PEPCK pathway as a decarboxylation mechanism, in addition to the NADP-ME pathway, under conditions of water deficit, particularly in the tolerant cultivar.

Functional annotation and distribution overview of RNA families in 27 Streptococcus agalactiae genomes

Background

Streptococcus agalactiae, also known as Group B Streptococcus (GBS), is a Gram-positive bacterium that colonizes the gastrointestinal and genitourinary tract of humans. This bacterium has also been isolated from various animals, such as fish and cattle. Non-coding RNAs (ncRNAs) can act as regulators of gene expression in bacteria, such as Streptococcus pneumoniae and Streptococcus pyogenes. However, little is known about the genomic distribution of ncRNAs and RNA families in S. agalactiae.

Results

Comparative genome analysis of 27 S. agalactiae strains showed more than 5 thousand genomic regions identified and classified as Core, Exclusive, and Shared genome sequences. We identified 27 to 89 RNA families per genome distributed over these regions, from these, 25 were in Core regions while Shared and Exclusive regions showed variations amongst strains. We propose that the amount and type of ncRNA present in each genome can provide a pattern to contribute in the identification of the clonal types.

Conclusions

The identification of RNA families provides an insight over ncRNAs, sRNAs and ribozymes function, that can be further explored as targets for antibiotic development or studied in gene regulation of cellular processes. RNA families could be considered as markers to determine infection capabilities of different strains. Lastly, pan-genome analysis of GBS including the full range of functional transcripts provides a broader approach in the understanding of this pathogen.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4951-z) contains supplementary material, which is available to authorized users.

PlaNC-TE: a comprehensive knowledgebase of non-coding RNAs and transposable elements in plants

Transposable elements (TEs) play an essential role in the genetic variability of eukaryotic species. In plants, they may comprise up to 90% of the total genome. Non-coding RNAs (ncRNAs) are known to control gene expression and regulation. Although the relationship between ncRNAs and TEs is known, obtaining the organized data for sequenced genomes is not straightforward. In this study, we describe the PlaNC-TE (http://planc-te.cp.utfpr.edu.br), a user-friendly portal harboring a knowledgebase created by integrating and analysing plant ncRNA-TE data. We identified a total of 14 350 overlaps between ncRNAs and TEs in 40 plant genomes. The database allows users to browse, search and download all ncRNA and TE data analysed. Overall, PlaNC-TE not only organizes data and provides insights about the relationship between ncRNA and TEs in plants but also helps improve genome annotation strategies. Moreover, this is the first database to provide resources to broadly investigate functions and mechanisms involving TEs and ncRNAs in plants.

Genome-wide identification, classification and transcriptional analysis of nitrate and ammonium transporters in Coffea

Nitrogen (N) is quantitatively the main nutrient required by coffee plants, with acquisition mainly by the roots and mostly exported to coffee beans. Nitrate (NO3-) and ammonium (NH4+) are the most important inorganic sources for N uptake. Several N transporters encoded by different gene families mediate the uptake of these compounds. They have an important role in source preference for N uptake in the root system. In this study, we performed a genome-wide analysis, including in silico expression and phylogenetic analyses of AMT1, AMT2, NRT1/PTR, and NRT2 transporters in the recently sequenced Coffea canephora genome. We analyzed the expression of six selected transporters in Coffea arabica roots submitted to N deficiency. N source preference was also analyzed in C. arabica using isotopes. C. canephora N transporters follow the patterns observed for most eudicots, where each member of the AMT and NRT families has a particular role in N mobilization, and where some of these are modulated by N deficiency. Despite the prevalence of putative nitrate transporters in the Coffea genome, ammonium was the preferential inorganic N source for N-starved C. arabica roots. This data provides an important basis for fundamental and applied studies to depict molecular mechanisms involved in N uptake in coffee trees.

Genome-wide analysis of EgEVE_1, a transcriptionally active endogenous viral element associated to small RNAs in Eucalyptus genomes

Endogenous viral elements (EVEs) are the result of heritable horizontal gene transfer from viruses to hosts. In the last years, several EVE integration events were reported in plants by the exponential availability of sequenced genomes. Eucalyptus grandis is a forest tree species with a sequenced genome that is poorly studied in terms of evolution and mobile genetic elements composition. Here we report the characterization of E. grandis endogenous viral element 1 (EgEVE_1), a transcriptionally active EVE with a size of 5,664 bp. Phylogenetic analysis and genomic distribution demonstrated that EgEVE_1 is a newly described member of the Caulimoviridae family, distinct from the recently characterized plant Florendoviruses. Genomic distribution of EgEVE_1 and Florendovirus is also distinct. EgEVE_1 qPCR quantification in Eucalyptus urophylla suggests that this genome has more EgEVE_1 copies than E. grandis. EgEVE_1 transcriptional activity was demonstrated by RT-qPCR in five Eucalyptus species and one intrageneric hybrid. We also identified that Eucalyptus EVEs can generate small RNAs (sRNAs),that might be involved in de novo DNA methylation and virus resistance. Our data suggest that EVE families in Eucalyptus have distinct properties, and we provide the first comparative analysis of EVEs in Eucalyptus genomes.

Differentially Accumulated Proteins in Coffea arabica Seeds during Perisperm Tissue Development and Their Relationship to Coffee Grain Size

Coffee is one of the most important crops for developing countries. Coffee classification for trading is related to several factors, including grain size. Larger grains have higher market value then smaller ones. Coffee grain size is determined by the development of the perisperm, a transient tissue with a highly active metabolism, which is replaced by the endosperm during seed development. In this study, a proteomics approach was used to identify differentially accumulated proteins during perisperm development in two genotypes with regular (IPR59) and large grain sizes (IPR59-Graudo) in three developmental stages. Twenty-four spots were identified by MALDI-TOF/TOF-MS, corresponding to 15 proteins. We grouped them into categories as follows: storage (11S), methionine metabolism, cell division and elongation, metabolic processes (mainly redox), and energy. Our data enabled us to show that perisperm metabolism in IPR59 occurs at a higher rate than in IPR59-Graudo, which is supported by the accumulation of energy and detoxification-related proteins. We hypothesized that grain and fruit size divergences between the two coffee genotypes may be due to the comparatively earlier triggering of seed development processes in IPR59. We also demonstrated for the first time that the 11S protein is accumulated in the coffee perisperm.

Transcriptionally active LTR retrotransposons in Eucalyptus genus are differentially expressed and insertionally polymorphic

BACKGROUND

In Eucalyptus genus, studies on genome composition and transposable elements (TEs) are particularly scarce. Nearly half of the recently released Eucalyptus grandis genome is composed by retrotransposons and this data provides an important opportunity to understand TE dynamics in Eucalyptus genome and transcriptome.

RESULTS

We characterized nine families of transcriptionally active LTR retrotransposons from Copia and Gypsy superfamilies in Eucalyptus grandis genome and we depicted genomic distribution and copy number in two Eucalyptus species. We also evaluated genomic polymorphism and transcriptional profile in three organs of five Eucalyptus species. We observed contrasting genomic and transcriptional behavior in the same family among different species. RLC_egMax_1 was the most prevalent family and RLC_egAngela_1 was the family with the lowest copy number. Most families of both superfamilies have their insertions occurring <3 million years, except one Copia family, RLC_egBianca_1. Protein theoretical models suggest different properties between Copia and Gypsy domains. IRAP and REMAP markers suggested genomic polymorphisms among Eucalyptus species. Using EST analysis and qRT-PCRs, we observed transcriptional activity in several tissues and in all evaluated species. In some families, osmotic stress increases transcript values.

CONCLUSION

Our strategy was successful in isolating transcriptionally active retrotransposons in Eucalyptus, and each family has a particular genomic and transcriptional pattern. Overall, our results show that retrotransposon activity have differentially affected genome and transcriptome among Eucalyptus species.

Galactinol synthase transcriptional profile in two genotypes of Coffea canephora with contrasting tolerance to drought

Increased synthesis of galactinol and raffinose family oligosaccharides (RFOs) has been reported in vegetative tissues in response to a range of abiotic stresses. In this work, we evaluated the transcriptional profile of a Coffea canephora galactinol synthase gene (CcGolS1) in two clones that differed in tolerance to water deficit in order to assess the contribution of this gene to drought tolerance. The expression of CcGolS1 in leaves was differentially regulated by water deficit, depending on the intensity of stress and the genotype. In clone 109A (drought-susceptible), the abundance of CcGolS1 transcripts decreased upon exposure to drought, reaching minimum values during recovery from severe water deficit and stress. In contrast, CcGolS1 gene expression in clone 14 (drought-tolerant) was stimulated by water deficit. Changes in galactinol and RFO content did not correlate with variation in the steady-state transcript level. However, the magnitude of increase in RFO accumulation was higher in the tolerant cultivar, mainly under severe water deficit. The finding that the drought-tolerant coffee clone showed enhanced accumulation of CcGolS1 transcripts and RFOs under water deficit suggests the possibility of using this gene to improve drought tolerance in this important crop.

Identification of the transcriptionally active cytochrome P450 repertoire in Coffea arabica

Cytochrome P450s (P450s) comprise a gene superfamily encoding enzymes that are involved in diverse plant metabolic pathways that produce primary and secondary metabolites such as phenylpropanoids, terpenoids, nitrogen-containing compounds, and plant hormones. They comprise one of the most diverse gene families in plant evolution. Although there are many studies that aim to characterize P450s in plants, there is no report on the characterization of this superfamily in Coffea arabica, where they might be related to plant tolerance to biotic and abiotic stresses, as well as aroma-related compounds. In this study, we report the characterization and annotation of 87 putative P450s from C. arabica obtained from the Brazilian Coffee Genome Project and describe their transcriptional pattern in different tissues and coffee organs. To validate our approach, we measured the transcriptional profile of the CaCYP81D8_1 gene by quantitative polymerase chain reaction in leaves, flowers, and fruits. This study is the first effort to present and analyze the P450 superfamily in C. arabica, which may assist in understanding the chemical diversity of coffee secondary metabolites.

Terminal-repeat retrotransposons with GAG domain in plant genomes: a new testimony on the complex world of transposable elements

A novel structure of nonautonomous long terminal repeat (LTR) retrotransposons called terminal repeat with GAG domain (TR-GAG) has been described in plants, both in monocotyledonous, dicotyledonous and basal angiosperm genomes. TR-GAGs are relatively short elements in length (<4 kb) showing the typical features of LTR-retrotransposons. However, they carry only one open reading frame coding for the GAG precursor protein involved for instance in transposition, the assembly, and the packaging of the element into the virus-like particle. GAG precursors show similarities with both Copia and Gypsy GAG proteins, suggesting evolutionary relationships of TR-GAG elements with both families. Despite the lack of the enzymatic machinery required for their mobility, strong evidences suggest that TR-GAGs are still active. TR-GAGs represent ubiquitous nonautonomous structures that could be involved in the molecular diversities of plant genomes.

CoffeebEST: an integrated resource for Coffea spp expressed sequence tags

Coffee is one of the most important commodities in the world, and its production relies mainly on two species, Coffea arabica and Coffea canephora. Although there are diverse transcriptome datasets available for coffee trees, few research groups have exploited the potential knowledge contained in these data, especially with respect to fruit and seed development. Here, we present a comparative analysis of the transcriptomes of Coffea arabica and Coffea canephora with a focus on fruit development using publicly available expressed sequence tags (ESTs). Most of the fruit and seed EST data has been obtained from C. canephora. Therefore, we performed a fruit EST analysis of the 5 developmental stages of this species (18, 22, 30, 42, and 46 weeks after flowering) comprising 29,009 sequences. We compared C. canephora fruit ESTs to reference unigenes of C. canephora (7710 contigs and 8955 singletons) and C. arabica (15,656 contigs and 16,351 singletons). Additional analyses included functional annotation based on Gene Onthology, as well as an annotation using PlantCyc, a curated plant protein database. The Coffee Bean EST (CoffeebEST) is a public database available at http://bioinfo-02.cp.utfpr.edu.br/. This database represents an additional resource for the coffee scientific community, offering a user-friendly collection of information for non-specialists in coffee molecular biology to support experimental research on comparative and functional genomics.

Nitrogen starvation, salt and heat stress in coffee (Coffea arabica L.): identification and validation of new genes for qPCR normalization

Abiotic stresses are among the most important factors that affect food production. One important step to face these environmental challenges is the transcriptional modulation. Quantitative real-time PCR is a rapid, sensitive, and reliable method for the detection of mRNAs and it has become a powerful tool to mitigate plant stress tolerance; however, suitable reference genes are required for data normalization. Reference genes for coffee plants during nitrogen starvation, salinity and heat stress have not yet been reported. We evaluated the expression stability of ten candidate reference genes using geNorm PLUS, NormFinder, and BestKeeper softwares, in plants submitted to nitrogen starvation, salt and heat stress. EF1, EF1α, GAPDH, MDH, and UBQ10 were ranked as the most stable genes in all stresses and software analyses, while RPL39 and RPII were classified as the less reliable references. For reference gene validation, the transcriptional pattern of a Coffea non-symbiotic hemoglobin (CaHb1) was analyzed using the two new recommended and the most unstable gene references for normalization. The most unstable gene may lead to incorrect interpretation of CaHb1 transcriptional analysis. Here, we recommend two new reference genes in Coffea for use in data normalization in abiotic stresses: MDH and EF1.

Construction and characterization of a BAC library from the Coffea arabica genotype Timor Hybrid CIFC 832/2

Most of the world's coffee production originates from Coffea arabica, an allotetraploid species with low genetic diversity and for which few genomic resources are available. Genomic libraries with large DNA fragment inserts are useful tools for the study of plant genomes, including the production of physical maps, integration studies of physical and genetic maps, genome structure analysis and gene isolation by positional cloning. Here, we report the construction and characterization of a Bacterial Artificial Chromosome (BAC) library from C. arabica Timor Hybrid CIFC 832/2, a parental genotype for several modern coffee cultivars. The BAC library consists of 56,832 clones with an average insert size of 118 kb, which represents a dihaploid genome coverage of five to sixfold. The content of organellar DNA was estimated at 1.04 and 0.5 % for chloroplast and mitochondrial DNA, respectively. The BAC library was screened for the NADPH-dependent mannose-6-phosphate reductase gene (CaM6PR) with markers positioned on four linkage groups of a partial C. arabica genetic map. A mixed approach using PCR and membrane hybridization of BAC pools allowed for the discovery of nine BAC clones with the CaM6PR gene and 53 BAC clones that were anchored to the genetic map with simple sequence repeat markers. This library will be a useful tool for future studies on comparative genomics and the identification of genes and regulatory elements controlling major traits in this economically important crop species.

The accumulation of endogenous proline induces changes in gene expression of several antioxidant enzymes in leaves of transgenic Swingle citrumelo

Plant exposure to abiotic stresses leads to an accumulation of reactive oxygen species with the concomitant increase in antioxidant defense mechanisms. Previous studies showed that exogenous application of proline mitigate the deleterious effects caused by oxidative stress due to its ability to increase the activity of antioxidant enzymes. However, there are no reports of the effects of high endogenous accumulation of proline in the transcriptional pattern of antioxidant enzymes genes under normal conditions of water supply or in response to water deficit. Here, we show that isoforms of four antioxidant enzymes genes (Ascorbate peroxidase-APX, Catalase-CAT, Superoxide dismutase-SOD and Glutathione reductase-GR) were differentially regulated in leaves of Swingle citrumelo transgenic plants with high endogenous proline accumulation submitted to water deficits and also under normal water supply condition. Proline per se caused a two-fold change in the transcription activity of APX1, APXcl, CAT2 and Cu/ZnSOD2, while during water deficit proline influenced mRNAs levels in APXs and Cu/ZnSODs isoforms, MnSODmit and GRcl. This study adds new information on the role of proline during drought conditions and, more important, without the potential confounding effects imposed by water deficiency. We showed that, in addition to its known effects on diverse plant physiological and biochemical processes, high endogenous proline can also acts as a regulatory/signalling molecule capable of altering the transcript levels of stress-related genes.

Gene expression and enzymatic activity of pectin methylesterase during fruit development and ripening in Coffea arabica L

Coffee quality is directly related to the harvest and post harvest conditions. Non-uniform maturation of coffee fruits, combined with inadequate harvest, negatively affects the final quality of the product. Pectin methylesterase (PME) plays an important role in fruit softening due to the hydrolysis of methylester groups in cell wall pectins. In order to characterize the changes occurring during coffee fruit maturation, the enzymatic activity of PME was measured during different stages of fruit ripening. PME activity progressively increased from the beginning of the ripening process to the cherry fruit stage. In silico analysis of expressed sequence tags of the Brazilian Coffee Genome Project database identified 5 isoforms of PME. We isolated and cloned a cDNA homolog of PME for further characterization. CaPME4 transcription was analyzed in pericarp, perisperm, and endosperm tissues during fruit development and ripening as well as in other plant tissues. Northern blot analysis revealed increased transcription of CaPME4 in the pericarp 300 days after flowering. Low levels of CaPME4 mRNAs were observed in the endosperm 270 days after flowering. Expression of CaPME4 transcripts was strong in the branches and lower in root and flower tissues. We showed that CaPME4 acts specifically during the later stages of fruit ripening and possibly contributes to the softening of coffee fruit, thus playing a significant role in pectin degradation in the fruit pericarp.

New mariner elements in Anastrepha species (Diptera: Tephritidae)

Mariner-like elements (MLE) are members from class II of transposable elements also known as DNA transposons. These elements have a wide distribution among different groups of organisms, including insects, which can be explained by horizontal and vertical gene-transfer. MLE families have been described in tephritid flies and other genera. During screening for Wolbachia bacteria in fruit flies of the genus Anastrepha, we discovered two sequences related to mariner-like elements. Based on these sequences, we designed primers that allowed us to isolate and characterize two new mariner-like elements (Anmar1 and Anmar2) in Anastrepha flies. These elements, which belong to the mellifera and rosa subfamilies have a low nucleotide diversity, and are probably inactive and acquired by vertical transfer. This is the first report of mariner-like transposons in flies found in South America.

Nitric oxide mediates the hormonal control of Crassulacean acid metabolism expression in young pineapple plants

Genotypic, developmental, and environmental factors converge to determine the degree of Crassulacean acid metabolism (CAM) expression. To characterize the signaling events controlling CAM expression in young pineapple (Ananas comosus) plants, this photosynthetic pathway was modulated through manipulations in water availability. Rapid, intense, and completely reversible up-regulation in CAM expression was triggered by water deficit, as indicated by the rise in nocturnal malate accumulation and in the expression and activity of important CAM enzymes. During both up- and down-regulation of CAM, the degree of CAM expression was positively and negatively correlated with the endogenous levels of abscisic acid (ABA) and cytokinins, respectively. When exogenously applied, ABA stimulated and cytokinins repressed the expression of CAM. However, inhibition of water deficit-induced ABA accumulation did not block the up-regulation of CAM, suggesting that a parallel, non-ABA-dependent signaling route was also operating. Moreover, strong evidence revealed that nitric oxide (NO) may fulfill an important role during CAM signaling. Up-regulation of CAM was clearly observed in NO-treated plants, and a conspicuous temporal and spatial correlation was also evident between NO production and CAM expression. Removal of NO from the tissues either by adding NO scavenger or by inhibiting NO production significantly impaired ABA-induced up-regulation of CAM, indicating that NO likely acts as a key downstream component in the ABA-dependent signaling pathway. Finally, tungstate or glutamine inhibition of the NO-generating enzyme nitrate reductase completely blocked NO production during ABA-induced up-regulation of CAM, characterizing this enzyme as responsible for NO synthesis during CAM signaling in pineapple plants.

Thermoperiod affects the diurnal cycle of nitrate reductase expression and activity in pineapple plants by modulating the endogenous levels of cytokinins

Nitrate reductase (NR, EC 1.6.6.1) activity in higher plants is regulated by a variety of environmental factors and oscillates with a characteristic diurnal rhythm. In this study, we have demonstrated that the diurnal cycle of NR expression and activity in pineapple (Ananas comosus, cv. Smooth Cayenne) can be strongly modified by changes in the day/night temperature regime. Plants grown under constant temperature (28 degrees C light/dark) showed a marked increase in the shoot NR activity (NRA) during the first half of the light period, whereas under thermoperiodic conditions (28 degrees C light/15 degrees C dark) significant elevations in the NRA were detected only in the root tissues at night. Under both conditions, increases in NR transcript levels occurred synchronically about 4 h prior to the corresponding elevation of the NRA. Diurnal analysis of endogenous cytokinins indicated that transitory increases in the levels of zeatin, zeatin riboside and isopentenyladenine riboside coincided with the accumulation of NR transcripts and preceded the rise of NRA in the shoot during the day and in the root at night, suggesting these hormones as mediators of the temperature-induced modifications of the NR cycle. Moreover, these cytokinins also induced NRA in pineapple when applied exogenously. Altogether, these results provide evidence that thermoperiodism can modify the diurnal cycle of NR expression and activity in pineapple both temporally and spatially, possibly by modulating the day/night changes in the cytokinin levels. A potential relationship between the day/night NR cycle and the photosynthetic pathway performed by the pineapple plants (C(3) or CAM) is also discussed.