Brassicaceae phylogeny and trichome evolution

Identification of shared single copy nuclear genes in Arabidopsis, Populus, Vitis and Oryza and their phylogenetic utility across various taxonomic levels

Background

Although the overwhelming majority of genes found in angiosperms are members of gene families, and both gene- and genome-duplication are pervasive forces in plant genomes, some genes are sufficiently distinct from all other genes in a genome that they can be operationally defined as 'single copy'. Using the gene clustering algorithm MCL-tribe, we have identified a set of 959 single copy genes that are shared single copy genes in the genomes of Arabidopsis thaliana, Populus trichocarpa, Vitis vinifera and Oryza sativa. To characterize these genes, we have performed a number of analyses examining GO annotations, coding sequence length, number of exons, number of domains, presence in distant lineages, such as Selaginella and Physcomitrella, and phylogenetic analysis to estimate copy number in other seed plants and to demonstrate their phylogenetic utility. We then provide examples of how these genes may be used in phylogenetic analyses to reconstruct organismal history, both by using extant coverage in EST databases for seed plants and de novo amplification via RT-PCR in the family Brassicaceae.

Results

There are 959 single copy nuclear genes shared in Arabidopsis, Populus, Vitis and Oryza ["APVO SSC genes"]. The majority of these genes are also present in the Selaginella and Physcomitrella genomes. Public EST sets for 197 species suggest that most of these genes are present across a diverse collection of seed plants, and appear to exist as single or very low copy genes, though exceptions are seen in recently polyploid taxa and in lineages where there is significant evidence for a shared large-scale duplication event. Genes encoding proteins localized in organelles are more commonly single copy than expected by chance, but the evolutionary forces responsible for this bias are unknown.

Regardless of the evolutionary mechanisms responsible for the large number of shared single copy genes in diverse flowering plant lineages, these genes are valuable for phylogenetic and comparative analyses. Eighteen of the APVO SSC single copy genes were amplified in the Brassicaceae using RT-PCR and directly sequenced. Alignments of these sequences provide improved resolution of Brassicaceae phylogeny compared to recent studies using plastid and ITS sequences. An analysis of sequences from 13 APVO SSC genes from 69 species of seed plants, derived mainly from public EST databases, yielded a phylogeny that was largely congruent with prior hypotheses based on multiple plastid sequences. Whereas single gene phylogenies that rely on EST sequences have limited bootstrap support as the result of limited sequence information, concatenated alignments result in phylogenetic trees with strong bootstrap support for already established relationships. Overall, these single copy nuclear genes are promising markers for phylogenetics, and contain a greater proportion of phylogenetically-informative sites than commonly used protein-coding sequences from the plastid or mitochondrial genomes.

Conclusions

Putatively orthologous, shared single copy nuclear genes provide a vast source of new evidence for plant phylogenetics, genome mapping, and other applications, as well as a substantial class of genes for which functional characterization is needed. Preliminary evidence indicates that many of the shared single copy nuclear genes identified in this study may be well suited as markers for addressing phylogenetic hypotheses at a variety of taxonomic levels.

Molecular phylogeny and systematics of the genus Draba (Brassicaceae) and identification of its most closely related genera

Phylogenetic analyses of sequence data from the internal transcribed spacer regions of the nuclear ribosomal DNA (ITS) and the plastid trnL-F region were conducted to reconstruct the phylogenetic history of Draba and related genera. Out of the approximately 370 Draba species, 169 geographically and morphologically representative species are sampled here, including such "controversial" segregates as Abdra, Arabis, Athysanus, Drabopsis, Erophila, Graellsia, Heterodraba, Schivereckia and Tomostima. Phylogenetic analyses of the combined ITS and trnL-F markers indicate that Draba represents a monophyletic evolutionary lineage within the tribe Arabideae, but previously proposed infrageneric classification systems are mostly artificial (i.e. do not reflect true evolutionary history). Draba muralis formed an independent genus, possibly between Draba and Arabis, whereas D. hystrix fell outside Draba and was closely associated with Arabis. The New World annual, lowland Draba, D. platycarpa, D. reptans, D. cuneifolia, D. australis, D. arabidoides, (section Abdra) and D. brachycarpa and D. aspera (section Tomostima), appear to be independent genera that fall outside Draba and are monophyletic with the endemic North American Heterodraba unilateralis and Athysanus pusillus. Graellsia hederifolia and Erophila verna appear to be earlier diverging Draba species, with weak evidence of ancient hybridization in G. hederifolia. Core Draba species were organized into three major groups that encompass the segregate genera Drabopsis and Schivereckia. The three groups have geographic significance: Group I-Europe to Iran; Group II-North and South American Cordillera; Group III-Asian, Arctic, and Beringian. These three groups also have significant petal color and base chromosome level similarities. We also found that Arabis rimarum is a synonym for Draba aucheri, and is well within Group I of Core Draba. Therefore, we propose that many of the lowland and coastal temperate species often associated with Draba should be excluded from that genus. The compilation of these data indicate that true Draba are highly migrating arctic and higher alpine species that are most often perennial with white petals and n=8 as plesiomorphic characters.

