A systems genomics and genetics approach to identify the genetic regulatory network for lignin content in Brassica napus seeds.

Seed quality traits of oilseed rape, Brassica napus (B. napus), exhibit quantitative inheritance determined by its genetic makeup and the environment via the mediation of a complex genetic architecture of hundreds to thousands of genes. Thus, instead of single gene analysis, network-based systems genomics and genetics approaches that combine genotype, phenotype, and molecular phenotypes offer a promising alternative to uncover this complex genetic architecture. In the current study, systems genetics approaches were used to explore the genetic regulation of lignin traits in B. napus seeds. Four QTL (qLignin_A09_1, qLignin_A09_2, qLignin_A09_3, and qLignin_C08) distributed on two chromosomes were identified for lignin content. The qLignin_A09_2 and qLignin_C08 loci were homologous QTL from the A and C subgenomes, respectively. Genome-wide gene regulatory network analysis identified eighty-three subnetworks (or modules); and three modules with 910 genes in total, were associated with lignin content, which was confirmed by network QTL analysis. eQTL (expression quantitative trait loci) analysis revealed four cis-eQTL genes including lignin and flavonoid pathway genes, cinnamoyl-CoA-reductase (CCR1), and TRANSPARENT TESTA genes TT4, TT6, TT8, as causal genes. The findings validated the power of systems genetics to identify causal regulatory networks and genes underlying complex traits. Moreover, this information may enable the research community to explore new breeding strategies, such as network selection or gene engineering, to rewire networks to develop climate resilience crops with better seed quality.

Practical Synthesis of 6-Amino-1-hydroxy-2,1-benzoxaborolane: A Key Intermediate of DNDI-6148.

Visceral leishmaniasis (VL), a parasitic, poverty-linked, neglected disease, is endemic across multiple regions of the world and fatal if untreated. There is an urgent need for a better and more affordable treatment for VL. DNDI-6148 is a promising drug candidate being evaluated for the treatment of VL; however, the current process for producing the key intermediate of DNDI-6148, 6-amino-1-hydroxy-2,1-benzoxaborolane, is expensive and difficult to scale up. Herein, we describe two practical approaches to synthesizing 6-amino-1-hydroxy-2,1-benzoxaborolane from inexpensive and readily available raw materials. Starting with 4-tolunitrile, the first approach is a five-step sequence involving a Hofmann rearrangement, resulting in an overall yield of 40%. The second approach utilizes 2-methyl-5-nitroaniline as the starting material and features borylation of aniline and continuous flow hydrogenation as the key steps, with an overall yield of 46%. Both routes bypass the nitration of 1-hydroxy-2,1-benzoxaborolane, which is challenging and expensive to scale. In particular, the second approach is more practical and scalable because of the mild operating conditions and facile isolation process.

Genomic asymmetry of the Brassica napus seed: epigenetic contributions of DNA methylation and small RNAs to subgenome bias.

Polyploidy is a persistent phenomenon in angiosperm genome evolution that is hypothesized to have contributed to the diversity of extant flowering plants. Brassica napus, one of the world's most important angiosperm oilseed species, originated from the interspecific hybridization of Brassica rapa (An ) and Brassica oleracea (Cn ). While the trends of genome dominance in transcriptomics are beginning to emerge, less is known about the epigenetic and small RNA landscapes in polyploids during reproductive development. The seed is the pivotal developmental transition into the new sporophytic generation, and experiences substantial epigenetic modifications over time. Here, we investigated the prevalence of bias in the contexts of DNA methylation and small interfering (si)RNA profiles in both subgenomes (An and Cn ), as well as the ancestral fractionated genomes across B. napus seed development. We report ubiquitous Cn subgenome bias of siRNA expression and cytosine methylation, with DNA methylation being particularly abundant on gene promoters in the Cn subgenome. Further, we provide evidence that siRNA transcriptional patterns were conserved within the ancestral triplicated subgenomes of B. napus, but not across the An and Cn subgenomes. We discuss how methylation patterns in the B. napus seed relate to genes, promoter regions, siRNA loci and transposable elements through the lens of genome fractionation and polyploidization. Taken together we provide evidence for epigenetic regulation selectively silencing the Cn subgenome during seed development, and explore the impact of genome fractionation on the epigenetic components of the B. napus seed.

