The domestication of wild emmer wheat led to the selection of modern durum wheat, grown mainly for pasta production. We describe the 10.45 gigabase (Gb) assembly of the genome of durum wheat cultivar Svevo. The assembly enabled genome-wide genetic diversity analyses revealing the changes imposed by thousands of years of empirical selection and breeding. Regions exhibiting strong signatures of genetic divergence associated with domestication and breeding were widespread in the genome with several major diversity losses in the pericentromeric regions. A locus on chromosome 5B carries a gene encoding a metal transporter (TdHMA3-B1) with a non-functional variant causing high accumulation of cadmium in grain. The high-cadmium allele, widespread among durum cultivars but undetected in wild emmer accessions, increased in frequency from domesticated emmer to modern durum wheat. The rapid cloning of TdHMA3-B1 rescues a wild beneficial allele and demonstrates the practical use of the Svevo genome for wheat improvement.

Fusarium fujikuroi is the causal agent of bakanae, the most significant seed-borne disease of rice. Molecular mechanisms regulating defence responses of rice towards this fungus are not yet fully known. To identify transcriptional mechanisms underpinning rice resistance, a RNA-seq comparative transcriptome profiling was conducted on infected seedlings of selected rice genotypes at one and three weeks post germination (wpg). Twelve rice genotypes were screened against bakanae disease leading to the identification of Selenio and Dorella as the most resistant and susceptible cultivars, respectively. Transcriptional changes were more appreciable at 3 wpg, suggesting that this infection stage is essential to study the resistance mechanisms: 3,119 DEGs were found in Selenio and 5,095 in Dorella. PR1, germin-like proteins, glycoside hydrolases, MAP kinases, and WRKY transcriptional factors were up-regulated in the resistant genotype upon infection with F. fujikuroi. Up-regulation of chitinases and down-regulation of MAP kinases and WRKY transcriptional factors were observed in the susceptible genotype. Gene ontology (GO) enrichment analyses detected in Selenio GO terms specific to response to F. fujikuroi: ‘response to chitin’, ‘jasmonic acid biosynthetic process’, and ‘plant-type hypersensitive response’, while Dorella activated different mechanisms, such as ‘response to salicylic acid stimulus’ and ‘gibberellin metabolic process’, which was in agreement with the production of gibberellin A3 in Dorella plants. RNA-seq profiling was performed for the first time to analyse response of rice to F. fujikuroi infection. Our findings allowed the identification of genes activated in one- and three- week-old rice seedlings of two genotypes infected with F. fujikuroi. Furthermore, we found the pathways involved in bakanae resistance, such as response to chitin, JA-dependent signalling and hypersensitive response. Collectively, this provides important information to elucidate the molecular and cellular processes occurring in rice during F. fujikuroi infection and to develop bakanae resistant rice germplasm.

A pan-transcriptome describes the transcriptional and post-transcriptional consequences of genome diversity from multiple individuals within a species. We developed a barley pan-transcriptome using 20 inbred genotypes representing domesticated barley diversity by generating and analyzing short- and long-read RNA-sequencing datasets from multiple tissues. To overcome single reference bias in transcript quantification, we constructed genotype-specific reference transcript datasets (RTDs) and integrated these into a linear pan-genome framework to create a pan-RTD, allowing transcript categorization as core, shell or cloud. Focusing on the core (expressed in all genotypes), we observed significant transcript abundance variation among tissues and between genotypes driven partly by RNA processing, gene copy number, structural rearrangements and conservation of promotor motifs. Network analyses revealed conserved co-expression module::tissue correlations and frequent functional diversification. To complement the pan-transcriptome, we constructed a comprehensive cultivar (cv.) Morex gene-expression atlas and illustrate how these combined datasets can be used to guide biological inquiry.

