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.

Genetic Diversity and Selection Footprints in the Genome of Brazilian Soybean Cultivars.

Although Brazil is currently the largest soybean producer in the world, only a small number of studies have analyzed the genetic diversity of Brazilian soybean. These studies have shown the existence of a narrow genetic base. The objectives of this work were to analyze the population structure and genetic diversity, and to identify selection signatures in the genome of soybean germplasms from different companies in Brazil. A panel consisting of 343 soybean lines from Brazil, North America, and Asia was genotyped using genotyping by sequencing (GBS). Population structure was assessed by Bayesian and multivariate approaches. Genetic diversity was analyzed using metrics such as the fixation index, nucleotide diversity, genetic dissimilarity, and linkage disequilibrium. The software BayeScan was used to detect selection signatures between Brazilian and Asian accessions as well as among Brazilian germplasms. Region of origin, company of origin, and relative maturity group (RMG) all had a significant influence on population structure. Varieties belonging to the same company and especially to the same RMG exhibited a high level of genetic similarity. This result was exacerbated among early maturing accessions. Brazilian soybean showed significantly lower genetic diversity when compared to Asian accessions. This was expected, because the crop's region of origin is its main genetic diversity reserve. We identified 7 genomic regions under selection between the Brazilian and Asian accessions, and 27 among Brazilian varieties developed by different companies. Associated with these genomic regions, we found 96 quantitative trait loci (QTLs) for important soybean breeding traits such as flowering, maturity, plant architecture, productivity components, pathogen resistance, and seed composition. Some of the QTLs associated with the markers under selection have genes of great importance to soybean's regional adaptation. The results reported herein allowed to expand the knowledge about the organization of the genetic variability of the Brazilian soybean germplasm. Furthermore, it was possible to identify genomic regions under selection possibly associated with the adaptation of soybean to Brazilian environments.

An 82 bp tandem repeat family typical of 3' non-coding end of Gypsy/TAT LTR retrotransposons is conserved in Coffea spp. pericentromeres.

Coffea spp. chromosomes are very small and accumulate a variety of repetitive DNA families around the centromeres. However, the proximal regions of Coffea chromosomes remain poorly understood, especially regarding the nature and organisation of the sequences. Taking advantage of the genome sequences of C. arabica (2n = 44), C. canephora, and C. eugenioides (C. arabica progenitors with 2n = 22) and good coverage genome sequencing of dozens of other wild Coffea spp., repetitive DNA sequences were identified, and the genomes were compared to decipher particularities of pericentromeric structures. The searches revealed a short tandem repeat (82 bp length) typical of Gypsy/TAT LTR retrotransposons, named Coffea_sat11. This repeat organises clusters with fragments of other transposable elements, comprising regions of non-coding RNA production. Cytogenomic analyses showed that Coffea_sat11 extends from the pericentromeres towards the middle of the chromosomal arms. This arrangement was observed in the allotetraploid C. arabica chromosomes, as well as in its progenitors. This study improves our understanding of the role of the Gypsy/TAT LTR retrotransposon lineage in the organisation of Coffea pericentromeres, as well as the conservation of Coffea_sat11 within the genus. The relationships between fragments of other transposable elements and the functional aspects of these sequences on the pericentromere chromatin were also evaluated. Highlights: A scattered short tandem repeat, typical of Gypsy/TAT LTR retrotransposons, associated with several fragments of other transposable elements, accumulates in the pericentromeres of Coffea chromosomes. This arrangement is preserved in all clades of the genus and appears to have a strong regulatory role in the organisation of chromatin around centromeres.

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.

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.

Population structure and genetic relationships between Ethiopian and Brazilian Coffea arabica genotypes revealed by SSR markers.

Information about population structure and genetic relationships within and among wild and brazilian Coffea arabica L. genotypes is highly relevant to optimize the use of genetic resources for breeding purposes. In this study, we evaluated genetic diversity, clustering analysis based on Jaccard's coefficient and population structure in 33 genotypes of C. arabica and of three diploid Coffea species (C. canephora, C. eugenioides and C. racemosa) using 30 SSR markers. A total of 206 alleles were identified, with a mean of 6.9 over all loci. The set of SSR markers was able to discriminate all genotypes and revealed that Ethiopian accessions presented higher genetic diversity than commercial varieties. Population structure analysis indicated two genetic groups, one corresponding to Ethiopian accessions and another corresponding predominantly to commercial cultivars. Thirty-four private alleles were detected in the group of accessions collected from West side of Great Rift Valley. We observed a lower average genetic distance of the C. arabica genotypes in relation to C. eugenioides than C. canephora. Interestingly, commercial cultivars were genetically closer to C. eugenioides than C. canephora and C. racemosa. The great allelic richness observed in Ethiopian Arabica coffee, especially in Western group showed that these accessions can be potential source of new alleles to be explored by coffee breeding programs.