Transcriptomic Analysis of Early Flowering Signals in 'Royal' Flax.

Canada is one of the world's leading producers and exporters of flax seed, with most production occurring in the Prairie Provinces. However, reduced season length and risk of frost restricts production in the northern grain belt of the Canadian Prairies. To expand the growing region of flax and increase production in Canada, flax breeders need to develop earlier-flowering varieties capable of avoiding the risk of abiotic stress. A thorough understanding of flowering control of flax is essential for the efficient breeding of such lines. We identified 722 putative flax flowering genes that span all major flowering-time pathways. Frequently, we found multiple flax homologues for a single Arabidopsis flowering gene. We used RNA sequencing to quantify the expression of genes in the shoot apical meristem (SAM) at 10, 15, 19, and 29 days after planting (dap) using the 'Royal' cultivar. We observed the expression of 80% of putative flax flowering genes and the differential expression of only 30%; these included homologues of major flowering regulators, such as SOC1, FUL, and AP1. We also found enrichment of differentially expressed genes (DEGs) in transcription factor (TF) families involved in flowering. Finally, we identified the candidates' novel flowering genes amongst the uncharacterized flax genes. Our transcriptomic dataset provides a useful resource for investigating the regulatory control of the transition to flowering in flax and for the breeding of northern-adapted varieties.

Gene expression profiling reveals transcription factor networks and subgenome bias during Brassica napus seed development.

We profiled the global gene expression landscape across the reproductive lifecycle of Brassica napus. Comparative analysis of this nascent amphidiploid revealed the contribution of each subgenome to plant reproduction. Whole-genome transcription factor networks identified BZIP11 as a transcriptional regulator of early B. napus seed development. Knockdown of BZIP11 using RNA interference resulted in a similar reduction in gene activity of predicted gene targets, and a reproductive-lethal phenotype. Global mRNA profiling revealed lower accumulation of Cn subgenome transcripts relative to the An subgenome. Subgenome-specific transcription factor networks identified distinct transcription factor families enriched in each of the An and Cn subgenomes early in seed development. Analysis of laser-microdissected seed subregions further reveal subgenome expression dynamics in the embryo, endosperm and seed coat of early stage seeds. Transcription factors predicted to be regulators encoded by the An subgenome are expressed primarily in the seed coat, whereas regulators encoded by the Cn subgenome were expressed primarily in the embryo. Data suggest subgenome bias are characteristic features of the B. napus seed throughout development, and that such bias might not be universal across the embryo, endosperm and seed coat of the developing seed. Transcriptional networks spanning both the An and Cn genomes of the whole B. napus seed can identify valuable targets for seed development research and that -omics level approaches to studying gene regulation in B. napus can benefit from both broad and high-resolution analyses.

3'-(Phenyl alkynyl) analogs of abscisic acid: synthesis and biological activity of potent ABA antagonists.

We report here the synthesis and biological testing of 3'-(phenyl alkynyl) abscisic ABA analogs, a new class of potent ABA antagonists. These ABA analogs incorporate a rigid framework of eight carbon atoms attached at the 3'-carbon atom of ABA that prevents folding of the ABA analog-bound receptor required for ABA signalling. The two-step synthesis is based upon the optimized conversion of natural (S)-ABA to 3'-iodo ABA which can be coupled to phenyl acetylenes using Sonogashira conditions, or to styryl compounds through Suzuki chemistry. The parent 3'-(phenyl alkynyl) ABA analog 7 was obtained in 29% yield, 74% yield based on recovered starting material. In a lentil seed germination assay, compound 7 was found to have more potent activity than other known 3'-substituted ABA antagonists to date. In a structure activity study parasubstituted phenyl alkynyl analogs had comparable activity to the analog 7 while the 3'-styryl ABA 18 was only slightly less active. Analog 7 overcame ABA inhibition of germination and seedling growth in a wide range of mono and dicot plant species, including canola, lentil, soybean, rice, wheat, barley, cannabis and canary seed. 3'-(Phenyl alkynyl) ABA analogs have numerous potential practical agricultural applications including promoting ripening of crops, dormancy breaking of seeds and woody perennials, as well as promoting seed germination, and growth under stress conditions as demonstrated in this report.

