Haitian coffee agroforestry systems harbor complex arabica variety mixtures and under-recognized genetic diversity

Though facing significant challenges, coffee (Coffea arabica) grown in Haitian agroforestry systems are important contributors to rural livelihoods and provide several ecosystem services. However, little is known about their genetic diversity and the variety mixtures used. In light of this, there is a need to characterize Haitian coffee diversity to help inform revitalization of this sector. We sampled 28 diverse farms in historically important coffee growing regions of northern and southern Haiti. We performed KASP-genotyping of SNP markers and HiPlex multiplex amplicon sequencing for haplotype calling on our samples, as well as several Ethiopian and commercial accessions from international collections. This allowed us to assign Haitian samples to varietal groups. Our analyses revealed considerable genetic diversity in Haitian farms, higher in fact than many farmers realized. Notably, genetic structure analyses revealed the presence of clusters related to Typica, Bourbon, and Catimor groups, another group that was not represented in our reference accession panel, and several admixed individuals. Across the study areas, we found both mixed-variety farms and monovarietal farms with the historical and traditional Typica variety. This study is, to our knowledge, the first to genetically characterize Haitian C. arabica variety mixtures, and report the limited cultivation of C. canephora (Robusta coffee) in the study area. Our results show that some coffee farms are repositories of historical, widely-abandoned varieties while others are generators of new diversity through genetic mixing.

Combined sensory, volatilome and transcriptome analyses identify a limonene terpene synthase as a major contributor to the characteristic aroma of a Coffea arabica L. specialty coffee

BACKGROUND

The fruity aromatic bouquet of coffee has attracted recent interest to differentiate high value market produce as specialty coffee. Although the volatile compounds present in green and roasted coffee beans have been extensively described, no study has yet linked varietal molecular differences to the greater abundance of specific substances and support the aroma specificity of specialty coffees.

RESULTS

This study compared four Arabica genotypes including one, Geisha Especial, suggested to generate specialty coffee. Formal sensory evaluations of coffee beverages stressed the importance of coffee genotype in aroma perception and that Geisha Especial-made coffee stood out by having fine fruity, and floral, aromas and a more balanced acidity. Comparative SPME-GC-MS analyses of green and roasted bean volatile compounds indicated that those of Geisha Especial differed by having greater amounts of limonene and 3-methylbutanoic acid in agreement with the coffee cup aroma perception. A search for gene ontology differences of ripening beans transcriptomes of the four varieties revealed that they differed by metabolic processes linked to terpene biosynthesis due to the greater gene expression of prenyl-pyrophosphate biosynthetic genes and terpene synthases. Only one terpene synthase (CaTPS10-like) had an expression pattern that paralleled limonene loss during the final stage of berry ripening and limonene content in the studied four varieties beans. Its functional expression in tobacco leaves confirmed its functioning as a limonene synthase.

CONCLUSIONS

Taken together, these data indicate that coffee variety genotypic specificities may influence ripe berry chemotype and final coffee aroma unicity. For the specialty coffee variety Geisha Especial, greater expression of terpene biosynthetic genes including CaTPS10-like, a limonene synthase, resulted in the greater abundance of limonene in green beans, roasted beans and a unique citrus note of the coffee drink.

Genetic diversity and structure of the coffee leaf rust fungus Hemileia vastatrix across different coffee management systems in Ethiopia

Although coffee leaf rust (CLR), caused by Hemileia vastatrix, poses an increasing threat to coffee production in Ethiopia, little is known regarding its genetic diversity and structure and how these are affected by coffee management. Here, we used genetic fingerprinting based on sequence-related amplified polymorphism (SRAP) markers to genotype H. vastatrix samples from different coffee shrubs, across 40 sites, covering four coffee production systems (forest coffee, semi plantation coffee, home garden coffee, and plantation coffee) and different altitudes in Ethiopia. In total, 96 H. vastatrix samples were successfully genotyped with three primer combinations, producing a total of 79 scorable bands. We found 35.44% of amplified bands to be polymorphic, and the polymorphic information content (PIC) was 0.45, suggesting high genetic diversity among our CLR isolates. We also found significant isolation-by-distance across the samples investigated and detected significant differences in fungal genetic composition among plantation coffee and home garden coffee and a marginally significant difference among plantation coffee and forest coffee. Furthermore, we found a significant effect of altitude on CLR genetic composition in the forest coffee and plantation systems. Our results suggest that both spore dispersal and different selection pressures in the different coffee management systems are likely responsible for the observed high genetic diversity and genetic structure of CLR isolates in Ethiopia. When selecting Ethiopian coffee genotypes for crop improvement, it is important that these genotypes carry some resistance against CLR. Because our study shows large variation in genetic composition across relatively short geographical distances, a broad selection of rust isolates must be used for coffee resistance screening.

The genome and population genomics of allopolyploid Coffea arabica reveal the diversification history of modern coffee cultivars

Coffea arabica, an allotetraploid hybrid of Coffea eugenioides and Coffea canephora, is the source of approximately 60% of coffee products worldwide, and its cultivated accessions have undergone several population bottlenecks. We present chromosome-level assemblies of a di-haploid C. arabica accession and modern representatives of its diploid progenitors, C. eugenioides and C. canephora. The three species exhibit largely conserved genome structures between diploid parents and descendant subgenomes, with no obvious global subgenome dominance. We find evidence for a founding polyploidy event 350,000-610,000 years ago, followed by several pre-domestication bottlenecks, resulting in narrow genetic variation. A split between wild accessions and cultivar progenitors occurred ~30.5 thousand years ago, followed by a period of migration between the two populations. Analysis of modern varieties, including lines historically introgressed with C. canephora, highlights their breeding histories and loci that may contribute to pathogen resistance, laying the groundwork for future genomics-based breeding of C. arabica.

Evaluation of suitable reference genes for expression profile analyses of target genes in the coffee berry borer, Hypothenemus hampei (Ferrari) (Coleoptera: Curculionidae)

The coffee berry borer, Hypothenemus hampei (Ferrari) (Coleoptera: Curculionidae), is a major destructive insect pest of coffee, which impacts the coffee crops negatively. As a draft genome has been completed for this insect, most molecular studies on gene transcriptional levels under different experimental conditions will be conducted using real-time reverse-transcription quantitative polymerase chain reactions (RT-qPCR). However, the lack of suitable internal reference genes will affect the accuracy of RT-qPCR results. In this study, the expression stability of nine candidate reference genes was evaluated under different developmental stages, temperature stress, and Beauveria bassiana infection. Data analyses were completed by four commonly used programs, BestKeeper, NormFinder, geNorm, and RefFinder. The result showed that RPL3 and EF1α combination were recommended as the most stable reference genes for developmental stages. EF1α and RPS3a combination were the top two stable reference genes for B. bassiana infection. RPS3a and RPL3 combination performed as the optimal reference genes both in temperature stress and all samples. Our results should provide a good foundation for the expression profile analyses of target genes in the future, especially for molecular studies on insect genetic development, temperature adaptability, and immune mechanism to entomogenous fungi in H. hampei.

The evolutionary history of three Baracoffea species from western Madagascar revealed by chloroplast and nuclear genomes

The wild species of the Coffea genus present a very wide morphological, genetic, and biochemical diversity. Wild species are recognized more resistant to diseases, pests, and environmental variations than the two species currently cultivated worldwide: C. arabica (Arabica) and C. canephora (Robusta). Consequently, wild species are now considered as a crucial resource for adapting cultivated coffee trees to climate change. Within the Coffea genus, 79 wild species are native to the Indian Ocean islands of Comoros, Mayotte, Mauritius, Réunion and Madagascar, out of a total of 141 taxa worldwide. Among them, a group of 9 species called "Baracoffea" are particularly atypical in their morphology and adaptation to the sandy soils of the dry deciduous forests of western Madagascar. Here, we have attempted to shed light on the evolutionary history of three Baracoffea species: C. ambongensis, C. boinensis and C. bissetiae by analyzing their chloroplast and nuclear genomes. We assembled the complete chloroplast genomes de novo and extracted 28,800 SNP (Single Nucleotide Polymorphism) markers from the nuclear genomes. These data were used for phylogenetic analysis of Baracoffea with Coffea species from Madagascar and Africa. Our new data support the monophyletic origin of Baracoffea within the Coffea of Madagascar, but also reveal a divergence with a sister clade of four species: C. augagneurii, C. ratsimamangae, C. pervilleana and C. Mcphersonii (also called C. vohemarensis), belonging to the Subterminal botanical series and living in dry or humid forests of northern Madagascar. Based on a bioclimatic analysis, our work suggests that Baracoffea may have diverged from a group of Malagasy Coffea from northern Madagascar and adapted to the specific dry climate and low rainfall of western Madagascar. The genomic data generated in the course of this work will contribute to the understanding of the adaptation mechanisms of these particularly singular species.

