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.

The absence of the caffeine synthase gene is involved in the naturally decaffeinated status of Coffea humblotiana, a wild species from Comoro archipelago.

Caffeine is the most consumed alkaloid stimulant in the world. It is synthesized through the activity of three known N-methyltransferase proteins. Here we are reporting on the 422-Mb chromosome-level assembly of the Coffea humblotiana genome, a wild and endangered, naturally caffeine-free, species from the Comoro archipelago. We predicted 32,874 genes and anchored 88.7% of the sequence onto the 11 chromosomes. Comparative analyses with the African Robusta coffee genome (C. canephora) revealed an extensive genome conservation, despite an estimated 11 million years of divergence and a broad diversity of genome sizes within the Coffea genus. In this genome, the absence of caffeine is likely due to the absence of the caffeine synthase gene which converts theobromine into caffeine through an illegitimate recombination mechanism. These findings pave the way for further characterization of caffeine-free species in the Coffea genus and will guide research towards naturally-decaffeinated coffee drinks for consumers.

Complex evolutionary history of coffees revealed by full plastid genomes and 28,800 nuclear SNP analyses, with particular emphasis on Coffea canephora (Robusta coffee).

For decades coffees were associated with the genus Coffea. In 2011, the closely related genus Psilanthus was subsumed into Coffea. However, results obtained in 2017-based on 28,800 nuclear SNPs-indicated that there is not substantial phylogenetic support for this incorporation. In addition, a recent study of 16 plastid full-genome sequences highlighted an incongruous placement of Coffea canephora (Robusta coffee) between maternal and nuclear trees. In this study, similar global features of the plastid genomes of Psilanthus and Coffea are observed. In agreement with morphological and physiological traits, the nuclear phylogenetic tree clearly separates Psilanthus from Coffea (with exception to C. rhamnifolia, closer to Psilanthus than to Coffea). In contrast, the maternal molecular tree was incongruent with both morphological and nuclear differentiation, with four main clades observed, two of which include both Psilanthus and Coffea species, and two with either Psilanthus or Coffea species. Interestingly, Coffea and Psilanthus taxa sampled in West and Central Africa are members of the same group. Several mechanisms such as the retention of ancestral polymorphisms due to incomplete lineage sorting, hybridization leading to homoploidy (without chromosome doubling) and alloploidy (for C. arabica) are involved in the evolutionary history of the coffee species. While sharing similar morphological characteristics, the genetic relationships within C. canephora have shown that some populations are well differentiated and genetically isolated. Given the position of its closely-related species, we may also consider C. canephora to be undergoing a long process of speciation with an intermediate step of (sub-)speciation.

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.

Inpactor, Integrated and Parallel Analyzer and Classifier of LTR Retrotransposons and Its Application for Pineapple LTR Retrotransposons Diversity and Dynamics.

One particular class of Transposable Elements (TEs), called Long Terminal Repeats (LTRs), retrotransposons, comprises the most abundant mobile elements in plant genomes. Their copy number can vary from several hundreds to up to a few million copies per genome, deeply affecting genome organization and function. The detailed classification of LTR retrotransposons is an essential step to precisely understand their effect at the genome level, but remains challenging in large-sized genomes, requiring the use of optimized bioinformatics tools that can take advantage of supercomputers. Here, we propose a new tool: Inpactor, a parallel and scalable pipeline designed to classify LTR retrotransposons, to identify autonomous and non-autonomous elements, to perform RT-based phylogenetic trees and to analyze their insertion times using High Performance Computing (HPC) techniques. Inpactor was tested on the classification and annotation of LTR retrotransposons in pineapple, a recently-sequenced genome. The pineapple genome assembly comprises 44% of transposable elements, of which 23% were classified as LTR retrotransposons. Exceptionally, 16.4% of the pineapple genome assembly corresponded to only one lineage of the Gypsy superfamily: Del, suggesting that this particular lineage has undergone a significant increase in its copy numbers. As demonstrated for the pineapple genome, Inpactor provides comprehensive data of LTR retrotransposons' classification and dynamics, allowing a fine understanding of their contribution to genome structure and evolution. Inpactor is available at https://github.com/simonorozcoarias/Inpactor.