Machine learning analyses of methylation profiles uncovers tissue-specific gene expression patterns in wheat.

DNA methylation is a mechanism of epigenetic modification in eukaryotic organisms. Generally, methylation within genes promoter inhibits regulatory protein binding and represses transcription, whereas gene body methylation is associated with actively transcribed genes. However, it remains unclear whether there is interaction between methylation levels across genic regions and which site has the biggest impact on gene regulation. We investigated and used the methylation patterns of the bread wheat cultivar Chinese Spring to uncover differentially expressed genes (DEGs) between roots and leaves, using six machine learning algorithms and a deep neural network. As anticipated, genes with higher expression in leaves were mainly involved in photosynthesis and pigment biosynthesis processes whereas genes that were not differentially expressed between roots and leaves were involved in protein processes and membrane structures. Methylation occurred preponderantly (60%) in the CG context, whereas 35 and 5% of methylation occurred in CHG and CHH contexts, respectively. Methylation levels were highly correlated (r = 0.7 to 0.9) between all genic regions, except within the promoter (r = 0.4 to 0.5). Machine learning models gave a high (0.81) prediction accuracy of DEGs. There was a strong correlation (p-value = 9.20×10-10 ) between all features and gene expression, suggesting that methylation across all genic regions contribute to gene regulation. However, the methylation of the promoter, the CDS and the exon in CG context was the most impactful. Our study provides more insights into the interplay between DNA methylation and gene expression and paves the way for identifying tissue-specific genes using methylation profiles.

Characterization of B-Genome Specific High Copy hAT MITE Families in Brassica nigra Genome.

Miniature inverted-repeat transposable elements (MITEs) are non-autonomous class II transposons which have been shown to influence genome evolution. Brassica nigra L. (B-genome) is one of three Brassica diploids cultivated primarily as an oil crop, which harbors novel alleles important for breeding. Two new high copy hAT MITE families (BniHAT-1 and BniHAT-2) from the B-genome were characterized and their prevalence assessed in the genomes of the related diploids, rapa L. (A) and Brassica oleracea L. (C). Both novel MITE families were present at high copy numbers in the B-genome with 434 and 331 copies of BniHAT-1 and BniHAT-2, respectively. Yet less than 20 elements were identified in the genome assemblies of the A, and C -genomes, supporting B-genome specific proliferation of these MITE families. Although apparently randomly distributed across the genome, 68 and 70% of the B-genome MITEs were present within 2 kb flanking regions of annotated genes suggesting they might influence gene expression and/or function. In addition, MITE derived microRNAs and transcription factor binding sites suggested a putative role in gene regulation. Age of insertion analysis revealed that the major proliferation of these elements occurred during 2-3 million years ago. Additionally, site-specific polymorphism analyses showed that 44% MITEs were undergoing active amplification into the B-genome. Overall, this study provides a comprehensive analysis of two high copy MITE families, which were specifically amplified in the B-genome, suggesting a potential role in shaping the Brassica B-genome.

A high-contiguity Brassica nigra genome localizes active centromeres and defines the ancestral Brassica genome.