Coffee Cell Suspensions as a Platform for Transient Gene Expression Analysis

Coffea arabica L. is a crucial crop globally, but its genetic homogeneity leads to its susceptibility to diseases and pests like the coffee berry borer (CBB). Chemical and cultural control methods are difficult due to the majority of the CBB life cycle taking place inside coffee beans. One potential solution is the use of the gene cyt1Aa from Bacillus thuringiensis as a biological insecticide. To validate candidate genes against CBB, a simple, rapid, and efficient transient expression system is necessary. This study uses cell suspensions as a platform for expressing the cyt1Aa gene in the coffee genome (C. arabica L. var. Catuaí) to control CBB. The Agrobacterium tumefaciens strain GV3101::pMP90 containing the bar and cyt1Aa genes are used to genetically transform embryogenic cell suspensions. PCR amplification of the cyt1Aa gene is observed 2, 5, and 7 weeks after infection. This chapter describes a protocol that can be used for the development of resistant varieties against biotic and abiotic stresses and CRISPR/Cas9-mediated genome editing.

Transcriptomic Data Analysis Using the Galaxy Platform: Coffee (Coffea arabica L.) Flowers as Example

Coffee, an important agricultural product for tropical producing countries, is facing challenges due to climate change, including periods of drought, irregular rain distribution, and high temperatures. These changes result in plant water stress, leading to significant losses in coffee productivity and quality. Understanding the processes that affect coffee flowering is crucial for improving productivity and quality. In this chapter, we describe a protocol for transcriptome analysis using available Internet software, mainly in the Galaxy Platform, using RNA-Seq data from flowers collected from different parts of the coffee tree. The methods presented in this chapter provide a comprehensive protocol for transcriptome analysis of differentially expressed genes from flowers of coffee plant. This knowledge can be utilized in coffee genetic improvement programs, particularly in the selection of cultivars that are tolerant to water deficit.

Agroinfiltration for Enhanced Transgene Expression in Coffee Leaves (Coffea arabica L.)

The Coffea spp. plant is a significant crop in Latin America, Africa, and Asia, and recent advances in genomics and transcriptomics have opened possibilities for studying candidate genes and introducing new desirable traits through genetic engineering. While stable transformation of coffee plants has been reported using various techniques, it is a time-consuming and laborious process. To overcome this, transient transformation methods have been developed, which avoid the limitations of stable transformation. This chapter describes an ex vitro protocol for transient expression using A. tumefaciens-mediated infiltration of coffee leaves, which could be used to produce coffee plants expressing desirable traits against biotic and abiotic stresses, genes controlling biochemical and physiological traits, as well as for gene editing through CRISPR/Cas9.

A General Protocol for Gene Expression Analysis in Chemically Mutant Coffee Plants in Response to Rust (Hemileia vastatrix Berk & Broome) Infection

Coffea spp. is the source of one of the most widely consumed beverages in the world. However, the cultivation of this crop is threatened by Hemileia vastatrix Berk & Broome, a fungal disease, which reduces the productivity and can cause significant economic losses. In this protocol, coffee leaf segment derived from a chemical mutagenesis process are inoculated with uredospores of the pathogen. Subsequently, the gene expression changes are analyzed over the time (0, 5, 24, 48, and 120 h) using quantitative real-time polymerase chain reaction (RT-qPCR). The procedures and example data are presented for expression analysis in the CaWRKY1 gene. This procedure can be applied for quantitative analysis of other genes of interest to coffee breeders and scientists for elucidating the molecular mechanisms involved in the interaction between the plant and pathogen, potentially leading to the development of more efficient approaches for managing this disease.

The terroir of Brazilian Coffea canephora: Characterization of the chemical composition

FTIR spectroscopy and multivariate analysis were used in the chemical study of the terroirs of Coffea canephora. Conilon coffees from Espírito Santo and Amazon robusta from Matas of Rondônia, were separated by PCA, with lipids and caffeine being the markers responsible for the separation. Coffees from Bahia, Minas Gerais, and São Paulo did not exhibit separation, indicating that the botanical variety had a greater effect on the terroir than geographic origin. Thus, the genetic factor was investigated considering the conilon and robusta botanical varieties. This last group was composed of hybrid robusta and apoatã. The DD-SIMCA favored the identification of the genetic predominance of the samples. PLS-DA had a high classification performance regarding the conilon, hybrid robusta, and apoatã genetic nature. Lipids, caffeine, chlorogenic acids, quinic acid, trigonelline, proteins, amino acids, and carbohydrates were identified as chemical markers that discriminated the genetic groups.

High-Density Chitosan Induces a Biochemical and Molecular Response in Coffea arabica during Infection with Hemileia vastatrix

The coffee industry faces coffee leaf rust caused by Hemileia vastratix, which is considered the most devastating disease of the crop, as it reduces the photosynthetic rate and limits productivity. The use of plant resistance inducers, such as chitosan, is an alternative for the control of the disease by inducing the synthesis of phytoalexins, as well as the activation of resistance genes. Previously, the effect of chitosan from different sources and physicochemical properties was studied; however, its mechanisms of action have not been fully elucidated. In this work, the ability of food-grade high-density chitosan (0.01% and 0.05%) to control the infection caused by the pathogen was evaluated. Subsequently, the effect of high-density chitosan (0.05%) on the induction of pathogenesis-related gene expression (GLUC, POX, PAL, NPR1, and CAT), the enzymatic activity of pathogenesis-related proteins (GLUC, POX, SOD, PPO, and APX), and phytoalexin production were evaluated. The results showed that 0.05% chitosan increased the activity and gene expression of ß-1,3 glucanases and induced a differentiated response in enzymes related to the antioxidant system of plants. In addition, a correlation was observed between the activities of polyphenol oxidase and the production of phytoalexin, which allowed an effective defense response in coffee plants.

New genetic markers for 100% arabica coffee demonstrate high discriminatory potential for InDel-HRM-based coffee authentication

Food authenticity is crucial in today's society, given the heightened consumer awareness and attention to the products they consume. Reliable and efficient techniques are needed to quickly detect potential food adulterations that can negatively impact product quality and economic value. Coffee, a globally traded agricultural product, holds immense economic importance, with an estimated value of USD 83 billion. It is widely consumed and recognized as a functional food that provides minerals (K, Mg, Mn, Cr), niacin, and antioxidants. However, the preferred coffee species, Coffea arabica, known for its superior drink quality, is often adulterated with Coffea canephora (Robusta and Conilon) beans, even in 100% Arabica coffee. To distinguish between these two coffee species, a comprehensive study was conducted using a robust approach to identify differences in Single-Ortholog Copy (SOC) based on InDel regions in these gene pairs. These differences were validated using a meticulous methodology that considered variations in amplicon size: electrophoretic profile, and high-resolution melting (HRM). The innovative combination of InDels and HRM resulted in highly distinctive HRM profiles, outperforming SNP-based methods previously used. The targeted InDel approach utilized in this study facilitated precise quantification of Coffea species beans with a detection sensitivity of 0.5%. The study's findings establish the reliability and accuracy in distinguishing between the two coffee species, showcasing the valuable application of InDels for quality control and ensuring the authenticity of coffee beans. This pioneering research contributes to the advancement of authenticity verification methods for both imported and exported coffee beans, as well as in future studies that require significant genetic differences between these species, such as C. arabica and C. canephora.

Potassium Phosphite Activates Components Associated with Constitutive Defense Responses in Coffea arabica Cultivars

