Andrew G. Sharpe

Early allopolyploid evolution in the post-neolithic Brassica napus oilseed genome

The genomic origins of rape oilseed Many domesticated plants arose through the meeting of multiple genomes through hybridization and genome doubling, known as polyploidy. Chalhoub et al. sequenced the polyploid genome of Brassica napus, which originated from a recent combination of two distinct genomes approximately 7500 years ago and gave rise to the crops of rape oilseed (canola), kale, and rutabaga. B. napus has undergone multiple events affecting differently sized genetic regions where a gene from one progenitor species has been converted to the copy from a second progenitor species. Some of these gene conversion events appear to have been selected by humans as part of the process of domestication and crop improvement. Science, this issue p. 950 The polyploid genome of oilseed rape exhibits evolution through homologous gene conversion. Oilseed rape (Brassica napus L.) was formed ~7500 years ago by hybridization between B. rapa and B. oleracea, followed by chromosome doubling, a process known as allopolyploidy. Together with more ancient polyploidizations, this conferred an aggregate 72× genome multiplication since the origin of angiosperms and high gene content. We examined the B. napus genome and the consequences of its recent duplication. The constituent An and Cn subgenomes are engaged in subtle structural, functional, and epigenetic cross-talk, with abundant homeologous exchanges. Incipient gene loss and expression divergence have begun. Selection in B. napus oilseed types has accelerated the loss of glucosinolate genes, while preserving expansion of oil biosynthesis genes. These processes provide insights into allopolyploid evolution and its relationship with crop domestication and improvement.

Draft genome sequence of chickpea (Cicer arietinum) provides a resource for trait improvement

Chickpea (Cicer arietinum) is the second most widely grown legume crop after soybean, accounting for a substantial proportion of human dietary nitrogen intake and playing a crucial role in food security in developing countries. We report the ∼738-Mb draft whole genome shotgun sequence of CDC Frontier, a kabuli chickpea variety, which contains an estimated 28,269 genes. Resequencing and analysis of 90 cultivated and wild genotypes from ten countries identifies targets of both breeding-associated genetic sweeps and breeding-associated balancing selection. Candidate genes for disease resistance and agronomic traits are highlighted, including traits that distinguish the two main market classes of cultivated chickpea—desi and kabuli. These data comprise a resource for chickpea improvement through molecular breeding and provide insights into both genome diversity and domestication.

The Brassica oleracea genome reveals the asymmetrical evolution of polyploid genomes

Polyploidization has provided much genetic variation for plant adaptive evolution, but the mechanisms by which the molecular evolution of polyploid genomes establishes genetic architecture underlying species differentiation are unclear. Brassica is an ideal model to increase knowledge of polyploid evolution. Here we describe a draft genome sequence of Brassica oleracea, comparing it with that of its sister species B. rapa to reveal numerous chromosome rearrangements and asymmetrical gene loss in duplicated genomic blocks, asymmetrical amplification of transposable elements, differential gene co-retention for specific pathways and variation in gene expression, including alternative splicing, among a large number of paralogous and orthologous genes. Genes related to the production of anticancer phytochemicals and morphological variations illustrate consequences of genome duplication and gene divergence, imparting biochemical and morphological variation to B. oleracea. This study provides insights into Brassica genome evolution and will underpin research into the many important crops in this genus.

The genome of flax (Linum usitatissimum) assembled de novo from short shotgun sequence reads

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(50) =694 kb, including contigs with N(50)=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 43384 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.

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.

Desaturase multigene families of Brassica napus arose through genome duplication

Abstract This paper reports the estimated gene copy number and restriction fragment length polymorphism (RFLP) map locations of five different desaturase cDNA clones from Brassica napus (oilseed rape). The desaturase enzymes encoded by four of these genes catalyze successive reactions that insert double bonds into lipid-linked fatty acid residues. Delta-12 (e2) and delta-15 (e3) desaturases are active in the endoplasmic reticulum, while omega-6 (p2) and omega-3 (p3) desaturases catalyze analogous desaturation reactions via a parallel pathway located in plastids. The fifth cDNA clone (b5) contains a desaturase-like domain bound to a cytochrome b5 segment. Estimates of gene copy number based on Southern blot analysis of 16 oilseed rape varieties and three different resynthesized Brassica napus lines indicated that e2 had 4–6 gene copies and e3, p2, p3 and b5 each had 6–8 gene copies per haploid genome. Estimates of the gene copy number for the two progenitor species, Brassica oleracea and Brassica rapa, supported the premise that all these genes were at least duplicated or triplicated in the two progenitor species before they combined to form B. napus. RFLP mapping results showed that the e2 probe detected 4 distinct loci, the e3 probe 6 loci and p2, p3 and b5 each detected 8 loci, with pairs of loci often mapping to homoeologous regions on 2 different linkage groups. The 28 mappable loci were distributed across 12 linkage groups of the B. napus map (Parkin et al. 1995) and were usually represented by single RFLP fragments. A collinear segment containing the e2 and p3 loci was positioned on B. napus linkage groups N1, N11, N3, N13, N5 and N15. This segment was collinear with a 30-cM region of Arabidopsis thaliana chromosome 3 that contains the homologous fad2 (e2) and fad7(p3) genes. This suggests that the desaturase multigene families arose as the result of duplication of large chromosome segments rather than duplication of individual genes.