It is only recently, with the advent of long-read sequencing technologies, that we are beginning to uncover previously uncharted regions of complex and inherently recursive plant genomes. To comprehensively study and exploit the genome of the neglected oilseed Brassica nigra, we generated two high-quality nanopore de novo genome assemblies. The N50 contig lengths for the two assemblies were 17.1 Mb (12 contigs), one of the best among 324 sequenced plant genomes, and 0.29 Mb (424 contigs), respectively, reflecting recent improvements in the technology. Comparison with a de novo short-read assembly corroborated genome integrity and quantified sequence-related error rates (0.2%). The contiguity and coverage allowed unprecedented access to low-complexity regions of the genome. Pericentromeric regions and coincidence of hypomethylation enabled localization of active centromeres and identified centromere-associated ALE family retro-elements that appear to have proliferated through relatively recent nested transposition events (<1 Ma). Genomic distances calculated based on synteny relationships were used to define a post-triplication Brassica-specific ancestral genome, and to calculate the extensive rearrangements that define the evolutionary distance separating B. nigra from its diploid relatives.

Core and Differentially Abundant Bacterial Taxa in the Rhizosphere of Field Grown Brassica napus Genotypes: Implications for Canola Breeding.

Modifying the rhizosphere microbiome through targeted plant breeding is key to harnessing positive plant-microbial interrelationships in cropping agroecosystems. Here, we examine the composition of rhizosphere bacterial communities of diverse Brassica napus genotypes to identify: (1) taxa that preferentially associate with genotypes, (2) core bacterial microbiota associated with B. napus, (3) heritable alpha diversity measures at flowering and whole growing season, and (4) correlation between microbial and plant genetic distance among canola genotypes at different growth stages. Our aim is to identify and describe signature microbiota with potential positive benefits that could be integrated in B. napus breeding and management strategies. Rhizosphere soils of 16 diverse genotypes sampled weekly over a 10-week period at single location as well as at three time points at two additional locations were analyzed using 16S rRNA gene amplicon sequencing. The B. napus rhizosphere microbiome was characterized by diverse bacterial communities with 32 named bacterial phyla. The most abundant phyla were Proteobacteria, Actinobacteria, and Acidobacteria. Overall microbial and plant genetic distances were highly correlated (R = 0.65). Alpha diversity heritability estimates were between 0.16 and 0.41 when evaluated across growth stage and between 0.24 and 0.59 at flowering. Compared with a reference B. napus genotype, a total of 81 genera were significantly more abundant and 71 were significantly less abundant in at least one B. napus genotype out of the total 558 bacterial genera. Most differentially abundant genera were Proteobacteria and Actinobacteria followed by Bacteroidetes and Firmicutes. Here, we also show that B. napus genotypes select an overall core bacterial microbiome with growth-stage-related patterns as to how taxa joined the core membership. In addition, we report that sets of B. napus core taxa were consistent across our three sites and 2 years. Both differential abundance and core analysis implicate numerous bacteria that have been reported to have beneficial effects on plant growth including disease suppression, antifungal properties, and plant growth promotion. Using a multi-site year, temporally intensive field sampling approach, we showed that small plant genetic differences cause predictable changes in canola microbiome and are potential target for direct and indirect selection within breeding programs.

Structural features of two major nucleolar organizer regions (NORs), Nor-B1 and Nor-B2, and chromosome-specific rRNA gene expression in wheat.