Phosphites have been used as inducers of resistance, activating the defense of plants and increasing its ability to respond to the invasion of the pathogen. However, the mode of action of phosphites in defense responses has not yet been fully elucidated. The objective of this study was to evaluate the effect of potassium phosphite (KPhi) in coffee cultivars with different levels of resistance to rust to clarify the mechanism by which KPhi activates the constitutive defense of plants. To this end, we studied the expression of genes and the activity of enzymes involved in the defense pathway of salicylic acid (SA) and reactive oxygen species (ROS), in addition to the levels of total soluble phenolic compounds and soluble lignin. Treatment with KPhi induced constitutive defense responses in cultivars resistant and susceptible to rust. The results suggest that KPhi acts in two parallel defense pathways, SA and ROS, which are essential for the induction of systemic acquired resistance (SAR) when activated simultaneously. The activation of the mechanisms associated with defense routes demonstrates that KPhi is a potential inducer of resistance in coffee plants.

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.

The chloroplast protein HCF164 is predicted to be associated with Coffea SH9 resistance factor against Hemileia vastatrix

To explore the connection between chloroplast and coffee resistance factors, designated as SH1 to SH9, whole genomic DNA of 42 coffee genotypes was sequenced, and entire chloroplast genomes were de novo assembled. The chloroplast phylogenetic haplotype network clustered individuals per species instead of SH factors. However, for the first time, it allowed the molecular validation of Coffea arabica as the maternal parent of the spontaneous hybrid "Híbrido de Timor". Individual reads were also aligned on the C. arabica reference genome to relate SH factors with chloroplast metabolism, and an in-silico analysis of selected nuclear-encoded chloroplast proteins (132 proteins) was performed. The nuclear-encoded thioredoxin-like membrane protein HCF164 enabled the discrimination of individuals with and without the SH9 factor, due to specific DNA variants linked to chromosome 7c (from C. canephora-derived sub-genome). The absence of both the thioredoxin domain and redox-active disulphide center in the HCF164 protein, observed in SH9 individuals, raises the possibility of potential implications on redox regulation. For the first time, the identification of specific DNA variants of chloroplast proteins allows discriminating individuals according to the SH profile. This study introduces an unexplored strategy for identifying protein/genes associated with SH factors and candidate targets of H. vastatrix effectors, thereby creating new perspectives for coffee breeding programs.

Identification of sources of male sterility in the Colombian Coffee Collection for the genetic improvement of Coffea arabica L

In coffee (Coffea arabica L.), male sterility is a prerequisite for the exploitation of heterosis since it provides an efficient and reliable method for the production of hybrid seeds. Given its relevance, the objective of this study was to identify male-sterile genotypes within the Colombian Coffee Collection that can be used in genetic improvement. For this purpose, Ethiopian germplasm and progenies derived from hybrids between C. arabica x C. canephora were explored between 2017 and 2021. In the first stage, genotypes without visual presence of pollen were preselected in the field, followed by selection through staining and verification of male sterility and female fertility through directed crosses (directed, reciprocal and selfing). In this stage, 9,753 trees were explored, preselecting 2.4% due to visual absence of pollen. The staining of structures allowed us to confirm the lack or sporadic production of pollen in 23 individuals of Ethiopian origin. The results of the directed crosses led to the identification of 11 male-sterile and 12 partially male-sterile genotypes belonging to 15 accessions. In all cases, the individuals were characterized by the presence of anthers but with an absence or low content of pollen, which is why the male sterility is possibly of the sporogenic type. The female receptivity values were between 2.9% and 72.6%, being higher than 30% in five genotypes. These genotypes are a valuable tool for the genetic improvement of C. arabica with the potential to facilitate the use of heterosis and to allow a deeper understanding the development of male gametophytes in the species.

An analysis of codon utilization patterns in the chloroplast genomes of three species of Coffea

BACKGROUND

The chloroplast genome of plants is known for its small size and low mutation and recombination rates, making it a valuable tool in plant phylogeny, molecular evolution, and population genetics studies. Codon usage bias, an important evolutionary feature, provides insights into species evolution, gene function, and the expression of exogenous genes. Coffee, a key crop in the global tropical agricultural economy, trade, and daily life, warrants investigation into its codon usage bias to guide future research, including the selection of efficient heterologous expression systems for coffee genetic transformation.

RESULTS

Analysis of the codon utilization patterns in the chloroplast genomes of three Coffea species revealed a high degree of similarity among them. All three species exhibited similar base compositions, with high A/T content and low G/C content and a preference for A/T-ending codons. Among the 30 high-frequency codons identified, 96.67% had A/T endings. Fourteen codons were identified as ideal. Multiple mechanisms, including natural selection, were found to influence the codon usage patterns in the three coffee species, as indicated by ENc-GC3s mapping, PR2 analysis, and neutral analysis. Nicotiana tabacum and Saccharomyces cerevisiae have potential value as the heterologous expression host for three species of coffee genes.

CONCLUSION

This study highlights the remarkable similarity in codon usage patterns among the three coffee genomes, primarily driven by natural selection. Understanding the gene expression characteristics of coffee and elucidating the laws governing its genetic evolution are facilitated by investigating the codon preferences in these species. The findings can enhance the efficacy of exogenous gene expression and serve as a basis for future studies on coffee evolution.

Embryonic abnormalities and genotoxicity induced by 2,4-dichlorophenoxyacetic acid during indirect somatic embryogenesis in Coffea

Indirect somatic embryogenesis (ISE) is a morphogenetic pathway in which somatic cells form callus and, later, somatic embryos (SE). 2,4-dichlorophenoxyacetic acid (2,4-D) is a synthetic auxin that promotes the proliferation and dedifferentiation of somatic cells, inducing the ISE. However, 2,4-D can cause genetic, epigenetic, physiological and morphological disorders, preventing the regeneration and/or resulting abnormal somatic embryos (ASE). We aimed to evaluate the toxic 2,4-D effect during the Coffea arabica and C. canephora ISE, assessing the SE morphology, global 5-methylcytosine levels (5-mC%) and DNA damage. Leaf explants were inoculated in media with different 2,4-D concentrations. After 90 days, the friable calli were transferred to the regeneration medium, and the number of normal and abnormal SE was monthly counted. The increase of the 2,4-D concentration increased the number of responsive explants in both Coffea. At 9.06, 18.08 and 36.24 μM 2,4-D, C. arabica presented the highest values of responsive explants, differing from C. canephora. Normal and abnormal SE regeneration increased in relation to the time and 2,4-D concentration. Global 5-mC% varied at different stages of the ISE in both Coffea. Furthermore, the 2,4-D concentration positively correlated with global 5-mC%, and with the mean number of ASE. All ASE of C. arabica and C. canephora exhibited DNA damage and showed higher global 5-mC%. The allotetraploid C. arabica exhibited greater tolerance to the toxic effect of 2,4-D than the diploid C. canephora. We conclude that synthetic 2,4-D auxin promotes genotoxic and phytotoxic disorders and promotes epigenetic changes during Coffea ISE.

Genetic diversity and structure in wild Robusta coffee (Coffea canephora A. Froehner) populations in Yangambi (DR Congo) and their relation to forest disturbance

Degradation and regeneration of tropical forests can strongly affect gene flow in understorey species, resulting in genetic erosion and changes in genetic structure. Yet, these processes remain poorly studied in tropical Africa. Coffea canephora is an economically important species, found in the understorey of tropical rainforests of Central and West Africa, and the genetic diversity harboured in its wild populations is vital for sustainable coffee production worldwide. Here, we aimed to quantify genetic diversity, genetic structure, and pedigree relations in wild C. canephora populations, and we investigated associations between these descriptors and forest disturbance and regeneration. Therefore, we sampled 256 C. canephora individuals within 24 plots across three forest categories in Yangambi (DR Congo), and used genotyping-by-sequencing to identify 18,894 SNPs. Overall, we found high genetic diversity, and no evidence of genetic erosion in C. canephora in disturbed old-growth forest, as compared to undisturbed old-growth forest. In addition, an overall heterozygosity excess was found in all populations, which was expected for a self-incompatible species. Genetic structure was mainly a result of isolation-by-distance, reflecting geographical location, with low to moderate relatedness at finer scales. Populations in regrowth forest had lower allelic richness than populations in old-growth forest and were characterised by a lower inter-individual relatedness and a lack of isolation-by-distance, suggesting that they originated from different neighbouring populations and were subject to founder effects. Wild Robusta coffee populations in the study area still harbour high levels of genetic diversity, yet careful monitoring of their response to ongoing forest degradation remains required.

A bitter future for coffee production? Physiological traits associated with yield reveal high vulnerability to hydraulic failure in Coffea canephora