Influence of Silver Nitrate on Somatic Embryogenesis Induction in Arabica Coffee (Coffea arabica L)

Abstract The influence of silver nitrate (AgNO3), benzyladenine (BAP), and indole-3-acetic acid (IAA) on low frequency somatic embryogenesis (LFSE) induction in Caturra and Catuaí arabica coffee was evaluated. For the Caturra cultivar, the production of somatic embryos was significantly increased by adding AgNO3 to the semisolid culture medium. The highest average number of somatic embryos for this cultivar was obtained using 6.6 μM BAP, 2.85 μM IAA, and 40 μM AgNO3. In contrast, for the Catuaí cultivar, the highest average number of somatic embryos was obtained using semisolid medium supplemented with 8.8 μM BAP, and 2.85 μM IAA. Using these protocols, somatic embryos were directly induced using leaf sections of in vitro plants of both coffee cultivars within 8 weeks. The somatic embryos developed into rooted plants with a 100% survival rate upon transfer to the greenhouse.

BASiNET-BiologicAl Sequences NETwork: a case study on coding and non-coding RNAs identification.

With the emergence of Next Generation Sequencing (NGS) technologies, a large volume of sequence data in particular de novo sequencing was rapidly produced at relatively low costs. In this context, computational tools are increasingly important to assist in the identification of relevant information to understand the functioning of organisms. This work introduces BASiNET, an alignment-free tool for classifying biological sequences based on the feature extraction from complex network measurements. The method initially transform the sequences and represents them as complex networks. Then it extracts topological measures and constructs a feature vector that is used to classify the sequences. The method was evaluated in the classification of coding and non-coding RNAs of 13 species and compared to the CNCI, PLEK and CPC2 methods. BASiNET outperformed all compared methods in all adopted organisms and datasets. BASiNET have classified sequences in all organisms with high accuracy and low standard deviation, showing that the method is robust and non-biased by the organism. The proposed methodology is implemented in open source in R language and freely available for download at https://cran.r-project.org/package=BASiNET.

Transcriptome Analysis of Leaves, Flowers and Fruits Perisperm of Coffea arabica L. Reveals the Differential Expression of Genes Involved in Raffinose Biosynthesis.

Coffea arabica L. is an important crop in several developing countries. Despite its economic importance, minimal transcriptome data are available for fruit tissues, especially during fruit development where several compounds related to coffee quality are produced. To understand the molecular aspects related to coffee fruit and grain development, we report a large-scale transcriptome analysis of leaf, flower and perisperm fruit tissue development. Illumina sequencing yielded 41,881,572 high-quality filtered reads. De novo assembly generated 65,364 unigenes with an average length of 1,264 bp. A total of 24,548 unigenes were annotated as protein coding genes, including 12,560 full-length sequences. In the annotation process, we identified nine candidate genes related to the biosynthesis of raffinose family oligossacarides (RFOs). These sugars confer osmoprotection and are accumulated during initial fruit development. Four genes from this pathway had their transcriptional pattern validated by quantitative reverse transcription polymerase chain reaction (qRT-PCR). Furthermore, we identified ~24,000 putative target sites for microRNAs (miRNAs) and 134 putative transcriptionally active transposable elements (TE) sequences in our dataset. This C. arabica transcriptomic atlas provides an important step for identifying candidate genes related to several coffee metabolic pathways, especially those related to fruit chemical composition and therefore beverage quality. Our results are the starting point for enhancing our knowledge about the coffee genes that are transcribed during the flowering and initial fruit development stages.

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.

Transcriptome analysis in Coffea eugenioides, an Arabica coffee ancestor, reveals differentially expressed genes in leaves and fruits.