The reference genome sequence of wheat 'Chinese Spring' (CS) is now available (IWGSC RefSeq v1.0), but the core sequences defining the nucleolar organizer regions (NORs) have not been characterized. We estimated that the total copy number of the rDNA units in the wheat genome is 11 160, of which 30.5%, 60.9% and 8.6% are located on Nor-B1 (1B), Nor-B2 (6B) and other NORs, respectively. The total length of the NORs is estimated to be 100 Mb, corresponding to approximately 10% of the unassembled portion of the genome not represented in RefSeq v1.0. Four subtypes (S1-S4) of the rDNA units were identified based on differences within the 3' external transcribed spacer regions in Nor-B1 and Nor-B2, and quantitative PCR indicated locus-specific variation in rDNA subtype contents. Expression analyses of rDNA subtypes revealed that S1 was predominantly expressed and S2 weakly expressed, in contrast to the relative abundance of rDNA subtypes in the wheat genome. These results suggest a regulation mechanism of differential rDNA expression based on sequence differences. S3 expression increased in the ditelosomic lines Dt1BL and Dt6BL, suggesting that S3 is subjected to chromosome-mediated silencing. Structural differences were detected in the regions surrounding the NOR among homoeologous chromosomes of groups 1 and 6. The adjacent regions distal to the major NORs were expanded compared with their homoeologous counterparts, and the gene density of these expanded regions was relatively low. We provide evidence that these regions are likely to be important for autoregulation of the associated major NORs as well as silencing of minor NORs.

Changes in gene expression in Camelina sativa roots and vegetative tissues in response to salinity stress.

The response of Camelina sativa to salt stress was examined. Salt reduced shoot, but not root length. Root and shoot weight were affected by salt, as was photosynthetic capacity. Salt did not alter micro-element concentration in shoots, but increased macro-element (Ca and Mg) levels. Gene expression patterns in shoots indicated that salt stress may have led to shuttling of Na+ from the cytoplasm to the tonoplast and to an increase in K+ and Ca+2 import into the cytoplasm. In roots, gene expression patterns indicated that Na+ was exported from the cytoplasm by the SOS pathway and that K+ was imported in response to salt. Genes involved in chelation and storage were up-regulated in shoots, while metal detoxification appeared to involve various export mechanisms in roots. In shoots, genes involved in secondary metabolism leading to lignin, anthocyanin and wax production were up-regulated. Partial genome partitioning was observed in roots and shoots based on the expression of homeologous genes from the three C. sativa sub-genomes. Sub-genome I and II were involved in the response to salinity stress to about the same degree, while about 10% more differentially-expressed genes were associated with sub-genome III.

Changes in the Sclerotinia sclerotiorum transcriptome during infection of Brassica napus.

BACKGROUND: Sclerotinia sclerotiorum causes stem rot in Brassica napus, which leads to lodging and severe yield losses. Although recent studies have explored significant progress in the characterization of individual S. sclerotiorum pathogenicity factors, a gap exists in profiling gene expression throughout the course of S. sclerotiorum infection on a host plant. In this study, RNA-Seq analysis was performed with focus on the events occurring through the early (1 h) to the middle (48 h) stages of infection. RESULTS: Transcript analysis revealed the temporal pattern and amplitude of the deployment of genes associated with aspects of pathogenicity or virulence during the course of S. sclerotiorum infection on Brassica napus. These genes were categorized into eight functional groups: hydrolytic enzymes, secondary metabolites, detoxification, signaling, development, secreted effectors, oxalic acid and reactive oxygen species production. The induction patterns of nearly all of these genes agreed with their predicted functions. Principal component analysis delineated gene expression patterns that signified transitions between pathogenic phases, namely host penetration, ramification and necrotic stages, and provided evidence for the occurrence of a brief biotrophic phase soon after host penetration. CONCLUSIONS: The current observations support the notion that S. sclerotiorum deploys an array of factors and complex strategies to facilitate host colonization and mitigate host defenses. This investigation provides a broad overview of the sequential expression of virulence/pathogenicity-associated genes during infection of B. napus by S. sclerotiorum and provides information for further characterization of genes involved in the S. sclerotiorum-host plant interactions.

Stable Epigenetic Variants Selected from an Induced Hypomethylated Fragaria vesca Population.