The increase in frequency and intensity of drought events have hampered coffee production in the already threatened Amazon region, yet little is known about key aspects underlying the variability in yield potential across genotypes, nor to what extent higher productivity is linked to reduced drought tolerance. Here we explored how variations in morphoanatomical and physiological leaf traits can explain differences in yield and vulnerability to embolism in 11 Coffea canephora genotypes cultivated in the Western Amazon. The remarkable variation in coffee yield across genotypes was tightly related to differences in their carbon assimilation and water transport capacities, revealing a diffusive limitation to photosynthesis linked by hydraulic constraints. Although a clear trade-off between water transport efficiency and safety was not detected, all the studied genotypes operated in a narrow and/or negative hydraulic safety margin, suggesting a high vulnerability to leaf hydraulic failure (HF), especially on the most productive genotypes. Modelling exercises revealed that variations in HF across genotypes were mainly associated with differences in leaf water vapour leakage when stomata are closed, reflecting contrasting growth strategies. Overall, our results provide a new perspective on the challenges of sustaining coffee production in the Amazon region under a drier and warmer climate.

Overexpression of Water-Responsive Genes Promoted by Elevated CO2 Reduces ROS and Enhances Drought Tolerance in Coffea Species

Drought is a major constraint to plant growth and productivity worldwide and will aggravate as water availability becomes scarcer. Although elevated air [CO₂] might mitigate some of these effects in plants, the mechanisms underlying the involved responses are poorly understood in woody economically important crops such as Coffea. This study analyzed transcriptome changes in Coffea canephora cv. CL153 and C. arabica cv. Icatu exposed to moderate (MWD) or severe water deficits (SWD) and grown under ambient (aCO₂) or elevated (eCO₂) air [CO₂]. We found that changes in expression levels and regulatory pathways were barely affected by MWD, while the SWD condition led to a down-regulation of most differentially expressed genes (DEGs). eCO₂ attenuated the impacts of drought in the transcripts of both genotypes but mostly in Icatu, in agreement with physiological and metabolic studies. A predominance of protective and reactive oxygen species (ROS)-scavenging-related genes, directly or indirectly associated with ABA signaling pathways, was found in Coffea responses, including genes involved in water deprivation and desiccation, such as protein phosphatases in Icatu, and aspartic proteases and dehydrins in CL153, whose expression was validated by qRT-PCR. The existence of a complex post-transcriptional regulatory mechanism appears to occur in Coffea explaining some apparent discrepancies between transcriptomic, proteomic, and physiological data in these genotypes.

Marker-Assisted Recurrent Selection for Pyramiding Leaf Rust and Coffee Berry Disease Resistance Alleles in Coffea arabica L

In this study, marker-assisted recurrent selection was evaluated for pyramiding resistance gene alleles against coffee leaf rust (CLR) and coffee berry diseases (CBD) in Coffea arabica. A total of 144 genotypes corresponding to 12 hybrid populations from crosses between eight parent plants with desired morphological and agronomic traits were evaluated. Molecular data were used for cross-certification, diversity study and resistance allele marker-assisted selection (MAS) against the causal agent of coffee leaf rust (Hemileia vastatrix) and coffee berry disease (Colletotrichum kahawae). In addition, nine morphological and agronomic traits were evaluated to determine the components of variance, select superior hybrids, and estimate genetic gain. From the genotypes evaluated, 134 were confirmed as hybrids. The genetic diversity between and within populations was 75.5% and 24.5%, respectively, and the cluster analysis revealed three primary groups. Pyramiding of CLR and CBD resistance genes was conducted in 11 genotypes using MAS. A selection intensity of 30% resulted in a gain of over 50% compared to the original population. Selected hybrids with increased gain also showed greater genetic divergence in addition to the pyramided resistance alleles. The strategies used were, therefore, efficient to select superior coffee hybrids for recurrent selection programs and could be used as a source of resistance in various crosses.

Ecological and genomic vulnerability to climate change across native populations of Robusta coffee (Coffea canephora)

The assessment of population vulnerability under climate change is crucial for planning conservation as well as for ensuring food security. Coffea canephora is, in its native habitat, an understorey tree that is mainly distributed in the lowland rainforests of tropical Africa. Also known as Robusta, its commercial value constitutes a significant revenue for many human populations in tropical countries. Comparing ecological and genomic vulnerabilities within the species' native range can provide valuable insights about habitat loss and the species' adaptive potential, allowing to identify genotypes that may act as a resource for varietal improvement. By applying species distribution models, we assessed ecological vulnerability as the decrease in climatic suitability under future climatic conditions from 492 occurrences. We then quantified genomic vulnerability (or risk of maladaptation) as the allelic composition change required to keep pace with predicted climate change. Genomic vulnerability was estimated from genomic environmental correlations throughout the native range. Suitable habitat was predicted to diminish to half its size by 2050, with populations near coastlines and around the Congo River being the most vulnerable. Whole-genome sequencing revealed 165 candidate SNPs associated with climatic adaptation in C. canephora, which were located in genes involved in plant response to biotic and abiotic stressors. Genomic vulnerability was higher for populations in West Africa and in the region at the border between DRC and Uganda. Despite an overall low correlation between genomic and ecological vulnerability at broad scale, these two components of vulnerability overlap spatially in ways that may become damaging. Genomic vulnerability was estimated to be 23% higher in populations where habitat will be lost in 2050 compared to regions where habitat will remain suitable. These results highlight how ecological and genomic vulnerabilities are relevant when planning on how to cope with climate change regarding an economically important species.

Coffea cytogenetics: from the first karyotypes to the meeting with genomics

Coffea karyotype organization and evolution has been uncovered by classical cytogenetics and cytogenomics. We revisit these discoveries and present new karyotype data. Coffea possesses ~ 124 species, including C. arabica and C. canephora responsible for commercial coffee production. We reviewed the Coffea cytogenetics, from the first chromosome counting, encompassing the karyotype characterization, chromosome DNA content, and mapping of chromosome portions and DNA sequences, until the integration with genomics. We also showed new data about Coffea karyotype. The 2n chromosome number evidenced the diploidy of almost all Coffea, and the C. arabica tetraploidy, as well as the polyploidy of other hybrids. Since then, other genomic similarities and divergences among the Coffea have been shown by karyotype morphology, nuclear and chromosomal C-value, AT and GC rich chromosome portions, and repetitive sequence and gene mapping. These cytogenomic data allowed us to know and understand the phylogenetic relations in Coffea, as well as their ploidy level and genomic origin, highlighting the relatively recent allopolyploidy. In addition to the euploidy, the role of the mobile elements in Coffea diversification is increasingly more evident, and the comparative analysis of their structure and distribution on the genome of different species is in the spotlight for future research. An integrative look at all these data is fundamental for a deeper understanding of Coffea karyotype evolution, including the key role of polyploidy in C. arabica origin. The 'Híbrido de Timor', a recent natural allotriploid, is also in the spotlight for its potential as a source of resistance genes and model for plant polyploidy research. Considering this, we also present some unprecedented results about the exciting evolutionary history of these polyploid Coffea.

Screening coffee genotypes for brown eye spot resistance in Brazil

Several researchers have attempted to develop coffee plants that are resistant to brown eye spot (BES); however, no coffee cultivars are resistant to the disease. In the present study, a blend of strains from Cercospora coffeicola was inoculated into 19 Brazilian commercial cultivars and 41 accessions from the Germplasm Collection of Minas Gerais to evaluate the genetic resistance ability within the population and select superior genotypes for the breeding program. After predicting the genotypic values of the estudied material, the evaluations number necessary for selecting genotypes with accuracy and efficiency was determined based on the data of severity to BES. The action of defense mechanisms plant was also investigated by assessing the levels of total soluble phenolic compounds and soluble lignin in contrasting genotypes for disease susceptibility. Based on the results, the accession MG 1207 Sumatra, had an intrinsic genetic capacity to maintain low levels of severity to BES. The genotype MG 1207 Sumatra can substantially contribute to the development of new cultivars, which may lead to the reduced use of pesticides. According to the accuracy and efficiency results obtained, four evaluations BES severity are sufficient to achieve accuracy, providing expressive genetic gains. Finally, the levels of lignin and phenolic compounds were not found to be associated with the resistance of coffee genotypes to BES.