The emerging biofuel crop Camelina sativa retains a highly undifferentiated hexaploid genome structure

Camelina sativa is an oilseed with desirable agronomic and oil-quality attributes for a viable industrial oil platform crop. Here we generate the first chromosome-scale high-quality reference genome sequence for C. sativa and annotated 89,418 protein-coding genes, representing a whole-genome triplication event relative to the crucifer model Arabidopsis thaliana. C. sativa represents the first crop species to be sequenced from lineage I of the Brassicaceae. The well-preserved hexaploid genome structure of C. sativa surprisingly mirrors those of economically important amphidiploid Brassica crop species from lineage II as well as wheat and cotton. The three genomes of C. sativa show no evidence of fractionation bias and limited expression-level bias, both characteristics commonly associated with polyploid evolution. The highly undifferentiated polyploid genome of C. sativa presents significant consequences for breeding and genetic manipulation of this industrial oil crop.

Alignment of the conserved C genomes of Brassica oleracea and Brassica napus

A population of 169 microspore-derived doubled-haploid lines was produced from a highly polymorphic Brassica oleracea cross. A dense genetic linkage map of B. oleracea was then developed based on the segregation of 303 RFLP-defined loci. It is hoped that these lines will be used by other geneticists to facilitate the construction of a unified genetic map of B. oleracea. When the B. oleracea map was compared to one ofB. napus (Parkin et al. 1995), based on the same RFLP probes (Sharpe et al. 1995), good collinearity between the C-genome linkage groups of the two species was observed.

Conserved structure and function of the Arabidopsis flowering time gene CONSTANS in Brassica napus.

The Arabidopsis thaliana CONSTANS (CO) gene which promotes flowering in long days was recently isolated by chromosome walking. The mapping of QTLs controlling flowering time in Brassica species has identified genomic regions that contain homologues of the CO gene. Four genes homologous to the Arabidopsis CO gene were isolated from a pair of homoeologous loci in each of two doubled-haploid Brassica napus lines displaying different flowering times, N-o-1 and N-o-9. The four genes, BnCOa1, BnCOa9, BnCOb1 and BnCOb9, are located on linkage groups N10 and N19, and are highly similar to each other and to the Arabidopsis CO gene. Two regions of the proteins are particularly well conserved, a N-terminal region with two putative zinc fingers and a C-terminal region which may contain a nuclear localization signal. All four genes appear to be expressed in B. napus. The BnCOa1 allele was shown to complement the co-2 mutation in Arabidopsis in a dosage-dependent manner causing earlier flowering than in wild type under both long- and short-day conditions.

A haplotype map of allohexaploid wheat reveals distinct patterns of selection on homoeologous genomes.

BackgroundBread wheat is an allopolyploid species with a large, highly repetitive genome. To investigate the impact of selection on variants distributed among homoeologous wheat genomes and to build a foundation for understanding genotype-phenotype relationships, we performed population-scale re-sequencing of a diverse panel of wheat lines.ResultsA sample of 62 diverse lines was re-sequenced using the whole exome capture and genotyping-by-sequencing approaches. We describe the allele frequency, functional significance, and chromosomal distribution of 1.57 million single nucleotide polymorphisms and 161,719 small indels. Our results suggest that duplicated homoeologous genes are under purifying selection. We find contrasting patterns of variation and inter-variant associations among wheat genomes; this, in addition to demographic factors, could be explained by differences in the effect of directional selection on duplicated homoeologs. Only a small fraction of the homoeologous regions harboring selected variants overlapped among the wheat genomes in any given wheat line. These selected regions are enriched for loci associated with agronomic traits detected in genome-wide association studies.ConclusionsEvidence suggests that directional selection in allopolyploids rarely acted on multiple parallel advantageous mutations across homoeologous regions, likely indicating that a fitness benefit could be obtained by a mutation at any one of the homoeologs. Additional advantageous variants in other homoelogs probably either contributed little benefit, or were unavailable in populations subjected to directional selection. We hypothesize that allopolyploidy may have increased the likelihood of beneficial allele recovery by broadening the set of possible selection targets.