Epigenetic inheritance was transmitted through selection over five generations of extreme early, but not late flowering time phenotypic lines in Fragaria vesca. Epigenetic variation was initially artificially induced using the DNA demethylation reagent 5-azacytidine (5-azaC). It is the first report to explore epigenetic variant selection and phenotypic trait inheritance in strawberry. Transmission frequency of these traits was determined across generations. The early flowering (EF4) and late stolon (LS) phenotypic traits were successfully transmitted across five and three generations through meiosis, respectively. Stable mitotic transmission of the early flowering phenotype was also demonstrated using clonal daughters derived from the 4th Generation (S4) mother plant. In order to further explore the DNA methylation patterns underlying the early flowering trait, the standard MSAP method using isoschizomers Hpa II/Msp I, and newly modified MSAP method using isoschizomers Tfi I/Pfe I which detected DNA methylation at CG, CHG, CHH sites were used in two early flowering lines, EF lines 1 (P2) and EF lines 2 (P3), and control lines (P1). A significant reduction in the number of fully-methylated bands was detected in P2 and P3 when compared to P1 using the novel MSAP method. In the standard MSAP, the symmetric CG and CHG methylation was maintained over generations in the early flowering lines based on the clustering in P2 and P3, the novel MSAP approach revealed the asymmetric CHH methylation pattern was not maintained over generations. This study provides evidence of stable selection of phenotypic traits, particularly early flowering through both meiosis and mitosis, which is meaningful to both breeding programs and commercial horticulture. The maintenance in CG and CHG methylation over generations suggests the early flowering phenotype might be related to DNA methylation alterations at the CG or CHG sites. Finally, this work provides a new approach for studying the role of epigenetics on complex quantitative trait improvement in strawberry, as well as providing a tool to expand phenotypic diversity and expedite potential new horticulture cultivar releases through either seed or vegetative propagation.

Quantitative trait variation is revealed in a novel hypomethylated population of woodland strawberry (Fragaria vesca).

BACKGROUND: Phenotypic variation is determined by a combination of genotype, environment and their interactions. The realization that allelic diversity can be both genetic and epigenetic allows the environmental component to be further separated. Partitioning phenotypic variation observed among inbred lines with an altered epigenome can allow the epigenetic component controlling quantitative traits to be estimated. To assess the contribution of epialleles on phenotypic variation and determine the fidelity with which epialleles are inherited, we have developed a novel hypomethylated population of strawberry (2n = 2x = 14) using 5-azacytidine from which individuals with altered phenotypes can be identified, selected and characterized. RESULTS: The hypomethylated population was generated using an inbred strawberry population in the F. vesca ssp. vesca accession Hawaii 4. Analysis of whole genome sequence data from control and hypomethylated lines indicate that 5-azacytidine exposure does not increase SNP above background levels. The populations contained only Hawaii 4 alleles, removing introgression of alternate F. vesca alleles as a potential source of variation. Although genome sequencing and genetic marker data are unable to rule out 5-azacytidine induced chromosomal rearrangements as a potential source of the trait variation observed, none were detected in our survey. Quantitative trait variation focusing on flowering time and rosette diameter was scored in control and treated populations where expanded levels of variation were observed among the hypomethylated lines. Methylation sensitive molecular markers indicated that 5-azacytidine induced alterations in DNA methylation patterns and inheritance of methylation patterns were confirmed by bisulfite sequencing of targeted regions. It is possible that methylation polymorphisms might underlie or have induced genetic changes underlying the observable differences in quantitative phenotypes. CONCLUSIONS: This population developed in a uniform genetic background provides a resource for the discovery of new variation controlling quantitative traits. Genome sequence analysis indicates that 5-azacytidine did not induce point mutations and the induced variation is largely restricted to DNA methylation. Using this resource, we have identified new variation and demonstrated the inheritance of both variant trait and methylation patterns. Although direct associations remain to be determined, these data suggest epigenetic variation might be subject to selection.

The developmental transcriptome atlas of the biofuel crop Camelina sativa.

