Guy C. Barker

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

BackgroundBrassica 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.ResultsWe 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.ConclusionsDifferential 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.

Mathematical modelling and the control of lymphatic filariasis

The current global initiative to eliminate lymphatic filariasis represents one of the largest mass drug administration programmes ever conceived for the control of a parasitic disease. Yet, it is still not known whether the WHO-recommended primary strategy of applying annual single-dose mass chemotherapy with a combination of two drugs for 4-6 years will effectively break parasite transmission from all endemic communities. Here we review recent work on the development and application of a deterministic mathematical model of filariasis transmission, to show how models of parasite transmission will help resolve the key currently debated questions regarding the ultimate effectiveness of the global strategy to control filariasis. These critical questions include the required duration of mass treatment in different endemic areas, the optimal drug coverage required to meet control targets within prescribed timeframes, the impact and importance of adding vector control to mass chemotherapy regimens, and the likelihood of the development of drug resistance by treated worm populations. The results demonstrate the vital role that integrating these models into control programming can have in providing effective decision-support frameworks for undertaking the optimal design and monitoring of regional and global filariasis-control programmes. Operationally, the models show that the effectiveness of the strategy to achieve filariasis control will be determined by successfully addressing two key factors: the need to maintain high community treatment coverages, and the need to include vector control measures especially in areas of high endemicity.

Pairing and recombination at meiosis of Brassica rapa (AA) × Brassica napus (AACC) hybrids

Interspecific crosses contribute significantly to plant evolution enabling gene exchanges between species. The efficiency of interspecific crosses depends on the similarity between the implicated genomes as high levels of genome similarity are required to ensure appropriate chromosome pairing and genetic recombination. Brassica napus (AACC) is an allopolyploid, resulting from natural hybridization between Brassica rapa (AA) and Brassica oleracea (CC), both being diploid species derived from a common ancestor. To study the relationships between genomes of these Brassica species, we have determined simultaneously the pairing and recombination pattern of A and C chromosomes during meiosis of AAC triploid hybrids, which result from the interspecific cross between natural B. napus and B. rapa. Different AAC triploid hybrids and their progenies have been analysed using cytogenetic, BAC-FISH, and molecular techniques. In 71% of the pollen mother cells, homologous A chromosomes paired regularly, and usually one chromosome of each pair was transmitted to the progeny. C chromosomes remained mainly univalent, but were involved in homoeologous pairing in 21.5% of the cells, and 13% of the transmitted C chromosomes were either recombined or broken. The rate of transmission of C chromosomes depended on the identity of the particular chromosome and on the way the hybrid was crossed, as the male or as the female parent, to B. napus or to B. rapa. Gene transfers in triploid hybrids are favoured between A genomes of B. rapa and B. napus, but also occur between A and C genomes though at lower rates.

Characterization of metabolite quantitative trait loci and metabolic networks that control glucosinolate concentration in the seeds and leaves of Brassica napus

• Glucosinolates are a major class of secondary metabolites found in the Brassicaceae, whose degradation products are proving to be increasingly important for human health and in crop protection. • The genetic and metabolic basis of glucosinolate accumulation was dissected through analysis of total glucosinolate concentration and its individual components in both leaves and seeds of a doubled-haploid (DH) mapping population of oilseed rape/canola (Brassica napus). • The quantitative trait loci (QTL) that had an effect on glucosinolate concentration in either or both of the organs were integrated, resulting in 105 metabolite QTL (mQTL). Pairwise correlations between individual glucosinolates and prior knowledge of the metabolic pathways involved in the biosynthesis of different glucosinolates allowed us to predict the function of genes underlying the mQTL. Moreover, this information allowed us to construct an advanced metabolic network and associated epistatic interactions responsible for the glucosinolate composition in both leaves and seeds of B. napus. • A number of previously unknown potential regulatory relationships involved in glucosinolate synthesis were identified and this study illustrates how genetic variation can affect a biochemical pathway.

Genetic regulation of glucoraphanin accumulation in Beneforté® broccoli

Diets rich in broccoli (Brassica oleracea var italica) have been associated with maintenance of cardiovascular health and reduction in risk of cancer. These health benefits have been attributed to glucoraphanin that specifically accumulates in broccoli. The development of broccoli with enhanced concentrations of glucoraphanin may deliver greater health benefits. Three high-glucoraphanin F1 broccoli hybrids were developed in independent programmes through genome introgression from the wild species Brassica villosa. Glucoraphanin and other metabolites were quantified in experimental field trials. Global SNP analyses quantified the differential extent of B. villosa introgression The high-glucoraphanin broccoli hybrids contained 2.5–3 times the glucoraphanin content of standard hybrids due to enhanced sulphate assimilation and modifications in sulphur partitioning between sulphur-containing metabolites. All of the high-glucoraphanin hybrids possessed an introgressed B. villosa segment which contained a B. villosa Myb28 allele. Myb28 expression was increased in all of the high-glucoraphanin hybrids. Two high-glucoraphanin hybrids have been commercialised as Beneforté® broccoli. The study illustrates the translation of research on glucosinolate genetics from Arabidopsis to broccoli, the use of wild Brassica species to develop cultivars with potential consumer benefits, and the development of cultivars with contrasting concentrations of glucoraphanin for use in blinded human intervention studies.

The CACTA transposon Bot1 played a major role in Brassica genome divergence and gene proliferation