Identification by the DArTseq method of the genetic origin of the Coffea canephora cultivated in Vietnam and Mexico.

BACKGROUND: The coffee species Coffea canephora is commercially identified as "Conilon" when produced in Brazil, or "Robusta" when produced elsewhere in the world. It represents approximately 40 % of coffee production worldwide. While the genetic diversity of wild C. canephora has been well studied in the past, only few studies have addressed the genetic diversity of currently cultivated varieties around the globe. Vietnam is the largest Robusta producer in the world, while Mexico is the only Latin American country, besides Brazil, that has a significant Robusta production. Knowledge of the genetic origin of Robusta cultivated varieties in countries as important as Vietnam and Mexico is therefore of high interest. RESULTS: Through the use of Sequencing-based diversity array technology-DArTseq method-on a collection of C. canephora composed of known accessions and accessions cultivated in Vietnam and Mexico, 4,021 polymorphic SNPs were identified. We used a multivariate analysis using SNP data from reference accessions in order to confirm and further fine-tune the genetic diversity of C. canephora. Also, by interpolating the data obtained for the varieties from Vietnam and Mexico, we determined that they are closely related to each other, and identified that their genetic origin is the Robusta Congo - Uganda group. CONCLUSIONS: The genetic characterization based on SNP markers of the varieties grown throughout the world, increased our knowledge on the genetic diversity of C. canephora, and contributed to the understanding of the genetic background of varieties from very important coffee producers. Given the common genetic origin of the Robusta varieties cultivated in Vietnam, Mexico and Uganda, and the similar characteristics of climatic areas and relatively high altitude where they are grown, we can state that the Vietnamese and the Mexican Robusta have the same genetic potential to produce good cup quality.

Major QTLs control resistance to rice hoja blanca virus and its vector Tagosodes orizicolus.

Rice hoja blanca (white leaf) disease can cause severe yield losses in rice in the Americas. The disease is caused by the rice hoja blanca virus (RHBV), which is transmitted by the planthopper vector Tagosodes orizicolus. Because classical breeding schemes for this disease rely on expensive, time-consuming screenings, there is a need for alternatives such as marker-aided selection. The varieties Fedearroz 2000 and Fedearroz 50, which are resistant to RHBV and to the feeding damage caused by T. orizicolus, were crossed with the susceptible line WC366 to produce segregating F2:3 populations. The F3 families were scored for their resistance level to RHBV and T. orizicolus. The F2:3 lines of both crosses were genotyped using microsatellite markers. One major QTL on the short arm of chromosome 4 was identified for resistance to RHBV in the two populations. Two major QTL on chromosomes 5 and 7 were identified for resistance to T. orizicolus in the Fd2000 × WC366 and Fd50 × WC366 crosses, respectively. This comparative study using two distinct rice populations allowed for a better understanding of how the resistance to RHBV and its vector are controlled genetically. Simple marker-aided breeding schemes based on QTL information can be designed to improve rice germplasm to reduce losses caused by this important disease.

Patterns of sequence divergence and evolution of the S orthologous regions between Asian and African cultivated rice species.