Camelina sativa is currently being embraced as a viable industrial bio-platform crop due to a number of desirable agronomic attributes and the unique fatty acid profile of the seed oil that has applications for food, feed and biofuel. The recent completion of the reference genome sequence of C. sativa identified a young hexaploid genome. To complement this work, we have generated a genome-wide developmental transcriptome map by RNA sequencing of 12 different tissues covering major developmental stages during the life cycle of C. sativa. We have generated a digital atlas of this comprehensive transcriptome resource that enables interactive visualization of expression data through a searchable database of electronic fluorescent pictographs (eFP browser). An analysis of this dataset supported expression of 88% of the annotated genes in C. sativa and provided a global overview of the complex architecture of temporal and spatial gene expression patterns active during development. Conventional differential gene expression analysis combined with weighted gene expression network analysis uncovered similarities as well as differences in gene expression patterns between different tissues and identified tissue-specific genes and network modules. A high-quality census of transcription factors, analysis of alternative splicing and tissue-specific genome dominance provided insight into the transcriptional dynamics and sub-genome interplay among the well-preserved triplicated repertoire of homeologous loci. The comprehensive transcriptome atlas in combination with the reference genome sequence provides a powerful resource for genomics research which can be leveraged to identify functional associations between genes and understand the regulatory networks underlying developmental processes.

Multi-environment QTL studies suggest a role for cysteine-rich protein kinase genes in quantitative resistance to blackleg disease in Brassica napus.

BACKGROUND: Resistance to the blackleg disease of Brassica napus (canola/oilseed rape), caused by the hemibiotrophic fungal pathogen Leptosphaeria maculans, is determined by both race-specific resistance (R) genes and quantitative resistance loci (QTL), or adult-plant resistance (APR). While the introgression of R genes into breeding material is relatively simple, QTL are often detected sporadically, making them harder to capture in breeding programs. For the effective deployment of APR in crop varieties, resistance QTL need to have a reliable influence on phenotype in multiple environments and be well defined genetically to enable marker-assisted selection (MAS). RESULTS: Doubled-haploid populations produced from the susceptible B. napus variety Topas and APR varieties AG-Castle and AV-Sapphire were analysed for resistance to blackleg in two locations over 3 and 4 years, respectively. Three stable QTL were detected in each population, with two loci appearing to be common to both APR varieties. Physical delineation of three QTL regions was sufficient to identify candidate defense-related genes, including a cluster of cysteine-rich receptor-like kinases contained within a 49 gene QTL interval on chromosome A01. Individual L. maculans isolates were used to define the physical intervals for the race-specific R genes Rlm3 and Rlm4 and to identify QTL common to both field studies and the cotyledon resistance response. CONCLUSION: Through multi-environment QTL analysis we have identified and delineated four significant and stable QTL suitable for MAS of quantitative blackleg resistance in B. napus, and identified candidate genes which potentially play a role in quantitative defense responses to L. maculans.

The compact genome of the plant pathogen Plasmodiophora brassicae is adapted to intracellular interactions with host Brassica spp.

BACKGROUND: The protist Plasmodiophora brassicae is a soil-borne pathogen of cruciferous species and the causal agent of clubroot disease of Brassicas including agriculturally important crops such as canola/rapeseed (Brassica napus). P. brassicae has remained an enigmatic plant pathogen and is a rare example of an obligate biotroph that resides entirely inside the host plant cell. The pathogen is the cause of severe yield losses and can render infested fields unsuitable for Brassica crop growth due to the persistence of resting spores in the soil for up to 20 years. RESULTS: To provide insight into the biology of the pathogen and its interaction with its primary host B. napus, we produced a draft genome of P. brassicae pathotypes 3 and 6 (Pb3 and Pb6) that differ in their host range. Pb3 is highly virulent on B. napus (but also infects other Brassica species) while Pb6 infects only vegetable Brassica crops. Both the Pb3 and Pb6 genomes are highly compact, each with a total size of 24.2 Mb, and contain less than 2 % repetitive DNA. Clustering of genome-wide single nucleotide polymorphisms (SNP) of Pb3, Pb6 and three additional re-sequenced pathotypes (Pb2, Pb5 and Pb8) shows a high degree of correlation of cluster grouping with host range. The Pb3 genome features significant reduction of intergenic space with multiple examples of overlapping untranslated regions (UTRs). Dependency on the host for essential nutrients is evident from the loss of genes for the biosynthesis of thiamine and some amino acids and the presence of a wide range of transport proteins, including some unique to P. brassicae. The annotated genes of Pb3 include those with a potential role in the regulation of the plant growth hormones cytokinin and auxin. The expression profile of Pb3 genes, including putative effectors, during infection and their potential role in manipulation of host defence is discussed. CONCLUSION: The P. brassicae genome sequence reveals a compact genome, a dependency of the pathogen on its host for some essential nutrients and a potential role in the regulation of host plant cytokinin and auxin. Genome annotation supported by RNA sequencing reveals significant reduction in intergenic space which, in addition to low repeat content, has likely contributed to the P. brassicae compact genome.

