Aaron J. Windsor

The draft genome of the transgenic tropical fruit tree papaya (Carica papaya Linnaeus)

Papaya, a fruit crop cultivated in tropical and subtropical regions, is known for its nutritional benefits and medicinal applications. Here we report a 3× draft genome sequence of ‘SunUp’ papaya, the first commercial virus-resistant transgenic fruit tree to be sequenced. The papaya genome is three times the size of the Arabidopsis genome, but contains fewer genes, including significantly fewer disease-resistance gene analogues. Comparison of the five sequenced genomes suggests a minimal angiosperm gene set of 13,311. A lack of recent genome duplication, atypical of other angiosperm genomes sequenced so far, may account for the smaller papaya gene number in most functional groups. Nonetheless, striking amplifications in gene number within particular functional groups suggest roles in the evolution of tree-like habit, deposition and remobilization of starch reserves, attraction of seed dispersal agents, and adaptation to tropical daylengths. Transgenesis at three locations is closely associated with chloroplast insertions into the nuclear genome, and with topoisomerase I recognition sites. Papaya offers numerous advantages as a system for fruit-tree functional genomics, and this draft genome sequence provides the foundation for revealing the basis of Carica’s distinguishing morpho-physiological, medicinal and nutritional properties.

Genome Wide Analyses Reveal Little Evidence for Adaptive Evolution in Many Plant Species

The relative contribution of advantageous and neutral mutations to the evolutionary process is a central problem in evolutionary biology. Current estimates suggest that whereas Drosophila, mice, and bacteria have undergone extensive adaptive evolution, hominids show little or no evidence of adaptive evolution in protein-coding sequences. This may be a consequence of differences in effective population size. To study the matter further, we have investigated whether plants show evidence of adaptive evolution using an extension of the McDonald-Kreitman test that explicitly models slightly deleterious mutations by estimating the distribution of fitness effects of new mutations. We apply this method to data from nine pairs of species. Altogether more than 2,400 loci with an average length of approximately 280 nucleotides were analyzed. We observe very similar results in all species; we find little evidence of adaptive amino acid substitution in any comparison except sunflowers. This may be because many plant species have modest effective population sizes.

A Gain-of-Function Polymorphism Controlling Complex Traits and Fitness in Nature

Natural Selection at Work Catching the evolution of a novel function and determining its selective parameters in nature remains an extremely difficult task. Prasad et al. (p. 1081) have undertaken this quest documenting the molecular basis of a natural allelic polymorphism and its effects on herbivory and survival in the Arabidopsis relative, Boechera stricta, living in the Rocky Mountains. Positive selection for a mutation that enhances resistance to herbivory in the model plant Boechera is described. Identification of the causal genes that control complex trait variation remains challenging, limiting our appreciation of the evolutionary processes that influence polymorphisms in nature. We cloned a quantitative trait locus that controls plant defensive chemistry, damage by insect herbivores, survival, and reproduction in the natural environments where this polymorphism evolved. These ecological effects are driven by duplications in the BCMA (branched-chain methionine allocation) loci controlling this variation and by two selectively favored amino acid changes in the glucosinolate-biosynthetic cytochrome P450 proteins that they encode. These changes cause a gain of novel enzyme function, modulated by allelic differences in catalytic rate and gene copy number. Ecological interactions in diverse environments likely contribute to the widespread polymorphism of this biochemical function.

Comparative genetic mapping in Boechera stricta, a close relative of Arabidopsis

The angiosperm family Brassicaceae contains both the research model Arabidopsis (Arabidopsis thaliana) and the agricultural genus Brassica. Comparative genomics in the Brassicaceae has largely focused on direct comparisons between Arabidopsis and the species of interest. However, the reduced genome size and chromosome number (n = 5) of Arabidopsis complicates comparisons. Arabidopsis shows extensive genome and chromosome reshuffling compared to its close relatives Arabidopsis lyrata and Capsella rubella, both with n = 8. To facilitate comparative genomics across the Brassicaceae we recently outlined a system of 24 conserved chromosomal blocks based on their positions in an ancestral karyotype of n = 8, rather than by their position in Arabidopsis. In this report we use this system as a tool to understand genome structure and evolution in Boechera stricta (n = 7). B. stricta is a diploid, sexual, and highly self-fertilizing species occurring in mostly montane regions of western North America. We have created an F2 genetic map of B. stricta based on 192 individuals scored at 196 microsatellite and candidate gene loci. Single-nucleotide polymorphism genotyping of 94 of the loci was done simultaneously using an Illumina bead array. The total map length is 725.8 cM, with an average marker spacing of 3.9 cM. There are no gaps greater than 19.3 cM. The chromosomal reduction from n = 8 to n = 7 and other genomic changes in B. stricta likely involved a pericentric inversion, a chromosomal fusion, and two reciprocal translocations that are easily visualized using the genomic blocks. Our genetic map will facilitate the analysis of ecologically relevant quantitative variation in Boechera.

Multilocus Patterns of Nucleotide Diversity, Population Structure and Linkage Disequilibrium in Boechera stricta, a Wild Relative of Arabidopsis