Genome-wide discovery of DNA polymorphism in Brassica rapa

Single nucleotide polymorphisms (SNPs) and/or insertion/deletions (InDels) are frequent sequence variations in the plant genome, which can be developed as molecular markers for genetic studies on crop improvement. The ongoing Brassica rapa genome sequencing project has generated vast amounts of sequence data useful in genetic research. Here, we report a genome-wide survey of DNA polymorphisms in the B. rapa genome based on the 557 bacterial artificial clone sequences of B. rapa ssp. pekinensis cv. Chiifu. We identified and characterized 21,311 SNPs and 6,753 InDels in the gene space of the B. rapa genome by re-sequencing 1,398 sequence-tagged sites (STSs) in eight genotypes. Comparison of our findings with a B. rapa genetic linkage map confirmed that STS loci were distributed randomly over the B. rapa whole genome. In the 1.4 Mb of aligned sequences, mean nucleotide polymorphism and diversity were theta = 0.00890 and pi = 0.00917, respectively. Additionally, the nucleotide diversity in introns was almost three times greater than that in exons, and the frequency of observed InDel was almost 17 times higher in introns than in exons. Information regarding SNPs/InDels obtained here will provide an important resource for genetic studies and breeding programs of B. rapa.

Genome-wide comparative analysis of the Brassica rapa gene space reveals genome shrinkage and differential loss of duplicated genes after whole genome triplication

Euchromatic regions of the Brassica rapa genome were sequenced and mapped onto the corresponding regions in the Arabidopsis thaliana genome.

Background

Brassica rapa is one of the most economically important vegetable crops worldwide. Owing to its agronomic importance and phylogenetic position, B. rapa provides a crucial reference to understand polyploidy-related crop genome evolution. The high degree of sequence identity and remarkably conserved genome structure between Arabidopsis and Brassica genomes enables comparative tiling sequencing using Arabidopsis sequences as references to select the counterpart regions in B. rapa, which is a strong challenge of structural and comparative crop genomics.

Results

We assembled 65.8 megabase-pairs of non-redundant euchromatic sequence of B. rapa and compared this sequence to the Arabidopsis genome to investigate chromosomal relationships, macrosynteny blocks, and microsynteny within blocks. The triplicated B. rapa genome contains only approximately twice the number of genes as in Arabidopsis because of genome shrinkage. Genome comparisons suggest that B. rapa has a distinct organization of ancestral genome blocks as a result of recent whole genome triplication followed by a unique diploidization process. A lack of the most recent whole genome duplication (3R) event in the B. rapa genome, atypical of other Brassica genomes, may account for the emergence of B. rapa from the Brassica progenitor around 8 million years ago.

Conclusions

This work demonstrates the potential of using comparative tiling sequencing for genome analysis of crop species. Based on a comparative analysis of the B. rapa sequences and the Arabidopsis genome, it appears that polyploidy and chromosomal diploidization are ongoing processes that collectively stabilize the B. rapa genome and facilitate its evolution.

Plant traits correlated with generation time directly affect inbreeding depression and mating system and indirectly genetic structure

Background

Understanding the mechanisms that control species genetic structure has always been a major objective in evolutionary studies. The association between genetic structure and species attributes has received special attention. As species attributes are highly taxonomically constrained, phylogenetically controlled methods are necessary to infer causal relationships. In plants, a previous study controlling for phylogenetic signal has demonstrated that Wright's F_ST, a measure of genetic differentiation among populations, is best predicted by the mating system (outcrossing, mixed-mating or selfing) and that plant traits such as perenniality and growth form have only an indirect influence on F_ST via their association with the mating system. The objective of this study is to further outline the determinants of plant genetic structure by distinguishing the effects of mating system on gene flow and on genetic drift. The association of biparental inbreeding and inbreeding depression with population genetic structure, mating system and plant traits are also investigated.

Results

Based on data from 263 plant species for which estimates of F_ST, inbreeding (F_IS) and outcrossing rate (t_m) are available, we confirm that mating system is the main influencing factor of F_ST. Moreover, using an alternative measure of F_ST unaffected by the impact of inbreeding on effective population size, we show that the influence of t_m on F_ST is due to its impact on gene flow (reduced pollen flow under selfing) and on genetic drift (higher drift under selfing due to inbreeding). Plant traits, in particular perenniality, influence F_ST mostly via their effect on the mating system but also via their association with the magnitude of selection against inbred individuals: the mean inbreeding depression increases from short-lived herbaceous to long-lived herbaceous and then to woody species. The influence of perenniality on mating system does not seem to be related to differences in stature, as proposed earlier, but rather to differences in generation time.