Polyploid evolution of the Brassicaceae during the Cenozoic era.

The Brassicaceae (Cruciferae) family, owing to its remarkable species, genetic, and physiological diversity as well as its significant economic potential, has become a model for polyploidy and evolutionary studies. Utilizing extensive transcriptome pyrosequencing of diverse taxa, we established a resolved phylogeny of a subset of crucifer species. We elucidated the frequency, age, and phylogenetic position of polyploidy and lineage separation events that have marked the evolutionary history of the Brassicaceae. Besides the well-known ancient α (47 million years ago [Mya]) and ß (124 Mya) paleopolyploidy events, several species were shown to have undergone a further more recent (∼7 to 12 Mya) round of genome multiplication. We identified eight whole-genome duplications corresponding to at least five independent neo/mesopolyploidy events. Although the Brassicaceae family evolved from other eudicots at the beginning of the Cenozoic era of the Earth (60 Mya), major diversification occurred only during the Neogene period (0 to 23 Mya). Remarkably, the widespread species divergence, major polyploidy, and lineage separation events during Brassicaceae evolution are clustered in time around epoch transitions characterized by prolonged unstable climatic conditions. The synchronized diversification of Brassicaceae species suggests that polyploid events may have conferred higher adaptability and increased tolerance toward the drastically changing global environment, thus facilitating species radiation.

Transcriptome and methylome profiling reveals relics of genome dominance in the mesopolyploid Brassica oleracea.

BACKGROUND: Brassica oleracea is a valuable vegetable species that has contributed to human health and nutrition for hundreds of years and comprises multiple distinct cultivar groups with diverse morphological and phytochemical attributes. In addition to this phenotypic wealth, B. oleracea offers unique insights into polyploid evolution, as it results from multiple ancestral polyploidy events and a final Brassiceae-specific triplication event. Further, B. oleracea represents one of the diploid genomes that formed the economically important allopolyploid oilseed, Brassica napus. A deeper understanding of B. oleracea genome architecture provides a foundation for crop improvement strategies throughout the Brassica genus. RESULTS: We generate an assembly representing 75% of the predicted B. oleracea genome using a hybrid Illumina/Roche 454 approach. Two dense genetic maps are generated to anchor almost 92% of the assembled scaffolds to nine pseudo-chromosomes. Over 50,000 genes are annotated and 40% of the genome predicted to be repetitive, thus contributing to the increased genome size of B. oleracea compared to its close relative B. rapa. A snapshot of both the leaf transcriptome and methylome allows comparisons to be made across the triplicated sub-genomes, which resulted from the most recent Brassiceae-specific polyploidy event. CONCLUSIONS: Differential expression of the triplicated syntelogs and cytosine methylation levels across the sub-genomes suggest residual marks of the genome dominance that led to the current genome architecture. Although cytosine methylation does not correlate with individual gene dominance, the independent methylation patterns of triplicated copies suggest epigenetic mechanisms play a role in the functional diversification of duplicate genes.