Summary: Flax (Linum usitatissimum) is an ancient crop that is widely cultivated as a source of fiber, oil and medicinally relevant compounds. To accelerate crop improvement, we performed whole‐genome shotgun sequencing of the nuclear genome of flax. Seven paired‐end libraries ranging in size from 300 bp to 10 kb were sequenced using an Illumina genome analyzer. A de novo assembly, comprised exclusively of deep‐coverage (approximately 94× raw, approximately 69× filtered) short‐sequence reads (44–100 bp), produced a set of scaffolds with N₅₀ = 694 kb, including contigs with N₅₀ = 20.1 kb. The contig assembly contained 302 Mb of non‐redundant sequence representing an estimated 81% genome coverage. Up to 96% of published flax ESTs aligned to the whole‐genome shotgun scaffolds. However, comparisons with independently sequenced BACs and fosmids showed some mis‐assembly of regions at the genome scale. A total of 43 384 protein‐coding genes were predicted in the whole‐genome shotgun assembly, and up to 93% of published flax ESTs, and 86% of A. thaliana genes aligned to these predicted genes, indicating excellent coverage and accuracy at the gene level. Analysis of the synonymous substitution rates (K_s) observed within duplicate gene pairs was consistent with a recent (5–9 MYA) whole‐genome duplication in flax. Within the predicted proteome, we observed enrichment of many conserved domains (Pfam‐A) that may contribute to the unique properties of this crop, including agglutinin proteins. Together these results show that de novo assembly, based solely on whole‐genome shotgun short‐sequence reads, is an efficient means of obtaining nearly complete genome sequence information for some plant species.

Fusarium oxysporum f. sp. lini is a hemibiotrophic fungus that causes wilt in flax. Along with rust, fusarium wilt has become an important factor in flax production worldwide. Resistant flax cultivars have been used to manage the disease, but the resistance varies, depending on the interactions between specific cultivars and isolates of the pathogen. This interaction has a strong molecular basis, but no genomic information is available on how the plant responds to attempted infection, to inform breeding programs on potential candidate genes to evaluate or improve resistance across cultivars. In the current study, disease progression in two flax cultivars (CDC Bethune and Lutea), showed earlier disease symptoms and higher susceptibility in the later cultivar. Chitinase gene expression was also divergent and demonstrated and earlier molecular response in Lutea. The most resistant cultivar (CDC Bethune) was used for a full RNA-seq transcriptome study through a time-course at 2, 4, 8 and 18 days post-inoculation (DPI). While over 100 genes were significantly differentially expressed at both 4 and 8 DPI, the broadest deployment of plant defense responses was evident at 18 DPI with transcripts of more than 1,000 genes responding to the treatment. These genes evidenced a reception and transduction of pathogen signals, a large transcriptional reprogramming, induction of hormone signalling, activation of pathogenesis-related (PR) genes, and changes in secondary metabolism. Among these several key genes, that consistently appear in studies of plant-pathogen interactions, had increased transcript abundance in our study, and constitute suitable candidates for resistance breeding programs. These included: an induced RPMI-induced protein kinase (RIPK); transcription factors WRKY3, WRKY70, WRKY75, MYB113 and MYB108; the ethylene response factors ERF1 and ERF14; two genes involved in auxin/glucosinolate precursor synthesis (CYP79B2 and CYP79B3); the flavonoid-related enzymes chalcone synthase, dihydroflavonol reductase and multiple anthocyanidin synthases; and a peroxidase implicated in lignin formation (PRX52). Additionally, regulation of some genes indicated potential pathogen manipulation to facilitate infection; these included four disease resistance proteins that were repressed, indole acetic acid amido/amino hydrolases which were upregulated, activated expansins and glucanases, amino acid transporters and aquaporins, and finally, repression of major latex proteins.

Clubroot, caused by the obligate parasite Plasmodiophora brassicae, is an important disease of the Brassicaceae and poses a significant threat to the $26.7 billion canola/oilseed rape (Brassica napus) industry in western Canada. While clubroot is managed most effectively by planting resistant host varieties, new pathotypes of P. brassicae have emerged recently that can overcome this resistance. Whole genome analyses provide both a toolbox and a systemic view of molecular mechanisms in host-pathogen interactions, which can be used to design new breeding strategies to increase P. brassicae resistance. We used RNA-seq to evaluate differential gene expression at 7, 14 and 21 days after inoculation (dai) of two B. napus genotypes with differential responses to P. brassicae pathotype 5X. Gall development was evident at 14 dai in the susceptible genotype (the oilseed rape ‘Brutor’), while gall development in the resistant genotype (the rutabaga (B. napus) ‘Laurentian’) was limited and not visible until 21 dai. Immune responses were better sustained through the time-course in ‘Laurentian’, and numerous genes from immune-related functional categories were associated with salicylic acid (SA)-mediated responses. Jasmonic acid (JA)-mediated responses seemed to be mostly inhibited, especially in the resistant genotype. The upregulation of standard defense-related proteins, like chitinases and thaumatins, was evident in ‘Laurentian’. The enrichment, in both host genotypes, of functional categories for syncytium formation and response to nematodes indicated that cell enlargement during P. brassicae infection, and the metabolic processes therein, share similarities with the response to infection by nematodes that produce similar anatomical symptoms. An analysis of shared genes between the two genotypes at different time-points, confirmed that the nematode-like responses occurred earlier for ‘Brutor’, along with cell metabolism and growth changes. Additionally, the susceptible cultivar turned off defense mechanisms earlier than ‘Laurentian’. Collectively, this study showed the importance of SA in triggering immune responses and suggested some key resistance and susceptibility factors that can be used in future studies for resistance breeding through gene-editing approaches.

Clubroot, caused by Plasmodiophora brassicae Woronin, is an important soilborne disease of Brassica napus L. and other crucifers. To improve understanding of the mechanisms of resistance and pathogenesis in the clubroot pathosystem, the rutabaga (B. napus subsp. rapifera Metzg) cultivars ‘Wilhelmsburger’ (resistant) and ‘Laurentian’ (susceptible) were inoculated with P. brassicae pathotype 3A and their transcriptomes were analyzed at 7, 14, and 21 days after inoculation (dai) by RNA sequencing (RNA-seq). Thousands of transcripts with significant changes in expression were identified in each host at each time-point in inoculated vs. non-inoculated plants. Molecular responses at 7 and 14 dai supported clear differences in the clubroot response mechanisms of the two genotypes. Both the resistant and the susceptible cultivars activated receptor-like protein (RLP) genes, resistance (R) genes, and genes involved in salicylic acid (SA) signaling as clubroot defense mechanisms. In addition, genes related to calcium signaling and genes encoding leucine-rich repeat (LRR) receptor kinases, the respiratory burst oxidase homolog (RBOH) protein, and transcription factors such as WRKYs, ethylene responsive factors, and basic leucine zippers (bZIPs), appeared to be upregulated in ‘Wilhelmsburger’ to restrict P. brassicae development. Some of these genes are essential components of molecular defenses, including ethylene (ET) signaling and the oxidative burst. Our study highlights the importance of activation of genes associated with SA- and ET-mediated responses in the resistant cultivar. A set of candidate genes showing contrasting patterns of expression between the resistant and susceptible cultivars was identified and includes potential targets for further study and validation through approaches such as gene editing.