Whole-genome resequencing of Coffea arabica L. (Rubiaceae) genotypes identify SNP and unravels distinct groups showing a strong geographical pattern

BACKGROUND

Coffea arabica L. is an economically important agricultural crop and the most popular beverage worldwide. As a perennial crop with recalcitrant seed, conservation of the genetic resources of coffee can be achieved through the complementary approach of in-situ and ex-situ field genebank. In Ethiopia, a large collection of C. arabica L. germplasm is preserved in field gene banks. Here, we report the whole-genome resequencing of 90 accessions from Choche germplasm bank representing garden and forest-based coffee production systems using Illumina sequencing technology.

RESULTS

The genome sequencing generated 6.41 billion paired-end reads, with a mean of 71.19 million reads per sample. More than 93% of the clean reads were mapped onto the C. arabica L. reference genome. A total of 11.08 million variants were identified, among which 9.74 million (87.9%) were SNPs (Single nucleotide polymorphisms) and 1.34 million (12.1%) were InDels. In all accessions, genomic variants were unevenly distributed across the coffee genome. The phylogenetic analysis using the SNP markers displayed distinct groups.

CONCLUSIONS

Resequencing of the coffee accessions has allowed identification of genetic markers, such as SNPs and InDels. The SNPs discovered in this study might contribute to the variation in important pathways of genes for important agronomic traits such as caffeine content, yield, disease, and pest in coffee. Moreover, the genome resequencing data and the genetic markers identified from 90 accessions provide insight into the genetic variation of the coffee germplasm and facilitate a broad range of genetic studies.

LysM receptors in Coffea arabica: Identification, characterization, and gene expression in response to Hemileia vastatrix

Pathogen‐associated molecular patterns (PAMPs) are recognized by pattern recognition receptors (PRRs) localized on the host plasma membrane. These receptors activate a broad-spectrum and durable defense, which are desired characteristics for disease resistance in plant breeding programs. In this study, candidate sequences for PRRs with lysin motifs (LysM) were investigated in the Coffea arabica genome. For this, approaches based on the principle of sequence similarity, conservation of motifs and domains, phylogenetic analysis, and modulation of gene expression in response to Hemileia vastatrix were used. The candidate sequences for PRRs in C. arabica (Ca1-LYP, Ca2-LYP, Ca1-CERK1, Ca2-CERK1, Ca-LYK4, Ca1-LYK5 and Ca2-LYK5) showed high similarity with the reference PRRs used: Os-CEBiP, At-CERK1, At-LYK4 and At-LYK5. Moreover, the ectodomains of these sequences showed high identity or similarity with the reference sequences, indicating structural and functional conservation. The studied sequences are also phylogenetically related to the reference PRRs described in Arabidopsis, rice, and other plant species. All candidates for receptors had their expression induced after the inoculation with H. vastatrix, since the first time of sampling at 6 hours post‐inoculation (hpi). At 24 hpi, there was a significant increase in expression, for most of the receptors evaluated, and at 48 hpi, a suppression. The results showed that the candidate sequences for PRRs in the C. arabica genome display high homology with fungal PRRs already described in the literature. Besides, they respond to pathogen inoculation and seem to be involved in the perception or signaling of fungal chitin, acting as receptors or co-receptors of this molecule. These findings represent an advance in the understanding of the basal immunity of this species.

Designing the best breeding strategy for Coffea canephora: Genetic evaluation of pure and hybrid individuals aiming to select for productivity and disease resistance traits

Breeding programs of the species Coffea canephora rely heavily on the significant genetic variability between and within its two varietal groups (conilon and robusta). The use of hybrid families and individuals has been less common. The objectives of this study were to evaluate parents and families from the populations of conilon, robusta, and its hybrids and to define the best breeding and selection strategies for productivity and disease resistance traits. As such, 71 conilon clones, 56 robusta clones, and 20 hybrid families were evaluated over several years for the following traits: vegetative vigor, incidence of rust and cercosporiosis, fruit ripening time, fruit size, plant height, canopy diameter, and yield per plant. Components of variance and genetic parameters were estimated via residual maximum likelihood (REML) and genotypic values were predicted via best linear unbiased prediction (BLUP). Genetic variability among parents (clones) and hybrid families was detected for most of the evaluated traits. The Mulamba-Rank index suggests potential gains up to 17% for the genotypic aggregate of traits in the hybrid population. An intrapopulation recurrent selection within the hybrid population would be the best breeding strategy because the genetic variability, narrow and broad senses heritabilities and selective accuracies for important traits were maximized in the crossed population. Besides, such strategy is simple, low cost and quicker than the concurrent reciprocal recurrent selection in the two parental populations, and this maximizes the genetic gain for unit of time.

Genetic diversity of wild and cultivated Coffea canephora in northeastern DR Congo and the implications for conservation

PREMISE

Many cultivated coffee varieties descend from Coffea canephora, commonly known as Robusta coffee. The Congo Basin has a century-long history of Robusta coffee cultivation and breeding, and is hypothesized to be the region of origin of many of the cultivated Robusta varieties. Since little is known about the genetic composition of C. canephora in this region, we assessed the genetic diversity of wild and cultivated C. canephora shrubs in the Democratic Republic of the Congo.

METHODS

Using 18 microsatellite markers, we studied the genetic composition of wild and backyard-grown C. canephora shrubs in the Tshopo and Ituri provinces and multiple accessions from the INERA Yangambi Coffee Collection. We assessed genetic clustering patterns, genetic diversity, and genetic differentiation between populations.

RESULTS

Genetic differentiation was relatively strong between wild and cultivated C. canephora shrubs, and both gene pools harbored multiple unique alleles. Strong genetic differentiation was also observed between wild populations. The level of genetic diversity in wild populations was similar to that of the INERA Yangambi Coffee Collection, but local wild genotypes were mostly missing from that collection. Shrubs grown in the backyards were genetically similar to the breeding material from INERA Yangambi.

CONCLUSIONS

Most C. canephora that is grown in local backyards originated from INERA breeding programs, while a few shrubs were obtained directly from surrounding forests. The INERA Yangambi Coffee Collection could benefit from an enrichment with local wild genotypes to increase the genetic resources available for breeding purposes and to support ex situ conservation.

Genomic and transcriptomic inventory of membrane transporters in coffee: Exploring molecular mechanisms of metabolite accumulation

The genus Coffea (Rubiaceae) encompasses a group of perennial plant species, including a commodity crop from which seeds are roasted, ground, and infused to make one of the most appreciated beverages in the world. As an important tropical crop restricted to specific regions of the world, coffee production is highly susceptible to the effects of environmental instabilities (i.e., local year-to-year weather fluctuations and global climate change) and threatening pest pressures, not to mention an increasing quality rigor by consumers in industrialized countries. Specialized metabolites are substances that largely affect plant-environment interactions as well as how consumers experience agricultural products. Membrane transporters are key targets, albeit understudied, for understanding and tailoring the spatiotemporal distribution of specialized metabolites as they mediate and control molecular trafficking and substance accumulation. Therefore, we analyzed the transportome of C. canephora encoded within the 25,574 protein-coding genes annotated in the genome of this species and identified 1847 putative membrane transporters. Following, we mined 152 transcriptional profiles of C. canephora and C. arabica and performed a comprehensive co-expression analysis to identify transporters potentially involved in the accumulation of specialized metabolites associated with beverage quality and bioactivity attributes. In toto, this report points to an avenue of possibilities on Coffea genomic and transcriptomic data mining for genetic breeding strategies, which can lead to the development of new, resilient varieties for more sustainable coffee production systems.

Multi-Approach Analysis Reveals Local Adaptation in a Widespread Forest Tree of Reunion Island

Detecting processes of local adaptation in forest trees and identifying environmental selective drivers are of primary importance for forest management and conservation. Transplant experiments, functional genomics and population genomics are complementary tools to efficiently characterize heritable phenotypic traits and to decipher the genetic bases of adaptive traits. Using an integrative approach combining phenotypic assessment in common garden, transcriptomics and landscape genomics, we investigated leaf adaptive traits in Coffea mauritiana, a forest tree endemic to Reunion Island. Eight populations of C. mauritiana originating from sites with contrasted environmental conditions were sampled in common garden to assess several leaf morphological traits, to analyze the leaf transcriptome and leaf cuticular wax composition. The relative alkane content of cuticular waxes was significantly correlated with major climatic gradients, paving the way for further transcriptome-based analyses. The expression pattern of cuticle biosynthetic genes was consistent with a modulation of alkane accumulation across the population studied, supporting the hypothesis that the composition of cuticular wax is involved in the local adaptation of C. mauritiana. Association tests in landscape genomics performed using RNA-seq-derived single-nucleotide polymorphisms revealed that genes associated with cell wall remodeling also likely play an adaptive role. By combining these different approaches, this study efficiently identified local adaptation processes in a non-model species. Our results provide the first evidence for local adaptation in trees endemic to Reunion Island and highlight the importance of cuticle composition for the adaptation of trees to the high evaporative demand in warm climates.