We isolated and characterized a Brassica C genome-specific CACTA element, which was designated Bot1 (Brassica oleracea transposon 1). After analysing phylogenetic relationships, copy numbers and sequence similarity of Bot1 and Bot1 analogues in B. oleracea (C genome) versus Brassica rapa (A genome), we concluded that Bot1 has encountered several rounds of amplification in the oleracea genome only, and has played a major role in the recent rapa and oleracea genome divergence. We performed in silico analyses of the genomic organization and internal structure of Bot1, and established which segment of Bot1 is C-genome specific. Our work reports a fully characterized Brassica repetitive sequence that can distinguish the Brassica A and C chromosomes in the allotetraploid Brassica napus, by fluorescent in situ hybridization. We demonstrated that Bot1 carries a host S locus-associated SLL3 gene copy. We speculate that Bot1 was involved in the proliferation of SLL3 around the Brassica genome. The present study reinforces the assumption that transposons are a major driver of genome and gene evolution in higher plants.

Novel insights into seed fatty acid synthesis and modification pathways from genetic diversity and quantitative trait loci analysis of the Brassica C genome

Natural genetic variation in fatty acid synthesis and modification pathways determine the composition of vegetable oils, which are major components of human diet and renewable products. Based on known pathways we combined diversity and genetic analysis of metabolites to infer the existence of enzymes encoded by distinct loci, and associated these with specific elongation steps or subpathways. A total of 107 lines representing different Brassica genepools revealed considerable variation for 18 seed fatty acid products. The effect of genetic variation within a single biochemical step on subsequent products was demonstrated using a correlation matrix of scatterplots, and by calculating relative step yields. Surprisingly, diploid Brassica oleracea segregating populations had a similar range of variation for individual fatty acids as across the whole genepool. This allowed identification of 22 quantitative trait loci (QTL) associated with activity in the plastid, early stages of synthesis, desaturation, and elongases. Four QTL were assigned to early stages of synthesis, seven to subpathway specific or general elongase activity, one to ketoacyl acyl-carrier protein synthetase, and two each to fatty acid desaturase and either desaturase or fatty acyl-carrier protein thioesterase. An additional 10 QTL had distinct effects but were not assigned specific functions. Where contrasting behavior in more than one subpathway was detected, we inferred QTL specificity for particular combinations of substrate and product. The assignment of enzyme function to QTL was consistent with the known position of some Brassicaeae candidate genes and collinear regions of the Arabidopsis (Arabidopsis thaliana) genome.

The pangenome of an agronomically important crop plant Brassica oleracea

There is an increasing awareness that as a result of structural variation, a reference sequence representing a genome of a single individual is unable to capture all of the gene repertoire found in the species. A large number of genes affected by presence/absence and copy number variation suggest that it may contribute to phenotypic and agronomic trait diversity. Here we show by analysis of the Brassica oleracea pangenome that nearly 20% of genes are affected by presence/absence variation. Several genes displaying presence/absence variation are annotated with functions related to major agronomic traits, including disease resistance, flowering time, glucosinolate metabolism and vitamin biosynthesis.

Turnip mosaic virus (TuMV) Is Able to Use Alleles of Both eIF4E and eIF(iso)4E from Multiple Loci of the Diploid Brassica rapa

Three copies of eIF4E and three copies of eIF(iso)4E have been identified and sequenced from a Turnip mosaic virus (TuMV)-susceptible, inbred, diploid Brassica rapa line, R-o-18. One of the copies of eIF4E lacked exons 2 and 3 and appeared to be a pseudogene. The two other copies of eIF4E and two of the three copies of eIF(iso)4E were isolated from a bacterial artificial chromosome library of R-o-18. Using an Arabidopsis line (Col-0::dSpm) with a transposon knock-out of the eIF(iso)4E gene which resulted in a change from complete susceptibility to complete resistance to TuMV, complementation experiments were carried out with the two versions of eIF4E and the two versions of eIF(iso)4E. When transformed into Col-0::dSpm, all four Brassica transgenes complemented the Arabidopsis eIF(iso)4E knock-out, conferring susceptibility to both mechanical and aphid challenge with TuMV. One of the copies of eIF4E did not appear to support viral replication as successfully as the other copy of eIF4E or the two copies of eIF(iso)4E. The results show that TuMV can use both eIF4E and eIF(iso)4E from B. rapa for replication and, for the first time, that a virus can use eIF4E and eIF(iso)4E from multiple loci of a single host plant.

A rapid method for the determination of artemisinin and its biosynthetic precursors in Artemisia annua L. crude extracts

A rapid high-pressure liquid chromatography (HPLC) tandem mass spectrometry (TQD) method for the determination of artemisinin, 9-epi-artemisinin, artemisitene, dihydroartemisinic acid, artemisinic acid and arteannuin B in Artemisia annua extracts is described. Detection and quantification of 9-epi-artemisinin in crude extracts are reported for the first time. In this method all six metabolites are resolved and eluted within 6 min with minimal sample preparation. A recovery of between 96.25% and 103.59% was obtained for all metabolites analysed and the standard curves were linear (r(2)>0.99) over the concentration range of 0.15-10 μg mL(-1) for artemisinin, 9-epi-artemisinin, artemisitene and arteannuin B, and the range of 3.75-120 μg mL(-1) for dihydroartemisinic acid and artemisinic acid. All validation indices were satisfactory, showing the method to be robust, quick, sensitive and adequate for a range of applications including high throughput (HTP) analysis.