Flavescence dorée is the most serious grapevine yellows disease in Europe. It is caused by phytoplasmas which are transmitted from grapevine to grapevine by the leafhopper Scaphoideus titanus. Differences in susceptibility among grapevine varieties suggest the existence of specific genetic features associated with resistance to the phytoplasma and/or possibly with its vector. In this work, RNA-Seq was used to compare early transcriptional changes occurring during the three-trophic interaction between the phytoplasma, its vector and the grapevine, represented by two different cultivars, one very susceptible to the disease and the other scarcely susceptible. The comparative analysis of the constitutive transcriptomic profiles suggests the existence of passive defense strategies against the insect and/or the phytoplasma in the scarcely-susceptible cultivar. Moreover, the attack by the infective vector on the scarcely-susceptible variety prompted immediate and substantial transcriptomic changes that led to the rapid erection of further active defenses. On the other hand, in the most susceptible variety the response was delayed and mainly consisted of the induction of phytoalexin synthesis. Surprisingly, the jasmonic acid- and ethylene-mediated defense reactions, activated by the susceptible cultivar following FD-free insect feeding, were not detected in the presence of the phytoplasma-infected vector. The comparison of the transcriptomic response in two grapevine varieties with different levels of susceptibility to Flavescence dorèe highlighted both passive and active defense mechanisms against the vector and/or the pathogen in the scarcely-susceptible variety, as well as the capacity of the phytoplasmas to repress the defense reaction against the insect in the susceptible variety.

Eggplant (Solanum melongena L.), similar to many other crops, suffers from soil-borne diseases, including Fusarium oxysporum f. sp. melongenae (Fom), causing wilting and heavy yield loss. To date, the genetic factors underlying plant responses to Fom are not well known. We previously developed a Recombinant Inbred Lines (RILs) population using as a female parent the fully resistant line ‘305E40’ and as a male parent the partially resistant line ‘67/3’. The fully resistant trait to Fom was introgressed from the allied species S. aethiopicum. In this work, the RIL population was assessed for the responses to Fom and by using a genomic mapping approach, two major QTLs on chromosomes CH02 and CH11 were identified, associated with the full and partial resistance trait to Fom, respectively. A targeted BSAseq procedure in which Illumina reads bulks of RILs grouped according to their resistance score was aligned to the appropriate reference genomes highlighted differentially enriched regions between resistant/susceptible progeny in the genomic regions underlying both QTLs. The characterization of such regions allowed us to identify the most reliable candidate genes for the two resistance traits. With the aim of revealing exclusive species-specific contigs and scaffolds inherited from the allied species and thus associated with the full resistance trait, a draft de-novo assembly of available Illumina sequences of the ‘305E40’ parent was developed to better resolve the non-recombining genomic region on its CH02 carrying the introgressed Fom resistance locus from S. aethiopicum.

Fusarium head blight (FHB) is a problem of great concern in small grain cereals, especially wheat. A quantitative trait locus (QTL) for FHB resistance (FHB_SFI) located on the long arm of chromosome 2D in the spring wheat genotype Wuhan 1 is a resistance locus which has potential to improve the FHB resistance of bread wheat since it confers effective resistance to wheat breeding lines. Recently, differentially expressed genes (DEG) have been identified by comparing near isogenic lines (NIL) carrying the susceptible and resistant alleles for the 2DL QTL, using RNA-Seq. In the present study, we aimed to identify candidate genes located within the genetic interval for the 2DL QTL for FHB resistance, as assessed by single floret inoculation (FHB_SFI), and possibly contributing to it. Combining previous and additional bioinformatics analyses, 26 DEG that were located on chromosome arm 2DL were selected for further characterization of their expression profile by RT-qPCR. Seven of those DEG showed a consistent differential expression profile between either three pairs of near isogenic lines or other genotypes carrying the R and S alleles for the 2DL QTL for FHB resistance. UN25696, which was identified in previous expression work using microarray was also confirmed to have a differential expression pattern. Those eight candidate genes were further characterized in 85 lines of a double haploid mapping population derived from the cross Wuhan 1/Nyubai, the population where the 2DL QTL was originally identified. The expression QTL for gene Traes_2DL_179570792 overlapped completely with the mapping interval for the 2DL QTL for FHB_SFI while the expression QTL for UN25696 mapped near the QTL, but did not overlap with it. None of the other genes had a significant eQTL on chromosome 2DL. Higher expression of Traes_2DL_179570792 and UN25696 was associated with the resistant allele at that locus. Of the 26 DEG from the 2DL chromosome further characterized in this study, only two had an expression QTL located in or near the interval for the 2DL QTL. Traes_2DL_179570792 is the first expression marker identified as associated with the 2DL QTL.