Lipidomic and metabolomic profiles of Coffea canephora L. beans cultivated in Southwestern Nigeria

Coffee (Coffea spp.) is one of the most popular refreshing beverages globally. Coffee lipid diversity has untapped potential for improving coffee marketability because lipids contribute significantly to both the health benefits and cup quality of coffee. However, in spite of its potential importance, there have not been extensive studies of lipids among C. canephora genotypes. In this study, ultra-performance liquid chromatography coupled with mass spectrometry (UPLC-MS) profiling of lipid molecules was performed for 30 genotypes consisting of 15 cultivated and 15 conserved genotypes of C. canephora in Southwestern Nigeria. We identified nine classes of lipids in the 30 genotypes which belong to the 'Niaouli', 'Kouillou' and 'Java Robusta' group: among these, the most abundant lipid class was the triacylglycerols, followed by the fatty acyls group. Although 'Niaouli' diverged from the 'Kouillou' and 'Java Robusta' genotypes when their lipid profiles were compared, there was greater similarity in their lipid composition by multivariate analysis, compared to that observed when their primary metabolites and especially their secondary metabolite profiles were examined. However, distinctions could be made among genotypes. Members of the fatty acyls group had the greatest power to discriminate among genotypes, however, lipids that were low in abundance e.g. a cholesterol ester (20:3), and phosphotidylethanolamine (34:0) were also helpful to understand the relationships among C. canephora genotypes. The divergent lipid profiles identified among the C. canephora genotypes, correlated with their Single Nucleotide Polymorphism grouping as assessed by genotype-by-sequencing, and will be exploited to improve coffee cup quality.

Photoperiod-dependent transcriptional modifications in key metabolic pathways in Coffea arabica

Photoperiod length induces in temperate plants major changes in growth rates, morphology and metabolism with, for example, modifications in the partitioning of photosynthates to avoid starvation at the end of long nights. However, this has never been studied for a tropical perennial species adapted to grow in a natural photoperiod close to 12 h/12 h all year long. We grew Coffea arabica L., an understorey perennial evergreen tropical species in its natural 12 h/12 h and in a short 8 h/16 h photoperiod, and we investigated its responses at the physiological, metabolic and transcriptomic levels. The expression pattern of rhythmic genes, including core clock genes, was affected by changes in photoperiod. Overall, we identified 2859 rhythmic genes, of which 89% were also rhythmic in Arabidopsis thaliana L. Under short-days, plant growth was reduced, and leaves were thinner with lower chlorophyll content. In addition, secondary metabolism was also affected with chlorogenic acid and epicatechin levels decreasing, and in agreement, the genes involved in lignin synthesis were overexpressed and those involved in the flavanol pathway were underexpressed. Our results show that the 8 h/16 h photoperiod induces drastic changes in morphology, metabolites and gene expression, and the responses for gene expression are similar to those observed in the temperate annual A. thaliana species. Short photoperiod induces drastic changes in gene expression, metabolites and leaf structure, some of these responses being similar to those observed in A. thaliana.

N-Methyltransferases of Caffeine Biosynthetic Pathway in Plants

Caffeine (Cf) is one of the important components of plant-derived drinks, such as tea, coffee, and cola. It can protect soft tissues from being infected by pathogens and is also medically beneficial for human health. In this review, we first introduced the Cf biosynthesis pathways in plants and the related N-methyltransferases (NMTs), with a focus on the current research status of the substrate specificity, structural basis for substrate recognition, and catalytic mechanism in members of the caffeine synthase gene family. In addition, we addressed the expression characteristics and potential regulatory mechanisms of NMTs and also projected the future research directions. The goal was to summarize the Cf biosynthetic pathway and related NMTs in plants and to provide the molecular basis for regulating the caffeine biosynthesis, so as to effectively guide future tea and coffee breeding.

Transcriptomic Leaf Profiling Reveals Differential Responses of the Two Most Traded Coffee Species to Elevated [CO2]

As atmospheric [CO₂] continues to rise to unprecedented levels, understanding its impact on plants is imperative to improve crop performance and sustainability under future climate conditions. In this context, transcriptional changes promoted by elevated CO₂ (eCO₂) were studied in genotypes from the two major traded coffee species: the allopolyploid Coffea arabica (Icatu) and its diploid parent, C. canephora (CL153). While Icatu expressed more genes than CL153, a higher number of differentially expressed genes were found in CL153 as a response to eCO₂. Although many genes were found to be commonly expressed by the two genotypes under eCO₂, unique genes and pathways differed between them, with CL153 showing more enriched GO terms and metabolic pathways than Icatu. Divergent functional categories and significantly enriched pathways were found in these genotypes, which altogether supports contrasting responses to eCO₂. A considerable number of genes linked to coffee physiological and biochemical responses were found to be affected by eCO₂ with the significant upregulation of photosynthetic, antioxidant, and lipidic genes. This supports the absence of photosynthesis down-regulation and, therefore, the maintenance of increased photosynthetic potential promoted by eCO₂ in these coffee genotypes.

Coffea arabica seedlings genotypes are tolerant to high induced selenium stress: Evidence from physiological plant responses and antioxidative performance

Selenium (Se) is considered a beneficial element to higher plants based on its regulation of antioxidative system under abiotic or biotic stresses. However, the limit of beneficial and toxic physiological effects of Se is very narrow. In the present study, the antioxidant performance, nutritional composition, long-distance transport of Se, photosynthetic pigments, and growth of Coffea arabica genotypes in response to Se concentration in solution were evaluated. Five Coffea arabica genotypes (Obatã, IPR99, IAC125, IPR100 and Catucaí) were used, which were grown in the absence and presence of Se (0 and 1.0 mmol L^-1) in nutrient solution. The application of 1 mmol L^-1 Se promoted root browning in all genotypes. There were no visual symptoms of leaf toxicity, but there was a reduction in the concentration of phosphorus and sulfur in the shoots of plants exposed to high Se concentration. Except for genotype Obatã, the coffee seedlings presented strategies for regulating Se uptake by reducing long-distance transport of Se from roots to shoots. The concentrations of total chlorophyll, total pheophytin, and carotenoids were negatively affected in genotypes Obatã, IPR99, and IAC125 upon exposure to Se at 1 mmol L^-1. H₂O₂ production was reduced in genotypes IPR99, IPR100, and IAC125 upon exposure to Se, resulting in lower activity of superoxide dismutase (SOD), and catalase (CAT). These results suggest that antioxidant metabolism was effective in regulating oxidative stress in plants treated with Se. The increase in sucrose, and decrease in SOD, CAT and ascorbate peroxidase (APX) activities, as well as Se compartmentalization in the roots, were the main biochemical and physiological modulatory effects of coffee seedlings under stress conditions due to excess of Se.

YUCCA-Mediated Biosynthesis of the Auxin IAA Is Required during the Somatic Embryogenic Induction Process in Coffea canephora

Despite the existence of considerable research on somatic embryogenesis (SE), the molecular mechanism that regulates the biosynthesis of auxins during the SE induction process remains unknown. Indole-3-acetic acid (IAA) is an auxin that is synthesized in plants through five pathways. The biosynthetic pathway most frequently used in this synthesis is the conversion of tryptophan to indol-3-pyruvic acid (IPA) by tryptophan aminotransferase of Arabidopsis (TAA) followed by the conversion of IPA to IAA by enzymes encoded by YUCCA (YUC) genes of the flavin monooxygenase family; however, it is unclear whether YUC-mediated IAA biosynthesis is involved in SE induction. In this study, we report that the increase of IAA observed during SE pre-treatment (plants in MS medium supplemented with 1-naphthaleneacetic acid (NAA) 0.54 µM and kinetin (Kin) 2.32 µM for 14 days) was due to its de novo biosynthesis. By qRT-PCR, we demonstrated that YUC gene expression was consistent with the free IAA signal found in the explants during the induction of SE. In addition, the use of yucasin to inhibit the activity of YUC enzymes reduced the signal of free IAA in the leaf explants and dramatically decreased the induction of SE. The exogenous addition of IAA restored the SE process in explants treated with yucasin. Our findings suggest that the biosynthesis and localization of IAA play an essential role during the induction process of SE in Coffea canephora.