Studies in diploid parental species of polyploid plants are important to understand their contributions to the formation of plant and species evolution. Coffea eugenioides is a diploid species that is considered to be an ancestor of allopolyploid Coffea arabica together with Coffea canephora. Despite its importance in the evolutionary history of the main economic species of coffee, no study has focused on C. eugenioides molecular genetics. RNA-seq creates the possibility to generate reference transcriptomes and identify coding genes and potential candidates related to important agronomic traits. Therefore, the main objectives were to obtain a global overview of transcriptionally active genes in this species using next-generation sequencing and to analyze specific genes that were highly expressed in leaves and fruits with potential exploratory characteristics for breeding and understanding the evolutionary biology of coffee. A de novo assembly generated 36,935 contigs that were annotated using eight databases. We observed a total of ~5000 differentially expressed genes between leaves and fruits. Several genes exclusively expressed in fruits did not exhibit similarities with sequences in any database. We selected ten differentially expressed unigenes in leaves and fruits to evaluate transcriptional profiles using qPCR. Our study provides the first gene catalog for C. eugenioides and enhances the knowledge concerning the mechanisms involved in the C. arabica homeologous. Furthermore, this work will open new avenues for studies into specific genes and pathways in this species, especially related to fruit, and our data have potential value in assisted breeding applications.

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.

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.

FISH using a gag-like fragment probe reveals a common Ty3-gypsy-like retrotransposon in genome of Coffea species.

The genus Coffea possesses about 100 species, and the most economically important are Coffea canephora and Coffea arabica. The latter is predominantly self-compatible with 2n = 4x = 44, while the others of the genus are diploid with 2n = 2x = 22 and mostly self-incompatible. Studies using molecular markers have been useful to detect differences between genomes in Coffea; however, molecular and cytogenetic studies have produced only limited information on the karyotypes organization. We used DOP-PCR to isolate repetitive elements from genome of Coffea arabica var. typica. The pCa06 clone, containing a fragment of 775 bp length, was characterized by sequencing and used as a probe in chromosomes of C. arabica and six other species: C. canephora, Coffea eugenioides, Coffea kapakata, Coffea liberica var. dewevrei, Coffea racemosa, and Coffea stenophylla. This insert shows similarities with a gag protein of the Ty3-gypsy-like super-family. Dot blot and FISH analyses demonstrated that pCa06 is differentially accumulated between species and chromosomes. Signals appeared scattered and clustered on the chromosomes and were also associated with heterochromatic regions. While the literature shows that there is a high karyotype similarity between Coffea species, our results point out differences in the accumulation and dispersion of this Ty3-gypsy-like retrotransposon during karyotype differentiation of Coffea.

Characterization and expression of two cDNA encoding 3-Hydroxy-3-methylglutaryl coenzyme A reductase isoforms in coffee (Coffea arabica L.).

In higher plants there are two independent pathways for isoprenoid biosynthesis, located in the cytosol (mevalonic acid or MVA pathway) or in the plastids [methylerythritol phosphate (MEP) pathway]. The 3-hydroxy-3-methyglutaryl-CoA reductase (HMGR) is the first committed step in the MVA pathway. Using the information available from the Brazilian Coffee Genome Project, we found 13 ESTs that originated two isoforms, CaHMGR1 and CaHMGR2, for the enzyme HMGR of Coffea arabica. A complementary DNA encoding the isoform CaHMGR1 was cloned, and its complete nucleotide sequence determined. The full-length cDNA of CaHMGR1 was 2,242 bp containing a 1,812-bp ORF encoding 604 amino acids. Bioinformatic analyses revealed that the deduced CaHMGR1 had extensive homology with other plant HMGRs and contained two transmembrane domains and two putative HMGR binding sites and two NADP(H)-binding sites. Under normal growth conditions, transcripts of isoform CaHMRG1 were detected in fruit tissues (pulp, perisperm, and endosperm) only at the initial stages of development, flower buds and leaves. CaHMRG2 was expressed in all tissues and during all fruit development stages examined. These results suggest a constitutive expression of isoform CaHMGR2, while the isoform CaHMGR1 shows temporal and tissue-specific transcriptional activation.

A high-throughput data mining of single nucleotide polymorphisms in Coffea species expressed sequence tags suggests differential homeologous gene expression in the allotetraploid Coffea arabica.