A strong postzygotic reproductive barrier separates the recently diverged Asian and African cultivated rice species, Oryza sativa and O. glaberrima. Recently a model of genetic incompatibilities between three adjacent loci: S(1)A, S(1) and S(1)B (called together the S(1) regions) interacting epistatically, was postulated to cause the allelic elimination of female gametes in interspecific hybrids. Two candidate factors for the S(1) locus (including a putative F-box gene) were proposed, but candidates for S(1)A and S(1)B remained undetermined. Here, to better understand the basis of the evolution of regions involved in reproductive isolation, we studied the genic and structural changes accumulated in the S(1) regions between orthologous sequences. First, we established an 813 kb genomic sequence in O. glaberrima, covering completely the S(1)A, S(1) and the majority of the S(1)B regions, and compared it with the orthologous regions of O. sativa. An overall strong structural conservation was observed, with the exception of three isolated regions of disturbed collinearity: (1) a local invasion of transposable elements around a putative F-box gene within S(1), (2) the multiple duplication and subsequent divergence of the same F-box gene within S(1)A, (3) an interspecific chromosomal inversion in S(1)B, which restricts recombination in our O. sativa×O. glaberrima crosses. Beside these few structural variations, a uniform conservative pattern of coding sequence divergence was found all along the S(1) regions. Hence, the S(1) regions have undergone no drastic variation in their recent divergence and evolution between O. sativa and O. glaberrima, suggesting that a small accumulation of genic changes, following a Bateson-Dobzhansky-Muller (BDM) model, might be involved in the establishment of the sterility barrier. In this context, genetic incompatibilities involving the duplicated F-box genes as putative candidates, and a possible strengthening step involving the chromosomal inversion might participate to the reproductive barrier between Asian and African rice species.

A genetic model for the female sterility barrier between Asian and African cultivated rice species.

S(1) is the most important locus acting as a reproductive barrier between Oryza sativa and O. glaberrima. It is a complex locus, with factors that may affect male and female fertility separately. Recently, the component causing the allelic elimination of pollen was fine mapped. However, the position and nature of the component causing female sterility remains unknown. To fine map the factor of the S(1) locus affecting female fertility, we developed a mapping approach based on the evaluation of the degree of female transmission ratio distortion (fTRD) of markers. Through implementing this methodology in four O. sativa x O. glaberrima crosses, the female component of the S(1) locus was mapped into a 27.8-kb (O. sativa) and 50.3-kb (O. glaberrima) region included within the interval bearing the male component of the locus. Moreover, evidence of additional factors interacting with S(1) was also found. In light of the available data, a model where incompatibilities in epistatic interactions between S(1) and the additional factors are the cause of the female sterility barrier between O. sativa and O. glaberrima was developed to explain the female sterility and the TRD mediated by S(1). According to our model, the recombination ratio and allelic combinations between these factors would determine the final allelic frequencies observed for a given cross.

A universal core genetic map for rice.

To facilitate the creation of easily comparable, low-resolution genetic maps with evenly distributed markers in rice (Oryza sativa L.), we conceived of and developed a Universal Core Genetic Map (UCGM). With this aim, we derived a set of 165 anchors, representing clusters of three microsatellite or simple sequence repeat (SSR) markers arranged into non-recombining groups. Each anchor consists of at least three, closely linked SSRs, located within a distance below the genetic resolution provided by common, segregating populations (<500 individuals). We chose anchors that were evenly distributed across the rice chromosomes, with spacing between 2 and 3.5 Mbp (except in the telomeric regions, where spacing was 1.5 Mbp). Anchor selection was performed using in silico tools and data: the O. sativa cv. Nipponbare rice genome sequence, the CHARM tool, information from the Gramene database and the OrygenesDB database. Sixteen AA-genome accessions of the Oryza genus were used to evaluate polymorphisms for the selected markers, including accessions from O. sativa, O. glaberrima, O. barthii, O. rufipogon, O. glumaepatula and O. meridionalis. High levels of polymorphism were found for the tested O. sativa x O. glaberrima or O. sativa x wild rice combinations. We developed Paddy Map, a simple database that is helpful in selecting optimal sets of polymorphic SSRs for any cross that involves the previously mentioned species. Validation of the UCGM was done by using it to develop three interspecific genetic maps and by comparing genetic SSR locations with their physical positions on the rice pseudomolecules. In this study, we demonstrate that the UCGM is a useful tool for the rice genetics and breeding community, especially in strategies based on interspecific hybridisation.