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

A chromosome-scale assembly reveals chromosomal aberrations and exchanges generating genetic diversity in Coffea arabica germplasm

In order to better understand the mechanisms generating genetic diversity in the recent allotetraploid species Coffea arabica, here we present a chromosome-level assembly obtained with long read technology. Two genomic compartments with different structural and functional properties are identified in the two homoeologous genomes. The resequencing data from a large set of accessions reveals low intraspecific diversity in the center of origin of the species. Across a limited number of genomic regions, diversity increases in some cultivated genotypes to levels similar to those observed within one of the progenitor species, Coffea canephora, presumably as a consequence of introgressions deriving from the so-called Timor hybrid. It also reveals that, in addition to few, early-occurring exchanges between homoeologous chromosomes, there are numerous recent chromosomal aberrations including aneuploidies, deletions, duplications and exchanges. These events are still polymorphic in the germplasm and could represent a fundamental source of genetic variation in such a lowly variable species.

Genotyping-by-sequencing provides new genetic and taxonomic insights in the critical group of Centaurea tenorei

Centaurea L. is one of the most widespread, differentiated, and critical genera of Asteraceae in the Euro-Mediterranean area, with more than 100 currently recognized species inhabiting the region. The controversial C. tenorei group, narrowly endemic to the Peninsula of Sorrento (Campania region, southern Italy), includes three weakly differentiated microspecies: C. tenorei Guss. ex Lacaita, C. montaltensis (Fiori) Peruzzi and C. lacaitae Peruzzi. However, their taxonomic distinctiveness and relationships with close or sympatric species are still unclear. In particular, the existence in several localities of individuals with intermediate morphology suggests inadequate taxonomic assessment within the group or hybridization and introgression with other species. In this study we aimed at defining population structure in this complex. With this objective, we sampled the three currently accepted species from their loci classici (i.e., the localities in which the taxa were originally described) and from other localities throughout the range, including populations of difficult identification occurring where the ranges of different taxa overlap. We employed a panel of SNPs obtained via genotyping-by-sequencing for investigations on genetic structure, admixture and ploidy inference, the latter also compared with chromosome counts. Our results showed that Centaurea tenorei s.l. is consistently tetraploid, contradicting the current taxonomy that was also based on ploidy level. Population structure analyses indicated the presence of four to seven clusters, most of which with clear evidence of admixture. Furthermore, contrarily to what previously supposed, we demonstrated a remarkable contribution of C. deusta, more that of C. cineraria in the genetic make-up of C. tenorei. However, we found a population of C. cineraria outside its ecological range, probably driven by climate change, which could be responsible in the future of further hybridization phenomena.

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.

Metabolite Profiling and Transcriptome Analysis Revealed the Conserved Transcriptional Regulation Mechanism of Caffeine Biosynthesis in Tea and Coffee Plants

Caffeine is a characteristic bioactive compound in tea and coffee plants, which is synthesized and accumulated extensively in leaves and seeds. However, little is known about the regulatory mechanism of caffeine synthesis in plants. This study compared the caffeine metabolite between tea and coffee plants. We found that tea leaves contained significantly higher caffeine than coffee leaves, which is perhaps due to more members of N-methyltransferase (NMT) genes as well as higher expression levels in tea plants. Substantial numbers of transcription factors were predicted to be involved in caffeine biosynthesis regulation, combining weighted gene co-expression network analysis and the cis-element of NMT promoter analysis in tea and coffee plants. Furthermore, analysis of the transcription factors from the caffeine-related modules suggested that the regulatory mechanism of caffeine biosynthesis was probably partly conservative in tea and coffee plants. This study provides an essential resource for the regulatory mechanism of caffeine biosynthesis in plants.

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.

TIP_finder: An HPC Software to Detect Transposable Element Insertion Polymorphisms in Large Genomic Datasets

Transposable elements (TEs) are non-static genomic units capable of moving indistinctly from one chromosomal location to another. Their insertion polymorphisms may cause beneficial mutations, such as the creation of new gene function, or deleterious in eukaryotes, e.g., different types of cancer in humans. A particular type of TE called LTR-retrotransposons comprises almost 8% of the human genome. Among LTR retrotransposons, human endogenous retroviruses (HERVs) bear structural and functional similarities to retroviruses. Several tools allow the detection of transposon insertion polymorphisms (TIPs) but fail to efficiently analyze large genomes or large datasets. Here, we developed a computational tool, named TIP_finder, able to detect mobile element insertions in very large genomes, through high-performance computing (HPC) and parallel programming, using the inference of discordant read pair analysis. TIP_finder inputs are (i) short pair reads such as those obtained by Illumina, (ii) a chromosome-level reference genome sequence, and (iii) a database of consensus TE sequences. The HPC strategy we propose adds scalability and provides a useful tool to analyze huge genomic datasets in a decent running time. TIP_finder accelerates the detection of transposon insertion polymorphisms (TIPs) by up to 55 times in breast cancer datasets and 46 times in cancer-free datasets compared to the fastest available algorithms. TIP_finder applies a validated strategy to find TIPs, accelerates the process through HPC, and addresses the issues of runtime for large-scale analyses in the post-genomic era. TIP_finder version 1.0 is available at https://github.com/simonorozcoarias/TIP_finder.