Conclusion

Plant traits correlated with generation time affect both inbreeding depression and mating system. These in turn modify genetic drift and gene flow and ultimately genetic structure.

Tests of adaptation: functional studies of pollen removal and estimates of natural selection on anther position in wild radish

BACKGROUND

There are a number of difficulties associated with the study of adaptation. One is a lack of variation in the trait, which is common in adaptations because past selection has removed unfit variants. This lack of variation makes it difficult to determine the relationship between trait variation and fitness. Another difficulty is proving causation in this trait-fitness relationship, because a correlated trait might be the actual adaptation. These difficulties can be ameliorated at least partially by combining studies of natural variation with studies of experimentally manipulated traits and traits whose variance has been augmented by artificial selection.

SCOPE

We review here a number of our studies on the adaptive value of two aspects of anther position in wild radish (Raphanus raphanistrum, Brassicaceae): anther exsertion, i.e. the degree to which anthers protrude from the mouth of the corolla tube, and anther height dimorphism, i.e. the difference in lengths of the filaments between the two short and four long stamens. We have used both functional analyses, in which the response variable is pollen removal, and measurements of selection, in which the response variable is lifetime male fitness estimated by molecular genetic paternity analyses. In these studies we use both the natural variation in populations as well as manipulated variation, the latter through both stamen removal and artificial selection, to re-create the ancestral trait conditions.

CONCLUSIONS

Our work provides convincing evidence that intermediate anther exsertion values are adaptive, and that this is probably an adaptation to a subset of the pollinator fauna, small bees. The picture for anther height dimorphism is much less clear, as the weight of current evidence suggests that current values of this trait might actually be maladaptive; however, if this is true it is difficult to understand how the dimorphism is maintained across the family Brassicaceae.

Identification of highly variable chloroplast sequences and development of cpDNA-based molecular markers that distinguish four cytoplasm types in radish (Raphanus sativus L.)

Four types of cytoplasms (Ogura, DCGMS, DBRMF1, and DBRMF2) were identified in the previous studies using molecular markers based on mitochondrial genome variations in radish (Raphanus sativus L.). However, mtDNA markers have limitations in obtaining clear results due to complexity of radish mitochondrial genomes. To improve fidelity, molecular markers based on variation of chloroplast genome sequences were developed in this study. We searched for the sequence variations of chloroplast genome among the four cytoplasm types in 11 noncoding intergenic regions of ~8.7 kb. Highly variable intergenic regions between trnK and rps16 were identified, and a couple of 4-34 bp indels were used to develop a simple PCR-based marker that distinguished the four cytoplasm types based on the PCR product length polymorphism. Two additional cpDNA markers were developed by using a single nucleotide polymorphism and 17-bp insertion. Analysis of 90 accessions using both mtDNA and cpDNA markers showed the perfect match of results of both the markers, suggesting strict co-transmission of mitochondria and chloroplast in radish. Phylogenetic trees showed that two male-sterility inducing cytoplasms, Ogura and DCGMS, were closely related to DBRMF1 and DBRMF2, respectively. Analysis of 120 radish germplasms introduced from diverse countries showed that the frequency of male-sterility inducing mitotypes of Ogura and DCGMS was very low, and DCGMS was predominately detected in eastern European countries. Majority of accessions from Europe and Asia were shown to contain DBRMF2 and DBRMF1 mitotypes, respectively.

Phylogenetic relationships among diploid species of Symphyotrichum (Asteraceae: Astereae) based on two nuclear markers, ITS and GAPDH