Resequencing 93 accessions of coffee unveils independent and parallel selection during Coffea species divergence

Coffea arabica, C. canephora and C. excelsa, with differentiated morphological traits and distinct agro-climatic conditions, compose the majority of the global coffee plantation. To comprehensively understand their genetic diversity and divergence for future genetic improvement requires high-density markers. Here, we sequenced 93 accessions encompassing these three Coffea species, uncovering 15,367,960 single-nucleotide polymorphisms (SNPs). These SNPs are unequally distributed across different genomic regions and gene families, with two disease-resistant gene families showing the highest SNP density, suggesting strong balancing selection. Meanwhile, the allotetraploid C. arabica exhibits greater nucleotide diversity, followed by C. canephora and C. excelsa. Population divergence (F_ST), population stratification and phylogeny all support strong divergence among species, with C. arabica and its parental species C. canephora being closer genetically. Scanning of genomic islands with elevated F_ST and structure-disruptive SNPs contributing to species divergence revealed that most of the selected genes in each lineage are independent, with a few being selected in parallel for two or three species, such as genes in root hair cell development, flavonols accumulation and disease-resistant genes. Moreover, some of the SNPs associated with coffee lipids exhibit significantly biased allele frequency among species, being valuable for interspecific breeding. Overall, our study not only uncovers the key population genomic patterns among species but also contributes a substantial genomic resource for coffee breeding.

Chloroplast genomes of Rubiaceae: Comparative genomics and molecular phylogeny in subfamily Ixoroideae

In Rubiaceae phylogenetics, the number of markers often proved a limitation with authors failing to provide well-supported trees at tribal and generic levels. A robust phylogeny is a prerequisite to study the evolutionary patterns of traits at different taxonomic levels. Advances in next-generation sequencing technologies have revolutionized biology by providing, at reduced cost, huge amounts of data for an increased number of species. Due to their highly conserved structure, generally recombination-free, and mostly uniparental inheritance, chloroplast DNA sequences have long been used as choice markers for plant phylogeny reconstruction. The main objectives of this study are: 1) to gain insight in chloroplast genome evolution in the Rubiaceae (Ixoroideae) through efficient methodology for de novo assembly of plastid genomes; and, 2) to test the efficiency of mining SNPs in the nuclear genome of Ixoroideae based on the use of a coffee reference genome to produce well-supported nuclear trees. We assembled whole chloroplast genome sequences for 27 species of the Rubiaceae subfamily Ixoroideae using next-generation sequences. Analysis of the plastid genome structure reveals a relatively good conservation of gene content and order. Generally, low variation was observed between taxa in the boundary regions with the exception of the inverted repeat at both the large and short single copy junctions for some taxa. An average of 79% of the SNP determined in the Coffea genus are transferable to Ixoroideae, with variation ranging from 35% to 96%. In general, the plastid and the nuclear genome phylogenies are congruent with each other. They are well-resolved with well-supported branches. Generally, the tribes form well-identified clades but the tribe Sherbournieae is shown to be polyphyletic. The results are discussed relative to the methodology used and the chloroplast genome features in Rubiaceae and compared to previous Rubiaceae phylogenies.

Seed comparative genomics in three coffee species identify desiccation tolerance mechanisms in intermediate seeds

In contrast to desiccation-tolerant 'orthodox' seeds, so-called 'intermediate' seeds cannot survive complete drying and are short-lived. All species of the genus Coffea produce intermediate seeds, but they show a considerable variability in seed desiccation tolerance (DT), which may help to decipher the molecular basis of seed DT in plants. We performed a comparative transcriptome analysis of developing seeds in three coffee species with contrasting desiccation tolerance. Seeds of all species shared a major transcriptional switch during late maturation that governs a general slow-down of metabolism. However, numerous key stress-related genes, including those coding for the late embryogenesis abundant protein EM6 and the osmosensitive calcium channel ERD4, were up-regulated during DT acquisition in the two species with high seed DT, C. arabica and C. eugenioides. By contrast, we detected up-regulation of numerous genes involved in the metabolism, transport, and perception of auxin in C. canephora seeds with low DT. Moreover, species with high DT showed a stronger down-regulation of the mitochondrial machinery dedicated to the tricarboxylic acid cycle and oxidative phosphorylation. Accordingly, respiration measurements during seed dehydration demonstrated that intermediate seeds with the highest DT are better prepared to cease respiration and avoid oxidative stresses.

Full moonlight-induced circadian clock entrainment in Coffea arabica

BACKGROUND

It is now well documented that moonlight affects the life cycle of invertebrates, birds, reptiles, and mammals. The lunisolar tide is also well-known to alter plant growth and development. However, although plants are known to be very photosensitive, few studies have been undertaken to explore the effect of moonlight on plant physiology.

RESULTS

Here for the first time we report a massive transcriptional modification in Coffea arabica genes under full moonlight conditions, particularly at full moon zenith and 3 h later. Among the 3387 deregulated genes found in our study, the main core clock genes were affected.

CONCLUSIONS

Moonlight also negatively influenced many genes involved in photosynthesis, chlorophyll biosynthesis and chloroplast machinery at the end of the night, suggesting that the full moon has a negative effect on primary photosynthetic machinery at dawn. Moreover, full moonlight promotes the transcription of major rhythmic redox genes and many heat shock proteins, suggesting that moonlight is perceived as stress. We confirmed this huge impact of weak light (less than 6 lx) on the transcription of circadian clock genes in controlled conditions mimicking full moonlight.

Development of an Efficient Protocol to Obtain Transgenic Coffee, Coffea arabica L., Expressing the Cry10Aa Toxin of Bacillus thuringiensis

This report presents an efficient protocol of the stable genetic transformation of coffee plants expressing the Cry10Aa protein of Bacillus thuringiensis. Embryogenic cell lines with a high potential of propagation, somatic embryo maturation, and germination were used. Gene expression analysis of cytokinin signaling, homedomains, auxin responsive factor, and the master regulators of somatic embryogenesis genes involved in somatic embryo maturation were evaluated. Plasmid pMDC85 containing the cry10Aa gene was introduced into a Typica cultivar of C. arabica L. by biobalistic transformation. Transformation efficiency of 16.7% was achieved, according to the number of embryogenic aggregates and transgenic lines developed. Stable transformation was proven by hygromycin-resistant embryogenic lines, green fluorescent protein (GFP) expression, quantitative analyses of Cry10Aa by mass spectrometry, Western blot, ELISA, and Southern blot analyses. Cry10Aa showed variable expression levels in somatic embryos and the leaf tissue of transgenic plants, ranging from 76% to 90% of coverage of the protein by mass spectrometry and from 3.25 to 13.88 μg/g fresh tissue, with ELISA. qPCR-based 2^−ΔΔCt trials revealed high transcription levels of cry10Aa in somatic embryos and leaf tissue. This is the first report about the stable transformation and expression of the Cry10Aa protein in coffee plants with the potential for controlling the coffee berry borer.

Plant population dynamics on oceanic islands during the Late Quaternary climate changes: genetic evidence from a tree species (Coffea mauritiana) in Reunion Island

Past climatic fluctuations have played a major role in shaping the current plant biodiversity. Although harbouring an exceptional biota, oceanic islands have received little attention in studies on species demographic history and past vegetation patterns. We investigated the impact of past climatic changes on the effective population size of a tree (Coffea mauritiana) that is endemic to Reunion Island, located in the south-western Indian Ocean (SWIO). Demographic changes were inferred using summary statistics calculated from genomic data. Using ecological niche modelling and the current distribution of genetic diversity, the paleodistribution of the species was also assessed. A reduction in the effective population size of C. mauritiana during the last glaciation maximum was inferred. The distribution of the species was reduced on the western side of the island, due to low rainfall. It appeared that a major reduction in rainfall and a slight temperature decrease prevailed in the SWIO. Our findings indicated that analyses on the current patterns of intraspecific genetic variations can efficiently contribute to past climatic changes characterisation in remote islands. Identifying area with higher resilience in oceanic islands could provide guidance in forest management and conservation faced to the global climate change.