Polyploidization constitutes a common mode of evolution in flowering plants. This event provides the raw material for the divergence of function in homeologous genes, leading to phenotypic novelty that can contribute to the success of polyploids in nature or their selection for use in agriculture. Mounting evidence underlined the existence of homeologous expression biases in polyploid genomes; however, strategies to analyze such transcriptome regulation remained scarce. Important factors regarding homeologous expression biases remain to be explored, such as whether this phenomenon influences specific genes, how paralogs are affected by genome doubling, and what is the importance of the variability of homeologous expression bias to genotype differences. This study reports the expressed sequence tag assembly of the allopolyploid Coffea arabica and one of its direct ancestors, Coffea canephora. The assembly was used for the discovery of single nucleotide polymorphisms through the identification of high-quality discrepancies in overlapped expressed sequence tags and for gene expression information indirectly estimated by the transcript redundancy. Sequence diversity profiles were evaluated within C. arabica (Ca) and C. canephora (Cc) and used to deduce the transcript contribution of the Coffea eugenioides (Ce) ancestor. The assignment of the C. arabica haplotypes to the C. canephora (CaCc) or C. eugenioides (CaCe) ancestral genomes allowed us to analyze gene expression contributions of each subgenome in C. arabica. In silico data were validated by the quantitative polymerase chain reaction and allele-specific combination TaqMAMA-based method. The presence of differential expression of C. arabica homeologous genes and its implications in coffee gene expression, ontology, and physiology are discussed.

Biochemical and genomic analysis of sucrose metabolism during coffee (Coffea arabica) fruit development.

Sucrose metabolism and the role of sucrose synthase were investigated in the fruit tissues (pericarp, perisperm, and endosperm) of Coffea arabica during development. Acid invertase, sucrose phosphate synthase, and sucrose synthase activities were monitored and compared with the levels of sucrose and reducing sugars. Among these enzymes, sucrose synthase showed the highest activities during the last stage of endosperm and pericarp development and this activity paralleled closely the accumulation of sucrose in these tissues at this stage. Carbon partitioning in fruits was studied by pulse-chase experiments with (14)C-sugars and revealed high rates of sucrose turnover in perisperm and endosperm tissues. Additional feeding experiments with (14)CO(2) showed that leaf photosynthesis contributed more to seed development than the pericarp in terms of photosynthate supply to the endosperm. Sugar analysis, feeding experiments, and histological studies indicated that the perisperm plays an important role in this downloading process. It was observed that the perisperm presents a transient accumulation of starch which is degraded as the seed develops. Two full-length cDNAs (CaSUS1 and CaSUS2) and the complete gene sequence of the latter were also isolated. They encode sucrose synthase isoforms that are phylogenetically distinct, indicating their involvement in different physiological functions during cherry development. Contrasting expression patterns were observed for CaSUS1 and CaSUS2 in perisperm, endosperm, and pericarp tissues: CaSUS1 mRNAs accumulated mainly during the early development of perisperm and endosperm, as well as during pericarp growing phases, whereas those of CaSUS2 paralleled sucrose synthase activity in the last weeks of pericarp and endosperm development. Taken together, these results indicate that sucrose synthase plays an important role in sugar metabolism during sucrose accumulation in the coffee fruit.

Production of herbicide-resistant coffee plants ( Coffea canephora P. ) via Agrobacterium tumefaciens-mediated transformation

Plantas transgênicas de Coffea canephora P resistentes ao herbicida glufosinato de amônio foram regeneradas a partir de explantes foliares co-cultivados com Agrobacterium tumefaciens EHA105 contendo o plasmídio pCambia3301 que contém os genes bar e uidA ambos sob controle do promotor 35S. Embriogênese somática direta foi induzida no meio contendo » da concentração de macro, metade da concentração de micronutrientes do meio MS, constituintes orgânicos do meio B5 e 30 g.L-1 de sacarose suplementado com 5µM N6 – (2-isopentenil)-adenina (2-iP) e 10 µM de glufosinato de amônio para seleção de embriões transgênicos putativos. A presença e a integração do gene bar foram confirmados pelas análises de PCR e Southern blot. As plantas transgênicas selecionadas de café, pulverizadas com 1600 mg.L-1 do herbicida FinaleÔ que contém glufosinato como ingrediente ativo, mantiveram a coloração e continuaram crescendo normalmente na aclimatação ex vitro.