The mostly North American subtribe Symphyotrichinae (Asteraceae: Astereae) comprises Canadanthus, Ampelaster, Psilactis, Almutaster, and Symphyotrichum. Intergeneric and interspecific relationships within the subtribe have been investigated in the past, particularly by Nesom [Nesom, G.L., 1994. Review of the taxonomy of Aster sensu lato (Asteraceae: Astereae), emphasizing the new world species, Phytologia 77, 141-297] and Semple [Semple, J.C., 2005. Classification of Symphyotrichum. Available from: <http://www.jcsemple.uwaterloo.ca/Symphyotrichumclassification.htm/>], using morphological and cytological approaches. Symphyotrichum is the largest and most complex genus within the subtribe and includes four subgenera: Symphyotrichum (x=7, 8), Virgulus (x=4, 5), Astropolium (x=5), and Chapmaniani (x=7). In this study we used two nuclear markers, the nrDNA internal transcribed spacer (ITS) and the low-copy nuclear gene glyceraldehyde 3-phosphate dehydrogenase (GAPDH), to resolve intergeneric and interspecific relationships within the subtribe at the diploid level, and to determine whether our phylogenies validate the classifications of Nesom or Semple. Our results confirm the distinct generic status of Canadanthus and Ampelaster, whereas Psilactis and Almutaster form a polytomy with Symphyotrichum. Within Symphyotrichum, subg. Virgulus is monophyletic based on ITS and appears polyphyletic based on GAPDH. Neither the ITS nor the GAPDH analyses support a distinct status for subg. Astropolium, which groups within subg. Symphyotrichum. In general, interspecific relationships within Symphyotrichum are unresolved. Lack of resolution may be interpreted as a case of recent and rapid evolutionary radiation.

Chapter 4. New model systems for the study of developmental evolution in plants

The number of genetically tractable plant model systems is rapidly increasing, thanks to the decreasing cost of sequencing and the wide amenability of plants to stable transformation and other functional approaches. In this chapter, I discuss emerging model systems from throughout the land plant phylogeny and consider how their unique attributes are contributing to our understanding of development, evolution, and ecology. These new models are being developed using two distinct strategies: in some cases, they are selected because of their close relationship to the established models, while in others, they are chosen with the explicit intention of exploring distantly related plant lineages. Such complementary approaches are yielding exciting new results that shed light on both micro- and macroevolutionary processes in the context of developmental evolution.

Lepidium as a model system for studying the evolution of fruit development in Brassicaceae

Fruits represent a key innovation of the flowering plants that facilitates seed dispersal. In many species of the plant family Brassicaceae dehiscent fruits develop in which seed dispersal occurs through a process termed 'pod-shatter'. In the case of dehiscence, the fruit opens during fruit maturation. Phylogeny reconstructions using molecular markers indicate that the development of dehiscent fruits is the ancestral condition within the genus Lepidium s.l., but that indehiscent fruits evolved independently several times from dehiscent fruits. With Lepidium campestre and Cardaria pubescens (also known as Lepidium appelianum), very closely related taxa with dehiscent and indehiscent fruits, respectively, were identified which constitute a well-suited model system to determine the molecular genetic basis of evolutionary changes in fruit dehiscence. Following the rationale of evolutionary developmental biology ('evo-devo') phylomimicking mutants with indehiscent fruits of the close relative Arabidopsis have been used to define the candidate genes ALC, FUL, IND, RPL, and SHP1/2 which might be involved in the origin of indehiscent fruits in Cardaria. Comparative expression studies in L. campestre and C. pubescens are used to identify differentially expressed genes and thus to narrow down the number of candidate genes. Reciprocal heterologous transformation experiments may help us to distinguish direct from indirect developmental genetic causes of fruit indehiscence, and to assess the contribution of cis- and trans-regulatory changes.

The dynamic ups and downs of genome size evolution in Brassicaceae

Crucifers (Brassicaceae, Cruciferae) are a large family comprising some 338 genera and c. 3,700 species. The family includes important crops as well as several model species in various fields of plant research. This paper reports new genome size (GS) data for more than 100 cruciferous species in addition to previously published C-values (the DNA amount in the unreplicated gametic nuclei) to give a data set comprising 185 Brassicaceae taxa, including all but 1 of the 25 tribes currently recognized. Evolution of GS was analyzed within a phylogenetic framework based on gene trees built from five data sets (matK, chs, adh, trnLF, and ITS). Despite the 16.2-fold variation across the family, most Brassicaceae species are characterized by very small genomes with a mean 1C-value of 0.63 pg. The ancestral genome size (ancGS) for Brassicaceae was reconstructed as (anc)1C=0.50 pg. Approximately 50% of crucifer taxa analyzed showed a decrease in GS compared with the ancGS. The remaining species showed an increase in GS although this was generally moderate, with significant increases in C-value found only in the tribes Anchonieae and Physarieae. Using statistical approaches to analyze GS, evolutionary gains or losses in GS were seen to have accumulated disproportionately faster within longer branches. However, we also found that GS has not changed substantially through time and most likely evolves passively (i.e., a tempo that cannot be distinguished between neutral evolution and weak forms of selection). The data reveal an apparent paradox between the narrow range of small GSs over long evolutionary time periods despite evidence of dynamic genomic processes that have the potential to lead to genome obesity (e.g., transposable element amplification and polyploidy). To resolve this, it is suggested that mechanisms to suppress amplification and to eliminate amplified DNA must be active in Brassicaceae although their control and mode of operation are still poorly understood.

Brassicaceae phylogeny inferred from phytochrome A and ndhF sequence data: tribes and trichomes revisited

The family Brassicaceae comprises 3710 species in 338 genera, 25 recently delimited tribes, and three major lineages based on phylogenetic results from the chloroplast gene ndhF. To assess the credibility of the lineages and newly delimited tribes, we sequenced an approximately 1.8-kb region of the nuclear phytochrome A (PHYA) gene for taxa previously sampled for the chloroplast gene ndhF. Using parsimony, likelihood, and Bayesian methods, we reconstructed the phylogeny of the gene and used the approximately unbiased (AU) test to compare phylogenetic results from PHYA with findings from ndhF. We also combined ndhF and PHYA data and used a Bayesian mixed model approach to infer phylogeny. PHYA and combined analyses recovered the same three large lineages as those recovered in ndhF trees, increasing confidence in these lineages. The combined tree confirms the monophyly of most of the recently delimited tribes (only Alysseae, Anchonieae, and Descurainieae are not monophyletic), while 13 of the 23 sampled tribes are monophyletic in PHYA trees. In addition to phylogenetic results, we documented the trichome branching morphology of species across the phylogeny and explored the evolution of different trichome morphologies using the AU test. Our results indicate that dendritic, medifixed, and stellate trichomes likely evolved independently several times in the Brassicaceae.

Chromosomal phylogeny and karyotype evolution in x=7 crucifer species (Brassicaceae)

Karyotype evolution in species with identical chromosome number but belonging to distinct phylogenetic clades is a long-standing question of plant biology, intractable by conventional cytogenetic techniques. Here, we apply comparative chromosome painting (CCP) to reconstruct karyotype evolution in eight species with x=7 (2n=14, 28) chromosomes from six Brassicaceae tribes. CCP data allowed us to reconstruct an ancestral Proto-Calepineae Karyotype (PCK; n=7) shared by all x=7 species analyzed. The PCK has been preserved in the tribes Calepineae, Conringieae, and Noccaeeae, whereas karyotypes of Eutremeae, Isatideae, and Sisymbrieae are characterized by an additional translocation. The inferred chromosomal phylogeny provided compelling evidence for a monophyletic origin of the x=7 tribes. Moreover, chromosomal data along with previously published gene phylogenies strongly suggest the PCK to represent an ancestral karyotype of the tribe Brassiceae prior to its tribe-specific whole-genome triplication. As the PCK shares five chromosomes and conserved associations of genomic blocks with the putative Ancestral Crucifer Karyotype (n=8) of crucifer Lineage I, we propose that both karyotypes descended from a common ancestor. A tentative origin of the PCK via chromosome number reduction from n=8 to n=7 is outlined. Comparative chromosome maps of two important model species, Noccaea caerulescens and Thellungiella halophila, and complete karyotypes of two purported autotetraploid Calepineae species (2n=4x=28) were reconstructed by CCP.

Evidence of neutral transcriptome evolution in plants

* The transcriptome of an organism is its set of gene transcripts (mRNAs) at a defined spatial and temporal locus. Because gene expression is affected markedly by environmental and developmental perturbations, it is widely assumed that transcriptome divergence among taxa represents adaptive phenotypic selection. This assumption has been challenged by neutral theories which propose that stochastic processes drive transcriptome evolution. * To test for evidence of neutral transcriptome evolution in plants, we quantified 18 494 gene transcripts in nonsenescent leaves of 14 taxa of Brassicaceae using robust cross-species transcriptomics which includes a two-step physical and in silico-based normalization procedure based on DNA similarity among taxa. * Transcriptome divergence correlates positively with evolutionary distance between taxa and with variation in gene expression among samples. Results are similar for pseudogenes and chloroplast genes evolving at different rates. Remarkably, variation in transcript abundance among root-cell samples correlates positively with transcriptome divergence among root tissues and among taxa. * Because neutral processes affect transcriptome evolution in plants, many differences in gene expression among or within taxa may be nonfunctional, reflecting ancestral plasticity and founder effects. Appropriate null models are required when comparing transcriptomes in space and time.

Molecular characterization of Lal2, an SRK-like gene linked to the S-locus in the wild mustard Leavenworthia alabamica

Single-locus sporophytic self-incompatibility inhibits inbreeding in many members of the mustard family (Brassicaceae). To investigate the genetics of self-incompatibility in the wild mustard Leavenworthia alabamica, diallel crosses were conducted between full siblings. Patterns of incompatibility were consistent with the action of single-locus sporophytic self-incompatibility. DNA sequences related to S-locus receptor kinase (SRK), the gene involved in self-pollen recognition in mustards, were cloned and sequenced. A single sequence with high identity to SRK and several other groups of sequences (Lal1, Lal2, Lal3, Lal8, and Lal14) were isolated from L. alabamica. We propose that either Lal2 sequences are divergent alleles of SRK or Lal2 is in tight linkage with SRK because (1) Lal2 alleles cosegregate with S-alleles inferred from dialleles in all 97 cases tested in five families; (2) Lal2 sequences are highly diverse at both synonymous and nonsynonymous sites and exhibit patterns of selective constraint similar to those observed at SRK in Brassica and Arabidopsis; and (3) transcripts of one Lal2 allele were detected in leaves and the styles of open flowers, but were most abundant in the stigmas of maturing buds. We discuss the utility of the S-linked polymorphism at Lal2 for studying the evolutionary forces acting on self-incompatibility in Leavenworthia.

Duplication and adaptive evolution of the COR15 genes within the highly cold-tolerant Draba lineage (Brassicaceae)

Plants have evolved diverse adaptive mechanisms that enable them to tolerate abiotic stresses, to varying degrees, and such stresses may have strongly influenced evolutionary changes at levels ranging from molecular to morphological. Previous studies on these phenomena have focused on the adaptive evolution of stress-related orthologous genes in specific lineages. However, heterogenetic evolution of the paralogous genes following duplication has only been examined in a very limited number of stress-response gene families. The COR15 gene encodes a low molecular weight protein that plays an important role in protecting plants from cold stresses. Although two different copies of this gene have been found in the model species, Arabidopsis thaliana, evolutionary patterns of this small gene family in plants have not been previously explored. In this study, we cloned COR15-like sequences and performed evolutionary analyses of these sequences (including those previously reported) in the highly cold-tolerant Draba lineage and related lineages of Brassicaceae. Our phylogenetic analyses indicate that all COR15-like sequences clustered into four clades that corresponded well to the morphological lineages. Gene conversions were found to have probably occurred before/during the divergence of Brassica and Draba lineage. However, repeated, independent duplications of this gene have occurred in different lineages of Brassicaceae. Further comparisons of all sequences suggest that there have been significant inter-lineage differences in evolutionary rates between the duplicated and original genes. We assessed the likelihood that the differences between two well-supported gene subfamilies that appear to have originated from a single duplication, COR15a and COR15b, within the Draba lineage have been driven by adaptive evolution. Comparisons of their non-synonymous/synonymous substitution ratios and rates of predicted amino acid changes indicate that these two gene groups are evolving under different selective pressures and may be functionally divergent. This functional divergence was confirmed by comparing site-specific shifts in evolution indexes of the two groups of predicted proteins. The evidence of differential selection and possible functional divergence suggests that the duplication may be of adaptive significance, with possible implications for the explosive diversification of the Draba lineage during the cooling Quaternary stages and the following worldwide colonization of arid alpine and artic regions.

Trade-offs between constitutive and induced resistance in wild crucifers shown by a natural, but not an artificial, elicitor

Constitutive and induced plant resistance against herbivores occurs throughout the plant kingdom, but little is known about the evolutionary relationship between these two types of resistances. We examined the relationships between constitutive and induced resistance to the diamondback moth, Plutella xylostella, in 11 wild species of crucifers, and analyzed the changes in volatiles associated with their expression in two species. We used larvae of P. xylostella and jasmonic acid (JA) as elicitors of the induced response. The level of resistance was estimated as the relative number of eggs laid on the plants by P. xylostella. Substantial variation in constitutive resistance was observed among the 11 crucifer species. When the plants were damaged by larvae, a negative correlation was found between constitutive and induced resistance. However, a positive correlation was detected between constitutive and induced resistance when the plants were treated by JA. The shift in resistance was associated with changes in the emission of volatiles. These results strongly suggest that (1) a trade-off occurs between constitutive and induced resistance in wild crucifers, and that (2) such a trade-off can be observed by treating the plants with a natural, but not an artificial, elicitor.

The recent evolution of a pseudogene: diversity and divergence of a mitochondria-localized small heat shock protein inArabidopsis thaliana

In this study we examined the evolution of the genes for three organelle-localized small heat shock proteins in Arabidopsis thaliana : the chloroplast-localized (CP) protein HSP21 and two mitochondria-localized (MT) proteins, HSP23.5 and HSP23.6. We found that the CP protein and one of the MT proteins, HSP23.6, are evolving under purifying selection to maintain function. In contrast, the gene for HSP23.5, the other MT protein, is highly variable within A. thaliana, and in some accessions or ecotypes this gene may be a pseudogene. HSP23.5 and HSP23.6 are related via a segmental duplication event, and the presence of orthologs of each gene in other species within the Brassicaceae indicates that the duplication generating HSP23.5 and HSP23.6 may have occurred as much as 20 million years ago. This is considerably longer than the 4 million year half-life of gene duplicates (functional genes as well as pseudogenes) reported by some studies. Our results are consistent with the prediction that after gene duplication one gene duplicate can be maintained for some time under relaxed selection while it accumulates random mutations. By capturing a pseudogene in the making our study provides important information on how pseudogenes are formed.

An exploration of LEAFY expression in independent evolutionary origins of rosette flowering in Brassicaceae

Whereas most Brassicaceae produce flowers on an elongated inflorescence, a few lineages produce flowers directly from the vegetative rosette on elongated pedicels. Knowing the extent to which independent origins of rosette flowering involve the same developmental and genetic mechanisms could clarify the constraints acting on plant architectural evolution. Prior work in Idahoa, Ionopsidium, and Leavenworthia suggested that changes in the activity or expression of the flower meristem identity gene, LEAFY (LFY), played a role in all three origins of rosette flowering. Here we studied the developmental morphology of L. crassa and immunolocalization of LFY protein in Leavenworthia and Ionopsidium to further compare independent origins of rosette flowering. Leavenworthia crassa differs from Ionopsidium and Idahoa in producing ebracteate flowers. Flowers are, however, associated with "squamules," here interpreted as stipules of a cryptic bract. LFY was detected in L. crassa flower primordia but not in inflorescence meristems. In contrast, the rosette flowering Io. acaule accumulated LFY protein in the inflorescence meristem, whereas its inflorescence-flowering close relative, Io. prolongoi, did not. Thus, although different cases of rosette flowering likely entailed modifications of the same meristem identity program, distinct developmental genetic mechanisms appear to be involved in each case.

Different flooding responses in Rorippa amphibia and Rorippa sylvestris, and their modes of expression in F1 hybrids

The river floodplain species Rorippa amphibia, Rorippa sylvestris, and their hybrid Rorippa x anceps were studied here, with the aim of identifying potential species differences with respect to flooding tolerance, and of assessing their expression in F1 hybrids. Parents and their F1 hybrids were subjected to three flooding treatments mimicking natural conditions, and growth-related and leaf morphological traits were compared. In contrast to R. sylvestris, R. amphibia responded to waterlogging by forming specialized roots, and its growth was not reduced. These traits were dominantly expressed in hybrids. Both species and the hybrids established shoot growth over 2 wk of complete submergence. Only in R. sylvestris was this not at the expense of root biomass, suggesting that R. sylvestris can photosynthesize underwater. Rorippa sylvestris also showed a hyponastic response. Hybrids were intermediate to the parents in this respect. This study shows that phenotypic expression of parental traits in F1 hybrids is mostly additive, but can also be dominant. This suggests that a large overlap in habitat use of parents and hybrids is likely. If such an overlap occurs, the main evolutionary consequences of hybridization in Rorippa will be the introgression of genes, as the hybrids are fully fertile.

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 F(2) 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. Sequence data from this article can be found in the GenBank/EMBL data libraries under accession numbers DU 667459 to DU 708532.

Characterization of terminal-repeat retrotransposon in miniature (TRIM) in Brassica relatives

We have newly identified five Terminal-repeat retrotransposon in miniature (TRIM) families, four from Brassica and one from Arabidopsis. A total of 146 elements, including three Arabidopsis families reported before, are extracted from genomics data of Brassica and Arabidopsis, and these are grouped into eight distinct lineages, Br1 to Br4 derived from Brassica and At1 to At4 derived from Arabidopsis. Based on the occurrence of TRIM elements in 434 Mb of B. oleracea shotgun sequences and 96 Mb of B. rapa BAC end sequences, total number of TRIM members of Br1, Br2, Br3, and Br4 families are roughly estimated to be present in 660 and 530 copies in B. oleracea and B. rapa genomes, respectively. Studies on insertion site polymorphisms of four elements across taxa in the tribe Brassiceae infer the taxonomic lineage and dating of the insertion time. Active roles of the TRIM elements for evolution of the duplicated genes are inferred in the highly replicated Brassica genome.

The ABC's of comparative genomics in the Brassicaceae: building blocks of crucifer genomes

In this review we summarize recent advances in our understanding of phylogenetics, polyploidization and comparative genomics in the family Brassicaceae. These findings pave the way for a unified comparative genomic framework. We integrate several of these findings into a simple system of 24 conserved chromosomal blocks (labeled A-X). The naming, order, orientation and color-coding of these blocks are based on their positions in a proposed ancestral karyotype (n=8), rather than by their position in the reduced genome of Arabidopsis thaliana (n=5). We show how these crucifer building blocks can be rearranged to model the genome structures of A. thaliana, Arabidopsis lyrata, Capsella rubella and Brassica rapa. A framework for comparison between species is timely because several crucifer genome-sequencing projects are underway.

Horizontal gene transfer of a plastid gene in the non-photosynthetic flowering plants Orobanche and Phelipanche (Orobanchaceae)

Plastid sequences are among the most widely used in phylogenetic and phylogeographic studies in flowering plants, where they are usually assumed to evolve like non-recombining, uniparentally transmitted, single-copy genes. Among others, this assumption can be violated by intracellular gene transfer (IGT) within cells or by the exchange of genes across mating barriers (horizontal gene transfer, HGT). We report on HGT of a plastid region including rps2, trnL-F, and rbcL in a group of non-photosynthetic flowering plants. Species of the parasitic broomrape genus Phelipanche harbor two copies of rps2, a plastid ribosomal gene, one corresponding to the phylogenetic position of the respective species, the other being horizontally acquired from the related broomrape genus Orobanche. While the vertically transmitted copies probably reside within the plastid genome, the localization of the horizontally acquired copies is not known. With both donor and recipient being parasitic plants, a possible pathway for the exchange of genetic material is via a commonly attacked host.

Supernetwork Identifies Multiple Events of Plastid trnF(GAA) Pseudogene Evolution in the Brassicaceae

The occurrence of nonfunctional trnF pseudogenes has been rarely described in flowering plants. However, we describe the first large-scale supernetwork for the Brassiccaeae built from gene trees for 5 loci (adh, chs, matK, trnL-F, and ITS) and report multiple independent origins for trnF pseudogenes in crucifers. The duplicated regions of the original trnF gene are comprised of its anticodon domain and several other highly structured motifs not related to the original gene. Length variation of the trnL-F intergenic spacer region in different taxa ranges from 219 to 900 bp as a result of differences in pseudocopy number (1-14). It is speculated that functional constraints favor 2-3 or 5-6 copies, as found in Arabidopsis and Boechera. The phylogenetic distribution of microstructural changes for the trnL-F region supports ancient patterns of divergence in crucifer evolution for some but not all gene loci.

Brassicaceae contain nortropane alkaloids

The report of cochlearine, the 3-hydroxybenzoate ester of tropine found in Cochlearia officinalis, Brassicaceae, initiated a screening for tropane alkaloids in Cochlearia species and for calystegines in further Brassicaceae. All ten Cochlearia species investigated contained cochlearine, tropine, and pseudotropine. Calystegines, nortropane alkaloids deriving from pseudotropine, were also identified in all Cochlearia species and accumulated up to 0.5% dry mass in leaves. Brassicaceae species of all major lineages of the family were analysed for calystegines. Of the 43 species included in the study, 18 accumulated calystegines of various structures. This is the first screening of Brassicaceae for products of the tropane alkaloid pathway, which is known as characteristic for plants of Solanaceae family. The identification of calystegines in all branches of the Brassicaceae family including Aethionema, a species at the basis of the family, suggests tropane alkaloids as secondary compound typical for Brassicaceae.

Toward a Global Phylogeny of the Brassicaceae

The Brassicaceae is a large plant family (338 genera and 3,700 species) of major scientific and economic importance. The taxonomy of this group has been plagued by convergent evolution in nearly every morphological feature used to define tribes and genera. Phylogenetic analysis of 746 nrDNA internal transcribed spacer (ITS) sequences, representing 24 of the 25 currently recognized tribes, 146 genera, and 461 species of Brassicaceae, produced the most comprehensive, single-locus-based phylogenetic analysis of the family published to date. Novel approaches to nrDNA ITS analysis and extensive taxonomic sampling offered a test of monophyly for a large complement of the currently recognized tribes and genera of Brassicaceae. In the most comprehensive analysis, tribes Alysseae, Anchonieae plus Hesperideae, Boechereae, Cardamineae, Eutremeae, Halimolobeae, Iberideae, Noccaeeae, Physarieae, Schizopetaleae, Smelowskieae, and Thlaspideae were all monophyletic. Several broadly defined genera (e.g., Draba and Smelowskia) were supported as monophyletic, whereas others (e.g., Sisymbrium and Alyssum) were clearly polyphyletic. Analyses of ITS data identified several problematic sequences attributable to errors in sample identification or database submission. Results from parsimony ratchet and Bayesian analyses recovered little support for the backbone of the phylogeny, suggesting that many lineages of Brassicaceae have undergone rapid radiations that may ultimately be difficult to resolve with any single locus. However, the development of a preliminary supermatrix including the combination of 10 loci for 65 species provides an initial estimate of intertribal relations and suggests that broad application of such a method will provide greater understanding of relationships in the family.

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 < or = 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.