Information about polymorphism, population structure, and linkage disequilibrium (LD) is crucial for association studies of complex trait variation. However, most genomewide studies have focused on model systems, with very few analyses of undisturbed natural populations. Here, we sequenced 86 mapped nuclear loci for a sample of 46 genotypes of Boechera stricta and two individuals of B. holboellii, both wild relatives of Arabidopsis. Isolation by distance was significant across the species range of B. stricta, and three geographic groups were identified by structure analysis, principal coordinates analysis, and distance-based phylogeny analyses. The allele frequency spectrum indicated a genomewide deviation from an equilibrium neutral model, with silent nucleotide diversity averaging 0.004. LD decayed rapidly, declining to background levels in ∼10 kb or less. For tightly linked SNPs separated by <1 kb, LD was dependent on the reference population. LD was lower in the specieswide sample than within populations, suggesting that low levels of LD found in inbreeding species such as B. stricta, Arabidopsis thaliana, and barley may result from broad geographic sampling that spans heterogeneous genetic groups. Finally, analyses also showed that inbreeding B. stricta and A. thaliana have ∼45% higher recombination per kilobase than outcrossing A. lyrata.

Ecological genomics of Boechera stricta: Identification of a QTL controlling the allocation of methionine- vs branched-chain amino acid-derived glucosinolates and levels of insect herbivory

In the Brassicaceae, glucosinolates influence the feeding, reproduction and development of many insect herbivores. Glucosinolate production and effects on herbivore feeding have been extensively studied in the model species, Arabidopsis thaliana and Brassica crops, both of which constitutively produce leaf glucosinolates mostly derived from the amino acid, methionine. Much less is known about the regulation or role in defense of glucosinolates derived from other aliphatic amino acids, such as the branched-chain amino acids (BCAA), valine and isoleucine. We have identified a glucosinolate polymorphism in Boechera stricta controlling the allocation to BCAA- vs methionine-derived glucosinolates in both leaves and seeds. B. stricta is a perennial species that grows in mostly undisturbed habitats of western North America. We have measured glucosinolate profiles and concentrations in 192 F2 lines that have earlier been used for genetic map construction. We also performed herbivory assays on six F3 replicates per F2 line using the generalist lepidopteran, Trichoplusia ni. Quantitative trait locus (QTL) analysis identified a single locus controlling both glucosinolate profile and levels of herbivory, the branched chain-methionine allocation or BCMA QTL. We have delimited this QTL to a small genomic region with a 1.0 LOD confidence interval just 1.9 cm wide, which, in A. thaliana, contains ∼100 genes. We also found that methionine-derived glucosinolates provided significantly greater defense than the BCAA-derived glucosinolates against feeding by this generalist insect herbivore. The future positional cloning of this locus will allow for testing various adaptive explanations.

Partial Shotgun Sequencing of the Boechera stricta Genome Reveals Extensive Microsynteny and Promoter Conservation with Arabidopsis

Comparative genomics provides insight into the evolutionary dynamics that shape discrete sequences as well as whole genomes. To advance comparative genomics within the Brassicaceae, we have end sequenced 23,136 medium-sized insert clones from Boechera stricta, a wild relative of Arabidopsis (Arabidopsis thaliana). A significant proportion of these sequences, 18,797, are nonredundant and display highly significant similarity (BLASTn e-value ≤ 10−30) to low copy number Arabidopsis genomic regions, including more than 9,000 annotated coding sequences. We have used this dataset to identify orthologous gene pairs in the two species and to perform a global comparison of DNA regions 5′ to annotated coding regions. On average, the 500 nucleotides upstream to coding sequences display 71.4% identity between the two species. In a similar analysis, 61.4% identity was observed between 5′ noncoding sequences of Brassica oleracea and Arabidopsis, indicating that regulatory regions are not as diverged among these lineages as previously anticipated. By mapping the B. stricta end sequences onto the Arabidopsis genome, we have identified nearly 2,000 conserved blocks of microsynteny (bracketing 26% of the Arabidopsis genome). A comparison of fully sequenced B. stricta inserts to their homologous Arabidopsis genomic regions indicates that indel polymorphisms >5 kb contribute substantially to the genome size difference observed between the two species. Further, we demonstrate that microsynteny inferred from end-sequence data can be applied to the rapid identification and cloning of genomic regions of interest from nonmodel species. These results suggest that among diploid relatives of Arabidopsis, small- to medium-scale shotgun sequencing approaches can provide rapid and cost-effective benefits to evolutionary and/or functional comparative genomic frameworks.

The shrunken genome of Arabidopsis thaliana

This paper examines macro and micro-level patterns of genome size evolution in the Brassicaceae. A phylogeny of 25 relatives of Arabidopsis thaliana was reconstructed using four molecular markers under both parsimony and Bayesian methods. Reconstruction of genome size (C value) evolution as a discrete character and as a continuous character was also performed. In addition, size dynamics in small chromosomal regions were assessed by comparing genomic clones generated for Arabidopsis lyrata and for Boechera stricta to the fully sequenced genome of A. thaliana. The results reveal a sevenfold variation in genome size among the taxa investigated and that the small genome size of A. thaliana is derived. Our results also indicate that the genome is free to increase or decrease in size across these evolutionary lineages without a directional bias. These changes are accomplished by insertions and deletions at both large and small-scales occurring mostly in intergenic regions, with repetitive sequences and transposable elements implicated in genome size increases. The focus upon taxa relatively closely related to the model organism A. thaliana, and the combination of complementary approaches, allows for unique insights into the processes driving genome size changes.