Fragmentation of Pooled PCR Products for Highly Multiplexed TILLING

Improvements to massively parallel sequencing have allowed the routine recovery of natural and induced sequence variants. A broad range of biological disciplines have benefited from this, ranging from plant breeding to cancer research. The need for high sequence coverage to accurately recover single nucleotide variants and small insertions and deletions limits the applicability of whole genome approaches. This is especially true in organisms with a large genome size or for applications requiring the screening of thousands of individuals, such as the reverse-genetic technique known as TILLING. Using PCR to target and sequence chosen genomic regions provides an attractive alternative as the vast reduction in interrogated bases means that sample size can be dramatically increased through amplicon multiplexing and multi-dimensional sample pooling while maintaining suitable coverage for recovery of small mutations. Direct sequencing of PCR products is limited, however, due to limitations in read lengths of many next generation sequencers. In the present study we show the optimization and use of ultrasonication for the simultaneous fragmentation of multiplexed PCR amplicons for TILLING highly pooled samples. Sequencing performance was evaluated in a total of 32 pooled PCR products produced from 4096 chemically mutagenized Hordeum vulgare DNAs pooled in three dimensions. Evaluation of read coverage and base quality across amplicons suggests this approach is suitable for high-throughput TILLING and other applications employing highly pooled complex sampling schemes. Induced mutations previously identified in a traditional TILLING screen were recovered in this dataset further supporting the efficacy of the approach.

Development and evaluation of a genome-wide Coffee 8.5K SNP array and its application for high-density genetic mapping and for investigating the origin of Coffea arabica L

Coffee species such as Coffea canephora P. (Robusta) and C. arabica L. (Arabica) are important cash crops in tropical regions around the world. C. arabica is an allotetraploid (2n = 4x = 44) originating from a hybridization event of the two diploid species C. canephora and C. eugenioides (2n = 2x = 22). Interestingly, these progenitor species harbour a greater level of genetic variability and are an important source of genes to broaden the narrow Arabica genetic base. Here, we describe the development, evaluation and use of a single-nucleotide polymorphism (SNP) array for coffee trees. A total of 8580 unique and informative SNPs were selected from C. canephora and C. arabica sequencing data, with 40% of the SNP located in annotated genes. In particular, this array contains 227 markers associated to 149 genes and traits of agronomic importance. Among these, 7065 SNPs (~82.3%) were scorable and evenly distributed over the genome with a mean distance of 54.4 Kb between markers. With this array, we improved the Robusta high-density genetic map by adding 1307 SNP markers, whereas 945 SNPs were found segregating in the Arabica mapping progeny. A panel of C. canephora accessions was successfully discriminated and over 70% of the SNP markers were transferable across the three species. Furthermore, the canephora-derived subgenome of C. arabica was shown to be more closely related to C. canephora accessions from northern Uganda than to other current populations. These validated SNP markers and high-density genetic maps will be useful to molecular genetics and for innovative approaches in coffee breeding.

Genetic diversity and population divergences of an indigenous tree (Coffea mauritiana) in Reunion Island: role of climatic and geographical factors

Oceanic islands are commonly considered as natural laboratories for studies on evolution and speciation. The evolutionary specificities of islands associated with species biology provide unique scenarios to study the role of geography and climate in driving population divergence. However, few studies have addressed this subject in small oceanic islands with heterogeneous climates. Being widely distributed in Reunion Island forest, Coffea mauritiana represents an interesting model case for investigating patterns of within-island differentiation at small spatial scale. In this study, we examined the genetic diversity and population divergences of C. mauritiana using SNP markers obtained from 323 individuals across 34 locations in Reunion Island. Using redundancy analysis, we further evaluated the contribution of geographic and climatic factors to shaping genetic divergence among populations. Genetic diversity analyses revealed that accessions clustered according to the source population, with further grouping in regional clusters. Genetic relationships among the regional clusters underlined a recent process of expansion in the form of step-by-step colonization on both sides of the island. Divergence among source populations was mostly driven by the joint effect of geographic distance and climatic heterogeneity. The pattern of isolation-by-geography was in accordance with the dispersal characteristics of the species, while isolation-by-environment was mostly explained by the heterogeneous rainfall patterns, probably associated with an asynchronous flowering among populations. These findings advance our knowledge on the patterns of genetic diversity and factors of population differentiation of species native to Reunion Island, and will also usefully guide forest management for conservation.

Genome-wide identification and expression analysis of HKT transcription factor under salt stress in nine plant species

High-affinity K+ (HKT) gene family is regulated the transport of Na+ and maintain the balance between Na+ and K+ in the process of plant growth, development and response to abiotic stress. Despite this fact, systemic and comprehensive studies on HKT in multiply plants remains unknown. A total of 29 HKT genes distributed on nine species were identified. Phylogenetic tree analysis results indicated that HKT genes were divided into five homology subfamilies. Combining structural analysis with protein contains five highly conservative motifs, HKT family has similar gene structures and special gene characteristics. Finally, the expression patterns of HKT showed two different dramatic changes in different organs and tissues under different salt stress in multiply plants. This study has many implications for research into the comparative genomics analysis of HTK gene family, which revealed regulation mechanism of HKT genes, is valuable for understanding development and response to abiotic stress in plant.

BDBM 1.0: A Desktop Application for Efficient Retrieval and Processing of High-Quality Sequence Data and Application to the Identification of the Putative Coffea S-Locus

Nowadays, bioinformatics is one of the most important areas in modern biology and the creation of high-quality scientific software supporting this recent research area is one of the core activities of many researchers. In this context, high-quality sequence datasets are needed to perform inferences on the evolution of species, genes, and gene families, or to get evidence for adaptive amino acid evolution, among others. Nevertheless, sequence data are very often spread over several databases, many useful genomes and transcriptomes are non-annotated, the available annotation is not for the desired coding sequence isoform, and/or is unlikely to be accurate. Moreover, although the FASTA text-based format is quite simple and usable by most software applications, there are a number of issues that may be critical depending on the software used to analyse such files. Therefore, researchers without training in informatics often use a fraction of all available data. The above issues can be addressed using already available software applications, but there is no easy-to-use single piece of software that allows performing all these tasks within the same graphical interface, such as the one here presented, named BDBM (Blast DataBase Manager). BDBM can be used to efficiently get gene sequences from annotated and non-annotated genomes and transcriptomes. Moreover, it can be used to look for alternatives to existing annotations and to easily create reliable custom databases. Such databases are essential to prepare high-quality datasets. The analyses that we have performed on the Coffea canephora genome using BDBM aimed at the identification of the S-locus region (that harbours the genes involved in gametophytic self-incompatibility) led to the conclusion that there are two likely regions, one on chromosome 2 (around region 6600000-6650000), and another on chromosome 5 (around 15830000-15930000). Such findings are discussed in the context of the Rubiaceae gametophytic self-incompatibility evolution.

Identification of Hydroxypyrazine O-Methyltransferase Genes in Coffea arabica: A Potential Source of Methoxypyrazines That Cause Potato Taste Defect

The goal of this study is to identify Coffea arabica O-methyltransferase (OMT) genes involved in the biosynthesis of methoxypyrazines. High levels of 2-isopropyl-3-methoxypyrazine (IPMP) and 2-isobutyl-3-methoxypyrazine (IBMP) in coffee beans are associated with the potato taste defect (PTD). Among the 34 putative O-methyltransferase genes identified in the published genome of C. canephora, three genes are highly homologous to known hydroxypyrazine OMT genes. Genes of interest were amplified and sequenced from genomic DNA of single C. arabica beans grown in eight different locations, including regions with endemic PTD. Although C. arabica OMT target sequences were almost identical regardless of source location, individual beans shared numerous polymorphisms in each of the target genes. Two of the predicted C. arabica OMT enzymes were successfully expressed in Escherichia coli and purified, and one enzyme shows slow yet measurable turnover of both 3-isobutyl-2-hydroxypyrazine (IBHP) and 3-isopropyl-2- hydroxypyrazine (IPHP), supporting a possible role of the coffee plant in PTD.